New quantification of symmetry energy from neutron skin thicknesses of ⁴⁸Ca and ²⁰⁸Pb

Nuclear symmetry energy (NSE), which characterizes the energy cost of converting isospin symmetric nuclear matter (SNM) into pure neutron matter (PNM), plays a vital role in determining the properties of finite nuclei and neutron stars PhysRevLett.86.5647 ; Li2019 ; Steiner:2004fi ; LATTIMER2007109 ; Xiachengjun . The density dependence of the NSE, that is, $E_{\mathrm{sym}}(\rho)$, can be expanded around the saturation density $\rho_{0}$ ($\simeq 0.16$ fm^-3). The slope parameter $L$ dominates the behavior of the equation of state (EoS) for asymmetric nuclear matter in the vicinity of $\rho_{0}$. Precise knowledge about the density dependence of the NSE is difficult to obtain owing to the uncertainties arising from the varying model-dependent slope parameter $L$. Fortunately, the characteristic behaviors of the NSE can be extracted indirectly from both extensive terrestrial nuclear experiments and observed astrophysical events Li:2008gp ; Li:2014oda ; Hu:2020ujf ; PhysRevC.90.064317 ; PhysRevC.101.034303 ; Liu_2018 ; PhysRevC.108.L021303 .

Nuclear symmetry energy has been extensively used to encode the implications of the degree of isospin-asymmetry in finite nuclei. This is especially useful in the formation of the neutron skin thickness (NST) or neutron halo structure PhysRevLett.77.3963 ; PhysRevC.82.011301 ; PhysRevC.102.044313 . The quantity of NST $\Delta{R_{\mathrm{np}}}=\sqrt{\langle{r_{n}^{2}}\rangle}-\sqrt{\langle{r_{p}^{2}}\rangle}$, which is defined as the difference between the root-mean-square (rms) radii of the neutrons and protons in a heavy nucleus and is strongly correlated to the slope parameter of the NSE, $L$ Brown20005296 ; PhysRevC.64.027302 ; PhysRevC.69.024318 ; PhysRevC.72.064309 ; PhysRevLett.106.252501 ; PhysRevC.80.024316 ; PhysRevC.84.034316 ; PhysRevLett.102.122502 ; PhysRevC.81.051303 ; PhysRevLett.109.262501 ; PhysRevC.87.051306 ; PhysRevC.87.034327 ; PhysRevC.93.051303 ; ZHANG2013234 ; PhysRevC.102.044316 ; particles6010003 ; PhysRevC.91.034315 ; PhysRevC.90.064310 ; chen2006 ; PhysRevC.97.034318 . Therefore, the NST of a heavy nucleus was undertaken to provide a constraint on the EoS of neutron-rich matter around $\rho_{0}$.

The neutron radius of ²⁰⁸Pb has been determined in a laboratory by measuring the parity violating asymmetry $A_{\mathrm{PV}}$ in polarized elastic electron scattering experiments such as PREX2 PhysRevLett.126.172502 . These efforts provided the latest value of NST with significantly improved precision: $\Delta{R}_{\mathrm{np}}^{208}=0.212\sim 0.354$ fm. Moreover, a precise measurement of the NST for ⁴⁸Ca was updated by the CREX group: $\Delta{R}_{\mathrm{np}}^{48}=0.071\sim 0.171$ fm PhysRevLett.129.042501 . The reported NST of ⁴⁸Ca is relatively thin compared to the measurement obtained by the high-resolution electric polarizability experiment ($\alpha_{D}$) in the RCNP collaboration ($\Delta{R}_{\mathrm{np}}^{48}=0.14\sim 0.20$ fm) PhysRevLett.118.252501 . In contrast, the NST of $\Delta{R}_{\mathrm{np}}^{208}$ obtained by the PREX2 Collaboration is larger than that measured by RCNP ($\Delta{R}_{\mathrm{np}}^{208}=0.135\sim 0.181$ fm) PhysRevLett.107.062502 . In Ref. PhysRevLett.127.192701 , the neutron skin thickness of ²⁰⁸Pb obtained by constraining the astrophysical observables favors smaller value; for example, $\Delta{R}_{\mathrm{np}}^{208}=0.17\pm 0.04$ fm. Likewise, the optimized new functionals obtained by calibrating the $A_{\mathrm{PV}}$ and $\alpha_{D}$ values of ²⁰⁸Pb predict an NST of $\Delta{R}_{\mathrm{np}}^{208}=0.19\pm 0.02$ fm and the symmetry-energy slope $L=54\pm 8$ MeV PhysRevLett.127.232501 . Recently theoretical studies have suggested that neutron star masses and radii are more sensitive to the NST of ²⁰⁸Pb than its dipole polarizability $\alpha_{D}$ PhysRevC.107.035802 . These results challenge our understanding of nuclear force and energy density functionals (EDFs).

