Axially deformed solution of the Skyrme-Hartree–Fock–Bogoliubov equations using the transformed harmonic oscillator basis (II) hfbtho v2.00d: A new version of the program

• Published in: Computer physics communications 2013-06, Vol.184 (6), p.1592-1604
• Main Authors: Stoitsov, M.V., Schunck, N., Kortelainen, M., Michel, N., Nam, H., Olsen, E., Sarich, J., Wild, S.
• Format: Article
• Language: English
• Subjects: Axial symmetry; Finite temperature HFB; Hartree–Fock–Bogoliubov; Nuclear density functional theory; Particle number projection; Self-consistent mean-field; Skyrme functionals; Transformed harmonic oscillator
• ISSN: 0010-4655; 1879-2944
• DOI: 10.1016/j.cpc.2013.01.013

We describe the new version 2.00d of the code hfbtho that solves the nuclear Skyrme-Hartree–Fock (HF) or Skyrme-Hartree–Fock–Bogoliubov (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. Program title: HFBTHO v2.00d Catalog identifier: ADUI_v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADUI_v2_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU General Public License version 3 No. of lines in distributed program, including test data, etc.: 167228 No. of bytes in distributed program, including test data, etc.: 2672156 Distribution format: tar.gz Programming language: FORTRAN-95. Computer: Intel Pentium-III, Intel Xeon, AMD-Athlon, AMD-Opteron, Cray XT5, Cray XE6. Operating system: UNIX, LINUX, WindowsXP. RAM: 200 Mwords Word size: 8 bits Classification: 17.22. Does the new version supercede the previous version?: Yes Catalog identifier of previous version: ADUI_v1_0 Journal reference of previous version: Comput. Phys. Comm. 167 (2005) 43 Nature of problem: The solution of self-consistent mean-field equations for weakly-bound paired nuclei requires a correct description of the asymptotic properties of nuclear quasi-particle wave functions. In the present implementation, this is achieved by using the single-particle wave functions of the transformed harmonic oscillator, which allows for an accurate description of deformation effects and pairing correlations in nuclei arbitrarily close to the particle drip lines. Solution method: The program uses the axial Transformed Harmonic Oscillator (THO) single- particle basis to expand quasi-particle wave functions. It iteratively diagonalizes the Hartree–Fock–Bogoliubov Hamiltonian based on generalized Skyrme-like energy densities and zero-range pairing interactions until a self-consistent soluti