Abstract:We present a new efficient way to perform hybrid density functional theory (DFT)-based electronic structure calculations. The new method uses an interpolative separable density fitting (ISDF) procedure to construct a set of numerical auxiliary basis vectors and a compact approximation of the matrix consisting of products of occupied orbitals represented in a large basis set such as the planewave basis. Such an approximation allows us to reduce the number of Poisson solves from [Formula: see text] to [Formula: see text] when we apply the exchange operator to occupied orbitals in an iterative method for solving the Kohn-Sham equations, where Ne is the number of electrons in the system to be studied. We show that the ISDF procedure can be carried out in [Formula: see text] operations, with a much smaller preconstant compared to methods used in existing approaches. When combined with the recently developed adaptively compressed exchange (ACE) operator formalism, which reduces the number of times the exchange operator needs to be updated, the resulting ACE-ISDF method significantly reduces the computational cost associated with the exchange operator by nearly 2 orders of magnitude compared to existing approaches for a large silicon system with 1000 atoms. We demonstrate that the ACE-ISDF method can produce accurate energies and forces for insulating and metallic systems and that it is possible to obtain converged hybrid functional calculation results for a 1000-atom bulk silicon within 10 min on 2000 computational cores. We also show that ACE-ISDF can scale to 8192 computational cores for a 4096-atom bulk silicon system. We use the ACE-ISDF method to geometrically optimize a 1000-atom silicon system with a vacancy defect using the HSE06 functional and computes its electronic structure. We find that that the computed energy gap from the HSE06 functional is much closer to the experimental value compared to that produced by semilocal functionals in the DFT calculations.

Machine Learning K-Means Clustering Algorithm for Interpolative Separable Density Fitting to Accelerate Hybrid Functional Calculations with Numerical Atomic Orbitals.

Machine Learning K-Means Clustering in Interpolative Separable Density Fitting Algorithm: Advancing Accurate and Efficient Cubic-Scaling Density Functional Perturbation Theory Calculations within Plane Waves

Complex-Valued K-Means Clustering of Interpolative Separable Density Fitting Algorithm for Large-Scale Hybrid Functional Enabled Ab Initio Molecular Dynamics Simulations within Plane Waves

Complex-valued K-means clustering of interpolative separable density fitting algorithm for large-scale hybrid functional enabled \textit{ab initio} molecular dynamics simulations within plane waves

Machine Learning K-Means Clustering of Interpolative Separable Density Fitting Algorithm for Accurate and Efficient Cubic-Scaling Exact Exchange Plus Random Phase Approximation within Plane Waves

Interpolative separable density fitting decomposition for accelerating Hartree-Fock exchange calculations within numerical atomic orbitals

Use of multigrids to reduce the cost of performing interpolative separable density fitting

Interpolative Separable Density Fitting for Accelerating Two-Electron Integrals: A Theoretical Perspective.

Accelerating Parallel First-Principles Excited-State Calculation by Low-Rank Approximation with K-Means Clustering

Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon

Accelerating Time-Dependent Density Functional Theory and GW Calculations for Molecules and Nanoclusters with Symmetry Adapted Interpolative Separable Density Fitting

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.

Low-Rank Approximations Accelerated Plane-Wave Hybrid Functional Calculations with K-Point Sampling

K-means Find Density Peaks in Molecular Conformation Clustering

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

The Dynamic Parallel Distribution Algorithm for Hybrid Density-Functional Calculations in HONPAS Package

Low-rank approximations for accelerating plane-wave hybrid functional calculations in unrestricted and noncollinear spin density functional theory

k-Means NANI: An Improved Clustering Algorithm for Molecular Dynamics Simulations

Hybrid MPI and OpenMP Parallel Implementation of Large-Scale Linear-Response Time-Dependent Density Functional Theory with Plane-Wave Basis Set

Statistically Optimal K-means Clustering via Nonnegative Low-rank Semidefinite Programming