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:2604.03137 · Journal of Chemical Theory and Computation

Low-Scaling Many-Body Green’s Function Calculations for Molecular Systems via Interacting-Bath Dynamical Embedding Theory

Christian Venturella, Jiachen Li, Tianyu Zhu

Published Jun 25, 2026• 1 citation
arXivHTMLPDF

Abstract

We present a molecular extension of our recently proposed Green’s function embedding method, interacting-bath dynamical embedding theory (ibDET), for computing charged excitation energies at the GW and EOM-CCSD levels. Starting from atom-centered impurities, we construct bath representations that capture the frequency-dependent entanglement between the impurity and its environment and can be systematically improved via the construction of cluster-specific natural orbitals. Utilizing a GW or coupled-cluster Green’s function solver, the self-energy of the full system is assembled from all embedding problems to obtain the interacting Green’s function. We show that ibDET provides accurate spectral properties with much reduced cost for a broad range of systems, including conjugated molecules and nanoclusters. Compared with full-system results, the errors in the predicted ionization potentials and electron affinities are around 0.1 eV or smaller, while each embedding problem includes only a small fraction of the total orbital space. This work provides an efficient and scalable framework for computing spectral properties of molecular systems.

DOI
10.1021/acs.jctc.6c00663

Related papers

Interacting-Bath Dynamical Embedding for Capturing Nonlocal Electron Correlation in Solids

Jiachen Li, Tianyu Zhu

Quantitative simulation of electronic structure of solids requires treating local and nonlocal electron correlations on an equal footing. We present a new ab initio formulation of Green's function embedding which, unlike dynamical mean-field theory that uses noninteracting bath, derives bath representation with general two-particle interactions in a systematically improvable manner. The resulting interacting-bath dynamical embedding theory (ibDET) utilizes an efficient real-axis coupled-cluster solver to compute the self-energy, approaching the full system limit at much reduced cost. When combined with the GW theory, GW+ibDET achieves good agreement with experimental spectral properties across a range of semiconducting, insulating, and metallic materials. Our approach also enables quantifying the role of nonlocal electron correlation in determining material properties and addressing the long-standing debate on the bandwidth narrowing of metallic sodium.

arXiv:2406.07531•Physical Review Letters•Nov 20, 2024•15 citations

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

Coupled Cluster as an Impurity Solver for Green’s Function Embedding Methods

Avijit Shee, Dominika Zgid

We investigate the performance of Green's function coupled cluster singles and doubles (CCSD) method as a solver for Green's function embedding methods. To develop an efficient CC solver, we construct the one-particle Green's function from the coupled cluster (CC) wave function based on the non-Hermitian Lanczos algorithm. The major advantage of this method is that its scaling does not depend on the number of frequency points. We have tested the applicability of the CC Green's function solver in the weakly to strongly correlated regimes by employing it for a half-filled one-dimensional (1D) Hubbard model projected onto a single-site impurity problem and a half-filled two-dimensional (2D) Hubbard model projected onto a four-site impurity problem. For the 1D Hubbard model, for all interaction strengths, we observe an excellent agreement with the full configuration interaction (FCI) technique, both for the self-energy and spectral function. For the 2D Hubbard, we have employed an open-shell version of the current implementation and observed some discrepancies from FCI in the strongly correlated regime. Finally, on an example of a small ammonia cluster, we analyze the performance of Green's function CCSD solver within the Self-Energy Embedding Theory (SEET) with Hartree-Fock (HF) and Green's Function second order (GF2) for the treatment of the environment.

arXiv:1906.04079•Journal of Chemical Theory and Computation•Sep 13, 2019•51 citations

Triple Excitations in Green’s Function Coupled Cluster Solver for Studies of Strongly Correlated Systems in the Framework of Self-Energy Embedding Theory

Avijit Shee, Chia-Nan Yeh, Bo Peng +2 more

Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handling a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green’s function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green’s function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO 3 solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.

arXiv:2211.12680•The Journal of Physical Chemistry Letters•Mar 1, 2023•5 citations

Spectral functions of strongly correlated extended systems via an exact quantum embedding

George H. Booth, Garnet Kin‐Lic Chan

