Pith. sign in

REVIEW 3 major objections 3 minor 2 cited by

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

T0 review · 3 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read HEOM auxiliaries encode phonon-assisted transport in the 1D Peierls model, giving numerically exact mobility.

desk verdict The formalism is the real contribution; the mobility numbers are honest but not literally 'numerically exact' because the hierarchy closing is uncontrolled. read the letter →

arxiv 2501.05054 v3 pith:4EI5UZ47 submitted 2025-01-09 cond-mat.str-el cond-mat.mes-hallphysics.chem-phquant-ph

Veljko Janković This is my paper · ORCID
classification cond-mat.str-elcond-mat.mes-hallphysics.chem-phquant-ph
keywords hierarchicalequationsofmotionnonlocalelectron-phononinteractionPeierlsmodelcarriermobilityphonon-assistedcurrentgeneralizedWick'stheoremorganicsemiconductorsquantumtransport
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper is about recovering the full quantum carrier–phonon dynamics from the hierarchical equations of motion (HEOM), even though phonons have been integrated out and only electronic auxiliary operators appear. It claims that, for harmonic undamped phonons with linear carrier–phonon coupling, each HEOM auxiliary operator has an explicit many-phonon form, so the auxiliaries describe genuine n-phonon-assisted carrier transitions with lower-order correlations subtracted. From that representation the paper proves the generalized Wick's theorem and uses it to compute the finite-temperature current–current correlation function including the phonon-assisted current, which previous HEOM treatments of Peierls-type models could not handle. Applying the framework to the one-dimensional Peierls model gives numerically exact carrier mobility over a parameter range: mobility decreases with interaction at moderate temperatures but increases at high temperatures, marking the crossover from band to phonon-assisted transport. The approach is practical for moderate interactions, moderate to high temperatures, and not too fast phonons, which is the range thought to be relevant for crystalline organic semiconductors.

What carries the argument

The load-bearing machinery is the explicit phononic representation of the HEOM auxiliary operators together with the generalized Wick's theorem. The auxiliary $F^{(n)}_{\mathbf{n}}$ is built from unit operators $f_{qm}\propto b_q$, $b_q^\dagger$ as a normal-ordered product with all lower-order equilibrium contractions subtracted (Eqs. 10 and 11), and the theorem $F^{(n)}_{\mathbf{n}} f_{qm} = F^{(n+1)}_{\mathbf{n}^+_{qm}} + \sum_{\mathbf{q}'m'} n_{\mathbf{q}'m'}\eta_{\mathbf{q}'\mathbf{q}m'} F^{(n-1)}_{\mathbf{n}^-_{\mathbf{q}'m'}}$ (and its mirror) is proved rather than assumed. That identity lets hybrid electron–phonon correlations be read off from purely electronic HEOM tiers. The computational machinery is the momentum-space HEOM with a Markovian–adiabatic hierarchy closing at maximum depth $D$, which is the main approximation of the framework; the random-phase approximation is used to simplify the closing term. The current–current correlation is assembled from the band, phonon-assisted, and cross contributions (Eqs. 30–33).

What would settle it

Take the one-dimensional Peierls model at $J=\omega_0=1$, $\lambda=0.25$, $T=1$, where the paper reports that its closing scheme struggles, and compute the zero-frequency mobility with an independent numerically exact method that does not rely on hierarchy truncation (for example, a small-chain calculation with converged phonon number per site, suitably extrapolated). If that mobility disagrees with the HEOM result beyond the stated uncertainty of about 10%, the closing approximation is the cause; agreement would confirm the auxiliary representation and the generalized Wick theorem.

Watch

Extended reading notes

Core claim

The central discovery is an operator identity for the phonon degrees of freedom. The HEOM auxiliary $F^{(n)}_{\mathbf{n}}$ equals a normally ordered product of phonon creation and annihilation operators minus all finite-temperature contractions down to lower orders; in equilibrium its average vanishes for $n>0$, so the auxiliary isolates genuine many-phonon correlations. Using this explicit expression (Eq. 10), the paper proves the generalized Wick's theorem: multiplying $F^{(n)}_{\mathbf{n}}$ by a phonon operator $f_{qm}$ produces the next-tier auxiliary plus lower-tier auxiliaries weighted by equilibrium phonon contractions. That identity is the missing link that lets the paper write the initial condition for the phonon-assisted part of the current operator in terms of HEOM auxiliaries (Eq. 28) and obtain the current–current correlation $C_{jj}(t)$ from the zeroth and first tiers (Eq. 29). With this, the paper computes numerically exact dynamical mobility of the one-dimensional Peierls model and identifies the regime in which the phonon-assisted current overtakes the band current.

Load-bearing premise

The load-bearing assumption is that the approximate formula used to terminate the hierarchy at maximum depth $D$ correctly captures what all deeper layers do to the mobility; the paper validates this only by comparing two closing schemes and by convergence checks.

Editorial extensions

