Effects of neutron excess in deformed $\Lambda$ hypernuclei

For each nucleus, deformation-constrained calculations were performed over the range $\beta_2 = -0.8$ to $+1.0$ with an interval of $\delta\beta_2 = 0.02$. In terms of spatial resolution, the calculations were carried out in cylindrical coordinates with a mesh spacing of $\Delta r = 0.1,\mathrm{fm}$, and the box radius was set to $15$–$20,\mathrm{fm}$ depending on the size of the nucleus. The main data file, HFe59.dat, contains the output results for all computed points.

Each row in the table corresponds to a specific combination of a nucleus and a deformation parameter $\beta_2$. The column labels include:inbeta2: Input quadrupole deformation $\beta_2$outbeta2: Converged quadrupole deformation $\beta_2$ after iterationq2_sum: Total quadrupole moment of the $\Lambda$ hypernucleus (or normal nucleus)Hf-e: Binding energy of the $\Lambda$ hypernucleus (or normal nucleus), in MeV<r_p>: Root-mean-square radius of the proton, in fm<r_n>: Root-mean-square radius of the neutron, in fmSome extreme deformation values, particularly near $\beta_2 \approx \pm 0.8$, may lead to non-convergent or unphysical solutions.

These points have been removed during data processing. As a result, there are no NaN entries in the final dataset. The main sources of numerical error stem from mesh discretization and the convergence threshold used in the self-consistent iteration.