The coordinate space formulation of the Hartree-Fock-Bogoliubov (HFB) method enables self-consistent treatment of mean-field and pairing in weakly bound systems whose properties are affected by the particle continuum space. Of particular interest are neutron-rich, deformed drip-line nuclei which can exhibit novel properties associated with neutron skin. To describe such systems theoretically, we developed an accurate 2D lattice Skyrme-HFB solver {hfbax} based on B-splines. Compared to previous implementations, we made a number of improvements aimed at boosting the solvers performance. These include: explicit imposition of axiality and space inversion, use of the modified Broydens method to solve self-consistent equations, and a partial parallelization of the code. {hfbax} has been benchmarked against other HFB solvers, both spherical and deformed, and the accuracy of the B-spline expansion was tested by employing the multiresolution wavelet method. Illustrative calculations are carried out for stable and weakly bound nuclei at spherical and very deformed shapes, including constrained fission pathways. In addition to providing new physics insights, {hfbax} can serve as a useful tool to assess the reliability and applicability of coordinate-space and configuration-space HFB solvers, both existing and in development.
Download