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:2208.06346 · Physical Review Letters

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization

Jon Lafuente-Bartolome, Chao Lian, Weng Hong Sio, Idoia G. Gurtubay, Asier Eiguren, Feliciano Giustino

Published Aug 9, 2022• 54 citations
arXivHTMLPDF

Abstract

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green's function theory of band structure renormalization that incorporates localization and self-trapping. We show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

DOI
10.1103/physrevlett.129.076402

Related papers

Renormalization group approach to the Fröhlich polaron model: application to impurity-BEC problem

Fabian Grusdt, Yulia E. Shchadilova, Alexey N. Rubtsov +1 more

When a mobile impurity interacts with a many-body system, such as a phonon bath, a polaron is formed. Despite the importance of the polaron problem for a wide range of physical systems, a unified theoretical description valid for arbitrary coupling strengths is still lacking. Here we develop a renormalization group approach for analyzing a paradigmatic model of polarons, the so-called Fröhlich model, and apply it to a problem of impurity atoms immersed in a Bose-Einstein condensate of ultra cold atoms. Polaron energies obtained by our method are in excellent agreement with recent diagrammatic Monte Carlo calculations for a wide range of interaction strengths. They are found to be logarithmically divergent with the ultra-violet cut-off, but physically meaningful regularized polaron energies are also presented. Moreover, we calculate the effective mass of polarons and find a smooth crossover from weak to strong coupling regimes. Possible experimental tests of our results in current experiments with ultra cold atoms are discussed.

arXiv:1410.2203•Scientific Reports•Jul 17, 2015•111 citations

Ab initio self-consistent many-body theory of polarons at all couplings

Jon Lafuente-Bartolome, Chao Lian, Weng Hong Sio +3 more

A self-consistent many-body theory of polarons is presented. This unifies weak- and strong-coupling approaches to the electron-phonon interaction from first principles, opening the door toward a full characterization of polaronic effects in materials at all coupling strengths.

arXiv:2208.06352•Physical review. B./Physical review. B•Aug 9, 2022•45 citations

Many-body Green's function approaches to the doped Fröhlich solid: Exact solutions and anomalous mass enhancement

Nikolaus Kandolf, Carla Verdi, Feliciano Giustino

In polar semiconductors and insulators, the Fr\"ohlich interaction between electrons and long-wavelength longitudinal optical phonons induces a many-body renormalization of the carrier effective masses and the appearance of characteristic phonon sidebands in the spectral function, commonly dubbed ``polaron satellites.'' The simplest model that captures these effects is the Fr\"ohlich model, whereby electrons in a parabolic band interact with a dispersionless longitudinal optical phonon. The Fr\"ohlich model has been employed in a number of seminal papers, from early perturbation-theory approaches to modern diagrammatic Monte Carlo calculations. One limitation of this model is that it focuses on undoped systems, thus ignoring carrier screening and Pauli blocking effects that are present in real experiments on doped samples. To overcome this limitation, we here extend the Fr\"ohlich model to the case of doped systems, and we provide exact solutions for the electron spectral function, mass enhancement, and polaron satellites. We perform the analysis using two approaches, namely, Dyson's equation with the Fan-Migdal self-energy, and the second-order cumulant expansion. We find that these two approaches provide qualitatively different results. In particular, Dyson's approach yields better quasiparticle masses and worse satellites, while the cumulant approach provides better satellite structures, at the price of worse quasiparticle masses. Both approaches yield an anomalous enhancement of the electron effective mass at finite doping levels, which in turn leads to a breakdown of the quasiparticle picture in a significant portion of the phase diagram.

arXiv:2202.05365•Physical review. B./Physical review. B•Feb 25, 2022•11 citations

Polarons from first principles

Zhenbang Dai, Jon Lafuente-Bartolome, Feliciano Giustino

This article reviews recent theoretical developments in the ab initio study of polarons in materials. The polaron is an emergent quasiparticle that arises from the interaction between electrons and phonons in solids, and consists of an electron or a hole accompanied by a distortion of the crystal lattice. Recent advances in experiments, theory, and computation have made it possible to investigate these quasiparticles with unprecedented detail, reigniting the interest in this classic problem of condensed matter physics. Recent theoretical and computational advances include ab initio calculations of polaron spectral functions, wavefunctions, lattice distortions, and transport and optical properties. These developments provide new insight into polaron physics, but they have evolved somewhat independently from the earlier effective Hamiltonian approaches that laid the foundation of the field. This article aims to bridge these complementary perspectives by placing them within a single unified conceptual framework. To this end, we start by reviewing effective Hamiltonians of historical significance in polaron theory, ab initio techniques based on density functional theory, and many-body first-principles approaches to polarons. After this survey, we outline a general field-theoretic framework that bridges between these diverse approaches to polaron physics. For completeness, we also review recent progress in the study of exciton polarons and self-trapped excitons and their relations to polarons. Beyond the methodology, we discuss recent applications to several classes of materials that attracted attention in the context of polaron physics.

