SciELO - Scientific Electronic Library Online

 
vol.68 issue6Half-metallic and thermodynamic properties of new d0 CsCaZ (Z= Ge, Sn and Pb) half Heusler alloys: A spin-based devices vision for the futureImpact of nonlinearity and wave dispersion parameters on the soliton pulses of the (2+1)-dimensional Kundu-Mukherjee-Naskar equation author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.68 n.6 México Nov./Dec. 2022  Epub July 31, 2023

https://doi.org/10.31349/revmexfis.68.061201 

Research

Nuclear Physics

Uncertainties in theoretical predictions for γdπ0d observables near threshold due to the use of different elementary amplitudes

H. M. Al-Ghamdia 

E. M. Darwishb  c  * 

A. A. Ibraheemd  e 

H. M. Abou-Elsebaab  c 

a Physics Department, College of Science, Princess Nourah bint, Abdulrahman University, P. O. Box 84428, Riyadh 11671, Saudi Arabia.

b Physics Department, Faculty of Science, Taibah University, Medina 41411, Saudi Arabia.

c Physics Department, Faculty of Science, Sohag University, Sohag 82524, Egypt.

d Physics Department, Faculty of Science, King Khalid University, Abha 62529, Saudi Arabia.

e Physics Department, Faculty of Science, Al-Azhar University, Assiut Branch, Assiut 71524, Egypt.


Abstract

We discuss possible uncertainties in theoretical predictions for γdπ0d observables near threshold due to the use of different elementary γNπN amplitudes using an approach which is based on time-ordered perturbation theory. Results are presented for unpolarized cross sections and all possible spin asymmetries of differential and total cross sections. Our results indicate that the estimations of the uncertainty on the γdπ0d observables show important sensitivity to the modeling of the elementary γNπN operator. A comparison to presently available experimental data is given. The results presented here are of particular interest for the evaluation of the systematic uncertainties caused by the use of different elementary operators in the analyses of γdπ0d measurements to extract information on the free neutron amplitude from deuteron data.

Keywords: Meson production; photoproduction reactions; few-body systems, deuteron; polarization phenomena in reactions; spin observables; polarized beams; polarized targets

1 Introduction

Understanding the internal structure and reactions of nuclear systems from first principles was, for a long time, an important goal of nuclear and particle physics. In this regard, few- and many-body systems provide a unique laboratory for studying nuclear interactions. Recent years have witnessed interesting developments in coherent and incoherent pion photoproduction on the deuteron and light nuclei, since relevant experimental studies have been performed (see Refs. [1-14] and references therein). These reactions can be utilized to investigate the internal structure of hadrons in the non-perturbative domain of Quantum Chromodynamics (QCD) and study the behavior of nucleon resonances in the nuclear medium.

In particular, the reaction γdπ0d may be used as an isospin filter and it is sensitive to the coherent sum of the proton and neutron amplitudes γpπ0p and γnπ0n, respectively. The γdπ0d reaction near π-threshold serves also as a test of our understanding of the chiral πN dynamics. The use of deuteron as an effective neutron target allows one to gain abundant information on the mechanisms of the elementary pion photoproduction off free neutron which otherwise is not possible due to the absence of any stable, dense, free-neutron targets. Furthermore, the deuteron represents also an ideal object for the study of NN interactions.

So far, many studies have been carried out dealing with the reaction γdπ0d in the photon energy region from π-threshold up to 1 GeV in order to obtain information on π0-photoproduction off free neutrons [15-33]. In spite of all these great theoretical efforts, a good description of the experimental data has not been yet found [3,4,6,7]. The disagreement between theoretical models and experimental data may be indicative of shortcomings in the elementary γNπN amplitudes. It was found in Refs. [34-39] that the γdπNN observables near 𝜋-threshold are considerably dependent on the elementary γNπN amplitude. During the last decades, more realistic models for the γNπN amplitude have been constructed [40-48]. Therefore, the construction of more realistic models for the investigation of the electromagnetic pion photo- and electroproduction processes on the deuteron and light nuclei remain challenging.

In addition, theoretical models which consider polarization observables for the reaction γdπ0d near π-threshold are very rare. Indeed, polarization observables constitute much more stringent tests of theoretical models. For instance, the helicity dependent total cross sections with respect to parallel (σP) and antiparallel (σA) spins of photon and deuteron provide important information on the nucleon spin structure. These cross sections are also required to test the validity of the Gerasimov-Drell-Hearn (GDH) sum rule [49]. Furthermore, the deuteron spin asymmetry σP — σP and the helicity E-asymmetry contain very interesting physics with respect to the internal hadron structure. Moreover, the double spin asymmetry T10c for circular polarized photons and vector polarized deuterons is also of great interest. It gives additional important information on the internal spin structure of the deuteron and provides a direct extraction of the GDH sum rule [49].

Therefore, the main aim of the present work is to discuss theoretical uncertainties in the analyses of γdπ0d observables near 𝜋-threshold due to the use of elementary γNπN amplitudes. Since the extraction of the elementary γpπ0p and γnπ0n amplitudes essentially relies on an interpretation in terms of the plane-wave impulse approximation (PWIA), it is necessary to study its validity, i.e., to have a quantitative estimate of all possible mechanisms going beyond the PWIA. However, in the present paper we want to restrict ourselves to the PWIA while a more realistic treatment including all possible reaction mechanisms will be reported in a forthcoming paper. Here, we want to explore the impact on the description of the γdπ0d observables from the choice of the pion photoproduction amplitudes off protons and neutrons. To the best of our knowledge, this has not been studied previously in the literature. We also investigate whether the uncertainties caused by the use of different elementary operators for the predictions of polarization observables in the γdπNN reaction channels [34-39] are also seen in the analyses of the coherent π0-photoproduction off deuterons, γdπ0d. It was demonstrated in Ref. [1] that the single- and double-polarization observables will allow us to select between different models for pion photoproduction on the nucleon.

Most recently, the sensitivity of the results for single-spin asymmetries and the helicity E-asymmetry in the γdπ0d reaction near threshold to the choice of the pion photoproduction amplitude has been investigated [33]. Unfortunately, the results presented in that work are based on an error in the definitions of initial and final relative momenta of the two nucleons in the deuteron process. Because of that the corresponding total energy of the photon-nucleon (γN) subsystem was not properly defined, but it strongly affects the predicted results especially at photon lab-energies close to π-threshold and forward pion angles. Thus, the importance of theoretical uncertainties in the γdπ0d observables still requires a more careful treatment as done in Ref. [33].

Therefore, the present work was motivated to present results for the unpolarized differential and total cross sections as well as for all possible beam, target, and beam-target spin asymmetries of the differential and total cross sections in the photon energy region from near π-threshold to 170 MeV. As elementary amplitudes, the one provided by the unitary isobar MAID-2007 model from [46] and those obtained using the Dubna-Mainz-Taipei dynamical model (DMT-2001) [47] and the chiral MAID model (χMAID-2013) [48] are used. Furthermore, we compare our results for the unpolarized differential cross section with the experimental data from TAPS [4]. In this comparison with experiment we concentrate our discussion on the unpolarized differential cross section, because data for polarization observables are still not available. The calculations presented in this paper are of particular interest for the evaluation of the systematic uncertainties caused by the use of different elementary operators in the analyses of γdπ0d measurements. These uncertainties should be kept in mind when extracting information on the free neutron amplitude from experimental data on the deuteron.

This paper is structured as follows. In Sec. 2, we briefly outline the theoretical framework used to describe the γdπ0d reaction. We give the explicit expressions for the γdπ0d observables in terms of the reduced t-matrix elements. Subsequently, we calculate the scattering amplitude for this reaction based on time-ordered perturbation theory. Section 3 is devoted to the main results together with a comparison with the available experimental data. Finally, we summarize our results in Sec. 4.

2 Theoretical framework

In this section we briefly describe the theoretical framework for the coherent π0-photoproduction reaction on the deuteron. We would like to mention that the general formalism for the γdπ0d reaction has been described in Ref. [29], and we refer the reader to that work for further details. Here, we briefly recall the necessary notation and definitions.

2.1 Cross section and observables

We consider the reaction γdπ0d in the photon-deuteron (γd) center-of-momentum (c.m.). We use a coordinate system with the z-axis along the photon momentum kez=k^=k/k, the y-axis parallel to k×q, and the x-axis in the direction of maximal linear photon polarization [50], in which the pion momentum q has spherical coordinates θ and ϕ. If the incoming photon beam is not linearly polarized, then the x-axis may be chosen arbitrary, and there is no dependence on the angle ϕ. For the deuteron states, a non-covariant normalization following to the conventions of Ref. [51] is used. The initial and final total three-momenta of the two nucleons in the deuteron are given in the c.m. system by P=p1+p2=-k and P'=p '1+p '2=-q, respectively. For the initial and final relative three-momenta of the two nucleons in the deuteron in the c.m. system, we obtain pr=p1-p2/2=p-k/2 and p 'r=p '1-p '2/2=p-q/2, respectively (see Fig. 1).

