Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation
arXiv:2111.02485 · doi:10.1016/j.cpc.2022.108344
Abstract
The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. The accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.
14pages, 1 figure
References in corpus (10)
- 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
- Corrections to nuclear energies and radii in finite oscillator spaces
- Surface Symmetry Energy of Nuclear Energy Density Functionals
- Universal properties of infrared oscillator basis extrapolations
- Deformed Coordinate-Space Hartree-Fock-Bogoliubov Approach to Weakly Bound Nuclei and Large Deformations
- Solution of the Skyrme HF+BCS equation on a 3D mesh. II. A new version of the Ev8 code
- Pairing renormalization and regularization within the local density approximation
- Gamow-Hartree-Fock-Bogoliubov Method: Representation of quasiparticles with Berggren sets of wave functions
- Adaptive Multi-resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure
- A New Efficient Method of Hartree-Fock-Bogoliubov Calculation for Weakly Bound Nuclei