HFB3: an axial HFB solver with Gogny forces using a 2-center HO basis (C++/Python)

TL;DR

HFB3 is a computational tool that solves axial nuclear HFB equations using a 2-center harmonic oscillator basis, effectively modeling elongated nuclear systems and fission processes with flexible observables and constraints.

Contribution

The paper introduces HFB3, a novel axial HFB solver utilizing a 2-center HO basis with Gogny forces, optimized for describing elongated nuclei and fission, with Python integration.

Findings

Accurately describes highly elongated nuclear systems.

Efficiently models nuclear fission processes.

Provides comprehensive observables and constraints.

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).