In Ref. TAGAMI2022106037 , 207 EoSs were employed to explore the systematic correlations between $\Delta{R}_{\mathrm{np}}^{48}$ and $L$(CREX) and between $\Delta{R}_{\mathrm{np}}^{208}$ and $L$(PREX2). The slope parameter of the NSE obtained by fitting $\Delta{R}_{\mathrm{np}}^{48}$ covers the interval range $L(\mathrm{CREX})=0\sim 50$ MeV; however, the calibrated correlation between the slope parameter $L$ and $\Delta{R}_{\mathrm{np}}^{208}$ yields $L$(PREX2)$=76\sim 165$ MeV. As mentioned in the literature, there is no overlap between $L$(CREX) and $L$(PREX2) at the one-$\sigma$ level. A combined analysis was also performed using a recent experimental determination of the parity violating asymmetry in ⁴⁸Ca and ²⁰⁸Pb PhysRevLett.129.232501 . The study demonstrated that the existing nuclear EDFs cannot simultaneously offer an accurate description of the skins of ⁴⁸Ca and ²⁰⁸Pb. The same scenario can also be encountered in Bayesian analysis, where the predicted $\Delta{R}_{\mathrm{np}}^{48}$ is closed to the CREX result, but considerably underestimates the result of $\Delta{R}_{\mathrm{np}}^{208}$ with respect to the PREX2 measurement Zhang:2022bni . Considering the isoscalar-isovector couplings in relativistic EDFs, the constraints from various high-density data cannot reconcile the recent results from PREX2 and CREX Collaboration measurements Miyatsu:2023lki . These investigations indicate that it is difficult to provide consistent constraints for the isovector components of the EoSs using existing nuclear EDFs, and further theoretical and experimental studies are urgently required YUKSEL2023137622 .

To reduce the discrepancies between the different measurements and observations, an extra term controlling the dominant gradient correction to the local functional in the isoscalar sector has been used to weaken the correlations between the properties of finite nuclei and the nuclear EoS PhysRevC.107.015801 . As demonstrated in Ref. PhysRevC.96.065805 , the influence of the isoscalar sector is nonnegligible in the analysis. Nuclear matter properties expressed in terms of their isoscalar and isovector counterparts are correlated PhysRevC.73.044320 . As noticed above, existing discussions focus on the isovector components in the EDFs model. Characteristic isoscalar quantities, such as the incompressibility of symmetric nuclear matter, are less considered when determining the slope parameter $L$ TAGAMI2022106037 . The nuclear incompressibility can be deduced from measurements of the isoscalar giant monopole resonance (ISGMR) in medium-heavy nuclei PhysRevC.70.024307 ; PhysRevC.86.054313 and multi-fragmentations of heavy ion collisions PhysRevC.70.041604 . The NSE obtained through the effective Skyrme-EDF is related to the isoscalar and isovector effective masses, which are also indirectly related to the incompressibility of symmetric nuclear matter PhysRevC.73.014313 . Although correlations between the incompressibility coefficients and isovector parameters are generally weaker than correlations between the slope parameter $L$ and NSE PhysRevC.104.054324 , quantification uncertainty due to nuclear matter incompressibility is inevitable in the evaluation. Therefore, the influence of the isoscalar nuclear matter properties is essential for evaluating slope parameter $L$.

The remainder of this paper is organized as follows. In Sec. II, we briefly describe our theoretical model. In Sec. III, we present the results and discussion. A short summary and outlook are provided in Sec. IV.

The sophisticated Skyrme-EDF, expressed as an effective zero-range force between nucleons with density- and momentum-dependent terms, has been succeeded in describing various physical phenomena RevModPhys.75.121 ; PhysRevC.95.014316 ; PhysRevC.102.054312 ; PhysRevC.90.024317 ; PhysRevC.102.014312 ; WU2022136886 ; Caotens11 ; Wenpw14 ; PhysRevC.87.064311 ; Caoquench . In this study, Skyrme-like effective interactions were calculated as follows CHABANAT1997710 ; CHABANAT1998231 :

	$\displaystyle V(\mathbf{r}_{1},\mathbf{r}_{2})$	$\displaystyle=$	$\displaystyle t_{0}(1+x_{0}\mathbf{P}_{\sigma})\delta(\mathbf{r})$		(1)
			$\displaystyle+\frac{1}{2}t_{1}(1+x_{1}\mathbf{P}_{\sigma})\left[\mathbf{P}^{\prime 2}\delta(\mathbf{r})+\delta(\mathbf{r})\mathbf{P}^{2}\right]$
			$\displaystyle+t_{2}(1+x_{2}\mathbf{P}_{\sigma})\mathbf{P}^{\prime}\cdot\delta(\mathbf{r})\mathbf{P}$
			$\displaystyle+\frac{1}{6}t_{3}(1+x_{3}\mathbf{P}_{\sigma})[\rho(\mathbf{R})]^{\alpha}\delta(\mathbf{r})$
			$\displaystyle+\mathrm{i}W_{0}\mathbf{\sigma}\cdot\left[\mathbf{P}^{\prime}\times\delta(\mathbf{r})\mathbf{P}\right],$

where $\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^{\prime}}$ 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, where $\vec{\sigma}_{1(2)}$ are the Pauli matrices. The final term denotes the spin-orbit force, where $\sigma=\vec{\sigma}_{1}+\vec{\sigma}_{2}$. Quantities $\alpha$, $t_{i}$ and $x_{i}$ ($i=0$-3) represent the effective interaction parameters of the Skyrme forces.

Generally, effective interaction parameter sets are calibrated by matching the properties of finite nuclei and nuclear matter at the saturation density. Notably,the Skyrme-EDF can provide an analytical expression of all variables characterizing infinite nuclear matter (see CHABANAT1997710 ; CHABANAT1998231 ; PhysRevC.85.035201 ; PhysRevC.82.024321 for details). The neutron skin of a heavy nucleus is regarded as the feasible indicator for probing the isovector interactions in the EoS of asymmetric nuclear matter. Thus, the neutron and proton density distributions can be self-consistently calculated using Skyrme DEFs with various parameter sets. To clarify this, we further inspected the correlations between the slope parameter $L$ and the NSTs of ⁴⁸Ca and ²⁰⁸Pb. The bulk properties were calculated using the standard Skyrme-type EDFs PhysRevC.70.024307 . The corresponding effective interactions were in accord with the calculated nuclear matter properties, such as binding energy per nucleon $E=\mathcal{E}/{\rho}$, symmetry energy $E_{\mathrm{sym}}(\rho)=\frac{1}{8}\partial^{2}(\mathcal{E}/\rho)/\partial\rho^{2}|_{\rho=\rho_{0}}$, slope parameter $L=3\rho_{0}\partial{E_{\mathrm{sym}}(\rho)}/{{\partial\rho}}|_{\rho=\rho_{0}}$, and the incompressibility coefficient $K=9\rho_{0}^{2}\partial^{2}(\mathcal{E}/\rho)/{\partial\rho^{2}}|_{\rho=\rho_{0}}$. The value of the isoscalar incompressibility $K$ from experimental data on giant monopole resonances covers a range of $230\pm 10$ MeV PhysRevLett.82.691 ; PhysRevC.69.051301 . In addition, the incompressibility of symmetric nuclear matter deduced from $\alpha$-decay properties is $K=241.28$ MeV PhysRevC.74.034302 .