If this is right

  • The phonon-assisted current, previously a serious obstacle for HEOM, is now computable on the same footing as the band current, so mobility, optical conductivity, and frequency-dependent response can be separated into their physical channels.
  • In the one-dimensional Peierls model, increasing carrier–phonon interaction lowers mobility at moderate temperatures but raises it at high temperatures because the phonon-assisted channel overtakes the band channel; the crossover sits near $T\simeq 2J$.
  • The method covers the moderate-interaction, moderate-to-high-temperature, not-too-fast-phonon regime, which is the regime believed relevant for charge transport in crystalline organic semiconductors.
  • For weak interaction ($\lambda=0.05$), HEOM mobility agrees with self-consistent semiclassical transport theory up to $T\simeq 5$, after which phonon-assisted and cross contributions become sizable and the semiclassical result underestimates the mobility.
  • Since the auxiliary representation only assumes harmonic undamped phonons and linear coupling, the same construction applies to local (Holstein-type) and multi-mode models, not only to the Peierls model.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the auxiliary representation and the generalized Wick theorem hold generally, the same HEOM machinery can produce other mixed correlation functions, such as $\langle j(t) B_q(0)\rangle$ or phonon-occupation response, opening a window into nonequilibrium phonon dynamics that this paper does not exploit.
  • The reported failure of the closing scheme at low temperature likely reflects the Markovian–adiabatic closure itself rather than the auxiliary representation, so a better closure could shift the practical boundary toward $T/\omega_0<2$ and slow phonons.
  • The paper's high-temperature comparison suggests a quantitative target: any approximate theory of organic-semiconductor transport should reproduce the power-law exponent $\alpha\simeq 0.5$ for the mobility's temperature decay and the weaker-than-linear dependence on interaction strength.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The manuscript develops a hierarchical equations of motion (HEOM) framework for computing real-time finite-temperature current–current correlation functions in models with nonlocal electron–phonon coupling and discrete undamped phonons. The central formal contributions are an explicit many-phonon representation of the HEOM auxiliary operators (Eq. 10) and a proof of the generalized Wick's theorem (Eqs. 12–15). These are used to express the phonon-assisted current and to compute the dynamical mobility in the one-dimensional Peierls model. Numerical results are reported for the mobility as a function of temperature and interaction strength, with a discussion of the hierarchy closing approximation and its limitations. The paper explicitly states that the method is applicable at moderate to high temperatures, moderate interactions, and not too fast phonons, and it presents results in that regime.

Significance. If the results are correct, the paper provides a much-needed reference method for transport in Peierls-type models, where fully quantum calculations have been largely inaccessible. The formal derivation of the generalized Wick's theorem is a genuine advance in connecting HEOM auxiliaries with many-phonon processes, and the open data and detailed convergence checks are commendable strengths. However, the 'numerically exact' claim is tempered by the uncontrolled hierarchy closing, which is the main approximation. The comparison with Boltzmann theory and the high-temperature formula (Eq. 59) provides useful physical context and partially validates the method in the weak-coupling regime. Overall, the formal results are valuable, while the numerical claims require more careful qualification or additional error control.

major comments (3)
  1. [Sec. IV B and Eq. (46)] The hierarchy closing is the main approximation, and its error is not controlled. The MA and DR schemes both rely on the random-phase approximation (Appendix D), so their agreement (Sec. V D) cannot exclude a common bias. The abstract and Sec. VI describe the results as 'numerically exact,' but this is not justified without an error estimate for the closing. I recommend either providing an independent validation (e.g., comparison with a different method for a limited parameter set) or rephrasing the claims to 'numerically converged with respect to hierarchy depth within the closing approximation.' This is load-bearing because the reported mobility values and trends depend on the closing.
  2. [Sec. V C and Fig. 2] At T=1, lambda=0.25, the time-dependent diffusion constant does not saturate for N=21 and D=4–6, and the mobility is obtained by fitting to an exponential saturation over a window 30 ≤ Jt ≤ 80 for N=45, with a stated ~10% uncertainty. The T=1 entries in Fig. 4(a) are flagged as possibly not fully converged. Given this, the use of these points to support the temperature and interaction trends (e.g., the crossover around T/J=2) is not fully robust. The paper should present these points with larger error bars or discuss them as preliminary estimates.
  3. [Sec. V D and Eq. (52)] The time-reversal check shows that the maximum of the LHS of Eq. (52) is about 10^-2 for the D values used in converged calculations. This is not negligible compared to the 10% uncertainty in mobility, and the paper should comment on whether this magnitude is small enough to ensure the accuracy of the individual contributions, especially the cross contribution, at the reported level.
minor comments (3)
  1. [Abstract and Sec. VI] The abstract and Sec. VI use 'numerically exact' without qualification, while Sec. V C states that T=1 results 'may not be fully converged.' This inconsistency should be resolved by adding a qualifier such as 'within the hierarchy-closing approximation' or by softening the claim.
  2. [Sec. I] The phrase 'numerically exact' is used in the abstract and elsewhere, but the paper later notes the results are exact only 'for all practical purposes' (Sec. V B). This imprecise usage should be harmonized.
  3. [Fig. 2 caption] The caption does not define the 'fit' curve clearly; it would be helpful to state explicitly that the fit is to an exponential saturation function with parameters given in the text.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the formal HEOM-phonon mapping and the generalized Wick theorem are derived from first principles, and the mobility computation is self-contained with external benchmarks.

full rationale

The derivation chain is self-contained. The central formal result, the explicit phonon-operator representation of HEOM auxiliaries (Eq. 10), is derived in Appendix B from the Feynman-Vernon influence-functional definition of the auxiliaries (Eqs. A4-A9) using the finite-temperature Wick theorem for free phonons. The generalized Wick theorem (Eqs. 12-15) is then proven in Appendix C by direct algebraic manipulation of Eq. 10; it is not assumed as an input. The transport framework in Sec. III uses this theorem to evaluate the phonon-assisted initial conditions (Eq. 28) and to extract the current-current correlation function (Eq. 29); no fitted parameter is renamed as a prediction. The only nontrivial approximation is the hierarchy closing scheme (Sec. IV B), which the paper explicitly identifies as 'the main approximation of our framework.' This is a numerical-convergence and accuracy risk, not a circularity: the closing is derived in Appendix D, compared against the TNL and DR schemes in Sec. V D, and benchmarked against Boltzmann theory in the weak-coupling regime (Fig. 5). Self-citations of the author's previous HEOM work (Refs. 60-62) supply the momentum-space HEOM machinery and the MA closing idea, but the closing is re-derived here and validated against independent schemes and external Boltzmann results, so those self-citations are not load-bearing. The high-temperature comparison to the analytic expression Eq. 59 is an independent approximate benchmark, not a fit target. No load-bearing step reduces by construction to its own input, so the correct finding is no significant circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are fitted to the target data; the model parameters J, omega0, g, T are inputs. The hierarchy closing uses a second-order scattering rate, not a fitted constant. No new physical entities are introduced; the F(n) operators are mathematical constructs within the existing HEOM/DEOM framework.

