We present solution of self-consistent equations for the N 3LO nuclear energy density functional. We derive general expressions for the mean fields expressed as differential operators depending on densities and for the densities expressed in terms of derivatives of wave functions. These expressions are then specified to the case of spherical symmetry. We also present the computer program hosphe (v1.02), which solves the self-consistent equations by using the expansion of single-particle wave functions on the spherical harmonic oscillator basis. Program summary Program title: HOSPHE (v1.02) Catalogue identifier: AEGK_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEGK_v1_0.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.: 45 809 No. of bytes in distributed program, including test data, etc.: 290 514 Distribution format: tar.gz Programming language: Fortran-90 Computer: PCs and workstations Operating system: Linux RAM: 50 MB Classification: 17.22 External routines: LAPACK ( http://www.netlib.org/lapack/), BLAS ( http://www.netlib.org/blas/) Nature of problem: The nuclear mean-field methods constitute principal tools of a description of nuclear states in heavy nuclei. Within the Local Density Approximation with gradient corrections up to N 3LO [1], the nuclear mean-field is local and contains derivative operators up to sixth order. The locality allows for an effective and fast solution of the self-consistent equations. Solution method: The program uses the spherical harmonic oscillator basis to expand single-particle wave functions of neutrons and protons for the nuclear state being described by the N 3LO nuclear energy density functional [1]. The expansion coefficients are determined by the iterative diagonalization of the mean-field Hamiltonian, which depends non-linearly on the local neutron and proton densities. Restrictions: Solutions are limited to spherical symmetry. The expansion on the harmonic-oscillator basis does not allow for a precise description of asymptotic properties of wave functions. Running time: 50 seconds of CPU time for the ground-state of 208Pb described by using N 0 = 50 as the maximum harmonic-oscillator shell included in the basis.
Read full abstract