Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory
- Univ. of California, Berkeley, CA (United States)
- Beijing Normal University, Beijing (China)
- Univ. of California, Berkeley, CA (United States); Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
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 [Math. Comp. 30 (1976), pp. 1–23] 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 (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, 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.
- Research Organization:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC); US Air Force Office of Scientific Research (AFOSR)
- Grant/Contract Number:
- SC0022198; SC0017867; FA9550-18-1-0095; 201906040071
- OSTI ID:
- 2281877
- Journal Information:
- Mathematics of Computation, Vol. 93, Issue 346; ISSN 0025-5718
- Publisher:
- American Mathematical SocietyCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Assessing electronic structure approaches for gas-ligand interactions in metal-organic frameworks: The CO{sub 2}-benzene complex
Finite-temperature second-order many-body perturbation and Hartree–Fock theories for one-dimensional solids: An application to Peierls and charge-density-wave transitions in conjugated polymers