The nuclear breathing model exhibits a moderate correlation with the slope of the NSE and a strong dependence on the isoscalar incompressibility coefficient $K$ of the symmetric nuclear matter Chen_2012 . The incompressibility of nuclear matter helps us understand the properties of neutron stars PhysRevC.94.052801 ; PhysRevC.104.055804 . Thus, it is essential to inspect the influence of isoscalar components on the slope parameter of the symmetry energy. To facilitate a quantitative discussion, a series of effective interaction sets classified by various nuclear incompressibility coefficients ($K=220$ MeV, 230 MeV, and 240 MeV)were employed, as shown in Table 1. Generally, analytical expressions at the saturation density $\rho_{0}$ have specific forms PhysRevC.85.035201 . Using these expressions, the density dependence of the symmetry energy can be expanded as a function of neutron excess. Under the corresponding $K$, the slope parameter $L$ and symmetry energy $E_{\mathrm{sym}}$ at the saturation density $\rho_{0}$ also cover a large range.

In Fig. 1, the NSTs of ⁴⁸Ca and ²⁰⁸Pb are determined under various effective interactions. The chosen parameter sets are classified by different incompressibility coefficients of symmetric nuclear matter, for example, $K=220$ MeV, $230$ MeV and $240$ MeV. The experimental constraint on the NST is indicated by a colored shadow. With increasing slope parameter $L$, the NST is increased monotonically, and strong linear correlations between $L$ and the NST of ⁴⁸Ca and ²⁰⁸Pb are observed. As shown in Fig. 1 (a), the linear correlations are almost similar, and the gradients for these three lines are in the ranges of $0.0008\sim 0.0009$.

Figure 1 (b) shows the related linear correlations between $\Delta{R_{\mathrm{np}}}$(²⁰⁸Pb) and the slope parameters $L$ for various nuclear matter incompressibility coefficients. However, with increasing incompressibility coefficient, the slopes of the fitted lines gradually decrease or a large deviation emerges at a high $L$. The nuclear matter EoS is conventionally defined as the binding energy per nucleon and can be expressed as Taylor-series expansion in terms of the isospin asymmetry. As suggested in Refs. PhysRevC.69.041301 ; Chen_2012 , the compression modulus of symmetric nuclear matter is sensitive to the density dependence of the NSE. With increasing neutron star mass, the correlation between $K$ and its slope $L$ increases PhysRevC.104.055804 . From this figure, we can see that the isoscalar quantity of the incompressibility coefficient has a significant influence on the determination of the slope parameter $L$ in ²⁰⁸Pb. However, for ⁴⁸Ca this influence can be ignored.

Herein, we assume that the value of $L$ is positive. Linear functions were fitted to the data classified by various nuclear matter incompressibility coefficients using the least-squares method. For $K=220$ MeV, we obtained the $L-\Delta{R_{\mathrm{np}}^{48}}$ relation as

$\displaystyle\Delta{R_{\mathrm{np}}^{48}}=0.0009L+0.1155>0.1155~{}~{}~{}~{}\mathrm{fm}.$

(2)

For $L-\Delta{R}_{\mathrm{np}}^{208}$, the linear function$K=220$ MeV is expressed as

$\displaystyle\Delta{R_{\mathrm{np}}^{208}}=0.0019L+0.0914~{}~{}~{}~{}\mathrm{fm},$

(3)

where a high correlation coefficient is located at $R=0.99$.

As suggested in Ref. TAGAMI2022106037 , the slope parameter $L$ ($0\sim 50$ MeV) deduced from $\Delta{R}_{\mathrm{np}}^{48}$ cannot overlap the interval range of the slope parameter $L$ ($76\sim 165$ MeV) deduced from $\Delta{R}_{\mathrm{np}}^{208}$. To facilitate a quantitative comparison of the experiments with these theoretical calculations, the slope parameter $L$ derived from the constraints of the NSTs of ⁴⁸Ca and ²⁰⁸Pb are presented for various nuclear matter incompressibility coefficients in Table 2. Remarkably, the gaps between $L-\Delta{R}_{\mathrm{np}}^{48}$ and $L-\Delta{R}_{\mathrm{np}}^{208}$ increase with increasing incompressibility coefficients from $K=220$ MeV to 240 MeV.

