Pith. sign in

REVIEW 2 major objections 4 minor 41 references

Coherent migration of the single excitation injected into finite-length segment of the biomolecular chain

T0 review · 2 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash

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

Pith's one-line read A single excitation injected off-center into a finite molecular segment propagates with a left-right asymmetric probability distribution, proving that initial-state geometry alone can break transport symmetry.

desk verdict Solid math, oversold novelty; the coherent regime may not hold for the paper's own parameters. read the letter →

arxiv 2608.02131 v1 pith:QSYXYI2T submitted 2026-08-03 cond-mat.other

D. Chevizovich , V. Matic , Z. Przulj , S. Galovic This is my paper · ORCID
classification cond-mat.other
keywords singleexcitationmigrationnonadiabaticpolaronfinitemolecularsegmentasymmetricprobabilitydistributioncontinuous-timequantumwalkChebyshevpolynomialsHolsteinmodelcoherentdephasing
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 studies the coherent migration of a single vibron excitation injected into a finite segment of a biomolecular chain, where the excitation couples to thermal vibrations and forms a self-trapped nonadiabatic polaron. The authors solve the mean-field Holstein model exactly, obtaining a closed-form expression for the time-dependent probability of finding the excitation at every site of the segment. Their central claim is that the probability distribution is generically asymmetric with respect to the initially excited site, even though the chain itself is homogeneous, and that this asymmetry comes solely from the off-center placement of the initial excitation. Only when the excitation starts at the exact center does the distribution become symmetric. They also show that non-equidistant mode frequencies prevent full revivals and cause the initially localized probability peak to fragment into a dominant maximum with secondary maxima.

What carries the argument

The key object is the exact analytical solution for the correlation function V_n(τ), built from modified Chebyshev polynomials of the second kind. After a Lang-Firsov transformation, a mean-field average over renormalized phonons, and a Laplace transform in time, each correlation function becomes a coherent superposition of K discrete modes with frequencies Ω_k = B e^{-S coth(1/2θ)} cos(kπ/(K+1)), weighted by Chebyshev products Ψ_{n,0}^{{M,N}}(φ_k). These weights carry the entire dependence on the initial position: M and N count the sites to the left and right of the injection point, and the asymmetric product sin[(M+n+1)φ_k] sin[(M+1)φ_k] on the right side versus its left-side counterpart i

What would settle it

Numerically diagonalize the full Holstein Hamiltonian including H_rest for a 21-site segment (S=0.3, B=0.1, θ=4) with an off-center initial site, and compare the resulting probability distribution p_n(t) with Eqs. (14)-(15). If the left-right imbalance deviates appreciably at times shorter than the dephasing time, or if a centered initial site develops asymmetry, the mean-field neglect is falsified. Alternatively, in an experimental realization of a finite tight-binding chain with phonon-like coupling, measure p_{n0+d}(t)/p_{n0-d}(t): it should differ from 1 for off-center injection and equal

Watch

Extended reading notes

Core claim

The central discovery is a symmetry-breaking effect in a closed, undriven, homogeneous quantum chain. For a finite segment with reflecting boundaries and a single initial excitation, the probability of finding the excitation at nodes on opposite sides of the injection point is generally different. The asymmetry is not introduced by the Hamiltonian—the transfer integrals are identical—but by the initial condition: the injection site's position relative to the two boundaries. In the exact analytical solution, this enters through Chebyshev weighting factors Ψ_{n,0}^{{M,N}} that differ between the left and right branches of the solution. The exception is the centered initial site, where left and

Load-bearing premise

The residual interaction between the dressed polaron and the renormalized phonons (H_rest) is neglected, so the excitation evolves coherently under the mean-field Hamiltonian; if that interaction is not weak, the predicted asymmetry and dephasing would be washed out by genuine decoherence.

Editorial extensions

If this is right

  • Off-center initial excitation in any finite, homogeneous, reflecting-boundary tight-binding chain produces a left-right asymmetric detection probability; only center injection stays symmetric.
  • Temperature enters the dynamics by renormalizing the hopping amplitude via the factor e^{-S coth(1/2θ)}, rather than by introducing decoherence, so the coherent quantum-walk description holds with a temperature-tuned hopping.
  • The non-equidistant mode spectrum means the system never fully revives; only partial rephasing occurs, so a localized excitation fragments into one dominant peak plus several smaller ones.
  • The residence time of the excitation near a given site is approximately T_min ∼ 2π/(B e^{S coth(1/2θ)}), which grows exponentially with coupling S and inverse temperature 1/θ.
  • The formal mapping to continuous-time quantum walks provides a direct route to use quantum-walk theory for polaron transport in biomolecular chains, with temperature-dependent hopping bridging idealized models and molecular environments.

