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

Diagrammatic quantum Monte Carlo toward the calculation of transport properties in disordered semiconductors

Yu‐Chen Wang, Yi Zhao

Published May 10, 2022• 8 citations
arXivHTMLPDF

Abstract

A new diagrammatic quantum Monte Carlo approach is proposed to deal with the imaginary time propagator involving both dynamic disorder (i.e., electron-phonon interactions) and static disorder of local or nonlocal nature in a unified and numerically exact way. The establishment of the whole framework relies on a general reciprocal-space expression and a generalized Wick's theorem for the static disorder. Since the numerical cost is independent of the system size, various physical quantities, such as the thermally averaged coherence, Matsubara one-particle Green's function, and current autocorrelation function, can be efficiently evaluated in the thermodynamic limit (infinite in the system size). The validity and performance of the proposed approach are systematically examined in a broad parameter regime. This approach, combined with proper numerical analytic continuation methods and first-principles calculations, is expected to be a versatile tool toward the calculation of various transport properties, such as mobilities in realistic semiconductors involving multiple electronic energy bands, high-frequency optical and low-frequency acoustic phonons, different forms of dynamic and static disorders, and anisotropy.

DOI
10.1063/5.0091124

Related papers

A stochastic Schrödinger equation and matrix product state approach to carrier transport in organic semiconductors with nonlocal electron–phonon interaction

Liqi Zhou, Xing Bao Gao, Zhigang Shuai

Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward-backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward-backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron-phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size.

The Journal of Chemical Physics•Aug 28, 2024•12 citations

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

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

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

Time correlation functions for quantum systems: the case of harmonic oscillators

Vladislav Efremkin, Jean‐Louis Barrat, Stefano Mossa +1 more

The quantum harmonic oscillator represents the most common fundamental building block to compute thermal properties of virtually any dielectric crystalline material at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo (PIMC) methods provide a powerful numerical tool to stochastically determine thermodynamic quantities without systematic bias, not relying on perturbative schemes. Addressing transport properties, for instance the calculation of thermal conductivity from PIMC, however, has turned out to be substantially more difficult. Although correlation functions of current operators can be determined by PIMC in the form of analytic continuation on the imaginary time axis, Bayesian methods are usually employed for the numerical inversion back to real-time response functions. This task not only relies strongly on the accuracy of the PIMC data, but also introduces noticeable dependence on the underlying model used for the inversion problem. Here, we address both difficulties with great care. In particular, we first devise improved estimators for current correlations which substantially reduce the variance of the PIMC data. Next, we provide a neat statistical approach to the inversion problem, blending into a fresh workflow the classical stochastic maximum entropy method together with more recent notions borrowed from statistical learning theory. We test our ideas on a single harmonic oscillator and a collection of oscillators with a continuous distribution of frequencies, establishing solid grounds for an unbiased, fully quantum mechanical calculation of transport properties in solids.

•May 20, 2021•0 citations

Real-frequency diagrammatic Monte Carlo at finite temperature

J. Vučičević, Michel Ferrero

Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [A. Taheridehkordi et al., Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $\mathrm{\ensuremath{\Sigma}}(\ensuremath{\omega})$ of the doped $32\ifmmode\times\else\texttimes\fi{}32$ cyclic square-lattice Hubbard model in a nontrivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.

arXiv:1908.11826•Physical review. B./Physical review. B•Feb 12, 2020•35 citations
A new diagrammatic quantum Monte Carlo approach is proposed to deal with the imaginary time propagator involving both dynamic disorder (i.e., electron-phonon interactions) and static disorder of local or nonlocal nature in a unified and numerically exact way. The establishment of the whole framework relies on a general reciprocal-space expression and a generalized Wick's theorem for the static disorder. Since the numerical cost is independent of the system size, various physical quantities, such as the thermally averaged coherence, Matsubara one-particle Green's function, and current autocorrelation function, can be efficiently evaluated in the thermodynamic limit (infinite in the system size). The validity and performance of the proposed approach are systematically examined in a broad parameter regime. This approach, combined with proper numerical analytic continuation methods and first-principles calculations, is expected to be a versatile tool toward the calculation of various transport properties, such as mobilities in realistic semiconductors involving multiple electronic energy bands, high-frequency optical and low-frequency acoustic phonons, different forms of dynamic and static disorders, and anisotropy.
Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward-backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward-backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron-phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size.
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.
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.
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 quantum harmonic oscillator represents the most common fundamental building block to compute thermal properties of virtually any dielectric crystalline material at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo (PIMC) methods provide a powerful numerical tool to stochastically determine thermodynamic quantities without systematic bias, not relying on perturbative schemes. Addressing transport properties, for instance the calculation of thermal conductivity from PIMC, however, has turned out to be substantially more difficult. Although correlation functions of current operators can be determined by PIMC in the form of analytic continuation on the imaginary time axis, Bayesian methods are usually employed for the numerical inversion back to real-time response functions. This task not only relies strongly on the accuracy of the PIMC data, but also introduces noticeable dependence on the underlying model used for the inversion problem. Here, we address both difficulties with great care. In particular, we first devise improved estimators for current correlations which substantially reduce the variance of the PIMC data. Next, we provide a neat statistical approach to the inversion problem, blending into a fresh workflow the classical stochastic maximum entropy method together with more recent notions borrowed from statistical learning theory. We test our ideas on a single harmonic oscillator and a collection of oscillators with a continuous distribution of frequencies, establishing solid grounds for an unbiased, fully quantum mechanical calculation of transport properties in solids.
Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [A. Taheridehkordi et al., Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $\mathrm{\ensuremath{\Sigma}}(\ensuremath{\omega})$ of the doped $32\ifmmode\times\else\texttimes\fi{}32$ cyclic square-lattice Hubbard model in a nontrivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.