Ideal regularization of the Coulomb singularity in exact exchange by Wigner-Seitz truncated interactions: towards chemical accuracy in non-trivial systems

Abstract

Hybrid density functionals show great promise for chemically-accurate first principles calculations, but their high computational cost limits their application in non-trivial studies, such as exploration of reaction pathways of adsorbents on periodic surfaces. One factor responsible for their increased cost is the dense Brillouin-zone sampling necessary to accurately resolve an integrable singularity in the exact exchange energy. We analyze this singularity within an intuitive formalism based on Wannier-function localization and analytically prove Wigner-Seitz truncation to be the ideal method for regularizing the Coulomb potential in the exchange kernel. We show that this method is limited only by Brillouin-zone discretization errors in the Kohn-Sham orbitals, and hence converges the exchange energy exponentially with the number of k-points used to sample the Brillouin zone for all but zero-temperature metallic systems. To facilitate the implementation of this method, we develop a general construction for the plane-wave Coulomb kernel truncated on the Wigner-Seitz cell in one, two or three lattice directions. We compare several regularization methods for the exchange kernel in a variety of real systems including low-symmetry crystals and low-dimensional materials. We find that our Wigner-Seitz truncation systematically yields the best k-point convergence for the exchange energy of all these systems and delivers an accuracy to hybrid functionals comparable to semi-local and screened-exchange functionals at identical k-point sets.

14 pages, 9 figures

References in corpus (3)

• Truncation of Periodic Image Interactions for Confined Systems
• Order-N implementation of exact exchange in extended systems
• A simple, efficient, and general treatment of the singularities in Hartree-Fock and exact-exchange Kohn-Sham methods for solids

Cited by in corpus (15)

• Computational 2D Materials Database: Electronic Structure of Transition-Metal Dichalcogenides and Oxides
• Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry
• JDFTx: software for joint density-functional theory
• The charge-asymmetric nonlocally-determined local-electric (CANDLE) solvation model
• Evaluating continuum solvation models for the electrode-electrolyte interface: challenges and strategies for improvement
• First-principles electrostatic potentials for reliable alignment at interfaces and defects
• All-electron periodic $G_0{W}_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks
• Plasmons on the edge of MoS2 nanostructures
• Electrochemical Capacitance of CO-terminated Pt(111) is Dominated by CO-Solvent Gap
• Fast periodic Gaussian density fitting by range separation
• Interfacial water asymmetry at ideal electrochemical interfaces
• Broken-symmetry self-consistent GW approach: degree of spin contamination and evaluation of effective exchange couplings in solid antiferromagnets
• Gap control in phosphorene/BN structures from first principles calculations
• Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2
• Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions