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

TL;DR

The EP v1.0 code provides an efficient, exact solution for nuclear pairing problems at zero temperature, capable of handling large configurations with improved speed and stability, suitable for future extensions.

Contribution

This work introduces the EP v1.0 code, a novel computational tool that efficiently computes exact nuclear pairing solutions for large-scale configurations using advanced matrix techniques.

Findings

Handles up to 26 nucleons and 26 levels on a desktop in ~100 seconds

Employs sparsity and symmetry for fast matrix construction

Uses Kahan algorithm for numerical stability

Abstract

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 odd 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$ seconds with double precision. With sufficient computational resources, the code can process up to 63 deformed single-particle levels, which can accomodate from 1 to 63 nucleon pairs. The EP v1.0 code is also designed for future extensions, including the finite-temperature and parallel computations.