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] aesn_v1_0.tar.gz(309 Kbytes)
Manuscript Title: DIRHB - a relativistic self-consistent mean-field framework for atomic nuclei.
Authors: T. Nikšić, N. Paar, D. Vretenar, P. Ring
Program title: DIRHB package (codes DIRHBS, DIRHBZ and DIRHBT)
Catalogue identifier: AESN_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 185(2014)1808
Programming language: Fortran 77.
Computer: All computers.
Operating system: All operating systems. The Makefiles are specific for a Unix OS and have to be modified for Windows.
RAM: Depends on the imposed symmetry and number of oscillator shells. For the triaxial test case it takes 300 Mb.
Keywords: Dirac-Hartree-Bogoliubov, Nuclear energy density functional, Relativistic self-consistent mean-field, Quadrupole deformation, Constrained calculation, Harmonic oscillator.
Classification: 17.22.

Nature of problem:
Ground-state properties of even-even open-shell nuclei can be calculated using the framework of self-consistent mean-field models based on relativistic energy density functionals. The structure of arbitrary heavy nuclei with spherical symmetry, axially symmetric quadrupole deformation, and triaxial quadrupole shapes, is modeled using the latest zero- and finite-range relativistic effective interactions. The particle-particle channel of the effective inter-nucleon interaction is described by a separable finite-range pairing force.

Solution method:
The current implementation of the model computes the mean- field solution of the nuclear many-body problem for even-even open-shell spherical and quadrupole deformed nuclei. The codes are used to solve the stationary relativistic Hartree-Bogoliubov equations in a self-consistent iteration scheme. At each iteration the matrix elements of the equations are updated using the modified Broyden method or the linear mixing method. The single-nucleon wave functions are expanded in a basis of spherical, axially symmetric or triaxial harmonic oscillator, depending on the assumed symmetry of the nuclear shape. For calculations that constrain the shape to specific values of the deformation parameters, the augmented Lagrangian method is used.

Time-reversal and reflection symmetries are assumed. Open-shell even-even spherical and quadrupole deformed nuclei are considered.

Running time:
Depends on the imposed symmetry and number of oscillator shells. For the test cases it runs from few seconds (spherical) up to few hours (triaxial).