Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0806.0957 · Physical Review B

Finite-size errors in continuum quantum Monte Carlo calculations

N. D. Drummond, Richard J. Needs, A. Sorouri, W. M. C. Foulkes

Published Sep 12, 2008• 233 citations
arXivHTMLPDF

Abstract

We analyze the problem of eliminating finite-size errors from quantum Monte Carlo (QMC) energy data. We demonstrate that both (i) adding a recently proposed [S. Chiesa et al., Phys. Rev. Lett. 97, 076404 (2006)] finite-size correction to the Ewald energy and (ii) using the model periodic Coulomb (MPC) interaction [L. M. Fraser et al., Phys. Rev. B 53, 1814 (1996); P. R. C. Kent et al., Phys. Rev. B 59, 1917 (1999); A. J. Williamson et al., Phys. Rev. B 55, R4851 (1997)] are good solutions to the problem of removing finite-size effects from the interaction energy in cubic systems provided the exchange-correlation (XC) hole has converged with respect to system size. However, we find that the MPC interaction distorts the XC hole in finite systems, implying that the Ewald interaction should be used to generate the configuration distribution. The finite-size correction of Chiesa et al. [Phys. Rev. Lett. 97, 076404 (2006)] is shown to be incomplete in systems of low symmetry. Beyond-leading-order corrections to the kinetic energy are found to be necessary at intermediate and high densities; we investigate the effect of adding such corrections to QMC data for the homogeneous electron gas. We analyze finite-size errors in two-dimensional systems and show that the leading-order behavior differs from that which has hitherto been supposed. We compare the efficiencies of different twist-averaging methods for reducing single-particle finite-size errors and we examine the performance of various finite-size extrapolation formulas. Finally, we investigate the system-size scaling of biases in diffusion QMC.

DOI
10.1103/physrevb.78.125106

Related papers

Finite-size errors in quantum many-body simulations of extended systems

Paul R. C. Kent, Randolph Q. Hood, Andrew Williamson +3 more

Further developments are introduced in the theory of finite-size errors in quantum many-body simulations of extended systems using periodic boundary conditions. We show that our recently introduced model periodic Coulomb interaction [A. J. Williamson et al., Phys. Rev. B 55, R4851 (1997)] can be applied consistently to all Coulomb interactions in the system. The model periodic Coulomb interaction greatly reduces the finite-size errors in quantum many-body simulations. We illustrate the practical application of our techniques with Hartree-Fock and variational and diffusion quantum Monte Carlo calculations for ground- and excited-state calculations. We demonstrate that the finite-size effects in electron promotion and electron addition/subtraction excitation energy calculations are very similar.

arXiv:cond-mat/9810243•Physical review. B, Condensed matter•Jan 15, 1999•113 citations

Finite-size effects and Coulomb interactions in quantum Monte Carlo calculations for homogeneous systems with periodic boundary conditions

Louisa M. Fraser, W. M. C. Foulkes, Gunaretnam Rajagopal +3 more

Quantum Monte Carlo (QMC) calculations are only possible in finite systems and so solids and liquids must be modeled using small simulation cells subject to periodic boundary conditions. The resulting finite-size errors are often corrected using data from local-density functional or Hartree-Fock calculations, but systematic errors remain after these corrections have been applied. The results of our jellium QMC calculations for simulation cells containing more than 600 electrons confirm that the residual errors are significant and decay very slowly as the system size increases. We show that they are sensitive to the form of the model Coulomb interaction used in the simulation cell Hamiltonian and that the usual choice, exemplified by the Ewald summation technique, is not the best. The finite-size errors can be greatly reduced and the speed of the calculations increased by a factor of 20 if a better choice is made. Finite-size effects plague most methods used for extended Coulomb systems and many of the ideas in this paper are quite general: they may be applied to any type of quantum or classical Monte Carlo simulation, to other many-body approaches such as the GW method, and to Hartree-Fock and density-functional calculations.

Physical review. B, Condensed matter•Jan 15, 1996•310 citations

Good-$Γ$ Numbers: A Shell-Optimized Route to the Thermodynamic Limit in Quantum Monte Carlo Calculations of the High-Density Two-Dimensional Electron Gas

S. Azadi, T. D. Kühne, A. Principi +1 more

