Numerical accuracy of mean-field calculations in coordinate space

Abstract

Background: Mean-field methods based on an energy density functional (EDF) are powerful tools used to describe many properties of nuclei in the entirety of the nuclear chart. The accuracy required on energies for nuclear physics and astrophysics applications is of the order of 500 keV and much effort is undertaken to build EDFs that meet this requirement. Purpose: The mean-field calculations have to be accurate enough in order to preserve the accuracy of the EDF. We study this numerical accuracy in detail for a specific numerical choice of representation for the mean-field equations that can accommodate any kind of symmetry breaking. Method: The method that we use is a particular implementation of 3-dimensional mesh calculations. Its numerical accuracy is governed by three main factors: the size of the box in which the nucleus is confined, the way numerical derivatives are calculated and the distance between the points on the mesh. Results: We have examined the dependence of the results on these three factors for spherical doubly-magic nuclei, neutron-rich ${}^{34}$Ne, the fission barrier of ${}^{240}$Pu and isotopic chains around Z = 50. Conclusions: Mesh calculations offer the user extensive control over the numerical accuracy of the solution scheme. By making appropriate choices for the numerical scheme the achievable accuracy is well below the model uncertainties of mean-field methods.

arXiv comment17 pages, 18 figures

References in corpus (28)

• CODATA Recommended Values of the Fundamental Physical Constants: 2010
• Nuclear Energy Density Optimization
• Nuclear energy density optimization: Large deformations
• Error Estimates of Theoretical Models: a Guide
• 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)
• Correlations beyond the mean field in Magnesium isotopes: Angular momentum projection and configuration mixing
• Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (II) HFBTHO v2.00d: a new version of the program
• The TDHF Code Sky3D
• Multidimensionally-constrained relativistic mean field models and potential energy surfaces of actinide nuclei
• Fission barriers and asymmetric ground states in the relativistic mean field theory
• Selfconsistent calculations of fission barriers in the Fm region
• Corrections to nuclear energies and radii in finite oscillator spaces
• 3D Solution of Hartree-Fock-Bogoliubov Equations for Drip-Line Nuclei
• Convergence properties of {\itshape ab initio} calculations of light nuclei in a harmonic oscillator basis
• Description of Induced Nuclear Fission with Skyrme Energy Functionals: I. Static Potential Energy Surfaces and Fission Fragment Properties
• Beyond the relativistic mean-field approximation: configuration mixing of angular momentum projected wave functions
• Microscopic calculation of 240Pu scission with a finite-range effective force
• Deformed Coordinate-Space Hartree-Fock-Bogoliubov Approach to Weakly Bound Nuclei and Large Deformations
• Solution of the Skyrme HF+BCS equation on a 3D mesh. II. A new version of the Ev8 code
• Propagation of uncertainties in the Skyrme energy-density-functional model
• Use of the Discrete Variable Representation Basis in Nuclear Physics
• Multidimensionally-constrained relativistic mean-field study of triple-humped barriers in actinides
• Coordinate Space HFB Calculations for the Zirconium Isotope Chain up to the Two-Neutron Dripline
• Error Analysis in Nuclear Density Functional Theory
• Infrared extrapolations for atomic nuclei
• Toward global beyond-mean-field calculations of nuclear masses and low-energy spectra
• Adaptive Multi-resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure

Cited by in corpus (25)

• Microscopic Theory of Nuclear Fission: A Review
• Future of Nuclear Fission Theory
• Fission Dynamics of 240Pu from Saddle-to-Scission and Beyond
• Nuclear Fission Dynamics: Past, Present, Needs, and Future
• Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider
• Projection on particle number and angular momentum: Example of triaxial Bogoliubov quasiparticle states
• 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
• Constraining the surface properties of effective Skyrme interactions
• The impact of the surface energy coefficient on the deformation properties of atomic nuclei as predicted by Skyrme energy density functionals
• Coordinate-Space Solver for Superfluid Many-Fermion Systems with Shifted Conjugate Orthogonal Conjugate Gradient Method
• Odd-odd neutron-rich rhodium isotopes studied with the double Penning trap JYFLTRAP
• Iterative approaches to the self-consistent nuclear energy density functional problem. Heavy ball dynamics and potential preconditioning
• Zero-pairing limit of Hartree-Fock-Bogoliubov reference states
• Structure of ${}^{128,129,130}$Xe through multi-reference energy density functional calculations
• Neck Rupture and Scission Neutrons in Nuclear Fission
• Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh. IIb. Fission properties of BSkG2
• Building Surrogate Models of Nuclear Density Functional Theory with Gaussian Processesand Autoencoders
• Systematic study of infrared energy corrections in truncated oscillator spaces
• Precision of finite-difference representation in 3D coordinate-space Hartree-Fock-Bogoliubov calculations
• Tensor force effect on the neutron shell closure in super-heavy elements
• Impact of choices for center-of-mass correction energy on the surface energy of Skyrme energy density functionals
• Skyrme pseudopotentials at next-to-next-to-leading order Construction of local densities and first symmetry-breaking calculations
• Framework for Polarized Superfluid Fermion Systems
• The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions