Iterative approaches to the self-consistent nuclear energy density functional problem. Heavy ball dynamics and potential preconditioning

Abstract

Large-scale applications of energy density functional (EDF) methods depend on fast and reliable algorithms to solve the associated non-linear self-consistency problem. When dealing with large single-particle variational spaces, existing solvers can become very slow, and their performance dependent on manual fine-tuning of numerical parameters. In addition, convergence can sensitively depend on particularities of the EDF's parametrisation under consideration. Using the widely-used Skyrme EDF as an example, we investigate the impact of the parametrisation of the EDF, both in terms of the operator structures present and the size of coupling constants, on the convergence of numerical solvers. We focus on two aspects of the self-consistency cycle, which are the diagonalisation of a fixed single-particle Hamiltonian on one hand and the evolution of the mean-field densities and potentials on the other. Throughout the article we use a coordinate-space representation, for which the behaviour of algorithms can be straightforwardly analysed. We propose two algorithmic improvements that are easily implementable in existing solvers, heavy-ball dynamics and potential preconditioning. We demonstrate that these methods can be made virtually parameter-free, requiring no manual fine-tuning to achieve near-optimal performance except for isolated cases. The combination of both methods decreases substantially the CPU time required to obtain converged results. The improvements are illustrated for the MOCCa code that solves the self-consistent HFB problem in a 3d coordinate space representation for parametrisations of the standard Skyrme EDF at next-to-leading order in gradients and its extension to next-to-next-to-leading order.

arXiv comment29 pages, 18 figures

References in corpus (34)

• Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. XII: Stiffness and stability of neutron-star matter
• Nuclear Energy Density Optimization
• Skyrme-Hartree-Fock-Bogoliubov nuclear mass formulas: Crossing the 0.6 MeV threshold with microscopically deduced pairing
• Nuclear energy density optimization: Large deformations
• Pairing gaps from nuclear mean-field models
• The tensor part of the Skyrme energy density functional. I. Spherical nuclei
• Axially Deformed Solution of the Skyrme-Hartree-Fock-Bogolyubov Equations using The Transformed Harmonic Oscillator Basis. The program HFBTHO (v1.66p)
• Local Density Approximation for proton-neutron pairing correlations. I. Formalism
• Nuclear energy density optimization: Shell structure
• The TDHF Code Sky3D
• Isovector splitting of nucleon effective masses, ab-initio benchmarks and extended stability criteria for Skyrme energy functionals
• Superdeformed rotational bands in the Mercury region; A Cranked Skyrme-Hartree-Fock-Bogoliubov study
• Tensor part of the Skyrme energy density functional. II: Deformation properties of magic and semi-magic nuclei
• 3D Solution of Hartree-Fock-Bogoliubov Equations for Drip-Line Nuclei
• Solution of the Skyrme-Hartree-Fock-Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VII) HFODD (v2.49t): a new version of the program
• Local nuclear energy density functional at next-to-next-to-next-to-leading order
• Broyden's Method in Nuclear Structure Calculations
• Application of the gradient method to Hartree-Fock-Bogoliubov theory
• Effective pseudopotential for energy density functionals with higher order derivatives
• Solution of the Skyrme HF+BCS equation on a 3D mesh. II. A new version of the Ev8 code
• Solving the Schrödinger eigenvalue problem by the imaginary time propagation technique using splitting methods with complex coefficients
• Solution of the Skyrme-Hartree-Fock equations in the Cartesian deformed harmonic oscillator basis. (II) The program HFODD
• The tensor part of the Skyrme energy density functional. III. Time-odd terms at high spin
• Spurious finite-size instabilities in nuclear energy density functionals
• Linear response of homogeneous nuclear matter with energy density functionals
• Skyrme functional from a three-body pseudo-potential of second-order in gradients. Formalism for central terms
• Skyrme pseudo-potential-based EDF parametrisation for spuriousity-free MR-EDF calculations
• Adaptive Multi-resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure
• Constraining the surface properties of effective Skyrme interactions
• Numerical accuracy of mean-field calculations in coordinate space
• Solution of Hartree-Fock-Bogoliubov equations and fitting procedure using N2LO Skyrme pseudo-potential in spherical symmetry
• The impact of the surface energy coefficient on the deformation properties of atomic nuclei as predicted by Skyrme energy density functionals
• Systematic study of infrared energy corrections in truncated oscillator spaces
• Tri-axial deformation in nuclei with realistic NN interactions

Cited by in corpus (14)

• Future of Nuclear Fission Theory
• Symmetry restoration in mean-field approaches
• Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D Mesh: Effect of triaxial shape
• Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: II. Time-reversal symmetry breaking
• Variational approximations to exact solutions in shell-model valence spaces: calcium isotopes in the pf-shell
• Odd-odd neutron-rich rhodium isotopes studied with the double Penning trap JYFLTRAP
• Symmetry-projected variational calculations with the numerical suite TAURUS I. Variation after particle-number projection
• Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh. IIb. Fission properties of BSkG2
• Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation
• Finite-temperature mean-field approximations for shell model Hamiltonians: the code HF-SHELL
• Generator coordinate method for nuclear octupole excitations: status and perspectives
• Impact of choices for center-of-mass correction energy on the surface energy of Skyrme energy density functionals
• Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method
• Skyrme pseudopotentials at next-to-next-to-leading order Construction of local densities and first symmetry-breaking calculations