Single-particle finite-size errors remain a central obstacle in quantum Monte Carlo (QMC) simulations of periodic systems, conventionally reduced via twist averaging (TA). Here we introduce Good-$Γ$ (GG) particle numbers as system sizes $N$ at which the $Γ$-point Hartree--Fock (HF) energy of the high-density two-dimensional uniform electron gas (2D UEG) lies anomalously close to its thermodynamic-limit (TL) value. All results in this work concern the 2D UEG in its high-density Fermi-liquid regime, $r_s=1$--$4$. We systematically investigate the kinetic- and exchange-energy errors, $ΔKE(N)$ and $ΔEX(N)$, that underlie this one-body finite-size-error cancellation across four densities and both spin polarizations. We then compare diffusion Monte Carlo (DMC) energies extrapolated to the TL from GG-$N$ sequences against twist-averaged extrapolations at the same densities and polarizations, finding agreement to within a fraction of a millihartree in every case, independent of whether a $1/N$ or an $N^{-5/4}$ extrapolation form is used.

arXiv:2610.27761•Oct 20, 2026•0 citations

Quantum Monte Carlo calculations of energy gaps from first principles

R. J. Hunt, Marcin Szyniszewski, Genki I. Prayogo +2 more

We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.

arXiv:1806.04750•Physical review. B./Physical review. B•Aug 13, 2018•50 citations

Exact special twist method for quantum Monte Carlo simulations

Mario Dagrada, Seher Karakuzu, Verónica Laura Vildosola +2 more

We present a systematic investigation of the special twist method introduced by Rajagopal et al. [Phys. Rev. B 51, 10591 (1995)] for reducing finite-size effects in correlated calculations of periodic extended systems with Coulomb interactions and Fermi statistics. We propose a procedure for finding special twist values which, at variance with previous applications of this method, reproduce the energy of the mean-field infinite-size limit solution within an adjustable (arbitrarily small) numerical error. This choice of the special twist is shown to be the most accurate single-twist solution for curing one-body finite-size effects in correlated calculations. For these reasons we dubbed our procedure ``exact special twist'' (EST). EST only needs a fully converged independent-particles or mean-field calculation within the primitive cell and a simple fit to find the special twist along a specific direction in the Brillouin zone. We first assess the performances of EST in a simple correlated model such as the three-dimensional electron gas. Afterwards, we test its efficiency within ab initio quantum Monte Carlo simulations of metallic elements of increasing complexity. We show that EST displays an overall good performance in reducing finite-size errors comparable to the widely used twist average technique but at a much lower computational cost since it involves the evaluation of just one wave function. We also demonstrate that the EST method shows similar performances in the calculation of correlation functions, such as the ionic forces for structural relaxation and the pair radial distribution function in liquid hydrogen. Our conclusions point to the usefulness of EST for correlated supercell calculations; our method will be particularly relevant when the physical problem under consideration requires large periodic cells.

arXiv:1606.06205•Physical review. B./Physical review. B•Dec 5, 2016•28 citations

Elimination of Coulomb finite-size effects in quantum many-body simulations

Andrew Williamson, Gunaretnam Rajagopal, Richard J. Needs +4 more

A model interaction is introduced for quantum many-body simulations of Coulomb systems using periodic boundary conditions. The interaction gives much smaller finite size effects than the standard Ewald interaction and is also much faster to compute. Variational quantum Monte Carlo simulations of diamond-structure silicon with up to 1000 electrons demonstrate the effectiveness of our method.