Figure 1 Diagrammatic representation of the γdπ0d reaction in the impulse approximation with definition of momenta in the γd c.m. system. The spectator nucleon is on-shell. 

The cross section for arbitrary polarized photons and initial deuterons can be calculated for a given transition M-matrix by applying the density matrix formalism [50]. All different cross sections for the γdπ0d reaction have the form:

O=Km'~dλ~m~dm'dλmdMm'~dλ~m~d*Ωm'~dm'dMm'dλmdρλλ~γ ρmdm~dd , (1)

where K is a kinematic factor, Mm'dλmd the scattering matrix, ργd) the density matrix for incoming photon (initial deuteron) polarization, mdm'd the spin projection of the initial (final) deuteron, and λ ± 1 the circular photon polarization. A polarimeter for the deuteron in the final state is described by the operator Ω. In the present work, we set Ω = 1, because we do not consider any polarization analysis of the final deuteron. The scattering Mm'dλmd-matrix is given by isolating the azimuthal dependence as follows [17,22]:

Mm'dλmd(θ,ϕ)=ei(λ+md)ϕ tm'dλmd(θ) , (2)

where the reduced t-matrix elements are defined by separating the ϕ-dependence from the M-matrix elements. Parity conservation gives for the t-matrix the symmetry:

t-m'd-λ-md=(-1)1+m'd+λ+md tm'dλmd . (3)

The kinematic factor K is given in the γd c.m.frame by:

K=116π2EdE'dWγd2 qk , (4)

where the deuteron energies in the initial and final states are given by Ed=Md2+k2 and E'd=Md2+q2, respectively, with Md as deuteron mass. The absolute values of the photon and pion three-momenta in the γd c.m. frame are given, respectively, by:

k=12Wγd (Wγd2-Md2)andq=12Wγd Wγd2-(Md-mπ)2Wγd2-(Md+mπ)2 . (5)

The invariant energy of the γd system is given as:

Wγd=Eγ+Ed=k+Md2+k2 ,=Eπ+E'd=mπ2+q2+Md2+q2 ,=Md2+2MdEγ , (6)

Where mπ is the neutral-pion mass.

Following the rules of Ref. [52], we write the differential cross section in terms of the unpolarized differential cross section dσ0/dΩ and the various spin asymmetries Σ, TIM, TIMc, and TIMl as follows:

dσdΩ=dσ0dΩ1+Plγ {Σ cos2ϕ+I=12PId M=-IITIMlcosψM-δI1 π2 dM0I(θd)}+I=12PId M=0I(TIMcosM(ϕd-ϕ)-δI1 π2+Pcγ TIMcsinM(ϕd-ϕ)+δI1 π2) dM0I(θd), (7)

where ψM=M (ϕd-ϕ)+2ϕ and dM0I(θd) is a small rotation matrix for which the convention of Ref. [53] is used. The photon polarization is characterized by the degree of circular Pcγ=Pzγ and linear Plγ=(Pxγ)2+(Pyγ)2 polarization, where the x-axis has been chosen in the direction of maximum linear polarization, i.e. Pxγ=-Plγ and Pyγ=0. The deuteron is characterized by the vector and tensor polarization components P1d and P2d, respectively, and the orientation angles (θd,ϕd) of the orientation axis of the deuteron with respect to which the density matrix of the deuteron has been assumed to be diagonal; the elements with mdmd' will therefore vanish.

In order to express the unpolarized differential cross section and various polarization observables in terms of the t-matrix, it is appropriate to define the two quantities [50]:

VIM=K3 2I+1 mdmd'(-)1-md11Imd'-mdMmd''tmd''1md'* tmd''1md , (8)

WIM=-K3 2I+1 mdmd'(-)1-md11Imd'-mdMmd''tmd''1md'* tmd''-1md , (9)

where the convention of Edmonds [53] is used for the Wigner 3j-symbol. The explicit formal expressions for unpolarized differential cross section and spin asymmetries can be derived as follows [28]:

(i) The unpolarized differential cross section:

dσ0dΩ=V00=K3 m'dmd|tm'd1md|2 . (10)

(ii) The linear photon asymmetry:

Σ dσ0dΩ=W00=-K3 m'dmdtm'd1md* tm'd-1md . (11)

(iii) The vector target asymmetry:

T11 dσ0dΩ=2 Im V11=23 K Im md(tmd1-1* tmd10+tmd10* tmd11) . (12)

(iv) The tensor target asymmetries:

T2M dσ0dΩ=(2-δM0) Re V2Mfor0M2,

with

T20 dσ0dΩ=K32 md(|tmd11|2+|tmd1-1|2-2 |tmd10|2) , (13)

T21 dσ0dΩ=23 K Re md (tmd1-1* tmd10-tmd10* tmd11) , (14)

T22 dσ0dΩ=2K3 Re md tmd1-1* tmd11 . (15)

(v) The spin asymmetries for circularly polarized photons and vector polarized deuterons:

T1Mc dσ0dΩ=-(2-δM0) Re V1M ,for0M1

with

T10c dσ0dΩ=K6 md(|tmd11|2-|tmd1-1|2) , (16)

T11c dσ0dΩ=-23 K Re md(tmd1-1* tmd10+tmd10* tmd11) . (17)

(vi) The spin asymmetries for circularly polarized photons and tensor polarized deuterons:

T2Mc dσ0dΩ=-2-δM0Im V2M ,for0M2

with

T21c dσ0dΩ=23 K Im md(tmd10* tmd11-tmd1-1* tmd10) , (18)

T22c dσ0dΩ=-2K3 Im mdtmd1-1* tmd11 . (19)

Because the quantity VI0 is real according to VIM*=(-)M VI-M under complex conjugation, the spin asymmetry T20c vanishes identically.

(vii) The spin asymmetries for linearly polarized photons and vector polarized deuterons:

T1Ml dσ0dΩ=i W1M,for-1M1,

with

T10l dσ0dΩ=23 K Im mdtmd11* tmd-11 , (20)

T11l dσ0dΩ=-23 K Im mdtmd1-1* tmd-10, (21)

T1-1l dσ0dΩ=23 K Im mdtmd11* tmd-10 . (22)

(viii) The spin asymmetries for linearly polarized photons and tensor polarized deuterons:

T2Ml dσ0dΩ=W2Mfor-2M2,

with

T20l dσ0dΩ=23 K Re md(tmd10* tmd-10-tmd11* tmd-11) , (23)

T21l dσ0dΩ=23 K Re mdtmd10* tmd-11 , (24)

T2-1l dσ0dΩ=23 K Re mdtmd10* tmd-1-1 , (25)

T22l dσ0dΩ=-K3 mdtmd1-1* tmd-11 , (26)

T2-2l dσ0dΩ=-K3 mdtmd11* tmd-1-1 . (27)

We would like to mention that for θ = 0 and π the spin asymmetries TIMc=0 for M0 and TIMl=0 for M2 because in that case the angle ϕ is undefined or arbitrary and, therefore, the differential cross section cannot depend on ϕ.

The deuteron spin asymmetry with respect to circularly polarized photons and the deuteron spin oriented parallel (P) and antiparallel (A) to the photon spin is related to the spin asymmetry T10c according to [54]:

d(σP-σA)dΩ=6 dσ0dΩ T10c . (28)

As mentioned in the Introduction, the asymmetry T10c is of special interest, because it is related to the spin asymmetry σP — σA which determines the GDH sum rule [49].

The general form of the total cross section with inclusion of photon and deuteron polarization effects is obtained from Eq. (7) by integrating dσ/dΩ over the pion spherical angle dΩ and reads [54]:

σ(Plγ,Pcγ,P1d,P2d)=σ0 1+P2d T~20 12(3cos2θd-1)+Pcγ P1d T~10 c cosθd+Plγ P2d T~22 lcos(2ϕd) 64sin2θd , (29)

where the unpolarized total cross section σ0 and the corresponding spin asymmetries σ0 T~20, σ0 T~10c and σ0 T~22l are given by:

σ0=dΩ dσ0dΩ , (30)

σ0 T~20=dΩ dσ0dΩ T20 , (31)

σ0 T~10c=dΩ dσ0dΩ T10c , (32)

σ0 T~22l=dΩ dσ0dΩ T22l . (33)

2.2 The γdπ0d amplitude

Next, the transition matrix elements Mm'dλmd are calculated in the frame of time-ordered perturbation theory. The impulse approximation (IA), in which the reaction will take place only on one of the two nucleons in the deuteron leaving the other as a pure spectator, usually serves as the starting point to calculate the amplitude for electromagnetic pion production on the deuteron [16] or, in general, on a nucleus. It corresponds to a direct embedding of the elementary γNπN amplitudes into the two-nucleon system. In our framework, we restrict the calculation of the transition matrix elements Mm'dλmd to IA in order to study the dependence of results for the γdπ0d observables on the elementary pion photoproduction operator.

The production operator tγπd for the deuteron process is obtained from the elementary operator tγπ by:

tγπd=tγπ(1)2+(1)tγπ(2) . (34)

The upper index refers to the nucleon on which the elementary operator acts. This means that tγπd contains pure single nucleon terms. The transition M-matrix of coherent π0-photoproduction on the deuteron has the following form:

Mm'dλmd(k,q)=2d3p(2π)3 ψm'dp-q2× p-q|tγπ(1)|p-k ψmdp-k2 , (35)

where ψmd(p) denotes the deuteron wave function and tγπ(1) the elementary γNπN operator.

For the intrinsic part of the deuteron wave function we use the ansatz:

ψmd(p )=L=0,2mLmS(LmL1mS|1md)uL(p) YLmL(p^) χmS ζ0 . (36)

The last two terms χmS and ζ0 denote spin and isospin wave functions, respectively. The S and D components of the deuteron wave function (DWF) are given by u0(P) and u2(P), respectively. In the present work, we compute the radial deuteron wave functions uL(p) using the realistic and high-precision Bonn full model [55].

For the elementary γNπN amplitude, we take the π-production operator from the unitary isobar MAID-2007 model introduced in Ref. [46]. This model is based on Born terms, ρ and ω vector-meson exchange contributions, and 13 four-star nucleon resonance excitations. It describes well the elementary γNπN amplitude and agrees with experimental observations very well. The MAID-2007 model is parameterized in terms of invariant amplitudes and therefore allows one to evaluate the transition deuteron amplitude in any frame of reference.

To study the uncertainties caused by the use of different elementary pion photoproduction operators in the analyses of γdπ0d observables, the dynamical DMT-2001 model [47] and the χMAID-2013 model [48] are used. The DMT-2001 model is a unitary dynamical model based on a non-resonant background described by Born terms and vector-meson exchange contributions in the t-channel (ρ and ω) and the following 8 four-star nucleon resonances in the s-channel: P33(1232), P11(1440), D13(1520), S11(1535), S31(1620), S11(1650), F15(1680) and D33(1700). The DMT-2001 model describes the pion photo- and electroproduction observables in terms of photon and nucleon degrees of freedom and provides a very good description of experimental data in the near-threshold region. The χMAID-2013 model for studying pion photo- and electroproduction on the nucleon in the near-threshold region has been constructed in relativistic chiral perturbation theory and included all multipole amplitudes up to and including l=4 or, in other words, G waves. This model provides a comprehensive description of the elementary process on the free nucleon up to photon energies of 170 MeV in the laboratory frame. The details of each of the used models (MAID-2007, DMT-2001, and χMAID-2013) can be found in their original works in Refs. [46-48] and therefore will not repeated here.

3 Results and discussions

In this section, we explore the sensitivity of the results for unpolarized cross sections as well as all various beam, target, and beam-target spin asymmetries of the γdπ0d reaction near threshold to the elementary pion photoproduction amplitude. For this purpose, we use for the elementary reaction amplitudes three different realistic models. These models are the MAID-2007 [46], DMT-2001 [47], and χMAID-2013 [48]. For the deuteron wave function, we use the realistic and high-precision Bonn full model [55]. We discuss the energy and angular dependences of the results for these observables and give a comparison of the obtained results with the available experimental data. In all the upcoming figures, the dashed, dotted, and solid curves represent the results using MAID-2007 [46], DMT-2001 [47], and χMAID-2013 [48] models for the elementary operator, respectively.

3.1 Differential and total cross sections

First, we show in Fig. 2 the energy and angular dependences of the results for unpolarized differential cross section, dσ0/dΩ, using different elementary γNπN amplitudes. We see that the dσ0/dΩ results with different elementary amplitudes are quite different especially at forward pion angles. The differences between the results with various elementary amplitudes decrease with increasing pion angle until they become very small at θ = 180o. When the dashed curve (MAID-2007) is compared with both the dotted (DMT-2001) and solid (χMAID-2013) curves, one can see that these differences are very obvious at forward pion angles and show up the discrepancies among elementary amplitudes. In this case, the calculated dσ0/dΩ within the MAID-2007 model is smaller than those within DMT-2001 and χMAID-2013. At extreme backward pion angles, we see that the curves represent the results of dσ0/dΩ using MAID-2007, DMT-2001, and χMAID-2013 are very close to each other and thus the influence of elementary operators on dσ0/dΩ is negligible in this case. This means that the dσ0/dΩ results are sensitive to the choice of the elementary amplitude at forward pion angles.

Figure 2 (Color online) The differential cross section for the reaction γdπ0d using different elementary operators and the DWF from Bonn full potential. The dashed, dotted, and solid curves represent the results using the MAID-2007, DMT-2001, and χMAID-2013 models for the elementary operator, respectively. Results at Eγ = 144 MeV are multiplied by the factor in the parentheses. 

Figure 3 shows the results for the unpolarized total cross section, σ0, for the reaction γdπ0d as a function of the photon energy in the laboratory system, Eγ, using different elementary operators. The importance of the choice of elementary amplitude is clearly addressed when the calculation with the MAID-2007 (dashed curve) is compared to the DMT-2001 (dotted curve) and χMAID-2013 (solid curve) curves. It is very clear that the results using MAID-2007 model differ significantly from the DMT-2001 and χMAID-2013 ones. However, the σ0 results using DMT-2001 and χMAID-2013 models are close to each other. We would like to point out that the MAID-2007 model is based on a single channel approach, whereas the DMT-2001 model is based on coupled channels. This is an important difference between both models.

Figure 3 (Color online) The unpolarized total cross section for the reaction γdπ0d as a function of photon lab-energy using different elementary operators and the DWF from Bonn full NN potential. Curve conventions as in Fig. 2

3.2 Single-spin asymmetries

Now, we focus our attention on the single-spin asymmetries of polarized photons (Σ) or polarized deuterons (T11, T20, T21 and T22) for γdπ0d and γdπ0d, respectively. We present in Fig. 4 the results for the linear photon asymmetry Σ at different photon lab-energies as a function of emission pion angle θ in the γd c.m. frame. In general, one can see the results for Σ remain same qualitatively but quantitatively slightly changed. Figure 4 shows that the photon Σ-asymmetry has dominately negative values in the photon energy domain of the present work. One displays also that the Σ-asymmetry decreases with increasing pion scattering angle until a minimum close to θ ≃ 130o at Eγ = 144 MeV is reached. This minimum is shifted towards lower pion angles with increasing the photon lab-energy. Then the photon asymmetry Σ increases with increasing pion angle until it reaches zero at θ = π. It is also noticeable from Fig. 4 that the photon asymmetry Σ vanishes at θ = 0o and 180o, because in that case the differential cross section cannot depend on the azimuthal angle ϕ, since at θ = 0o and π the angle ϕ is undefined or can be arbitrary. At extreme forward and backward emission pion angles, one notes that the photon Σ-asymmetry is relatively small in comparison to the results when θ varies from 60o to 150o.

Figure 4 (Color online) Same as in Fig. 2 but for the photon Σ-asymmetry with linearly polarized photons and unpolarized deuterons. 

We found that the asymmetry Σ is sensitive to the elementary amplitude, especially in the peak region where sizeable differences are obtained in the pion angle range from 60o to 150o. It has a minimal value around θ ≃ 135o in the case of χMAID-2013 (solid curve) and it is shifted towards lower pion angles in the case of DMT-2001 (dotted curve) and MAID-2007 (dashed curve). It is also very obvious that the computations with different elementary amplitudes are quite different with, in absolute size, a larger Σ-asymmetry predicted using χMAID than the ones obtained with DMT and MAID models. This discrepancy shows up the differences among elementary pion photoproduction operators and means that Σ is very sensitive to the choice of the elementary amplitude, in particular in the vicinity of the peak.

For polarization observables with polarized deuteron targets and unpolarized photon beams, we present in Fig. 5 the sensitivity of the results for the vector T11 and tensor T2M (M = 0,1,2) deuteron spin asymmetries for the reaction γdπ0d to the choice of elementary pion photoproduction amplitude. We see that the results for T11 asymmetry exhibit qualitatively similar behaviors for different elementary pion photoproduction operators. The T11 asymmetry is sensitive to the imaginary parts of the scattering amplitudes and its values vanish identically at θ = 0 and π. This asymmetry depends on the relative phase of the matrix elements as can be seen from Eq. (12). It would vanish for a constant overall phase of the t-matrix. One can also see that the results using different elementary operators are quantitatively rather different even at forward and backward pion angles. This discrepancy displays the differences among elementary pion photoproduction operators which means that the T11-asymmetry is sensitive to the choice of the elementary amplitude.

Figure 5 (Color online) Same as in Fig. 2 but for the deuteron asymmetries with vector T11 and tensor T2M (M = 0, 1, 2) polarized deuterons and unpolarized photons. 

