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:2011.08113 · The Journal of Chemical Physics

Accelerating convergence to the thermodynamic limit with twist angle selection applied to methods beyond many-body perturbation theory

Tina N. Mihm, William Z. Van Benschoten, James J. Shepherd

Published Jan 12, 2021• 8 citations
arXivHTMLPDF

Abstract

We recently developed a scheme to use low-cost calculations to find a single twist angle where the coupled cluster doubles energy of a single calculation matches the twist-averaged coupled cluster doubles energy in a finite unit cell. We used initiator full configuration interaction quantum Monte Carlo as an example of an exact method beyond coupled cluster doubles theory to show that this selected twist angle approach had comparable accuracy in methods beyond coupled cluster. Furthermore, at least for small system sizes, we show that the same twist angle can also be found by comparing the energy directly (at the level of second-order Moller-Plesset theory), suggesting a route toward twist angle selection, which requires minimal modification to existing codes that can perform twist averaging.

DOI
10.1063/5.0033408

Related papers

Power Laws Used to Extrapolate the Coupled Cluster Correlation Energy to the Thermodynamic Limit

Tina N. Mihm, Bingdi Yang, James J. Shepherd

Recent calculations using coupled cluster on solids have raised the discussion of using a N –1/3 power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used N –1 power law, which is, for example, often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of 922 electrons for an extremely wide range of densities (0.1 < r s < 100.0) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for the basis set incompleteness error and use a selected twist angle approach to mitigate the finite size error from shell filling effects. Analyzing these data, we initially find that a power law of N –1/3 appears to fit the data better than a N –1 power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that N –1 still applies to large system sizes and that the apparent N –1/3 power law occurs only at low N .

arXiv:2007.11696•Journal of Chemical Theory and Computation•Apr 8, 2021•14 citations

Single- and two-particle finite size effects in interacting lattice systems

Sergei Iskakov, Hanna Terletska, Emanuel Gull

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

arXiv:2208.08521•Physical review. B./Physical review. B•Dec 5, 2022•3 citations

A shortcut to the thermodynamic limit for quantum many-body calculations of metals

Tina N. Mihm, Tobias Schäfer, Sai Kumar Ramadugu +3 more

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.

arXiv:2103.06318•Nature Computational Science•Dec 16, 2021•34 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

Size Extensive Auxiliary-Field Quantum Monte Carlo with Perturbative Coupled Cluster Trial Wavefunction

Yichi Zhang, Ankit Mahajan, Yann Damour +1 more

In this work, we develop a size extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as $O(N^5)$ for local energy evaluation by treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions perturbatively. Comprehensive numerical examinations, spanning from main-group molecules to $3d$ transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles (CISD) trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas (UEG) provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit (TDL).

arXiv:2605.11257•May 11, 2026•0 citations

Size-extensive auxiliary-field quantum Monte Carlo with perturbative coupled cluster trial wavefunctions

Yichi Zhang, Ankit Mahajan, Yann Damour +1 more

In this work, we develop a size-extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5) for local energy evaluation by perturbatively treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions. Comprehensive numerical examinations, spanning from main-group molecules to 3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit.

The Journal of Chemical Physics•Sep 1, 2026•0 citations
We recently developed a scheme to use low-cost calculations to find a single twist angle where the coupled cluster doubles energy of a single calculation matches the twist-averaged coupled cluster doubles energy in a finite unit cell. We used initiator full configuration interaction quantum Monte Carlo as an example of an exact method beyond coupled cluster doubles theory to show that this selected twist angle approach had comparable accuracy in methods beyond coupled cluster. Furthermore, at least for small system sizes, we show that the same twist angle can also be found by comparing the energy directly (at the level of second-order Moller-Plesset theory), suggesting a route toward twist angle selection, which requires minimal modification to existing codes that can perform twist averaging.
Recent calculations using coupled cluster on solids have raised the discussion of using a N –1/3 power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used N –1 power law, which is, for example, often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of 922 electrons for an extremely wide range of densities (0.1 < r s < 100.0) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for the basis set incompleteness error and use a selected twist angle approach to mitigate the finite size error from shell filling effects. Analyzing these data, we initially find that a power law of N –1/3 appears to fit the data better than a N –1 power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that N –1 still applies to large system sizes and that the apparent N –1/3 power law occurs only at low N .
Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.
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.
In this work, we develop a size extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5)O(N^5)O(N5) for local energy evaluation by treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions perturbatively. Comprehensive numerical examinations, spanning from main-group molecules to 3d3d3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles (CISD) trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas (UEG) provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit (TDL).
In this work, we develop a size-extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5) for local energy evaluation by perturbatively treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions. Comprehensive numerical examinations, spanning from main-group molecules to 3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit.