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.

Local Hybrid Divide-and-Conquer Method for the Computation of Medium and Large Molecules.

A Localized Molecular-Orbital Assembler Approach For Hartree-Fock Calculations Of Large Molecules

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

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

Dcmb That Combines Divide-And-Conquer and Mixed-Basis Set Methods for Accurate Geometry Optimizations, Total Energies, and Vibrational Frequencies of Large Molecules

Accelerating Linear-Response Time-Dependent Hybrid Density Functional Theory with Low-Rank Decomposition Techniques in the Plane-Wave Basis

Accurate Localized Resolution Of Identity Approach For Linear-Scaling Hybrid Density Functionals And For Many-Body Perturbation Theory

Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections

GW with hybrid functionals for large molecular systems

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

Low-rank Approximations to Accelerate Hybrid Functional Enabled Real-Time Time-Dependent Density Functional Theory Within Plane Waves

Extended Energy Divide-and-Conquer Method Based on Charge Conservation.

Deterministic/Fragmented-Stochastic Exchange for Large Scale Hybrid DFT Calculations

Efficient real space formalism for hybrid density functionals

Short- and long-range corrected hybrid density functionals with the D3 dispersion corrections

Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations

A General and Transferable Local Hybrid Functional for Electronic Structure Theory and Many-Fermion Approaches

Efficient HF exchange evaluation through Fourier convolution in Cartesian grid for orbital-dependent density functionals

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

Consistent Improvement on Both Exchange and Correlation within Hybrid Multiconfiguration Pair-density Functional Theory: tB4LYP

Hybrid Coupled Cluster Methods: Combining Active Space Coupled Cluster Methods with Coupled Cluster Singles, Doubles, and Perturbative Triples