The results for tensor deuteron spin asymmetries T2M (M = 0,1,2) are also shown in Fig. 4. In general, we see that the results for these asymmetries at the lowest photon energy, Eγ = 144 MeV, are rather different than the results at higher photon energies, especially in the case of T22 asymmetry. The results for T20 and T21 asymmetries using various elementary amplitudes exhibit qualitatively, but not quantitatively, similar behaviors. One can see that the results using different elementary operators are quantitatively rather different, which means that these spin asymmetries are slightly sensitive to the choice of the elementary amplitude. When the photon energy increases, we see that the curves represent the results of T20 and T21 using various elementary amplitudes are close to each other and thus the influence of elementary operators on T20 is small in this case. The T22 asymmetry indicates at photon energies greater than 144 MeV that it has an oscillatory shape. In contrast to the T20 and T21 cases, we find that the results for T22 asymmetry using different elementary amplitudes exhibit qualitatively and quantitatively different behaviors and show up the sensitivity of its results to the elementary amplitude.

A serial system of four columns packed with Eichhornia crassipes and 16 cm radio and 150 cm length is capable of keeping an outlet concentration below 10 ppb

3.3 Beam-target double-spin asymmetries

Here, we report the numerical results for the beam-target double spin asymmetries for γdπ0d reaction near threshold. We start with the results for T10c and T11c asymmetries with circular polarized photons and vector polarized deuterons as shown in Fig. 6 as functions of pion angle in the γd c.m. frame at the same photon lab-energies as mentioned before in the case of single-spin asymmetries. We see that the results for T10c and T11c at different photon lab-energies remain same qualitatively but quantitatively slightly changed. The T10c and T11c asymmetries behave the same as T20 and T21 ones, respectively. Figure 6 displays that the T10c asymmetry has values between —1 and 1. It begins negative at θ = 0o and increases with increasing pion angle until a maximum value at θ ≃ 90o. Then, it rapidly falls down to negative values when increasing the pion angles. The maximum value of T10c is shifted towards higher pion angles at Eγ = 144 MeV. As mentioned in the Introduction, the spin asymmetry T10c determines the GDH sum rule [49]. The T11c results vanish at θ = 0o and π and indicate that it has an oscillatory shape. The T10c and T11c results using different elementary operators are quantitatively rather different at the lowest photon energy, which means that these asymmetries are sensitive to the choice of the elementary amplitude close to threshold. By increasing photon energies, the curves represent the results of T10c and T11c using different pion production amplitudes are close to each other and thus a small dependence of these asymmetries on the elementary operators is obtained.

Figure 6 (Color online) Same as in Fig. 2 but for the beam-target double-spin asymmetries with circularly polarized photons and vector polarized deuterons. 

Figure 7 shows the results for the beam-target double spin asymmetries T21c and T22c for circularly polarized photons and tensor polarized deuterons. A quick glance reveals that the results for T21c and T22c are negative and vanish at θ = 0o and π. In contrast to the vector spin asymmetries T10c and T11c, it is obvious from Fig. 7 that the tensor spin asymmetries T21c and T22c are much more sensitive to the choice of elementary amplitude. Figure 7 shows also that the sensitivity of T21c and T22c results to the elementary pion photoproduction amplitude is very important, especially in the pion angle range between 30o and 120o.

Figure 7 (Color online) Same as in Fig. 2 but for the beam-target double-spin asymmetries with circularly polarized photons and tensor polarized deuterons. Results for the T22c asymmetry are multiplied by the factor in the parentheses. 

In Fig. 8 we present our results for T10l, T1+1l, and T1-1l double spin asymmetries with longitudinal polarized photons and vector polarized deuterons. In general, we see that these double spin asymmetries are sensitive to the elementary amplitude, especially in the pion angle range between 30o and 150o. The asymmetry T10l is smaller than the other ones and it differs in size between the results with different elementary amplitudes. This emphasizes the very important dependence of the T10l results on the elementary amplitude. As discussed in Refs. [31,38], the vector asymmetries T1Ml are considerably small for π0 production near threshold. Figure 8 displays also that the difference between the results for T10l and T1-1l asymmetries using χMAID-2013 model and both DMT-2001 and MAID-2007 increases with increasing photon lab-energy, which indicates that these spin asymmetries are very sensitive the choice of the elementary γNπN amplitude.

Figure 8 (Color online) Same as in Fig. 2 but for the beam-target double-spin asymmetries with linearly polarized photons and vector polarized deuterons. Results for the T1-1l asymmetry are multiplied by the factor in the parentheses. 

The results for the beam-target double spin asymmetries T20l, T2+1l, T2+2l, T2-1l, and T2-2l with longitudinal polarized photons and tensor polarized deuterons are shown in Fig. 9. We see that these spin asymmetries are also affected by the choice of various elementary amplitudes, especially at photon energies close to threshold. The sensitivity of the results for these asymmetries to the elementary amplitude is very obvious at the lowest photon energy, Eγ = 144 MeV, and decreases with increasing photon energy. This sensitivity is much more obvious in the case of T20l and T2-1l asymmetries since the differences between the results with various elementary operators are significant, in particular in the peak region. The T2±2l asymmetries are sensitive to both the real and the imaginary parts of the scattering amplitudes.

Figure 9 (Color online) Same as in Fig. 2 but for the beam-target double-spin asymmetries with linearly polarized photons and tensor polarized deuterons. Results for the T2-1l and T2-2l asymmetries are multiplied by the factor in the parentheses. 

Figure 10 shows the results for the beam-target double-spin asymmetries σ0 T~20, σ0 T~10c and σ0 T~22l of the total cross section as functions of the photon lab-energy. We see that the results for these asymmetries with different elementary amplitudes are quite different. The differences between the results with various elementary amplitudes increase with increasing photon lab-energy. When the dashed curve (MAID-2007) is compared with both the dotted (DMT-2001) and solid (χMAID-2013) curves, one can see that these differences are obvious and show up the discrepancies among elementary amplitudes. The calculated σ0 T~20, σ0 T~10c, and σ0 T~22l within MAID-2007 model is, in absolute size, smaller than those within DMT-2001 and χMAID-2013. At photon energies close to threshold, we see that the curves represent the results of σ0 T~20, σ0 T~10c, and σ0 T~22l using MAID-2007, DMT-2001, and χMAID-2013 models are close to each other and thus the influence of elementary operator on these spin asymmetries is negligible in this case.

Figure 10 (Color online) Same as in Fig. 3 but for the beam-target double-spin asymmetries σ0 T~20, σ0 T~10c and σ0 T~22l of the total cross section. 

3.4 Helicity-dependent cross sections and E-asymmetry

Next, we present in Fig. 11 the results for the doubly polarized differential cross sections for parallel (upper part) and antiparallel (middle part) spins of photon and deuteron as functions of pion angle in the c.m. frame at various photon lab-energies using different elementary amplitudes. These polarized differential cross sections are given by

dσPdΩ=dσ0dΩ 1+12T20+32T10c , (37)

dσAdΩ=dσ0dΩ 1+12T20-32T10c . (38)

Figure 11 (Color online) Same as in Fig. 2 but for the doubly polarized differential cross sections for parallel (upper part) and antiparallel (middle part) spins of photon and deuteron and their difference (lower part). Results at Eγ = 144 MeV are multiplied by the factor in the parentheses. 

In Fig. 11 we also present an important physical observable which is the difference dP — σA)/dΩ (lower part) that has a direct relation with the GDH sum rule [49]. The results displayed show the contribution of each dσP/dΩ and dσA/dΩ on the difference. Since their values have opposite behavior with increasing θ, we can see that the difference has negative values at small pion angles 0o < θ < 30o because of dσA/dΩ has the larger contribution in this range of pion angles. Thereafter, for θ ranges between 30o and 90o, dP — σA)/dΩ starts to increase and behaves the same as dσP/dΩ. For θ = 90o, the values of dP — σA)/dΩ start to decrease due to dσA/dΩ which has larger contribution again.

We would like to mention that the threshold region is dominated by the pion production with spin 1/2, which makes the antiparallel term dσA/dΩ larger than the parallel one dσP/dΩ. The estimated values of dP — σA)/dΩ close to threshold at Eγ = 144MeV are rather small compared to their values at higher energies under consideration.

As for the dσP/dΩ results (upper part in Fig. 11), we observe that the results are sensitive to the choice of elementary amplitude at the peak region since we obtain smaller values using MAID-2007 than using DMT-2001 and χMAID-2013. At extreme forward and backward pion angles, similar results are obtained and the sensitivity of dσP/dΩ results to the elementary amplitude is negligible. In the case of dσA/dΩ results (middle part in Fig. 11), we see that the influence of dσA/dΩ on the elementary amplitude is similar to the case of unpolarized differential cross section. One can see that the differences among dσA/dΩ results using different elementary amplitudes are very obvious at forward pion angles and the calculation within the MAID-2007 is smaller than those within DMT-2001 and χMAID-2013 but provides similar results at extreme pion backward angles. This discrepancy shows up the differences among elementary amplitudes. The computations of the differential spin asymmetry with respect to circularly polarized photons and oriented deuterons, dP — σA)/dΩ, (lower part in Fig. 11) shows that the sensitivity to the choice of elementary amplitude is also important in the peak position. This means that the difference dP — σA)/dΩ is also sensitive to the choice of the elementary amplitude.

