Pith. sign in

REVIEW 3 major objections 6 minor 66 references

Fast, accurate simulation of polaron dynamics and multidimensional spectroscopy by multiple Davydov trial states

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

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

Pith's one-line read A superposition of multiple Davydov trial states, evolved under the Dirac-Frenkel variational principle, reproduces exact Holstein polaron dynamics and yields the first computed 2D spectra for off-diagonal exciton-phonon coupling.

desk verdict A solid diagonal-coupling validation of multi-Davydov dynamics, but the first-ever 2D spectra claim for off-diagonal coupling rests on an unbenchmarked propagator and an under-specified Appendix D. read the letter →

arxiv 1908.09243 v1 pith:NM7O6MED submitted 2019-08-25 cond-mat.stat-mech cond-mat.str-elphysics.chem-phquant-ph

Nengji Zhou , Lipeng Chen , Zhongkai Huang , Kewei Sun , Yoshitaka Tanimura , Yang Zhao This is my paper · ORCID
classification cond-mat.stat-mechcond-mat.str-elphysics.chem-phquant-ph
keywords HolsteinpolaronDavydovtrialstatesDirac-Frenkelvariationalprincipleoff-diagonalexciton-phononcouplinghierarchyequationsofmotiontwo-dimensionalelectronicspectroscopylinearabsorptionspectrumexcitonself-trapping
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

The paper aims to show that a linear combination of Davydov trial states, called the multi-D1 and multi-D2 Ansätze, gives fast and accurate variational dynamics for the Holstein molecular crystal when exciton-phonon coupling is both diagonal and off-diagonal. It argues that with large enough multiplicity M the variational results become numerically exact, supporting this with near-overlap of multi-D1 exciton probabilities against hierarchy-equations-of-motion results for a diagonal-coupling case. It then uses the multi-D2 Ansatz to compute, for the first time, two-dimensional spectra of a system with off-diagonal coupling. The reported spectra show a single peak for weak off-diagonal coupling and a vibronic multi-peak structure for strong coupling. A sympathetic reader cares because this offers a practical variational route to multidimensional spectroscopy of molecular aggregates without exact wavefunction propagation.

What carries the argument

The load-bearing object is the multiple Davydov trial state: |D^M_1⟩ = Σ_{i=1}^M Σ_{n=1}^N ψ_{i,n}|n⟩|λ_{i,n}⟩ for multi-D1 (phonon displacements depend on the exciton site) and |D^M_2⟩ = Σ_{i=1}^M Σ_{n=1}^N ψ_{i,n}|n⟩|λ_i⟩ for multi-D2 (site-independent displacements). Time evolution of the variational parameters ψ and λ is generated by the Dirac-Frenkel time-dependent variational principle from the Lagrangian L = ⟨D|(iℏ/2)∂↔/∂t − H|D⟩. The multiplicity M controls accuracy: M=1 restores the standard Davydov D1 or D2 Ansatz, and increasing M systematically reduces the deviation vector δ(t) = $iℏ^{{-1}}$H|D⟩ − ∂_t|D⟩. For the 2D spectra, system propagators are approximated by the multi-D2 Ansatz and the solvent bath is folded into analytic lineshape factors via a second-order cumulant expansion, justified by the assumption that the system-bath coupling commutes with the system Hamiltonian.

What would settle it

Run the multi-D2 calculation for an off-diagonal case such as φ=0.4, J=g=0, N=10, and compare the exciton probability and the 2D spectra with HEOM results for the same parameters; disagreement at the level of the diagonal-case benchmark (differences two orders of magnitude smaller than Pex) would refute the claim that the method is fast and accurate for off-diagonal spectroscopy.

Watch

Extended reading notes

Core claim

