Accurate Frozen Core Approximation for All-Electron Density-Functional Theory

Abstract

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over two-fold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the periodic table, and a large-scale simulation of CsPbBr3 with 2,560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below -200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

References in corpus (13)

• Canonical sampling through velocity-rescaling
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Quantum ESPRESSO toward the exascale
• QuantumATK: An integrated platform of electronic and atomic-scale modelling tools
• All-electron magnetic response with pseudopotentials: NMR chemical shifts
• Ab initio pseudopotentials for electronic structure calculations of poly-atomic systems using density-functional theory
• NWChem: Past, Present, and Future
• SIESTA: recent developments and applications
• Questaal: a package of electronic structure methods based on the linear muffin-tin orbital technique
• 103-Compound Band Structure Benchmark of Post-SCF Spin-Orbit Coupling Treatments in Density-Functional Theory
• ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
• ELSI -- An Open Infrastructure for Electronic Structure Solvers
• GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

Cited by in corpus (2)

• Equivariant Electronic Hamiltonian Prediction with Many-Body Message Passing
• Extending Quantum Computing through Subspace, Embedding and Classical Molecular Dynamics Techniques