Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory
From MaRDI portal
Publication:6140844
Perturbations in context of PDEs (35B20) Discrete version of topics in analysis (39A12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15) Yang-Mills and other gauge theories in mechanics of particles and systems (70S15) Computational stability and error-correcting codes for quantum computation and communication processing (81P73) Fixed points and periodic points of finite-dimensional Hamiltonian and Lagrangian systems (37J12)
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 rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order M{o}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 method 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, 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.
Recommendations
- Finite-size effects in periodic coupled cluster calculations
- Error bars and quadratically convergent methods for the numerical simulation of the Hartree-Fock equations
- Discretization error cancellation in the plane-wave approximation of periodic Hamiltonians with Coulomb singularities
- Discretization error cancellation in electronic structure calculation: toward a quantitative study
- \(h-P\) finite element approximation for full-potential electronic structure calculations
Cites work
- A scalable and accurate algorithm for the computation of Hartree-Fock exchange
- A trapezoidal rule error bound unifying the Euler–Maclaurin formula and geometric convergence for periodic functions
- An Error Functional Expansion for N-Dimensional Quadrature with an Integrand Function Singular at a Point
- Asymptotic Expansions for Integration Formulas in One or More Dimensions
- Correlation Energy of an Electron Gas at High Density
- Electronic Structure
- Numerical quadrature in the Brillouin zone for periodic Schrödinger operators
- On the remainder term in the N-dimensional Euler Maclaurin expansion
- Optimal decay of Wannier functions in Chern and quantum Hall insulators
- The exponentially convergent trapezoidal rule
This page was built for publication: Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6140844)
