Effect of Incompressibility and Symmetry Energy Density on Charge Distribution and Radii of Closed-Shell Nuclei | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Kirkuk Journal of Science | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Article 2, Volume 17, Issue 3, September 2022, Pages 17-28 PDF (949.75 K) | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Document Type: Research Paper | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
DOI: 10.32894/kujss.2022.135889.1073 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Authors | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Shaymaa H. Amin* 1; Ahmed Aziz Al-Rubaiee2; Ali H. Taqi3 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Department of Physics, College of Science, Al-Mustansiriyah University, Baghdad, Iraq | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2Department of Physics, College of Science, Al-Mustansiriyah University, Baghdad, Iraq. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3Department of Physics, College of Science, Kirkuk University, Kirkuk, Iraq. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Abstract | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
In this work, the effect of incompressibility modulus KNM and symmetry energy density J on charge distribution and root-mean-square radii of neutron Rn and proton Rp has been investigated for light, medium and heavy closed-shell nuclei 40, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb within the framework of self-consistent Hartree-Fock (HF) with 20 types of Skyrme interactions having wide range of nuclear properties such as incompressibility modulus KNM and symmetry energy density J. The nuclear charge densities have been obtained and compared with the experiment data to give us a picture about the internal structure of the investigated nuclei. Also, the sensitivity of proton and neutron root-mean-square radii to nuclear matter properties has been examined by determining and discussion the statistical Pearson linear correlation coefficient. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Keywords | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Hartree-Fock HF; Skyrme interaction; Charge distribution; Root-mean-square radii | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Full Text | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1. Introduction: The main problem in physics is many-body problem [1], In nuclear physics, the many-body systems and the nature of the nuclear force [2] are the most important challenges facing physics the theoretical nuclear physics, therefore, theoretical researches in nuclear physics aim to describe the structure of the nuclei with a universal nuclear model. The properties of ground-state of atomic nuclei can be studied by Hartree-Fock (HF) [3, 3] or Hartree-Fock-Bogoliubov [5, 6]. There are many models used to study the excited levels, and the Random Phase Approximation (RPA) is considered successful in studying the closed-shell nuclei [7]. Charge density distribution and radii are fundamental nuclear properties applied to probe the nucleus structure and describes the effect of effective interactions on nuclear structure [8, 9] which were measured by electron elastic scattering and muonic X-ray data [10] by methods based on the electromagnetic interaction between the nucleus and electrons or muons. It has also been indicated by many research groups that the nuclear properties are closely related to the symmetry term of the equation of state (EOS) [11, 12]. Charge radii is a good indicator of nuclear properties and the nuclear symmetry energy that characterizes the equation of state of matter [13], also, it found correlates strongly with other observables in nuclei, such as nuclear charge and radii [14-16], or neutron skin thickness [17-19], and their relationship to the nuclear matter properties [20, 21]. The present work aims to study the effect of incompressibility modulus KNM and symmetry energy density J on charge distribution and root-mean-square radii of neutron (Rn)and proton (Rp) for light, medium and heavy closed-shell nuclei 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb using the self-consistent Hartree-Fock (HF) with 20 types of Skyrme interactions. Having a large number of Skyrme-force parameterizations requires continuous search for the best to describe the experimental data, therefore, the nuclear charge densities to be obtained and compared with the experiment data to give us an imagination about the internal structure of the investigated nuclei, as well as, the dependence of proton and neutron root-mean-square radii to nuclear matter properties such as incompressibility modulus KNM and symmetry energy density J to be examined by determining the statistical Pearson linear correlation coefficient.
2. Theoretical Framework: The second quantization Hamiltonian of a many-body fermion system is a summation of the kinetic energy t is and the anti-symmetrized two body matrix-elements as given below [22],
Where are annihilation and creation operators, respectively with The indices: covers all degrees of freedom of all available single particle states. The Hartree-Fock energy can be expressed as,
where single-particle density is diagonal which has the eigenvalues 1 or 0 for occupied and non-occupied states, and the energy should be minimized to determine the HF-basis [22],
with the self-consistent field.
The Skyrme effective interaction [23, 24] is one of the most convenient forces used in HF calculations to describe the nuclear probertite in ground-state, commonly contains central, non-local and density-dependent terms [25, 26]:
where , is the spin-exchange operator and is the Pauli spin operator. The Skyrme energy-density functional contains many terms and given by [27],
Where is a kinetic-energyterm, a zero-range term, is a density-dependent term, is an effective-mass term, is a finite-range term, is a spin-orbit term, is a tensor coupling with spin and gradient term and is a Coulomb interaction term. The charge density distribution in ground-state (in terms of Skyrme HF radial wave function of the state and occupation probability ) can be obtained by [28, 29],
The root-mean-square (rms) radii is defined as,
3- Result and Discussion: In this work, The HF calculations have been investigated for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei by solving the equations numerically with 20 Skyrme-type interaction: SkM* [30], SQMC650 [31], SV-K218 [32], SQMC700 [31], SkO [33], NRAPR [34], KDE0 [35], T31 [36], Skχm∗ [37], v075 [38], SkSC10 [39], RATP [40], SAMi [41], SK255 [42], SkD [43], GS2 [44], SIV [45], QMC2 [46], SIII [44], SVII [47]) the meaning of the abbreviation for each type of interaction and its parameters can be found in the cited references). Nuclear matter properties of the used Skyrme interaction such as the incompressibility modulus KNM (MeV) and the symmetry energy density J (MeV), are presented in Table 1. Our results are consistent with those of HFB calculations [5. 6]. The nuclear charge densities have been obtained and compared with the experiment data to give us a picture about the internal structure of the investigated. Figs 1-4 show the calculated HF charge distribution for 40Ca and 48Ca. Good agreements were obtained between our HF results for all of the used types of Skyrme interaction with the experimental data [8, 49] at the surface and interior regions, while our results underestimated the experimental data line inside the investigated nuclei SVII, SIII, QMC2 and GS2 sets. For 90Zr, our results of charge distribution for all interactions were in a good agreement with the experimental data [8, 50] in the center and at the surface of the nucleus but not so good inside the nucleus as GS2, QMC2 and SVII interactions underestimated the experimental data line inside the nucleus and SQMC650, SQMC700, NRAPR, Skχm∗, SK255, T31, SAMi, RATP parameters in less overestimated the experimental data line inside the nucleus. There are many factors that control the success of the theoretical model, in Table 1, we notice that the value of KMN for the first group is higher than that of the second group. For 116Sn, our results of charge distribution for all interactions were in a good agreement with the experimental data [8] in the center and at the surface of the nucleus but not so good inside the nucleus as in SQMC650, SKO, Skχm∗ parameters in less overestimated the experimental data line inside the nucleus and GS2, QMC2, SIII and SVII interactions underestimated the experimental data line inside the nucleus. The differences in describing data may due to differences in nuclear properties such as KNM, where its value the first group is lower than of the second group. For 144Sm our results of charge distribution for all interactions were in a good agreement with the experimental data [8] at the surface of the nucleus but not so good inside the nucleus. For 208Pb our results of charge distribution for all interactions were in a good agreement with the experimental data [8] in the center and at the surface of the nucleus but not so good inside the nucleus as GS2, QMC2 and SVII and less for SIV and SIII interactions underestimated the experimental data line inside the nucleus and SQMC650, KDE0 parameters in less overestimated the experimental data line inside the nucleus. This may be due to differences in the nuclear matter properties of the investigated types of Skyrme interaction.
Table 1: Nuclear matter properties of the used Skyrme interactions [48].
Fig.1 Our obtained charge density distribution (fm-3) of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei with Skyrme interactions Skm*, SQMC650, SV-K218, SQMC700, SKO. Charge densities are plotted by black solid line. The red dotted line shows experimental data taken from Ref. [8].
Fig. 2 Our obtained charge density distribution (fm-3) of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei with Skyrme interactions NRAPER, KDE0, T31 Skχm∗ and v075. Charge densities are plotted by black solid line. The red dotted line shows experimental data taken from Ref. [8].
Fig. 3 Our obtained charge density distribution (fm-3) of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei with Skyrme interactions SKSC10, RATP, SAMi, SK255 and SKD. Charge densities are plotted by black solid line. The red dotted line shows experimental data taken from Ref. [8].
Fig. 4 Our obtained charge density distribution (fm-3) of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei with Skyrme interactions GS2, SIV, QMC2, SIII and SVII. Charge densities are plotted by black solid line. The red dotted line shows experimental data taken from Ref. [8].
By considering a light, medium and heavy nuclei, we can illustrate the impact of the effect of incompressibility modulus KNM and symmetry energy density J on root-mean-square radii of neutron Rn and proton Rp for light, medium and heavy closed-shell nuclei 40, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb. In Fig. 5. the incompressibility KNM are shown as function of Rn and Rp for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei respectively. the linear correlation coefficients of KNM versus Rn, within the Skyrme forces were 0.395, 0.369, 0.512*, 0.515*, 0.605** and 0.544* for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei, respectively, where it is small (weak) in 40Ca and 48Ca. The linear correlation coefficients of KNM versus Rp within the Skyrme forces reaches 0.378, 0.489*, 0.552*, 0.562*, 0.591* and 0.606** for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei, respectively, also it is small (weak) in 40Ca and 48Ca. From Fig. 5, we conclude there is an impact of the incompressibility KNM on Rn and Rp. In Fig. 6, the symmetry energy density J are shown as function of Rn and Rp for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei respectively within the Skyrme forces, the linear correlation coefficients of J versus Rn were -0.433-, -0.179-, -0.308-, -0.254-, -0.212- and -0.032- for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei, respectively and it is very small correlated except for 40Ca it is rather better than other. The linear correlation coefficients of J versus Rp were -0.462*-, -0.497*-, -0.496*-, -0.499*-, -0.475*- and -0.476*- for 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei, respectively are small (weak). From Fig. 6, It can be concluded that the impact of the symmetry energy density J on neutron radii (J ↔ Rn) is weaker than that of Rp (J ↔ Rp) or nonexistent. In general, from Figs. 5 and 6., there is an influence of two nuclear matter properties (J and KNM) on proton and neutron radii in a space of Skyrme functional, we demonstrate the existence of a relation between (KNM↔Rn andRp) for all the investigated nuclei except for 40Ca and 48Ca (KNM↔Rn) it is weak. This illustrates that these quantities are coupled by the skyrme force and higher than (J↔Rn , Rp), where the relation (J↔Rn) is smaller (weaker) than (J↔Rp).
Fig. 5. Incompressibility KNM of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei as a function of (a) Rn and (b) Rp.
Fig. 6. Symmetry energy J of 40Ca, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb nuclei as a function of (a) Rn and (b) Rp.
4. Conclusions: In this study, the effect of incompressibility modulus KNM and symmetry energy density J on charge distribution and root-mean-square radii of neutron Rn and proton Rp has been investigated for light, medium and heavy closed-shell nuclei 40, 48Ca, 90Zr, 116Sn, 144Sm and 208Pb within the framework of self-consistent Hartree-Fock (HF) with 20 types of Skyrme interaction. Good agreements were obtained between our HF results of charge density distribution for all of the used types of Skyrme interaction with the experimental data at the surface and interior regions, while our results underestimated the experimental data line inside the investigated nuclei for SVII, SIII, QMC2 and GS2 sets. Concerning, the sensitivity of proton and neutron root-mean-square radii to nuclear matter properties, we conclude there is an impact of the incompressibility KNM on Rn and Rp. Also, it can be concluded that the impact of the symmetry energy density J on neutron radii (J ↔ Rn) is weaker than that of Rp (J ↔ Rp) or nonexistent.
Funding: None. Data Availability Statement: All of the data supporting the findings of the presented study are available from corresponding author on request. Declarations: Conflict of interest: The authors declare that they have no conflict of interest. Ethical approval: The manuscript has not been published or submitted to another journal, nor is it under review. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
References | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
[1] P Ring and P Schuck "The Nuclear Many-Body Problem" (New York: Springer-Verlag) (1980).
[2] K. S. Krane “Introductory Nuclear Physics” 2nd Ed., Wiley, United State of America (1988).
[3] Ali H. Taqi and E. G. Khidher, “Nuclear multipole excitations in the framework of self-consistent Hartree–Fock random phase approximation for Skyrme forces”, Pramana-J. Phys. 93, 60 (2019).
[4] Ali H. Taqi and G. I. Alawi, “Isoscalar giant resonance in 100, 116, 132Sn isotopes using Skyrme HF-RPA”, Nuc. Phy. A983, 103 (2019).
[5] Ali H. Taqi, Malik A. Hasan, “Ground‑State Properties of Even–Even Nuclei from He (Z = 2) to Ds (Z = 110) in the Framework of Skyrme–Hartree–Fock–Bogoliubov Theory”, Arabian Journal for Science and Engineering”, 47, pages761–775 (2022).
[6] Ali H. Taqi and Safa M. Qatal, “Nuclear Structure of Samarium Isotopes Using Skyrme and Gogny Hartree–Fock–Bogoliubov Method”, Iran J Sci Technol Trans A Sci, 46, 967(2020).
[7] Ali H. Taqi and Abdullah H. Ibrahim, “Collective Excitations of 14N and 10B Nuclei”, Kirkuk University Journal-Scientific Studies, 11(1), 202 (2016).
[8] I. Angeli, K.P. Marinova, "Table of experimental nuclear ground state charge radii: an update", At. Data Nucl. Data Tables 99(1), 69 (2013).
[9] Junjie Yang, J. Piekarewicz, "Difference in proton radii of mirror nuclei as a possible surrogate for the neutron skin", Phys. Rev. C97, 014314 (2018).
[10] H. De Vries, C. W. de Jager, and C. De Vries, "Nuclear charge-density-distribution parameters from elastic electron scattering", At. Data Nucl. Data Tables 36, 495 (1987).
[11] S. Typel and B. A. Brown, " Neutron radii and the neutron equation of state in relativistic models", Phys. Rev. C64, 027302 (2001).
[12] P. Danielewicz, "Surface symmetry energy", Nucl. Phys. A727, 233 (2003).
[13] M. Kortelainen, J. Erler, W. Nazarewicz, N. Birge, Y. Gao, and E. Olsen, " Neutron-skin uncertainties of Skyrme energy density functionals", Phys. Rev. C 88, 031305(R) (2013).
[14] P.-G. Reinhard and W. Nazarewicz, " Nuclear charge and neutron radii and nuclear matter: Trend analysis in Skyrme density-functional-theory approach", Phys. Rev. C 93, 051303(R) (2016).
[15] B. A. Brown, " Neutron Radii in Nuclei and the Neutron Equation of State", Phys. Rev. Lett. 85, 5296 (2000).
[16] R. Furnstahl, " Neutron radii in mean-field models", Nucl. Phys. A 706, 85 (2002).
[17] P.-G. Reinhard and W. Nazarewicz, "Information content of a new observable: The case of the nuclear neutron skin", Phys. Rev. C 81, 051303 (2010).
[18] F. J. Fattoyev and J. Piekarewicz, " Has a thick neutron skin in 208Pb been ruled out?", Phys. Rev. Lett. 111, 162501 (2013).
[19] J. Piekarewicz, B. K. Agrawal, G. Colò, W. Nazarewicz, N. Paar, P.-G. Reinhard, X. Roca-Maza, and D. Vretenar, " Electric dipole polarizability and the neutron skin", Phys. Rev. C85, 041302 (2012).
[20] X. Roca-Maza, M. Centelles, X. Viñas, and M. Warda, " Neutron Skin of 208Pb, Nuclear Symmetry Energy, and the Parity Radius Experiment", Phys. Rev. Lett. 106, 252501 (2011).
[21] W. Nazarewicz, P.-G. Reinhard, W. Satuła, and D. Vretenar, "Symmetry energy in nuclear density functional theory", EPJ A volume 50, Article number: 20 (2014).
[22] Ali H. Taqi and M. S. Ali, “Self-consistent Hartree-Fock RPA calculations in 208Pb” Indian J. Phys. 92, 69 (2017).
[23] T. H. R. Skyrme, "The effective nuclear potential", Nucl. Phys., 9, 615 (1959).
[24] D. Vautherin and D. M. Brink, " Hartree-Fock Calculations with Skyrme's Interaction. I. Spherical Nuclei", Phys. Rev., C5 (1972) 626 (1959).
[25] M. Bender, P.-H. Heenen, and P.-G. Reinhard, "Self-consistent mean-field models for nuclear structure", Rev. Mod. Phys.75, 121 (2003).
[26] J. R. Stone and P. -G. Reinhard, "The Skyrme interaction in finite nuclei and nuclear matter", Prog. Part. Nucl. Phys.58(2), 587 (2007).
[27] Ali H. Taqi, Zainab Q. Mosa, “Nuclear Structure of Zirconium Isotopes 90,96,106Zr Using Skyrme HF-RPA”, Nucl. Egypt. J. Phy, 48, 1 (2020).
[28] Ali H. Taqi, Ebtihal G. Khider, “Ground and transition properties of 40Ca and 48Ca nuclei”, Nucl. Phys. At. Energy, 19, 326 (2018).
[29] Ali H. Taqi. A visual Fortran 90 program for the two-particle or two-hole excitations of nuclei: The PPRPA program. SoftwareX 5, 51 (2016).
[30] J. Bartel, P. Quentin, M. Brack, C. Guet and H.B. Hakansson, "Toward a Better Parametrisation of Skyrme-Like Effective Forces: A Critical Study of the SKM* Force", Nucl. Phys. A386, 79 (1982).
[31] S. Shlomo, V. M. Kolomietz, and G. Colo, " Deducing the nuclear-matter incompressibility coefficient from data on isoscalar compression modes", Eur. Phys. J. A30, 23 (2006).
[32] P. Klupfel, P.-G. Reinhard, T. J. Burvenich, and J. A. Maruhn, " Variations on a theme by Skyrme: A systematic study of adjustments of model parameters", Phys. Rev. C79, 034310 (2009).
[33] P.-G. Reinhard, D. J. Dean, W. Nazarewicz, J. Dobaczewski, J. A. Maruhn, and M. R. Strayer, "Shape coexistence and the effective nucleon-nucleon interaction", Phys. Rev. C 60, 014316 (1999).
[34] A.W. Steiner, M. Prakash, J. M. Lattimer, and P. J. Ellis, " Isospin asymmetry in nuclei and neutron stars", Phys. Rep. 411, 325 (2005).
[35] M. N. Harakeh and A. van der Woude, “Giant Resonances”, Oxford Studies on Nuclear Physics. Oxford Science Publications, 2001.
[36] T. Lesinski, M. Bender, K. Bennaceur, T. Duguet, and J. Meyer, " Tensor part of the Skyrme energy density functional: Spherical nuclei", Phys. Rev. C76, 014312 (2007).
[37] Zhen Zhang, Yeunhwan Lim, Jeremy W. Holt, Che Ming Ko, " Nuclear dipole polarizability from mean-field modeling constrained by chiral effective field theory", Phys Lett. B777, 73 (2018).
[38] J. M. Pearson and S. Goriely, "Isovector effective mass in the Skyrme-Hartree-Fock method", Phys. Rev. C64, 027301 (2001).
[39] M. Onsi, H. Przysiezniak, and J. M. Pearson, " Equation of state of homogeneous nuclear matter and the symmetry coefficient", Phys. Rev. C 50, 460 (1994).
[40] M. Rayet, M. Amould, E Tondeur and G. Paulus, "Nuclear force and the properties of matter at high temperature and density", Astron. Astrophys. 116, 183 (1982).
[41] X. Roca-Maza, G. Colo, and H. Sagawa, " New Skyrme interaction with improved spin-isospin properties", Phys. Rev. C 86, 031306 (2012).
[42] B. K. Argawal, S. Shlomo and V. Kim Au, "Nuclear matter incompressibility coefficient in relativistic and nonrelativistic microscopic models", Phys. Rev. C68 (2003), 031304 (2003).
[43] Yongli Xu, Hairui Guo, Yinlu Han, and Qingbiao Shen, “New Skyrme interaction parameters for a unified description of the nuclear properties”, J. Phys. G: Nucl. Part. Phys. 41, 015101 (2014)
[44] Q. B. Shen, Y. L. Han, and H. R. Guo, “Isospin dependent nucleon-nucleus optical potential with Skyrme interactions”, Phys. Rev. C80, 024604 (2009).
[45] B. A. Brown, G. Shen, G. C. Hillhouse, J. Meng, and A. Trzcińska, “Neutron skin deduced from antiprotonic atom data”, Phys. Rev. C76, 034305 (2007).
[46] P. A. M. Guichon and A. W. Thomas, “Quark Structure and Nuclear Effective Forces”, Phys. Rev. Lett. 93, 13 (2004).
[47] M. J. Giannoni and P. Quentin, " The mass parameter in the generator coordinate method", Phys. Rev. C21, 2076 (1980).
[48] M. Dutra, O. Lourenco, J. S. Sa Martins, and A. Delfino, “Skyrme interaction and nuclear matter constraints”, Phys. Rev. C85, 035201 (2012).
[49] J. Zenihiro, H. Sakaguchi, S. Terashima, T. Uesaka, G. Hagen, M. Itoh, T. Murakami, Y. Nakatsugawa, T. Ohnishi, H. Sagawa, H. Takeda, M. Uchida, H.P. Yoshida, S. Yoshida, M. Yosoi, "Direct determination of the neutron skin thicknesses in 40,48Ca from proton elastic scattering at Ep=295 MeV", arXiv:1810.11796 (2018).
[50] H. de Vries, C.W. de Jager, C. de Vries, "Nuclear Charge-Density-Distribution Parameters from Elastic Electron Scattering ", At. Data. Nucl. Data Tables 36, 495 (1987). | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Statistics Article View: 313 PDF Download: 172 |