×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

Correlation between the charge radii difference in mirror partner nuclei and thesymmetry energy slope

Figures(4) / Tables(2)

Get Citation
Xiao-Rong Ma, Shuai Sun, Rong An and Li-Gang Cao. Correlation between the charge radii difference in mirror partner nuclei...[J]. Chinese Physics C. doi: 10.1088/1674-1137/ad47a8
Xiao-Rong Ma, Shuai Sun, Rong An and Li-Gang Cao. Correlation between the charge radii difference in mirror partner nuclei...[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ad47a8 shu
Milestone
Received: 2024-01-10
Article Metric

Article Views(1178)
PDF Downloads(33)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Correlation between the charge radii difference in mirror partner nuclei and thesymmetry energy slope

    Corresponding author: Rong An, rongan@nxu.edu.cn
    Corresponding author: Li-Gang Cao, caolg@bnu.edu.cn
  • 1. School of Physics, Ningxia University, Yinchuan 750021, China
  • 2. Key Laboratory of Beam Technology of Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China
  • 3. Key Laboratory of Beam Technology of Ministry of Education, Institute of Radiation Technology, Beijing Academy of Science and Technology, Beijing 100875, China

Abstract: A correlation between the charge radii difference of mirror partner nuclei $ \Delta{R_{\mathrm{ch}}} $ and the slope parameter L of symmetry energy has been built to ascertain the equation of state of isospin asymmetric nuclear matter. In this work, the influences of pairing correlations and isoscalar compression modulus on the $ \Delta{R_{\mathrm{ch}}} $ are systematically investigated based on the Skyrme energy density functional theory. The calculated results suggest that the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L is decreased by the surface pairing correlations. The slope parameter deduced from the difference of charge radii of mirror-pair nuclei 32Ar-32Si, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe falls into the range of $ L=42.57 $$ 50.64 $ MeV, that is, the rather soft equation of state of asymmetric nuclear matter. Besides, the range of the slope parameter can also be influenced by the effective forces classified by various isoscalar incompressibility coefficients.

    HTML

    I.   INTRODUCTION
    • Nuclear symmetry energy is generally employed to comprehend the underlying physical mechanism from terrestrial nuclei to dense astrophysical events [16]. So far, plenty of methods have been applied to constrain the nuclear symmetry energy, such as the extracted neutron skin thickness (NST) [712], properties of giant resonances [1317], and quantities of the heavy-ion collision [1823]. However, a unified density-dependence of symmetry energy is hardly determined owing to the uncertainties of the different theoretical models. Thus, alternative observables are urgently required for constraining the equation of state (EoS) of asymmetric nuclear matter. The difference in charge radii of mirror-pair nuclei has been proposed as a tentative probe to evaluate the slope parameter of symmetry energy L [24, 25]. Therefore, the precisely measured charge radii differences in mirror partner nuclei 32Ar-32Si, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe have been used to validate the range of the slope parameter L [2628].

      Charge radii differences in mirror pair nuclei are strongly associated with the neutron-skin thickness of neutron-rich nuclei and with the slope parameter L [24, 25]. As shown in Ref. [29], the calculations performed using relativistic energy density functionals (EDFs) further inspect the validity of this assumption. It points out that the difference in proton radii of mirror nuclei can be indeed employed to constrain the EoS of isospin asymmetry nuclear matter. This may provide an alternative approach to extract content information about the neutron skin thickness [30]. A recent study suggests that the mirror-difference in nuclear charge radii is proportional to the isospin asymmetry $ I=(N-Z)/A $, where N and Z represent the neutron and proton number of a nucleus, respectively, and $ A=N+Z $ is the corresponding mass number [31]. Here, the same scenario can also be observed between the difference in binding energy per nucleon of mirror nuclei and the Coulomb asymmetry. Furthermore, combing the accurately detected charge radii differences in mirror partner nuclei, the upper limit range of the symmetry energy slope has been determined to be $ L\approx100 $ MeV [32].

      This method requires more accurate charge radii data, which should be obtained through different nuclear models. With the improvements in the experimental technique, high-precision charge radii of exotic nuclei can be detected [33]. Thus, more charge radii data of nuclei far away from the β-stability line can be compiled [34, 35]. This provides a fundamental guideline for theoretical studies to access the charge radii of mirror partner nuclei. A reliable description of nuclear charge radii is influenced by various mechanisms, such as shape deformation [3639], pairing correlations [40, 4142], and shell closure effects [43, 44]. As is well known, pairing correlations play an important role in describing the ground state properties of finite nuclei [4547] or the quasiparticle resonant states [48]. In particular, it has been mentioned that the radius of neutron-rich nuclei can also be influenced by the pairing correlations [49].

      As demonstrated in Ref. [50], the pairing correlations should be taken into account in calculating the nuclear charge radii. As presented in Ref. [51], the density energy functionals consisting of Skyrme and covariant models are used to systematically investigate the correlations between the charge radii differences in mirror partner nuclei and slope parameter L. The calculated results suggest that the linear correlation between the mirror difference in charge radii and slope parameter L is further decreased when considering the pairing effects. In contrast, the correlation between the neutron-skin thickness of a neutron-rich nucleus and the proton radii difference of the corresponding mirror nuclei is mostly enhanced by pairing effects [30]. Meanwhile, it has also been mentioned that the compression modulus of symmetric nuclear matter has an influence on determining the slope parameter L [52].

      The density-dependence of nuclear symmetry energy, which is mostly determined by the slope parameter L, is associated with the isovector-sensitive indicators in the EoS of isospin asymmetry systems. In Ref. [53], it is stated that the correlation between the isoscalar incompressibility K and the NST is not clear. Further research demonstrated that the correlations between K and the isovector parameters are generally weaker in comparison with the strong correlation between NST and symmetry energy coefficients [54, 55]. Actually, the nuclear matter properties recognizing the isoscalar and isovector counterparts are mutually correlated with the effective model parameters. This means that the influence coming from the isoscalar sector is non-negligible in discussing the isovector properties [56, 57]. As suggested in Refs. [5860], the compression modulus of symmetric nuclear matter is sensitive to the density dependence of symmetry energy. Therefore, it is also essential to investigate the correlation between the difference in charge radii of mirror-partner nuclei and the slope parameter L under various incompressibility coefficients K. In this work, the influences of pairing correlations and isoscalar compression modulus on the charge radii differences of mirror partner nuclei 32Ar-32Si, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe are investigated using the spherical Skyrme EDFs. Furthermore, the range of the slope parameter L derived from the charge radii differences of mirror-pair nuclei is discussed.

      The contents of this paper are organized as follows. Sec. II summarizes the theoretical framework. In Sec. III, the numerical results and discussion are provided. Finally, a summary and outlook are presented in Sec. IV.

    II.   THEORETICAL FRAMEWORK
    • The sophisticated Skyrme-type EDF, which is expressed as an effective zero-range force between nucleons with density-dependent and momentum-dependent terms, has achieved remarkable success in describing various physical phenomena [48, 49, 6175]. In general, the effective force sets are determined by calibrating the properties of the finite nuclei and symmetric nuclear matter at saturation density $ \rho_{0} $($ \approx0.16\; \mathrm{fm}^{-3} $). Notably, it has been mentioned that some interactions can be used to characterize the bulk properties of finite nuclei and infinite nuclear matter (see, e.g., [76, 77] for details). In this work, the Skyrme-like effective interaction has been recalled as follows [7880]:

      $ \begin{aligned}[b] V(\mathbf{r}_{1},\mathbf{r}_{2})=\;&t_{0}(1+x_{0}\mathbf{P}_{\sigma})\delta(\mathbf{r})\\ &+\frac{1}{2}t_{1}(1+x_{1}\mathbf{P}_{\sigma})\left[\mathbf{P}'^{2}\delta(\mathbf{r})+\delta(\mathbf{r})\mathbf{P}^{2}\right]\\ &+t_{2}(1+x_{2}\mathbf{P}_{\sigma})\mathbf{P}'\cdot\delta(\mathbf{r})\mathbf{P}\\ &+\frac{1}{6}t_{3}(1+x_{3}\mathbf{P}_{\sigma})[\rho(\mathbf{R})]^{\alpha}\delta(\mathbf{r})\\ &+\mathrm{i}W_{0}\mathbf{\sigma}\cdot\left[\mathbf{P}'\times\delta(\mathbf{r})\mathbf{P}\right]. \end{aligned} $

      (1)

      Here, $ \mathbf{r}=\mathbf{r}_{1}-\mathbf{r}_{2} $ and $ \mathbf{R}=(\mathbf{r}_{1}+\mathbf{r}_{2})/2 $ are related to the positions of two nucleons $ \mathbf{r}_{1} $ and $ \mathbf{r}_{2} $; $ \mathbf{P}=(\nabla_{1}-\nabla_{2})/2\mathrm{i} $ is the relative momentum operator and $ \mathbf{P'} $ is its complex conjugate acting on the left; and $ \mathbf{P_{\sigma}}=(1+\vec{\sigma}_{1}\cdot\vec{\sigma}_{2})/2 $ is the spin exchange operator that controls the relative strength of the $ S=0 $ and $ S=1 $ channels for a given term in the two-body interactions, with $ \vec{\sigma}_{1(2)} $ being the Pauli matrices. The last term represents the spin-orbit force where $ \sigma=\vec{\sigma}_{1}+\vec{\sigma}_{2} $. The quantities α, $ t_{i} $, and $ x_{i} $ ($ i=0 $-3) represent the parameters of the Skyrme forces used in this work.

      Solving the Schrödinger-like equations by self-consistent iteration leads to the eigenenergies and wave functions of the constituent nucleons. The pairing correlations can be generally treated either by the BCS method or by the Bogoliubov transformation [45, 46, 81]. In this work, the spherical Skyrme Hartree-Fock-Bogoliubov (HFB) approach is utilized for all the calculations [82]. The density-dependent zero-range pairing force is employed as follows [83]:

      $ V_{\mathrm{pair}}(\mathbf{r}_{1},\mathbf{r}_{2})=V_{0}\left[1-\eta\left(\frac{\rho(\mathbf{r})}{\rho_{0}}\right)\right]\delta(\mathbf{r}_{1}-\mathbf{r}_{2}), $

      (2)

      where $ \rho(\mathbf{r}) $ is the baryon density in the coordinate space and $ \rho_{0}=0.16\; \mathrm{fm}^{-3} $ represents the nuclear saturation density. The value of η is taken as $ 0.0 $, $ 0.5 $, or $ 1.0 $ for volume-, mixed-, or surface-type pairing interactions, respectively. The quantity $ V_{0} $ is adjusted by calibrating the empirical pairing gaps with the three-point formula [47, 83], and the values of $ V_{0} $ are 188.2, 370.2, and 509.6 MeV fm3 for the corresponding volume-, mixed-, and surface-type pairing interactions, respectively.

      After achieving the convergence of the total binding energy, the root-mean square (rms) radii of the neutron and proton matter, and the neutron-skin thickness can be obtained naturally. The quantity of nuclear charge radius ($ R_{\mathrm{ch}} $) can be calculated as follows (in units of fm2):

      $ \begin{array}{*{20}{l}} R_{\mathrm{ch}}^{2}=\langle{r_{\mathrm{p}}^{2}}\rangle+0.64\; \mathrm{fm^{2}}. \end{array} $

      (3)

      The first term represents the charge density distributions of point-like protons and the second term is due to the finite size of protons [84]. For mirror pair nuclei, the difference in charge radii ($ \Delta{R_{\mathrm{ch}}} $) can be obtained through the formula mentioned above.

      The density dependence of symmetry energy can be expanded around the saturation density $ \rho_{0} $ as follows:

      $ \begin{equation} E_{\mathrm{sym}}(\rho)\approx{E_{\mathrm{sym}}}(\rho_{0})+\frac{L}{3}\left(\frac{\rho-\rho_{0}}{\rho_{0}}\right)+\frac{K_{\mathrm{sym}}}{18}\left(\frac{\rho-\rho_{0}}{\rho_{0}}\right)^{2}+\dots, \end{equation} $

      (4)

      where L and $ K_{\mathrm{sym}} $ are the slope and curvature of symmetry energy at nuclear saturation density $ \rho_{0} $, respectively. These two quantities are defined as

      $ L=3\rho_{0}\frac{\partial{E_{\mathrm{sym}}(\rho)}}{{{\partial\rho}}}\Large{|}_{\rho=\rho_{0}}, $

      (5)

      $ K_{\mathrm{sym}}=9\rho_{0}^{2}\frac{\partial^{2}E_{\mathrm{sym}}(\rho)}{{\partial\rho^{2}}}\Large{|}_{\rho=\rho_{0}}. $

      (6)

      Further, the isoscalar incompressibility K of symmetric nuclear matter is recalled as

      $ K=9 \rho_0^2 \frac{\partial^2E_0(\rho)}{\partial \rho^2}|_{\rho=\rho_0}.$

      (7)

      Here, the quantity $ E_{0}(\rho) $ is the energy per particle of symmetric nuclear matter [30]. The nuclear incompressibility can be extracted from the measurements of the isoscalar giant monopole resonance (ISGMR) in medium-heavy nuclei [83, 85]. In addition, it is also mentioned that the incompressibility of symmetric nuclear matter can be deduced from the α-decay properties [86]. Table 1 presents the parameters of the Skyrme forces used in this work. The corresponding saturation properties in infinite nuclear matter derived from the Skyrme forces, such as isoscalar incompressibility K, slope parameter L, and symmetry energy $ E_{\mathrm{sym}} $, are also listed explicitly. Both the slope parameter L and symmetry energy $ E_{\mathrm{sym}} $ at saturation density cover a wide range.

      K Sets $ t_{0} $ $ t_{1} $ $ t_{2} $ $ t_{3} $ $ x_{0} $ $ x_{1} $ $ x_{2} $ $ x_{3} $ $ W_{0} $ L $ E_{\mathrm{sym}} $
      $ K\approx230 $s3028$ -2461.903 $472.082$ -530.239 $13599.531.6830$ -0.3349 $$ -1.0 $2.6602124.766$ -11.2262 $28
      s3030$ -2477.013 $475.604$ -533.657 $13709.091.1408$ -0.3365 $$ -1.0 $1.8302117.48222.871530
      s3032$ -2491.849 $486.614$ -579.661 $13842.350.9962$ -0.2791 $$ -1.0 $1.5763128.02836.224632
      s3034$ -2503.455 $489.112$ -591.648 $13939.400.7308$ -0.2622 $$ -1.0 $1.1527126.84656.144234
      s3036$ -2513.951 $491.018$ -615.662 $14045.310.5457$ -0.2210 $$ -1.0 $0.8437128.52871.542836
      s3038$ -2524.594 $494.060$ -635.510 $14142.730.3550$ -0.1916 $$ -1.0 $0.5315129.20087.615538
      s3040$ -2531.688 $489.598$ -625.976 $14203.480.1249$ -0.1940 $$ -1.0 $0.1662122.045106.086240
      $ K\approx240 $s4028$ -2296.534 $515.586$ -336.464 $11786.331.7714$ -0.8491 $$ -1.0 $3.1703118.8033.977428
      s4030$ -2317.546 $532.898$ -355.019 $11907.071.2924$ -0.8449 $$ -1.0 $2.3824118.26334.073530
      s4032$ -2321.087 $520.235$ -491.504 $12164.331.2626$ -0.5755 $$ -1.0 $2.1916136.98934.428332
      s4034$ -2329.794 $510.932$ -423.708 $12179.540.8546$ -0.6759 $$ -1.0 $1.5447117.93262.588434
      s4036$ -2337.516 $507.834$ -424.568 $12253.220.7156$ -0.6664 $$ -1.0 $1.2882118.41475.667936
      s4038$ -2355.765 $523.051$ -465.145 $12389.160.3790$ -0.6244 $$ -1.0 $0.7299118.23898.652238
      s4040$ -2359.001 $513.325$ -475.563 $12466.460.2959$ -0.5827 $$ -1.0 $0.5558119.406108.174140
      $ K\approx250 $s5028$ -2157.179 $599.608$ -396.271 $10352.381.4017$ -0.9193 $$ -1.0 $2.9016138.87433.003728
      s5030$ -2155.199 $570.378$ -299.172 $10324.291.6366$ -1.0334 $$ -1.0 $3.2711131.21830.024830
      s5032$ -2172.510 $583.512$ -348.992 $10472.921.4649$ -0.9688 $$ -1.0 $2.9140141.58943.587132
      s5034$ -2179.610 $584.069$ -331.053 $10501.791.2593$ -0.9977 $$ -1.0 $2.5371141.09460.320234
      s5036$ -2195.402 $599.442$ -365.556 $10609.250.9683$ -0.9629 $$ -1.0 $2.0020145.63880.176236
      s5038$ -2205.358 $602.093$ -375.998 $10690.190.7284$ -0.9493 $$ -1.0 $1.5517145.12497.492538
      s5040$ -2214.249 $600.714$ -404.383 $10802.960.5278$ -0.9019 $$ -1.0 $1.1569146.358112.207940

      Table 1.  Parameters of the forces used in this work, where the parameter $ \alpha=1/6 $. Saturation properties of the Skyrme parametrization sets used in this work, such as the incompressibility coefficients K (MeV), slope parameter L (MeV), and symmetry energy $ E_{\mathrm{sym}} $ (MeV) at saturation density $ \rho_{0} $ (fm$ ^{-3} $), are also listed.

    III.   RESULTS AND DISCUSSION
    • The linear correlation between the charge radii differences in the mirror partner nuclei 32Ar-32Si, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe and the slope parameter of symmetry energy have been sequentially employed to constrain the EoS of isospin asymmetric nuclear matter [2628]. In Table 2, the experimental data of $ R_{\mathrm{ch}} $ and $ \Delta{R_{\mathrm{ch}}} $ for $ A = 32 $, $ 36 $, $ 38 $, and $ 54 $ mirror-pair nuclei are given. As mentioned above, the nuclear charge radius is influenced by the pairing correlations [4042]. Thus, the influence on the correlations between the charge radii difference in mirror partner nuclei and the slope parameter L should be systematically investigated.

      A $ R_{\mathrm{ch}} $/fm $ \Delta{R_{\mathrm{ch}}} $/fm
      32Ar3.3468(62)
      Si3.153(12) 0.194(14)
      36Ca3.4484(27)
      S3.2982(12)0.150(4)
      38Ca3.4652(17)
      Ar3.4022(15)0.063(3)
      54Ni3.7370(30)
      Fe3.6880(17)0.049(4)

      Table 2.  $ R_{\mathrm{ch}} $ and $ \Delta{R_{\mathrm{ch}}} $ database for the $ A = 32 $, $ 36 $, $ 38 $, and $ 54 $ mirror-pair nuclei. The parentheses on the values of charge radii and the difference of charge radii show the systematic uncertainties [2628, 34, 87].

      In Fig. 1, plots for the charge radii difference $ \Delta{R_{\mathrm{ch}}} $ in 54Ni-54Fe against symmetry energy slope L are shown for various pairing forces. For the convenience of discussion, the Pearson coefficient R, which is used to measure the degree of correlation between two variables, is employed in our study. Without considering the pairing correlations, the Pearson coefficient R between the $ \Delta{R_{\mathrm{ch}}} $ of mirror partner nuclei 54Ni-54Fe and L is determined to be $ 0.993 $. The same scenario can also be found in mixed-type pairing forces. For the volume-type pairing force, the Pearson coefficient R is almost equivalent to those for the no pairing and $ \eta=0.0 $ cases. As demonstrated in Ref. [51], the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L is decreased by the pairing effects. From Fig. 1, one can find that the surface-type pairing force actually decreases the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L in mirror pair nuclei 54Ni-54Fe, where the Pearson coefficient R falls to $ 0.882 $. Actually, more particles are scattered into the higher single levels under the surface pairing force. This leads to the rather poor correlation between the difference of charge radii of mirror partner nuclei and the slope parameter L.

      Figure 1.  (color online) Charge radii difference $ \Delta{R_{\mathrm{ch}}} $ of mirror partner nuclei 54Ni-54Fe as a function of slope parameter L is depicted under various pairing forces: no pairing (a), volume type $ \eta=0.0 $ (b), mixed type $ \eta=0.5 $ (c), and surface type $ \eta=1.0 $ (d) pairing forces. The Pearson coefficients R are presented in parentheses. The database is taken from Refs. [28, 34, 35], and the systematic error bands are represented in light blue.

      For further focusing on the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L, the charge radii differences $ \Delta{R_{\mathrm{ch}}} $ of mirror partner nuclei 32Ar-32Si, 36Ca-36S, and 38Ca-38Ar against slope parameter L are shown in Fig. 2. In this figure, clear linear correlations can be observed for these mirror-pair nuclei (32Ar-32Si, 36Ca-36S, and 38Ca-38Ar). For the volume-type pairing force, the Pearson coefficients show almost no change in the mirror partner nuclei 32Ar-32Si, 36Ca-36S, and 38Ca-38Ar, with R values of 0.926, 0.881, and 0.836, respectively. Considering the mixed-type pairing force, the Pearson coefficients R are slightly improved with respect to the no pairing case. Actually, the Pearson coefficients R are rather poor if the surface pairing force is taken into account. The corresponding figures are not shown here, but the Pearson coefficients R for the surface-type pairing force are $ 0.887 $, $ 0.791 $, and $ 0.743 $ for the mirror partner nuclei 32Ar-32Si, 36Ca-36S, and 38Ca-38Ar, respectively. This means that the pairing effect can decrease the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L if the pairing interactions are only addressed by the surface-type pairing force.

      Figure 2.  (color online) Charge radii differences $ \Delta{R_{\mathrm{ch}}} $ of mirror partner nuclei 32Ar-32Si, 36Ca-36S, and 38Ca-38Ar as a function of slope parameter L depicted under no pairing and mixed-type $ \eta=0.5 $ pairing forces. The Pearson coefficients R are presented in parentheses. The database is taken from Refs. [26, 27, 34, 35], and the systematic error bands are represented in light blue.

      As mentioned in Refs. [2428], charge radii differences in mirror partner nuclei provide an alternative approach to pin down the isovector components in the EoS of asymmetric nuclear matter. To analyze the relatively clear linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L, the mixed-type pairing force is employed in the following discussion. In Fig. 3, the deduced slope parameters L from the differences of charge radii of mirror-pair nuclei 54Ni-54Fe, 38Ca-38Ar, 36Ca-36S, and 32Ar-32Si are drawn by shaping the range covered in Fig. 2. Combing the experimental data, the extracted slope parameter falls into the range of $ L=42.57 $$ 50.64 $ MeV. From this figure, one can find that the deduced value of L covers the uncertainty range derived from the mirror charge radii differences [24, 2628, 32, 88]. Here, the upper limit range of L coming from 32Ar-32Si is 96.65 MeV. This result is overestimated in comparison with that in Ref. [28], where the upper range is confined to $ L\leq60 $ MeV. This may be due to the fact that shape deformation is ignored in our study. Besides, the pairing correlations and the limited Skyrme parameter sets used here may also have an influence on evaluating the interval range of slope parameter L.

      Figure 3.  (color online) Plots of the values of slope parameter L obtained in this work (54Ni-54Fe, 38Ca-38Ar, 36Ca-36S, and 32Ar-32Si) and in previous studies. The partially enumerated values of L from Wang $ et\; al $. [24], Brown $ et\; al $. [26], Pineda $ et\; al $. [28], An $ et\; al $. [88], Bano $ et\; al $. [32], König $ et\; al $. [28], pygmy dipole resonance (PDR) [89], electric dipole polarizability $ \alpha_{\mathrm{D}} $ [15], GW170817 [90], pion production ($ \pi^{-}/\pi^{+} $) [91], and charge exchange spin-dipole (SD) excitations [17], are also presented. The color band represents the slope parameter $ L=42.57 $$ 50.64 $ MeV obtained in this study.

      The value of the symmetry energy slope L extracted from the pygmy dipole resonance (PDR) locates at the interval range of $L=64.8\pm15.7\;\;{\rm MeV}$ [89]. The covered range is very narrow compared with that of the other cases depicted in this figure. Furthermore, the L obtained in this work can cover the uncertainty ranges derived from the nuclear electric dipole polarizability $ \alpha_{\mathrm{D}} $ [17], dense astrophysical event GW170817 [90], pion production ratios in heavy-ion collisions [91], and charge exchange spin-dipole (SD) excitations [17]. Besides, this deduced result covers most of the range of L presented in Ref. [92], in which the slope parameter captures the range of $L=47.3\pm7.8\;\;{\rm MeV}$. In summary, the extracted range of slope parameter $ L=42.57 $$ 50.64 $ MeV means that the rather soft EoS of the isospin asymmetric nuclear matter can be obtained from the implication of charge radii difference in mirror-pair nuclei.

      As shown in Fig. 2, the Pearson coefficients become lower for the mirror-pair nuclei 36Ca-36S and 38Ca-38Ar. This suggests that the profoundly theoretical uncertainties are encountered in the calibrating protocol. To clarify this phenomenon, the charge radii differences $ \Delta{R_{\mathrm{ch}}} $ of the mirror partner nuclei 36Ca-36S and 38Ca-38Ar against the slope parameter L classified by various isoscalar incompressibility coefficients K are shown in Fig. 4. From both panels of Fig. 4, one can find that the isoscalar incompressibility has an influence on the charge radii difference of mirror partner nuclei. Even from $ K\approx230 $ to $ K\approx250 $ MeV, the deviated range for L is approximately $ 38.0 $ MeV. Here, as shown in Fig. 3, we mention that the influence of the isoscalar incompressibility coefficient cannot be taken into account in evaluating the slope parameter L. Significantly, the isoscalar incompressibility may be indispensable in the calibrating protocol. As is well known, the accurate description of the density-dependent behavior of symmetry energy plays an important role in recognizing various physical mechanisms through nuclear models, such as the production of superheavy elements [93] and cluster radioactivity [94, 95]. In particular, a recent study suggests that the effective proton-neutron chemical potential difference of neutron-rich nuclei is strongly sensitive to the symmetry energy [96]. Therefore, more aspects should be taken into account in constraining the slope parameter L.

      Figure 4.  (color online) Charge radii differences $ \Delta{R_{\mathrm{ch}}} $ of mirror partner nuclei 36Ca-36S and 38Ca-38Ar plotted as a function of slope parameter L classified by various isoscalar incompressibility coefficients K. The corresponding experimental data are taken from Ref. [26].

    IV.   SUMMARY AND OUTLOOK
    • In this work, the correlation between the charge radii differences ($ \Delta{R_{\mathrm{ch}}} $) in mirror-pair nuclei 32Ar-32Si, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe and slope parameter L at saturation density is investigated based on the spherical Skyrme EDFs. The calculated results suggest that the surface pairing force can further decrease the linear correlation between $ \Delta{R_{\mathrm{ch}}} $ and L. The deduced slope parameter locates at the interval range of $ L=42.57 $$ 50.64 $ MeV under the mixed pairing interaction, that is, the rather soft EoS. This is in agreement with most of the previous studies. Meanwhile, the influence of the isoscalar nuclear matter property on $ \Delta{R_{\mathrm{ch}}} $ inevitably constrains the nuclear matter EoS. This may suggest that most of the parametrization forces are ruled out in the calibrated procedure. As shown in Fig. 4, the correlation between the charge radii difference of mirror nuclei and the slope parameter of symmetry energy is influenced by nuclear incompressibility. This seems to suggest that the Skyrme parameter sets adjusted to different slope parameters under a specific incompressibility coefficient can result in a highly linear correlation between the charge radii difference of mirror partner nuclei and the slope parameter L. Meanwhile, this may provide an alternative approach to evaluate the possible influential quantities in determining the EoS of asymmetric nuclear matter. As shown in Ref. [51], the calculated results obtained by density energy functionals consisting of Skyrme and covariant models suggest that more mirror-pair nuclei are ruled out in constraining the slope parameter L, but the superior candidates of mirror-pair nuclei 44Cr-44Ca and 46Fe-46Ca may be used to ascertain the limit range of L. Thus, more discussion should be performed in the proceeding work.

      Recent studies have shown that the charge-changing cross section measurements can be used to determine the EoS of nuclear matter with mirror nuclei [97]. This method can also be employed to investigate the charge radii of exotic nuclei [98] and the isospin interactions [99]. This provides an accessable approach to ascertain the nuclear symmetry energy. From the aspect of the charge radii difference in mirror-pair nuclei, the underlying mechanism should be taken into account in determining the nuclear symmetry energy, such as the local relations among the adjacent neighbouring nuclei [100103] and α-cluster formation [104, 105]. Thus, more high-precision charge radii data are urgently required in experimental and theoretical studies.

Reference (105)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return