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.

Large-scale efficient Langevin dynamics, and why it works

Self-consistent field theory based molecular dynamics with linear system-size scaling

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Semi-classical generalized Langevin equation for equilibrium and nonequilibrium molecular dynamics simulation

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

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

Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Large Scale and Quantum Accurate Molecular Dynamics Simulation:Liquid Iron under Extreme Condition

A modified Ehrenfest formalism for efficient large-scale ab initio molecular dynamics

An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics

Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

Gaussian Statistics of the Cosmic Microwave Background: Correlation of Temperature Extrema In the COBE DMR Two-Year Sky Maps

A Concurrent Multiscale Method Based on Smoothed Molecular Dynamics for Large-Scale Parallel Computation at Finite Temperature

Effective Langevin Equation Approach to the Molecular Diffusion on Optical Lattices

Born-Oppenheimer Molecular Dynamics with a Linear Combination of Atomic Orbitals and Hybrid Functionals for Condensed Matter Simulations Made Possible. Theory and Performance for the Microcanonical and Canonical Ensembles

Contraction and Convergence Rates for Discretized Kinetic Langevin Dynamics

Global Langevin model of multidimensional biomolecular dynamics

Molecular dynamics in a grand ensemble: Bergmann–Lebowitz model and adaptive resolution simulation

Exploring A Multi-Scale Method for Molecular Simulations in Continuum Solvent Model: Explicit Simulation of Continuum Solvent As An Incompressible Fluid

Next generation extended Lagrangian first principles molecular dynamics

The velocity jump Langevin process and its splitting scheme: long time convergence and numerical accuracy