Neutron skin of $^{48}$Ca deduced from reaction and interaction cross sections

In our previous paper, we predicted $sigma_{rm R}$ for $^{40-60,62,64}$Ca+ $^{12}$C scattering at 280 MeV/u, using the Kyushu (chiral) $g$-matrix folding model with the densities calculated with D1S-GHFB with and without the AMP. Interaction cross sections $sigma_{rm I}$ are available for $^{42-51}$Ca + $^{12}$C scattering, whereas $sigma_{rm R}$ are available for p+$^{48}$Ca scattering. As for $^{48}$Ca, the high-resolution $E1$ polarizability experiment ($E1$pE) yields $r_{rm skin}^{48}(E1{rm pE}) =0.14 sim 0.20~{rm fm}$. We determine $r_{rm skin}^{48}({rm exp})$ from the data on $sigma_{rm R}$ for p+$^{48}$Ca scattering and from the data on $sigma_{rm I}$ for $^{48}$Ca+$^{12}$C scattering. We use the chiral (Kyushu) $g$-matrix folding model with the densities calculated with the Gogny-D1M Hartree-Fock-Bogoliubov with the AMP. The D1M-GHFB+AMP proton and neutron densities are scaled so as to reproduce the data under the condition that the radius $r_{rm p}$ of the scaled proton density equals the data $r_{rm p}({rm exp})$ of the electron scattering. The neutron radius $r_{rm n}$ thus obtained is an experimental value. Our results are $r_{rm skin}^{48}({rm exp})=-0.031sim 0.183$fm for p+$^{48}$Ca and $0.100 sim 0.218$fm for $^{48}$Ca + $^{12}$C scattering. Using the $r_{rm skin}^{48}$-$r_{rm skin}^{208}$ relation with a high correlation coefficient $R=0.99$, we have transformed $r_{rm skin}^{208}({rm PREXII})$ and $r_{rm skin}^{208}(E1{rm pE})$ to the corresponding values $r_{rm skin}^{48}({rm tPREXII})$ and $r_{rm skin}^{48}({rm t}E1{rm pE})$. The transformed data $r_{rm skin}^{48}({rm tPREXII})=0.190 sim 0.268$fm is consistent with $r_{rm skin}^{48}=0.102 sim 0.218$fm for $^{48}$Ca + $^{12}$C. Our final result is $r_{rm skin}^{48}=0.102 sim 0.218$fm determined from $^{48}$Ca + $^{12}$C scattering.

Reaction cross section of proton scattering consistent with PREX-II

Background: The neutron skin thickness $R_{rm skin}^{rm PV}$ of PREX-II is presented in Phys. Rev. Lett. {bf 126}, 172502 (2021). The reaction cross section $sigma_R$ is useful to determine the matter radius $R_m$ and $R_{rm skin}$. For proton scattering, the reaction cross section $sigma_R$ are available for $E_{rm in} > 400$ MeV. Method and results: We determine $R_n^{rm exp}=5.727 pm 0.071$ fm and $R_m^{rm exp}=5.617 pm 0.044$ fm from $R_p^{rm exp}$ = 5.444 fm and $R_{rm skin}^{rm PV}$. The $R_p^{rm GHFB}$ calculated with Gongny-D1S HFB (GHFB) with the angular momentum projection (AMP). agrees with $R_p^{rm exp}$. The neutron density calculated with GHFB+AMP is scaled so as to $R_n^{rm scaling}=5.727$ fm. The Love-Franey $t$-matrix model with the scaled densities reproduces the data on $sigma_R$. Aim: Our aim is to find the $sigma_R$ of proton scattering consistent with $R_{rm skin}^{rm PV}$. Conclusion: The $sigma_R$ of proton scattering consistent with $R_{rm skin}^{rm PV}$ are $sigma_R^{rm exp}$ at $E_{rm in} = 534.1, 549, 806$ MeV.

Determination of matter radius and neutron skin of $^{58}$Ni from reaction cross section of proton+$^{58}$Ni scattering based on chiral $g$-matrix model

Background: Using the chiral (Kyushu) $g$-matrix folding model with the densities calculated with Gogny-HFB (GHFB) with the angular momentum projection (AMP), we determined the central values of matter radius and neutron skin from the central values of reaction cross sections $sigma_{rm R}({rm EXP})$ of p+$^{40,48}$Ca and p+$^{208}$Pb scattering. As for p+$^{58}$Ni scattering, $sigma_{rm R}({rm EXP})$ are available as a function of incident energy $E_{rm in}$. Aim: Our aim is to determine matter radius $r_{m}$ and skin $r_{rm skin}$ for $^{58}$Ni from the $sigma_{rm R}({rm EXP})$ of p+$^{58}$Ni scattering by using the Kyushu $g$-matrix folding model with the GHFB+AMP densities. Results: For p+$^{58}$Ni scattering, the Kyushu $g$-matrix folding model with the GHFB+AMP densities reproduces $sigma_{rm R}({rm EXP})$ in $8.8 leq E_{rm in} leq 81$MeV. For $E_{rm in}=81$MeV, we define the factor $F$ as $F=sigma_{rm R}({rm EXP})/sigma_{rm R}({rm AMP})=0.9775$. The $Fsigma_{rm R}({rm AMP})$ be much the same as the center values of $sigma_{rm R}({rm EXP})$ in $8.8 leq E_{rm in} leq 81$MeV. We then determine $r_{rm m}({rm EXP})$ from the center values of $sigma_{rm R}({rm EXP})$, using $sigma_{rm R}({rm EXP})=C r_{m}^{2}({rm EXP})$ with $C=r_{m}^{2}({rm AMP})/ (Fsigma_{rm R}({rm AMP}))$. The $r_{m}({rm EXP})$ thus obtained are averaged over $E_{rm in}$. The averaged value is $r_{m}({rm EXP})=3.697$fm. Eventually, we obtain $r_{rm skin}({rm EXP})=0.023$fm from $r_{rm m}=3.697$fm and $r_p({rm EXP})=3.685$fm of electron scattering.

Chiral $g$-matrix model for reaction cross section of proton+$^{40}$Ca scattering: Matter radius $r_{rm m}$ and skin $r_{rm skin}$ for $^{40}$Ca

For $^{48}$Ca, we determined $r_{m}$fm and $r_{rm skin}$fm from the central values of $sigma_{rm R}({rm EXP})$ of p+$^{48}$Ca scattering, using the chiral (Kyushu) $g$-matrix folding model with the GHFB+AMP densities. For $^{40}$Ca, Zenihiro {it et al.} determined $r_n({rm RCNP})=3.375$~fm and $r_{rm skin}({rm RCNP})=-0.01 pm 0.023$fm from the differential cross section and the analyzing powers for p+$^{40}$Ca scattering. For $^{40}$Ca, $sigma_{rm R}({rm EXP})$ are available with high accuracy. Our aim is to determine matter radius $r_{m}^{40}$ and skin $r_{rm skin}^{40}$ from $sigma_{rm R}({rm EXP})$ by using the Kyushu $g$-matrix folding model with the GHFB+AMP densities. We first determine $r_m({rm RCNP})=3.380$fm from the central value -0.01~fm of $r_{rm skin}({rm RCNP})$ and $r_p({rm RCNP})=3.385$fm. The folding model with the GHFB+AMP densities reproduces $sigma_{rm R}({rm EXP})$ in $30 leq E_{rm in} leq 180$MeV, in 2-$sigma$ level. We scale the GHFB+AMP densities so as to $r_p({rm AMP})=r_p({rm RCNP})$ and $r_n({rm AMP})=r_n({rm RCNP})$. The $sigma_{rm R}({rm RCNP})$ thus obtained agrees with the original one $sigma_{rm R}({rm AMP})$ for each $E_{rm in}$. For $E_{rm in}=180$MeV, we define $F$ as $F=sigma_{rm R}({rm EXP})/sigma_{rm R}({rm AMP})=0.929$. The $Fsigma_{rm R}({rm AMP})$ be much the same as the center values of $sigma_{rm R}({rm EXP})$ in $30 leq E_{rm in} leq 180$MeV. We then determine $r_{rm m}^{40}({rm EXP})$ from the center values of $sigma_{rm R}({rm EXP})$, using $sigma_{rm R}({rm EXP})=C r_{m}^{2}({rm EXP})$ with $C=r_{m}^{2}({rm AMP})/(Fsigma_{rm R}({rm AMP}))$. The $r_{m}({rm EXP})$ are averaged over $E_{rm in}$. The averaged value is $r_{m}({rm EXP})=3.380$fm. Eventually, we obtain $r_{rm skin}({rm EXP})=-0.01$fm from the averaged $r_{rm m}({rm EXP})$~fm and $r_p({rm PCNP})=3.385$fm.