assumptions (5)
  • domain assumption Phonons are undamped harmonic oscillators with Gaussian statistics, so the finite-temperature Wick theorem applies.
    Invoked in Sec. II B and Appendix B to express auxiliaries through phonon operators and to derive Eqs. (10)-(11).
  • domain assumption The carrier-phonon interaction is linear in both phonon displacements and single-carrier densities (Eqs. 1-3).
    The HEOM and the auxiliary representation rely on this linear coupling; the author notes generalization but derives for this case.
  • domain assumption Phonons remain in thermal equilibrium, so equilibrium expectation values subtract lower-order correlations in Eq. (11).
    The paper states this 'effectively assuming that phonons are in thermal equilibrium' in Sec. II B.
  • ad hoc to paper The Markovian and adiabatic approximations plus random-phase approximation for the hierarchy closing (Eqs. 46-49, Appendix D) are valid.
    The closing scheme is the main approximation of the framework; it is validated only by comparison of MA and DR schemes and by convergence checks.
  • domain assumption The 1D Peierls model with a single carrier, periodic boundary conditions, and dispersion epsilon_k = -2J cos k is the physical model.
    Defines the model in Sec. III; all transport results pertain to this model.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Charge transport limited by nonlocal electron-phonon interaction. I. Hierarchical equations of motion approach." pith.science (2026). https://pith.science/paper/4EI5UZ47

@misc{pith2026250105054,
  author       = {Pith},
  title        = {Pith review of: Charge transport limited by nonlocal electron-phonon interaction. I. Hierarchical equations of motion approach},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4EI5UZ47}},
  note         = {Machine review of arXiv:2501.05054}
}
read the original abstract

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 general 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 appropriate 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).

Figures

