Exact nuclear pairing solution for large-scale configurations: I. The EP (v1.0) program at zero temperature

In this work, we present the “EP code” (version 1.0), a user-friendly and robust computational tool. It computes the exact pairing eigenvalues and eigenvectors directly from the general nuclear pairing Hamiltonian, represented using SU(2) quasi-spin algebra with basis vectors in binary representation, at zero temperature for both odd and even deformed nucleon systems. In this initial release, the sparsity and symmetry of the pairing matrix are exploited for the first time to quickly construct the pairing matrix. The ARPACK and LAPACK packages are employed for the diagonalization of large- and small-scale sparse matrices, respectively. In addition, the calculation speed for odd-nucleon systems is significantly improved by employing a novel technique to accurately identify the block containing the ground state in such odd-nucleon configurations. To ensure the high numerical stability, the Kahan compensation algorithm is employed. The current version of the EP code can efficiently expand the computational space to handle up to 26 doubly folded (deformed) single-particle levels and 26 nucleons on a standard desktop computer in approximately 10^2 s with double precision. With sufficient computational resources, the code can process up to 63 deformed single-particle levels, which can accommodate from 1 to 63 nucleon pairs. The EP v1.0 code is also designed for future extensions, including the finite-temperature and parallel computations.