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:2306.13328 · The Journal of Chemical Physics

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

Veljko Janković

Published Sep 6, 2023• 17 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/5.0165532

Related papers

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

Charge transport limited by nonlocal electron-phonon interaction. I. Hierarchical equations of motion approach

Veljko Janković

Studying charge transport in models with nonlocal carrier-phonon interaction is difficult because it requires finite-temperature real-time correlation functions of mixed carrier-phonon operators. Focusing on models with discrete undamped phonon modes, we show that such correlation functions can be retrieved from the hierarchical equations of motion (HEOM), although phonons have been integrated out. Our procedure relies on the explicit expression of HEOM auxiliaries in terms of phonon creation and annihilation operators. It reveals that the auxiliaries describe multiphonon-assisted carrier transitions induced by genuine many-phonon correlations, from which lower-order correlations are subtracted according to the finite-temperature Wick's theorem. Applying the procedure to our recently developed momentum-space HEOM method featuring a specific hierarchy closing, we compute the numerically exact dynamical mobility of a carrier within the one-dimensional Peierls model. The carrier mobility at moderate temperatures decreases with increasing interaction, whereas high temperatures see the opposite trend, reflecting the prevalence of the phonon-assisted current over the purely electronic band current. The pronounced finite-size effects and HEOM instabilities delimit the range of applicability of our approach to moderate interactions, moderate to high temperatures, and not too fast phonons. Importantly, this range comprises the values relevant for charge transport in crystalline organic semiconductors, and we present and discuss the corresponding numerically exact results in a companion paper (arXiv:2501.05055).

arXiv:2501.05054•Jan 9, 2025•1 citation

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

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

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

Carrier mobility in Holstein–Peierls models of organic materials: A tensor-train HEOM approach

H. Takahashi, Raffaele Borrelli

We present a comprehensive theoretical study of charge transport in organic materials on the basis of the one-dimensional dispersive Holstein-Peierls model. The hierarchical equations of motion in tensor-train format (HEOM-TT) method is employed to accurately describe both local Holstein-type and nonlocal Peierls-type electron-phonon couplings in the presence of low-frequency dispersive phonons. Linear response formulation of carrier mobility is obtained by employing the generating functional method. We demonstrate the capability of the method by calculating current correlation functions, diffusivities, and mobilities for various transport regimes. Our results provide new insights into the mechanisms underlying charge transport in organic semiconductors and establish HEOM-TT as a powerful tool for studying complex vibronic quantum dynamics and transport properties in condensed matter systems.

The Journal of Chemical Physics•Nov 17, 2025•2 citations
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.
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.
Studying charge transport in models with nonlocal carrier-phonon interaction is difficult because it requires finite-temperature real-time correlation functions of mixed carrier-phonon operators. Focusing on models with discrete undamped phonon modes, we show that such correlation functions can be retrieved from the hierarchical equations of motion (HEOM), although phonons have been integrated out. Our procedure relies on the explicit expression of HEOM auxiliaries in terms of phonon creation and annihilation operators. It reveals that the auxiliaries describe multiphonon-assisted carrier transitions induced by genuine many-phonon correlations, from which lower-order correlations are subtracted according to the finite-temperature Wick's theorem. Applying the procedure to our recently developed momentum-space HEOM method featuring a specific hierarchy closing, we compute the numerically exact dynamical mobility of a carrier within the one-dimensional Peierls model. The carrier mobility at moderate temperatures decreases with increasing interaction, whereas high temperatures see the opposite trend, reflecting the prevalence of the phonon-assisted current over the purely electronic band current. The pronounced finite-size effects and HEOM instabilities delimit the range of applicability of our approach to moderate interactions, moderate to high temperatures, and not too fast phonons. Importantly, this range comprises the values relevant for charge transport in crystalline organic semiconductors, and we present and discuss the corresponding numerically exact results in a companion paper (arXiv:2501.05055).
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.
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.
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 present a comprehensive theoretical study of charge transport in organic materials on the basis of the one-dimensional dispersive Holstein-Peierls model. The hierarchical equations of motion in tensor-train format (HEOM-TT) method is employed to accurately describe both local Holstein-type and nonlocal Peierls-type electron-phonon couplings in the presence of low-frequency dispersive phonons. Linear response formulation of carrier mobility is obtained by employing the generating functional method. We demonstrate the capability of the method by calculating current correlation functions, diffusivities, and mobilities for various transport regimes. Our results provide new insights into the mechanisms underlying charge transport in organic semiconductors and establish HEOM-TT as a powerful tool for studying complex vibronic quantum dynamics and transport properties in condensed matter systems.