Neutron skin $r_{rm skin}^{48}$ determined from reaction cross section of proton+$^{48}$Ca scattering

{bf Background:} Using the chiral (Kyushu) $g$-matrix folding model with the densities calculated with GHFB+AMP, we determined $r_{rm skin}^{208}=0.25$fm from the central values of $sigma_{rm R}$ of p+$^{208}$Pb scattering in $E_{rm in}=40-81$MeV. The high-resolution $E1$ polarizability experiment ($E1$pE) yields $r_{rm skin}^{48}(E1{rm pE}) =0.14-0.20$fm. The data on $sigma_{rm R}$ are available as a function of $E_{rm in}$ for $p$+$^{48}$Ca scattering. {bf Aim:} Our aim is to determine $r_{rm skin}^{48}$ from the central values of $sigma_{rm R}$ for $p$+$^{48}$Ca scattering by using the folding model. {bf Results:} As for $^{48}$Ca, we determine $r_n(E1{rm pE})=3.56$fm from the central value 0.17fm of $r_{rm skin}^{48}(E1{rm pE})$ and $r_p({rm EXP})=3.385$fm of electron scattering, and evaluate $r_m(E1{rm pE})=3.485$fm from the $r_n(E1{rm pE})$ and the $r_p({rm EXP})$ of electron scattering. The folding model with GHFB+AMP densities reproduces $sigma_{rm R}$ in $23 leq E_{rm in} leq 25.3$ MeV in one-$sigma$ level, but slightly overestimates the central values of $sigma_{rm R}$ there. In $23 leq E_{rm in} leq 25.3$MeV, the small deviation allows us to scale the GHFB+AMP densities to the central values of $r_p({rm EXP})$ and $r_n(E1{rm pE})$. The $sigma_{rm R}(E1{rm pE})$ obtained with the scaled densities almost reproduce the central values of $sigma_{rm R}$ when $E_{rm in}=23-25.3$MeV, so that the $sigma_{rm R}({rm GHFB+AMP})$ and the $sigma_{rm R}(E1{rm pE})$ are in 1-$sigma$ of $sigma_{rm R}$ there. In $E_{rm in}=23-25.3$MeV, we determine the $r_{m}({rm EXP})$ from the central values of $sigma_{rm R}$ and take the average for the $r_{m}({rm EXP})$. The averaged value is $r_{m}({rm EXP})=3.471$fm. Eventually, we obtain $r_{rm skin}^{48}({rm EXP})=0.146$fm from $r_{m}({rm EXP})=3.471$fm and $r_p({rm EXP})=3.385$fm.

Neutron skin thickness of ${}^{208}$Pb determined from reaction cross section for proton scattering

The reaction cross section $sigma_R$ is useful to determine the neutron radius $R_n$ as well as the matter radius $R_m$. The chiral (Kyushu) $g$-matrix folding model for $^{12}$C scattering on $^{9}$Be, $^{12}$C, $^{27}$Al targets was tested in the incident energy range of $30 lsim E_{rm in} lsim 400 $ MeV, and it is found that the model reliably reproduces the $sigma_R$ in $30 lsim E_{rm in} lsim 100 $ MeV and $250 lsim E_{rm in} lsim 400$ MeV. item[Aim] We determine $R_n$ and the neutron skin thickness $R_{rm skin}$ of ${}^{208}{rm Pb}$ by using high-quality $sigma_R$ data for the $p+{}^{208}{rm Pb}$ scattering in $30 leq E_{rm in} leq 100$ MeV. The theoretical model is the Kyushu $g$-matrix folding model with the densities calculated with Gongny-D1S HFB (GHFB) with the angular momentum projection (AMP). item[Results] The Kyushu $g$-matrix folding model with the GHFB+AMP densities underestimates $sigma_{rm R}$ in $30 leq E_{rm in} leq 100$~MeV only by a factor of 0.97. Since the proton radius $R_p$ calculated with GHFB+AMP agrees with the precise experimental data of 5.444 fm, the small deviation of the theoretical result from the data on $sigma_R$ allows us to scale the GHFB+AMP neutron density so as to reproduce the $sigma_R$ data. In $E_{rm in}$ = 30--100 MeV, the experimental $sigma_R$ data can be reproduced by assuming the neutron radius of ${}^{208}{rm Pb}$ as $R_n$ = $5.722 pm 0.035$ fm. item[Conclusion] The present result $R_{rm skin}$ = $0.278 pm 0.035$ fm is in good agreement with the recent PREX-II result of $r_{rm skin}$ = $0.283pm 0.071$ fm.