Reading between the lines

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

  • Because the asymmetry arises purely from the initial condition relative to boundaries, the prediction should also hold for any coherent tight-binding simulator—for example, arrays of coupled waveguides or cold atoms in optical lattices—provided the effective Hamiltonian has reflecting boundaries and negligible dissipation; measuring the imbalance would directly test the mean-field approximation.
  • The model's neglect of H_rest means the long-time dephasing is an interference effect, not true environmental decoherence; in a real molecular chain, polaron-phonon scattering would likely wash out the fine interference pattern after a dissipative timescale, so the predicted fragmented maxima may only be observable at short to intermediate times.
  • The same Chebyshev decomposition could be extended to absorbing or semi-infinite boundary conditions to compute first-passage and transfer probabilities to an active site, turning the asymmetry into a design principle for directing excitation flow in finite molecular wires.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The paper studies single-excitation transport in a finite segment of a molecular chain modeled by the Holstein Hamiltonian with optical phonons. After a Lang-Firsov transformation and thermal mean-field averaging, the dynamics reduces to a K-site tight-binding chain with renormalized hopping J_eff = (B/2)e^{-S coth(1/2θ)} and reflecting boundaries (Eqs. (1)-(5)). The authors solve the polaron correlation function V_n(τ) via Laplace transform and Chebyshev polynomials, obtaining the closed-form superposition Eq. (14), and compute p_n(τ)=|V_n(τ)|^2. They report that for an off-center initial site, symmetrically placed sites have different temporal probability spectra, with exact symmetry only for a central initial site; the non-equidistant spectrum Ω_k causes dephasing and fragmentation of the probability peak. The model is connected to continuous-time quantum walks. Numerical results use S=0.3, B=0.1, θ=4.

Significance. The central formal result is a strength: Eq. (14) is an explicit, easily checkable exact solution of the effective tight-binding model, and the claimed symmetry properties follow directly from the weighting coefficients in Eq. (13) rather than from an assumed ansatz. The connection to CTQW is appropriate and gives the result a useful conceptual context. The paper is also explicit that dissipative processes are not captured. However, the physical claim that this describes coherent migration in real biomolecular segments at the stated biological parameters is not yet supported because the validity of the neglect of H_rest is not quantified. The paper's main value is therefore as an exact treatment of a clearly defined mean-field model; its applicability to the motivating biomolecular scenario remains conditional.

major comments (2)
  1. [Sec. II, Eqs. (5)-(6); Sec. III parameters] The central approximation is that the residual polaron-phonon coupling H_rest is weak and can be neglected, but no quantitative validity condition is given. For the parameter set used in all figures (S=0.3, B=0.1, θ=4), S/B=3 is only marginally above the stated S/B≫1 criterion. Moreover, W=S coth(1/(2θ))≈2.4, so e^{-W}≈0.09 and J_eff=(B/2)e^{-W}ℏω0≈0.0045ℏω0, while the residual couplings in H_rest scale as the bare J0=0.05ℏω0 times O(1) fluctuations of T± at the thermal occupation n_q≈3.5. A rough estimate gives H_rest/J_eff=O(e^W)≈10, suggesting that the polaron damping rate may be comparable to or larger than the mode spacings Ω_k−Ω_l that produce the interference asymmetry. The manuscript needs a quantitative estimate of the incoherent rate, or an explicit restricted regime, before the physical asymmetry claim for biomolecular segments can be accepted.
  2. [Sec. II, Eqs. (7)-(9)] The initial state and the observable are defined in the polaron representation: after the Lang-Firsov transformation, a†_n creates a dressed polaron, V_n(τ) is the polaron propagator, and p_n(τ) is the polaron occupation probability. The text, however, repeatedly refers to an 'injected excitation' and to 'the probability of finding the excitation'. A bare exciton injected at n0 would have a phonon coherent component in the polaron representation, so Eq. (9) would not describe it, and the bare-exciton occupation would differ from p_n(τ). Please state explicitly that the initial state is assumed to be a pre-formed polaron and justify the polaron-formation timescale, or compute the appropriate observable for the injection scenario.
