The nucleon single-particle energies (SPEs) of the selected nuclei, that is, O16, Ca40, and Ni56, are obtained by using the diagonal matrix elements of two-body effective interaction, which generated through the lowest-order constrained variational (LOCV) calculations for the symmetric nuclear matter with the Aυ18 phenomenological nucleon-nucleon potential. The SPEs at the major levels of nuclei are calculated by employing a Hartree-Fock inspired scheme in the spherical harmonic oscillator basis. In the scheme, the correlation influences are taken into account by imposing the nucleon effective mass factor on the radial wave functions of the major levels. Replacing the density-dependent one-body momentum distribution functions of nucleons, n(k,ρ), with the Heaviside functions, the role of n(k,ρ) in the nucleon SPEs at the major levels of the selected closed shell nuclei is investigated. The best fit of spin-orbit splitting is taken into account when correcting the major levels of the nuclei by using the parameterized Wood-Saxon potential and the Aυ18 density-dependent mean field potential which is constructed by the LOCV method. Considering the point-like protons in the spherical Coulomb potential well, the single-proton energies are corrected. The results show the importance of including n(k,ρ), instead of the Heaviside functions, in the calculation of nucleon SPEs at the different levels, particularly the valence levels, of the closed shell nuclei.
1. Introduction
The study of shell structure and relevant properties, for example, single-particle energies (SPEs), of closed shell nuclei plays a crucial role in the description of various astrophysical scenarios. For many years, the impact of nuclear shell structure on the astrophysical nucleosynthesis processes has been recognized by observing the so-called equilibrium-, s-, and r-process peaks in the galactic abundance distribution curves [1]. SPEs have frequently been used in a few mass models where the total ground-state energies of finite nuclei are calculated as the sum of a microscopic term, determined by the calculated SPEs, and a macroscopic term, based on the finite range droplet model [2, 3]. Nowadays, the benchmark role of closed shell nuclei properties on a reliable description of astrophysical events becomes more obvious by daily progress in getting the observational data, for example, from the r-process sources such as binary compact stars (white dwarf and neutron star mergers) and explosive supernovae (Ia and II), or from the weak interaction processes in the stellar burning, whether via neutron/electron capture or via neutrino absorption on nucleons/nuclei (see, e.g., [4, 5] and the relevant references therein). Hence, investigating shell structure properties of closed shell nuclei such as SPEs is necessary for the state-of-the-art of the nuclear astrophysics calculations.
Recently, many works have been presented to describe the shell structure properties of exotic nuclei based on the closed shell nuclei properties. We can point, for instance, to the performing systematic microscopic calculations of medium-mass nuclei [6, 7], the calculations based on microscopic valence-space interactions [8–16], and the calculations by using Self-Consistent Green’s Function [17, 18] or In-Medium Similarity Renormalization Group [19, 20]. An important step to study nuclei far from stability is the use of nuclear interactions predicted by chiral effective field theory, which provides consistent two- and three-nucleon forces without any fit to the particular region in the nuclear chart [21, 22]. It is worth mentioning that shell evolutions of various neutron-rich isotopes have been studied by many authors recently, for example, for Si isotopes [23], for Sn isotopes with A=103–110 [24], for Ca isotopes [7], and for single-particle-like levels in Sb and Cu isotopes [25]. Lately, the authors of [26] studied the nuclear states with single-particle and collective characters in the framework of beyond mean field approaches and in the study of semimagic S44. The structure obtained is compatible with results from the state-of-the-art of shell-model calculations and experimental data [27–29]. More recently, by applying the newly developed theory of effective nucleon-nucleon interactions on the study of exotic nuclei, the authors of [30] performed shell-model calculations with several major shells for the exotic neutron-rich Ne, Mg, and Si isotopes and predicted the drip lines and showed the effective single-particle energies by exhibiting the shell evolution. Thus, accounting for the shell structure of exotic nuclei, the nuclei far from stability, is essential to both nuclear physics and astrophysics.
Although most of the attempts, mentioned above, lead to results which are compatible with the relevant experimental data, they are not the end of the story in the study of the shell structure of nuclei because of the highly complicated nature of the strong nuclear force. During the past decade, many authors investigated [15, 31, 33–45] how the presence of strong short-range repulsive force in finite nuclei leads to the presence of a high-momentum component in the ground-state wave function, as predicted by scattering experiments. Recently, the high-momentum transfer measurements have shown that nucleons in nuclear ground states can form short-range correlated pairs with large relative momentum [45–50]. Furthermore, high-energy electron scattering shows that internucleon short-range repulsive force forms correlated high-momentum neutron-proton pairs and thus, in exotic nuclei, protons are more probable to have a momentum greater than the Fermi momentum in comparison to neutrons [51]. Nevertheless, the effect of the high-momentum component of nuclear force on the properties of nucleons such as nucleon momentum distribution, n(k), and SPE is not well-accommodated in the shell-model theory since this model does not consider the nucleon-nucleon interaction at small distances with precision. For many years, some authors have theoretically predicted the experimental features of n(k), such as the long-tail component in the high-momentum region and the gap at the Fermi surface compared to the fully occupied case of the noninteracting Fermi gas model [52, 53]. However, there are still differences, particularly in the high-momentum region, between experimental n(k) data and those on the basis of theoretical models whether they are microscopic or phenomenological. It is worth mentioning that in the framework of the lowest-order constrained variational (LOCV) method [54] and based on the Ristig-Clark formalism [55] the authors of [56] constructed one-body density- and momentum-dependent distribution functions n(k,ρ). Their calculated n(k,ρ) have been in overall agreement with those of other methods. By applying n(k,ρ) on the calculation of single-nucleon properties, the authors obtained comparable results with other methods [57]. The model-dependent nature of n(k,ρ), particularly in the high-momentum part, leads to disagreements between different methods in the calculation of single-nucleon properties. Hence, the investigation of the influence of high-momentum component on the shell structure properties of nuclei, for example, SPEs, is an open problem in nuclear physics.
In this work, we use the LOCV method to generate n(k,ρ) and two-body effective interaction matrix elements. LOCV is a fully self-consistent technique without any free parameter and with state-dependent correlation functions. The method uses the cluster expansion formalism with an adequate convergence up to two-body interactions. The insignificant effect of higher-order terms, for example, three-body, on the nuclear matter ground-state energy, and the convergence of the cluster expansion in the LOCV formalism have been discussed by the authors of [58] and the formalism has been validated. The method has been employed for asymmetric/symmetric nuclear matter (A/SNM) and neutron matter calculations at zero and nonzero temperatures, for different types of neutron stars [59–67], and for finite nuclei calculations from which results have been compatible with experimental data [68–75]. Lately, in order to reduce the effect of the high-momentum component of strong nuclear force on finite nuclei properties, a Fermi momentum cut-off has been imposed on both average effective interactions and density- and channel-dependent ones [76, 77]. More recently, we investigated the effect of the gap of n(k,ρ) at the Fermi surface on the binding energy of the selected medium-mass nuclei [78]. The results of our last work motivate us to investigate the effect of n(k,ρ) on the calculation of the nucleon SPEs of the closed shell nuclei.
There is a problem in the definition and the calculation of single-particle levels (SPLs) for atomic nuclei comprising strongly interacting nucleons. Many years ago, in order to describe experimental measurements, some authors (see, e.g., [79]) attempted to calculate the SPLs of nuclei by presenting a reasonable definition of SPLs. Furthermore, the presence of correlations between interacting nucleons has not been taken into account in the Hartree-Fock (HF) mean field levels [80]. In order to calculate the nucleon SPEs of nuclei, we present a formalized scheme, called here the HF inspired scheme, in which the two-body correlations are considered. We include nk,ρ in the calculation of density-dependent SPEs and impose a factor of nucleon effective mass, M∗ at the Fermi surface, on the radial parts of density functions when extracting SPLs. Thus, we control high-momentum components of two-nucleon interactions by nk,ρ, rather than using the usual Heaviside functions, and account for the presence of correlations in HF mean field by M∗. Both nk,ρ and M∗ are generated by the LOCV calculations for SNM with Aυ18(Jmax=2) phenomenological nucleon-nucleon potential [75]. In order to do this, we organize the paper as follows. In Section 2, we formalize a variational calculation of nucleon SPEs at the major levels of nuclei in the spherical harmonic oscillator (SHO) basis. During the calculation, we investigate the role of nk,ρ on the calculations at the minimum point of the total energy of a nucleus, denoted by γmin (γ is the SHO parameter). In Section 3, we correct the major levels of nuclei by applying the phenomenological parameterized Wood-Saxon (WS) potential and the Aυ18 mean field potential on the calculation of spin-orbit interaction at γmin. In Section 4, by including the Coulomb interactions, we improve the single-proton energies. In Section 5, we shall show the effect of nk,ρ on the calculation of SPEs in the ground-state of selected closed shell nuclei, that is, O16, Ca40, and Ni56. Finally, in Section 6, we discuss the results.
2. Calculation of Nucleon SPEs at Major Levels
Although one can obtain the nucleon SPEs of a nucleus by solving the HF self-consistent equations, we calculate them by using a variational method and introducing the HF inspired scheme. In our point of view, applying a variational method for the present calculations is a more straightforward approach to show the effect of nk,ρ on SPEs. On the other hand, we believe that including a factor of M∗/M with M as the nucleon mass within the radial parts of SHO functions covers the correlation effects on the mean field calculations. Now, the general form of the Hamiltonian of a nucleus comprising A nucleons is given by(1)H=∑i=1ATi+12∑ijAVij,where T(i), the single-particle kinetic energy (SPKE) operator for ith nucleon, is given by ℏ2ki2/2M and V(ij) denotes the two-nucleon interaction potential. By denoting the occupied single-particle states by c,c′,…, the total energy of occupied states (the core energy), Ec, in diagonal configuration, is given by(2)Ec=∑cocccTc+12∑cc′occcc′Vcc′,where “occ” shows that the summation topped to occupied states. Now, we define the nucleon single-particle interaction energy (SPIE) as the energy which a specific nucleon obtains via the mean field interaction of all other nucleons in the core. By denoting configuration c as the state of specific nucleon and SPIE as Uc=1/2∑c′occcc′Vcc′ and rewriting interaction part of the Hamiltonian, the total energy of core gets the following form:(3)Ec=∑cocccTc+∑cocccUcc.Now, we can present an expression to calculate the SPE of a nucleon in the core as follows:(4)εc=cTc+cUcc.The total energy of a closed shell nucleus with configuration c plus one nucleon (one absent-nucleon, so-called, one hole) in the state p above (below) the Fermi surface is as follows:(5)Ec±p=Ec±pTp±∑cocccpVcp.By defining Up=∑cocccpVcp, we can rewrite the total energy as follows:(6)Ec±p=Ec±pTp±pUpp.The 1p-1h SPE (near the Fermi surface) reads(7)εph=pTp+pUpp.In this work, we set the SHO basis; that is, c=α;γ in which α=nlml with n, l, and ml as the principle, the orbital angular momentum, and the projection of orbital angular momentum quantum numbers of the nucleon in the core, respectively. The SHO parameter, which plays the role of variation parameter, is given by γ=Mω/ℏ in which ℏω is the oscillator energy to measure the size of the specific nucleus.
In order to calculate SPEs, our strategy is as follows: first, we obtain the density- and momentum-dependent diagonal matrix elements of the nucleon-nucleon effective interactions, Heff12, by using the LOCV calculations for SNM with Aυ18(Jmax=2) potential. Then we employ the elements in order to construct SPIE. Next, while we keep the quantum numbers of one orbit, αc(p), we do summations over the relevant quantum numbers of other orbits in all over occupied states. In order to obtain the expectation values of SPKE and SPIE in configuration space, we use the correlated wave functions by imposing the nucleon effective mass on the radial SHO functions. Finally, we find the minimum value of the total energy of the nucleus, (3), via variation of γ and thus the ground-state energy of the nucleus is obtained as we did in our previous works [69–78]. We assign the corresponding γ at the minimum point of total energy of the nucleus by γmin and calculate the nucleon SPEs of the nucleus at γmin. During the calculation of SPEs, we shall investigate how the one-body density-dependent momentum distribution functions are important in the evaluation of SPEs by replacing n(k,ρ) with Heaviside functions.
According to (3) and above explanations, the ground-state energy of the core reads(8)Eg.scγmin=∑αctαc;γmin+∑αcuαc;γmin,where tαc;γmin=αc;γminTαc;γmin0020and uαc;γmin=αc;γminUcαc;γmin stand for SPKE and SPIE, respectively. Thus, the SPE of nucleon in a specific single-particle state of a closed shell nucleus at γmin is given by(9)εαc=tαc;γmin+uαc;γmin,and the 1p-1h SPE is as follows:(10)εαph=tαp;γmin+uαp;γmin,where (11a)tαp;γmin=αp;γminTαp;γmin,(11b)uαp;γmin=αp;γminUpαp;γmin.In the SHO basis, one can obtain SPIE from (4) and (7) as follows:(12)ucpγ=αcp;γUcpαcp;γ=12∑α2αcpα2;γHeff12αcpα2;γ,where 1/2 means that we need a factor of 1/2 only in the calculation of SPIE for αc (cf. (2)). Regarding the LOCV method, the two-body effective interaction has the following form [54]:(13)Heff12=-ℏ22MF12,∇122,F12+F12V12F12,where F12, the density- and channel-dependent two-body correlation functions, are operator type-like Aυ18 [73]. The operator V(12) for the Aυ18 potential has the following form [81]:(14)V12=∑λυλ12Oλ12,where the forms of eighteen components of the two-body operator Oλ=1-18(12) are as follows:(15)Oλ=1-1812=1,σ1·σ2,τ1·τ2,σ1·σ2τ1·τ2,S12,τ1·τ2S12,L·S,τ1·τ2L·S,L2,σ1·σ2L2,σ1·σ2τ1·τ2L2,L·S2,τ1·τ2L·S2,T12,S12T12,τz1·τz2.In (14), all components υλ are central functions of distance between two interacting nucleons without any dependency on quantum numbers.
Now, the matrix elements of two-body interactions are obtained by sandwiching Heff12 (a) between ∑k→1k→2,η1η2k→1k→2,η1η2 and ∑k→1′k→2′,η1′η2′k→1′k→2′,η1′η2′, two complete basis sets of plane waves together with spin and isospin parts, and (b) between ∫d3rr→R→r→R→ and ∫d3r′r→′R→′r→′R→′, two orthogonal bases in the configuration space, and (c) between ∑λλλ and ∑λ′λ′λ′, two complete sets in which λ denotes L, S, T, MT, and J where L, S, T, and MT are the relative orbital momentum, the total spin, isospin, and isospin projection of two interacting nucleons, respectively, and J comes from J→=L→+S→. Thus, a straightforward calculation, by taking just diagonal matrix elements and considering orthogonality conditions, (12) is simplified by(16)uαcpρ,γ=12Ωγ2π3∫dk→1nk1,ρϕαcpk→1,γ-12Uαcpk1,ρ,γ,where(17)Uαcpk1,ρ,γ=12π2∑α2,mτ2∑λCmτ1mτ2TMT22J+11--1L+S+T2∫dk→2nk2,ρϕα2k→2,γ-12∫Heffλr,ρjLrk2r2dr.In (17), Cmτ1mτ2TMT are the isospin Clebsch-Gordan coefficients and jL(rk) are the well-known spherical Bessel functions in which the relative momentum and distance of two interacting nucleons are obtained by the familiar relations k=k→1-k→2/2 and r=r→1-r→2, respectively. In the above equations, we converted the summations over k→1 and k→2 to integrals by the well-known relation as follows:(18)1Ω∑k→⟶12π3∫nkdk→,where n(k) are the one-body momentum distribution functions. In this work, n(k), generated by using the LOCV calculations for SNM with Aυ18(Jmax=2) potential based on the Ristig-Clark formalism [55], are density-dependent [78]. We suppose that the volume of box, Ω, is obtained by Ωγ=4π[Rrms(γ)]3/3 where Rrms(γ) calls(19)Rrms2γ=∫ρR,γR4dR.In (16) and (17), we can replace ϕαk→,γ-1 by 2l+1Rnlk,γ-12/4π, since SPKE and SPIE depend only on the magnitude of the relevant single-particle momentum and also the interest of this work is to consider the closed shell nuclei. It should be noted that k’s are topped to the Fermi momentum, KFρ=6π2ρ/ν1/3, with ν being degeneracy parameter of nucleons. Now, by use of (9), (10), (16), and (17), and recalling O=trρO, we find SPE as follows:(20)εnlcpγ=∫dRR2R~nlR,γ2tnlcpρ,γ+unlcpρ,γ,where we replaced α with nl. In the above equation, R~nlR,γ reads(21)R~nlR,γ=CnlM∗R/MRnlR,γ,where Cnl is a normalization coefficient for each natural orbit and M∗(R), the effective mass of nucleon, is calculated at the Fermi surface by the LOCV method with considering Aυ18(Jmax=2) potential [78]. In this work, we calculate M∗ at a fixed energy, so-called k-mass definition, which is a measure of the nonlocality of single-particle potential, due to nonlocality in space. This definition is clearly different from the E-mass definition which shows the nonlocality of single-particle potential, due to nonlocality in time. By imposing this factor we mean the shape of density distribution function is altered by correlations in the HF inspired scheme. After varying γ and finding the minimum total energy of nucleus, we obtain the SPE of a specific single-particle state at γmin(εnlcp(γmin)). Table 1 shows εnl(γmin) for the single-nucleon major levels of O16, Ca40, and Ni56. In order to clarify the importance of including n(k,ρ) in calculations, the values of εnl are considered in two different cases: (i) including the constructed n(k,ρ) and (ii) including the Heaviside function. In Table 1, the magnitudes of SPEs at the major levels, especially near the valence levels, indicate the importance of including n(k,ρ), as discussed in Section 6.
The single-nucleon major levels of energies (in MeV) of the sample nuclei at γmin (in fm−1) by including two distribution functions: (i) n(k,ρ) and (ii) Heaviside (see text).
Nuclei
Distribution
γmin
ε0s
ε0p
ε0d
ε1s
ε0f
ε1p
16O
(i)
0.310
−23.003
−11.100
−5.404
−14.203
—
—
(ii)
0.330
−26.483
−35.336
−34.465
−21.434
—
—
40Ca
(i)
0.330
−22.273
−11.939
−7.299
−15.946
+0.674
−1.127
(ii)
0.360
−19.140
−30.205
−36.423
−23.177
−29.531
−18.033
56Ni
(i)
0.340
−19.094
−13.321
−7.644
−13.661
−0.581
−4.993
(ii)
0.360
−17.219
−31.566
−36.567
−19.694
−29.535
−20.029
3. Spin-Orbit Interaction
In order to complete the calculations of SPLs and the evaluation of the n(k,ρ) effect on them, we include the spin-orbit effect in our calculations. There are many calculations of spin-orbit corrections, from the old calculations to new relativistic ones [82–87]. Considering spin-orbit interaction, the spin orientation degeneracy is removed in the major SPLs and they will be split. Indeed, the SHO base kets, α,γ, are not appropriate in the presence of spin-orbit interaction. Including total angular momentum quantum number, j, the appropriate base kets are nlsj,mj,γ, where mj=ml+ms and ms is the spin projection. Thus, SPEs at γmin, corresponding γ to the minimum point of total energy of the nucleus, are given by(22)εnljγmin=εnlγmin+Δεnljγmin,where(23)Δεnljγmin=nlsj,mj,γminV^lsnlsj,mj,γmin.Sandwiching the spin-orbit operator, V^ls~l→·s→, between two complete orthogonal sets of position base kets, imposing effective mass factor on the radial SHO functions, and performing straightforward mathematics, the spin-orbit correction is written as(24)Δεnljγmin=υSO∫dfRdRR~nlR,γmin2RdR-l+12;j=l-12,l2;j=l+12.In (13), the forms of constant factor υSO and the WS type function fR are as follows [82]:(25a)υSO=λV0ℏ/μc22,(25b)fR=1+expR-R0a-1,where λ is a fit parameter, V0 presents the total strength, μ is the reduced mass of the single particle and residual nucleus mass, R0 is given by R0=r0A1/3, and a is the skin depth which shows diffuseness of the surface. In order to achieve the best agreement with the experimental spin-orbit splitting of nuclei the parameters take the following values at γmin:(26)V0,best=55.00MeV;r0,best=1.25fm;abest=0.66fm;λbest=95.00.
As another approach, we use our density-dependent single-particle potential, uαcpρ,γ, as fR. In this way, we fit the spin-orbit correction with the experimental data only by one parameter, that is, λ, through υSO=λℏ/μc2/2. In order to obtain the best agreement with the experimental spin-orbit splitting of nuclei, the parameter λbest takes the value 0.528 at γmin. Table 2 shows the results of spin-orbit interaction corrections on the major SPLs of O16, Ca40, and Ni56, calculated at γmin, according to the best-fit values for the parameters of WS and u(ρ,γ) potentials. We put the values of spin-orbit splitting corresponding to the Aυ18 mean field potential, (16), inside the parentheses. In order to characterize the accuracy of the best-adjustment data in the description of the calculated spin-orbit splitting, δnlcal=Δεnl,j=l-1/2-Δεnl,j=l+1/2, the average theoretical error of predictions is given by(27)ADεnl=1N∑i=1Nδnl,ical-δnl,iexp2,where δnl,iexp is the empirical spin-orbit splitting of each nucleus. Regarding the lack of empirical spin-orbit splitting data for some nuclei, as seen in Table 2, the value of N is 14 in the above relation.
Spin-orbit corrections (in MeV) on the single-neutron and -proton levels of 16O, 40Ca, and 56Ni at the corresponding γmin (in fm-1). The corrections have been obtained by employing WS and u(ρ,γ) potentials, where the relevant values of latter potential have been put inside parentheses. The average error has been calculated by considering the empirical spin-orbit splitting [31, 32].
Nuclei
γmin
δε0pcal
δε0dcal
δε0fcal
δε1pcal
δ0pexp
δε0dexp
δε0fexp
δε1pexp
ADεnl
16O-n
0.31
6.05(2.69)
4.89(4.61)
—
—
6.18
5.08
—
—
0.37(0.75)
16O-p
0.31
6.07(2.70)
4.90(4.63)
—
—
6.32
5.04
—
—
40Ca-n
0.33
5.10(2.55)
7.22(4.40)
5.92(6.24)
2.32(2.49)
—
5.63
5.71
2.00
40Ca-p
0.33
5.11(2.56)
7.24(4.42)
5.94(6.27)
2.32(2.50)
—
5.40
5.69
1.75
56Ni-n
0.36
3.38(2.56)
6.42(4.42)
6.50(6.30)
1.26(2.54)
—
—
7.17
1.11
56Ni-p
0.36
3.39(2.56)
6.44(4.44)
6.52(7.45)
1.27(2.55)
—
—
7.45
1.11
4. Coulomb Interaction
To impose the electromagnetic effects on the single-proton levels, one may suppose that each proton lies in the repulsive Coulomb potential well, made by the other protons. In fact, the potential is determined by a certain nuclear charge distribution, ρCR, as follows:(28)υCR=eR∫0RρCR′R′2dR′;R′≤R,e∫R∞ρCR′R′dR′;R′>R,where e=1.6×10-19 C. In this work, the protons are assumed to be point-like particles and to have the uniform spherical charge distribution with the total value Ze. Then the Coulomb potential is given by(29)υCR=Ze23-R/RC22RC;R≤RC,Ze2R;R>RC.In the above equation, RC, the spherical charge distribution radius is calculated by RC=5/3Rrms(γmin).
Comparing to nuclear interactions, there is no significant difference if one supposes that the protons are point-like or they have certain radii. In the latter case, it is sufficient to use RC2=Rrmsγmin2+Rint2 in which Rint≈0.8 fm [88]. Table 3 shows the results of imposing Coulomb interaction corrections on the single-proton levels, εC,nl, of the selected nuclei at the corresponding γmin by including n(k,ρ) and the Heaviside functions, characterized with (i) and (ii), respectively.
The Coulomb interaction corrections (in MeV) on the single-proton levels of 16O, 40Ca, and 56Ni at the corresponding γmin (in fm−1) by including two functions: (i) n(k,ρ) and (ii) Heaviside.
Nuclei
γmin
εC,0s
εC,0p
εC,0d
εC,1s
εC,0f
εC,1p
16O (i)
0.310
0.338
0.714
0.737
0.166
—
—
16O (ii)
0.330
0.365
0.794
0.920
0.171
—
—
40Ca (i)
0.330
0.367
0.785
0.920
0.175
0.741
0.547
40Ca (i)
0.360
0.402
0.885
1.135
0.198
1.121
0.524
56Ni (i)
0.34
0.387
0.807
0.988
0.186
0.864
0.533
56Ni (ii)
0.36
0.409
0.877
1.133
0.202
1.121
0.516
5. The Effect of nk,ρ on SPEs
Including the spin-orbit corrections which have been calculated by WS potential, ΔεnljWS, and Coulomb corrections, εC,nl, the nucleon SPEs, εnlj, read(30)εnljγmin=εnlγmin+ΔεnljWSγmin+εC,nlγmin.Table 4 shows the calculated SPLs of O16, Ca40, and Ni56 at the corresponding γmin, obtained by the HF inspired scheme and the LOCV calculations for SNM with Aυ18(Jmax=2) by including nk,ρ and the Heaviside functions, identified by (i) and (ii), respectively, besides the experimental data. The letters p, n and the abbreviated word “exp” specify the single-proton, single-neutron, and the experimental levels of each nucleus, respectively. As shown in Table 4, including n(k,ρ) leads to a reduction in the magnitude values of all SPEs with respect to the considering the Heaviside functions. This reduction is more significant for the valence levels. We discuss in Section 6 the results.
The nucleon SPEs of 16O, 40Ca, and 56Ni by considering (i) n(k,ρ) and (ii) Heaviside functions.
Nucleus scheme
ε0s1/2
ε0p3/2
ε0p1/2
ε0d5/2
ε0d3/2
ε1s1/2
ε0f7/2
ε0f5/2
ε1p3/2
ε1p1/2
16O-p (i)
−22.7
−12.4
−6.3
−6.6
−1.7
−14.0
—
—
—
—
16O-p (ii)
−26.1
−36.9
−29.9
−36.1
−29.7
−21.2
—
—
—
—
16O-n (i)
−23.0
−13.1
−7.1
−7.4
−2.5
−14.2
—
—
—
—
16O-n (ii)
−26.5
−37.7
−30.7
−37.0
−30.7
−21.4
—
—
—
—
40Ca-p (i)
−21.9
−12.8
−7.7
−9.3
−2.0
−15.7
−1.1
+4.8
−1.3
+0.9
40Ca-p (ii)
−18.7
−31.0
−25.9
−38.8
−30.0
−23.0
−32.2
−23.3
−18.2
−16.1
40Ca-n (i)
−22.3
−13.6
−8.5
−10.2
−3.0
−15.9
−1.8
+4.1
−1.9
+0.4
40Ca-n (ii)
−19.1
−31.9
−26.8
−39.9
−31.2
−23.2
−33.4
−24.4
−18.7
−16.7
56Ni-p (i)
−18.7
−13.6
−10.2
−9.2
−2.8
−13.5
−2.5
+4.0
−4.9
−3.6
56Ni-p (ii)
−16.8
−31.7
−28.6
−38.2
−31.2
−19.5
−32.9
−23.6
−19.9
−18.8
56Ni-n (i)
−19.1
−14.4
−11.1
−10.2
−3.8
−13.7
−3.4
+3.1
−5.4
−4.1
56Ni-n (ii)
−17.2
−32.6
−29.4
−39.3
−32.4
−19.7
−33.1
−24.8
−20.4
−19.3
6. Discussion and Conclusion
The momentum and energy distributions of nucleons are affected by the strong short-range correlations of nucleons in finite nuclei. This fact becomes more important when we calculate the nuclei properties, for example, nucleon SPEs, which play a main role in the description of some astrophysical processes (cf. Section 1). Thus, it would be useful if we measure how much n(k,ρ) can influence the prediction of nucleon SPEs of closed shell nuclei. In this work, we have investigated the issue for the selected closed shell nuclei, that is, O16, Ca40, and Ni56, by replacing n(k,ρ), generated by the LOCV calculations for SNM with Aυ18(Jmax=2), with Heaviside functions. By constructing a density- and momentum-dependent mean field potential in the framework of LOCV method and by employing the HF inspired scheme in which the nucleon effective mass has been imposed on the natural orbits of nuclei, we have calculated the nucleon SPEs at the major levels of the selected nuclei and measured the effect of nk,ρ on the calculations. Table 1 illustrates the nucleon SPEs at major levels of the selected nuclei by considering nk,ρ and Heaviside functions. We see that there is an absolute increase around 70%, 30%, and 100% in the energy of valence levels of O16, Ca40, and Ni56, respectively, when we replace nk,ρ with the Heaviside functions. In order to obtain the correct values of the subshell energies, we have included the spin-orbit splitting corrections by using the parameterized WS potential and our density-dependent single-particle potential, uρ,γ, as the Thomas type mean field potential. The scheme incorporating this effect has been considered at γmin, according to the best fit to the empirical spin-orbit splitting data, and the results have been presented in Table 2. In the case of WS best-fit, the calculated average deviation, ADεnl, from empirical data is surprising with respect to the available data of empirical spin-orbit splitting for both neutron and proton levels of all selected nuclei. Although there is an advantage in the use of the Aυ18 mean field potential relative to the WS potential because of the existing only one adjustment parameter (λ), the spin-orbit splitting energies are not compatible with the experimental data for 0p levels of the selected nuclei. In order to correct single-proton levels, we assumed that single point-like proton lies in a uniformly charged spherical potential well with a radius determined by the corresponding RMS radius, Rrmsγmin. As shown in Table 3, whereas these corrections do not lead to a significant change in the energy levels of the nuclei, there is a maximum about 20% and 30% decrease for the energy levels 0d (all selected nuclei) and 0f (Ca40), respectively, by replacing the Heaviside functions with nk,ρ. Table 4 compares the results of two cases of the SPEs calculations by considering two momentum distribution functions: (i) nk,ρ and (ii) Heaviside at the corresponding γmin. Regarding the magnitude of SPEs, considering n(k,ρ) leads to a significant reduction with respect to the Heaviside case for the valence levels of the selected closed shell nuclei. The values of SPEs depend on the calculation scheme, which encompasses the way the interactions are regularized, and details of the many-body method. However, the results of the employed scheme in the present work show how the selection of momentum distribution function is important in the SPE calculations. As a result of the above discussions, by taking into account the depletion at the Fermi surface and the tail of the one-body density- dependent momentum distribution in the calculations, through using n(k,ρ), we find a large difference in the nucleon SPEs at the valence levels and other levels of the selected closed shell nuclei.
In conclusion, n(k,ρ) has been employed to calculate the SPEs of closed shell nuclei, that is, O16, Ca40, and Ni56, by using the matrix elements of nucleon-nucleon effective interaction and n(k,ρ) from the LOCV calculations for SNM with Aυ18 and also employing an HF inspired scheme to consider the correlation effects through imposing the nucleon effective mass on the density distribution functions of nucleons. The surprising results show that the use of n(k,ρ) has a significant effect on the calculation of SPEs. The results of the present formalism might be improved by including all interaction matrix elements and constructing the diagonalized interaction matrix, by including three-body force effects, as well as the density-dependent Coulomb interactions. Although the calculations have been limited by the LOCV method (to generate the two-body effective interaction matrix elements), the Ristig-Clark formalism (to construct n(k,ρ) functions), and the HF inspired scheme (to include correlation effects), the results show how much the inclusion of n(k,ρ) can be important in the prediction of closed shell properties such as single-nucleon energies.
Competing Interests
The author declares that he has no competing interests.
Acknowledgments
The author wishes to thank the Laboratory for Advanced Computing at the University of Coimbra for providing high performance computing resources in order to calculate the two-body effective interactions matrix elements. The author would also like to thank Professor Constança Providência for her hospitality at the Centro de Física Computacional, Departamento de Física, Universidade de Coimbra, Portugal, Professor Orfeu Bertolami for his hospitality at the Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Portugal, Professor I. Vidaña for their good discussions, and Dr. S. Chiacchiera for her useful help on the use of the cluster.
BurbidgeE. M.BurbidgeG. R.FowlerW. A.HoyleF.Synthesis of the elements in stars195729454765010.1103/revmodphys.29.5472-s2.0-36149014927von GrooteH.HilfE. R.TakahashiK.A new semiempirical shell correction to the droplet model: gross theory of nuclear magics1976175-641842710.1016/0092-640x(76)90031-02-s2.0-0017123907MöllerP.NixJ. R.MyersW. D.SwiateckiW. J.Nuclear ground-state masses and deformations199559218538110.1006/adnd.1995.10022-s2.0-1842850439GraweH.LangankeK.Martínez-PinedoG.Nuclear structure and astrophysics2007709, article R021525158210.1088/0034-4885/70/9/r022-s2.0-34548427906LangankeK.ThielemannF.-K.WiescherM.Al-KhaliliJ.RoecklE.Nuclear astrophysicsand nuclei far from stability2004651383467Lecture Notes in Physics10.1007/978-3-540-44490-9_11OtsukaT.SuzukiT.HoltJ. D.SchwenkA.AkaishiY.Three-body forces and the limit of oxygen isotopes2010105303250110.1103/physrevlett.105.032501HoltJ. D.OtsukaT.SchwenkA.SuzukiT.Three-body forces and shell structure in calcium isotopes201239808511110.1088/0954-3899/39/8/085111GallantA. T.BaleJ. C.BrunnerT.ChowdhuryU.EttenauerS.LennarzA.RobertsonD.SimonV. V.ChaudhuriA.New precision mass measurements of neutron-rich calcium and potassium isotopes and three-nucleon forces201210903250610.1103/physrevlett.109.032506WienholtzF.BeckD.BlaumK.BorgmannC.BreitenfeldtM.CakirliR. B.GeorgeS.HerfurthF.HoltJ. D.KowalskaM.KreimS.LunneyD.ManeaV.MenéndezJ.NeidherrD.RosenbuschM.SchweikhardL.SchwenkA.SimonisJ.StanjaJ.WolfR. N.ZuberK.Masses of exotic calcium isotopes pin down nuclear forces2013498745434634910.1038/nature122262-s2.0-84879396564HoltJ. D.MenéndezJ.SchwenkA.Chiral three-nucleon forces and bound excited states in neutron-rich oxygen isotopes201349, article 3910.1140/epja/i2013-13039-2CaesarC.SimonisJ.AdachiT.Beyond the neutron drip line: the unbound oxygen isotopes 25O and 26O201388303431310.1103/physrevc.88.034313HoltJ. D.MenéndezJ.SchwenkA.Three-body forces and proton-rich nuclei20131102502250210.1103/PhysRevLett.110.022502HoltJ. D.MenéndezJ.SchwenkA.The role of three-nucleon forces and many-body processes in nuclear pairing201340707510510.1088/0954-3899/40/7/075105HoltJ. D.MenéndezJ.SimonisJ.SchwenkA.Three-nucleon forces and spectroscopy of neutron-rich calcium isotopes201490202431210.1103/physrevc.90.024312HagenG.Hjorth-JensenM.JansenG. R.MacHleidtR.PapenbrockT.Evolution of shell structure in neutron-rich calcium isotopes2012109303250210.1103/PhysRevLett.109.0325022-s2.0-84863892602RothR.BinderS.VobigK.CalciA.LanghammerJ.NavrátilP.Medium-mass nuclei with normal-ordered chiral NN+3N interactions2012109505250110.1103/physrevlett.109.052501SomàV.BarbieriC.DuguetT.Ab initio Gorkov-Green's function calculations of open-shell nuclei2013871501130310.1103/PhysRevC.87.011303SomàV.CipolloneA.BarbieriC.NavrátilP.DuguetT.Chiral two-and three-nucleon forces along medium-mass isotope chains201489606130110.1103/physrevc.89.061301TsukiyamaK.BognerS. K.SchwenkA.In-medium similarity renormalization group for nuclei20111062222250210.1103/physrevlett.106.222502HergertH.BognerS. K.BinderS.CalciA.LanghammerJ.RothR.SchwenkA.In-medium similarity renormalization group with chiral two-plus three-nucleon interactions201387303430710.1103/physrevc.87.034307EpelbaumE.HammerH.-W.MeißnerU.-G.Modern theory of nuclear forces20098141773182510.1103/revmodphys.81.17732-s2.0-75249090464HammerH.-W.NoggaA.SchwenkA.Colloquium: three-body forces: from cold atoms to nuclei201385119710.1103/revmodphys.85.197StrobergS. R.GadeA.TostevinJ. A.Neutron single-particle strength in silicon isotopes: constraining the driving forces of shell evolution201591404130210.1103/PhysRevC.91.041302DikmenE.ÖztürkO.CengizY.Fitting single particle energies in sdgh major shell201563222222810.1088/0253-6102/63/2/152-s2.0-84923204996UtsunoY.OtsukaT.TsunodaY.Recent advances in shell evolution with shell-model calculations20156010007EgidoJ. L.BorrajoM.RodŕiguezT.Collective and single-particle motion in beyond mean field approaches201611605250210.1103/PhysRevLett.116.052502Santiago-GonzalezD.WiedenhöverI.AbramkinaV.Triple configuration coexistence in 44S20118306130510.1103/physrevc.83.061305UtsunoY.ShimizuN.OtsukaT.YoshidaT.TsunodaY.Nature of isomerism in exotic sulfur isotopes2015114303250110.1103/physrevlett.114.032501NowackiF.PovesA.New effective interaction for 0ℏω shell-model calculations in the sd-pf valence space200979101431010.1103/physrevc.79.014310TsunodaN.OtsukaT.ShimizuN.Hjorth-JensenM.TakayanagiK.SuzukiT.Exotic neutron-rich medium-mass nuclei with realistic nuclear forceshttps://arxiv.org/abs/1601.06442CoraggioL.ItacoN.CovelloA.GarganoA.KuoT. T. S.Ground-state properties of closed-shell nuclei with low-momentum realistic interactions200368303432010.1103/PhysRevC.68.0343202-s2.0-0344927980GnezdilovN. V.BorzovI. N.SapersteinE. E.TolokonnikovS. V.Self-consistent description of single-particle levels of magic nuclei201489303430410.1103/physrevc.89.0343042-s2.0-84898796931NavrátilP.VaryJ. P.BarrettB. R.Properties of 12C in the Ab initio nuclear shell model200084255728573110.1103/physrevlett.84.5728Arias de SaavedraF.CóG.FabrociniA.Correlated model for Λ hypernuclei200163606430810.1103/physrevc.63.064308PieperS. C.WiringaR. B.CarlsonJ.Quantum Monte Carlo calculations of excited states in A=6–8 nuclei20047051105432510.1103/PhysRevC.70.054325DeanD. J.Hjorth-JensenM.Coupled-cluster approach to nuclear physics200469505432010.1103/physrevc.69.0543202-s2.0-15844390031CaurierE.Martínez-PinedoG.NowackiF.PovesA.ZukerP. A.The shell model as a unified view of nuclear structure200577242748810.1103/RevModPhys.77.4272-s2.0-23944445856EpelbaumE.Few-nucleon forces and systems in chiral effective field theory200657265474110.1016/j.ppnp.2005.09.002HagenG.PapenbrockT.DeanD. J.Hjorth-JensenM.Medium-mass nuclei from chiral nucleon-nucleon interactions2008101909250210.1103/PhysRevLett.101.0925022-s2.0-50849095226FujiiS.OkamotoR.SuzukiK.Publisher's note: ground-state and single-particle energies of nuclei around 16O, 40Ca, and 56Ni from realistic nucleon-nucleon forces [Phys. Rev. Lett. 103, 182501 (2009)]20091031918250110.1103/physrevlett.103.199903OtsukaT.SuzukiT.HoltJ. D.SchwenkA.AkaishiY.Three-body forces and the limit of oxygen isotopes20101053403250110.1103/PhysRevLett.105.032501HergertH.BinderS.CalciA.LanghammerJ.RothR.Ab Initio calculations of even oxygen isotopes with chiral two-plus-three-nucleon interactions20131102424250110.1103/physrevlett.110.242501EpelbaumE.KrebsH.LähdeT. A.LeeD.MeißnerG.RupakG.Ab Initio calculation of the spectrum and structure of 16O20141121010250110.1103/physrevlett.112.102501CipolloneA.BarbieriC.NavrátilP.Chiral three-nucleon forces and the evolution of correlations along the oxygen isotopic chain201592101430610.1103/PhysRevC.92.014306KellyJ. J.NegeleJ. W.VogtE.Nucleon knockout by intermediate energy electrons199623New York, NY, USASpringer75294DickhoffW. H.BarbieriC.Self-consistent Green's function method for nuclei and nuclear matter200452237749610.1016/j.ppnp.2004.02.038BenmokhtarF.Measurement of the 3Hee,ep´pn reaction at high missing energies and momenta200594808230510.1103/PhysRevLett.94.082305EgiyanK. S.DashyanN. B.SargsianM. M.Measurement of two- and three-nucleon short-range correlation probabilities in nuclei200696808250110.1103/physrevlett.96.082501PiasetzkyE.SargsianM.FrankfurtL.StrikmanM.WatsonJ. W.Evidence for strong dominance of proton-neutron correlations in nuclei20069716416250410.1103/physrevlett.97.1625042-s2.0-33750179982ShneorR.MonaghanP.SubediR.Investigation of proton-proton short-range correlations via the 12Ce,e´pppn200799707250110.1103/physrevlett.99.072501HenO.SargsianM.WeinsteinL. B.Momentum sharing in imbalanced Fermi systems2014346620961461710.1126/science.1256785AntonovA. N.HodgsonP. E.PetkovI. Zh.1988Oxford, UKClarendon PressPandharipandeV. R.SickI.DeWitt HubertsP. K. A.Independent particle motion and correlations in fermion systems199769398199210.1103/RevModPhys.69.9812-s2.0-0031480081BishopR. F.HowesC.IrvineJ. M.ModarresM.Constrained variational calculations for finite nuclei19784111709172310.1088/0305-4616/4/11/0052-s2.0-0009970183RistingM. L.ClarkJ. W.Density matrix of quantum fluids1976142875288710.1103/physrevb.14.2875ModarresM.RajabiA.MoshfeghH. R.The one-body momentum distribution of nuclear matter at finite temperature20088081–4607210.1016/j.nuclphysa.2008.05.0132-s2.0-45849092943ModarresM.RajabiA.The single-particle potential of nuclear matter in the LOCV framework2011867111110.1016/j.nuclphysa.2011.07.0092-s2.0-80052263867ModarresM.RajabiA.MoshfeghH. R.State-dependent calculation of three-body cluster energy for nuclear matter and the validity of the lowest order constrained variational formalism200776606431110.1103/physrevc.76.064311ModarresM.IrvineJ. M.LOCV calculations with a self-consistent treatment of isobars19795451152610.1088/0305-4616/5/4/016MoshfeghH. R.ModarresM.The properties of hot nuclear matter in LOCV formalism20057491–413013310.1016/j.nuclphysa.2004.12.0212-s2.0-13444251536ModarresM.BordbarG. H.Incompressibility of hot asymmetrical nuclear matter: lowest order constrained variational approach19985852781278610.1103/physrevc.58.27812-s2.0-0032327117ModarresM.Lowest-order constrained variational calculation for hot asymmetric nuclear matter199723892393710.1088/0954-3899/23/8/0072-s2.0-0031526514ModarresM.MoshfeghH. R.Lowest order constrained variational calculation for nuclear and neutron matter with a new charge-dependent reid potential20041121213610.1143/ptp.112.21MoshfeghH. R.ModarresM.The effect of three-body cluster energy on LOCV calculation for hot nuclear and neutron matter199824482182910.1088/0954-3899/24/4/0122-s2.0-0040120223GoudarziS.MoshfeghH. R.Effects of three-body forces on the maximum mass of neutron stars in the lowest-order constrained variational formalism20159151005432010.1103/PhysRevC.91.054320BordbarG. H.BigdeliM.YazdizadehT.Structure of neutron star with a quark core20062128-295991600110.1142/S0217751X060338542-s2.0-33845346763BordbarG. H.ZebarjadS. M.ZahediniaR.Lowest order constrained variational calculation of structure properties of protoneutron star2009481617010.1007/s10773-008-9782-22-s2.0-60449115355ModarresM.Local density approximation for alpha-particle binding energy19841011L251L25510.1088/0305-4616/10/11/003ModarresM.MoshfeghH. R.MarijiH.Lowest order constrained variational and local density approximation approach to the hot alpha particle200280891192310.1139/p02-0412-s2.0-67049086341ModarresM.RasekhinejadN.Ground state of heavy closed shell nuclei: an effective interaction and local density approximation approach200572601430110.1103/physrevc.72.064306ModarresM.RasekhinejadN.Ground state of heavy closed shell nuclei: an effective interaction and local density approximation approach20057261006430610.1103/physrevc.72.064306ModarresM.MarijiH.RasekhinejadN.The ground state binding energy of the closed shell nuclei with the density dependent Av18 effective interaction in LOCV method2011312909204310.1088/1742-6596/312/9/0920432-s2.0-80455158433ModarresM.RasekhinejadN.MarijiH.The density-dependent Av18 effective interaction and ground state of closed shell nuclei201120367970310.1142/S02183013110181622-s2.0-79955823075ModarresM.MarijiH.RasekhinejadN.The effect of density dependent Av18 effective interaction on the ground state properties of heavy closed shell nuclei20118591162810.1016/j.nuclphysa.2011.04.0132-s2.0-79957872952ModarresM.MarijiH.Effect of nondiagonal lowest order constrained variational effective two-body matrix elements on the binding energy of closed shell nuclei201286505432410.1103/physrevc.86.054324MarijiH.Imposing Fermi momentum cut-off on the channel- and density-dependent effective interaction and the ground-state properties of closed shell nuclei2014503, article 5610.1140/epja/i2014-14056-3MarijiH.ModarresM.The effect of Fermi momentum cutoff on the binding energy of closed shell nuclei in the LOCV framework201411324525110.1134/S15474771140301332-s2.0-84900796870MarijiH.The effect of gap in n(k,ρ) on the single-particle properties of nucleons and the ground-state binding energy of closed-shell nuclei2016524, article 10911010.1140/epja/i2016-16109-yMahauxC.SartorR.1991New York, NY, USASpringer US; Plenum PressEdited by J. W. Negele and E. VogtFroisP. B.FroisP. B.SickI.Nucleon distributions and the nuclear many-body problem1991London, UKWorld ScientificWiringaR. B.StoksV. G. J.SchiavillaR.Accurate nucleon-nucleon potential with charge-independence breaking1995511385110.1103/physrevc.51.382-s2.0-35949007980BohrN.Mottelson19691SchwierzN.WiedenhöverI.VolyaA.Parameterization of the woods-saxon potential for shell-model calculationshttps://arxiv.org/abs/0709.3525KoepfW.RingP.The spin-orbit field in superdeformed nuclei: a relativistic investigation19913391819010.1007/bf012829362-s2.0-0002471216MairleG.GrabmayrP.A phenomenological spin-orbit and nuclear potential for 208Pb20009331332610.1007/s1005000700162-s2.0-0034361514GraweH.BlazhevA.GórskaM.GrzywaczR.MachH.MukhaI.Nuclear structure far off stability—implications for nuclear astrophysics200627125726710.1140/epja/i2006-08-040-72-s2.0-33646421644PieperS. C.PandharipandeV. R.Origins of spin-orbit splitting in 15N19937017254110.1103/physrevlett.70.2541ThompsonI. J.NunesF. M.2009London, UKCambridge University Press