Abstract:By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semi-local density functional theory (DFT), and thereby furnish a more accurate and reliable description of the electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we have devised a linear-scaling yet formally exact approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions, MLWFs) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems at the hybrid DFT level. We provide a comprehensive description of the exx algorithm, which is currently implemented in the Quantum ESPRESSO program and employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize high-performance computing (HPC) resources. This is followed by a critical assessment of the accuracy and parallel performance of this approach when performing AIMD simulations of liquid water in the canonical ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT based AIMD simulations of condensed-phase systems containing 500-1000 atoms with a walltime cost that is comparable to semi-local DFT. In doing so, exx takes us closer to routinely performing AIMD simulations of large-scale condensed-phase systems for sufficiently long timescales at the hybrid DFT level of theory.

Generalized Linear-Scaling Localized-Density-Matrix Method

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System

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

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

Density functional calculation of many-electron systems in cartesian coordinate grid

Linear Scaling Local Correlation Approach for Solving the Coupled Cluster Equations of Large Systems.

Extending Density Functional Embedding Theory for Covalently Bonded Systems.

Linear Scaling Calculation of Maximally Localized Wannier Functions with Atomic Basis Set.

Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections

Semilocal density functional obeying a strongly tightened bound for exchange

Diagonalization-free self-consistent field approach with localized molecular orbitals

Local Self-Energies for V and Pd Emergent from a Non-Local LDA+FLEX Implementation

Global Hybrid Multiconfiguration Pair-Density Functional Theory

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

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

Orbital densities functional

Linear Scaling Time-Dependent Density-Functional Tight-Binding Method for Absorption Spectra of Large Systems

Basis-Set Limit CCSD(T) Energies for Large Molecules with Local Natural Orbitals and Reduced-Scaling Basis-Set Corrections

Density functional theory of electron transfer beyond the Born-Oppenheimer approximation: Case study of LiF

Enhancing robustness and efficiency of density matrix embedding theory via semidefinite programming and local correlation potential fitting

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