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.

The decoupled DFT-$\frac{1}{2}$ method for defect excitation energies

Numerical Fitting of Transient Decays in the High Defect Density Limit

Comparison of charged-defect finite-size supercell correction methods in a general framework

DFT-1/2 and shell DFT-1/2 methods: Electronic structure calculation for semiconductors at LDA complexity

Accurate and convergent energetics of color centers by wavefunction theory

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

The role of the exact Hartree-Fock exchange in the investigation of defects in crystalline systems.

Optical Properties of Defects in Solids via Quantum Embedding with Good Active Space Orbitals

Linear Scaling DFT for defects in metals

Defect Energy Levels in Density Functional Calculations: Alignment and Band Gap Problem

Quantum embedding methods for correlated excited states of point defects: Case studies and challenges

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Comparative study of defect transition energy calculation methods: The case of oxygen vacancy in In 2 O 3 and ZnO

Energy density in density functional theory: Application to crystalline defects and surfaces

Analysis of nearly planar defects using the Thomas--Fermi--von Weiszacker model

Fundamental Principles for Calculating Charged Defect Ionization Energies in Ultrathin Two-Dimensional Materials

Density fitting in periodic systems: application to TDHF in diamond and oxides

Fingerprinting defects in diamond: Partitioning the vibrational spectrum

Aperiodic defects in periodic solids

Dispersion with Fixed Diagonal Matrices: Exchange energy correction and an assessment of the Becke-Roussel exchange hole

A New View on Density Corrected DFT: Can One Get a Better Answer for a Good Reason?