Figure 12 shows the results for the helicity-dependent total cross sections for the γdπ0d reaction as functions of photon lab-energies using different elementary amplitudes. Results are displayed as follows: (left panel) total cross section σP for circularly polarized photons on a deuteron target with spin parallel to the photon spin; (middle panel) σA, the same for antiparallel spins of photon and deuteron; (right panel) spin asymmetry σP — σA. We see that the cross sections σP and σA as well as the deuteron spin asymmetry σP — σA present similar behaviors for all elementary amplitudes used in the present work. Compared the calculation with the MAID-2007 (dashed curve) to the DMT-2001 (dotted curve) and χMAID-2013 (solid curve) ones, we clearly address the importance of the choice of elementary amplitude. We see that the results using MAID-2007 model differ from the DMT-2001 and χMAID-2013 ones. We also find that the results for σP — σA starts out negative due to the E0+ multipole, which is dominant in the threshold region. The influence of elementary amplitude is important in σP, σA, and σP — σA. We would like to emphasize here that several experiments to measure the deuteron spin asymmetry σP — σA are presently underway [3,8,9].

Figure 12 (Color online) Same as in Fig. 3 but for the doubly polarized total cross sections for parallel (left panel) and antiparallel (middle panel) spins of photon and deuteron and their difference (right panel). 

During the recent years, there is a considerable interest in experiments [1-7] to measure the double polarization observable E for the γdπ0d reaction. This asymmetry is given by

E(θ)=d(σA-σP)/dΩd(σA+σP)/dΩ=d(σA-σP)/dΩ2dσ0/dΩ . (39)

In Fig. 13 we present the results for E-asymmetry as a function of the emission pion angle θ in the γd c.m.frame at four fixed values of the incident photon lab-energy. We see that the E-asymmetry has qualitatively a similar behavior for all incident photon lab-energies considered. Its maximum equals unity at θ = 0o and 180o. The curves begin with unity and decrease as the pion angle increases until a minimum value at θ ≃ 120o is reached. Then, it increases again to unity. The minimum value is shifted towards lower pion angles with increasing Eγ. The negative values in the E-asymmetry come mainly from higher positive contribution in dσP/dΩ. As mentioned in Refs. [22,34], the beam-target double polarization E-asymmetry is an excellent observable to test any weakness in the underlying elementary pion photoproduction model. Figure 13 shows that the helicity E-asymmetry is sensitive to the choice of the elementary γNπN amplitude at incident photon lab-energies close to threshold. When the photon energy increases, we see that the sensitivity of the results for E-asymmetry to the choice of the pion photoproduction operator is rather small.

Figure 13 (Color online) Same as in Fig. 3 but for the double polarization E-asymmetry. 

3.5 Comparison with experimental data

We close this section by comparing our results with the available experimental data. In fact, several experiments to measure both single and double polarization observables for coherent π0-photoproduction on the deuteron are presently underway. Unfortunately, no data are currently available for these polarization observables in the kinematic region under consideration in the present work and therefore we cannot make a comparison to experimental data for polarization observables. Thus, in the present work we compare the calculated results for the unpolarized differential cross section with the available data.

Figure 14 shows a comparison between our results for the unpolarized differential cross section dσ0/dΩ at two values of photon lab-energies (Eγ = 151.4 and 171.8 MeV) and the experimental data from TAPS [4]. One readily sees, that the estimated results for dσ0/dΩ using different elementary amplitudes underestimate the experimental data at backward pion angles. The results with the χMAID-2013 model (solid curve) are the nearest one to the experimental data in this case. At forward pion angles and Eγ = 151.4 MeV, an overestimation of the results using χMAID-2013 and DMT-2001 elementary amplitudes is found. But the results for dσ0/dΩ using MAID-2007 model underestimate the last three data points at high angles. At Eγ =171.8 MeV and for forward pion angles, a good agreement with the experimental data is also obtained. In this case, the shape given by MAID model is consistent with data lacking only in strength. We would like to mention that the elementary operators considered in the present work use some values of resonance photocouplings which especially in case of neutron are not yet well understood, and fitting these values to the data could improve the agreement between theory and experiment.

Figure 14 (Color online) The unpolarized differential cross section for γdπ0d in comparison with the experimental data from TAPS [4]. Curve conventions as in Fig. 2

We would like to point out that even at the PWIA level studied in the present work, there are also uncertainties due to the deuteron wave function, since one is involving large momentum components of the deuteron. Therefore, in Fig. 15 we present results for dσ0/dΩ using different deuteron wave functions in comparison with the experimental data from TAPS [4]. For this purpose, we use the deuteron wave functions of the Bonn full [55] (solid curve), CD-Bonn [56] (dashed curve), and Paris [57] (dotted curve) NN potentials. These potentials are exceedingly employed for numerical estimations of electromagnetic reactions on the deuteron, give a precise characterization of the NN scattering data and phase shifts, and used to characterize the range of the NN interaction.

Figure 15 (Color online) The unpolarized differential cross section for γdπ0d using χMAID-2013 elementary amplitude in comparison with the experimental data from TAPS [4]. The dotted, dashed, and solid curves represent the results using the Paris, CD-Bonn, and Bonn-full NN potentials for the deuteron wave function, respectively. 

It is clear from Fig. 15 that the results using the deuteron wave function of the Paris potential is smaller than those using CD-Bonn and Bonn full potentials. At Eγ = 151.4 MeV, we note an overestimation of the results using CD-Bonn and Bonn full potentials at θ < 60o, whereas a good agreement with the experimental data using the Paris potential is obtained in this case. On the contrary, a good agreement between the results using CD-Bonn and Bonn full potentials and the experimental data is obtained at Eγ = 171.8 MeV and forward pion angles. The results using the Paris potential underestimate the experimental data in this case. At backward pion angles, we see that the results using various NN potentials underestimate the experimental data. The origin of the differences obtained using various realistic deuteron wave functions maybe due to the tensor force between two nucleons. It is well-known that a practical measure for the strength of the tensor-force component contained in a nuclear potential is the predicted D-state probability of the deuteron, PD (see, for example, Refs. [58,59]). The PD values for the NN potentials used in this work are 4.85% for CD-Bonn, 4.25% for Bonn full, and 5.77% for Paris. The difference in these values of PD is related to the different tensor part of the NN potentials. The dependence of the γdπ0d and ede'd' observables on the D-state component of the deuteron wave function was investigated in Refs. [20] and [60], respectively. It was found that the D-wave contribution becomes visible at backward scattering angles.

The obtained discrepancies between our estimations for dσ0/dΩ and the experimental data can be attributed to the neglected contributions from two-body effects in the transition M-matrix. It was found in Refs. [16,61] that the influence of intermediate first-order pion rescattering and the two-loop diagram which includes πN and NN rescattering in the intermediate state are very important at photon energies near threshold. In addition, the meson-exchange current effects were found to be quite significant for coherent π+ photoproduction on 3He in Ref. [62].

From the preceding discussions it is apparent that the choice of the elementary amplitude has a visible effect on unpolarized cross sections as well as on various beam, target, and beam-target spin asymmetries in the γdπ0d reaction near threshold. Summarizing, we can say that the MAID-2007 model provides different predictions for all possible observables in γdπ0d than the χMAID-2013 and DMT-2001 models and that these observables provide excellent test object for different elementary operators.

4 Summary and outlook

The main topic of this article was to discuss theoretical uncertainties in the analyses of γdπ0d observables near threshold due to the use of elementary γNπN amplitudes using a model which is based on time-ordered perturbation theory. As input we have used the realistic MAID-2007 model [46] for the elementary amplitude and the high-precision Bonn full NN potential [55] for the deuteron wave function. We have presented results for the unpolarized differential and total cross sections as well as for all possible beam, target, and beam-target spin asymmetries of the differential and total cross sections in the photon energy region from near π-threshold to 170 MeV. In particular, we have studied the sensitivity of the calculated results to the choice of elementary pion photoproduction amplitude. For this purpose, we have used in addition to the MAID-2007 model [46] the realistic elementary amplitudes from DMT-2001 [47] and χMAID-2013 [48] models. We have also compared our calculations with the available experimental data.

We have found that the estimations of the uncertainty on the γdπ0d observables show important sensitivity to the modeling of the elementary γNπN operator. In many cases the deviation among results obtained using different elementary amplitudes is very large. It may be possible that the background contribution from the crossed nucleon pole amplitude plays a significant role in the results presented in this work, because the relative importance of direct and crossed nucleon pole amplitudes grows below the Δ-resonance region. The direct nucleon pole amplitude always interferes constructively with the Δ -resonance one, whereas the crossed nucleon pole amplitude interferes constructively for energies below the Δ -resonance region. The differences found in the predictions of γdπNN observables [34-39] using different elementary operators are also seen here in the results for γdπ0d observables. The calculated results for the unpolarized differential cross section using different elementary amplitudes are compared with the experimental data from TAPS [4]. A satisfactory agreement between theory and experiment is obtained only at forward pion angles. At backward pion angles, the calculations underestimate the experimental data and a disagreement is obtained.

