Effects of neutron excess in deformed $\Lambda$ hypernuclei

This dataset is generated from deformation-constrained Skyrme–Hartree–Fock (SHF) calculations in axially symmetric coordinates, aiming to study the structural evolution of $\Lambda$ hypernuclei under various deformation conditions. The calculations were performed using a self-developed Fortran program running on a Linux system, with large-scale batch computations carried out on a high-performance computing platform. A quadrupole deformation constraint ($\beta_2$) was introduced in the calculations, and the energy surfaces and related physical quantities were obtained through point-by-point scanning. The adopted Skyrme interaction parameter set includes central and spin–orbit terms, and all results were output after convergence of both the single-particle energies and the total energy. The dataset covers a series of $\Lambda$ hypernuclei with mass numbers ranging from $A = 12$ to $A = 208$. 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.