dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations

Abstract

A robust and general solver for the radial Schrödinger, Dirac, and Kohn--Sham equations is presented. The formulation admits general potentials and meshes: uniform, exponential, or other defined by nodal distribution and derivative functions. For a given mesh type, convergence can be controlled systematically by increasing the number of grid points. Radial integrations are carried out using a combination of asymptotic forms, Runge-Kutta, and implicit Adams methods. Eigenfunctions are determined by a combination of bisection and perturbation methods for robustness and speed. An outward Poisson integration is employed to increase accuracy in the core region, allowing absolute accuracies of $10^{-8}$ Hartree to be attained for total energies of heavy atoms such as uranium. Detailed convergence studies are presented and computational parameters are provided to achieve accuracies commonly required in practice. Comparisons to analytic and current-benchmark density-functional results for atomic number $Z$ = 1--92 are presented, verifying and providing a refinement to current benchmarks. An efficient, modular Fortran 95 implementation, \ttt{dftatom}, is provided as open source, including examples, tests, and wrappers for interface to other languages; wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

Submitted to Computer Physics Communication on August 27, 2012, revised February 1, 2013

References in corpus (2)

• A B-spline Galerkin method for the Dirac equation
• High-order finite element method for atomic structure calculations

Cited by in corpus (12)

• 103-Compound Band Structure Benchmark of Post-SCF Spin-Orbit Coupling Treatments in Density-Functional Theory
• Accurate absolute and relative core-level binding energies from $GW$
• Wide Ranging Equation of State with Tartarus: a Hybrid Green's Function/Orbital based Average Atom Code
• Relativistic correction scheme for core-level binding energies from $GW$
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Electronic friction coefficients from the atom-in-jellium model for $Z=1-92$
• High-order finite element method for atomic structure calculations
• Relativistic EELS scattering cross-sections for microanalysis based on Dirac solutions
• Average Atom Model with Siegert States
• Influence of atomic modeling on electron capture and shaking processes
• SClib, a hack for straightforward embedded C functions in Python
• Toward Modern Fortran Tooling and a Thriving Developer Community