Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory
arXiv:2108.00206 · doi:10.1090/mcom/3877
Abstract
Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of [Lyness, 1976] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations [Xing, Li, Lin, J. Chem. Theory Comput. 2021]. Our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.
References in corpus (11)
- Exponential localization of Wannier functions in insulators
- The Finite Size Error in Many-body Simulations with long-Ranged Interactions
- Ideal regularization of the Coulomb singularity in exact exchange by Wigner-Seitz truncated interactions: towards chemical accuracy in non-trivial systems
- Electrostatics in Periodic Boundary Conditions and Real-space Corrections
- Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit
- Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller-Plesset Perturbation Theory
- Quartic scaling MP2 for solids: A highly parallelized algorithm in the plane wave basis
- A simple, efficient, and general treatment of the singularities in Hartree-Fock and exact-exchange Kohn-Sham methods for solids
- An optimized twist angle to find the twist-averaged correlation energy applied to the uniform electron gas
- Staggered mesh method for correlation energy calculations of solids: Random phase approximation in direct ring coupled cluster doubles and adiabatic connection formalisms
- Staggered mesh method for correlation energy calculations of solids: Second order Møller-Plesset perturbation theory