HFB3: an axial HFB solver with Gogny forces using a 2-center HO basis (C++/Python)
Abstract
The HFB3 program solves the axial nuclear Hartree-Fock-Bogoliubov (HFB) equations using bases formed by either one or two sets of deformed Harmonic Oscillator (HO) solutions with D1-type and D2-type Gogny effective nucleon-nucleon interactions. Using two sets of HO solutions shifted along the z-axis (2-center basis) allows to accurately describe highly elongated nuclear systems while keeping a moderate basis size, making this type of basis very convenient for the description of the nuclear fission process. For the description of odd-even and odd-odd systems, the equal-filling-approximation is used. Several observables can be calculated by the program, including the mean values of the multipole moments, nuclear radii, inertia tensors following Adiabatic Time-Dependent Hartree-Fock-Bogoliubov (ATDHFB) or Generator Coordinate Method (GCM) prescriptions, local and non-local one-body densities, local and non-local pairing densities, some fission fragment properties, etc. The program can ensure that the mean values associated with some specific operators take pre-defined values (constraints). Such constraints can be set on the usual multipole moments (for protons, neutrons or total mass). This program can be used as a monoprocess and monothreaded CLI executable, or through full-featured Python bindings (available through the Python Package Index PyPI).
Keywords
Cite
@article{arxiv.2506.10745,
title = {HFB3: an axial HFB solver with Gogny forces using a 2-center HO basis (C++/Python)},
author = {N. Dubray and J. P. Ebran and P. Carpentier and M. Frosini and A. Zdeb and N. Pillet and J. Newsome and M. Verrière and G. Accorto and D. Regnier},
journal= {arXiv preprint arXiv:2506.10745},
year = {2025}
}