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:2608.18267

Finite-temperature Green's function cluster expansion from thermofield doubles: Breakdown of the polaron picture

M. R. Carbone, S. Fomichev, B. Kloss, Andrew J. Millis, Mona Berciu, David R. Reichman, John Sous

Published Aug 18, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce a method, numerically exact in principle, for computing the momentum- and frequency-resolved single-particle Green's function of a polaron at finite temperature. The method, which we refer to as the finite-temperature Green's function cluster expansion, combines two ingredients: the generalized Green's function cluster expansion, a numerically exact extension of the momentum average family of methods that solves the polaron problem at zero temperature through a hierarchy of equations of motion for restricted phonon cloud configurations; and the thermofield double formalism, which maps the thermal trace onto a pure-state expectation value over a doubled Hilbert space. The resulting equations of motion have the same algebraic structure as those of the multi-boson zero-temperature theory, with the temperature entering through a Bogoliubov-type mixing angle that controls the coupling to a set of fictitious bath bosons. We implement the method in our open-source software package and benchmark it on the one-dimensional Holstein polaron, computing spectral functions, dispersions, lifetimes, and effective masses across coupling regimes and temperatures up to $T/Ω\sim 1$. Where finite-temperature density matrix renormalization group results are available, we find quantitative agreement at affordable computational cost. The method recovers momentum-resolved spectra directly in frequency space, with no time evolution or analytic continuation. We also discuss the practical costs of the approach. In particular, since the doubled phonon Hilbert space has a non-trivial configuration structure in which real and fictitious clouds compete, convergence in the corresponding cloud parameters requires care.

DOI
10.48550/arxiv.2608.18267

Related papers

Thermofield Theory for Finite-Temperature Electronic Structure

Gaurav Harsha, Thomas M. Henderson, Gustavo E. Scuseria

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. They also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

arXiv:2303.12318•The Journal of Physical Chemistry A•Apr 5, 2023•10 citations

Linked-cluster expansion for the Green's function of the infinite- U Hubbard model

Ehsan Khatami, Edward Perepelitsky, Marcos Rigol +1 more

We implement a highly efficient strong-coupling expansion for the Green's function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green's function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Padé approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green's function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order.

arXiv:1310.8029•Physical Review E•Jun 2, 2014•10 citations

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

Matthew R. Carbone, David R. Reichman, John Sous

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein$+$Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. The GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

arXiv:2103.01972•Physical review. B./Physical review. B•Jul 6, 2021•18 citations

Wave function theories for finite-temperature electronic structure

Gaurav Harsha

Wave function methods have offered a robust, systematically improvable means to study the ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogs of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space. Approximating every state in the Hilbert space is an impossible task. Besides, excited-state theories are not as well studied as ground-state ones. In this thesis, we overcome these difficulties by employing thermofield dynamics, a theory that allows the purification of the ensemble density matrix and constructs a pure wave function that encodes the equilibrium thermal behavior of the system. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we develop a framework to extend ground-state wave function theories to non-zero temperatures. We discuss explicit formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of fermions in the grand canonical ensemble. For fermions in the canonical ensemble, we present mean-field, configuration interaction, and perturbation theories. We also discuss thermal mean-field, configuration interaction, and coupled cluster formalisms for \SU systems. To assess the quality of these approximations, we study both model and small atomic and molecular as benchmark electronic problems while comparing against exact results. For spins, we use the Lipkin-Meshkov-Glick and the one-dimensional transverse field Ising models as our reference systems. Indeed, the thermal methods perform reasonably well while merely adding a prefactor to the computational cost. They also inherit all the properties, good or bad, from their ground-state counterparts, signifying the robustness of our formalism.

•Nov 30, 2021•0 citations

Delocalized Coupled-Cluster Theory for Polaron Structure and Dynamics

Hamlin Wu, Moritz K. A. Baumgarten, Tong Jiang +1 more

Polaron ground states and finite-temperature dynamics remain challenging to simulate because existing methods struggle to combine nonperturbative accuracy, systematic improvability, and scalability from models to materials-specific Hamiltonians. We introduce a translationally invariant variational coupled-cluster (CC) theory for polarons, termed delocalized CC (dCC), with closed-form energies at cost as low as $\mathcal{O}(N^3)$ and no phonon-number cutoff. dCC accurately describes the ground states of the one- and two-dimensional Holstein and Su--Schrieffer--Heeger (optical and bond) models and the Fr{ö}hlich model, in close agreement with density matrix renormalization group (DMRG) and diagrammatic Monte Carlo benchmarks. A projected tangent-space response formalism built on the same ansatz yields electron-addition spectral functions and optical conductivities at zero and finite temperature. The resulting spectra agree well with DMRG, Lanczos, and neural-network quantum-state benchmarks while retaining a physically interpretable excitation hierarchy, and extend to two-dimensional lattices at finite temperature beyond the practical reach of these methods. The same framework applies directly to \textit{ab initio} electron--phonon matrix elements, yielding LiF electron- and hole-polaron binding energies that match state-of-the-art many-body calculations. These results establish dCC as a unified variational framework for polaron ground states and dynamics, from model systems to real materials.

arXiv:2608.04979•Aug 5, 2026•0 citations

Universal thermal breakdown of polaron coherence in one, two, and three dimensions

Jeet Shannigrahi, J. Bonča, Mona Berciu

How a polaron loses its quasiparticle coherence with increasing temperature is a long-standing open problem. Holstein addressed it in 1959 only in the extreme antiadiabatic, strong-coupling limit, while more recent numerically exact approaches are largely restricted to one dimension. Here we solve this problem on square and simple cubic lattices across the weak-, intermediate-, and strong-coupling regimes. We show that the polaron effective mass $m^*$ and inverse lifetime $1/τ$ increase monotonically with temperature until, at $T \sim 0.5\,Ω$, the quasiparticle peak dissolves into a broad incoherent thermal continuum. The phonon frequency $Ω$ therefore defines a universal coherence scale, independent of dimensionality and coupling strength, validating Holstein's prediction far beyond the regime in which it was derived. These results follow from a finite-temperature generalization of the Momentum Average (MA) approximation, yielding a closed-form, diagrammatically derived self-energy that is asymptotically exact in the strong-coupling limit at all temperatures. Benchmark comparisons demonstrate excellent quantitative agreement of the resulting 1D spectral functions with the numerically exact Variational Exact Diagonalization-Finite-Temperature Lanczos Method (VED-FTLM) and finite-$T$ Density Matrix Renormalization Group (DMRG).

