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

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

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

Published Jul 6, 2021• 18 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.104.035106

Related papers

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

M. R. Carbone, S. Fomichev, B. Kloss +4 more

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.

arXiv:2608.18267•Aug 18, 2026•0 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

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

Niko Säkkinen, Yang Peng, Heiko Appel +1 more

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

arXiv:1507.04726•The Journal of Chemical Physics•Dec 15, 2015•36 citations

Exact single and two particle non-equilibrium Green's functions for the time-dependent Holstein-Hubbard model in the atomic limit

R. D. Nesselrodt, James K. Freericks

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales -- for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections -- but strongly correlated electrons require a more complex description and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We derive nonequilibrium spectral moment sum rules for the first four moments of the single-particle Green's function. We then analytically confirm that this Green's function exactly satisfies the zeroth, first, and second moment sum rules, and provide numerical evidence that the third spectral moment sum rule also holds. Finally we calculate the two-particle Green's function in the atomic limit. These two functions are required inputs to create a self-consistent strong-coupling expansion for the lattice Green's functions of the Hubbard-Holstein model.

•Mar 26, 2021•0 citations

Green’s function of the Holstein polaron

Glen L. Goodvin, Mona Berciu, G. A. Sawatzky

We present a highly efficient yet accurate analytical approximation for the Green's function of a Holstein polaron. It is obtained by summing all the self-energy diagrams, but with each self-energy diagram averaged over the momenta of its free propagators. The result becomes exact for both zero bandwidth and for zero electron-phonon coupling and is accurate everywhere in the parameter space. The resulting Green's function satisfies exactly the first six spectral weight sum rules. All higher sum rules are satisfied with great accuracy, becoming asymptotically exact for coupling both much larger and much smaller than the free particle bandwidth. Comparison with existing numerical data also confirms this accuracy. We use this approximation to analyze in detail the redistribution of the spectral weight as the coupling strength varies.

arXiv:cond-mat/0609597•Physical Review B•Dec 5, 2006•128 citations

Green’s Function of a Dressed Particle

Mona Berciu

We present a new, highly efficient yet accurate approximation for the Green's functions of dressed particles, using the Holstein polaron as an example. Instead of summing a subclass of self-energy diagrams (e.g., the noncrossed ones, in the self-consistent Born approximation), we sum all the diagrams, but with each diagram averaged over its free propagators' momenta. The resulting Green's function satisfies exactly the first six spectral weight sum rules. All higher sum rules are satisfied with great accuracy, becoming asymptotically exact for coupling both much larger and much smaller than the free particle bandwidth. Possible generalizations to other models are also discussed.

arXiv:cond-mat/0602195•Physical Review Letters•Jul 19, 2006•145 citations
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.
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.
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 present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.
One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales -- for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections -- but strongly correlated electrons require a more complex description and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We derive nonequilibrium spectral moment sum rules for the first four moments of the single-particle Green's function. We then analytically confirm that this Green's function exactly satisfies the zeroth, first, and second moment sum rules, and provide numerical evidence that the third spectral moment sum rule also holds. Finally we calculate the two-particle Green's function in the atomic limit. These two functions are required inputs to create a self-consistent strong-coupling expansion for the lattice Green's functions of the Hubbard-Holstein model.
We present a highly efficient yet accurate analytical approximation for the Green's function of a Holstein polaron. It is obtained by summing all the self-energy diagrams, but with each self-energy diagram averaged over the momenta of its free propagators. The result becomes exact for both zero bandwidth and for zero electron-phonon coupling and is accurate everywhere in the parameter space. The resulting Green's function satisfies exactly the first six spectral weight sum rules. All higher sum rules are satisfied with great accuracy, becoming asymptotically exact for coupling both much larger and much smaller than the free particle bandwidth. Comparison with existing numerical data also confirms this accuracy. We use this approximation to analyze in detail the redistribution of the spectral weight as the coupling strength varies.
We present a new, highly efficient yet accurate approximation for the Green's functions of dressed particles, using the Holstein polaron as an example. Instead of summing a subclass of self-energy diagrams (e.g., the noncrossed ones, in the self-consistent Born approximation), we sum all the diagrams, but with each diagram averaged over its free propagators' momenta. The resulting Green's function satisfies exactly the first six spectral weight sum rules. All higher sum rules are satisfied with great accuracy, becoming asymptotically exact for coupling both much larger and much smaller than the free particle bandwidth. Possible generalizations to other models are also discussed.