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

Spectral Functions of the Holstein Polaron: Exact and Approximate Solutions

Petar Mitrić, Veljko Janković, Nenad Vukmirović, D. Tanasković

Published Aug 22, 2022• 30 citations
arXivHTMLPDF

Abstract

It is generally accepted that the dynamical mean field theory gives a good solution of the Holstein model, but only in dimensions greater than two. Here, we show that this theory, which becomes exact in the weak coupling and in the atomic limit, provides an excellent, numerically cheap, approximate solution for the spectral function of the Holstein model in the whole range of parameters, even in one dimension. To establish this, we make a detailed comparison with the spectral functions that we obtain using the newly developed momentum-space numerically exact hierarchical equations of motion method, which yields electronic correlation functions directly in real time. We crosscheck these conclusions with our path integral quantum Monte Carlo and exact diagonalization results, as well as with the available numerically exact results from the literature.

DOI
10.1103/physrevlett.129.096401

Related papers

Cumulant expansion in the Holstein model: Spectral functions and mobility

Petar Mitrić, Veljko Janković, Nenad Vukmirović +1 more

We examine the range of validity of the second-order cumulant expansion (CE) for the calculation of spectral functions, quasiparticle properties, and mobility of the Holstein polaron. We devise an efficient numerical implementation that allows us to make comparisons in a broad interval of temperature, electron-phonon coupling, and phonon frequency. For a benchmark, we use the dynamical mean-field theory which gives, as we have recently shown, rather accurate spectral functions in the whole parameter space, even in low dimensions. We find that in one dimension, the CE resolves well both the quasiparticle and the first satellite peak in a regime of intermediate coupling. At high temperatures, the charge mobility assumes a power law $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}2}$ in the limit of weak coupling and $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}3/2}$ for stronger coupling. We find that, for stronger coupling, the CE gives slightly better results than the self-consistent Migdal approximation (SCMA), while the one-shot Migdal approximation is appropriate only for a very weak electron-phonon interaction. We also analyze the atomic limit and the spectral sum rules. We derive an analytical expression for the moments in CE and find that they are exact up to the fourth order, as opposed to the SCMA where they are exact to the third order. Finally, we analyze the results in higher dimensions.

arXiv:2212.13846•Physical review. B./Physical review. B•Mar 30, 2023•16 citations

Phonon spectral function of the Holstein polaron

Jan Loos, Martin Hohenadler, Andreas Alvermann +1 more

The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations.The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion.Particular emphasis is laid on the different physics occurring in the adiabatic and antiadiabatic regimes, respectively.Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit.This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.

•Jan 1, 2006•12 citations

Phonon spectral function of the Holstein polaron

J Loos, M Hohenadler, A Alvermann +1 more

The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations. The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion. Particular emphasis is laid on the different physics occurring in the adiabatic and anti-adiabatic regimes, respectively. Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit. This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.

arXiv:cond-mat/0604040•Journal of Physics Condensed Matter•Jul 21, 2006•20 citations

Spectral and thermodynamic properties of the Holstein polaron: Hierarchical equations of motion approach

Veljko Janković, Nenad Vukmirović

We develop a hierarchical equations of motion (HEOM) approach to compute real-time single-particle correlation functions and thermodynamic properties of the Holstein model at finite temperature. We exploit the conservation of the total momentum of the system to formulate the momentum-space HEOM whose dynamical variables explicitly keep track of momentum exchanges between the electron and phonons. Our symmetry-adapted HEOM enable us to overcome the numerical instabilities inherent to the commonly used real-space HEOM. The HEOM method is then used to study the spectral function and thermodynamic quantities of chains containing up to ten sites. The HEOM results compare favorably to existing literature. To provide an independent assessment of the HEOM approach and to gain insight into the importance of finite-size effects, we devise a quantum Monte Carlo (QMC) procedure to evaluate finite-temperature single-particle correlation functions in imaginary time and apply it to chains containing up to twenty sites. QMC results reveal that finite-size effects are quite weak, so that the results on 5 to 10-site chains, depending on the parameter regime, are representative of larger systems. A detailed comparison between the HEOM and QMC data place our HEOM method among reliable methods to compute real-time finite-temperature correlation functions in parameter regimes ranging from low- to high-temperature, and weak- to strong-coupling regime.

arXiv:2111.06430•Physical review. B./Physical review. B•Feb 18, 2022•24 citations

Spectral function of electron-phonon models by cluster perturbation theory

Martin Hohenadler, Markus Aichhorn, Wolfgang von der Linden

Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.

arXiv:cond-mat/0306740•Physical review. B, Condensed matter•Nov 19, 2003•56 citations

Identification of the transport regimes of the one-dimensional Holstein model

Suzana Miladić‐Tešić, Nenad Vukmirović

