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:2304.03330 · Physical Review X

Inverse Volume Scaling of Finite-Size Error in Periodic Coupled Cluster Theory

Xin Xing, Lin Lin

Published Mar 28, 2024• 13 citations
arXivHTMLPDF

Abstract

Coupled cluster theory is one of the most popular post-Hartree-Fock methods for molecular quantum chemistry. The finite-size error of the correlation energy in periodic coupled cluster calculations for three-dimensional insulating systems has been observed to satisfy the inverse volume scaling, even in the absence of any correction schemes. This is surprising, as simpler theories that utilize only a subset of the coupled cluster diagrams exhibit much slower decay of the finite-size error, which scales inversely with the length of the system. In this study, we review the current understanding of finite-size error in quantum chemistry methods for periodic systems. We introduce new tools that elucidate the mechanisms behind this phenomenon in the context of coupled cluster doubles calculations. This reconciles some seemingly paradoxical statements related to finite-size scaling. Our findings also show that singularity subtraction can be a powerful method to effectively reduce finite-size errors in practical quantum chemistry calculations for periodic systems. Published by the American Physical Society 2024

DOI
10.1103/physrevx.14.011059

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

Lowering of the complexity of quantum chemistry methods by choice of representation

Narbe Mardirossian, James D. McClain, Garnet Kin‐Lic Chan

The complexity of the standard hierarchy of quantum chemistry methods is not invariant to the choice of representation. This work explores how the scaling of common quantum chemistry methods can be reduced using real-space, momentum-space, and time-dependent intermediate representations without introducing approximations. We find the scalings of exact Gaussian basis Hartree–Fock theory, second-order Møller-Plesset perturbation theory, and coupled cluster theory (specifically, linearized coupled cluster doubles and the distinguishable cluster approximation with doubles) to be O(N3), O(N3), and O(N5), respectively, where N denotes the system size. These scalings are not asymptotic and hold over all ranges of N.

arXiv:1710.01004•The Journal of Chemical Physics•Jan 23, 2018•32 citations

Coupled Cluster Theory in Materials Science

Igor Ying Zhang, Andreas Grüneis

The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.

arXiv:2004.06424•Frontiers in Materials•Jun 5, 2019•138 citations

How the Exchange Energy Can Affect the Power Laws Used to Extrapolate the Coupled Cluster Correlation Energy to the Thermodynamic Limit

Tina N. Mihm, Laura Weiler, James J. Shepherd

High Resolution Image Download MS PowerPoint Slide Finite size error is commonly removed from coupled cluster theory calculations by N –1 extrapolations over correlation energy calculations of different system sizes ( N ), where the N –1 scaling comes from the total energy rather than the correlation energy. However, previous studies in the quantum Monte Carlo community suggest an exchange-energy-like power law of N –2/3 should also be present in the correlation energy when using the conventional Coulomb interaction. The rationale for this is that the total energy goes as N –1 and the exchange energy goes as N –2/3; thus, the correlation energy should be a combination of these two power laws. Further, in coupled cluster theory, these power laws are related to the low G scaling of the transition structure factor, S ( G ), which is a property of the coupled cluster wave function calculated from the amplitudes. We show here that data from coupled cluster doubles calculations on the uniform electron gas fit a function with a low G behavior of S ( G ) ∼ G . The prefactor for this linear term is derived from the exchange energy to be consistent with an N –2/3 power law at large N . Incorporating the exchange structure factor into the transition structure factor results in a combined structure factor of S ( G ) ∼ G 2, consistent with an N –1 scaling of the exchange-correlation energy. We then look for the presence of an N –2/3 power law in the energy. To do so, we first develop a plane-wave cutoff scheme with less noise than the traditional basis set used for the uniform electron gas. Then, we collect data from a wide range of electron numbers and densities to systematically test five methods using N –1 scaling, N –2/3 scaling, or combinations of both scaling behaviors. We find that power laws that incorporate both N –1 and N –2/3 scaling perform better than either alone, especially when the prefactor for N –2/3 scaling can be found from exchange energy calculations.

arXiv:2302.05051•Journal of Chemical Theory and Computation•Mar 14, 2023•11 citations

Coupled cluster method tailored with quantum computing

Luca Erhart, Yuichiro Yoshida, Viktor E. Khinevich +1 more

