We describe the program HFBTHO for axially deformed configurational Hartree–Fock–Bogolyubov calculations with Skyrme-forces and zero-range pairing interaction using Harmonic-Oscillator and/or ...Transformed Harmonic-Oscillator states. The particle-number symmetry is approximately restored using the Lipkin–Nogami prescription, followed by an exact particle number projection after the variation. The program can be used in a variety of applications, including systematic studies of wide ranges of nuclei, both spherical and axially deformed, extending all the way out to nucleon drip lines.
Title of the program: HFBTHO (v1.66p)
Catalogue number: ADUI
Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland
Program summary URL:
http://cpc.cs.qub.ac.uk/summaries/ADUI
Licensing provisions: none
Computers on which the program has been tested: Pentium-III, Pentium-IV, AMD-Athlon, IBM Power 3, IBM Power 4, Intel Xeon
Operating systems: LINUX, Windows
Programming language used: FORTRAN-95
Memory required to execute with typical data: 59 MB when using
N
sh
=
20
No. of bits in a word: 64
No. of processors used: 1
Has the code been vectorized?: No
No. of bytes in distributed program, including test data, etc.: 195 285
No. of lines in distributed program: 12 058
Distribution format: tar.gz
Nature of physical problem: The solution of self-consistent mean-field equations for weakly bound paired nuclei requires a correct description of the asymptotic properties of nuclear quasiparticle 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.
Method of solution: The program uses the axially Transformed Harmonic Oscillator (THO) single-particle basis to expand quasiparticle wave functions. It iteratively diagonalizes the Hartree–Fock–Bogolyubov Hamiltonian based on the Skyrme-forces and zero-range pairing interaction until a self-consistent solution is found.
Restrictions on the complexity of the problem: Axial-, time-reversal-, and space-inversion symmetries are assumed. Only quasiparticle vacua of even–even nuclei can be calculated.
Typical running time: 4 s per iteration on an Intel Xeon 2.8 GHz processor when using
N
sh
=
20
Unusual features of the program: none
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
.
The augmented Lagrangiam method (ALM), widely used in quantum chemistry constrained optimization problems, is applied in the context of the nuclear Density Functional Theory (DFT) in the ...self-consistent constrained Skyrme Hartree-Fock-Bogoliubov (CHFB) variant. The ALM allows precise calculations of multi-dimensional energy surfaces in the space of collective coordinates that are needed to,
e.g.
, determine fission pathways and saddle points; it improves the accuracy of computed derivatives with respect to collective variables that are used to determine collective inertia; and is well adapted to supercomputer applications.
Full text
Available for:
DOBA, EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, IZUM, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, SIK, UILJ, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Calculations of nuclear masses, using nuclear density functional theory, are presented for even-even nuclei spanning the nuclear chart. The resulting binding energy differences can be interpreted in ...terms of valence proton-neutron interactions. These are compared globally, regionally, and locally with empirical values. Overall, excellent agreement is obtained. Discrepancies highlight neglected degrees of freedom and can point to improved density functionals.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
The program of systematic large-scale self-consistent nuclear mass calculations that is based on the nuclear density functional theory represents a rich scientific agenda that is closely aligned with ...the main research directions in modern nuclear structure and astrophysics, especially the radioactive nuclear beam physics. The quest for the microscopic understanding of the phenomenon of nuclear binding represents, in fact, a number of fundamental and crucial questions of the quantum many-body problem, including the proper treatment of correlations and dynamics in the presence of symmetry breaking. Recent advances and open problems in the field of nuclear mass calculations are presented and discussed.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
The Hartree-Fock-Bogoliubov equations for deformed weakly bound nuclei are solved using a basis of Poschl-Teller-Ginocchio and Bessel/Coulomb wave functions, which possess proper asymptotic behavior. ...Results compare well with standard Hartree-Fock- Bogoliubov calculations using box boundary condition or transformed harmonic oscillator basis.
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 solution is found. A previous version of the program was presented in: M.V. Stoitsov, J. Dobaczewski, W. Nazarewicz, P. Ring, Comput. Phys. Commun. 167 (2005) 43–63.
Reasons for new version:
Version 2.00d of HFBTHO provides a number of new options such as the optional breaking of reflection symmetry, the calculation of axial multipole moments, the finite temperature formalism for the HFB method, optimized multi-constraint calculations, the treatment of odd–even and odd–odd nuclei in the blocking approximation, and the framework for generalized energy density with arbitrary density-dependences. It is also the first version of HFBTHO to contain threading capabilities.
Summary of revisions:1.The modified Broyden method has been implemented,2.Optional breaking of reflection symmetry has been implemented,3.The calculation of all axial multipole moments up to λ=8 has been implemented,4.The finite temperature formalism for the HFB method has been implemented,5.The linear constraint method based on the approximation of the Random Phase Approximation (RPA) matrix for multi-constraint calculations has been implemented,6.The blocking of quasi-particles in the Equal Filling Approximation (EFA) has been implemented,7.The framework for generalized energy density functionals with arbitrary density-dependence has been implemented,8.Shared memory parallelism via OpenMP pragmas has been implemented.
Restrictions:
Axial- and time-reversal symmetries are assumed.
Unusual features:
The user must have access to (i)the LAPACK subroutines DSYEVD, DSYTRF and DSYTRI, and their dependences, which compute eigenvalues and eigenfunctions of real symmetric matrices,(ii)the LAPACK subroutines DGETRI and DGETRF, which invert arbitrary real matrices, and(iii)the BLAS routines DCOPY, DSCAL, DGEMM and DGEMV for double-precision linear algebra (or provide another set of subroutines that can perform such tasks). The BLAS and LAPACK subroutines can be obtained from the Netlib Repository at the University of Tennessee, Knoxville: http://netlib2.cs.utk.edu/.
Running time:
Highly variable, as it depends on the nucleus, size of the basis, requested accuracy, requested configuration, compiler and libraries, and hardware architecture. An order of magnitude would be a few seconds for ground-state configurations in small bases Nmax≈8−12, to a few minutes in very deformed configuration of a heavy nucleus with a large basis Nmax>20.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
Within the nuclear density functional theory (DFT) we study the effect of reflection-asymmetric shapes on ground-state binding energies and binding energy differences. To this end, we developed the ...new DFT solver AXIALHFB that uses an approximate second-order gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even-even isotopes of radium and thorium.
We present an efficient method for calculating strength functions using the finite amplitude method (FAM) for deformed superfluid heavy nuclei within the framework of the nuclear density functional ...theory. We demonstrate that FAM reproduces strength functions obtained with the fully self-consistent quasi-particle random-phase approximation (QRPA) at a fraction of computational cost. As a demonstration, we compute the isoscalar and isovector monopole strength for strongly deformed configurations in ^{240}Pu by considering huge quasi-particle QRPA spaces. Our approach to FAM, based on Broyden's iterative procedure, opens the possibility for large-scale calculations of strength distributions in well-bound and weakly bound nuclei across the nuclear landscape.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM