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:2501.18667 · Physical review. B./Physical review. B

Exploring the accuracy of the equation-of-motion coupled-cluster band gap of solids

Evgeny Moerman, Henrique Pereira Coutada Miranda, Alejandro Gallo, Andreas Irmler, Tobias Schäfer, Felix Hummel, Manuel Engel, Georg Kresse, Matthias Scheffler, Andreas Grüneis

Published Mar 26, 2025• 9 citations
arXivHTMLPDF

Abstract

While the periodic equation-of-motion coupled-cluster (EOM-CC) method promises systematic improvement of electronic band gap calculations in solids, its practical application at the singles and doubles level (EOM-CCSD) is hindered by severe finite-size errors in feasible simulation cells. We present a hybrid approach combining EOM-CCSD with the computationally less demanding G W approximation to estimate thermodynamic limit band gaps for several insulators and semiconductors. Our method substantially reduces required cell sizes while maintaining accuracy. Comparisons with experimental gaps and self-consistent G W calculations reveal that deviations in EOM-CCSD predictions correlate with reduced single excitation character of the excited many-electron states. Our work not only provides a computationally tractable approach to EOM-CC calculations in solids but also reveals fundamental insights into the role of single excitations in electronic-structure theory.

DOI
10.1103/physrevb.111.l121202

Related papers

Resolving finite-size errors in EOM-CCSD band gaps of solids with interacting-bath dynamical embedding theory

Jiachen Li, Christopher Hillenbrand, Christian Venturella +2 more

Periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) has shown promise for quantitative calculations of band structures in solids. However, its steep computational scaling has limited calculations to relatively coarse k-point meshes, leading to sizable finite-size errors and discrepant estimates of thermodynamic-limit band gaps in recent benchmarks. In this work, we revisit EOM-CCSD band gaps for ten semiconductors and insulators using interacting-bath dynamical embedding theory (ibDET), a systematically improvable Green's function embedding framework that enables dense Brillouin-zone sampling at modest computational cost. By pushing the k-point sampling up to 10 × 10 × 10, well beyond the system sizes accessible in canonical periodic EOM-CCSD calculations, we significantly reduce finite-size errors and obtain stable thermodynamic-limit extrapolations. We further compare G0W0@PBE, G0W0@HF, and EOM-CCSD on an equal footing using the same numerical settings in PySCF. We find that EOM-CCSD yields a mean absolute error of 0.32 eV relative to experimental band gaps for a test set of ten semiconductors and insulators, which is lower than that of G0W0@PBE. For ZnO, EOM-CCSD also accurately describes the Zn 3d-band binding energy, despite overestimating the band gap. These results demonstrate that ibDET offers a practical route to high-accuracy many-body electronic structure calculations in periodic systems.

arXiv:2606.12621•The Journal of Chemical Physics•Sep 1, 2026•0 citations

Resolving molecular biexcited states using classical leading-order triples corrections to quantum equation-of-motion UCCSD

Zachary W. Windom, Prince Frederick Kwao, Luke W. Bertels +3 more

Excited electronic states are difficult to accurately and tractably model, particularly those exhibiting multiexcitonic character. Scalable single-reference excited state methods typically fail to describe such states as they suffer from the defects of the mean-field approximation and/or neglect higher-rank excitation operators. To address this, we derive the [T] perturbative correction to the quantum self-consistent equation-of-motion unitary coupled cluster singles and doubles (q-sc-EOM-UCCSD) method. The method captures triple excitation effects through a classical step of second-order in many-body perturbation theory following a q-sc-EOM-UCCSD calculation on a quantum computer. The perturbative correction developed has a $\mathcal{O}(N^7)$ classical computational time scaling. We incorporate [T] into a hybrid compute strategy wherein eigenvectors of the effective Hamiltonian from the hybrid quantum-classical q-sc-EOM-UCCSD algorithm are relayed to a classical computer which performs the postprocessing required for the [T] correction. The benefits of [T] are quantified in a benchmark covering a variety of electronically excited states of the isoelectronic CH$^+$ and BH molecules. We find that the addition of the [T] correction to both trotterized and full operator variants of q-sc-EOM-CCSD can offer dramatic improvements over baseline q-sc-EOM-UCCSD. Our assessment clearly demonstrates the importance of considering higher-rank excitation operators in q-sc-EOM-UCCSD, quantifies the overall success of the [T] correction, and discusses some of its limitations.

