Advances in High Energy Physics

Volume 2016 (2016), Article ID 2398198, 10 pages

http://dx.doi.org/10.1155/2016/2398198

## The Calculation of Single-Nucleon Energies of Nuclei by Considering Two-Body Effective Interaction, , and a Hartree-Fock* Inspired* Scheme

^{1}Centro de Física Computacional, Department of Physics, University of Coimbra, 3004-516 Coimbra, Portugal^{2}Centro de Física do Porto, Departamento de Física e Astronomia, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre 687, 4169-007 Porto, Portugal

Received 25 June 2016; Revised 2 September 2016; Accepted 20 September 2016

Academic Editor: Nasser Kalantar-Nayestanaki

Copyright © 2016 H. Mariji. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

The nucleon single-particle energies (SPEs) of the selected nuclei, that is, , , and , 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 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, , with the Heaviside functions, the role of 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 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 , 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 [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 . 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, , 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 , 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 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 . Their calculated have been in overall agreement with those of other methods. By applying on the calculation of single-nucleon properties, the authors obtained comparable results with other methods [57]. The model-dependent nature of , 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 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 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 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 in the calculation of density-dependent SPEs and impose a factor of nucleon effective mass, 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 , rather than using the usual Heaviside functions, and account for the presence of correlations in HF mean field by . Both and are generated by the LOCV calculations for SNM with 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 on the calculations at the minimum point of the total energy of a nucleus, denoted by ( 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 mean field potential on the calculation of spin-orbit interaction at . In Section 4, by including the Coulomb interactions, we improve the single-proton energies. In Section 5, we shall show the effect of on the calculation of SPEs in the ground-state of selected closed shell nuclei, that is, , , and . 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 on SPEs. On the other hand, we believe that including a factor of with 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 bywhere , the single-particle kinetic energy (SPKE) operator for* i*th nucleon, is given by and denotes the two-nucleon interaction potential. By denoting the occupied single-particle states by , the total energy of occupied states (the core energy), , in diagonal configuration, is given bywhere “” 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 as the state of specific nucleon and SPIE as and rewriting interaction part of the Hamiltonian, the total energy of core gets the following form:Now, we can present an expression to calculate the SPE of a nucleon in the core as follows:The total energy of a closed shell nucleus with configuration plus one nucleon (one absent-nucleon, so-called, one hole) in the state above (below) the Fermi surface is as follows:By defining , we can rewrite the total energy as follows:The 1*p*-1*h* SPE (near the Fermi surface) readsIn this work, we set the SHO basis; that is, in which with , , and 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 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, , by using the LOCV calculations for SNM with potential. Then we employ the elements in order to construct SPIE. Next, while we keep the quantum numbers of one orbit, , 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 and calculate the nucleon SPEs of the nucleus at . 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 with Heaviside functions.

According to (3) and above explanations, the ground-state energy of the core readswhere 0020and stand for SPKE and SPIE, respectively. Thus, the SPE of nucleon in a specific single-particle state of a closed shell nucleus at is given byand the 1*p*-1*h* SPE is as follows:where In the SHO basis, one can obtain SPIE from (4) and (7) as follows:where means that we need a factor of only in the calculation of SPIE for (cf. (2)). Regarding the LOCV method, the two-body effective interaction has the following form [54]:where , the density- and channel-dependent two-body correlation functions, are operator type-like [73]. The operator for the potential has the following form [81]:where the forms of eighteen components of the two-body operator are as follows: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 (a) between and , two complete basis sets of plane waves together with spin and isospin parts, and (b) between and , two orthogonal bases in the configuration space, and (c) between and , two complete sets in which denotes , , , , and where , , , and are the relative orbital momentum, the total spin, isospin, and isospin projection of two interacting nucleons, respectively, and comes from . Thus, a straightforward calculation, by taking just diagonal matrix elements and considering orthogonality conditions, (12) is simplified bywhereIn (17), are the isospin Clebsch-Gordan coefficients and are the well-known spherical Bessel functions in which the relative momentum and distance of two interacting nucleons are obtained by the familiar relations and , respectively. In the above equations, we converted the summations over and to integrals by the well-known relation as follows:where are the one-body momentum distribution functions. In this work, , generated by using the LOCV calculations for SNM with potential based on the Ristig-Clark formalism [55], are density-dependent [78]. We suppose that the volume of box, , is obtained by where callsIn (16) and (17), we can replace by , 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 ’s are topped to the Fermi momentum, , with being degeneracy parameter of nucleons. Now, by use of (9), (10), (16), and (17), and recalling , we find SPE as follows:where we replaced with . In the above equation, readswhere is a normalization coefficient for each natural orbit and , the effective mass of nucleon, is calculated at the Fermi surface by the LOCV method with considering potential [78]. In this work, we calculate 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 . Table 1 shows for the single-nucleon major levels of , , and . In order to clarify the importance of including in calculations, the values of are considered in two different cases: (i) including the constructed 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 , as discussed in Section 6.