Solution of self-consistent equations for the N3LO nuclear energy density functional in spherical symmetry. The program hosphe (v1.02)

We present solution of self-consistent equations for the N3LO 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 summaryProgram title: HOSPHE (v1.02)Catalogue identifier: AEGK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGK_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 45 809No. of bytes in distributed program, including test data, etc.: 290 514Distribution format: tar.gzProgramming language: Fortran-90Computer: PCs and workstationsOperating system: LinuxRAM: 50 MBClassification: 17.22External 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 N3LO [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 N3LO 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 N0=50N0=50 as the maximum harmonic-oscillator shell included in the basis.References:[1]B.G. Carlsson, J. Dobaczewski, M. Kortelainen, Phys. Rev. C 78 (2008) 044326.