Figures reproduced from arXiv: 2501.05054 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Time-dependent diffusion constant [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) and (b): Time-dependent diffusion constant for [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Short-time HEOM dynamics of Re [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Temperature dependence of (a) the carrier mobility [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Temperature-dependent charge mobility computed [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Temperature dependence of (a) the carrier mobility [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Purely electronic (dashed lines, label ”e”), phonon [PITH_FULL_IMAGE:figures/full_fig_p015_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Temperature-dependent mobility computed using the [PITH_FULL_IMAGE:figures/full_fig_p025_8.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Charge transport limited by nonlocal electron-phonon interaction. II. Numerically exact quantum dynamics in the slow-phonon regime

    cond-mat.str-el 2025-01 conditional novelty 6.0 of 10

    In the slow-phonon one-dimensional Peierls model, carrier diffusion is achieved from the superdiffusive side after a transient subdiffusive phase, and the mobility shows both a zero-frequency peak and a dip at the pho...

  2. Dynamical quantum typicality: Simple method for investigating transport properties applied to the Holstein model

    cond-mat.str-el 2024-12 conditional novelty 6.0 of 10

    Quantum typicality produces well-converged frequency-dependent electron mobility for the Holstein model at strong coupling, confirmed by quantum Monte Carlo and used to quantify vertex corrections.

Reference graph

Works this paper leans on

114 extracted references · 73 canonical work pages · cited by 2 Pith papers

  1. [1]

    Markovian and adiabatic scheme We developed and tested this scheme on the one-dimensional Holste in model in Ref. 60. It transforms Eq. ( D2) by neglecting the hierarchical couplings to auxiliaries at depth D + 2 [the second term on the RHS of Eq. ( D2)] and 22 retaining only the coupling to ρ(D) D (t) for which Eq. ( D1) is formulated [the summand with q...

  2. [2]

    Derivative-resum scheme In Eq. ( D2), we assume that [91] ∂tρ(D+1) D+ qm (t) j 2i ∑ q2m2 √ 1 + Dq2m2 + δq2qδm2m √ |cm2 |V × q2 ρ(D+2) D+, + qm,q 2 m2 (t), (D7) 23 and retain only the coupling to ρ(D) D (t) [the summand with q2 = q and m2 = m in the third term on the RHS of Eq. ( D2)] [92]. We thus obtain ρ(D+1) D+ qm (t) = 2 √ 1 + Dqm √ |cm| [Le 2 i(µ D +...

  3. [3]

    Nematiaram and A

    T. Nematiaram and A. Troisi, Modeling charge transport in high-mobility molecular semiconduc- tors: Balancing electronic structure and quantum dynamics methods with the help of experiments, J. Chem. Phys. 152, 190902 (2020)

    2020

  4. [4]

    Giannini and J

    S. Giannini and J. Blumberger, Charge trans- port in organic semiconductors: The perspec- tive from nonadiabatic molecular dynamics, Acc. Chem. Res. 55, 819 (2022)

    2022

  5. [5]

    Therefore, we give preference to the MA over the DR closing scheme at lower temperatures. Overall, we conclude that the MA closing scheme sta- bilizes HEOM dynamics in a manner that does not in- troduce spurious high-frequency features and does not compromise low-frequency features of the dynamical mo- bility, and, in particular, the magnitude of the carr...

  6. [6]

    M. J. Schilcher, P. J. Robinson, D. J. Abramovitch, L. Z. Tan, A. M. Rappe, D. R. Reichman, and D. A. Egger, The significance of polarons and dynamic disorder in halide perovskites, ACS Energy Lett. 6, 2162 (2021)

    2021

  7. [7]

    and ( 10)]. The auxiliaries are found to describe many phonon-assisted transitions between free-carrier states that are mediated by genuine many-phonon cor- relations, from which the redundant information already present in lower-order auxiliaries is eliminated [Eq. ( 11)]. We then rigorously prove the generalized Wick’s theo- rem [Eqs. (

  8. [8]

    Ponc´ e, W

    S. Ponc´ e, W. Li, S. Reichardt, and F. Giustino, First-principles calculations of charge car- rier mobility and conductivity in bulk semi- conductors and two-dimensional materials, Rep. Prog. Phys. 83, 036501 (2020)

    2020

Show all 114 references
  1. [9]

    Franchini, M

    C. Franchini, M. Reticcioli, M. Setvin, and U. Diebold, Polarons in materials, Nat. Rev. Mater. 6, 560 (2021)

    2021

  2. [10]

    Holstein, Studies of polaron motion: Part I

    T. Holstein, Studies of polaron motion: Part I. The molecular-crystal model, Ann. Phys. 8, 325 (1959)

    1959

  3. [11]

    Fratini, M

    S. Fratini, M. Nikolka, A. Salleo, G. Schweicher, and H. Sirringhaus, Charge transport in high-mobility conjugated polymers and molecular semiconductors, Nat. Mater. 19, 491 (2020)

    2020

  4. [12]

    A. S. Alexandrov and J. T. Devreese, Advances in Polaron Physics (Springer-Verlag Berlin Heidelberg, 2010)

    2010

  5. [13]

    L. R. V. Buizza and L. M. Herz, Polarons and charge localization in metal-halide semicon- ductors for photovoltaic and light-emitting devices, Adv. Mater. 33, 2007057 (2021)

    2021

  6. [14]

    D. A. Micha and I. Burghardt, eds., Quantum Dy- namics of Complex Molecular Systems , Springer Series in Chemical Physics (CHEMICAL, vol. 83) (Springer- Verlag Berlin Heidelberg, 2007)

    2007

  7. [15]

    III) that handles finite-temperature correlation functions of mixed carrier– phonon operators

    and ( 14)], which is the essential ingredient of the computational framework (Sec. III) that handles finite-temperature correlation functions of mixed carrier– phonon operators. This framework is then used to obtain numerically ex- act transport properties of the one-dimensiona...

  8. [16]

    Nematiaram, D

    T. Nematiaram, D. Padula, A. Landi, and A. Troisi, On the largest possible mobility of molecular semiconductors and how to achieve it, Adv. Funct. Mater. 30, 2001906 (2020)

    2020

  9. [17]

    Mahan, Many-Particle Physics (Kluwer Academic, New York, 2000)

    G. Mahan, Many-Particle Physics (Kluwer Academic, New York, 2000)

    2000

  10. [18]

    L. Wang, O. V. Prezhdo, and D. Beljonne, Mixed quantum-classical dynamics for charge transport in or- ganics, Phys. Chem. Chem. Phys. 17, 12395 (2015)

    2015

  11. [19]

    W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in polyacetylene, Phys. Rev. Lett. 42, 1698 (1979)

    1979

  12. [20]

    Capone, W

    M. Capone, W. Stephan, and M. Grilli, Small- 26 polaron formation and optical absorption in Su-Schrieffer-Heeger and Holstein models, Phys. Rev. B 56, 4484 (1997)

    1997

  13. [21]

    D. J. J. Marchand, G. De Filippis, V. Cataudella, M. Berciu, N. Nagaosa, N. V. Prokof’ev, A. S. Mishchenko, and P. C. E. Stamp, Sharp Transition for Single Polarons in the One-Dimensional Su-Schrieffer- Heeger Model, Phys. Rev. Lett. 105, 266605 (2010)

    2010

  14. [22]

    Troisi and G

    A. Troisi and G. Orlandi, Charge-Transport Regime of Crystalline Organic Semiconductors: Diffusion Lim- ited by Thermal Off-Diagonal Electronic Disorder, Phys. Rev. Lett. 96, 086601 (2006)

    2006

  15. [23]

    L. Wang, A. V. Akimov, L. Chen, and O. V. Prezhdo, Quantized Hamiltonian dynamics captures the low- temperature regime of charge transport in molecular crystals, J. Chem. Phys. 139, 174109 (2013)

    2013

  16. [24]

    Rammal, A

    H. Rammal, A. Ralko, S. Ciuchi, and S. Fratini, Tran- sient localization from the interaction with quantum bosons, Phys. Rev. Lett. 132, 266502 (2024)

    2024

  17. [25]

    J. E. Runeson, T. J. G. Drayton, and D. E. Manolopoulos, Charge transport in organic semiconduc- tors from the mapping approach to surface hopping, J. Chem. Phys. 161, 144102 (2024)

    2024

  18. [26]

    Ciuchi, S

    S. Ciuchi, S. Fratini, and D. Mayou, Transient localization in crystalline organic semiconductors, Phys. Rev. B 83, 081202 (2011)

    2011

  19. [27]

    Ciuchi and S

    S. Ciuchi and S. Fratini, Electronic transport and quantum localization effects in organic semiconductors, Phys. Rev. B 86, 245201 (2012)

    2012

  20. [28]

    Fratini, D

    S. Fratini, D. Mayou, and S. Ciuchi, The transient lo- calization scenario for charge transport in crystalline or - ganic materials, Adv. Funct. Mater. 26, 2292 (2016)

    2016

  21. [29]

    Fratini and S

    S. Fratini and S. Ciuchi, Dynamical lo- calization corrections to band transport, Phys. Rev. Res. 2, 013001 (2020)

    2020

  22. [30]

    Hutsch, M

    S. Hutsch, M. Panhans, and F. Ortmann, Time- consistent hopping transport with vibration- mode-resolved electron-phonon couplings, Phys. Rev. B 104, 054306 (2021)

    2021

  23. [31]

    Prodanovi´ c and N

    N. Prodanovi´ c and N. Vukmirovi´ c, Charge carrier mo- bility in systems with local electron-phonon interaction, Phys. Rev. B 99, 104304 (2019)

    2019

  24. [32]

    J. H. Fetherolf, D. Goleˇ z, and T. C. Berkelbach, A unification of the Holstein polaron and dynamic dis- order pictures of charge transport in organic crystals, Phys. Rev. X 10, 021062 (2020)

    2020

  25. [33]

    Hannewald and P

    K. Hannewald and P. A. Bobbert, Anisotropy effects in phonon-assisted charge-carrier transport in organic molecular crystals, Phys. Rev. B 69, 075212 (2004)

    2004

  26. [34]

    Cheng and R

    Y.-C. Cheng and R. J. Silbey, A unified the- ory for charge-carrier transport in organic crystals, J. Chem. Phys. 128, 114713 (2008)

    2008

  27. [35]

    Ortmann, F

    F. Ortmann, F. Bechstedt, and K. Hannewald, Theory of charge transport in organic crys- tals: Beyond Holstein’s small-polaron model, Phys. Rev. B 79, 235206 (2009)

    2009

  28. [36]

    Mitri´ c, V

    P. Mitri´ c, V. Jankovi´ c, N. Vukmirovi´ c, and D. Tanaskovi´ c, Cumulant expansion in the Hol- stein model: Spectral functions and mobility, Phys. Rev. B 107, 125165 (2023)

    2023

  29. [37]

    Merkel, M

    K. Merkel, M. Panhans, S. Hutsch, and F. Ort- mann, Interplay of band occupation, localization, and polaron renormalization for electron transport in molecular crystals: Naphthalene as a case study, Phys. Rev. B 105, 165136 (2022)

    2022

  30. [38]

    Kornjaˇ ca and N

    M. Kornjaˇ ca and N. Vukmirovi´ c, Polaron mobility ob- tained by a variational approach for lattice Fr¨ ohlich models, Ann. Phys. 391, 183 (2018)

    2018

  31. [39]

    D. Chen, J. Ye, H. Zhang, and Y. Zhao, On the Munn– Silbey approach to polaron transport with off-diagonal coupling and temperature-dependent canonical trans- formations, J. Phys. Chem. B 115, 5312 (2011)

    2011

  32. [40]

    Mitri´ c, V

    P. Mitri´ c, V. Dobrosavljevi´ c, and D. Tanaskovi´ c, Pre- cursors to Anderson Localization in the Holstein Model: Quantum and Quantum-Classical Solutions (2024), arXiv:2412.17770 [cond-mat.str-el]

    2024 arXiv

  33. [41]

    Fratini and S

    S. Fratini and S. Ciuchi, Bandlike motion and mo- bility saturation in organic molecular semiconductors, Phys. Rev. Lett. 103, 266601 (2009)

    2009

  34. [42]

    Song and Q

    L. Song and Q. Shi, A new approach to calculate charge carrier transport mobility in organic molecular crystals from imaginary time path integral simulations, J. Chem. Phys. 142, 174103 (2015)

    2015

  35. [43]

    Fratini and S

    S. Fratini and S. Ciuchi, Dynamical mean- field theory of transport of small polarons, Phys. Rev. Lett. 91, 256403 (2003)

    2003

  36. [44]

    Fratini and S

    S. Fratini and S. Ciuchi, Optical properties of small polarons from dynamical mean-field theory, Phys. Rev. B 74, 075101 (2006)

    2006

  37. [45]

    Landi, R

    A. Landi, R. Borrelli, A. Capobianco, A. Ve- lardo, and A. Peluso, Hole hopping rates in organic semiconductors: A second-order cumulant approach, J. Chem. Theory Comput. 14, 1594 (2018)

    2018

  38. [46]

    Landi, R

    A. Landi, R. Borrelli, A. Capobianco, A. Velardo, and A. Peluso, Second-order cumulant approach for the eval- uation of anisotropic hole mobility in organic semicon- ductors, J. Phys. Chem. C 122, 25849 (2018)

    2018

  39. [47]

    Landi, Charge mobility prediction in organic semiconductors: Comparison of second-order cumu- lant approximation and transient localization theory, J

    A. Landi, Charge mobility prediction in organic semiconductors: Comparison of second-order cumu- lant approximation and transient localization theory, J. Phys. Chem. C 123, 18804 (2019)

    2019

  40. [48]

    Feynman and F

    R. Feynman and F. Vernon, The theory of a general quantum system interacting with a linear dissipative system, Ann. Phys. 24, 118 (1963)

    1963

  41. [49]

    Balzer, T

    D. Balzer, T. J. A. M. Smolders, D. Blyth, S. N. Hood, and I. Kassal, Delocalised kinetic Monte Carlo for simulating delocalisation-enhanced charge and exciton transport in disordered materials, Chem. Sci. 12, 2276 (2021)

    2021

  42. [50]

    J. T. Willson, W. Liu, D. Balzer, and I. Kassal, Jumping kinetic Monte Carlo: Fast and accurate simulations of partially delocalized charge transport in organic semi- conductors, J. Phys. Chem. Lett. 14, 3757 (2023)

    2023

  43. [51]

    Tanimura, Stochastic Liouville, Langevin, Fokker–Planck, and master equation ap- proaches to quantum dissipative systems, J

    Y. Tanimura, Stochastic Liouville, Langevin, Fokker–Planck, and master equation ap- proaches to quantum dissipative systems, J. Phys. Soc. Jpn. 75, 082001 (2006)

    2006

  44. [52]

    Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), J

    Y. Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), J. Chem. Phys. 153, 020901 (2020)

    2020

  45. [53]

    Xu and Y

    R.-X. Xu and Y. J. Yan, Dynamics of quantum dissipation systems interacting with bosonic canoni- cal bath: Hierarchical equations of motion approach, Phys. Rev. E 75, 031107 (2007)

    2007

  46. [54]

    I. S. Dunn, R. Tempelaar, and D. R. Reichman, Re- moving instabilities in the hierarchical equations of mo- tion: Exact and approximate projection approaches, J. Chem. Phys. 150, 184109 (2019)

    2019

  47. [55]

    Ishizaki and G

    A. Ishizaki and G. R. Fleming, Theoretical ex- 27 amination of quantum coherence in a photo- synthetic system at physiological temperature, Proc. Natl. Acad. Sci. 106, 17255 (2009)

    2009

  48. [56]

    Kreisbeck, T

    C. Kreisbeck, T. Kramer, M. Rodr ´ ıguez, and B. Hein, High-performance solution of hierarchical equations of motion for studying energy transfer in light-harvesting complexes, J. Chem. Theory Comput. 7, 2166 (2011)

    2011

  49. [57]

    Str¨ umpfer and K

    J. Str¨ umpfer and K. Schulten, Open quan- tum dynamics calculations with the hierar- chy equations of motion on parallel computers, J. Chem. Theory Comput. 8, 2808 (2012)

    2012

  50. [58]

    D. M. Wilkins and N. S. Dattani, Why Quantum Co- herence Is Not Important in the Fenna–Matthews–Olsen Complex, J. Chem. Theory Comput. 11, 3411 (2015)

    2015

  51. [59]

    For λ = 0

    [22]. For λ = 0 . 05, we expect that the mobility within the Boltzmann transport theory [11, 100, 101], which consid- ers only the purely electronic contribution to µ dc, should closely follow HEOM results at least at lower tempera- tures featuring small phonon-assisted and cr...

  52. [60]

    Kramer, M

    T. Kramer, M. Noack, A. Reinefeld, M. Rodr ´ ıguez, and Y. Zelinskyy, Efficient calculation of open quantum system dynamics and time-resolved spec- troscopy with distributed memory HEOM (DM- HEOM), J. Comput. Chem. 39, 1779 (2018)

    2018

  53. [61]

    Y. Yan, T. Xing, and Q. Shi, A new method to im- prove the numerical stability of the hierarchical equa- tions of motion for discrete harmonic oscillator modes, J. Chem. Phys. 153, 204109 (2020)

    2020

  54. [62]

    T. Li, Y. Yan, and Q. Shi, Is there a finite mobility for the one vibrational mode Holstein model? Implications from real time simulations, J. Chem. Phys. 160, 111102 (2024)

    2024

  55. [63]

    R. Kubo, M. Toda, and N. Hashitsume, Statisti- cal Physics II: Nonequilibrium Statistical Mechanics (Springer-Verlag Berlin Heidelberg, 1985)

    1985

  56. [64]

    Gosar and S

    P. Gosar and S. Choi, Linear-response theory of the electron mobility in molecular crystals, Phys. Rev. 150, 529 (1966)

    1966

  57. [65]

    Gosar and I

    P. Gosar and I. Vilfan, Phonon-assisted current in or- ganic molecular crystals, Mol. Phys. 18, 49 (1970)

    1970

  58. [66]

    Jankovi´ c, Holstein polaron transport from numeri- cally “exact” real-time quantum dynamics simulations, J

    V. Jankovi´ c, Holstein polaron transport from numeri- cally “exact” real-time quantum dynamics simulations, J. Chem. Phys. 159, 094113 (2023)

    2023

  59. [67]

    Jankovi´ c, P

    V. Jankovi´ c, P. Mitri´ c, D. Tanaskovi´ c, and N. Vuk- mirovi´ c, Vertex corrections to conductivity in the Holstein model: A numerical-analytical study, Phys. Rev. B 109, 214312 (2024)

    2024

  60. [68]

    Jankovi´ c and N

    V. Jankovi´ c and N. Vukmirovi´ c, Spectral and thermodynamic properties of the Holstein po- laron: Hierarchical equations of motion approach, Phys. Rev. B 105, 054311 (2022)

    2022

  61. [69]

    Song and Q

    L. Song and Q. Shi, Calculation of correlated ini- tial state in the hierarchical equations of motion method using an imaginary time path integral approach, J. Chem. Phys. 143, 194106 (2015)

    2015

  62. [70]

    T. Xing, T. Li, Y. Yan, S. Bai, and Q. Shi, Applica- tion of the imaginary time hierarchical equations of mo- tion method to calculate real time correlation functions, J. Chem. Phys. 156, 244102 (2022)

    2022

  63. [71]

    Z. H. Li, N. H. Tong, X. Zheng, D. Hou, J. H. Wei, J. Hu, and Y. J. Yan, Hierarchical Liouville-space approach for accurate and universal characterization of quantum im- purity systems, Phys. Rev. Lett. 109, 266403 (2012)

    2012

  64. [72]

    Bhattacharyya, T

    S. Bhattacharyya, T. Sayer, and A. Montoya-Castillo, Anomalous transport of small polarons arises from transient lattice relaxation or immovable boundaries, J. Phys. Chem. Lett. 15, 1382 (2024)

    2024

  65. [73]

    Song and Q

    L. Song and Q. Shi, Hierarchical equations of mo- tion method applied to nonequilibrium heat transport in model molecular junctions: Transient heat cur- rent and high-order moments of the current operator, Phys. Rev. B 95, 064308 (2017)

    2017

  66. [74]

    Schinabeck, R

    C. Schinabeck, R. H¨ artle, and M. Thoss, Hier- archical quantum master equation approach to electronic-vibrational coupling in nonequilibrium transport through nanosystems: Reservoir formu- lation and application to vibrational instabilities, Phys. Rev. B 97, 235429 (2018)

    2018

  67. [75]

    Seibt and T

    J. Seibt and T. Manˇ cal, Treatment of Herzberg- Teller and non-Condon effects in optical spec- tra with Hierarchical Equations of Motion, Chem. Phys. 515, 129 (2018)

    2018

  68. [76]

    L. Zhu, H. Liu, W. Xie, and Q. Shi, Explicit system-bath correlation calculated using the hierarchical equations o f motion method, J. Chem. Phys. 137, 194106 (2012)

    2012

  69. [77]

    Kato and Y

    A. Kato and Y. Tanimura, Quantum heat cur- rent under non-perturbative and non-Markovian conditions: Applications to heat machines, J. Chem. Phys. 145, 224105 (2016)

    2016

  70. [78]

    Yan, Theory of open quantum systems with bath of electrons and phonons and spins: Many-dissipaton density matrixes approach, J

    Y. Yan, Theory of open quantum systems with bath of electrons and phonons and spins: Many-dissipaton density matrixes approach, J. Chem. Phys. 140, 054105 (2014)

    2014

  71. [79]

    Y. Yan, J. Jin, R.-X. Xu, and X. Zheng, Dissipation equation of motion approach to open quantum systems, Front. Phys. 11, 110306 (2016)

    2016

  72. [80]

    Zhang, R.-X

    H.-D. Zhang, R.-X. Xu, X. Zheng, and Y. Yan, Statis- tical quasi-particle theory for open quantum systems, Mol. Phys. 116, 780 (2018)

    2018

  73. [81]

    Wang and Y

    Y. Wang and Y. Yan, Quantum mechan- ics of open systems: Dissipaton theories, J. Chem. Phys. 157, 170901 (2022)

    2022

  74. [82]

    Jankovi´ c, Charge transport limited by nonlo- cal electron–phonon interaction

    V. Jankovi´ c, Charge transport limited by nonlo- cal electron–phonon interaction. II. Numerically exact quantum dynamics in the slow-phonon regime (2025), arXiv:2501.05055 [cond-mat.str-el]

    2025 arXiv

  75. [83]

    Shun-jin and W

    W. Shun-jin and W. Cassing, Explicit treatment of n- body correlations within a density-matrix formalism, Ann. Phys. 159, 328 (1985)

    1985

  76. [84]

    Kira and S

    M. Kira and S. Koch, Many-body correlations and excitonic effects in semiconductor spectroscopy, Prog. Quantum Electron. 30, 155 (2006)

    2006

  77. [85]

    Kira and S

    M. Kira and S. W. Koch, Quantum-optical spectroscopy of semiconductors, Phys. Rev. A 73, 013813 (2006)

    2006

  78. [86]

    Kira and S

    M. Kira and S. W. Koch, Semiconductor Quantum Op- tics (Cambridge University Press, New York, 2012)

    2012

  79. [87]

    Schilp, T

    J. Schilp, T. Kuhn, and G. Mahler, Electron- phonon quantum kinetics in pulse-excited semi- conductors: Memory and renormalization effects, Phys. Rev. B 50, 5435 (1994)

    1994

  80. [88]

    V. M. Axt and S. Mukamel, Nonlinear optics of semicon- ductor and molecular nanostructures; a common per- spective, Rev. Mod. Phys. 70, 145 (1998)

    1998

  81. [89]

    Jankovi´ c and N

    V. Jankovi´ c and N. Vukmirovi´ c, Dynamics of exciton formation and relaxation in photoexcited semiconduc- tors, Phys. Rev. B 92, 235208 (2015)

    2015

  82. [90]

    See the Supplemental Material for (i) the details on inferring the generalized Wick’s theorem from the dy- namical equations of the HEOM method, (ii) the im- 28 plications of the time-reversal symmetry for the cross contribution to Cjj (t), (iii) the details of HEOM com- putat...

  83. [91]

    De Filippis, V

    G. De Filippis, V. Cataudella, A. S. Mishchenko, N. Nagaosa, A. Fierro, and A. de Candia, Crossover from super- to subdiffusive motion and mem- ory effects in crystalline organic semiconductors, Phys. Rev. Lett. 114, 086601 (2015)

    2015

  84. [92]

    Y. Zhao, D. W. Brown, and K. Lindenberg, On the Munn–Silbey approach to nonlocal exciton–phonon cou- pling, J. Chem. Phys. 100, 2335 (1994)

    1994

  85. [93]

    Grabert, P

    H. Grabert, P. Schramm, and G.-L. Ingold, Quantum brownian motion: The functional integral approach, Phys. Rep. 168, 115 (1988)

    1988

  86. [94]

    Q. Shi, L. Chen, G. Nan, R.-X. Xu, and Y. Yan, Efficient hierarchical Liouville space propagator to quantum dis- sipative dynamics, J. Chem. Phys. 130, 084105 (2009)

    2009

  87. [95]

    Kuhn, Density matrix theory of coherent ultrafast dynamics, in Theory of Transport Properties of Semi- conductor Nanostructures , edited by E

    T. Kuhn, Density matrix theory of coherent ultrafast dynamics, in Theory of Transport Properties of Semi- conductor Nanostructures , edited by E. Sch¨ oll (Chap- man and Hall, London, 1998)

    1998

  88. [96]

    Rossi and T

    F. Rossi and T. Kuhn, Theory of ultrafast phenomena in photoexcited semiconductors, Rev. Mod. Phys. 74, 895 (2002)

    2002

  89. [97]

    D. Hou, S. Wang, R. Wang, L. Z. Ye, R. X. Xu, X. Zheng, and Y. J. Yan, Improving the efficiency of hierarchical equations of motion approach and applica- tion to coherent dynamics in Aharonov–Bohm interfer- ometers, J. Chem. Phys. 142, 104112 (2015)

    2015

  90. [98]

    T. P. Fay, A simple improved low temperature cor- rection for the hierarchical equations of motion, J. Chem. Phys. 157, 054108 (2022)

    2022

  91. [99]

    Bonˇ ca, S

    J. Bonˇ ca, S. A. Trugman, and M. Berciu, Spectral function of the Holstein polaron at finite temperature, Phys. Rev. B 100, 094307 (2019)

    2019

  92. [100]

    Jansen, J

    D. Jansen, J. Bonˇ ca, and F. Heidrich-Meisner, Finite-temperature density-matrix renormalization group method for electron-phonon systems: Ther- modynamics and Holstein-polaron spectral functions, Phys. Rev. B 102, 165155 (2020)

    2020

  93. [101]

    Jansen, J

    D. Jansen, J. Bonˇ ca, and F. Heidrich-Meisner, Finite-temperature optical conductivity with density- matrix renormalization group methods for the Hol- stein polaron and bipolaron with dispersive phonons, Phys. Rev. B 106, 155129 (2022)

    2022

  94. [102]

    Miladi´ c and N

    S. Miladi´ c and N. Vukmirovi´ c, Method for ob- taining polaron mobility using real and imag- inary time path-integral quantum Monte Carlo, Phys. Rev. B 107, 184315 (2023)

    2023

  95. [103]

    Mitri´ c, Dynamical quantum typicality: A simple method for investigating transport prop- erties applied to the Holstein model (2024), arXiv:2412.17436 [cond-mat.str-el]

    P. Mitri´ c, Dynamical quantum typicality: A simple method for investigating transport prop- erties applied to the Holstein model (2024), arXiv:2412.17436 [cond-mat.str-el]

    2024 arXiv

  96. [105]

    Sumi, Theory of electrical conduction in organic molecular crystals: Temperature-independent mobili- ties, J

    H. Sumi, Theory of electrical conduction in organic molecular crystals: Temperature-independent mobili- ties, J. Chem. Phys. 70, 3775 (1979)

    1979

  97. [106]

    Allen, Chapter 6 electron transport, in Conceptual Foundations of Materials , Contempo- rary Concepts of Condensed Matter Science, Vol

    P. Allen, Chapter 6 electron transport, in Conceptual Foundations of Materials , Contempo- rary Concepts of Condensed Matter Science, Vol. 2, edited by S. G. Louie and M. L. Cohen (Elsevier, 2006) pp. 165–218

    2006

  98. [107]

    Jacoboni, Theory of Electron Transport in Semicon- ductors (Springer, Berlin, 2010)

    C. Jacoboni, Theory of Electron Transport in Semicon- ductors (Springer, Berlin, 2010)

    2010

  99. [108]

    Ishizaki and G

    A. Ishizaki and G. R. Fleming, Unified treatment of quantum coherent and incoherent hopping dynamics in electronic energy transfer: Reduced hierarchy equation approach, J. Chem. Phys. 130, 234111 (2009)

    2009

  100. [109]

    Jankovi´ c and T

    V. Jankovi´ c and T. Manˇ cal, Exact description of exci- tonic dynamics in molecular aggregates weakly driven by light, J. Chem. Phys. 153, 244122 (2020)

    2020

  101. [110]

    Tanimura, Reduced hierarchical equations of motion in real and imaginary time: Corre- lated initial states and thermodynamic quantities, J

    Y. Tanimura, Reduced hierarchical equations of motion in real and imaginary time: Corre- lated initial states and thermodynamic quantities, J. Chem. Phys. 141, 044114 (2014)

    2014

  102. [111]

    A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (McGraw-Hill Book Company, 1971)

    1971

  103. [112]

    Y. Su, Y. Wang, and W. Dou, Nonperturbative open quantum dynamics bypass influence functional (2025), arXiv:2503.00297 [quant-ph]

    2025 arXiv

  104. [113]

    Lee, J.-J

    N.-E. Lee, J.-J. Zhou, L. A. Agapito, and M. Bernardi, Charge transport in organic molecular semiconductors from first principles: The bandlike hole mobility in a naphthalene crystal, Phys. Rev. B 97, 115203 (2018)

    2018

  105. [114]

    Vukmirovi´ c, Calculations of electron mobility in II- VI semiconductors, Phys

    N. Vukmirovi´ c, Calculations of electron mobility in II- VI semiconductors, Phys. Rev. B 104, 085203 (2021) . Supplemental Material for Charge transport limited by nonlocal electron–phonon inte raction. I. Hierarchical equations of motion approach Veljko Jankovi´ c∗ Institute...

    2021

  106. [115]

    Jankovi´ c, Numerical investigation of transport pro perties of the one-dimensional Peierls model based on the hi erarchical equations of motion, 10.5281/zenodo.14637019 (2025)

    V. Jankovi´ c, Numerical investigation of transport pro perties of the one-dimensional Peierls model based on the hi erarchical equations of motion, 10.5281/zenodo.14637019 (2025)

    2025 doi

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.