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

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

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

Published Aug 12, 2025• 3 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/t1g7-r95d

Related papers

Method for obtaining polaron mobility using real and imaginary time path-integral quantum Monte Carlo

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

We developed a path-integral quantum Monte Carlo-based methodology for calculation of polaron mobility in systems with electron-phonon interaction. Within the method, the current-current correlation function in both the imaginary and real time is calculated in a numerically exact way. The choice of basis for representation of the path integral enabled us to reduce the sign problem and perform real-time calculations for longer times. The DC polaron mobility was extracted by performing analytic continuation that makes use of both the real and imaginary-time data. The method was applied to the Holstein polaron model in one dimension. We obtained reliable results for the temperature dependence of the Holstein polaron mobility for interactions ranging from weak to strong and temperatures that are not too low.

Physical review. B./Physical review. B•May 31, 2023•10 citations

Holstein polaron transport from numerically “exact” real-time quantum dynamics simulations

Veljko Janković

Numerically "exact" methods addressing the dynamics of coupled electron-phonon systems have been intensively developed. Nevertheless, the corresponding results for the electron mobility μdc are scarce, even for the one-dimensional (1d) Holstein model. Building on our recent progress on single-particle properties, here we develop the momentum-space hierarchical equations of motion (HEOM) method to evaluate real-time two-particle correlation functions of the 1d Holstein model at a finite temperature. We compute numerically "exact" dynamics of the current-current correlation function up to real times sufficiently long to capture the electron's diffusive motion and provide reliable results for μdc in a wide range of model parameters. In contrast to the smooth ballistic-to-diffusive crossover in the weak-coupling regime, we observe a temporally limited slow-down of the electron on intermediate time scales already in the intermediate-coupling regime, which translates to a finite-frequency peak in the optical response. Our momentum-space formulation lowers the numerical effort with respect to existing HEOM-method implementations, while we remove the numerical instabilities inherent to the undamped-mode HEOM by devising an appropriate hierarchy closing scheme. Still, our HEOM remains unstable at too low temperatures, for too strong electron-phonon coupling, and for too fast phonons.

arXiv:2306.13328•The Journal of Chemical Physics•Sep 6, 2023•17 citations

Spectral Functions of the Holstein Polaron: Exact and Approximate Solutions

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

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.

arXiv:2112.15542•Physical Review Letters•Aug 22, 2022•30 citations

Vertex corrections to conductivity in the Holstein model: A numerical-analytical study

Veljko Janković, Petar Mitrić, D. Tanasković +1 more

The optical-conductivity profile carries information on electronic dynamics in interacting quantum many-body systems. Its computation is a formidable task that is usually approached by invoking the single-particle (bubble) approximation and neglecting vertex corrections, the importance of which remains elusive even in model Hamiltonian calculations. Here, we combine analytical arguments with our recent breakthroughs in numerically exact and approximate calculations of finite-temperature real-time correlation functions to thoroughly assess the importance of vertex corrections in the one-dimensional Holstein polaron model. We find, both analytically and numerically, vanishing vertex corrections to optical conductivity in the limits of zero electron-phonon interaction, zero electronic bandwidth, and infinite temperature. Furthermore, our numerical results show that vertex corrections to the electron mobility also vanish in many parameter regimes between these limits. In some of these cases, the vertex corrections still introduce important qualitative changes to the optical-conductivity profile in comparison to the bubble approximation even though the self-energy remains approximately local. We trace these changes back to the bubble approximation not fully capturing a time-limited slowdown of the electron on intermediate timescales between ballistic and diffusive transport. We find that the vertex corrections are overall most pronounced for intermediate electron-phonon interaction and may increase or decrease the bubble-approximation mobility depending on the values of the model parameters.

arXiv:2403.18394•Physical review. B./Physical review. B•Jun 26, 2024•12 citations

Optical absorption and activated transport in polaronic systems

G. Schubert, Gerhard Wellein, Alexander Weiße +2 more

We present exact results for the optical response in the one-dimensional Holstein model. In particular, by means of a refined kernel polynomial method, we calculate the ac and dc electrical conductivities at finite temperatures for a wide parameter range of electron-phonon interaction. We analyze the deviations from the results of standard small polaron theory in the intermediate coupling regime and discuss nonadiabaticity effects in detail.

