High-order finite element method for atomic structure calculations

DOI10.1016/J.CPC.2023.109051arXiv2307.05856OpenAlexW4389778506MaRDI QIDQ6124635FDOQ6124635

Authors: Ondřej Čertík, John E. Pask, Isuru Fernando, Rohit Goswami, N. Sukumar, Lee A. Collins, Gianmarco Manzini, Jiří Vackář 

Publication date: 28 March 2024

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: We introduce featom, an open source code that implements a high-order finite element solver for the radial Schr"odinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number of, or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. To address the slow convergence of the $kappa=pm1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy ($10^{-8}$ Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schr"odinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers ($Z$) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples, tests, and wrappers to interface to other languages; wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

Full work available at URL: https://arxiv.org/abs/2307.05856

Recommendations

• \texttt{Dftatom}: a robust and general Schrödinger and Dirac solver for atomic structure calculations
• Accurate finite element method for atomic calculations based on density functional theory and Hartree-Fock method
• Finite-element solution of the coupled-channel Schrödinger equation using high-order accuracy approximations
• On the adaptive finite element analysis of the Kohn-Sham equations: methods, algorithms, and implementation
• Finite-element methods in electronic-structure theory

zbMATH Keywords

finite element methodDirac equationdensity functional theoryelectronic structureatomic structureKohn-Sham equationsSchrödinger equationFortran 2008

Mathematics Subject Classification ID

Statistical mechanics, structure of matter (82-XX) Quantum theory (81-XX)

Cites Work

• \texttt{Dftatom}: a robust and general Schrödinger and Dirac solver for atomic structure calculations
• DBSR\_HF: a B-spline Dirac-Hartree-Fock program
• A spectral element method for fluid dynamics: Laminar flow in a channel expansion
• Electronic Structure
• A B-spline Galerkin method for the Dirac equation
• Stabilized finite element method for the radial Dirac equation
• Variation and series approach to the Thomas-Fermi equation
• Application of the dual-kinetic-balance sets in the relativistic many-body problem of atomic structure
• Finite element method for solving the Dirac eigenvalue problem with linear basis functions

Cited In (1)

• Spectral scheme for atomic structure calculations in density functional theory

This page was built for publication: High-order finite element method for atomic structure calculations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6124635)