Introducing an active space approximation is inevitable for the quantum computations of chemical systems. However, this approximation ignores the electron correlations related to nonactive orbitals. Here, we propose a computational method for correcting quantum computing results using a well-established classical theory called coupled cluster theory. Our approach efficiently extracts the quantum state from a quantum device by computational basis tomography. The extracted expansion coefficients of the quantum state are embedded into the coupled cluster ansatz within the framework of the tailored coupled cluster method. We demonstrate the performance of our method by verifying the potential energy curves of LiH, H 2 O , and N 2 with a correlation-energy correction scheme. Our method demonstrates reasonable potential energy curves even when the standard coupled cluster fails. The sufficient numbers of measurements for tomography were also investigated. Furthermore, this method successfully estimated the activation energy of the Cope rearrangement reaction of 1,5-hexadiene together with perturbative triples correction. These demonstrations suggest that our approach has the potential for practical quantum chemical calculations using quantum computers. Published by the American Physical Society 2024

arXiv:2312.11012•Physical Review Research•Jun 3, 2024•8 citations

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
Coupled cluster theory is one of the most popular post-Hartree-Fock methods for molecular quantum chemistry. The finite-size error of the correlation energy in periodic coupled cluster calculations for three-dimensional insulating systems has been observed to satisfy the inverse volume scaling, even in the absence of any correction schemes. This is surprising, as simpler theories that utilize only a subset of the coupled cluster diagrams exhibit much slower decay of the finite-size error, which scales inversely with the length of the system. In this study, we review the current understanding of finite-size error in quantum chemistry methods for periodic systems. We introduce new tools that elucidate the mechanisms behind this phenomenon in the context of coupled cluster doubles calculations. This reconciles some seemingly paradoxical statements related to finite-size scaling. Our findings also show that singularity subtraction can be a powerful method to effectively reduce finite-size errors in practical quantum chemistry calculations for periodic systems. Published by the American Physical Society 2024
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.
The complexity of the standard hierarchy of quantum chemistry methods is not invariant to the choice of representation. This work explores how the scaling of common quantum chemistry methods can be reduced using real-space, momentum-space, and time-dependent intermediate representations without introducing approximations. We find the scalings of exact Gaussian basis Hartree–Fock theory, second-order Møller-Plesset perturbation theory, and coupled cluster theory (specifically, linearized coupled cluster doubles and the distinguishable cluster approximation with doubles) to be O(N3), O(N3), and O(N5), respectively, where N denotes the system size. These scalings are not asymptotic and hold over all ranges of N.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
High Resolution Image Download MS PowerPoint Slide Finite size error is commonly removed from coupled cluster theory calculations by N –1 extrapolations over correlation energy calculations of different system sizes ( N ), where the N –1 scaling comes from the total energy rather than the correlation energy. However, previous studies in the quantum Monte Carlo community suggest an exchange-energy-like power law of N –2/3 should also be present in the correlation energy when using the conventional Coulomb interaction. The rationale for this is that the total energy goes as N –1 and the exchange energy goes as N –2/3; thus, the correlation energy should be a combination of these two power laws. Further, in coupled cluster theory, these power laws are related to the low G scaling of the transition structure factor, S ( G ), which is a property of the coupled cluster wave function calculated from the amplitudes. We show here that data from coupled cluster doubles calculations on the uniform electron gas fit a function with a low G behavior of S ( G ) ∼ G . The prefactor for this linear term is derived from the exchange energy to be consistent with an N –2/3 power law at large N . Incorporating the exchange structure factor into the transition structure factor results in a combined structure factor of S ( G ) ∼ G 2, consistent with an N –1 scaling of the exchange-correlation energy. We then look for the presence of an N –2/3 power law in the energy. To do so, we first develop a plane-wave cutoff scheme with less noise than the traditional basis set used for the uniform electron gas. Then, we collect data from a wide range of electron numbers and densities to systematically test five methods using N –1 scaling, N –2/3 scaling, or combinations of both scaling behaviors. We find that power laws that incorporate both N –1 and N –2/3 scaling perform better than either alone, especially when the prefactor for N –2/3 scaling can be found from exchange energy calculations.
Introducing an active space approximation is inevitable for the quantum computations of chemical systems. However, this approximation ignores the electron correlations related to nonactive orbitals. Here, we propose a computational method for correcting quantum computing results using a well-established classical theory called coupled cluster theory. Our approach efficiently extracts the quantum state from a quantum device by computational basis tomography. The extracted expansion coefficients of the quantum state are embedded into the coupled cluster ansatz within the framework of the tailored coupled cluster method. We demonstrate the performance of our method by verifying the potential energy curves of LiH, H 2 O , and N 2 with a correlation-energy correction scheme. Our method demonstrates reasonable potential energy curves even when the standard coupled cluster fails. The sufficient numbers of measurements for tomography were also investigated. Furthermore, this method successfully estimated the activation energy of the Cope rearrangement reaction of 1,5-hexadiene together with perturbative triples correction. These demonstrations suggest that our approach has the potential for practical quantum chemical calculations using quantum computers. Published by the American Physical Society 2024
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 .