Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Abstract

We present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham DFT calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block-diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit-cells and supercells. The benchmark studies show a staggering $130\times$ speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the LAPW+lo basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of $\sim$10,000 electrons.

References in corpus (20)

• Optimized norm-conserving Vanderbilt pseudopotentials
• Stable Generalized Finite Element Method (SGFEM)
• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• The Elephant in the Room of Density Functional Theory Calculations
• Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element approach
• All-electron density functional theory and time-dependent density functional theory with high-order finite elements
• Large-scale all-electron density functional theory calculations using an enriched finite element basis
• How Large is the Elephant in the Density Functional Theory Room?
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations
• Higher-order adaptive finite-element methods for orbital-free density functional theory
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids
• Real-time time-dependent density functional theory using higher order finite element methods
• Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations
• Orbital-enriched Flat-top Partition of Unity Method for the Schrödinger Eigenproblem
• Structural (B1 to B8) Phase Transition in MnO under Pressure: Comparison of All-electron and Pseudopotential Approaches
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations

Cited by in corpus (4)

• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Simulation of 24,000 Electrons Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory