Portal de Periódicos da CAPES

Solution of the Skyrme–Hartree–Fock–Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis.

2009; Elsevier BV; Volume: 180; Issue: 11 Linguagem: Inglês

10.1016/j.cpc.2009.08.009

1879-2944

J. Dobaczewski, W. Satuła, B. G. Carlsson, J. Engel, P. Olbratowski, P. Powałowski, M. Sadziak, Jason Sarich, N. Schunck, A. Staszczak, M. V. Stoitsov, M. Zalewski, H. Zduńczuk,

Quantum Chromodynamics and Particle Interactions

We describe the new version (v2.40h) of the code hfodd which solves the nuclear Skyrme–Hartree–Fock or Skyrme–Hartree–Fock–Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented: (i) projection on good angular momentum (for the Hartree–Fock states), (ii) calculation of the GCM kernels, (iii) calculation of matrix elements of the Yukawa interaction, (iv) the BCS solutions for state-dependent pairing gaps, (v) the HFB solutions for broken simplex symmetry, (vi) calculation of Bohr deformation parameters, (vii) constraints on the Schiff moments and scalar multipole moments, (viii) the DT2h transformations and rotations of wave functions, (ix) quasiparticle blocking for the HFB solutions in odd and odd–odd nuclei, (x) the Broyden method to accelerate the convergence, (xi) the Lipkin–Nogami method to treat pairing correlations, (xii) the exact Coulomb exchange term, (xiii) several utility options, and we have corrected three insignificant errors. Program title: HFODD (v2.40h) Catalogue identifier: ADFL_v2_2 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADFL_v2_2.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 79 618 No. of bytes in distributed program, including test data, etc.: 372 548 Distribution format: tar.gz Programming language: FORTRAN-77 and Fortran-90 Computer: Pentium-III, AMD-Athlon, AMD-Opteron Operating system: UNIX, LINUX, Windows XP Has the code been vectorised or parallelized?: Yes, vectorised RAM: 10 Mwords Word size: The code is written in single-precision for use on a 64-bit processor. The compiler option -r8 or +autodblpad (or equivalent) has to be used to promote all real and complex single-precision floating-point items to double precision when the code is used on a 32-bit machine. Classification: 17.22 Catalogue identifier of previous version: ADFL_v2_1 Journal reference of previous version: Comput. Phys. Commun. 167 (2005) 214 External routines: Lapack (http://www.netlib.org/lapack/), Blas (http://www.netlib.org), linpack (http://www.netlib/linpack/) Does the new version supersede the previous version?: Yes Nature of problem: The nuclear mean-field and an analysis of its symmetries in realistic cases are the main ingredients of a description of nuclear states. Within the Local Density Approximation, or for a zero-range velocity-dependent Skyrme interaction, the nuclear mean-field is local and velocity dependent. The locality allows for an effective and fast solution of the self-consistent Hartree–Fock equations, even for heavy nuclei, and for various nucleonic (n-particle n-hole) configurations, deformations, excitation energies, or angular momenta. Similar Local Density Approximation in the particle–particle channel, which is equivalent to using a zero-range interaction, allows for a simple implementation of pairing effects within the Hartree–Fock–Bogolyubov method. Solution method: The program uses the Cartesian harmonic oscillator basis to expand single-particle or single-quasiparticle wave functions of neutrons and protons interacting by means of the Skyrme effective interaction and zero-range pairing interaction. The expansion coefficients are determined by the iterative diagonalization of the mean field Hamiltonians or Routhians which depend non-linearly on the local neutron and proton densities. Suitable constraints are used to obtain states corresponding to a given configuration, deformation or angular momentum. The method of solution has been presented in [1]. Summary of revisions: Projection on good angular momentum (for the Hartree–Fock states) has been implemented. Calculation of the GCM kernels has been implemented. Calculation of matrix elements of the Yukawa interaction has been implemented. The BCS solutions for state-dependent pairing gaps have been implemented. The HFB solutions for broken simplex symmetry have been implemented. Calculation of Bohr deformation parameters has been implemented. Constraints on the Schiff moments and scalar multipole moments have been implemented. The DT2h transformations and rotations of wave functions have been implemented. The quasiparticle blocking for the HFB solutions in odd and odd–odd nuclei has been implemented. The Broyden method to accelerate the convergence has been implemented. The Lipkin–Nogami method to treat pairing correlations has been implemented. The exact Coulomb exchange term has been implemented. Several utility options have been implemented. Three insignificant errors have been corrected. Restrictions: The main restriction is the CPU time required for calculations of heavy deformed nuclei and for a given precision required. Unusual features: The user must have access to an implementation of the BLAS (Basic Linear Algebra Subroutines), the NAGLIB subroutine F02AXE, or LAPACK subroutines ZHPEV, ZHPEVX, or ZHEEVR, which diagonalize complex Hermitian matrices, and the LINPACK subroutines ZGEDI and ZGECO, which invert arbitrary complex matrices and calculate determinants Running time: One Hartree–Fock iteration for the superdeformed, rotating, parity conserving state of 15266Dy86 takes about six seconds on the AMD-Athlon 1600+ processor. Starting from the Woods–Saxon wave functions, about fifty iterations are required to obtain the energy converged within the precision of about 0.1 keV. In the case where every value of the angular velocity is converged separately, the complete superdeformed band with precisely determined dynamical moments J(2) can be obtained in forty minutes of CPU time on the AMD-Athlon 1600+ processor. This time can be often reduced by a factor of three when a self-consistent solution for a given rotational frequency is used as a starting point for a neighboring rotational frequency. References: J. Dobaczewski, J. Dudek, Comput. Phys. Commun. 102 (1997) 166.