The results presented above highlight the sensitivity of γdπ0d observables to the choice of pion photoproduction amplitude. The obtained differences in the numerical results for γdπ0d observables with various elementary operators can be attributed to the numerical account for the off-shell effects when the elementary nucleon amplitude is embedded into the deuteron process. In fact, it is difficult to give precise quantitative account for these off-shell effects, because the equivalence existing in the on-shell γNπN operator is broken in deuteron calculations when one nucleon or both of them are off their mass shells.

In summary, we conclude that the calculated results are of particular interest for the evaluation of the systematic uncertainties caused by the use of different elementary operators in the analyses of γdπ0d measurements. The calculated differential cross section does not described the TAPS data [4] at backward angles. This discrepancy between theory and experiment maybe resolved by considering the neglected effects from two-body mechanisms in the transition M-matrix. To study the influence of these reaction mechanisms on the γdπ0d observables, a more realistic treatment including all possible reaction mechanisms will be reported in a forthcoming paper [63]. An independent evaluation in the framework of effective field theory would also be very interesting. On the experimental side, precise measurements of γdπ0d observables can check these predictions and also provide a rigorous test of theoretical models.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through research groups program under grant number R.G.P.1/1/42. E.M. Darwish would like to thank M.I. Levchuk, V.V. Gauzshtein, and M.N. Nevmerzhitsky for many interesting and helpful discussions.

References

1. B. Krusche and S. Schadmand, Study of non-strange baryon resonances with meson photoproduction, Prog. Part. Nucl. Phys. 51 (2003) 399, https://doi.org/10.1016/S0146-6410(03)90005-6; V. D. Burkert and T.-S. H. Lee, Electromagnetic meson production in the nucleon resonance region, Int. J. Mod. Phys. E 13 (2004) 1035, https://doi.org/10.1142/S0218301304002545; K. Helbing, The Gerasimov-Drell-Hearn sum rule, Prog. Part. Nucl. Phys. 57 (2006) 405, https://doi.org/10.1016/j.ppnp.2005.09.003; A. Thomas, The Gerasimov-Drell-Hearn sum rule at MAMI, Eur. Phys. J. A 28 (2006) 161, https://doi.org/10.1140/epja/i2006-09-017-2; A. M. Sandorfi et al., Recent Polarization Experiments and the GDH Sum Rule, AIP Conf. Proc. 1155 (2009) 73, https://doi.org/10.1063/1.3203304; H. R. Weller et al., Research opportunities at the upgraded HI𝛾S facility, Prog. Part. Nucl. Phys. 62 (2009) 257, https://doi.org/10.1016/j.ppnp.2008.07.001. [ Links ]

2. B. Krusche, Photoproduction of mesons off nuclei, Eur. Phys. J. Spec. Topics 198 (2011) 199, https://doi.org/10.1140/epjst/e2011-01491-2; I. Jaegle (CBELSA/TAPS Collaboration), Quasi-free photoproduction of η-mesons off the deuteron, Eur. Phys. J. A 47 (2011) 89, https://doi.org/10.1140/epja/i2011-11089-0; D. K. Hasell et al., Spin- Dependent Electron Scattering from Polarized Protons and Deuterons with the BLAST Experiment at MIT-Bates, Ann. Rev. Nucl. Part. Sci. 61 (2011) 409, https://doi.org/10.1146/annurev-nucl-100809-131956; B. Strandberg et al., Threshold π production on the deuteron, EPJ Web Conf. 130 (2016) 05019, https://doi.org/10.1051/epjconf/201613005019. [ Links ]

3. B. Strandberg, Ph.D. thesis, University of Glasgow, 2017; S. E. Lukonin et al., Measurement of the tensor analyzing power for the reaction γd→pnπ0, Int. J. Mod. Phys. E 28 (2019) 1950010, https://doi.org/10.1142/S0218301319500101; B. Strandberg et al. (PIONS@Max-lab Collaboration), Near-threshold π—photoproduction on the deuteron, Phys. Rev. C 101 (2020) 035207, https://doi.org/10.1103/PhysRevC.101.035207. [ Links ]

4. U. Siodlaczek, Ph.D. thesis, Tübingen University, 2000; U. Siodlaczek et al., Coherent and incoherent π0 photoproduction from the deuteron, Eur. Phys. J. A 10 (2001) 365, https://doi.org/10.1007/s100500170120. [ Links ]

5. D.K. Hasell et al., Ann. Rev. Nucl. Part. Sci. 61 (2011) 409. [ Links ]

6. S. A. Zevakov et al., Measuring tensor analyzing power component T20 of the coherent photoproduction of a neutral pion on a tensor-polarized deuteron in the VEPP-3 storage ring, Bull. Russ. Acad. Sci. Phys. 79 (2015) 864, https://doi.org/10.3103/S1062873815070266; S. A. Zevakov et al., Neutral pion photoproduction on tensor-polarized deuterium on the VEPP-3 storage ring, Bull. Russ. Acad. Sci. Phys. 78 (2014) 611, https://doi.org/10.3103/S1062873814070260; D. M. Nikolenko et al., Measurement of the tensor analyzing power components in the coherent photoproduction of a π0 meson on a deuteron, JETP Lett. 89 (2009) 432, https://doi.org/10.1134/S0021364009090021; D. M. Nikolenko et al., Tensor observables in electro- and photoreactions on the deuteron, Phys. Part. Nucl. 48 (2017) 102, https://doi.org/10.1134/S1063779617010154; I. Rachek et al., Measurement of Tensor Analyzing Power T20 in Coherent π0 Photoproduction on Deuteron, Few Body Syst. 58 (2017) 29, https://doi.org/10.1007/s00601-016-1191-0; V. V. Gauzshtein et al., Measurements of the tensor analyzing power T20 of the reaction γd→dπ0, Int. J. Mod. Phys. E 29 (2020) 2050011, https://doi.org/10.1142/S0218301320500111; V. V. Gauzshtein et al., Measurement of the tensor analyzing power T20 for the reaction γd→dπ0, Eur. Phys. J. A 56 (2020) 169, https://doi.org/10.1140/epja/s10050-020-00175-z. [ Links ]

7. Y. Ilieva et al. (CLAS Collaboration), Coherent Pion Photoproduction on Deuterium, arXiv:nucl-ex/0309017; Y. Ilieva, Evidence for a backward peak in the γd→π0d cross section near the η threshold, Eur. Phys. J. A 43 (2010) 261, https://doi.org/10.1140/epja/i2010-10918-x. [ Links ]

8. W. J. Briscoe et al., (A2 Collaboration at MAMI), Cross section for γn→π0n measured at the Mainz A2 experiment, Phys. Rev. C 100 (2019) 065205, https://doi.org/10.1103/PhysRevC.100.065205; W. J. Briscoe, A. E. Kudryavtsev, I. I. Strakovsky, V. E. Tarasov, and R. L. Workman, Threshold π photoproduction on the neutron, Eur. Phys. J. A 56 (2020) 218, https://doi.org/10.1140/epja/s10050-020-00221-w; W. J. Briscoe, A. E. Kudryavtsev, I. I. Strakovsky, V. E. Tarasov, and R. L. Workman, On the photoproduction reactions γd→πNN, Eur. Phys. J. A 58 (2022) 23, https://doi.org/10.1140/epja/s10050-022-00671-4. [ Links ]

9. D.G. Ireland, E. Pasyuk, and I. Strakovsky, Prog. Part. Nucl. Phys. 111 (2020) 103752. [ Links ]

10. V.V. Gauzshtein, E.M. Darwish et al., Eur. Phys. J. A 56 (2020) 169. [ Links ]

11. V. V. Gauzshtein et al., Measurement of the tensor analyzing power T20 for the reaction γd→pnπ0 , Mod. Phys. Lett. A 36 (2021) 2150199, https://doi.org/10.1142/S0217732321501996. [ Links ]

12. O. Kolesnikov and A. Fix, Three-body calculation of coherent π0 photoproduction on the deuteron in the Δ region, J. Phys. G 48 (2021) 085101, https://doi.org/10.1088/1361-6471/ac010b. [ Links ]

13. T. Ishikawa et al (A2 Collaboration at MAMI)., Coherent photoproduction of the neutral-pion and η-meson on the deuteron at incident energies below 1.15 GeV, Phys. Rev. C 105 (2022) 045201, https://doi.org/10.1103/PhysRevC.105.045201. [ Links ]

14. F. Cividini et al. (A2 Collaboration at MAMI), Measurement of the helicity dependence for single π0 photoproduction from the deuteron, arXiv:2203.00535. [ Links ]

15. J. H. Koch and R. M. Woloshyn, Near threshold photoproduction of neutral pions from the deuteron, Phys. Rev. C 16 (1977) 1968, https://doi.org/10.1103/PhysRevC.16.1968. [ Links ]

16. P. Bosted and J. M. Laget, Electromagnetic properties of the πNN system: (II). The γD→Dπ0 reaction, Nucl. Phys. A 296 (1978) 413, https://doi.org/10.1016/0375-9474(78)90082-9; J. M. Laget, Electromagnetic properties of the πNN system: (I). The reaction γD→NNπ, Nucl. Phys. A 296 (1978) 388, https://doi.org/10.1016/0375-9474(78)90081-7; J. M. Laget, Pion photoproduction on few body systems, Phys. Rep. 69 (1981) 1, https://doi.org/10.1016/0370-1573(81)90164-2. [ Links ]

17. P. Wilhelm and H. Arenhövel, Dynamical Treatment of the NΔ- Interaction in Photo Reactions on the Deuteron, Few Body Syst. Suppl. 7 (1994) 235; P. Wilhelm and H. Arenhövel, Coherent pion photoproduction on the deuteron in the Δ resonance region, Nucl. Phys. A 593 (1995) 435, https://doi.org/10.1016/0375-9474(95)00348-5; P. Wilhelm and H. Arenhövel, Rescattering effects in coherent pion photoproduction on the deuteron in the Δ resonance region, Nucl. Phys. A 609 (1996) 469, https://doi.org/10.1016/S0375-9474(96)00277-1. [ Links ]

18. H. Garcilazo and E. Moya de Guerra, Pion photoproduction of the deuteron, The reaction γd→π0d , Phys. Rev. C 52 (1995) 49, https://doi.org/10.1103/PhysRevC.52.49; H. Garcilazo and E. Moya de Guerra, Reaction γd→π0d and the small components of the deuteron wave function, Phys. Rev. C 49 (1994) R601(R), https://doi.org/10.1103/PhysRevC.49.R601. [ Links ]

19. F. Blaazer, B. L. G. Bakker, and H. J. Boersma, Coherent pion production on the deuteron spin observables, Nucl. Phys. A 568 (1994) 681, https://doi.org/10.1016/0375-9474(94)90355-7; F. Blaazer, B. L. G. Bakker, and H. J. Boersma, Rescattering effects in coherent pion production on the deuteron, Nucl. Phys. A 590 (1995) 750, https://doi.org/10.1016/0375- 9474(95)00218-P; F. Blaazer, Ph.D. thesis, Free University of Amsterdam, 1995. [ Links ]

20. S. S. Kamalov, L. Tiator, and C. Bennhold, Coherent π0 and η photoproduction on the deuteron, Phys. Rev. C 55 (1997) 98, https://doi.org/10.1103/PhysRevC.55.98; S. S. Kamalov, L. Tiator, and C. Bennhold, Polarization observables in pion photoproduction on 3He, Nucl. Phys. A 547 (1992) 599, https://doi.org/10.1016/0375-9474(92)90653-2; S. S. Kamalov, L. Tiator, and C. Bennhold, Pion photoproduction on 3He including final-state interaction, Few Body Syst. 10 (1991) 143, https://doi.org/10.1007/BF01076291. [ Links ]

21. E. M. Darwish and S. S. Al-Thoyaib, Invariant amplitudes for electromagnetic pion production from the nucleon and its implication for separated structure functions, Arab. J. Sci. Eng. 33 (2008) 401; E. M. Darwish and S. S. Al-Thoyaib, Talk presented at the 7th NUPPAC’09 Conference on Nuclear and Particle Physics, November 11-15, 2009, Egypt. [ Links ]

22. E. M. Darwish and M. Y. Hussein, Contribution of coherent pion photoproduction to the Gerasimov-Drell-Hearn sum rule for the deuteron, Appl. Math. Inf. Sci. 3 (2009) 321. [ Links ]

23. E. M. Darwish, N. Akopov, and M. El-Zohry, Coherent π0- photoproduction on the deuteron including polarization observables, AIP Conf. Proc. 1370 (2011) 242, https://doi.org/10.1063/1.3638108; E. M. Darwish, N. Akopov, and M. El-Zohry, Coherent π0-Photoproduction on the Deuteron, Proceeding of the 35th International Conference of High Energy Physics, July 22- 28, Paris, https://doi.org/10.22323/1.120.0185. [ Links ]

24. E. M. Darwish and S. S. Al-Thoyaib, Coherent π0-photoproduction on the deuteron near the η-production threshold including polarization observables, Ann. Phys. 351 (2014) 35, https://doi.org/10.1016/j.aop.2014.08.008; E. M. Darwish and A. Hemmdan, Influence of intermediate ηNN interaction on spin asymmetries for γd→π0d reaction near the η-threshold within a three-body approach, Ann. Phys. 356 (2015) 128, https://doi.org/10.1016/j.aop.2015.02.035. [ Links ]

25. E. M. Darwish and H. Mansour, Coherent π-Production off Deuteron Near η-Threshold (Lambert Academic Publishing, Saarbrücken, 2015). [ Links ]

26. E. M. Darwish, Influence of the Choice NN Potential Model on γd→π0d Observables near η-Threshold, Chin. Phys. Lett. 33 (2016) 041301, https://doi.org/10.1088/0256-307X/33/4/041301. [ Links ]

27. E. M. Darwish, H. M. Abou-Elsebaa, and K. S. A. Hassaneen, Tensor Target Spin Asymmetries in Coherent π0-Photoproduction on the Deuteron Including Intermediate ηNN Interaction Within a Three-Body Approach, Braz. J. Phys. 48 (2018) 168, https://doi.org/10.1007/s13538-018-0559-7; E. M. Darwish, E. M. Mahrous, and M. E. Alsehli, Coherent π0-photoproduction from the deuteron near threshold, AIP Conf. Proc. 1976 (2018) 020035, https://doi.org/10.1063/1.5042402. [ Links ]

28. E. M. Darwish, Coherent Photoproduction of π0-meson from the Deuteron Including Polarization Effects, Q. Phys. Rev. 4 (2018) 1729. [ Links ]

29. E. M. Darwish and M. Saleh Yousef, Coherent π0-Photoproduction on the Deuteron near Threshold and the Role of D-Wave Component of the Deuteron Wave Function, Moscow Univ. Phys. Bull. 74 (2019) 595, https://doi.org/10.3103/S0027134919060122. [ Links ]

30. E. M. Darwish, H. M. Abou-Elsebaa, Kh. S. Alsadi, and M. Saleh Yousef, The Spin Response of the γd→π0d Reaction Near Threshold and Its Implication to the GDH Sum Rule and the Double Polarization E-Asymmetry, Moscow Univ. Phys. Bull. 75 (2020) 198, https://doi.org/10.3103/S002713492003008X. [ Links ]

31. H. M. Al-Ghamdi, E. S. Almogai, E. M. Darwish, and S. Abdel-Khalek, Beam-target double spin asymmetries in the reaction γd→π0d near threshold and the role of D-wave component of the deuteron wave function, Braz. J. Phys. 50 (2020) 615, https://doi.org/10.1007/s13538-020- 00771-w. [ Links ]

32. H. Al-Ghamdi, E. S. Almogait, E. M. Darwish, and S. Abdel-Khalek, Tensor analyzing power component T20 of the γd→π0d process in the photon energy range from 200 to 400 MeV, Res. Phys. 18 (2020) 103238, https://doi.org/10.1016/j.rinp.2020.103238. [ Links ]

33. E. M. Darwish and H. M. Al-Ghamdi, Sensitivity of Polarization Observables in γd→π0d Reaction Near Threshold to the Choice of Elementary γN→πN Amplitude and Deuteron Wave Function, Moscow Univ. Phys. Bull. 76 (2021) 136, https://doi.org/10.3103/S0027134921030036. [ Links ]

34. E. M. Darwish, C. Fernández-Ramírez, E. Moya de Guerra, and J. M. Udías, Helicity dependence and contribution to the Gerasimov-Drell-Hearn sum rule of the γd→πNN reaction channels in the energy region from threshold up to the Δ(1232) resonance, Phys. Rev. C 76 (2007) 044005, https://doi.org/10.1103/PhysRevC.76.044005. [ Links ]

35. M. I. Levchuk, Helicity-dependent reaction γd→πNN and its contribution to the Gerasimov- Drell-Hearn sum rule for the deuteron, Phys. Rev. C 82 (2010) 044002, https://doi.org/10.1103/PhysRevC.82.044002. [ Links ]

36. E. M. Darwish and S. S. Al-Thoyaib, Incoherent pion photoproduction on the deuteron including polarization effects, Ann. Phys. 326 (2011) 604, https://doi.org/10.1016/j.aop.2010.09.014. [ Links ]

37. E. M. Darwish, M. M. Almarashi, E. M. Mahrous, M. A. Hassanain, and M. Saleh Yousef, Near-threshold incoherent pion photoproduction on the deuteron with final-state interaction effects, Ann. Phys. 411 (2019) 167990, https://doi.org/10.1016/j.aop.2019.167990. [ Links ]

