-
The interaction between nucleons is very important in understanding the nuclear experiments and the properties of nuclei and neutron stars. However, the interaction between nucleons is still poorly understood, particularly at high densities. The main reason for this is that the interaction is not a directly observable quantity in experiments. The interaction between nucleons as residual interactions of the QCD cannot be solved via perturbation expansion. Therefore, the interaction between nucleons leads to the uncertainty of the symmetry energy. Considerable effort has been directed toward understanding the interaction between nucleons. The global Dirac phenomenological optical potentials, which describe the elastic proton scattering observables of
12 C and208 Pb well for incident energies between 20 and 1040 MeV [1, 2], are powerful tools in understanding the interaction between nucleons. They have been widely used to extract the symmetry energy directly or indirectly [3, 4]. The relativistic impulse approximation (RIA) [5–7] can also be used to describe the elastic proton scattering observables. Theoretically, the symmetry energy can be extracted from the microscopic potentials of the RIA.The original RIA with scalar and vector potentials can reproduce the analyzing power and spin-rotation parameter in proton-nucleus scattering for incident energies above 500 MeV [5, 8]. Because the direct and exchange parts of the t-matrix are not separated explicitly in the original RIA, the optical potential of the original RIA diverges at energies below 500 MeV [8]. By separating the direct and exchange parts of the t-matrix, it is found that the pion with pseudoscalar leads to the divergent optical potential. This can be cured by introducing the pion with pseudovector coupling, in a process known as the relativistic Love-Franey (RLF) model [6, 9,10]. The RLF model has been used to describe the elastic proton scattering observables successfully for incident energies of 135, 200, 300, and 400 MeV [6]. The shortcoming of the original RLF model is that the scattering observables were fitted separately at each energy. The separate fits of the scattering observables cannot reflect the energy-dependent trend of the scattering observables. This motivates the proposal of the energy-dependent RLF model [11, 12]. The energy-dependent RLF model has been used to fit empirical amplitudes successfully at incident energies of 50–200 MeV [13], 200–500 MeV [11], and 500–800 MeV [12]. It should be noted that these fits are separated from each other. In this work, we try to construct the energy-dependent RLF model over the laboratory energy range of 20 to 800 MeV within a unified fit.
The RLF model has been used to study both elastic proton-nucleus scattering [6, 7] and the (p, 2p) reaction [14]. The RLF model without Pauli blocking (PB) corrections cannot well describe the differential cross sections
dΩ/dσ , analyzing powersAy , and spinrotation functions Q of elastic proton-nucleus scattering. With the PB corrections, the RLF model can describedΩ/dσ ,Ay , and Q well [6, 7]. The RLF model is also helpful in understanding the in-medium effect. In-medium nucleon-nucleon (NN) cross sections, which are key inputs of various transport models, can be obtained from RLF. It has been recognized that the theoretical uncertainty in extracting the density dependence of the symmetry energy from data of heavy-ion reactions can be largely ascribed to the poor knowledge and premature treatment of the isospin-dependent in-medium NN cross sections [15, 16]. Although the NN cross section in free space is well-determined, the in-medium NN cross section is still model dependent [17–28]. Under the original RIA framework, the in-medium NN cross sections have been studied [24, 28]. However, no such work has been done in the framework of the RLF model. At a low nucleon kinetic energy, it would be interesting to study the in-medium NN cross sections with the RLF model. Additionally, the energy and density dependent symmetry potential of nuclear matter has been investigated with the original RIA and RLF models [24, 28, 29]. The energy and density dependence symmetry potential of the energy-dependent RLF model may have some differences from the results of the original RIA and RLF models.The remainder of this paper is organized as follows. In Sec. II, we briefly introduce formulas and approaches used in this work. The formulas include the fits between RIA scattering amplitudes and the empirical amplitudes, as well as the elastic proton scattering observables of
208 Pb. Results and discussions are presented in Sec. III. Finally, a summary is given in Sec. IV. -
In order to obtain the on-shell NN matrix elements of RLF scattering amplitudes, one needs to transform the NN matrix elements of RLF scattering amplitudes into those of the non-relativistic NN scattering amplitudes. The relation between the matrix elements of the RLF scattering amplitudes and the matrix elements of the non-relativistic NN scattering amplitudes in the center-of-mass frame can be written as [5, 6, 9, 10, 13]
(2ikc)−1χ†s′1χ†s′2ˆfcmχs1χs2=ˉu(k′1,s′1)ˉu(k′2,s′2)ˆFu(k1,s1)u(k2,s2),
(1) where
k1≡(Ec,kc) andk2≡(Ec,−kc) are the incident 4-momenta in the center-of-mass frame.χ′ s represent the usual Pauli spinors.kc represents the momentum in the center of mass frame. The Dirac spinors are normalized (ˉuu =1).ˆfcm andˆF are the non-relativistic scattering matrix and scattering matrix of RLF respectively.ˆfcm contains five real and imaginary amplitudes [30]:(2ikc)−1ˆfcm=12[a+b+(a−b)σ1nσ2n+(c+d)σ1mσ2m+(c−d)σ1lσ2l+e(σ1n+σ2n)],
(2) where a, b, c, d, and e are the non-relativistic real and imaginary amplitudes. The scattering matrix of RLF
ˆF also contains five real and imaginary amplitudes:ˆF=T∑i=SFi(s,t,u)λ(1)iλ(2)i,
(3) where i includes the scalar (S), vector (V), pseudoscalar (P), axial vector (A), and tensor (T).
FS ,FV ,FP ,FA , andFT are the relativistic amplitudes of RLF. The Dirac matrices of S, V, P, A, and T areI4×4 ,γμ ,γ5 ,γ5γμ , andσμν respectively. With these Dirac matrices and Eq. (1), the relation between the non-relativistic amplitudes and the relativistic amplitudes is obtained as follows [5, 9, 10, 13]:(abcde)=ikcO5×5(FSFVFPFAFT),
(4) The
O5×5 matrix is written as [10, 13]O5×5=(cosθcosθcosθ−cosθ−4sinθ1−1110−1111011−110−isinθ−isinθ−isinθisinθ−4icosθ)×(αhα+pβ0pα+hβ−2β−hγ−γ−pγ−δ2(p+h)δα(p+h)α0(h−p)α2αhγγ−pγ−γ2γ−ϵ−ϵ0ϵ−2ϵ),
(5) with
α=cos2θ2,β=1+sin2θ2,γ=sin2θ2,δ=1+cos2θ2,ϵ=EcMsinθ2cosθ2,Ec=(k2c+M2)1/2,h=E2cM2,p=k2cM2,
(6) where θ represents the scattering angle in the two-body center-of-mass, and M represents the free nucleon mass. Now, the relation between the matrix elements of the RLF scattering amplitudes and the matrix elements of the non-relativistic NN scattering amplitudes is clear. For comparison with the empirical amplitudes, the relativistic amplitudes
Fi(s,t,u) need to be solved. TheFi(s,t,u) contains direct and exchange contributions:Fi(s,t,u)=iM22Eckc[FDi(s,t)+FEi(s,u)],
(7) with
FDi(s,t)=N∑j=1δi,i(j)⟨τ1⋅τ2⟩Tjfj(Ec,|q|),FEi(s,u)=(−1)TNNN∑j=1Ci(j),i⟨τ1⋅τ2⟩Tjfj(Ec,|Q|),
(8) where N represents the number of mesons in this work,
Ci(j),i is the Fierz matrix [6, 9, 31],Tj is the isospin of the jth meson, andTNN is the total isospin of the two-nucleon state. WhenTNN=1 ,⟨τ1⋅τ2⟩Tj=1 for bothTj=0 andTj=1 . WhenTNN=0 ,⟨τ1⋅τ2⟩Tj=1 forTj=0 , and⟨τ1⋅τ2⟩Tj=−3 forTj=1 . For instance, whenTNN=0 , the scalar amplitudes are written asFDS(s,t)=fS(Tj=0)(Ec,|q|)−3fS(Tj=1)(Ec,|q|),FES(s,u)=∑j=S,TCj,S[fj(Tj=0)(Ec,|Q|)−3fj(Tj=1)(Ec,|Q|)],
(9) where
fj(Ec,|q|) andfj(Ec,|Q|) contain the real and imaginary parts. For instance:fj(Ec,|q|)=fjR(Ec,|q|)−ifjI(Ec,|q|),fjR(Ec,|q|)=g2j(Ec)q2+m2j(1+q2/Λ2j)−2,fjI(Ec,|q|)=ˉg2j(Ec)q2+ˉm2j(1+q2/ˉΛ2j)−2,
(10) where
g2j(Ec) ,mj , andΛj are the coupling constant, mass, and momentum cutoff of the jth real meson, respectively.ˉg2j(Ec) ,ˉmj , andˉΛj are the coupling constant, mass, and cutoff parameter of the jth imaginary meson, respectively. The direct momentum transfer|q| and the exchange momentum transfer|Q| are written as [9]|q|=|kc−k′c|=2kcsinθ2,|Q|=|kc+k′c|=2kcsin(π−θ2).
(11) For a given projectile's kinetic energy
Tlab , the momentum and energy of the center of mass frame arekc=√2MTlab2 andEc=√k2c+M2 , respectively. With these fitting formulas, the parameterization of the RLF model can be obtained. Then, the Dirac optical potential is generated with the fitting parameters of the RLF model.The Dirac optical potential of the RLF model is the same as that of the original RIA model [6, 32, 33]:
Uopt(q,E)=−4πiplabM⟨ˉΨ|A∑n=1ˆF(q,E;n)eiq⋅r(n)|Ψ⟩,
(12) where
plab represents the laboratory momentum of incident nucleons. q and E represent the momentum transfer and collision energy, respectively.|Ψ⟩ represents the A-particle ground state of the target nucleus.In proton-nucleus scattering, the scattering process can be approximately treated as incident protons scattered by each of the nucleons in the target nucleus. Assuming the incident proton wave function as
U0(r) and the relativistic Hartree wave function of the nucleus ground state asU0(r) , the optical potential can act on the incident proton wave function and project into the coordinate space. In this work, we only list the formulas after Fourier transforms [33]:⟨r|Uopt(q,E)|U0(r)⟩=−4iπplabM∑L[∫d3r′ρL(r′)tLD(|r′−r|;E)]λLU0(r)−4iπplabM∑L[∫d3r′ρL(r′,r)tLE(|r′−r|;E)]λLU0(r′),
(13) where
tLD(|r|;E)]≡∫d3q(2π)3tLD(q,E)e−iq⋅r,
(14) tLE(|r|;E)]≡∫d3q(2π)3tLE(Q,E)e−iQ⋅r,
(15) with
tLD(q,E)≡(iM2/2Eckc)FLD(q) andtLE(Q,E)≡(iM2/2Eckc)FLE(Q) , respectively. The first term of the above equation represents the direct optical potential in the coordinate space:ULD(r,E)=−4iπplabM∫d3r′ρL(r′)tLD(|r′−r|;E),
(16) and the second term of the above equation represents the exchange optical potential. With the local-density approximation, the exchange optical potential is written as
ULE(r,E)=−4iπplabM∫d3r′ρL(r′,r)×tLE(|r′−r|;E)j0(plab|r′−r|),
(17) j0 is a spherical Bessel function. The nuclear densities of Eqs. (13), (16), and (17) are defined asρL(r′,r)≡occ′∑αˉUα(r′)λLUα(r),ρL(r)≡ρL(r,r).
(18) The occupied states indicate that one sums over target protons for the density with pp amplitudes and over target neutrons for the density with pn amplitudes. The densities of target protons and neutrons can be obtained from relativistic mean field theory (RMF) [7, 33]. When the densities of target protons and neutrons predicted by various models are substituted into Eq. (18), the results are almost the same in p+
208 Pb elastic scattering. In the following sections, we take the densities predicted by the code of Ref. [33] as the input.With above formulas, the Dirac optical potential can be obtained. For the spin-saturated nucleus, the Dirac optical potential only contains scalar, vector, and tensor terms. Because the tensor contribution is small, the tensor contribution has been neglected in the calculation. The Dirac optical potential is simplified as
Uopt(r;E)=US(r;E)+γ0UV(r;E),UL(r;E)=ULD(r;E)+ULE(r;E).
(19) For the laboratory energy near 200 MeV, the contribution of Pauli blocking (PB) plays an important role in the calculations. The modification of PB with a local-density approximation is [6, 13, 33]
ULPB(r;E)=[1−aL(E)(ρB(r)ρ0)2/3]UL(r;E).
(20) Here,
ρB(r) represents the local baryon density of the target nuclus, andρ0=0.1934 fm−3 .With the Dirac optical potential, the Dirac equation of the projectile is written as
hU0(r′)=EU0(r′)={−iα⋅∇+UV(r;E)+β[M+US(r;E)]}U0(r′).
(21) By solving the Dirac equation, i.e., Eq. (21), the scattering observables (the differential cross section
dσ/dΩ , analyzing powerAy , and spin rotation function Q) are determined. The detailed solution processes are presented in Ref. [33]. In this work, we concentrate ondσ/dΩ ,Ay , and Q of p+208 Pb elastic scattering. -
The RLF pp and pn scattering amplitudes are written as [9]
Fi(pp)=Fi(TNN=1),Fi(pn)=12[Fi(TNN=1)+Fi(TNN=0)],
(22) where the relativistic amplitudes
FS ,FV ,FP , andFA of pp and pn and the non-relativistic amplitudes a, b, c, d, and e of pp and pn are generated from Eq. (4). However, the non-relativistic amplitudes are no longer available on the on-line Scattering Analysis Interactive Dial-in (SAID) facility. The scattering observables are still available on SAID. There are several models available for computing the scattering observables on SAID. We choose the scattering observables obtained from the weighted fits (WF16) on SAID [34]. In contrast to previous studies [11–13], we fit 10 scattering observables of the weighted fits (WF16) rather than the real and imaginary amplitudes of pp and pn. These 10 scattering observables are written as [30]dσdΩ≡σ=12(|a|2+|b|2+|c|2+|d|2+|e|2),D=12(|a|2+|b|2−|c|2−|d|2+|e|2)/σ,DT=12(|a|2−|b|2+|c|2−|d|2+|e|2)/σ,Ayy=12(|a|2−|b|2−|c|2+|d|2+|e|2)/σ,
P=ℜ(a∗e)/σ,CKP=ℑ(d∗e)/σ,CKK=ℜ(a∗d+b∗c)/σ,CPP=−ℜ(a∗d−b∗c)/σ,R=[ℜ(a∗b)cos(α+θ2)+ℜ(c∗d)cos(α−θ2)−ℑ(b∗e)sin(α+θ2)],A=[−ℜ(a∗b)sin(α+θ2)+ℜ(c∗d)sin(α−θ2)−ℑ(b∗e)cos(α+θ2)],
(23) These 10 scattering observables are the unpolarized differential cross section (
dσdΩ ), depolarization tensor for polarized beam (D), polarization transfer from beam to recoil particle (DT), polarization correlation for initially unpolarized particlesAyy , CKP, CKK, CPP, polarization of scattered particle (P), and triple scattering parameters A and R. α is a function of the angle:α=θ2−θlab,
(24) where
θlab represents the laboratory scattering angle. Theoretically, the fitting quality with these 10 scattering observables of pp and pn should be as good as that in the previous case with the non-relativistic amplitudes of pp and pn.The real isovector pseudoscalar (π), isoscalar vector (ω), and isovector vector (ρ) masses are fixed to their physical values
mπ=138 MeV,mω=782 MeV, andmρ= 770 MeV, respectively. The remaining meson masses are set in accordance with Ref. [13]. Because the energy boundary between low energy and high energy is around 200 MeV, we first fit the scattering observables to obtain the initial coupling strengths of RLF. The initial coupling constantg20 (ˉg20 ) and the cutoff parameter Λ (ˉΛ ) are obtained by fitting the scattering observables forTlab= 200 MeV alone. In order to obtain the best fit to the scattering observables forTlab= 200 MeV, the Levenberg-Marquadt method has been employed to minimize the value ofχ2 .χ2 is defined as follows:χ2=∑data(χempirical−χfit)2⟨χ2empirical⟩,
(25) where
⟨χ2empirical⟩ represents an angle-averaged value and is defined as⟨χ2empirical⟩=175∘∑θ=5∘χ2empirical(θ)Nempirical,
(26) where
Nempirical=35 represents the number of angles per energy. For a single incident energy, the number of the data points is 700: 35 angles× 10 scattering observables× 2 pp and pn. As a result, the real (imaginary) initial coupling constantg20 (ˉg20 ), the real (imaginary) meson massm (ˉm ), and the real (imaginary) cutoff parameter Λ (ˉΛ ) are presented in Table 1. It should be noted thatΛ>m has been employed for a given meson.Real parameters Meson Isospin Coupling type m g20 a1 a2 a3 a4 Λ σ 0 Scalar (S) 600 −8.348 6.527 ×10−1 2.773 −2.741 ×10−2 −7.639 962.7 ω 0 Vector (V) 782 1.019 ×101 −6.246 ×10−1 −1.400 ×101 8.492 ×10−1 2.677 ×101 1155 t0 0 Tensor (T) 550 2.794 ×10−1 3.303 1.172 ×101 −6.676 ×10−1 −3.008 ×101 1945 a0 0 Axial vector (A) 500 4.921 ×10−1 6.066 1.953 ×101 −1.089 −5.133 ×101 1597 η 0 Pseudoscalar (P) 950 1.180 ×101 −1.694 ×102 −1.350 ×102 5.417 5.550 ×102 957.3 δ 1 Scalar (S) 500 6.095 ×10−3 5.635 ×10−1 1.287 −6.040 ×10−2 −4.101 2945 ρ 1 Vector (V) 770 −1.667 ×10−1 9.420 ×10−5 1.644 ×10−1 −2.905 ×10−2 2.915 ×10−1 3001 t1 1 Tensor (T) 600 −2.500 ×10−1 −1.893 −7.848 4.448 ×10−1 1.989 ×101 1280 a1 1 Axial vector (A) 650 −1.356 −4.678 −1.529 ×101 8.102 ×10−1 4.146 ×101 745.1 π 1 Pseudoscalar (P) 138 1.095 ×101 −3.834 ×10−9 −4.889 4.005 ×10−1 1.082 ×101 677.4 Imaginary parameters Meson Isospin Coupling type ˉm ˉg20 a1 a2 a3 a4 ˉΛ σ 0 Scalar (S) 600 −2.849 −4.637 −1.475 ×101 6.301 ×10−1 4.377 ×101 800.0 ω 0 Vector (V) 700 4.405 4.825 2.077 ×101 −9.616 ×10−1 −5.672 ×101 716.0 t0 0 Tensor (T) 750 −9.457 ×10−1 −1.474 −7.535 4.140 ×10−1 1.914 ×101 1340 a0 0 Axial vector (A) 750 −2.198 −3.583 −1.374 ×101 6.840 ×10−1 3.762 ×101 1437 η 0 Pseudoscalar (P) 1000 1.420 ×101 −1.619 ×102 −2.011 ×102 9.779 6.591 ×102 2758 δ 1 Scalar (S) 650 3.094 1.593 2.945 1.614 ×10−2 −1.328 ×101 751.9 ρ 1 Vector (V) 600 −2.421 −9.543 ×10−1 −3.644 9.074 ×10−2 1.292 ×101 766.6 t1 1 Tensor (T) 750 6.140 ×10−1 1.513 5.775 −3.011 ×10−1 −1.557 ×101 1662 a1 1 Axial vector (A) 1000 2.265 5.071 1.997 ×101 −1.011 −5.450 ×101 1233 π 1 Pseudoscalar (P) 500 −5.891 2.029 ×101 2.363 ×101 −9.593 ×10−1 −7.574 ×101 1342 Table 1. Parameter sets for the energy-dependent RLF over the laboratory energy range of 20 to 800 MeV. The masses m (
ˉm ) and cutoff parameters Λ (ˉΛ ) are in MeV, whereas the remaining parameters are dimensionless.The energy-dependent RLF model has been used to fit the empirical NN data for separate incident energies (50–200 MeV [13], 200–500 MeV [11], and 500–800 MeV [12]). To relate these energy domains, we first fit the scattering observables of various energies (
Tlab= 20, 50, 75, 100, 150, 200, 300, 400, 500, 600, and 800 MeV) separately with the same cutoff parameters and meson masses as 200 MeV. As shown in Fig. 1, we find that the coupling strengths, such asg2π ,g2η ,g2a1 , andg2t1 , approximately follow a logarithmic or parabolic curve opening to the right. Combining the linear and square terms [11, 12], the energy-dependent coupling strengths can be written asFigure 1. (color online) The coupling constants (
g2π ,g2η ,g2a1 , andg2t1 ) as a function of the laboratory energiesTlab . The coupling constants are obtained from separate fits withTlab= 20, 50, 75, 100, 150, 200, 300, 400, 500, 600, and 800 MeV.g2(Ec)=g20+a1ln(TlabT0)+a2(TlabT0−1)+a3[(TlabT0)2−1]+a4(√TlabT0−1),
(27) with
T0=200 MeV. This equation is also valid for the imaginary coupling strengths. When the incident energy is 200 MeV, the energy-dependent coupling strengths should equal the initial coupling strengths. The meson masses and cutoff parameters of RLF are assumed to be independent of the energy in this work.For various incident energies, the initial coupling constant, the meson mass, and the cutoff parameter remain the same. However, the coupling constants
g2j(Ec) andˉg2j(Ec) change with respect to the incident energy. The parametersa1 ,a2a3 , anda4 of Eq. (27) are obtained by fitting the scattering observables for various incident energies (Tlab= 20, 50, 75, 100, 150, 200, 300, 400, 500, 600, and 800 MeV) uniformly. The number of the data points is 7700: 11 energies× 35 angles× 10 scattering observables× 2 pp and pn. The final results fora1 ,a2a3 , anda4 are presented in Table 1.To illustrate the fitting quality, we show the unpolarized differential cross section (
dσdΩ ), depolarization tensor for the polarized beam (D), polarization transfer from beam to recoil particle (DT), polarization correlation for initially unpolarized particles (Ayy ), and polarization of scattered particles (P) of pp and pp atTlab= 20, 100, 200, 400, and 800 MeV in Figs. 2–6. AtTlab= 200 MeV, as shown in Fig. 4, because the scattering observables ofTlab= 200 MeV are fitted alone, the scattering observables are described well for both np and pp. As the Coulomb corrections dominate the empirical pp scattering observables of forward and backward angles, the Coulomb corrections are included in the pp scattering observables in this work. The Coulomb corrections are obtained by subtracting theTNN=1 amplitudes of SAID from the pp amplitudes of SAID.Figure 3. (color online) Same as Fig. 2 but for
Tlab = 100 MeV.Figure 4. (color online) Same as Fig. 2 but for
Tlab = 200 MeV.Figure 5. (color online) Same as Fig. 2 but for
Tlab = 400 MeV.Figure 6. (color online) Same as Fig. 2 but for
Tlab = 800 MeV.For energies above and below 200 MeV, because the fits cover a wide range of energies, it is expected that the scattering observables are not described as well as those at 200 MeV. As shown in Figs. 2 and 3, for energies far from 200 MeV (
Tlab= 20 MeV), the unpolarized differential cross section (dσdΩ ) and depolarization tensor for the polarized beam (D) of pn are not perfectly described. However, because the Coulomb corrections of all the pp scattering observables are very significant at energies below 200 MeV, the pp scattering observables are well described. For incident energies above 200 MeV, the results are similar to those of previous studies [11, 12]. As shown in Figs. 5 and 6, the polarization transfer from beam to recoil particle (DT) and polarization correlation for initially unpolarized particles (Ayy ) of pn are not perfectly described. When the Coulomb corrections are included, the pp scattering observables are better described than those of pn. For energies above and below 200 MeV, although some fits are not perfectly described, they are not far from the experiment observables. In conclusion, the RLF model of the NN interaction has been constructed over the laboratory energy range of 20 to 800 MeV. -
To examine the validity of the above fit, we investigate the p+
208 Pb elastic scattering for a wide energy range (Tlab= 30, 65, 200, 295, 500, and 800 MeV). The differential cross sectionsdΩ/dσ , analyzing powersAy , and spinrotation functions Q have been investigated with RLF [6, 7]. However, the RLF model has not been used to describedΩ/dσ ,Ay , and Q for energies above 400 MeV. Moreover, because the parameters of this work are fitted for a wide energy range, the physical quantities may differ somewhat from those in previous studies. For our calculations, we have taken pseudovector coupling for the pion of RLF. Additionally, because the tensor potential is small, it has been neglected in the calculations.The PB corrections have been proved to be significant in RLF at low energy [6, 7]. However, the PB correction factors are still model-dependent. The initial PB parameters are obtained via a relativistic Dirac-Brueckner approach [6]. However, the initial PB parameters cannot well describe the
dΩ/dσ ,Ay , and Q in the energy-dependent RLF model [7]. Therefore, the PB factors of the energy-dependent RLF model are obtained by acquiring the best fits ofdΩ/dσ ,Ay , and Q. Because the fits of this work cover a wide range of energies, the PB factors of the previous energy-dependent RLF model may not be able to describe the p+208 Pb elastic scattering. In order to obtain the PB factors and evaluate our fits, we first compare the optical potentials of this work with the Dirac global optical potentials (GOPs) [1] and the optical potentials of Ref. [7] at 65 MeV. As shown in upper panel of Fig. 7, without considering the PB corrections, the real (Re) values of the optical potentials of this work are similar to those of the previous energy-dependent RLF model, and the imaginary (Im) values of the optical potentials of this work are slightly higher than those of the previous energy-dependent RLF model. These results indicate that the fits of this work are reasonable. The Im values of the optical potentials of both this work and the previous energy-dependent RLF model are higher than those of the GOPs; however, the Re values of the optical potentials of both this work and the previous energy-dependent RLF model are lower than those of the GOPs. The higher Im and lower Re optical potentials imply that the signs of the PB factors of the Im and Re optical potentials in this work should be opposite. At low energy, the extrapolation of Dirac-Brueckner PB factors are too small to properly change the optical potentials. The Im values with PB factors of the phenomenological approach [7] can change the Im values of the optical potentials properly; however, the Re values with PB factors of the phenomenological approach are too low and positive. As shown in the lower panel of Fig. 7, the optical potentials are obtained with the fit of this work and various PB factors. With the extrapolation of Dirac-Brueckner PB factors, the optical potentials do not change much compared with those without the PB corrections. With the PB factors of the phenomenological approach [7], the Im values of the optical potentials are consistent with those of the GOPs, and the Re values of optical potentials are lower than those of GOPs. Therefore, with phenomenological approach, we re-fit the optical potentials by adjusting the PB factors. The PB factors of this work are presented in Table 2. The PB factors of previous studies are also listed, for comparison. The optical potentials with the PB factors of this work are close to the GOPs. When the optical potentials become reasonable, the energy-dependent RLF can well describe the p+208 Pb scattering.Figure 7. (color online) The real (Re) and imaginary (Im) parts of the scalar S and vector V optical potentials for p
+208Pb at 65 MeV. The black solid curves of the upper panel are obtained via the fit of this work (TW) without the Pauli blocking (PB) corrections. The violet dot curves of the upper panel were obtained in Li's work (LI) without PB corrections [7]. The blue dashed curves of the lower panel are obtained by considering the PB corrections of this work. The orange short-dash curves of the lower panel are obtained by considering the PB corrections of Li's work. The cyan dash-dot-dot curves of lower panel are obtained by considering the extrapolation of Dirac-Brueckner (DB) PB corrections [6]. The Dirac GOPs (red dash-dot curves) come from Ref. [1].The Pauli blocking correction factors of this work Energy/MeV Real scalar Imaginary scalar Real vector Imaginary vector 30 −1.30 0.95 −0.85 1.05 65 −0.47 0.85 −0.41 0.88 200 −0.05 0.65 − 0.65 295 − 0.65 − 0.50 500 − − −0.05 − 650 − − −0.1 − 800 − − −0.1 − The phenomenological Pauli blocking factors Energy/MeV Real scalar Imaginary scalar Real vector Imaginary vector 65 0.020 0.68 0.14 0.85 100 0.007 0.60 0.11 0.73 135 0.005 0.53 0.09 0.60 170 0.008 0.42 0.07 0.50 200 0.010 0.35 0.605 0.42 The Pauli blocking correction factors of Dirac-Brueckner approach Energy/MeV Real scalar Imaginary scalar Real vector Imaginary vector 135 0.00377 0.11236 0.08403 0.24535 200 −0.0078 0.098 0.0605 0.207 300 −0.0256 0.0759 0.0243 0.148 400 −0.0434 0.061 −0.0119 0.089 As shown in Fig. 8, the
dΩ/dσ values ofTlab= 65 MeV without PB contributions are lower than the experimental data, andAy and Q fail to describe the experimental data over the entire angular range. With the PB corrections of this work, the experimentaldΩ/dσ ,Ay , and Q ofTlab= 65 MeV can be well described for the entire angular range. Our predictions of 65 MeV are as good as the previous results obtained by changingg2σ andg2ω [7]. ThedΩ/dσ ,Ay , and Q with the PB corrections of previous studies are also shown, for comparison. In this work, the PB factors of the phenomenological approach [7] and the extrapolation of Dirac-Brueckner fail to describe thedΩ/dσ ,Ay , and Q. Generally, different fitting methods need different PB corrections to describe elastic scattering. It should be noted that the PB factors of energies below 200 MeV are far larger than the initial PB parameters generated from the relativistic Dirac-Brueckner approach [6]. In the following, we only consider the PB corrections of this work.Figure 8. (color online) Differential cross sections (
dΩ/dσ ), analyzing powers (Ay ), and spinrotation functions (Q) as a function of the center-of-mass scattering angle for 65 MeV p+208 Pb scattering. The black solid curves do not include the Pauli blocking (PB) corrections. The blue dashed curves include the PB corrections of this work. The orange short-dash curves of the lower panel are obtained by considering the PB corrections of Li's work (LI) [7]. The cyan dash-dot-dot curves are obtained by considering the extrapolation of PB corrections of Dirac-Brueckner (DB) [6]. The open circles indicate the experimental data taken from Ref. [35].For energies far below 200 MeV, as shown in Fig. 9, the
dΩ/dσ andAy ofTlab= 30 MeV without PB contributions lack the property of fluctuating with respect to the angle in the p+208 Pb elastic scattering and cannot describe the experimental data well. When the PB corrections are taken into account, the property of fluctuating with respect to the angle is generated, and thedΩ/dσ andAy of 30 MeV can be consistent with the experimental data. For scattering angles smaller than 120 degrees, our result of 30 MeV shows some improvement over the previous results based on changing the realσN andωN coupling constants (g2σ andg2ω ) [7]. Because the spin observable Q of 30 MeV lacks experimental data, there may be uncertainties in the predictions of 30 MeV.Figure 9. (color online) Differential cross sections (
dΩ/dσ ), analyzing powers (Ay ), and spinrotation functions (Q) as a function of the center-of-mass scattering angle for 30 MeV p+208 Pb scattering. The solid curves do not consider the Pauli blocking (PB) corrections, and the dashed curves consider the PB corrections of this work. The open circles indicate the experimental data taken from Refs. [36, 37].When
Tlab= 200 MeV, as shown in Figs. 4 and 10, even though the np and pp scattering observables are the best fit, the experimentaldΩ/dσ ,Ay , and Q without PB contributions cannot fit the experimental data well. In particular, the spin observable Q is poorly fitted for the entire angular range. When the PB corrections are included, the experimentaldΩ/dσ ,Ay , and Q are fitted well for scattering angles smaller than 40 degrees. The result ofTlab= 200 MeV is similar to that of Ref. [7].For energies above 200 MeV, the situation of PB corrections is complex. When the energy equals 295 MeV, as shown in Fig. 11, the
dΩ/dσ andAy without PB contributions are not well fitted. When a large imaginary factor of PB contributions is used forTlab= 295 MeV, thedΩ/dσ andAy are only slightly corrected for scattering angles smaller than 30 degrees. As the energy increases, the factors of PB contributions become small; however, the effect of PB corrections cannot be neglected. For instance, as shown in Figs. 12–14, whenTlab= 500, 650, and 800 MeV, thedΩ/dσ ,Ay , and Q without PB contributions are not described well. When a small factor of the real vector of PB contributions is included, thedΩ/dσ ,Ay , and Q are fitted better than those in the case without PB contributions. -
The energy-dependent RLF model has been used to fit empirical amplitudes for various energy domains: 50–200, 200–500, and 500–800 MeV. In this work, we try to relate these energy domains within a unified fit. In contrast to previous studies, we fit the scattering observables of pp and pn of WF16 instead of the non-relativistic amplitudes of pp and pn. Because an energy of 200 MeV is chosen as the energy boundary between low and high energy, we first fit the scattering observables of pp and pn of 200 MeV alone. When the scattering observables of pp and pn of 200 MeV are well fitted, the scattering observables of pp and pn of other energies are fitted by using the energy-dependent RLF model. Even though the scattering observables of other energies are not described as well as those of 200 MeV, their trends are consistent with the experimental value, and their values are not far from the experimental values. In conclusion, we constructed the RLF model of the NN interaction over the laboratory energy range of 20 to 800 MeV.
To examine the validity of our fit, we use the energy-dependent RLF model in the range of 20 to 800 MeV to investigate p+
208 Pb elastic scattering forTlab= 30, 65, 200, 295, 500, and 800 MeV. The differential cross sectionsdΩ/dσ , analyzing powersAy , and spin rotation functions Q of p+208 Pb elastic scattering without PB corrections are not well fitted at all the energies considered. Although the scattering observables of pp and pn of 200 MeV best fit the experiment values, the RLF model of 200 MeV without PB corrections fails to describe the experimentaldΩ/dσ ,Ay , and Q. For a better description, the PB corrections must be taken into account. We vary the PB factors to obtain the best fits ofdΩ/dσ ,Ay , and Q. In our calculations, the PB factors are larger than those obtained via a relativistic Dirac-Brueckner approach at energies below 295 MeV but similar to those of phenomenological PB effects [7]. ForTlab≥ 500 MeV, the RLF model with a small factor of the real vector of PB contributions can better describe thedΩ/dσ ,Ay , and Q than those without PB corrections. Our results clearly indicate the importance of PB corrections in the RLF model. -
We thank Ying Kuang and Profs. Zhi-Pan Li and Wei-Zhou Jiang for helpful discussions.
Real parameters | |||||||||
Meson | Isospin | Coupling type | m | ![]() | ![]() | ![]() | ![]() | ![]() | Λ |
σ | 0 | Scalar (S) | 600 | −8.348 | 6.527![]() | 2.773 | −2.741![]() | −7.639 | 962.7 |
ω | 0 | Vector (V) | 782 | 1.019![]() | −6.246![]() | −1.400![]() | 8.492![]() | 2.677![]() | 1155 |
![]() | 0 | Tensor (T) | 550 | 2.794![]() | 3.303 | 1.172![]() | −6.676![]() | −3.008![]() | 1945 |
![]() | 0 | Axial vector (A) | 500 | 4.921![]() | 6.066 | 1.953![]() | −1.089 | −5.133![]() | 1597 |
η | 0 | Pseudoscalar (P) | 950 | 1.180![]() | −1.694![]() | −1.350![]() | 5.417 | 5.550![]() | 957.3 |
δ | 1 | Scalar (S) | 500 | 6.095![]() | 5.635![]() | 1.287 | −6.040![]() | −4.101 | 2945 |
ρ | 1 | Vector (V) | 770 | −1.667![]() | 9.420![]() | 1.644![]() | −2.905![]() | 2.915![]() | 3001 |
![]() | 1 | Tensor (T) | 600 | −2.500![]() | −1.893 | −7.848 | 4.448![]() | 1.989![]() | 1280 |
![]() | 1 | Axial vector (A) | 650 | −1.356 | −4.678 | −1.529![]() | 8.102![]() | 4.146![]() | 745.1 |
π | 1 | Pseudoscalar (P) | 138 | 1.095![]() | −3.834![]() | −4.889 | 4.005![]() | 1.082![]() | 677.4 |
Imaginary parameters | |||||||||
Meson | Isospin | Coupling type | ![]() | ![]() | ![]() | ![]() | ![]() | ![]() | ![]() |
σ | 0 | Scalar (S) | 600 | −2.849 | −4.637 | −1.475![]() | 6.301![]() | 4.377![]() | 800.0 |
ω | 0 | Vector (V) | 700 | 4.405 | 4.825 | 2.077![]() | −9.616![]() | −5.672 ![]() | 716.0 |
![]() | 0 | Tensor (T) | 750 | −9.457![]() | −1.474 | −7.535 | 4.140![]() | 1.914![]() | 1340 |
![]() | 0 | Axial vector (A) | 750 | −2.198 | −3.583 | −1.374![]() | 6.840![]() | 3.762![]() | 1437 |
η | 0 | Pseudoscalar (P) | 1000 | 1.420![]() | −1.619![]() | −2.011![]() | 9.779 | 6.591![]() | 2758 |
δ | 1 | Scalar (S) | 650 | 3.094 | 1.593 | 2.945 | 1.614![]() | −1.328![]() | 751.9 |
ρ | 1 | Vector (V) | 600 | −2.421 | −9.543![]() | −3.644 | 9.074![]() | 1.292![]() | 766.6 |
![]() | 1 | Tensor (T) | 750 | 6.140![]() | 1.513 | 5.775 | −3.011![]() | −1.557![]() | 1662 |
![]() | 1 | Axial vector (A) | 1000 | 2.265 | 5.071 | 1.997![]() | −1.011 | −5.450![]() | 1233 |
π | 1 | Pseudoscalar (P) | 500 | −5.891 | 2.029![]() | 2.363![]() | −9.593![]() | −7.574![]() | 1342 |