The Holstein model is a benchmark model of systems with electron-phonon interaction. However, its electrical transport properties are not yet fully understood. In this work, we performed numerically exact calculations of the imaginary-time current-current correlation function of the Holstein Hamiltonian for a broad range of model parameters. These calculations were performed using a path-integral-based quantum Monte Carlo method. We compared these results with the results obtained under the assumption of conventional band transport, small polaron hopping, and polaron band transport. From this comparison, we identified the regions in parameter space where each of these transport regimes is valid. In some cases when an imaginary-time comparison of current-current correlation functions could not give conclusive results, we complemented it with real-time comparisons or a comparison with the literature data on numerically exact dc mobilities. Overall, we found that the parameter space is almost completely covered by the three mentioned transport regimes.

arXiv:2509.01335•Physical review. B./Physical review. B•Aug 12, 2025•3 citations
It is generally accepted that the dynamical mean field theory gives a good solution of the Holstein model, but only in dimensions greater than two. Here, we show that this theory, which becomes exact in the weak coupling and in the atomic limit, provides an excellent, numerically cheap, approximate solution for the spectral function of the Holstein model in the whole range of parameters, even in one dimension. To establish this, we make a detailed comparison with the spectral functions that we obtain using the newly developed momentum-space numerically exact hierarchical equations of motion method, which yields electronic correlation functions directly in real time. We crosscheck these conclusions with our path integral quantum Monte Carlo and exact diagonalization results, as well as with the available numerically exact results from the literature.
We examine the range of validity of the second-order cumulant expansion (CE) for the calculation of spectral functions, quasiparticle properties, and mobility of the Holstein polaron. We devise an efficient numerical implementation that allows us to make comparisons in a broad interval of temperature, electron-phonon coupling, and phonon frequency. For a benchmark, we use the dynamical mean-field theory which gives, as we have recently shown, rather accurate spectral functions in the whole parameter space, even in low dimensions. We find that in one dimension, the CE resolves well both the quasiparticle and the first satellite peak in a regime of intermediate coupling. At high temperatures, the charge mobility assumes a power law $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}2}$ in the limit of weak coupling and $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}3/2}$ for stronger coupling. We find that, for stronger coupling, the CE gives slightly better results than the self-consistent Migdal approximation (SCMA), while the one-shot Migdal approximation is appropriate only for a very weak electron-phonon interaction. We also analyze the atomic limit and the spectral sum rules. We derive an analytical expression for the moments in CE and find that they are exact up to the fourth order, as opposed to the SCMA where they are exact to the third order. Finally, we analyze the results in higher dimensions.
The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations.The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion.Particular emphasis is laid on the different physics occurring in the adiabatic and antiadiabatic regimes, respectively.Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit.This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.
The phonon spectral function of the one-dimensional Holstein model is obtained within weak and strong-coupling approximations based on analytical self-energy calculations. The characteristic excitations found in the limit of small charge-carrier density are related to the known (electronic) spectral properties of Holstein polarons such as the polaron band dispersion. Particular emphasis is laid on the different physics occurring in the adiabatic and anti-adiabatic regimes, respectively. Comparison is made with a cluster approach exploiting exact numerical results on small systems to yield an approximation for the thermodynamic limit. This method, similar to cluster perturbation theory, confirms the analytical findings, and yields accurate results also in the intermediate-coupling regime.
We develop a hierarchical equations of motion (HEOM) approach to compute real-time single-particle correlation functions and thermodynamic properties of the Holstein model at finite temperature. We exploit the conservation of the total momentum of the system to formulate the momentum-space HEOM whose dynamical variables explicitly keep track of momentum exchanges between the electron and phonons. Our symmetry-adapted HEOM enable us to overcome the numerical instabilities inherent to the commonly used real-space HEOM. The HEOM method is then used to study the spectral function and thermodynamic quantities of chains containing up to ten sites. The HEOM results compare favorably to existing literature. To provide an independent assessment of the HEOM approach and to gain insight into the importance of finite-size effects, we devise a quantum Monte Carlo (QMC) procedure to evaluate finite-temperature single-particle correlation functions in imaginary time and apply it to chains containing up to twenty sites. QMC results reveal that finite-size effects are quite weak, so that the results on 5 to 10-site chains, depending on the parameter regime, are representative of larger systems. A detailed comparison between the HEOM and QMC data place our HEOM method among reliable methods to compute real-time finite-temperature correlation functions in parameter regimes ranging from low- to high-temperature, and weak- to strong-coupling regime.
Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.
The Holstein model is a benchmark model of systems with electron-phonon interaction. However, its electrical transport properties are not yet fully understood. In this work, we performed numerically exact calculations of the imaginary-time current-current correlation function of the Holstein Hamiltonian for a broad range of model parameters. These calculations were performed using a path-integral-based quantum Monte Carlo method. We compared these results with the results obtained under the assumption of conventional band transport, small polaron hopping, and polaron band transport. From this comparison, we identified the regions in parameter space where each of these transport regimes is valid. In some cases when an imaginary-time comparison of current-current correlation functions could not give conclusive results, we complemented it with real-time comparisons or a comparison with the literature data on numerically exact dc mobilities. Overall, we found that the parameter space is almost completely covered by the three mentioned transport regimes.