Physical review. B, Condensed matter•Feb 15, 1997•111 citations
We analyze the problem of eliminating finite-size errors from quantum Monte Carlo (QMC) energy data. We demonstrate that both (i) adding a recently proposed [S. Chiesa et al., Phys. Rev. Lett. 97, 076404 (2006)] finite-size correction to the Ewald energy and (ii) using the model periodic Coulomb (MPC) interaction [L. M. Fraser et al., Phys. Rev. B 53, 1814 (1996); P. R. C. Kent et al., Phys. Rev. B 59, 1917 (1999); A. J. Williamson et al., Phys. Rev. B 55, R4851 (1997)] are good solutions to the problem of removing finite-size effects from the interaction energy in cubic systems provided the exchange-correlation (XC) hole has converged with respect to system size. However, we find that the MPC interaction distorts the XC hole in finite systems, implying that the Ewald interaction should be used to generate the configuration distribution. The finite-size correction of Chiesa et al. [Phys. Rev. Lett. 97, 076404 (2006)] is shown to be incomplete in systems of low symmetry. Beyond-leading-order corrections to the kinetic energy are found to be necessary at intermediate and high densities; we investigate the effect of adding such corrections to QMC data for the homogeneous electron gas. We analyze finite-size errors in two-dimensional systems and show that the leading-order behavior differs from that which has hitherto been supposed. We compare the efficiencies of different twist-averaging methods for reducing single-particle finite-size errors and we examine the performance of various finite-size extrapolation formulas. Finally, we investigate the system-size scaling of biases in diffusion QMC.
Further developments are introduced in the theory of finite-size errors in quantum many-body simulations of extended systems using periodic boundary conditions. We show that our recently introduced model periodic Coulomb interaction [A. J. Williamson et al., Phys. Rev. B 55, R4851 (1997)] can be applied consistently to all Coulomb interactions in the system. The model periodic Coulomb interaction greatly reduces the finite-size errors in quantum many-body simulations. We illustrate the practical application of our techniques with Hartree-Fock and variational and diffusion quantum Monte Carlo calculations for ground- and excited-state calculations. We demonstrate that the finite-size effects in electron promotion and electron addition/subtraction excitation energy calculations are very similar.
Quantum Monte Carlo (QMC) calculations are only possible in finite systems and so solids and liquids must be modeled using small simulation cells subject to periodic boundary conditions. The resulting finite-size errors are often corrected using data from local-density functional or Hartree-Fock calculations, but systematic errors remain after these corrections have been applied. The results of our jellium QMC calculations for simulation cells containing more than 600 electrons confirm that the residual errors are significant and decay very slowly as the system size increases. We show that they are sensitive to the form of the model Coulomb interaction used in the simulation cell Hamiltonian and that the usual choice, exemplified by the Ewald summation technique, is not the best. The finite-size errors can be greatly reduced and the speed of the calculations increased by a factor of 20 if a better choice is made. Finite-size effects plague most methods used for extended Coulomb systems and many of the ideas in this paper are quite general: they may be applied to any type of quantum or classical Monte Carlo simulation, to other many-body approaches such as the GW method, and to Hartree-Fock and density-functional calculations.
Good-ΓΓΓ Numbers: A Shell-Optimized Route to the Thermodynamic Limit in Quantum Monte Carlo Calculations of the High-Density Two-Dimensional Electron Gas
Single-particle finite-size errors remain a central obstacle in quantum Monte Carlo (QMC) simulations of periodic systems, conventionally reduced via twist averaging (TA). Here we introduce Good-ΓΓΓ (GG) particle numbers as system sizes NNN at which the ΓΓΓ-point Hartree--Fock (HF) energy of the high-density two-dimensional uniform electron gas (2D UEG) lies anomalously close to its thermodynamic-limit (TL) value. All results in this work concern the 2D UEG in its high-density Fermi-liquid regime, rs=1r_s=1rs​=1--444. We systematically investigate the kinetic- and exchange-energy errors, ΔKE(N)ΔKE(N)ΔKE(N) and ΔEX(N)ΔEX(N)ΔEX(N), that underlie this one-body finite-size-error cancellation across four densities and both spin polarizations. We then compare diffusion Monte Carlo (DMC) energies extrapolated to the TL from GG-NNN sequences against twist-averaged extrapolations at the same densities and polarizations, finding agreement to within a fraction of a millihartree in every case, independent of whether a 1/N1/N1/N or an N−5/4N^{-5/4}N−5/4 extrapolation form is used.
We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.
We present a systematic investigation of the special twist method introduced by Rajagopal et al. [Phys. Rev. B 51, 10591 (1995)] for reducing finite-size effects in correlated calculations of periodic extended systems with Coulomb interactions and Fermi statistics. We propose a procedure for finding special twist values which, at variance with previous applications of this method, reproduce the energy of the mean-field infinite-size limit solution within an adjustable (arbitrarily small) numerical error. This choice of the special twist is shown to be the most accurate single-twist solution for curing one-body finite-size effects in correlated calculations. For these reasons we dubbed our procedure ``exact special twist'' (EST). EST only needs a fully converged independent-particles or mean-field calculation within the primitive cell and a simple fit to find the special twist along a specific direction in the Brillouin zone. We first assess the performances of EST in a simple correlated model such as the three-dimensional electron gas. Afterwards, we test its efficiency within ab initio quantum Monte Carlo simulations of metallic elements of increasing complexity. We show that EST displays an overall good performance in reducing finite-size errors comparable to the widely used twist average technique but at a much lower computational cost since it involves the evaluation of just one wave function. We also demonstrate that the EST method shows similar performances in the calculation of correlation functions, such as the ionic forces for structural relaxation and the pair radial distribution function in liquid hydrogen. Our conclusions point to the usefulness of EST for correlated supercell calculations; our method will be particularly relevant when the physical problem under consideration requires large periodic cells.