Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation

Abstract

Kohn-Sham density functional theory is one of the most widely used electronic structure theories. In the pseudopotential framework, uniform discretization of the Kohn-Sham Hamiltonian generally results in a large number of basis functions per atom in order to resolve the rapid oscillations of the Kohn-Sham orbitals around the nuclei. Previous attempts to reduce the number of basis functions per atom include the usage of atomic orbitals and similar objects, but the atomic orbitals generally require fine tuning in order to reach high accuracy. We present a novel discretization scheme that adaptively and systematically builds the rapid oscillations of the Kohn-Sham orbitals around the nuclei as well as environmental effects into the basis functions. The resulting basis functions are localized in the real space, and are discontinuous in the global domain. The continuous Kohn-Sham orbitals and the electron density are evaluated from the discontinuous basis functions using the discontinuous Galerkin (DG) framework. Our method is implemented in parallel and the current implementation is able to handle systems with at least thousands of atoms. Numerical examples indicate that our method can reach very high accuracy (less than 1meV) with a very small number ($4\sim 40$) of basis functions per atom.

References in corpus (3)

• Daubechies wavelets as a basis set for density functional pseudopotential calculations
• Pole-based approximation of Fermi-Dirac function
• Optimized local basis set for Kohn-Sham density functional theory

Cited by in corpus (51)

• Quantum algorithms for quantum chemistry and quantum materials science
• Adaptively Compressed Exchange Operator
• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
• The Elephant in the Room of Density Functional Theory Calculations
• DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
• Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• Model reduction methods for nuclear emulators
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Symmetry-adapted real-space density functional theory for cylindrical geometries: application to large X (X=C, Si, Ge, Sn) nanotubes
• Training and Projecting: A Reduced Basis Method Emulator for Many-Body Physics
• Chebyshev polynomial filtered subspace iteration in the Discontinuous Galerkin method for large-scale electronic structure calculations
• Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory
• Higher-order adaptive finite-element methods for orbital-free density functional theory
• Multi-sliced Gausslet Basis Sets for Electronic Structure
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• Efficient Calculations with Multisite Local Orbitals in the large-scale DFT Code CONQUEST
• SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization
• Fast real-time time-dependent hybrid functional calculations with the parallel transport gauge and the adaptively compressed exchange formulation
• Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations
• Optimized local basis set for Kohn-Sham density functional theory
• GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code
• Reproducibility of Hybrid Density Functional Calculations for Equation-of-State Properties and Band Gaps
• Orbital-enriched Flat-top Partition of Unity Method for the Schrödinger Eigenproblem
• PSelInv -- A Distributed Memory Parallel Algorithm for Selected Inversion : the Symmetric Case
• Stochastic Density Functional Theory: Real- and Energy-Space Fragmentation for Noise Reduction
• Low-cost orbital-based linear-scaling \emph{ab initio} molecular dynamics for weakly-interacting systems
• Element orbitals for Kohn-Sham density functional theory
• Localized density matrix minimization and linear scaling algorithms
• Enhancing the scalability and load balancing of the parallel selected inversion algorithm via tree-based asynchronous communication
• Robust Determination of the Chemical Potential in the Pole Expansion and Selected Inversion Method for Solving Kohn-Sham density functional theory
• The Wigner Function of Ground State and One-Dimensional Numerics
• SPARC: Simulation Package for Ab-initio Real-space Calculations
• A discontinuous Galerkin scheme for full-potential electronic structure calculations
• Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators
• Noise Tolerant Force Calculations in Density Functional Theory: A Surface Integral Approach for Wavelet-Based Methods
• Localized spectrum slicing
• Quantum Dynamics with the Parallel Transport Gauge
• Randomized sampling for basis functions construction in generalized finite element methods
• Discontinuous Galerkin discretization for quantum simulation of chemistry
• A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations
• Partition of unity finite element method for quantum mechanical materials calculations
• A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part I: Second order linear PDE
• PEXSI-$\Sigma$: A Green's function embedding method for Kohn-Sham density functional theory
• Discontinuous Galerkin method with Voronoi partitioning for Quantum Simulation of Chemistry
• Interpolative Separable Density Fitting through Centroidal Voronoi Tessellation With Applications to Hybrid Functional Electronic Structure Calculations
• A Left-Looking Selected Inversion Algorithm and Task Parallelism on Shared Memory Systems
• Parallel Transport Time-Dependent Density Functional Theory Calculations with Hybrid Functional on Summit
• A Unified Heterogeneous Implementation of Numerical Atomic Orbitals-Based Real-Time TDDFT within the ABACUS Package