Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] adui_v1_0.tar.gz(190 Kbytes)
Manuscript Title: Axially Deformed Solution of the Skyrme-Hartree-Fock-Bogolyubov Equations using The Transformed Harmonic Oscillator Basis. The program HFBTHO
Authors: M.V. Stoitsov, J. Dobaczewski, W. Nazarewicz, P. Ring
Program title: HFBTHO (v1.66p)
Catalogue identifier: ADUI_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 167(2005)43
Programming language: FORTRAN-95.
Computer: Pentium-III, Pentium-IV, AMD-Athlon, IBM Power 3, IBM Power 4, Intel Xeon.
Operating system: LINUX, Window.
RAM: 59 MB when using Nsh = 20
Word size: 64 bits
Keywords: Hartree-Fock, Hartree-Fock-Bogolyubov, Nuclear many-body problem, Skyrme interaction, Self-consistent mean-field, Quadrupole deformation, Constrained calculations, Energy surface, Pairing, Particle number projection, Nuclear radii, Quasiparticle spectra, Harmonic oscillator, Coulomb field.
PACS: 07.05.T, 21.60.-n, 21.60.Jz.
Classification: 17.22.

Nature of problem:
Solution of self-consistent mean-field equations for weakly bound paired nuclei requires a correct description of 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.

Solution method:
The program uses axial 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 the self-consistent solution is found.

Axial-, time-reversal-, and space-inversion symmetries are assumed. Only quasiparticle vacua of even-even nuclei can be calculated.

Running time:
4 seconds per iteration on an Intel Xeon 2.8 GHz processor when using Nsh = 20