arXiv:2610.09260•Oct 7, 2026•0 citations

Periodic Coupled-Cluster Green’s Function for Photoemission Spectra of Realistic Solids

Katelyn Laughon, Jason M. Yu, Tianyu Zhu

We present an efficient implementation of the coupled-cluster Green’s function (CCGF) method for simulating photoemission spectra of periodic systems. We formulate the periodic CCGF approach with Brillouin zone sampling in the Gaussian basis at the coupled-cluster singles and doubles (CCSD) level. To enable CCGF calculations of realistic solids, we propose an active-space self-energy correction scheme by combining CCGF with the cheaper many-body perturbation theory ( GW ) and implement the model order reduction (MOR) frequency interpolation technique. We find that the active-space self-energy correction and MOR techniques significantly reduce the computational cost of CCGF while maintaining the high accuracy. We apply the developed CCGF approaches to compute spectral properties and band structure of silicon (Si) and zinc oxide (ZnO) crystals using triple-ζ Gaussian basis sets and medium-size k -point sampling and find good agreement with experimental measurements.

arXiv:2208.07478•The Journal of Physical Chemistry Letters•Sep 26, 2022•15 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

N-representability violations in truncated equation-of-motion coupled-cluster methods

Stephen H. Yuwono, A. Eugene DePrince

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. These observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

arXiv:2305.06457•The Journal of Chemical Physics•Aug 7, 2023•3 citations

Accelerating Periodic Coupled Cluster and Algebraic Diagrammatic Construction Theories With Frozen Natural Orbitals

Ning-Yuan Chen, James D. Serna, Alexander Yu. Sokolov

We present a unified frozen natural orbital (FNO) framework for periodic Gaussian-orbital post-Hartree-Fock calculations with explicit k-point sampling. The approach includes conventional ground-state FNOs together with state-specific (SS-FNO) and state-averaged (SA-FNO) variants constructed from perturbative one-particle density matrices. We implement these approximations for periodic Moller-Plesset perturbation theory, algebraic diagrammatic construction, and equation-of-motion coupled cluster methods in the PySCF software package and benchmark them for correlation energies, equations of state, fundamental band gaps, and quasiparticle band structures of semiconductors and insulators. FNO truncation reproduces canonical results with small errors while substantially reducing computational cost. SS-FNO enables accurate large-basis band-gap calculations, including quadruple-zeta results, whereas SA-FNO efficiently compresses the virtual space for multi-state band-structure calculations. Combined with basis-set and thermodynamic-limit extrapolations, the framework provides a practical route to high-accuracy correlated calculations in periodic systems.

