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 (19)

• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• Predicting Band Gaps with Hybrid Density Functionals
• Order-N implementation of exact exchange in extended systems
• Accurate and efficient linear scaling DFT calculations with universal applicability
• Adaptively Compressed Exchange Operator
• Tunable Semiconductors: Control over Carrier States and Excitations in Layered Hybrid Organic-Inorganic Perovskites
• Large scale and linear scaling DFT with the CONQUEST code
• 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
• Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance
• Shortcomings of meta-GGA functionals when describing magnetism
• GPU-Acceleration of the ELPA2 Distributed Eigensolver for Dense Symmetric and Hermitian Eigenproblems
• The random phase approximation applied to ice
• Thermodynamic equilibrium conditions of graphene films on SiC
• Optimizations of the Eigensolvers in the ELPA Library
• Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library
• SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization
• A Koopmans-compliant screened exchange potential with correct asymptotic behavior for semiconductors
• 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