cond-mat.mtrl-sciCross-listedarXiv:2512.06176•Dec 9, 2025

The effects of electron-phonon interactions on bandgaps

J. P. Hague

I compute the spectral functions of a 1D Holstein polaron moving in a modulated potential, to examine the effects of electron-phonon interactions on band gaps. The imaginary time Green function is computed using diagrammatic quantum Monte Carlo, which exactly sums the diagrammatic series. From the imaginary time Green function, spectral functions are computed. The electron-phonon interaction flattens the electronic dispersion and leads to an increase in the gap at momentum k = π/2. At strong coupling, polaron sidebands form in the gap. These results demonstrate the strong effect that electron-phonon interactions can have on band gaps.

Journal of Physics Conference Series•Mar 1, 2011•3 citations

Ab initioelectron propagators in molecules with strong electron-phonon interaction: II. Electron Green’s function

Yuri Dahnovsky

Ab initio electron propagator methods are developed to study electronic properties of molecular systems with strong electron-electron and electron-phonon interactions. For the calculation of electron Green's functions we apply a canonical small polaron transformation that intrinsically contains strong electron-phonon effects. In the transformed Hamiltonian, the energy levels for the noninteracting particles are shifted down by the relaxation (solvation) energies. The Coulomb integrals are also renormalized by the electron-phonon interaction. For certain values of the electron-phonon coupling constants, the renormalized Coulomb integrals can be negative which implies the attraction between two electrons. Within the small polaron transformation we develop a diagrammatic technique for the calculation of electron Green's function in which the electron-phonon interaction is already included into the multiple phonon correlation functions. Since the decoupling of the phonon correlation functions is impossible, and therefore, a Wick's theorem for such correlation functions is invalid, there is no Dyson equation for the electron Green's function. To find the electron Green's function, we use different approximations. One of them is a link-cluster approximation that includes diagonal transitions for the renormalized zeroth Green's function. In the linked-cluster approach the Dyson equation is derived in the most general case, where the self-energy operator is an arbitrary functional (not only in the Hartree-Fock approximation). It is shown that even a Hartree-Fock electron (hole) is not a particle any longer. It is a quasiparticle with a finite lifetime that depends on energy of particle and hole states in different ways. As a consequence of this, a standard description of a Hartree-Fock approximation in terms of wave functions becomes inappropriate in this problem. To challenge the linked-cluster approximation we develop a different approach: a sequential propagation approximation where scattering events occur only for sequential transitions. A self-consistent Hartree-Fock equation for a four-index Green's function matrix is derived. In conclusion, the proposed schemes can be considered for future method developments for quantum chemical calculations for large molecules with strong nonadiabatic effects, e-e correlated electron transfer reactions, and electron transport in molecular transport junctions.