arXiv:2608.18529•Aug 19, 2026•0 citations
We introduce a method, numerically exact in principle, for computing the momentum- and frequency-resolved single-particle Green's function of a polaron at finite temperature. The method, which we refer to as the finite-temperature Green's function cluster expansion, combines two ingredients: the generalized Green's function cluster expansion, a numerically exact extension of the momentum average family of methods that solves the polaron problem at zero temperature through a hierarchy of equations of motion for restricted phonon cloud configurations; and the thermofield double formalism, which maps the thermal trace onto a pure-state expectation value over a doubled Hilbert space. The resulting equations of motion have the same algebraic structure as those of the multi-boson zero-temperature theory, with the temperature entering through a Bogoliubov-type mixing angle that controls the coupling to a set of fictitious bath bosons. We implement the method in our open-source software package and benchmark it on the one-dimensional Holstein polaron, computing spectral functions, dispersions, lifetimes, and effective masses across coupling regimes and temperatures up to T/Ω∼1T/Ω\sim 1T/Ω∼1. Where finite-temperature density matrix renormalization group results are available, we find quantitative agreement at affordable computational cost. The method recovers momentum-resolved spectra directly in frequency space, with no time evolution or analytic continuation. We also discuss the practical costs of the approach. In particular, since the doubled phonon Hilbert space has a non-trivial configuration structure in which real and fictitious clouds compete, convergence in the corresponding cloud parameters requires care.
Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. They also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.
We implement a highly efficient strong-coupling expansion for the Green's function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green's function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Padé approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green's function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order.
We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein+++Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. The GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.
Wave function methods have offered a robust, systematically improvable means to study the ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogs of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space. Approximating every state in the Hilbert space is an impossible task. Besides, excited-state theories are not as well studied as ground-state ones. In this thesis, we overcome these difficulties by employing thermofield dynamics, a theory that allows the purification of the ensemble density matrix and constructs a pure wave function that encodes the equilibrium thermal behavior of the system. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we develop a framework to extend ground-state wave function theories to non-zero temperatures. We discuss explicit formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of fermions in the grand canonical ensemble. For fermions in the canonical ensemble, we present mean-field, configuration interaction, and perturbation theories. We also discuss thermal mean-field, configuration interaction, and coupled cluster formalisms for \SU systems. To assess the quality of these approximations, we study both model and small atomic and molecular as benchmark electronic problems while comparing against exact results. For spins, we use the Lipkin-Meshkov-Glick and the one-dimensional transverse field Ising models as our reference systems. Indeed, the thermal methods perform reasonably well while merely adding a prefactor to the computational cost. They also inherit all the properties, good or bad, from their ground-state counterparts, signifying the robustness of our formalism.
Polaron ground states and finite-temperature dynamics remain challenging to simulate because existing methods struggle to combine nonperturbative accuracy, systematic improvability, and scalability from models to materials-specific Hamiltonians. We introduce a translationally invariant variational coupled-cluster (CC) theory for polarons, termed delocalized CC (dCC), with closed-form energies at cost as low as O(N3)\mathcal{O}(N^3)O(N3) and no phonon-number cutoff. dCC accurately describes the ground states of the one- and two-dimensional Holstein and Su--Schrieffer--Heeger (optical and bond) models and the Fr{ö}hlich model, in close agreement with density matrix renormalization group (DMRG) and diagrammatic Monte Carlo benchmarks. A projected tangent-space response formalism built on the same ansatz yields electron-addition spectral functions and optical conductivities at zero and finite temperature. The resulting spectra agree well with DMRG, Lanczos, and neural-network quantum-state benchmarks while retaining a physically interpretable excitation hierarchy, and extend to two-dimensional lattices at finite temperature beyond the practical reach of these methods. The same framework applies directly to \textit{ab initio} electron--phonon matrix elements, yielding LiF electron- and hole-polaron binding energies that match state-of-the-art many-body calculations. These results establish dCC as a unified variational framework for polaron ground states and dynamics, from model systems to real materials.
How a polaron loses its quasiparticle coherence with increasing temperature is a long-standing open problem. Holstein addressed it in 1959 only in the extreme antiadiabatic, strong-coupling limit, while more recent numerically exact approaches are largely restricted to one dimension. Here we solve this problem on square and simple cubic lattices across the weak-, intermediate-, and strong-coupling regimes. We show that the polaron effective mass m∗m^*m∗ and inverse lifetime 1/τ1/τ1/τ increase monotonically with temperature until, at T∼0.5 ΩT \sim 0.5\,ΩT∼0.5Ω, the quasiparticle peak dissolves into a broad incoherent thermal continuum. The phonon frequency ΩΩΩ therefore defines a universal coherence scale, independent of dimensionality and coupling strength, validating Holstein's prediction far beyond the regime in which it was derived. These results follow from a finite-temperature generalization of the Momentum Average (MA) approximation, yielding a closed-form, diagrammatically derived self-energy that is asymptotically exact in the strong-coupling limit at all temperatures. Benchmark comparisons demonstrate excellent quantitative agreement of the resulting 1D spectral functions with the numerically exact Variational Exact Diagonalization-Finite-Temperature Lanczos Method (VED-FTLM) and finite-TTT Density Matrix Renormalization Group (DMRG).