Nuclear matter properties consisting of isovector and isoscalar components are correlated with each other. Ref. PhysRevC.73.044320 suggests that there is no clear correlation between the incompressibility $K$ and NSE, and between the slope of the NSE and incompressibility $K$. The correlations between $K$ and the isovector parameters are generally weaker than those between the NST and NSE coefficients PhysRevC.69.024318 ; PhysRevC.104.054324 . As seen in Fig. 1 (b), the increasing incompressibility coefficient $K$ influences the determination of the covered range of the slope parameter $L$. Table 2 shows that the gap between $L-\Delta{R}_{\mathrm{np}}^{48}$ and $L-\Delta{R}_{\mathrm{np}}^{208}$ is smaller than the theoretical uncertainty when the nuclear incompressibility is $K=220$ MeV. This is instructive for calibrating new sets of Skyrme parameters for reproducing various nuclear matter properties as auxiliary conditions.

To facilitate the influence of incompressibility coefficient on determining the slope parameter $L$, the “data-to-data” relations between the NST of ²⁰⁸Pb and the incompressibility coefficients $K$ are presented in Fig. 2. Here, the slope parameters of the NSE were chosen to be approximately $L=34$ MeV and $L=73$ MeV. From this figure, in can be seen that the NST of ²⁰⁸Pb decreases with increasing incompressibility coefficient. This further demonstrates that the isoscalar compression modulus should be appropriately considered in the calibration protocol.

In our calculations, the upper limits of $L$ are gradually overestimated as the increasing incompressibility coefficients $K$ increased. Combined with the latest PREX2 experiment, the result extracted from the relativistic EDFs leads to a covered range of $L=106\pm 37$ MeV PhysRevLett.126.172503 . The induced slope parameter $L$ is more consistent with that obtained when the incompressibility coefficient is $K=220$ MeV.

In Refs. PhysRevC.88.011301 ; PhysRevLett.119.122502 ; PhysRevResearch.2.022035 , the highly linear correlation between the slope parameter $L$ and the differences of charge radii of mirror-partner nuclei $\Delta{R_{\mathrm{ch}}}$ was demonstrated. The nuclear charge radius of ⁵⁴Ni has been determined using the collinear laser spectroscopy PhysRevLett.127.182503 . By combining the charge radii of the mirror-pair nuclei ⁵⁴Fe, the deduced slope parameter covers the interval range $21\leq{L}\leq 88$ MeV. A recent study suggested that the upper or lower limits of $L$ may be constrained if precise data on the mirror charge radii of ⁴⁴Cr-⁴⁴Ca and ⁴⁶Fe-⁴⁶Ca are selected PhysRevC.107.034319 . In all of these studies, isoscalar nuclear matter properties were not considered. In fact, the value deduced from the relativistic and non-relativistic Skyrme EDFs with identical incompressibility coefficients $K=230$ MeV gives a narrow range of $22.50\leq{L}\leq 51.55$ MeV An:2023ahu . This is in agreement with the results in Ref. Konig:2023rwe where a soft EoS is obtained, for example, $L\leq 60$ MeV.

In atomic nuclei, the NST is regarded as a perfect signal for describing the isovector property, and is highly correlated with the slope parameter of the NSE. The difference in the charge radii of the mirror-pair nuclei and the slope of the NSE exhibits a highly linear relationship PhysRevC.97.014314 ; PhysRevLett.130.032501 ; PhysRevC.108.015802 . To facilitate the influence of the isoscalar properties on determining the EoS of nuclear matter, the data-to-data relations between the difference in charge radii $\Delta{R}_{\mathrm{ch}}$ of the mirror-pair nuclei ⁵⁴Ni-⁵⁴Fe and the NSTs of ⁴⁸Ca and ²⁰⁸Pb are shown in Fig. 3. Notably, highly linear correlations between $\Delta{R}_{\mathrm{ch}}$ and the NSTs of ⁴⁸Ca and ²⁰⁸Pb are observed.

In Fig. 3 (a), the linear functions fit the experimental data well across various incompressibility coefficients $K$, that is, the slope parameter can be constrained concurrently through the calculated NST of ⁴⁸Ca and the $\Delta{R}_{\mathrm{ch}}$ of mirror-pair nuclei ⁵⁴Ni-⁵⁴Fe. However, as shown in Fig. 3 (b), the fitting lines deviate from the cross-over region between $\Delta{R}_{\mathrm{ch}}$ and the NST of ²⁰⁸Pb except for $K=220$ MeV. Although the linear function captures a relatively narrow region, this further demonstrates the need to extract valid information about the nuclear EoS by considering the isoscalar components in the calibration procedure.