The Journal of Chemical Physics•Jul 2, 2007•16 citations
Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green's function theory of band structure renormalization that incorporates localization and self-trapping. We show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.
When a mobile impurity interacts with a many-body system, such as a phonon bath, a polaron is formed. Despite the importance of the polaron problem for a wide range of physical systems, a unified theoretical description valid for arbitrary coupling strengths is still lacking. Here we develop a renormalization group approach for analyzing a paradigmatic model of polarons, the so-called Fröhlich model, and apply it to a problem of impurity atoms immersed in a Bose-Einstein condensate of ultra cold atoms. Polaron energies obtained by our method are in excellent agreement with recent diagrammatic Monte Carlo calculations for a wide range of interaction strengths. They are found to be logarithmically divergent with the ultra-violet cut-off, but physically meaningful regularized polaron energies are also presented. Moreover, we calculate the effective mass of polarons and find a smooth crossover from weak to strong coupling regimes. Possible experimental tests of our results in current experiments with ultra cold atoms are discussed.
In polar semiconductors and insulators, the Fr\"ohlich interaction between electrons and long-wavelength longitudinal optical phonons induces a many-body renormalization of the carrier effective masses and the appearance of characteristic phonon sidebands in the spectral function, commonly dubbed ``polaron satellites.'' The simplest model that captures these effects is the Fr\"ohlich model, whereby electrons in a parabolic band interact with a dispersionless longitudinal optical phonon. The Fr\"ohlich model has been employed in a number of seminal papers, from early perturbation-theory approaches to modern diagrammatic Monte Carlo calculations. One limitation of this model is that it focuses on undoped systems, thus ignoring carrier screening and Pauli blocking effects that are present in real experiments on doped samples. To overcome this limitation, we here extend the Fr\"ohlich model to the case of doped systems, and we provide exact solutions for the electron spectral function, mass enhancement, and polaron satellites. We perform the analysis using two approaches, namely, Dyson's equation with the Fan-Migdal self-energy, and the second-order cumulant expansion. We find that these two approaches provide qualitatively different results. In particular, Dyson's approach yields better quasiparticle masses and worse satellites, while the cumulant approach provides better satellite structures, at the price of worse quasiparticle masses. Both approaches yield an anomalous enhancement of the electron effective mass at finite doping levels, which in turn leads to a breakdown of the quasiparticle picture in a significant portion of the phase diagram.
This article reviews recent theoretical developments in the ab initio study of polarons in materials. The polaron is an emergent quasiparticle that arises from the interaction between electrons and phonons in solids, and consists of an electron or a hole accompanied by a distortion of the crystal lattice. Recent advances in experiments, theory, and computation have made it possible to investigate these quasiparticles with unprecedented detail, reigniting the interest in this classic problem of condensed matter physics. Recent theoretical and computational advances include ab initio calculations of polaron spectral functions, wavefunctions, lattice distortions, and transport and optical properties. These developments provide new insight into polaron physics, but they have evolved somewhat independently from the earlier effective Hamiltonian approaches that laid the foundation of the field. This article aims to bridge these complementary perspectives by placing them within a single unified conceptual framework. To this end, we start by reviewing effective Hamiltonians of historical significance in polaron theory, ab initio techniques based on density functional theory, and many-body first-principles approaches to polarons. After this survey, we outline a general field-theoretic framework that bridges between these diverse approaches to polaron physics. For completeness, we also review recent progress in the study of exciton polarons and self-trapped excitons and their relations to polarons. Beyond the methodology, we discuss recent applications to several classes of materials that attracted attention in the context of polaron physics.
I compute the spectral functions of a 1D Holstein polaron moving in a modulated potential, to examine the effects of electron-phonon interactions on band gaps. The imaginary time Green function is computed using diagrammatic quantum Monte Carlo, which exactly sums the diagrammatic series. From the imaginary time Green function, spectral functions are computed. The electron-phonon interaction flattens the electronic dispersion and leads to an increase in the gap at momentum k = π/2. At strong coupling, polaron sidebands form in the gap. These results demonstrate the strong effect that electron-phonon interactions can have on band gaps.
Ab initio electron propagator methods are developed to study electronic properties of molecular systems with strong electron-electron and electron-phonon interactions. For the calculation of electron Green's functions we apply a canonical small polaron transformation that intrinsically contains strong electron-phonon effects. In the transformed Hamiltonian, the energy levels for the noninteracting particles are shifted down by the relaxation (solvation) energies. The Coulomb integrals are also renormalized by the electron-phonon interaction. For certain values of the electron-phonon coupling constants, the renormalized Coulomb integrals can be negative which implies the attraction between two electrons. Within the small polaron transformation we develop a diagrammatic technique for the calculation of electron Green's function in which the electron-phonon interaction is already included into the multiple phonon correlation functions. Since the decoupling of the phonon correlation functions is impossible, and therefore, a Wick's theorem for such correlation functions is invalid, there is no Dyson equation for the electron Green's function. To find the electron Green's function, we use different approximations. One of them is a link-cluster approximation that includes diagonal transitions for the renormalized zeroth Green's function. In the linked-cluster approach the Dyson equation is derived in the most general case, where the self-energy operator is an arbitrary functional (not only in the Hartree-Fock approximation). It is shown that even a Hartree-Fock electron (hole) is not a particle any longer. It is a quasiparticle with a finite lifetime that depends on energy of particle and hole states in different ways. As a consequence of this, a standard description of a Hartree-Fock approximation in terms of wave functions becomes inappropriate in this problem. To challenge the linked-cluster approximation we develop a different approach: a sequential propagation approximation where scattering events occur only for sequential transitions. A self-consistent Hartree-Fock equation for a four-index Green's function matrix is derived. In conclusion, the proposed schemes can be considered for future method developments for quantum chemical calculations for large molecules with strong nonadiabatic effects, e-e correlated electron transfer reactions, and electron transport in molecular transport junctions.