arXiv:cond-mat/0505447•Physical Review B•Sep 27, 2005•32 citations

Mobility of Holstein Polaron at Finite Temperature: An Unbiased Approach

Andrei S. Mishchenko, Naoto Nagaosa, G. De Filippis +2 more

We present the first unbiased results for the mobility μ of a one-dimensional Holstein polaron obtained by numerical analytic continuation combined with diagrammatic and worldline Monte Carlo methods in the thermodynamic limit. We have identified for the first time several distinct regimes in the λ-T plane including a band conduction region, incoherent metallic region, an activated hopping region, and a high-temperature saturation region. We observe that although mobilities and mean free paths at different values of λ differ by many orders of magnitude at small temperatures, their values at T larger than the bandwidth become very close to each other.

arXiv:1406.6486•Physical Review Letters•Apr 7, 2015•76 citations
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.
We developed a path-integral quantum Monte Carlo-based methodology for calculation of polaron mobility in systems with electron-phonon interaction. Within the method, the current-current correlation function in both the imaginary and real time is calculated in a numerically exact way. The choice of basis for representation of the path integral enabled us to reduce the sign problem and perform real-time calculations for longer times. The DC polaron mobility was extracted by performing analytic continuation that makes use of both the real and imaginary-time data. The method was applied to the Holstein polaron model in one dimension. We obtained reliable results for the temperature dependence of the Holstein polaron mobility for interactions ranging from weak to strong and temperatures that are not too low.
Numerically "exact" methods addressing the dynamics of coupled electron-phonon systems have been intensively developed. Nevertheless, the corresponding results for the electron mobility μdc are scarce, even for the one-dimensional (1d) Holstein model. Building on our recent progress on single-particle properties, here we develop the momentum-space hierarchical equations of motion (HEOM) method to evaluate real-time two-particle correlation functions of the 1d Holstein model at a finite temperature. We compute numerically "exact" dynamics of the current-current correlation function up to real times sufficiently long to capture the electron's diffusive motion and provide reliable results for μdc in a wide range of model parameters. In contrast to the smooth ballistic-to-diffusive crossover in the weak-coupling regime, we observe a temporally limited slow-down of the electron on intermediate time scales already in the intermediate-coupling regime, which translates to a finite-frequency peak in the optical response. Our momentum-space formulation lowers the numerical effort with respect to existing HEOM-method implementations, while we remove the numerical instabilities inherent to the undamped-mode HEOM by devising an appropriate hierarchy closing scheme. Still, our HEOM remains unstable at too low temperatures, for too strong electron-phonon coupling, and for too fast phonons.
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.
The optical-conductivity profile carries information on electronic dynamics in interacting quantum many-body systems. Its computation is a formidable task that is usually approached by invoking the single-particle (bubble) approximation and neglecting vertex corrections, the importance of which remains elusive even in model Hamiltonian calculations. Here, we combine analytical arguments with our recent breakthroughs in numerically exact and approximate calculations of finite-temperature real-time correlation functions to thoroughly assess the importance of vertex corrections in the one-dimensional Holstein polaron model. We find, both analytically and numerically, vanishing vertex corrections to optical conductivity in the limits of zero electron-phonon interaction, zero electronic bandwidth, and infinite temperature. Furthermore, our numerical results show that vertex corrections to the electron mobility also vanish in many parameter regimes between these limits. In some of these cases, the vertex corrections still introduce important qualitative changes to the optical-conductivity profile in comparison to the bubble approximation even though the self-energy remains approximately local. We trace these changes back to the bubble approximation not fully capturing a time-limited slowdown of the electron on intermediate timescales between ballistic and diffusive transport. We find that the vertex corrections are overall most pronounced for intermediate electron-phonon interaction and may increase or decrease the bubble-approximation mobility depending on the values of the model parameters.
We present the first unbiased results for the mobility μ of a one-dimensional Holstein polaron obtained by numerical analytic continuation combined with diagrammatic and worldline Monte Carlo methods in the thermodynamic limit. We have identified for the first time several distinct regimes in the λ-T plane including a band conduction region, incoherent metallic region, an activated hopping region, and a high-temperature saturation region. We observe that although mobilities and mean free paths at different values of λ differ by many orders of magnitude at small temperatures, their values at T larger than the bandwidth become very close to each other.