Reanalyses for $^{42-51}$Ca scattering on a $^{12}$C target at $280$ MeV/nucleon based on chiral $g$ folding mode with Gogny-D1S Hartree-Fock-Bogoliubov densities

In our previous paper, we predicted $r_{rm skin}$, $r_{rm p}$, $r_{rm n}$, $r_{rm m}$ for $^{40-60,62,64}$Ca after determining the neutron dripline, using the Gogny-D1S HFB (GHFB) with and without the angular momentum projection (AMP). Using the chiral $g$-matrix folding model, we predicted $sigma_{rm R}$ for Ca scattering on a $^{12}$C target at 280 MeV/nucleon, since Tanaka {it el al.} measured interaction cross sections $sigma_{rm I} (approx sigma_{rm R})$ for $^{42-51}$Ca in RIKEN. After our prediction, they determine $r_{rm m}({rm RIKEN})$, $r_{rm skin}({rm RIKEN})$, $r_{rm n}({rm RIKEN})$. In this paper, we reanalyses the $sigma_{rm I}$, since they assumed the Wood-Saxon densities for $^{42-51}$Ca. The $sigma_{rm R}$ calculated with the folding model with GHFB and GHFB+AMP densities almost reproduce the $sigma_{rm I}$. We then scale proton and neutron densities so that $r_{rm p}$ and $r_{rm n}$ may agree with the central values of $r_{rm p}(rm exp)$ and $r_{rm n}({rm RIKEN})$, respectively. The $sigma_{rm R}$ calculated with the scaled densities do not reproduce the central values of $sigma_{rm I}$ perfectly. We then determine the $r_{rm m}$ that agree with the central values of $sigma_{rm I}$, using the chiral $g$-matrix folding model. The fitted $r_{rm m}$ do not reproduce the central values of $r_{rm m}({rm RIKEN})$ perfectly, but are in one $sigma$ level. Finally, we show the $r_{rm skin}$, $r_{rm n}$ determined from the fitted $r_{rm m}$ are close to the original ones except for $r_{rm skin}^{48}$. The fitted $r_{rm skin}^{48}$ is 0.105 fm, while the central value of $r_{rm m}^{48}({rm RIKEN})$ is 0.146 fm. When we fit $r_{rm m}$ to the upper bound of $sigma_{rm I}$, the fitted $r_{rm skin}^{48}$ is 0.164~fm and near the central vale 0.17 fm of the high-resolution $E1$ polarizability experiment.

Neutron skin of $^{48}$Ca consistent with experimental data on skins

[Background]: In our previous paper, we predicted $r_{rm skin}$, $r_{rm p}$, $r_{rm n}$, $r_{rm m}$ for $^{40-60,62,64}$Ca after determining the neutron dripline, using the Gogny-D1S HFB with and without the angular momentum projection (AMP). We found that effects of the AMP are small. Very lately, Tanaka {it et al.} measured interaction cross sections $sigma_{rm I}$ for $^{42-51}$Ca, determined $r_{rm m}$ from the $sigma_{rm I}$, and deduced skin $r_{rm skin}$ and $r_{rm n}$ from the $r_{rm m}$ and the $r_{rm p}(rm {exp})$ evaluated from the electron scattering. Comparing our results with the data, we find for $^{42-48}$Ca that GHFB and GHFB+AMP reproduce the data on $r_{rm skin}$, $r_{rm n}$, $r_{rm m}$, but not for $r_{rm p}(rm {exp})$. [Aim]: Our purpose is to determine a value of $r_{rm skin}^{48}$ by using GHFB+AMP and the constrained GHFB (cGHFB) in which the calculated value is fitted to $r_{rm p}(rm {exp})$. [Results]: For $^{42,44,46,48}$Ca, cGHFB hardly changes $r_{rm skin}$, $r_{rm m}$, $r_{rm n}$ calculated with GHFB+AMP, except for $r_{rm skin}^{48}$. For $r_{rm skin}^{48}$, the cGHFB result is $r_{rm skin}^{48}=0.190$fm, while $r_{rm skin}^{48}=0.159$fm for GHFB+AMP. We should take the upper and the lower bound of GHFB+AMP and cGHFB. The result $r_{rm skin}^{48}=0.159-0.190$fm consists with the $r_{rm skin}^{48}(sigma_{rm I})$ and the data $r_{rm skin}^{48}(rm $E1$pE)$ obtained from high-resolution $E1$ polarizability experiment ($E1$pE). Using the $r_{rm skin}^{48}$-$r_{rm skin}^{208}$ relation with strong correlation of Ref.[3], we transform the data $r_{rm skin}^{208}$ determined by PREX and $E1$pE to the corresponding values, $r_{rm skin}^{48}(rm tPREX)$ and $r_{rm skin}^{48}(rm t$E1$pE)$. Our result is consistent also for $r_{rm skin}^{48}(rm tPREX)$ and $r_{rm skin}^{48}(rm t$E1$pE)$.

