Implementation of the hybrid exchange-correlation functionals in the SIESTA code

Abstract

We present an efficient and accurate implementation of hybrid exchange-correlation (XC) functionals in the SIESTA code, enabling large-scale simulations based on Hartree-Fock-type exact exchange combined with strictly localized numerical atomic orbitals (NAOs). Our approach exploits a fitted representation of the NAOs in terms of Gaussian-type orbitals (GTOs), which allows for the analytical evaluation of four-center electron repulsion integrals (ERIs) via the LIBINT library. This framework is seamlessly integrated with SIESTA's real-space grid and sparse-matrix infrastructure, and is combined with multiple screening techniques to control the computational complexity. We also introduce a fully analytical formulation of hybrid-functional forces and a dynamic parallel distribution scheme that ensures excellent scalability. We validate our implementation through benchmark calculations on a broad set of systems (including semiconductors, insulators, and two-dimensional materials) and demonstrate that the HSE06 functional significantly improves the prediction of band gaps compared to PBE, in close agreement with G0W0 and experimental data. We analyze in detail the trade-offs between accuracy and computational efficiency as a function of the number of Gaussians, basis set range, and integral screening thresholds. Our results confirm that hybrid functional calculations in SIESTA are now feasible for large extended systems, making accurate first-principles predictions of electronic and structural properties accessible at scale.

27 pages, 5 figures, 5 tables

References in corpus (30)

• Restoring the density-gradient expansion for exchange in solids and surfaces
• Generalized gradient approximation for solids and their surfaces
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Electronic transport in two dimensional graphene
• Optimized norm-conserving Vanderbilt pseudopotentials
• The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table
• Hexagonal boron nitride is an indirect bandgap semiconductor
• Localization and delocalization errors in density functional theory and implications for band-gap prediction
• Recent developments in the PySCF program package
• Numerical atomic orbitals for linear scaling
• 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
• The GW compendium: A practical guide to theoretical photoemission spectroscopy
• Polaronic hole localization and multiple hole binding of acceptors in oxide wide-gap semiconductors
• The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water
• Thermal physics of the lead chalcogenides PbS, PbSe, and PbTe from first principles
• Toward a realistic description of multilayer black phosphorus: from $GW$ approximation to large-scale tight-binding simulations
• Order-N implementation of exact exchange in extended systems
• Adaptively Compressed Exchange Operator
• Large scale and linear scaling DFT with the CONQUEST code
• Efficient many-body calculations of 2D materials using exact limits for the screened potential: Band gaps of MoS${}_2$, hBN, and phosphorene
• The PSML format and library for norm-conserving pseudopotential data curation and interoperability
• Hybrid functional pseudopotentials
• Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms
• The CECAM Electronic Structure Library and the modular software development paradigm
• The dynamic parallel distribution algorithm for hybrid density-functional calculations in HONPAS package
• Benchmarking exchange-correlation potentials with the mstar60 dataset: Importance of the nonlocal exchange potential for effective mass calculations in semiconductors
• The reliability of hybrid functionals for accurate fundamental and optical gap prediction of bulk solids and surfaces
• The static parallel distribution algorithms for hybrid density-functional calculations in HONPAS package
• Predicting Fundamental Gaps of Chromium-Based 2D Materials Using GW Methods