Coordinate-space solver for superfluid many-fermion systems with the shifted conjugate-orthogonal conjugate-gradient method

Abstract

Self-consistent approaches to superfluid many-fermion systems in 3-dimensions (and subsequent time-dependent approaches) require a large number of diagonalizations of very large dimension hermitian matrices, which results in enormous computational costs. We present an approach based on the shifted conjugate-orthogonal conjugate-gradient (COCG) method for the evaluation of the Green's function, from which we subsequently extract various densities (particle number, spin, current, kinetic energy, etc.) of a nuclear system needed in self-consistent approaches. The approach eschews the construction of the quasiparticle wavefunctions and their corresponding quasiparticle energies, which are never explicitly needed in any density functional approaches. As benchmarks we present calculations for nuclei with axial symmetry, including the ground state of spherical (magic or semi-magic) and axially deformed nuclei, the saddle-point in the $^{240}$Pu constrained fission path, and a vortex in the neutron star crust.