The central discovery is that the multi-D1 and multi-D2 Ansätze converge toward the exact Schrödinger dynamics of the Holstein polaron as the multiplicity M grows, and that the multi-D2 Ansatz can be used to evaluate nonlinear response functions. For the diagonal-coupling case, the multi-D1 Ansatz with M=8 reproduces HEOM exciton probabilities to within about two orders of magnitude, and energy components from multi-D1 with M=5 agree with multi-D2 with M=32. For off-diagonal coupling, the multi-D2 Ansatz reveals an increased effective mass when both transfer and off-diagonal coupling act, and yields linear absorption spectra whose zero-phonon line matches the k=0 polaron energy band. The claimed first 2D spectra for off-diagonal coupling show a crossover from a single diagonal peak at φ=0.1 to a vibronic multi-peak structure at φ=0.4.

Load-bearing premise

The central claim rests on the assumption that the multi-D2 trial state faithfully captures off-diagonal exciton-phonon dynamics and that the solvent bath can be folded into analytic lineshape factors, because the paper's exact comparison against HEOM covers only a diagonal-coupling case.

Editorial extensions

If this is right

  • Multi-D1 with modest M (about 4–8) can replace single Davydov trial states in the diagonal-coupling Holstein model, eliminating artifacts such as spurious self-trapping at long times.
  • Linear absorption spectra computed variationally locate the zero-phonon line at the k=0 polaron energy, giving a direct way to read ground-state band energies from dynamics.
  • The multi-D2 calculations predict that off-diagonal coupling can localize an exciton even when the bare transfer integral would delocalize it, which should be observable as reduced mobility in molecular crystals.
  • The first 2D spectra for off-diagonal coupling imply a spectroscopic fingerprint: weak off-diagonal coupling gives one peak, strong coupling gives vibronic multi-peaks, with population cascading to lower energies as the waiting time grows.
  • The same variational pipeline should extend to larger lattices and to simultaneous diagonal-plus-off-diagonal coupling without changing the method.

Reading between the lines

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

  • A natural next test, not performed in the paper, is a HEOM calculation of off-diagonal-coupling dynamics: the only exact benchmark shown is diagonal, and the paper's own deviation measure for off-diagonal multi-D1 is larger (σ ≈ 0.54 at M=6), so an exact off-diagonal comparison would sharpen or bound the claim.
  • Because the 2D calculation separates the bath through cumulant lineshape factors that require the system-bath coupling to commute with the system Hamiltonian, the same approach may need modification when the bath couples to off-diagonal transfer degrees of freedom; a direct nonlinear-response HEOM check would test whether the multi-D2 2D peaks survive.
  • If the single-peak to multi-peak transition is robust, the ratio of peak splittings at fixed population time could be used experimentally to estimate the off-diagonal coupling strength φ in J-aggregate-like systems.
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 / 6 minor

Summary. The paper develops time-dependent variational dynamics for the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling, using superpositions of Davydov D1 and D2 trial states called multi-D1 and multi-D2 Ansätze. It derives equations of motion from the Dirac-Frenkel variational principle, introduces a relative-deviation error measure, and benchmarks multi-D1 dynamics against HEOM for a diagonal-coupling case. It then computes linear absorption spectra and presents 2D spectra for off-diagonal coupling using the multi-D2 Ansatz, claiming a transition from a single peak for weak off-diagonal coupling to a vibronic multi-peak structure for strong coupling.

Significance. If the central computational claims are correct, the multi-D trial states provide an efficient route to polaron dynamics and nonlinear spectroscopy in a regime where single Davydov Ansätze are known to fail. The paper's diagonal-coupling validation is a genuine strength: exciton probabilities agree with HEOM to roughly two orders of magnitude below the signal, and the relative deviation decreases with multiplicity. The linear-absorption zero-phonon line is also cross-checked against variational energy-band calculations. However, the distinctive new result, the off-diagonal 2D spectra, lacks an independent benchmark and rests on response-function formulas that are not fully specified, so the significance is currently conditional.

major comments (3)
  1. [Sec. III D; Fig. 4; Sec. III B] The 2D spectra in Fig. 13, which are the paper's central new claim, are not validated by any independent method or by a convergence study for the multi-D2 multiplicity at the parameters used (J=g=W=0, N=10, φ=0.1 and 0.4). The only exact benchmark in the paper is for diagonal coupling (Sec. III B, Figs. 5-6). For off-diagonal coupling, the reported relative deviation remains as large as σ=0.54 at the largest multiplicity shown in Fig. 4, with J=W=g=0 and φ=0.4, and the text's assertion that multi-D2 gives 'considerable improvements' is delegated to Ref. [47] rather than demonstrated in this parameter regime. The multiplicity M used for Fig. 13 is not stated. Without a convergence test in M or an external reference solution for the off-diagonal dynamics, the quantitative content of the 2D spectra is unsupported.
  2. [Appendix D, Eqs. (D5), (D10), (D11)] The derivation of the nonlinear response functions is not self-contained. Eq. (D10) approximates e^{-iH_S t}|n>|0>_ph, but the initial-site label n is not carried into the right-hand side; Eq. (D11) then introduces superscripts n on ψ and λ (e.g., ψ_{j n'}^{n*}(T), λ_{jq}^{n}(T)) without defining them. The phase factors e^{iω_q t}, e^{iω_q(t+T)}, and e^{-iω_q T} in the coherent-state overlaps do not follow from Eq. (D10) as written, and Eq. (D5) displays system propagators without the imaginary unit in several terms (R1: e^{-H_S(t+T+τ)}; R2 and R3 second factors; R4 second factor). These formula-level inconsistencies affect the central first-time 2D claim and must be corrected with a full derivation or with defining expressions for all labels and phases.
  3. [Sec. III B, Fig. 8] The off-diagonal dynamics that motivate the use of multi-D2 are not independently benchmarked. In Fig. 8 the multi-D2 and single-D2 results are compared with each other, which only establishes that the two variational approximations differ; it does not show that the multi-D2 result, such as the reported localization at φ=0.1, is closer to the exact dynamics. Since this off-diagonal accuracy is the premise for the 2D calculation, an independent test for at least one off-diagonal parameter set (e.g., off-diagonal HEOM or a converged wave-function/TD-DMRG calculation) is required.
minor comments (6)
  1. [Fig. 3(a) caption; Sec. III A] The caption of Fig. 3(a) states g=1, while the text in Sec. III A states g=0.1 for the same panel; please reconcile the value of the diagonal coupling strength.
  2. [Sec. III D, Fig. 13] The multiplicity M used for the 2D spectra in Fig. 13 is not reported; please state it and include a convergence check in M for at least one panel.
  3. [Abstract; Sec. III B] The phrase 'perfect agreement' with HEOM is stronger than the data support; the text reports a difference two orders of magnitude below the signal, so a quantitative wording would be more accurate.
  4. [Fig. 4; Sec. III A] The text states that σ=0.54 corresponds to M=6, but the x-axis 1/M=0.2 corresponds to M=5; please correct the multiplicity labeling in the figure or in the text.
  5. [Eq. (14), Sec. II B] Please clarify whether the denominator in σ is a time average of N_err(t) over [0,t_max], and specify the units of Δ(t) and N_err(t) so that the dimensionless character of σ is transparent.
  6. [General editorial] Minor typographical and reference issues: 'Lorenz' in the Fig. 12 caption should be 'Lorentzian'; 'the the' appears in the Fig. 4 caption; Ref. [53] lists the year as '2012)'; and the second line of the abstract contains a misplaced space in 'multidim ensional'.

Circularity Check

1 steps flagged · score 4.0 of 10

Off-diagonal 2D-spectroscopy claim leans on same-group Ref. [47] for multi-D2 accuracy; diagonal benchmarks are independent, so overall circularity is moderate.

  1. self citation load bearing [Section III A, paragraph after Fig. 4 (off-diagonal validity of multi-D2)]
    "For off-diagonal coupling, considerable improvements in accuracy can be achieved by utilizing multi-D2 with the increase of multiplicity M (see discussions in Ref. [47])."

    This sentence is the only support offered in this paper for the accuracy of the multi-D2 propagator in the off-diagonal regime used for the headline 2D spectra (Sec. III D, phi = 0.1 and 0.4). Ref. [47] is authored by the same group (N. J. Zhou, Z. K. Huang, J. F. Zhu, V. Chernyak, and Y. Zhao; current authors include N. Zhou, Z. Huang, and Y. Zhao). The paper's own off-diagonal error metric for multi-D1 gives sigma = 0.54 at M = 6 (Fig. 4), and no multi-D2 convergence data or HEOM comparison for off-diagonal dynamics is shown here. Thus the premise that the Appendix D response functions are computed with an accurate propagator is not established in this paper but is inherited from a self-citation, making the first-time off-diagonal 2D claim load-bearing on that citation.

full rationale

There is no equation-level circularity of the self-definitional or fitted-input kind: the multi-D1 and multi-D2 equations of motion are derived from the Dirac-Frenkel variational principle, and no fitted parameter is renamed as a prediction. The diagonal-coupling exciton probabilities are checked against an independent HEOM benchmark (Figs. 5 and 6), and the linear absorption zero-phonon line is compared with variationally computed energy bands, which provides independent content. However, the paper's flagship novelty, the first 2D spectra for off-diagonal exciton-phonon coupling, relies on the accuracy of the multi-D2 propagator in the off-diagonal regime. For that regime the paper provides no HEOM or other exact nonlinear-response benchmark; instead, the only justification cited is Ref. [47], a prior paper by the same authors. Since the 2D peak-structure distinction (single peak for weak coupling, vibronic multi-peak for strong coupling) is produced by that unbenchmarked propagator, the central claim is partly load-bearing on a self-citation. This warrants a moderate circularity score of 4 rather than a higher one, because the diagonal dynamics and the numerical construction of the 2D spectra are independently substantive.

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

Most numbers in the paper are physical model parameters (J, g, phi, W, omega_0, N) that are inputs rather than fitted quantities. The free-parameter burden rests on the multiplicity M, the spectral broadening width, the bath parameters, and the artificial noise amplitude; none is determined by data. The main axioms are the variational approximation, the Holstein model, the exactness of truncated HEOM, and the cumulant factorization used for 2D spectra. No new physical entities are introduced.

free parameters (4)
  • Multiplicity M = M = 1, 4, 5, 8, 16, 32 in different runs
    M controls the flexibility of the trial state; accuracy claims depend on M being large, and no automated convergence criterion is given for new parameter regimes.
  • Spectral damping factor = 0.08 omega_0
    Used to broaden the linear absorption lines in Sec. III C; chosen by hand and affects line widths, not stated peak positions.
  • Drude-Lorentz bath parameters = eta = 0.1, beta = 5, gamma = 0.02
    Chosen for the 2D spectra in Appendix D; these control dephasing and peak shapes and are not derived from experimental data.
  • Initial-state noise amplitude = 1e-5
    Uniform noise added to variational parameters to avoid singularities; chosen ad hoc and not justified by a convergence study.
assumptions (4)
  • domain assumption The Dirac-Frenkel time-dependent variational principle with a finite coherent-state manifold gives a controlled approximation to the exact Schrodinger dynamics.
    Used to derive equations of motion in Sec. II B; accuracy is monitored by the relative deviation sigma, but sigma itself is not an external check.
  • domain assumption The extended Holstein Hamiltonian with linear phonon dispersion and N = 10 or 16 periodic sites captures the essential physics of the target molecular crystals.
    The model is introduced in Sec. II A and all numerical results are obtained within it; no experimental data are used.
  • domain assumption Truncated hierarchy equations of motion provide a numerically exact reference for the diagonal-coupling dynamics.
    HEOM is used as the benchmark in Sec. III B; truncation depth parameters are not reported fully.
  • domain assumption For the 2D spectra, the bath can be factored out through second-order cumulant lineshape functions because the system-bath coupling commutes with the system Hamiltonian.
    Appendix D factorizes the nonlinear response into variational system propagators and analytic lineshape factors; this standard approximation is not tested against an exact nonlinear-response calculation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Fast, accurate simulation of polaron dynamics and multidimensional spectroscopy by multiple Davydov trial states." pith.science (2026). https://pith.science/paper/NM7O6MED

@misc{pith2026190809243,
  author       = {Pith},
  title        = {Pith review of: Fast, accurate simulation of polaron dynamics and multidimensional spectroscopy by multiple Davydov trial states},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NM7O6MED}},
  note         = {Machine review of arXiv:1908.09243}
}
abstract