The Coulomb term does not contribute to infinite nuclear matter calculations, in which the NSE plays an essential role in determining the evolution of isospin asymmetry components. However, in atomic nuclei, the actual proton and neutron density distributions are mostly dominated by the degree of isospin asymmetry and Coulomb forces. It is evident that the competition between the Coulomb interaction and the NSE is related to the stability of dripline nuclei against nucleon emission RocaMaza:2018ujj ; PhysRevC.82.027301 . The NST is associated with the symmetry energy and significantly influenced by the NSE, which corresponds to the EoS of neutron-rich matter. Meanwhile, a strong high liner correlation between the slope parameter $L$ and the difference in the charge radii of mirror-pair nuclei is evident PhysRevC.88.011301 ; PhysRevLett.119.122502 ; PhysRevResearch.2.022035 ; PhysRevLett.127.182503 ; PhysRevC.107.034319 ; An:2023ahu . As shown in Fig. 3, this highly linear correlation extends to the NST and the difference in charge radii of mirror-pair nuclei, owing to the isospin-symmetry breaking PhysRevC.108.015802 .

As is well known, the Skyrme parameters can be characterized analytically by the isoscalar and isovector nuclear matter properties of the Hamiltonian density. More effective statistical methods have also been used to discuss the theoretical uncertainties PhysRevC.93.051303 ; Balliet_2021 ; PhysRevC.105.L021301 . In this study, we reviewed the influence of nuclear matter incompressibility on the determination of the slope parameter of symmetry energy $L$. The NSTs of ⁴⁸Ca and ²⁰⁸Pb were calculated using Skyrme EDFs. The slope parameter $L$ deduced from ²⁰⁸Pb is sensitive to the incompressibility coefficients, whereas that for ⁴⁸Ca is not. A continuous range of $L$ can be obtained if the nuclear matter is incompressible at $K=220$ MeV. This is in agreement with that in Ref. PhysRevC.104.054324 where the nuclear matter incompressibility covers the interval range of $K=223^{+7}_{-8}$ MeV. This implies that the isoscalar components should be considered when determining the slope parameter $L$. In addition, it is desirable to review the influence of the incompressibility coefficient $K$ on the determination of the slope parameter $L$ within the framework of relativistic EDFs.

The nuclear symmetry energy can be obtained using different methods and models Giuliani2023 ; FURNSTAHL200285 ; PhysRevC.93.064303 ; GAIDAROV2020122061 ; Cheng_2023 ; Caoantigdr ; Caoantigdr1 ; Cao08 ; Cao041 ; LiuMin2011 ; Fangdeqing ; GAUTAM2024122832 ; fang2023 ; Gao2023 ; Xu_2021 ; He2023 . The precise determination of the slope parameter $L$ is related to various quantities such as the charge changing cross section XU2022137333 ; ZHAO2023138269 , sub-barrier fusion cross-section, and astrophysical $S$-factor in asymmetric nuclei ghosh2023neutron . Generally, the proton and neutron density distributions are mutually determined by the isospin asymmetry and Coulomb force. The isospin symmetry breaking effect influences the determination of the charge density distributions DONG2019133 ; PhysRevC.105.L021304 ; PhysRevC.107.064302 ; SENG2023137654 . Thus, more accurate descriptions of NST and charge radii are required. In addition, the curvature of the symmetry energy $K_{\mathrm{sym}}$ PhysRevC.69.041301 and three-body interactions in the Skyrme forces PhysRevC.94.064326 may also influence the determination of the neutron skin.

This work is supported partly by the National Key R$\&$D Program of China under Grant No. 2023YFA1606401 and the National Natural Science Foundation of China under Grants No. 12135004, No. 11635003, No. 11961141004 and No. 12047513. L.-G. C. is grateful for the support of the National Natural Science Foundation of China under Grants No. 12275025, No. 11975096 and the Fundamental Research Funds for the Central Universities (2020NTST06)¡£