{"id":"b3a77d29-5e75-460d-aebb-0e395413ab45","arxiv_id":"2501.05054","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"HEOM auxiliaries are expressed through phonon operators, the generalized Wick theorem is proved, and exact mobility is computed for the 1D Peierls model.","lead":"This paper devises a way to compute how an electron moves through a vibrating crystal when the vibrations change the electron's hopping strength. The method yields numerically exact mobilities for a standard model and shows that at high temperature the vibrations themselves transport the charge.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hierarchy closing (Sec. IV B) is an uncontrolled approximation; the 'numerically exact' mobility claim rests on MA/DR agreement that cannot exclude a common RPA bias.","rationale":"The reader's weakest_assumption is exactly the hierarchy closing scheme of Sec. IV B and Appendix D. My stress-test analysis independently reaches the same conclusion: the closing scheme is the main uncontrolled approximation in the numerical pipeline, and the paper's own validation (two closing schemes, D-convergence at high temperature, and the sum-rule accuracy) does not provide a rigorous error bound. The concern is load-bearing because every reported mobility value and trend depends on the accuracy of Eq. (46) with the MA or DR form of Γ(k,n). A concrete exact-diagonalization cross-check at a representative parameter point would settle whether the closing bias is within the claimed uncertainty. I do not identify a separate fatal flaw in the formal derivation of Eq. (10) or the generalized Wick theorem; the proof in Appendix C appears internally consistent, though the derivation of Eq. (10) in Appendix B contains a step (the assertion that each term in Eq. (B12) vanishes separately) that is not fully justified, but the exact-diagonalization test would also implicitly verify the overall framework. Therefore the conditional verdict is appropriate; no change is needed.","tokens_in":40812,"tokens_out":8483,"duration_ms":90405,"concrete_test":"Perform exact diagonalization of the full electron-phonon Hamiltonian for N=7, ω0=J=1, λ=0.5, T=5, with a phonon truncation of, say, up to 10 quanta per mode, and compute μ_dc via the Kubo formula. Compare with the HEOM result at the same parameters. If the difference exceeds the claimed ~10% relative uncertainty, the closing scheme introduces uncontrolled error and the 'numerically exact' label must be revised; if the difference is below this threshold, the closing is validated in a representative case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim is that HEOM yields the numerically exact dynamical mobility of the 1D Peierls model. The paper explicitly identifies the hierarchy closing scheme as the main approximation (Sec. IV B): the MA closing (Eq. 47) is derived by neglecting depth-D+2 couplings, applying Markovian and adiabatic approximations, and invoking the random-phase approximation to discard off-diagonal (in momentum) terms in Eq. (D6). The DR closing (Eq. D10) shares the same RPA. Neither closing has a controlled error estimate, and the validation in Sec. V D compares only these two schemes; agreement between them does not rule out a common bias introduced by RPA. Moreover, at moderate temperatures and interactions the method does not converge with D or N: Fig. 2 shows D(t) failing to saturate for λ=0.25, T=1, requiring a fitting procedure and yielding a stated ~10% uncertainty. The T=1 entries in Fig. 4(a) are explicitly flagged as possibly not fully converged. Since the abstract and Sec. VI describe the results as numerically exact, the uncontrolled closing error is load-bearing: if the closing is inaccurate, the reported mobility values and the temperature/interaction trends could be systematically wrong even in the claimed regime of applicability.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":41051,"tokens_out":3948,"duration_ms":38596,"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":[{"comment":"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.","section":"Sec. IV B and Eq. (46)"},{"comment":"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.","section":"Sec. V C and Fig. 2"},{"comment":"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.","section":"Sec. V D and Eq. (52)"}],"minor_comments":[{"comment":"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.","section":"Abstract and Sec. VI"},{"comment":"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.","section":"Sec. I"},{"comment":"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.","section":"Fig. 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the journal's readership given the importance of nonlocal electron–phonon interactions in organic semiconductors. The main concern is the overstatement of 'numerically exact' in light of the uncontrolled hierarchy closing. I believe the formal results and the method are significant, but the numerical claims need to be made more cautious or supported by additional validation. The self-citations are appropriate in the context of prior HEOM implementations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news is the formalism, not the mobility curves. Explicitly expressing HEOM auxiliaries as products of phonon operators with lower-order correlations subtracted, then proving the generalized Wick theorem for undamped discrete modes, is a solid and useful step. It clears up something DEOM treatments usually leave opaque. The derivations in Sec. II and the appendices are detailed, self-contained, and look internally consistent. The decomposition of the current–current correlation function into band, phonon-assisted, and cross terms is clean, and the companion-paper plan makes sense.\n\nThe soft spot is the closing, and the paper is more honest in the body than in the abstract. The abstract says 'numerically exact,' but the body concedes the hierarchy closing is the main approximation. The MA and DR closures both invoke an RPA, and neither has a controlled error bound. Agreement between them is weak evidence: a common RPA bias would survive. At lambda=0.25, T=1, D(t) does not saturate; the mobility is obtained by fitting and comes with a stated ~10% uncertainty. The T=1 points in Fig. 4(a) are flagged as possibly not fully converged. The fast-phonon results at omega0=3 are explicitly not entirely reliable. So the 'numerically exact' label overstates the numerical part. The formal part is not affected by the closing, and that is the part worth publishing.\n\nMinor points: no per-point error bars, and code is not provided (data are). The self-citations are fine—they cite prior HEOM implementations and do not force the conclusion.\n\nWho is this for? People building HEOM/DEOM methods and those comparing approximate transport schemes to a benchmark. The paper deserves a serious referee, not a desk rejection. I would send it out and tell the referee to focus on Sec. IVB/D and to decide whether the abstract's 'numerically exact' should be tempered. With the wording fixed and a bit more transparency on closing bias, the paper is publishable.","headline":"The formalism is the real contribution; the mobility numbers are honest but not literally 'numerically exact' because the hierarchy closing is uncontrolled.","tokens_in":41554,"tokens_out":2937,"would_cite":true,"duration_ms":29097,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"HEOM auxiliaries encode phonon-assisted transport in the 1D Peierls model, giving numerically exact mobility.","keywords":["hierarchical equations of motion","nonlocal electron-phonon interaction","Peierls model","carrier mobility","phonon-assisted current","generalized Wick's theorem","organic semiconductors","quantum transport"],"falsifier":"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.","tokens_in":40599,"feed_emoji":"⚡","tokens_out":8031,"duration_ms":75997,"temperature":0.7,"pith_summary":"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.","feed_headline":"Exact phonon-assisted mobility computed in 1D Peierls model","feed_subtitle":"HEOM recovers the full current correlation, so mobility flips from decreasing to increasing with interaction as temperature rises.","key_machinery":"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).","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the hierarchical equations of motion formalism that the paper extends to mixed electron–phonon correlation functions.","marker":"[46]"},{"why":"Provides the HEOM formulation for bosonic canonical baths that is the formal starting point for the auxiliary hierarchy.","marker":"[47]"},{"why":"The dissipaton theory in which the generalized Wick's theorem originally appears; the paper proves it here rather than postulating it.","marker":"[72]"},{"why":"Systematizes the dissipaton algebra and the generalized Wick's theorem notation that Eqs. (14) and (15) of the paper adopt.","marker":"[74]"},{"why":"A prior HEOM treatment of the Holstein model contributes the momentum-space hierarchy and the Markovian–adiabatic closing scheme adapted in this paper.","marker":"[60]"},{"why":"Introduces the phonon-assisted current in linear-response theory and provides the early high-temperature formula that the paper compares against.","marker":"[58]"},{"why":"Defines the Su–Schrieffer–Heeger-type model, its mobility, and the optical sum rule used here as a convergence check.","marker":"[14]"},{"why":"Defines the organic-semiconductor parameter regime of interest and the approximate high-temperature mobility formula, Eq. (59), used for comparison.","marker":"[22]"}],"fun_headline_variants":["Mobility trend flips as temperature rises in 1D Peierls model","Exact mobility: phonon-assisted current overtakes band at high T","HEOM yields exact mobility for nonlocal phonon coupling in 1D","Phonon-assisted current flips mobility trend at high temperatures","Exact 1D mobility from HEOM: nonlocal phonons solved"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Mobility trend flips as temperature rises in 1D Peierls model","Exact mobility: phonon-assisted current overtakes band at high T","HEOM yields exact mobility for nonlocal phonon coupling in 1D","Phonon-assisted current flips mobility trend at high temperatures","Exact 1D mobility from HEOM: nonlocal phonons solved"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001346,"raw_usage":{"total_tokens":5510,"prompt_tokens":1030,"completion_tokens":4480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":4382}},"tokens_in":646,"tokens_out":4480,"duration_ms":30431,"temperature":1.0,"reasoning_tokens":4382,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:19:22.879994+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Landi, Charge mobility prediction in organic semiconductors: Comparison of second-order cumu- lant approximation and transient localization theory, J","cited_arxiv_id":null,"evidence_quote":"Provides the HEOM formulation for bosonic canonical baths that is the formal starting point for the auxiliary hierarchy."},{"cited_title":"Landi, R","cited_arxiv_id":null,"evidence_quote":"Supplies the hierarchical equations of motion formalism that the paper extends to mixed electron–phonon correlation functions."},{"cited_title":"Bhattacharyya, T","cited_arxiv_id":null,"evidence_quote":"The dissipaton theory in which the generalized Wick's theorem originally appears; the paper proves it here rather than postulating it."},{"cited_title":"Schinabeck, R","cited_arxiv_id":null,"evidence_quote":"Systematizes the dissipaton algebra and the generalized Wick's theorem notation that Eqs. (14) and (15) of the paper adopt."},{"cited_title":"Kramer, M","cited_arxiv_id":null,"evidence_quote":"A prior HEOM treatment of the Holstein model contributes the momentum-space hierarchy and the Markovian–adiabatic closing scheme adapted in this paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the phonon-assisted current in linear-response theory and provides the early high-temperature formula that the paper compares against."}],"review_version":1}