Determination of symmetry energy from experimental and observational constraints; prediction on CREX

We determine the symmetry energy $S_{rm sym}(rho)$ at the saturation density $rho=rho_0$, particularly for three quantities of $J equiv S_{rm sym}(rho_0)$, the sloop $L$, the curvature $K_{rm sym}$. The values are evaluated from the observational constraint $M_{rm max} ge 2{rm M}_{rm sun}$ of neutron star and the experimental constraint on the neutron skin thickness $r_{rm skin}^{208}$. The experimental constraint is obtained by taking the weighted mean and its error for three reliable measurements of PREX, electric-dipole polarizability, pygmy strength. One result is $J=28-34$MeV, $L= 43-67$MeV, $K_{rm sym}=(-150)-(-3)$MeV based on both the observational constraint and the experimental constraint. The other is $J=27-35$MeV, $L= 43-67$MeV, $K_{rm sym}=(-270)-(29)$MeV based on the experimental constraint only. We also determine the relation between $r_{rm skin}^{48}$ and $r_{rm skin}^{208}$ with the correlation coefficient $R=0.99$. Using the $r_{rm skin}^{48}$--$r_{rm skin}^{208}$ relation, we transform the central value of PREX to the corresponding value for $r_{rm skin}^{48}$ and predict a result of CREX.

$g$-matrix folding-model approach to reaction cross sections for scattering of Ca isotopes on a C target

We first predict the ground-state properties of Ca isotopes, using the Gogny-D1S Hartree-Fock-Bogoliubov (GHFB) with and without the angular momentum projection (AMP). We find that $^{64}$Ca is an even-dripline nucleus and $^{59}$Ca is an odd-dripline nucleus, using $A$ dependence of the one-neutron separation energy $S_{1}$ and the two-neutron separation energy, $S_{2}$. As for $S_{1}$, $S_{2}$ and the binding energies $E_{rm B}$, our results agree with the experimental data in $^{40-58}$Ca. As other ground-state properties of $^{40-60,62,64}$Ca, we predict charge, proton, neutron, matter radii, neutron skin and deformation. As for charge radii, our results are consistent with the experimental data in $^{40-52}$Ca. For $^{48}$Ca, our results on proton, neutron, matter radii agree with the experimental data. Very lately, Tanaka et. al. measured interaction cross sections for $^{42-51}$Ca scattering on a $^{12}$C target at an incident energy per nucleon of $E_{rm lab}=280$MeV. Secondly, we predict reaction cross sections $sigma_{rm R}$ for $^{40-60,62,64}$Ca, using a chiral $g$-matrix double-folding model (DFM). To show the reliability of the present DFM for $sigma_{rm R}$, we apply the DFM for the data on $^{12}$C scattering on $^{9}$Be, $^{12}$C, $^{27}$Al targets in $30 < E_{rm lab} < 400 $MeV, and show that the present DFM is good in $30 < E_{rm lab} < 100 $MeV and $250 < E_{rm lab} < 400 $MeV. For $110 < E_{rm lab} < 240 $MeV, our results have small errors. To improve the present DFM for $sigma_{rm R}$, we propose two prescriptions.