minor comments (4)
  1. [Sec. II, Eqs. (11)-(12)] The definition D_n(x)=U_n(2x) is inconsistent with the trigonometric representation D_n(x)=sin((n+1)φ)/sinφ, x=2cosφ, and with the zero locations x_k=2cos(kπ/(n+1)). The correct relation is D_n(x)=U_n(x/2). Please correct the definition or the substitution.
  2. [Sec. III.B] The printed formula T_min∼2π/[B e^{S coth(1/2θ)}] is inconsistent with Ω_k=B e^{-S coth(1/2θ)}cosφ. The minimal period should be (2π/B)e^{S coth(1/2θ)}≈700, matching the numerical value given in the text; as printed the formula gives ≈5.6.
  3. [Sec. II, Eqs. (10)-(13)] The derivation from Eq. (10) to Eq. (13) is sketched only by reference to Ref. 9. An appendix containing the partial-fraction identity and the relevant Chebyshev identities would make the paper self-contained.
  4. [Fig. 7 caption] The caption contains typos ('ant the fifith SE') and 'Rectangled area' should be 'Rectangular area'. The definition of 'above/below' relative to the panel layout should also be clarified.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the asymmetry and dephasing claims are derived consequences of the explicitly stated equations of motion; the lone self-citation (Ref. 9) is a solution-method citation that is independently checkable from Eqs. (8)-(14), and the paper's central limitation (neglect of H_rest) is stated openly rather than disguised as a prediction.

full rationale

The paper's derivation chain is self-contained and non-circular. Starting from the Holstein-type Hamiltonian (1)-(4), the Lang-Firsov transformation and a mean-field average produce the effective Hamiltonian (5), whose only temperature dependence is the narrowing factor e^{-W(T)} derived from Eq. (6) (for Einstein phonons, W = S coth(1/2θ)). The equations of motion (8) for the correlation functions are written explicitly, the Laplace substitution is given, and the solution (10), (14) is displayed in full; the asymmetry claim in Sec. III (observation B) is read directly from the differing Chebyshev-numerator products D_{N-n}D_M versus D_{M-n}D_N in Eq. (10), and the center-of-chain symmetry is the M = N case. Thus the central results (asymmetry, non-equidistant spectrum dephasing, fragmented probability) are derived, not assumed, and no parameter is fitted to any predicted output (S = 0.3, B = 0.1, θ = 4 are literature values for amide-I systems). The self-citation 'following the procedure of 9' (Sec. II) is not load-bearing: the same linear system and solution are stated in the present paper and can be verified independently (Chebyshev finite-chain Green's functions are standard); the citation is not used to forbid alternatives or to import an unverifiable premise. The CTQW analogy is explicitly labeled 'formal analogy', and the temperature-independence of coherence is an honest consequence of Hamiltonian (5), with the paper itself adding the qualifier that this 'is valid only for systems where the residual interaction between the polaron and the renormalized phonons is weak' (Conclusion). The real qualification is the explicit mean-field neglect of H_rest: Sec. II states 'dissipative processes associated with incoherent polaron-phonon scattering are not captured by the present model.' Whether S/B = 3 and the θ = 4 narrowing (e^{-0.3·coth(1/8)} ≈ 0.09) place the system in that regime is an unquantified validity concern (correctness risk), not a circular step; the predictions follow from the stated model and do not presuppose the conclusions. Score 2 reflects the single non-load-bearing self-citation; otherwise no circularity.

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

The paper introduces no new entities; the nonadiabatic polaron is a known quasiparticle. The free parameters are physical inputs chosen from literature ranges, not fitted to the central claim. The key assumptions are the reflecting boundaries and the neglect of residual polaron-phonon scattering, both standard in this modeling tradition.

free parameters (4)
  • S (polaron-phonon coupling) = 0.3
    Dimensionless coupling, chosen from the literature range for amide-I vibrons (chi ~ 32-62 pN). Not fitted to the central result, but sets the magnitude of the narrowing factor.
  • B (adiabatic parameter) = 0.1
    Chosen from J0 ~ 0.5-0.97 meV. Sets the hopping amplitude and the time scale of the dynamics.
  • theta (normalized temperature) = 4
    Room-temperature value, chosen as biologically relevant. Controls the exponential narrowing of the transfer integral.
  • M, N (segment geometry) = M=4, N=16, K=21
    Illustrative segment geometry used in the figures. The asymmetry claim holds for any M != N.
