Large-scale all-electron density functional theory calculations using an enriched finite element basis

Abstract

We present a computationally efficient approach to perform large-scale all-electron density functional theory calculations by enriching the classical finite element basis with compactly supported atom-centered numerical basis functions that are constructed from the solution of the Kohn-Sham (KS) problem for single atoms. We term these numerical basis functions as enrichment functions, and the resultant basis as the enriched finite element basis. The enrichment functions are compactly supported through the use of smooth cutoff functions, which enhances the conditioning and maintains the locality of the basis. The integrals involved in the evaluation of the discrete KS Hamiltonian and overlap matrix in the enriched finite element basis are computed using an adaptive quadrature grid based on the characteristics of enrichment functions. Further, we propose an efficient scheme to invert the overlap matrix by using a block-wise matrix inversion in conjunction with special reduced-order quadrature rules to transform the discrete Kohn-Sham problem to a standard eigenvalue problem. Finally, we solve the resulting standard eigenvalue problem by using a Chebyshev polynomial based filtering technique to compute the relevant eigenspectrum. We demonstrate the accuracy, efficiency and parallel scalability of the proposed method on semiconducting and heavy-metallic systems of various sizes, with the largest system containing 8694 electrons. We obtain accuracies in the ground-state energies that are within $\sim 1$mHa with ground-state energies obtained from classical finite element as well as gaussian basis sets. We observe a $50-300$ and $\sim 8$ fold reduction in the overall computational time when compared to classical finite element and gaussian basis, respectively. We also observe good parallel scalability up to 384 processors for a benchmark system comprising of 280-atom silicon nano-cluster.

References in corpus (10)

• Numerical atomic orbitals for linear scaling
• Density functional study of Au${}_n$ (n=2-20) clusters: lowest-energy structures and electronic properties
• Systematic generation of finite-range atomic basis sets for linear-scaling calculations
• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• 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
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Higher-order adaptive finite-element methods for orbital-free density functional theory
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations
• Partition of unity finite element method for quantum mechanical materials calculations

Cited by in corpus (19)

• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• Inverse Kohn-Sham Density Functional Theory: Progress and Challenges
• A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules
• Roadmap on Electronic Structure Codes in the Exascale Era
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Frontiers, challenges, and solutions in modeling of swift heavy ion effects in materials
• 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
• Electronic structure study of screw dislocation core energetics in Aluminum and core energetics informed forces in a dislocation aggregate
• Real-time time-dependent density functional theory using higher order finite element methods
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Orbital-enriched Flat-top Partition of Unity Method for the Schrödinger Eigenproblem
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations
• Solving the electronic structure problem for over 100,000 atoms in real-space
• Treecode-accelerated Green Iteration for Kohn-Sham Density Functional Theory
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Discontinuous Galerkin discretization for quantum simulation of chemistry
• Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory
• Many Electrons and the Photon Field -- The many-body structure of nonrelativistic quantum electrodynamics