arXiv:2609.13098•Sep 11, 2026•0 citations
While the periodic equation-of-motion coupled-cluster (EOM-CC) method promises systematic improvement of electronic band gap calculations in solids, its practical application at the singles and doubles level (EOM-CCSD) is hindered by severe finite-size errors in feasible simulation cells. We present a hybrid approach combining EOM-CCSD with the computationally less demanding G W approximation to estimate thermodynamic limit band gaps for several insulators and semiconductors. Our method substantially reduces required cell sizes while maintaining accuracy. Comparisons with experimental gaps and self-consistent G W calculations reveal that deviations in EOM-CCSD predictions correlate with reduced single excitation character of the excited many-electron states. Our work not only provides a computationally tractable approach to EOM-CC calculations in solids but also reveals fundamental insights into the role of single excitations in electronic-structure theory.
Periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) has shown promise for quantitative calculations of band structures in solids. However, its steep computational scaling has limited calculations to relatively coarse k-point meshes, leading to sizable finite-size errors and discrepant estimates of thermodynamic-limit band gaps in recent benchmarks. In this work, we revisit EOM-CCSD band gaps for ten semiconductors and insulators using interacting-bath dynamical embedding theory (ibDET), a systematically improvable Green's function embedding framework that enables dense Brillouin-zone sampling at modest computational cost. By pushing the k-point sampling up to 10 × 10 × 10, well beyond the system sizes accessible in canonical periodic EOM-CCSD calculations, we significantly reduce finite-size errors and obtain stable thermodynamic-limit extrapolations. We further compare G0W0@PBE, G0W0@HF, and EOM-CCSD on an equal footing using the same numerical settings in PySCF. We find that EOM-CCSD yields a mean absolute error of 0.32 eV relative to experimental band gaps for a test set of ten semiconductors and insulators, which is lower than that of G0W0@PBE. For ZnO, EOM-CCSD also accurately describes the Zn 3d-band binding energy, despite overestimating the band gap. These results demonstrate that ibDET offers a practical route to high-accuracy many-body electronic structure calculations in periodic systems.
Excited electronic states are difficult to accurately and tractably model, particularly those exhibiting multiexcitonic character. Scalable single-reference excited state methods typically fail to describe such states as they suffer from the defects of the mean-field approximation and/or neglect higher-rank excitation operators. To address this, we derive the [T] perturbative correction to the quantum self-consistent equation-of-motion unitary coupled cluster singles and doubles (q-sc-EOM-UCCSD) method. The method captures triple excitation effects through a classical step of second-order in many-body perturbation theory following a q-sc-EOM-UCCSD calculation on a quantum computer. The perturbative correction developed has a O(N7)\mathcal{O}(N^7)O(N7) classical computational time scaling. We incorporate [T] into a hybrid compute strategy wherein eigenvectors of the effective Hamiltonian from the hybrid quantum-classical q-sc-EOM-UCCSD algorithm are relayed to a classical computer which performs the postprocessing required for the [T] correction. The benefits of [T] are quantified in a benchmark covering a variety of electronically excited states of the isoelectronic CH+^++ and BH molecules. We find that the addition of the [T] correction to both trotterized and full operator variants of q-sc-EOM-CCSD can offer dramatic improvements over baseline q-sc-EOM-UCCSD. Our assessment clearly demonstrates the importance of considering higher-rank excitation operators in q-sc-EOM-UCCSD, quantifies the overall success of the [T] correction, and discusses some of its limitations.
We present an efficient implementation of the coupled-cluster Green’s function (CCGF) method for simulating photoemission spectra of periodic systems. We formulate the periodic CCGF approach with Brillouin zone sampling in the Gaussian basis at the coupled-cluster singles and doubles (CCSD) level. To enable CCGF calculations of realistic solids, we propose an active-space self-energy correction scheme by combining CCGF with the cheaper many-body perturbation theory ( GW ) and implement the model order reduction (MOR) frequency interpolation technique. We find that the active-space self-energy correction and MOR techniques significantly reduce the computational cost of CCGF while maintaining the high accuracy. We apply the developed CCGF approaches to compute spectral properties and band structure of silicon (Si) and zinc oxide (ZnO) crystals using triple-ζ Gaussian basis sets and medium-size k -point sampling and find good agreement with experimental measurements.
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.
One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. These observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.
We present a unified frozen natural orbital (FNO) framework for periodic Gaussian-orbital post-Hartree-Fock calculations with explicit k-point sampling. The approach includes conventional ground-state FNOs together with state-specific (SS-FNO) and state-averaged (SA-FNO) variants constructed from perturbative one-particle density matrices. We implement these approximations for periodic Moller-Plesset perturbation theory, algebraic diagrammatic construction, and equation-of-motion coupled cluster methods in the PySCF software package and benchmark them for correlation energies, equations of state, fundamental band gaps, and quasiparticle band structures of semiconductors and insulators. FNO truncation reproduces canonical results with small errors while substantially reducing computational cost. SS-FNO enables accurate large-basis band-gap calculations, including quadruple-zeta results, whereas SA-FNO efficiently compresses the virtual space for multi-state band-structure calculations. Combined with basis-set and thermodynamic-limit extrapolations, the framework provides a practical route to high-accuracy correlated calculations in periodic systems.