38. E. M. Darwish, M. M. Almarashi, and M. Saleh Yousef, Beam-targeet double spin asymmetries in the reaction γd→πNN near threshold with final-state rescattering effects, Ann. Phys. 420 (2020) 168254, https://doi.org/10.1016/j.aop.2020.168254. [ Links ]

39. E. M. Darwish, M. I. Levchuk, M. N. Nevmerzhitsky, M. M. Almarashi, and M. Saleh Yousef, Influence of rescattering effects on the helicity dependence of γd→πNN near threshold and its implication to the E-asymmetry and the GDH sum rule, Chin. J. Phys. 71 (2021) 319, https://doi.org/10.1016/j.cjph.2020.07.015. [ Links ]

40. R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, and R. L. Workman, Analysis of pion photoproduction data, Phys. Rev. C 66 (2002) 055213, https://doi.org/10.1103/PhysRevC.66.055213. [ Links ]

41. D. Drechsel, O. Hanstein, S. S. Kamalov, and L. Tiator, A unitary isobar model for pion photo- and electroproduction on the proton up to 1 GeV, Nucl. Phys. A 645 (1999) 145, https://doi.org/10.1016/S0375-9474(98)00572-7. [ Links ]

42. A. V. Anisovich et al., Photoproduction of baryons decaying into and , Eur. Phys. J. A 25 (2005) 427, https://doi.org/10.1140/epja/i2005- 10120-5; A. V. Anisovich, E. Klempt, A. V. Sarantsev, and U. Thoma, Partialwave decomposition of pion and photoproduction amplitudes, Eur. Phys. J. A 24 (2005) 111, https://doi.org/10.1140/epja/i2004-10125-6; A. V. Anisovich et al., Properties of baryon resonances from a multichannel partial wave analysis, Eur. Phys. J. A 48 (2012) 15, https://doi.org/10.1140/epja/i2012-12015-8. [ Links ]

43. T. Sato and T.-S. H. Lee, Meson-exchange model for πN scattering and γN→πN reaction, Phys. Rev. C 54 (1996) 2660, https://doi.org/10.1103/PhysRevC.54.2660; T. Sato and T.-S. H. Lee, Dynamical study of the Δ excitation in N(e,e'π) reactions, Phys. Rev. C 63 (2001) 055201, https://doi.org/10.1103/PhysRevC.63.055201; H. Kamano, S. X. Nakamura, T.-S. H. Lee, and T. Sato, Nucleon resonances within a dynamical coupled-channels model of πN and γN reactions, Phys. Rev. C 88 (2013) 035209, https://doi.org/10.1103/PhysRevC.88.035209; H. Kamano, T.-S. H. Lee, S. X. Nakamura, and T. Sato, The ANL-Osaka Partial-Wave Amplitudes of and πN and γN reactions, arXiv:1909.11935. [ Links ]

44. T. Feuster and U. Mosel, Unitary model for meson-nucleon scattering, Phys. Rev. C 58 (1998) 457, https://doi.org/10.1103/PhysRevC.58.457; V. Shklyar, H. Lenske, and U. Mosel, Analysis of baryon resonances in the coupled-channel Giessen model, Proceedings of the 11th Workshop on The Physics of Excited Nucleons, September 5-8, Bonn, https://doi.org/10.1007/978-3-540-85144-8_16; V. Shklyar, H. Lenske, and U. Mosel, η photoproduction in the resonance energy region, Phys. Lett. B 650 (2007) 172, https://doi.org/10.1016/j.physletb.2007.05.005; V. Shklyar, H. Lenske, and U. Mosel, Giessen coupled-channel model for pion- and photon-induced reactions, AIP Conf. Proc. 1432 (2012) 349, https://doi.org/10.1063/1.3701246. [ Links ]

45. C. Fernández-Ramírez, E. Moya de Guerra, and J. M. Udías, Crossing symmetry and phenomenological widths in effective Lagrangian models of the pion photoproduction process, Phys. Lett. B 660 (2008) 188, https://doi.org/10.1016/j.physletb.2007.11.099; C. Fernández-Ramírez, Ph.D. tesis, Universidad Complutense de Madrid, 2006. [ Links ]

46. D. Drechsel, S. S. Kamalov, and L. Tiator, Unitary isobar model- MAID2007, Eur. Phys. J. A 34 (2007) 69, https://doi.org/10.1140/epja/i2007-10490-6. [ Links ]

47. S. S. Kamalov and S. Nan Yang, Pion Cloud and the Q2 Dependence of γ*N↔Δ Transition Form Factors, Phys. Rev. Lett. 83 (1999) 4494, https://doi.org/10.1103/PhysRevLett.83.4494; S. S. Kamalov, S. Nan Yang, D. Drechsel, O. Hanstein, and L. Tiator, γ*N→Δ transition form factors: A new analysis of data on p(e,e'p)π0 at Q2 = 208 and 4.0 (GeV/c)2, Phys. Rev. C 64 (2001) 032201(R), https://doi.org/10.1103/PhysRevC.64.032201. [ Links ]

48. M. Hilt, S. Scherer, and L. Tiator, Threshold π0 photoproduction in relativistic chiral perturbation theory, Phys. Rev. C 87 (2013) 045204, https://doi.org/10.1103/PhysRevC.87.045204; M. Hilt, B. C. Lehnhart, S. Scherer, and L. Tiator, Pion photo-and electroproduction in relativistic baryon chiral perturbation theory and the chiral MAID interface, Phys. Rev. C 88 (2013) 055027, https://doi.org/10.1103/PhysRevC.88.055207. [ Links ]

49. S. B. Gerasimov, A sum rule for magnetic moments and the damping of the nucleon magnetic moment in nuclei, Yad. Fiz. 2 (1965) 598 [Sov. J. Nucl. Phys. 2 (1966) 430]; S. D. Drell and A. C. Hearn, Exact Sum Rule for Nucleon Magnetic Moments, Phys. Rev. Lett. 16 (1966) 908, https://doi.org/10.1103/PhysRevLett.16.908. [ Links ]

50. H. Arenhövel, General formulae for polarization observables in two-body break-up of deuteron photodisintegration, Few Body Syst. 4 (1988) 55, https://doi.org/10.1007/BF01076329. [ Links ]

51. J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics (McGraw- Hill, New York, 1964). [ Links ]

52. H. Arenhövel, Invariant Amplitudes for Coherent Electromagnetic Pseudoscalar Production from a Spin-One Target, II: Crossing, Multipoles, and Observables, Few Body Syst. 27 (1999) 141, https://doi.org/10.1007/s006010050127. [ Links ]

53. A. R. Edmonds, Angular Momentum in Quantum Mechanics (Princeton University Press, New Jersey, 1957). [ Links ]

54. H. Arenhövel, Spin degrees and polarization observables in electromagnetic reactions, Int. J. Mod. Phys. E 18 (2009) 1226, https://doi.org/10.1142/S0218301309013476. [ Links ]

55. R. Machleidt, K. Holinde, and Ch. Elster, The Bonn meson-exchange model for the nucleon-nucleon interaction, Phys. Rep. 149 (1987) 1, https://doi.org/10.1016/S0370-1573(87)80002-9. [ Links ]

56. R. Machleidt, High-precision, charge-dependent Bonn nucleon-nucleon potential, Phys. Rev. C 63 (2001) 024001, https://doi.org/10.1103/PhysRevC.63.024001. [ Links ]

57. M. Lacombe et al., Parametrization of the deuteron wave function of the Paris N-N potential, Phys. Lett. B 101 (1981) 139, https://doi.org/10.1016/0370-2693(81)90659-6. [ Links ]

58. B. D. Day, Nuclear Saturation from Two-Nucleon Potentials, Phys. Rev. Lett. 47 (1981) 226, https://doi.org/10.1103/PhysRevLett.47.226. [ Links ]

59. M. Hjorth-Jensen, T. T. S. Kuo, and E. Osnes, Realistic effective interactions for nuclear systems, Phys. Rep. 261 (1995) 125, https://doi.org/10.1016/0370- 1573(95)00012-6. [ Links ]

60. E. M. Darwish, H. M. Abou-Elsebaa, E. M. Mahrous, and S. S. Al-Thoyaib, Sensitivity of tensor and vector analyzing powers in elastic ed scattering to modern local and nonlocal NN potentials, Ind. J. Phys. 94 (2020) 1025, https://doi.org/10.1007/s12648-019-01537-0. [ Links ]

61. M. Benmerrouche and E. Tomusiak, γD→π0D reaction in the threshold region, Phys. Rev. C 58 (1998) 1777, https://doi.org/10.1103/PhysRevC.58.1777. [ Links ]

62. J. A. Gómez Tejedor, S. S. Kamalov, and E. Oset, Meson Exchange currents in the 3He(γ,π+)3H reaction, Phys. Rev. C 54 (1996) 3160, https://doi.org/10.1103/PhysRevC.54.3160. [ Links ]

63. M.I. Levchuk, V.V. Gauzshtein, E.M. Darwish, M.N. Nevmerzhitsky, (in preparation). [ Links ]

Received: March 16, 2022; Accepted: March 31, 2022

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License