assumptions (4)
  • domain assumption The excitation is confined to the finite segment: J0(n)=0 at the segment edges (reflecting boundaries).
    Stated in Section II. Required for the central asymmetry claim, which depends on finite boundaries.
  • domain assumption The mean-field approximation is valid and the residual polaron-phonon interaction H_rest is negligible.
    Stated in Section II. If this fails, the coherent tight-binding picture and the calculated probabilities break down.
  • standard math The Lang-Firsov transformation and thermal averaging yield the standard narrowing factor e^{-W} with W = S coth(1/(2θ)).
    Standard polaron theory; used in Eq. (5) and Eq. (14b).
  • domain assumption Only optical phonons at frequency ω0 couple to the excitation.
    Justified by neutron and pump-probe experiments on ACN (Refs 25-26); restricts the model to optical phonons.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Coherent migration of the single excitation injected into finite-length segment of the biomolecular chain." pith.science (2026). https://pith.science/paper/QSYXYI2T

@misc{pith2026260802131,
  author       = {Pith},
  title        = {Pith review of: Coherent migration of the single excitation injected into finite-length segment of the biomolecular chain},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QSYXYI2T}},
  note         = {Machine review of arXiv:2608.02131}
}
read the original abstract

We study the migration of a single excitation excited at a structural element of a finite molecular segment, which is a part of a long biomolecular chain. The excitation cannot leave the segment and is locally coupled to thermal vibrations of the lattice, forming a self-trapped state corresponding to a nonadiabatic polaron. The time-dependent probability distribution of finding the excitation at the nodes of the segment is calculated, with particular emphasis on the role of the initial excitation position. A formal analogy is observed between the present model and continuous-time quantum walk models on finite chains with reflecting boundaries. The results reveal an asymmetry in the probability distribution for nodes symmetrically positioned with respect to the initially excited site, which arises solely from the asymmetric placement of the initial excitation within the finite segment. The only exception occurs when the initially excited node is located at the center of the segment, where the probability distribution becomes symmetric. The complex interference pattern and the absence of well-defined revivals stem from the non-equidistant spectrum of mode frequencies, leading to progressive dephasing of the constituent modes. As a result, the initially well-localized probability maximum fragments into one dominant maximum accompanied by several secondary maxima of lower intensity. These findings highlight the importance of boundary conditions and initial-state geometry in controlling quantum transport in finite molecular systems.

Figures

Figures reproduced from arXiv: 2608.02131 by the authors.

