Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (II) HFBTHO v2.00d: a new version of the program

TL;DR

The paper introduces version 2.00d of HFBTHO, a computational code for nuclear structure calculations, with new features like finite temperature, symmetry breaking, and parallelism, enhancing its accuracy and efficiency.

Contribution

It presents significant updates to the HFBTHO code, including new numerical methods, symmetry options, and parallel computing capabilities, improving nuclear modeling tools.

Findings

Enhanced code stability and convergence with the modified Broyden method.

Ability to perform multi-constraint and finite temperature calculations.

Improved computational efficiency through shared memory parallelism.

Abstract

We describe the new version 2.00d of the code HFBTHO that solves the nuclear Skyrme Hartree-Fock (HF) or Skyrme Hartree-Fock-Bogolyubov (HFB) problem by using the cylindrical transformed deformed harmonic-oscillator basis. In the new version, we have implemented the following features: (i) the modified Broyden method for non-linear problems, (ii) optional breaking of reflection symmetry, (iii) calculation of axial multipole moments, (iv) finite temperature formalism for the HFB method, (v) linear constraint method based on the approximation of the Random Phase Approximation (RPA) matrix for multi-constraint calculations, (vi) blocking of quasi-particles in the Equal Filling Approximation (EFA), (vii) framework for generalized energy density with arbitrary density-dependences, and (viii) shared memory parallelism via OpenMP pragmas.