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

Abstract

Hybrid density functional approximations (DFAs) offer compelling accuracy for ab initio electronic-structure simulations of molecules, nanosystems, and bulk materials, addressing some deficiencies of computationally cheaper, frequently used semilocal DFAs. However, the computational bottleneck of hybrid DFAs is the evaluation of the non-local exact exchange contribution, which is the limiting factor for the application of the method for large-scale simulations. In this work, we present a drastically optimized resolution-of-identity-based real-space implementation of the exact exchange evaluation for both non-periodic and periodic boundary conditions in the all-electron code FHI-aims, targeting high-performance CPU compute clusters. The introduction of several new refined Message Passing Interface (MPI) parallelization layers and shared memory arrays according to the MPI-3 standard were the key components of the optimization. We demonstrate significant improvements of memory and performance efficiency, scalability, and workload distribution, extending the reach of hybrid DFAs to simulation sizes beyond ten thousand atoms. As a necessary byproduct of this work, other code parts in FHI-aims have been optimized as well, e.g., the computation of the Hartree potential and the evaluation of the force and stress components. We benchmark the performance and scaling of the hybrid DFA based simulations for a broad range of chemical systems, including hybrid organic-inorganic perovskites, organic crystals and ice crystals with up to 30,576 atoms (101,920 electrons described by 244,608 basis functions).

References in corpus (8)

• Predicting Band Gaps with Hybrid Density Functionals
• Order-N implementation of exact exchange in extended systems
• Roadmap on Electronic Structure Codes in the Exascale Era
• All-electron periodic $G_0{W}_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks
• Shortcomings of meta-GGA functionals when describing magnetism
• The random phase approximation applied to ice
• SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization
• High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach

Cited by in corpus (6)

• Implementation of the hybrid exchange-correlation functionals in the SIESTA code
• An atomic-scale perspective on individual thiol-terminated molecules anchored to single S vacancies in MoS${}_2$
• Efficient Computation of the Long-Range Exact Exchange using an Extended Screening Function
• Atomic Confinement Potentials and the Generation of Numerical Atomic Orbitals
• Low-scaling \textit{GW} calculations of quasi-particle energies for extended systems within the numerical atomic orbital framework
• Consistent GMTKN55 and molecular-crystal accuracy using minimally empirical DFT with XDM(Z) dispersion