Figure 1
Figure 1. FIG. 1. Schematic representation of the considered molecular [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. A schematic representation of the renumbered MS. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. TPS corresponding to the two nearest-neighbor sites [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: FIG. 5. TPS on the edge nodes of the MS, for the configura [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Excitation probability distribution, in dependence on the position of the initially excited node. Parameters: [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Time distribution of the excitation probability over [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references

  1. [1]

    Voet , title =

    Donald Voet and Judith G. Voet , title =

  2. [2]

    Nelson and Michael M

    David L. Nelson and Michael M. Cox , title =

  3. [3]

    Thierry Dauxois and Michel Peyrard , title =

  4. [4]

    Coherent Excitations in Biological Systems , publisher =

  5. [5]

    Lambert and Y

    N. Lambert and Y. N. Chen and Y. C. Cheng and C. M. Li and G. Y. Chen and F. Nori , title =. Nature Physics , volume =. 2013 , doi =

    2013

  6. [6]

    Angewandte Chemie International Edition , volume =

    Xi Chen and Xue Zhang and Xiao Xiao and Zhijia Wang and Jianzhang Zhao , title =. Angewandte Chemie International Edition , volume =. 2023 , doi =

    2023

  7. [7]

    Cruzeiro , title =

    L. Cruzeiro , title =. Low Temperature Physics , volume =. 2022 , doi =

    2022

  8. [8]

    P. H. Alvarez and L. Gerhards and I. A. Solov'yov and M. C. de Oliveira , title =. Frontiers in Quantum Science and Technology , volume =. 2024 , doi =

    2024

Show all 41 references
  1. [9]

    Chevizovich and S

    D. Chevizovich and S. Galovic and V. Matic and Z. Ivic and Z. Przulj , title =. AIP Advances , volume =. 2026 , doi =

    2026

  2. [10]

    Chevizovich , editor =

    D. Chevizovich , editor =. Vibron Self-trapping in Quasi-One-Dimensional Biomolecules: Non-adiabatic Polaron Approach , booktitle =. 2022 , doi =

    2022

  3. [11]

    Holstein , title =

    T. Holstein , title =. Annals of Physics , volume =. 1959 , doi =

    1959

  4. [12]

    I. G. Lang and Yu. A. Firsov , title =. Soviet Physics JETP , volume =

  5. [13]

    D. M. Alexander and J. A. Krumhansl , title =. Physical Review B , volume =. 1986 , doi =

    1986

  6. [14]

    Pouthier , title =

    V. Pouthier , title =. Journal of Chemical Physics , volume =. 2010 , doi =

    2010

  7. [15]

    Edler and R

    J. Edler and R. Pfister and V. Pouthier and C. Falvo and P. Hamm , title =. Physical Review Letters , volume =. 2004 , doi =

    2004

  8. [16]

    Careri and U

    G. Careri and U. Buontempo and F. Carta and E. Gratton and A. C. Scott , title =. Physical Review Letters , volume =. 1983 , doi =

    1983

  9. [17]

    G. B. Blanchet and C. R. Fincher Jr. , title =. Physical Review Letters , volume =. 1985 , doi =

    1985

  10. [18]

    Hamm and G

    P. Hamm and G. P. Tsironis , title =. European Physical Journal Special Topics , volume =. 2007 , doi =

    2007

  11. [19]

    Hamm and G

    P. Hamm and G. P. Tsironis , title =. Physical Review B , volume =. 2008 , doi =

    2008

  12. [20]

    Pouthier , title =

    V. Pouthier , title =. Physical Review E , volume =. 2008 , doi =

    2008

  13. [21]

    A. H. Castro Neto and A. O. Caldeira , title =. Physical Review B , volume =. 1992 , doi =

    1992

  14. [22]

    J. T. Devreese and A. S. Alexandrov , title =. Reports on Progress in Physics , volume =. 2009 , doi =

    2009

  15. [23]

    Kalosakas and S

    G. Kalosakas and S. Aubry and G. P. Tsironis , title =. Physical Review B , volume =. 1998 , doi =

    1998

  16. [24]

    Barthes , title =

    M. Barthes , title =. Journal of Molecular Liquids , volume =. 1989 , doi =

    1989

  17. [25]

    Edler and P

    J. Edler and P. Hamm and A. C. Scott , title =. Physical Review Letters , volume =. 2002 , doi =

    2002

  18. [26]

    Yarkony and R

    D. Yarkony and R. Silbey , title =. Journal of Chemical Physics , volume =. 1976 , doi =

    1976

  19. [27]

    Hennig , title =

    D. Hennig , title =. Physical Review B , volume =. 2002 , doi =

    2002

  20. [28]

    E. I. Rashba , editor =. Excitons , publisher =

  21. [29]

    G. D. Mahan , title =

  22. [30]

    N. A. Nevskaya and Yu. N. Chirgadze , title =. Biopolymers , volume =. 1976 , doi =

    1976

  23. [31]

    Falvo and V

    C. Falvo and V. Pouthier , title =. Journal of Chemical Physics , volume =. 2005 , doi =

    2005

  24. [32]

    Hamm and J

    P. Hamm and J. Edler , title =. Physical Review B , volume =. 2006 , doi =

    2006

  25. [33]

    Kalosakas , title =

    G. Kalosakas , title =. Physical Review E , volume =. 2011 , doi =

    2011

  26. [34]

    Kempe , title =

    J. Kempe , title =. Contemporary Physics , volume =. 2003 , doi =

    2003

  27. [35]

    O. M. Continuous-time quantum walks: Models for coherent transport on complex networks , journal =. 2011 , doi =

    2011

  28. [36]

    Fuss and L

    I. Fuss and L. White and P. Sherman and S. Naguleswaran , title =. 2007 , eprint =

    2007

  29. [37]

    Apers and L

    S. Apers and L. Miclo , title =. 2024 , eprint =

    2024

  30. [38]

    Reimann , title =

    P. Reimann , title =. Physics Reports , volume =. 2002 , doi =

    2002

  31. [39]

    S. V. Kozyrev and A. N. Pechen , title =. Physica Scripta , volume =. 2023 , doi =

    2023

  32. [40]

    Lau and O

    B. Lau and O. Kedem and J. Schwabacher and D. Kwasnieski and E. A. Weiss , title =. Materials Horizons , volume =. 2017 , doi =

    2017

  33. [41]

    The Journal of Chemical Physics , volume =

    Mohseni, Masoud and Rebentrost, Patrick and Lloyd, Seth and Aspuru-Guzik, Alán , title =. The Journal of Chemical Physics , volume =. 2008 , month =

    2008

Pith tools

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