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.

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

Towards adiabatic-connection interpolation model with broader applicability

High-Level Correlated Approach to the Jellium Surface Energy, Without Uniform-Electron-Gas Input

Assessment of interaction-strength interpolation formulas for gold and silver clusters

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

Møller-Plesset and density-fixed adiabatic connections for a model diatomic system at different correlation regimes

Energy densities in the strong-interaction limit of density functional theory

A similarity renormalization group approach to Green's function methods

Semilocal Meta-GGA Exchange-Correlation Approximation From Adiabatic Connection Formalism: Extent and Limitations

Energy Density Functionals From the Strong-Coupling Limit Applied to the Anions of the He Isoelectronic Series

Efficient anisotropic Migdal-Eliashberg calculations with the Intermediate Representation basis and Wannier interpolation

A quantitative theory and atomistic simulation study on the soft-sphere crystal-melt interfacial properties. I. Kinetic coefficients

Ab initio path integral Monte Carlo simulations of the uniform electron gas on large length scales

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

Effective interaction between guest charges immersed in 2D jellium

Interacting-bath dynamical embedding for capturing non-local electron correlation in solids

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

Optimized numerical Density Functional Theory calculation of rotationally symmetric jellium model systems

Strongly correlated metal interfaces in the Gutzwiller approximation

Multiscale simulation of rarefied gas dynamics via direct intermittent GSIS-DSMC coupling