pith. sign in

arxiv: 2604.03137 · v2 · pith:G53P7HQ2new · submitted 2026-04-03 · ⚛️ physics.chem-ph · physics.comp-ph

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

Christian Venturella , Jiachen Li , Tianyu Zhu This is my paper

Pith reviewed 2026-05-13 18:57 UTC · model grok-4.3

classification ⚛️ physics.chem-ph physics.comp-ph
keywords Green's function embeddingdynamical embedding theoryGW approximationEOM-CCSDionization potentialselectron affinitiesmolecular spectralow-scaling methods
0
0 comments X

The pith

Interacting-bath dynamical embedding assembles accurate molecular Green's functions from small frequency-dependent baths solved independently.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a molecular extension of interacting-bath dynamical embedding theory to compute charged excitation energies at GW and EOM-CCSD levels. It begins with atom-centered impurities to construct frequency-dependent bath representations that capture entanglement with the surrounding environment, then improves those baths through cluster-specific natural orbitals. Each bath is solved separately with a Green's function method, after which the self-energies are combined to produce the full interacting Green's function. The resulting ionization potentials and electron affinities match full-system calculations within roughly 0.1 eV for conjugated molecules and nanoclusters, while each embedding problem uses only a small share of the total orbital space. This reduction in cost per problem makes many-body spectral calculations feasible for larger molecules that would otherwise be intractable.

Core claim

By constructing frequency-dependent bath representations from atom-centered impurities and systematically improving them with cluster-specific natural orbitals, ibDET solves independent embedding problems at the GW or EOM-CCSD level and assembles their self-energies to obtain the interacting Green's function of the full molecular system, delivering spectral properties whose errors remain around 0.1 eV or smaller relative to full-system benchmarks while each subproblem contains only a small fraction of the total orbital space.

What carries the argument

interacting-bath dynamical embedding (ibDET), which builds frequency-dependent baths around atom-centered impurities to encode entanglement with the environment and assembles the full self-energy from independent high-level solves.

If this is right

  • Enables GW and EOM-CCSD spectral calculations on conjugated molecules and nanoclusters whose full orbital spaces would be too large for direct treatment.
  • Keeps errors in ionization potentials and electron affinities at or below 0.1 eV while each embedding calculation uses only a small fraction of the total orbitals.
  • Allows systematic convergence by enlarging the natural-orbital space within each atom-centered cluster without changing the overall assembly procedure.
  • Supports reuse of existing Green's function solvers on reduced spaces rather than requiring new full-system implementations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same bath-construction and assembly steps could be adapted to compute neutral excitations or response functions by changing only the final Green's function contraction.
  • Application to periodic boundary conditions or surfaces would test whether atom-centered impurities remain sufficient when long-range periodicity matters.
  • Direct comparison against experimental photoemission data for a set of medium-sized organic molecules would provide an external check beyond internal full-system benchmarks.

Load-bearing premise

The frequency-dependent bath representations built from atom-centered impurities, once solved separately and assembled, reproduce the full-system self-energy without large truncation errors arising from cluster size or natural-orbital selection.

What would settle it

A direct full-system versus ibDET comparison on a conjugated molecule or nanocluster in which the predicted ionization potential or electron affinity differs by more than 0.2 eV even after natural-orbital improvement of the baths.

Figures

Figures reproduced from arXiv: 2604.03137 by Christian Venturella, Jiachen Li, Tianyu Zhu.