Density matrix embedding theory (DMET) [Phys. Rev. Lett. 109, 186404 (2012)], introduced an approach to quantum cluster embedding methods whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost but was inherently limited by the construction of a bath designed to reproduce ground-state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency-dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the one- (1D) and two-dimensional (2D) Hubbard model, where we obtain zero temperature and thermodynamic limit spectral functions, and show the trivial extension to two-particle Green's functions. This advance therefore extends the scope and applicability of DMET in condensed-matter problems as a computationally tractable route to correlated spectral functions of extended systems and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.

arXiv:1309.2320•Physical Review B•Apr 7, 2015•60 citations

Machine Learning Many-Body Green’s Functions for Molecular Excitation Spectra

Christian Venturella, Christopher Hillenbrand, Jiachen Li +1 more

We present a machine learning (ML) framework for predicting Green’s functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green’s functions through analytic continuation and Dyson’s equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory ( GW ) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW @PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.

arXiv:2310.09911•Journal of Chemical Theory and Computation•Dec 27, 2023•17 citations
We present a molecular extension of our recently proposed Green’s function embedding method, interacting-bath dynamical embedding theory (ibDET), for computing charged excitation energies at the GW and EOM-CCSD levels. Starting from atom-centered impurities, we construct bath representations that capture the frequency-dependent entanglement between the impurity and its environment and can be systematically improved via the construction of cluster-specific natural orbitals. Utilizing a GW or coupled-cluster Green’s function solver, the self-energy of the full system is assembled from all embedding problems to obtain the interacting Green’s function. We show that ibDET provides accurate spectral properties with much reduced cost for a broad range of systems, including conjugated molecules and nanoclusters. Compared with full-system results, the errors in the predicted ionization potentials and electron affinities are around 0.1 eV or smaller, while each embedding problem includes only a small fraction of the total orbital space. This work provides an efficient and scalable framework for computing spectral properties of molecular systems.
Quantitative simulation of electronic structure of solids requires treating local and nonlocal electron correlations on an equal footing. We present a new ab initio formulation of Green's function embedding which, unlike dynamical mean-field theory that uses noninteracting bath, derives bath representation with general two-particle interactions in a systematically improvable manner. The resulting interacting-bath dynamical embedding theory (ibDET) utilizes an efficient real-axis coupled-cluster solver to compute the self-energy, approaching the full system limit at much reduced cost. When combined with the GW theory, GW+ibDET achieves good agreement with experimental spectral properties across a range of semiconducting, insulating, and metallic materials. Our approach also enables quantifying the role of nonlocal electron correlation in determining material properties and addressing the long-standing debate on the bandwidth narrowing of metallic sodium.
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.
We investigate the performance of Green's function coupled cluster singles and doubles (CCSD) method as a solver for Green's function embedding methods. To develop an efficient CC solver, we construct the one-particle Green's function from the coupled cluster (CC) wave function based on the non-Hermitian Lanczos algorithm. The major advantage of this method is that its scaling does not depend on the number of frequency points. We have tested the applicability of the CC Green's function solver in the weakly to strongly correlated regimes by employing it for a half-filled one-dimensional (1D) Hubbard model projected onto a single-site impurity problem and a half-filled two-dimensional (2D) Hubbard model projected onto a four-site impurity problem. For the 1D Hubbard model, for all interaction strengths, we observe an excellent agreement with the full configuration interaction (FCI) technique, both for the self-energy and spectral function. For the 2D Hubbard, we have employed an open-shell version of the current implementation and observed some discrepancies from FCI in the strongly correlated regime. Finally, on an example of a small ammonia cluster, we analyze the performance of Green's function CCSD solver within the Self-Energy Embedding Theory (SEET) with Hartree-Fock (HF) and Green's Function second order (GF2) for the treatment of the environment.
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handling a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green’s function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green’s function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO 3 solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.
Density matrix embedding theory (DMET) [Phys. Rev. Lett. 109, 186404 (2012)], introduced an approach to quantum cluster embedding methods whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost but was inherently limited by the construction of a bath designed to reproduce ground-state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency-dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the one- (1D) and two-dimensional (2D) Hubbard model, where we obtain zero temperature and thermodynamic limit spectral functions, and show the trivial extension to two-particle Green's functions. This advance therefore extends the scope and applicability of DMET in condensed-matter problems as a computationally tractable route to correlated spectral functions of extended systems and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.
We present a machine learning (ML) framework for predicting Green’s functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green’s functions through analytic continuation and Dyson’s equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory ( GW ) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW @PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.