By employing the Dirac-Frenkel time-dependent variational principle, we study the dynamical properties of the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling. A linear combination of the Davydov D$_1$ (D$_2$) Anstaz, referred to as the multi-D$_1$ Ansatz (multi-D$_2$ Ansatz), is used as the trial state with enhanced accuracy but without sacrificing efficiency. The time evolution of the exciton probability is found to be in perfect agreement with that of the hierarchy equations of motion, demonstrating the promise the multiple Davydov trial states hold as an efficient, robust description of dynamics of complex quantum systems. In addition to the linear absorption spectra computed for both diagonal and off-diagonal cases, for the first time, $2$D spectra have been calculated for systems with off-diagonal exciton-phonon coupling by employing the multiple $D_2$ Ansatz to compute the nonlinear response function, testifying to the great potential of the multiple $D_2$ Ansatz for fast, accurate implementation of multidimensional spectroscopy. It is found that the signal exhibits a single peak for weak off-diagonal coupling, while a vibronic multi-peak structure appears for strong off-diagonal coupling.

Figures

Figures reproduced from arXiv: 1908.09243 by the authors.

Figure 1
Figure 1. (a) Schematic of the Holstein ring. A simplified [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. The energies of the exciton, phonon, and exciton [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. (a) The relative deviation σ of the multi-D1 Ansatz in a 16-site molecular ring is displayed as a function of 1/M representing the inverse of the multiplicity. The set of param￾eters J = 0.1, g = 1, W = 0.5 and φ = 0 is used. Moreover, the relative deviation σ for the diagonal coupling case is also plotted as a function of the transfer integral J in (b) and di￾agonal coupling strength g in (c). In both of them, the … view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: The relative deviation σ from the multi-D1 Ansatz is displayed as a function of 1/M for the the off-diagonal cou￾pling case with the strength φ = 0.4, and other parameters J = W = g = 0 are set. B. Exciton probabilities and phonon displacements Dynamical properties of …
Figure 5
Figure 5. Figure 5: Time evolution of the exciton probability [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: Time evolution of the exciton probability [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 8
Figure 8. Figure 8: Time evolution of the exciton probability [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 10
Figure 10. Figure 10: Polaron energy bands Ek/ω0 are calculated varia￾tionally using the Delocalized D1 Ansatz (solid line) and the Toyozawa Ansatz (open circles) for the case of g = 0.2, J = 0.1, W = 0.1 and φ = 0. The position of zero-phonon line ωm/ω0 is marked by the dashed line, consi…
Figure 11
Figure 11. Figure 11: denote the positions of the zero-phonon lines (ωm/ω0 = −0.08, 0.369, −0.956 and −1.93). For strong off-diagonal coupling, such as the case of φ = 1, the linear absorption spectra, shown in [PITH_FULL_IMAGE:figures/full_fig_p009_11.png]
Figure 13
Figure 13. Figure 13: 2D spectra of the molecular ring for off-diagonal [PITH_FULL_IMAGE:figures/full_fig_p010_13.png]
Figure 12
Figure 12. Figure 12: Linear absorption spectrum F(ω) of the DM=16 2 Ansatz is displayed for the off-diagonal coupling case with φ = 1. In the inset, the power-law and Lorenz fittings are given in the log-log scale with the dashed and dotted lines, respectively. ωτ ωt T=0 −3 −2 −1 0 1 2 3 …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

66 extracted references · 63 canonical work pages

  1. [47]

    M. J. ˇSkrinjar, D. V. Kapor and S. D. Stojanovi´ c , Phys. Rev. A 38, 6402 (1988), and references therein

    work page 1988

  2. [1]

    obtained by the multi-D 1 Ansatz with M = 4 are pre- sented in Figs

    5 and φ = 0. obtained by the multi-D 1 Ansatz with M = 4 are pre- sented in Figs. 7(b) and 7(d), respectively. Quite obvi- ous difference is found in the excitonic behavior for the two cases when t/(2π/ω0)> 3. To be specific, the exci- ton probability calculated by the single D 1 Ansatz stag- gers around two sites in the ring before being eventually trapped...

    work page

  3. [2]

    1 and φ = 0

    1, W = 0 . 1 and φ = 0. The position of zero-phonon line ω m/ω 0 is marked by the dashed line, consistent with the val- ues of Ek=0/ω 0. A lattice of N = 16 sites is used in calcula- tions. In the inset, the transition moment Pk is plotted as a function of the crystal momentum k. responding polaron energy band Ek =⟨Ψ k| ˆH|Ψ k⟩ cal- culated. Variations ca...

    work page

  4. [3]

    multi-D 1 Ansatz

    1, 0. 2, 0. 3 and 0 . 4, respectively. Other parameters g = W = J = 0 and N = 16 are set. The vertical dash-dotted lines indicate locations of zero-phonon lines. 10 panel). For simplicity, we set J = g = W = 0, and adopt the toy model of J-aggregates with the tangen- tial (head-to-tail) orientations of the transition dipoles. We first consider weak off-diag...

    work page 2009

  5. [4]

    Z. An, C. Q. Wu, and X. Sun, Phys. Rev. Lett. 93, 216407 (2004)

    work page 2004

  6. [5]

    Zheng, J

    B. Zheng, J. Wu, W. Sun, and C. Liu, Chem. Phys. Lett. 425, 123 (2006)

    work page 2006

  7. [6]

    X. Liu, K. Gao, J. Fu, Y. Li, J. Wei, and S. Xie, Phys. Rev. B 74, 172301 (2006)

    work page 2006

  8. [7]

    Bronner and P

    C. Bronner and P. Tegeder, Phys. Rev. B 89, 115105 (2014)

    work page 2014

Show all 66 references
  1. [8]

    Timrov, T

    I. Timrov, T. Kampfrath, J. Faure, N. Vast, C. R. Ast, C. Frischkorn, M. Wolf, P. Gava, and L. Perfetti, Phys. Rev. B 85, 155139 (2012)

    2012

  2. [9]

    Chikan and D

    V. Chikan and D. F. Kelleya, J. Phys. Chem. 117, 8944 (2002)

    2002

  3. [10]

    V. I. Klimov, A. A. Mikhailovsky, D. W. McBranch, C. A. Leatherdale, and M. G. Bawendi, Phys. Rev. B 61, R13349 (2000)

    2000

  4. [11]

    S. H. Kim, R. H. Wolters, and J. R. Heath, J. Chem. Phys. 105, 7957 (1996)

    1996

  5. [12]

    Sauer, Annu

    K. Sauer, Annu. Rev. Phys. Chem. 30, 155 (1979)

    1979

  6. [13]

    Renger, V

    T. Renger, V. May, and Oliver K¨ uKhn, Phys. Rep. 343, 137 (2001)

    2001

  7. [14]

    R. E. Blankenship, Molecular Mechanisms of Photosyn- thesis (Blackwell Science, Oxford/Malden, 2002)

    2002

  8. [15]

    van Grondelle and V

    R. van Grondelle and V. I. Novoderezhkin, Phys. Chem. Chem. Phys. 8, 793 (2006)

    2006

  9. [16]

    L. P. Chen, P. Shenai, F. L. Zheng, A. Somoza, and Y. Zhao, Molecules. 20, 15224, (2015)

    2015

  10. [17]

    Tomimoto, H

    S. Tomimoto, H. Nansei, S. Saito, T. Suemoto, J. Takeda, and S. Kurita, Phys. Rev. Lett. 81, 417 (2000)

    2000

  11. [18]

    S. L. Dexheimer, A. D. Van Pelt, J. A. Brozik, and B. I. Swanson,Phys. Rev. Lett. 84, 4425 (2000)

    2000

  12. [19]

    Sugita, T

    A. Sugita, T. Saito, H. Kano, M. Yamashita, and T. Kobayashi, Phys. Rev. Lett. 86, 2158 (2001)

    2001

  13. [20]

    F. X. Morrissey and S. L. Dexheimer, Phys. Rev. B 81, 094302 (2010)

    2010

  14. [21]

    A. S. Alexandrov and Sir N. Mott, Polarons and Bipo- larons (World Scientific, London, 1995)

    1995

  15. [22]

    F. M. Peeters and J. T. Devreese, Solid State Phys. 38, 81 (1984)

    1984

  16. [23]

    Ranninger, in Polarons in Bulk Materials and Sys- tems with Reduced Dimensionality , edited by G

    J. Ranninger, in Polarons in Bulk Materials and Sys- tems with Reduced Dimensionality , edited by G. Iadonisi, J. Ranninger, and G. De Filippis, International School of Physics Enrico Fermi, (IOS Press, Amsterdam) 161, 1(2006)

    2006

  17. [24]

    Holstein, Ann

    T. Holstein, Ann. Phys. (N.Y.) 8, 325 (1959)

    1959

  18. [25]

    Holstein, Ann

    T. Holstein, Ann. Phys. (N.Y.) 8, 343 (1959)

    1959

  19. [26]

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

    1979

  20. [27]

    L. A. Dissado and S. H. Walmsley, Chem. Phys. 86, 375 (1984)

    1984

  21. [28]

    Sumi, Chem

    H. Sumi, Chem. Phys. 130, 433 (1989)

    1989

  22. [29]

    Y. Zhao, D. W. Brown, and K. Lindenberg, J. Chem. Phys. 100, 2335 (1994)

    1994

  23. [30]

    D. M. Chen, J. Ye, H. J. Zhang, and Y. Zhao, J. Phys. Chem. B 115, 5312 (2011)

    2011

  24. [31]

    B. Luo, J. Ye, C. B. Guan and Y. Zhao, Phys. Chem. Chem. Phys. 12, 6045 (2010)

    2010

  25. [32]

    J. Sun, B. Luo, and Y. Zhao, Phys. Rev. B 82, 014305 (2010)

    2010

  26. [33]

    R. W. Munn and R. Silbey, J. Chem. Phys. 83, 1843(1985); 83, 1854 (1985)

    1985

  27. [34]

    Zhao, Doctoral thesis, University of California, Sa n Diego, (1994)

    Y. Zhao, Doctoral thesis, University of California, Sa n Diego, (1994)

    1994

  28. [35]

    Y. Zhao, G. Q. Li, J. Sun, and W. H. Wang, J. Chem. Phys. 129, 124114 (2008)

    2008

  29. [36]

    Q. Liu, Y. Zhao, W. Wang, and T. Kato, Phys. Rev. B 79, 165105 (2009)

    2009

  30. [37]

    L. C. Ku, A. Trugman, Phys. Rev. B 75, 014307 (2007)

    2007

  31. [38]

    Meier, Y

    T. Meier, Y. Zhao, V. Chernyak, and S. Mukamel, J. Chem. Phys. 107, 3876 (1997)

    1997

  32. [39]

    Y. Zhao, D. W. Brown, and K. Lindenberg, J. Chem. Phys. 106, 2728 (1997)

    1997

  33. [40]

    Y. Zhao, D. W. Brown, and K. Lindenberg, J. Chem. Phys. 107, 3159 (1997); D. Brown, K. Lindenberg, and Y. Zhao, ibid. 107, 3179 (1997)

    1997

  34. [41]

    J. Sun, L. W. Duan, and Y. Zhao, J. Chem. Phys. 138, 174116 (2013)

    2013

  35. [42]

    N. J. Zhou, L.P. Chen, Y. Zhao, D. Mozyrsky, V. Chernyak, and Y. Zhao, Phys. Rev. B 90, 155135 (2014)

    2014

  36. [43]

    P. A. M. Dirac, Proc. Cambridge Philos. Soc. 26, 376 (1930); J. Frenkel, Wave Mechanics (Oxford University Press, 1934)

    1930

  37. [44]

    Y. Zhao, B. Luo, Y. Y. Zhang, and J. Ye, J. Chem. Phys. 137, 084113 (2012)

    2012

  38. [45]

    A. S. Davydov and N. I. Kislukha, Zh. Eksp. Teor. Fiz,. 71,1090 (1976) [Sov. Phys. JETP. 44, 571 (1976)]

    1976

  39. [46]

    A. S. Davydov, Solitons in Molecular Systems (Reidel, Dordrecht, 1985)

    1985

  40. [48]

    F¨ orner, J

    W. F¨ orner, J. Phys.: Condens. Matter 5, 3897 (1993); Phys. Rev. B 53, 6291 (1996)

    1993

  41. [49]

    L. C. Hansson, Phys. Rev. Lett. 73, 2927 (1994)

    1994

  42. [50]

    N. J. Zhou, Z. K. Huang, J. F. Zhu, V. Chernyak, and Y. Zhao, J. Chem. Phys. 143, 014113 (2015)

    2015

  43. [51]

    Brixner, J

    T. Brixner, J. Stenger, H. M. Vaswani, M. Cho, R. E. Blankenship, G. R. Fleming, Nature 434, 625 (2005)

    2005

  44. [52]

    G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, Y. C. Chung, R. E. Blankenship, G. R. Fleming. Nature 446, 782 (2007)

    2007

  45. [53]

    Collini, C

    E. Collini, C. Y. Wong, K. E. Wilk, P. M. G. Curmi, P. Brumer, G. D. Scholes. Nature, 463, 644 (2010)

    2010

  46. [54]

    Panitchayangkoon, D

    G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. R. Caram, E. Harel, J. Z. Wen, R. E. Blankenship, G. S. Engel. Proc. Natl. Acad. Sci. USA 107, 12766 (2010)

    2010

  47. [55]

    J. A. Myers, K. L. M. Lewis, F. D. Fuller, P. F. Tekavec, C. F. Yocum, J. P. Ogilvie. J. Phys. Chem. Lett, 1, 2774 (2010)

    2010

  48. [56]

    K. L. M. Lewis, J. P. Ogilvie. J. Phys. Chem. Lett, 3, 503 (2012)?

    2012

  49. [57]

    Romero, R

    E. Romero, R. Augulis, V. I. Novoderezhkin, M. Ferretti , J. Thieme, D. Zigmantas, R. van Grondelle. Nat, Phys. 10,676 (2014)

    2014

  50. [58]

    Mukamel, Principles of Nonlinear Optical Spectroscopy (Oxford University Press, New York, 1995)

    S. Mukamel, Principles of Nonlinear Optical Spectroscopy (Oxford University Press, New York, 1995)

    1995

  51. [59]

    T. D. Huynh, K. W. Sun, M. F. Gelin, and Y. Zhao, J. Chem. Phys. 139, 104103 (2013)

    2013

  52. [60]

    K. W. Sun, M. F. Gelin, V. Y. Chernyak, and Y. Zhao, J. Chem. Phys. 142, 212448 (2015)

    2015

  53. [61]

    Tanimura, R

    Y. Tanimura, R. Kubo, J. Phys. Soc. Jpn, 58, 101 (1989)

    1989

  54. [62]

    Tanimura, Phys

    Y. Tanimura, Phys. Rev. A, 41, 6676 (1990)

    1990

  55. [63]

    Tanimura, J

    Y. Tanimura, J. Phys. Soc. Jpn, 75, 082001 (2006). 16

    2006

  56. [64]

    Ishizaki, Y

    A. Ishizaki, Y. Tanimura, J. Phys. Soc. Jpn, 74, 3131 (2005)

    2005

  57. [65]

    L. P. Chen, Y. Zhao, and Y. Tanimura, J. Phys. Chem. Lett. 6, 3110 (2015)

    2015

  58. [66]

    L. P. Chen, M. F. Gelin, W. Domcke, and Y. Zhao, J. Chem. Phys. 142, 164106 (2015)

    2015

Pith tools

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