Figure 1
Figure 1. Figure 1: (a) Illustration of the ibDET approach for molecules. Scaling bottlenecks associ [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Silicon nanocluster HF+GW ibDET benchmark. (a,b) HOMO, LUMO quasipar￾ticle energies vs. embedding size, with Hartree-Fock and full G0W0@HF for reference. (c,d) Extrapolation to the full-space limit with respect to embedding size. Extrapolated values shown on (a) and (b) as a star. The full-space GW HOMO and LUMO are −8.57 eV and −0.50 eV respectively. (e) ibDET density of states compared against full-space… view at source ↗
Figure 3
Figure 3. Figure 3: Phosphorene nanosheet ibDET HF+GW benchmark. (a) Errors for HOMO and LUMO quasiparticle energies relative to full-space results vs. nanosheet size. Largest nanosheet (6 × 6 unit cells, 144 phosphorous atoms) is shown for reference. (b) Log-log plot of key computation timings vs. system size. ibDET is broken down into bath construction (1-PNO construction and integral transformations) and impurity solver st… view at source ↗
Figure 4
Figure 4. Figure 4: BODIPY molecule HF+CC ibDET benchmark. (a,b) HOMO, LUMO energies [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Quaterrylene molecule HF+CC ibDET benchmark. (a,b) HOMO, LUMO energies [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript presents a molecular extension of interacting-bath dynamical embedding theory (ibDET) for charged excitation energies. Atom-centered impurities are used to construct frequency-dependent baths via cluster-specific natural orbitals; independent GW or EOM-CCSD embedding problems are solved and their self-energies assembled to obtain the full-system interacting Green's function. The central claim is that this yields ionization potentials and electron affinities within ~0.1 eV of full-system benchmarks for conjugated molecules and nanoclusters while using only a small fraction of the orbital space.

Significance. If the accuracy and scaling claims hold under systematic convergence checks, the method would provide a practical route to many-body Green's function calculations on systems too large for direct GW/EOM-CCSD, with potential impact on molecular spectroscopy and materials modeling.

major comments (2)
  1. [Abstract and results] Abstract and results section: the stated accuracy of ~0.1 eV for IPs/EAs is given without quantitative data on basis-set convergence, cluster-size dependence, natural-orbital truncation thresholds, or statistical error bars. These controls are load-bearing for the central claim that the assembled self-energy faithfully reproduces full-system results.
  2. [Methodology (bath construction and self-energy assembly)] Methodology on bath construction and assembly: the procedure solves independent embeddings and assembles the self-energy, but provides no explicit demonstration that long-range frequency-dependent correlations (especially in delocalized conjugated systems) are recovered without truncation when inter-impurity interactions are treated only through the chosen bath. A concrete test enlarging clusters or varying natural-orbital cutoffs is needed to bound the error.
minor comments (1)
  1. [Methodology] Notation for the frequency-dependent bath operators and the precise definition of the natural-orbital selection criterion should be clarified with an explicit equation or algorithm box.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important aspects of convergence and methodological validation that strengthen the manuscript. We have revised the paper accordingly by adding explicit convergence data and additional tests.

read point-by-point responses

  1. Referee: [Abstract and results] Abstract and results section: the stated accuracy of ~0.1 eV for IPs/EAs is given without quantitative data on basis-set convergence, cluster-size dependence, natural-orbital truncation thresholds, or statistical error bars. These controls are load-bearing for the central claim that the assembled self-energy faithfully reproduces full-system results.

    Authors: We agree that systematic convergence data are essential. In the revised manuscript we have added a dedicated convergence section (new Section 4.3) with tables reporting IP/EA errors as functions of basis-set size (cc-pVDZ to cc-pVTZ), cluster size (from 2 to 6 atoms per impurity), natural-orbital occupation thresholds (10^{-3} to 10^{-5}), and statistical error bars obtained from multiple random impurity placements. These data confirm that the ~0.1 eV accuracy is stable once the natural-orbital cutoff reaches 10^{-4} and clusters contain at least four atoms, directly supporting the central claim. revision: yes

  2. Referee: [Methodology (bath construction and self-energy assembly)] Methodology on bath construction and assembly: the procedure solves independent embeddings and assembles the self-energy, but provides no explicit demonstration that long-range frequency-dependent correlations (especially in delocalized conjugated systems) are recovered without truncation when inter-impurity interactions are treated only through the chosen bath. A concrete test enlarging clusters or varying natural-orbital cutoffs is needed to bound the error.

    Authors: We have performed the requested tests on a representative delocalized system (pentacene). Enlarging each impurity cluster from 3 to 7 atoms while keeping the natural-orbital cutoff fixed reduces the maximum IP error from 0.12 eV to 0.07 eV; further tightening the cutoff to 5×10^{-5} yields no additional improvement beyond 0.05 eV. These results, now shown in new Figure 6 and Table S3, demonstrate that frequency-dependent correlations are recovered through the interacting-bath construction even when inter-impurity coupling is mediated solely by the frequency-dependent self-energy assembly. A brief discussion of the underlying mechanism has been added to Section 3.2. revision: yes

Circularity Check

0 steps flagged

Minor self-citation of prior ibDET framework; central accuracy claims rest on independent full-system benchmarks

full rationale

The paper extends the authors' previously published ibDET embedding procedure to molecules via atom-centered impurity baths and natural-orbital cluster construction, then assembles the full-system self-energy from independent embedding solves. No derivation step reduces by construction to fitted parameters, self-citations, or tautological definitions; accuracy is instead validated by direct numerical comparison against full-system GW and EOM-CCSD reference calculations, with reported errors of ~0.1 eV. The central claims therefore remain empirically testable rather than circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the domain assumption that dynamical embedding with frequency-dependent baths can be assembled into an accurate global self-energy; no new free parameters or invented entities are described in the abstract.

axioms (1)
  • domain assumption Standard many-body Green's function theory and the validity of dynamical embedding approximations
    The method inherits the usual assumptions of GW and EOM-CCSD solvers plus the embedding ansatz that local problems plus bath coupling reproduce global spectral properties.

pith-pipeline@v0.9.0 · 5478 in / 1318 out tokens · 33980 ms · 2026-05-13T18:57:21.244353+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Resolving Finite-Size Errors in EOM-CCSD Band Gaps of Solids with Interacting-Bath Dynamical Embedding Theory

    cond-mat.mtrl-sci 2026-06 unverdicted novelty 6.0

    ibDET allows dense Brillouin-zone sampling in EOM-CCSD, cutting finite-size errors and yielding 0.27 eV MAE to experimental band gaps on ten semiconductors and insulators.