{"id":"eea73475-5106-4d6c-a7d4-c3c6278911dd","arxiv_id":"2411.17021","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Polaron diffusion in the dispersive Holstein model is claimed to onset only exponentially slowly with lattice size, a scaling that the paper's own generalized master equation cannot produce.","lead":"This paper uses a space-and-time-local memory-kernel method to simulate polaron motion in the Holstein model on lattices far larger than exact methods can reach, and reports that the onset of diffusive transport is exponentially delayed with system size. The headline scaling claim is not derived and appears inconsistent with the method's own Markovian long-time propagation, while the large-lattice predictions are not checked against any exact calculation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central exponential-delay claim is contradicted by the paper's own STL-GME propagation: a fixed translation-invariant rate matrix cannot produce 1 - alpha_max decaying exponentially in N, since the slowest finite-size mode gives an algebraic ~1/N^2 deficit.","rationale":"I read the paper as claiming a physically surprising result: nonequilibrium polaron transport in the homogeneous dispersive Holstein model remains subdiffusive until exponentially large times and system sizes, with Fig. 2d/e as the direct evidence. The most load-bearing premise is that the observed log-linear decline of 1 - alpha_max with N and tau_R reflects the true model. My concern is that this premise contradicts the paper's own STL-GME algorithm. Once the generator has reached its lifetime, propagation is a fixed finite-range stochastic matrix; such matrices have a quadratic dispersion near q=0 and hence a boundary mode with rate ~D/N^2. Consequently the maximum diffusion exponent before boundary reflection can only miss unity by an algebraic amount ~1/N^2, and before that 1 - alpha(t) ~ 1/t. An exponential deficit would require a mechanism to push the boundary onset to exponentially long times; none is present in Eq. (M2). The paper does not derive the exponential form, give error bars on the fits, or validate against exact dynamics for N > 30. This is not a disagreement with consensus; it is an internal consistency issue for the central claim. I therefore affirm the reader's REJECT verdict, while recognizing that the STL-GME idea, the 1D-to-2D outer-product observation, and the localization results may still be useful. The concrete check is to recompute alpha_max(N) from the same generator, either by direct propagation or by diagonalizing the N-site rate matrix; an algebraic scaling would confirm that the exponential claim is an artifact of the construction or analysis rather than a physical property of the model.","tokens_in":17706,"tokens_out":13847,"duration_ms":149727,"concrete_test":"Reconstruct U(tauU=90 fs) from the 10-site HEOM reference exactly as in Methods B-C. Then, for each N in {100, 200, 400, 800, 1500}, diagonalize the padded, renormalized periodic N-site generator and compute alpha_max(N) from the exact mode expansion, propagating at least to a few times N^2/D, where D is obtained from the q -> 0 Fourier eigenvalue of the same generator. Plot 1 - alpha_max(N) on both log-linear and log-log axes. If the log-log slope is approximately -2, the exponential claim in Fig. 2d is an artifact; a truly exponential fall-off would require a mechanism absent from Eq. (M2). As a cross-check, compare STL-GME alpha(t) against exact HEOM for N=30 and N=50 over accessible times to verify that the padded generator reproduces the finite-size rise time and not just the small-system plateau.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result is the exponential delay of diffusive onset (Fig. 2d/e and Introduction). This claim is internally inconsistent with the STL-GME construction used to produce it. For t > tauU = 90 fs, Eq. (M2) propagates C by repeated application of the fixed, translation-invariant generator U(tauU). Any finite-range translation-invariant generator with a nonzero variance has a Fourier eigenvalue lambda(q) = 1 - D q^2 delta_t + O(q^4), so on the infinite lattice MSD(t) = 2Dt + c and 1 - alpha(t) ~ c/(2Dt). On a ring of N sites the slowest nonconstant mode has rate D(2*pi/N)^2; boundary effects therefore set in at t ~ N^2/D, giving 1 - alpha_max ~ C/N^2 (algebraic), not exp(-cN). An exponentially small 1 - alpha_max would require the maximum to occur at t ~ exp(cN), which cannot happen for a local rate matrix with finite speed and diffusive spreading. The paper offers no derivation of the exponential form, reports no error bars on the log-linear fits, and provides no exact large-N check. The observed exponential scaling is therefore most plausibly an artifact of the 10-site reference used to build U(tauU), the zero-padding/renormalization in Methods C, or the choice of alpha_max sampling window, rather than a property of the dispersive Holstein model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a space- and time-local generalized master equation (STL-GME) method for simulating polaron transport in the dispersive Holstein model. The method constructs a memory kernel U(t) from small-lattice HEOM reference dynamics, truncates it in space and time, and then extends the generator to much larger lattices by zero padding and translational invariance, propagating with the fixed matrix U(τU) for long times (Eq. M2). The central claim is that the maximum of the scaling exponent α (defined through MSD ∝ t^α) approaches unity exponentially slowly as a function of lattice size, so that the onset of diffusion is delayed to exponentially large times and system sizes. The paper also reports that 1D transport, via an outer-product construction, quantitatively reproduces 2D transport coefficients, and that anisotropic energy landscapes in periodic 2D lattices direct polaron motion.","tokens_in":18029,"tokens_out":7769,"duration_ms":76027,"significance":"If the exponential-delay claim were correct, it would challenge the standard extrapolation of finite-lattice simulations to the thermodynamic limit and would have broad implications for the interpretation of nonequilibrium polaron transport. The STL-GME idea of exploiting joint temporal and spatial memory locality is attractive and the small-system agreement between GME and HEOM (errors below 1e-7 for τU and 5e-6 for dU) is credible. However, the central claim is not merely unsupported; it is mathematically incompatible with the paper's own propagation rule, Eq. (M2), as detailed below. The 2D transport and localization results are of interest but do not rescue the headline claim, which is the basis for the paper's stated significance.","major_comments":[{"comment":"The exponential decay of 1−α_max with N is contradicted by the STL-GME propagation rule itself. For t > τU = 90 fs, the dynamics are generated by repeated application of the fixed, finite-range, translationally invariant matrix U(τU) via C(t+nδt) = [U(τU)]^n C(t=τU). Any such matrix has Fourier eigenvalues of the form λ(q) = 1 − D q^2 δt + O(q^4), where D > 0 is the diffusion constant. Consequently, on an infinite lattice MSD(t) = 2Dt + c, which gives 1−α(t) = c/(2Dt) + O(t^−2); on a ring of N sites the slowest nonconstant mode has relaxation time τ_R ∼ N^2/D, so α_max is reached at t ∼ N^2/D and 1−α_max ∼ const/N^2. This is an algebraic, not exponential, dependence on N. The authors provide no derivation of the exponential form and no independent large-N check; the log-linear behavior in Fig. 2(d)–(e) is therefore most plausibly an artifact of the 10-site reference, the zero-padding/renormalization in Methods C, or the choice of the α_max sampling window. Please provide a direct test: plot log(1−α_max) versus log N for the STL-GME with the fixed generator U(τU) and report the slope; if it is not −2, explain explicitly how the generator violates the small-q expansion.","section":"Polaron transport; Methods, Eq. (M2)"},{"comment":"The memory cutoffs τU = 90 fs and dU = 4r0 are determined by comparing GME predictions against the same small-system HEOM reference from which U was constructed. This establishes self-consistency on the 10-site system, but it does not validate the zero-padded extension to M = 1500 sites propagated to 1.5 ns. The claim that the kernel is exactly zero beyond dU and exactly constant after τU is a strong truncation assumption; if the true kernel has slowly decaying tails, or if truncation error accumulates under repeated powers [U(τU)]^n, the extrapolated α_max(N) is uncontrolled. The paper repeatedly labels the STL-GME results as 'exact' (Abstract, Introduction, Outlook), yet no exact or independent large-N benchmark is provided. This is load-bearing because the central quantitative claim rests entirely on the extrapolated large-N dynamics.","section":"Methods B and Figs. S1–S2"},{"comment":"The exponential fits lack error bars, fit ranges, and residuals. The text states that for γ = 400 cm^−1 the asymptotic behavior is evident for N = 600–1000, while for γ = 1000 cm^−1 it only arises for N = 1000–1500, but each curve contains only a handful of points and no statistical uncertainty is reported. Given the algebraic prediction of Eq. (M2) described above, the authors should show whether the data actually follow 1−α_max ∼ N^−2 or demonstrate, with a concrete mathematical argument, that U(τU) produces an exponentially small spectral gap that depends on N. Without such evidence, the exponential-delay claim is unsupported.","section":"Fig. 2(d)–(e)"}],"minor_comments":[{"comment":"In the paragraph on periodic 2D lattices, 'another form blue site' should read 'another from blue site'.","section":"Methods D"},{"comment":"The caption begins with 'PPolaron', which should be 'Polaron'.","section":"Fig. S5 caption"},{"comment":"The phrase 'logarithmic behavior of 1−αmax ... scales linearly' is ambiguous; please specify the axes of the inset and state explicitly that the fit is log(1−αmax) versus N.","section":"Fig. 2(d) caption and main text"},{"comment":"The error threshold is given as 1e-7 per element for τU in Methods B, but Fig. S2 and the surrounding text state a threshold of 5e-6 per element for dU. Please clarify whether these are different thresholds and why.","section":"Methods B and Fig. S2"}],"recommendation":"reject","confidential_remarks":"The central claim of the paper is not merely weakly supported; it is mathematically incompatible with the authors' own propagation formula, Eq. (M2). Because the STL-GME after τU is a fixed finite-range rate matrix, the scaling of 1−α_max with N must be algebraic (∼1/N^2), not exponential. This is a fundamental issue that cannot be remedied by adding benchmarks or error bars within the current framework; it requires either a different propagation scheme or a reformulation of the main claim. I therefore recommend rejection, with the caveat that the 2D transport and localization results might be publishable separately if the exponential-delay claim is removed and the method is validated against an independent large-system reference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: don't take the \"exponentially delayed diffusion\" claim at face value. It is contradicted by the paper's own propagation scheme, and the surprising numerical trend is very likely an artifact. What is genuinely useful is the space-time-local memory construction, and the observation that 1D transport determines 2D transport; those deserve attention even if the headline result falls.\n\nThe small-system validation is the strongest part. The STL-GME reproduces HEOM to about 1e-7 (tauU) and 5e-6 (dU), with a clear protocol for choosing cutoffs, and the 10-site reference extended to 30 and 100 sites looks faithful. That gives real ground under the method, within its tested regime. The convergence checks on HEOM depth, Matsubara terms, and timestep are reported, which is more than many papers do. The 1D-outer-product equivalence for homogeneous 2D transport is empirically useful, and the supplementary difference plots support it.\n\nThe soft spot is the central claim. For t > tauU, Eq. M2 propagates C by repeated application of a fixed, translation-invariant, finite-range generator. On a ring of N sites the slowest nonconstant Fourier mode decays at rate D(2pi/N)^2, so finite-size reflection sets in at tauR ~ N^2/D and 1 - alpha_max should scale as C/N^2, not exp(-cN). An exponentially small 1 - alpha_max would require the maximum of alpha to occur at t ~ exp(cN), which a local generator cannot produce in this setup. The paper offers no derivation of the exponential form and no error bars on the fits; the text says the log-linear scaling is asymptotic, but the mechanism that would produce it is absent.\n\nThe 2D \"exact\" simulations are also not exact: they use the high-temperature approximation K=0, acknowledged in Methods but not in the abstract. That overstatement matters because the abstract says \"exact\" repeatedly. I also note no code or data are shipped; \"available upon request\" is not reproducible for a method paper.\n\nWho is this for? People working on GME memory kernels and polaron transport will find the STL-GME method and the 1D/2D comparison worth reading. The exponential-delay conclusion should not be cited until it is derived or checked against an exact large-system reference. Editor should send it to referees, but with the expectation that the central claim needs either proof or withdrawal. My vote: engage, but don't believe the scaling.","headline":"The exponential finite-size scaling the paper sells is an artifact of its own STL-GME propagator; the useful parts are the method itself and the 1D/2D transport comparison.","tokens_in":18592,"tokens_out":3833,"would_cite":false,"duration_ms":38713,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Nonequilibrium polaron relaxation reaches diffusion only exponentially slowly in time and system size.","keywords":["dispersive Holstein model","polaron transport","space- and time-local GME","nonequilibrium relaxation","subdiffusion","thermodynamic limit","hierarchical equations of motion"],"falsifier":"Directly compute the memory kernel $U(r,t)$ from a larger exact reference, say 30 sites run well past $\\tau_U=90$ fs, and check whether transitions beyond $r=4r_0$ or times beyond $\\tau_U$ are genuinely negligible; alternatively, compare the STL-GME populations on sites farther than $d_U$ from the injection site against exact dynamics on the same lattice. If non-negligible tails appear or the error grows with propagation time, the exponential $1-\\alpha_{\\max}$ scaling is an artifact of the truncation.","tokens_in":17419,"feed_emoji":"⏳","tokens_out":13821,"duration_ms":112443,"temperature":0.7,"pith_summary":"Polarons—charges dragging a local lattice deformation behind them—are central to energy transport in many materials. This paper claims that after a sudden excitation, a polaron in the dispersive Holstein model does not settle into ordinary diffusive motion until exponentially large times and system sizes: the transport exponent $\\alpha$ in $\\mathrm{MSD}\\propto t^{\\alpha}$ approaches 1 only asymptotically, with $1-\\alpha_{\\max}$ falling linearly on a log scale as the lattice grows to 1500 sites and the dynamics run to 1.5 ns. To see this, the authors use a generalized master equation whose memory kernel is local in both time and space, so that a short, small-lattice numerically exact simulation can be zero-padded into a thermodynamically large one. They also find that dilute-limit transport in 1D quantitatively determines transport in 2D after dimensional renormalization, and that anisotropic 2D energy landscapes localize polaron motion along specific directions. If the claims hold, finite-size simulations cannot be casually extrapolated to the diffusive regime, and experimentally observed polaron motion may be probing a very long transient rather than asymptotic transport.","feed_headline":"Polaron diffusion waits exponentially long in large lattices","feed_subtitle":"New simulations out to 1500 sites and 1.5 ns show the approach to diffusion slows exponentially with system size.","key_machinery":"The load-bearing object is the space- and time-local generalized master equation (STL-GME), a generator $U(r,t)$ that updates site populations according to $C(r,t+\\delta t)=\\int dr'\\, U(r-r',t)\\,C(r',t)$, where $r$ is the separation between the initial injection site and the measurement site. The generator is taken to be finite in time—with a lifetime $\\tau_U=90$ fs—and finite in space—with a cutoff $d_U=4r_0$—both measured against a 10-site hierarchical equations of motion (HEOM) reference. Because $d_U<N/2$, the small-lattice generator contains no finite-size contamination, and translational invariance lets one zero-pad it to 1500 sites in 1D and $30\\times30$ sites in 2D; after $\\tau_U$ the generator is a constant matrix, so propagation beyond that time is simple repeated matrix multiplication. The spatial locality is justified by a light-cone picture: by the time the memory dies, a quasiparticle moving at the characteristic speed has traveled at most $d_U\\sim\\tau_U V$, so transitions beyond that distance can be set to zero.","core_discovery":"On its own terms, the paper's central discovery is that $\\alpha_{\\max}$, the maximum of the scaling exponent $\\alpha$ before finite-size reflection sets in, approaches unity exponentially slowly as a function of system size $N$ and of the finite-size reflection time $\\tau_R$, for every phonon decorrelation rate $\\gamma$ studied. A 10-site chain gives $\\alpha_{\\max}\\approx0.78$ near 200 fs, a 30-site chain $\\alpha_{\\max}\\approx0.95$ near 2 ps, and a 100-site chain $\\alpha_{\\max}\\approx0.99$ near 20 ps; extending to 1500 sites and 1.5 ns shows $\\log(1-\\alpha_{\\max})$ decreasing linearly, i.e., diffusion is reached only asymptotically in the thermodynamic limit. The same space-time-local GME applied to homogeneous 2D lattices up to $30\\times30$ shows that the dimensionality-renormalized $d\\mathrm{MSD}/dt$ and $\\alpha$ agree quantitatively with the 1D results, although the full 2D density is not exactly the outer product of the 1D density; the difference is orders of magnitude smaller than the populations and decays over time. In 2D periodic lattices with two site energies, the paper finds that the relative hopping strengths direct polaron localization: a small inter-column hopping confines the polaron along the vertical axis, while a larger cross-hopping restores an elliptical spread.","pith_inferences":["The exponential onset suggests that measured polaron mobility will depend on the observation window, so two experiments on the same material could report different effective diffusion constants simply by probing different times after excitation.","If the space-local memory picture is correct, other short-range carrier-phonon models with similar local couplings may also show exponentially slow onsets of diffusion in their nonequilibrium relaxation.","The cutoff $d_U$ could shift with temperature, reorganization energy, or phonon frequency; mapping where the exponential regime begins would give a testable boundary between practically diffusive and asymptotically diffusive behavior.","The predicted direction-dependent polaron densities in anisotropic 2D lattices could be tested by preparing polarons at specific sites and imaging the spread, since the two hopping arrangements give sharply different shapes (a localized stripe versus an ellipse)."],"forward_implications":["Nonequilibrium polaron transport in the dispersive Holstein model remains subdiffusive for times that grow exponentially with system size, so a finite-lattice simulation underestimates the difficulty of reaching diffusion.","Equilibrium linear-response diffusion constants describe the asymptotic limit, not the nonequilibrium relaxation: the system can appear subdiffusive on every experimentally practical timescale even though its equilibrium transport is diffusive.","In the dilute limit, 1D relaxation quantitatively fixes the dimensionality-renormalized $\\frac{1}{2N_d}\\frac{d\\mathrm{MSD}}{dt}$ and $\\alpha$ in 2D homogeneous lattices, so 2D transport can be built from a 1D outer product without loss of accuracy in the transport metrics.","Anisotropic 2D energy landscapes direct polaron flow: lowering one hopping pathway confines the polaron along one axis, while raising the cross-hopping restores an elliptical distribution—predictions accessible to spatiotemporal microscopy.","The STL-GME construction extends beyond polarons to other short-range lattice models, including spin relaxation, qubit crosstalk, topological magnons, and thermal conduction."],"supporting_citations":[{"why":"supplies the space-local memory construction that lets small-lattice dynamics generate thermodynamic-limit dynamics","marker":"[43]"},{"why":"defines the time-local GME generator whose finite lifetime enables arbitrary-time propagation","marker":"[54]"},{"why":"provides the HEOM solver used to generate the numerically exact small-lattice reference dynamics","marker":"[52]"},{"why":"supplies the n-particle approximation and convergence checks that make the HEOM reference affordable","marker":"[53]"},{"why":"prior observation of long-lived subdiffusion on ~40-site lattices, the finite-size-limited result this paper extends","marker":"[50]"},{"why":"GME treatment of Holstein polaron transport that the paper builds on for long-time dynamics","marker":"[42]"},{"why":"equilibrium linear-response calculation expecting diffusive transport, which the nonequilibrium result contrasts with","marker":"[58]"},{"why":"equilibrium linear-response diffusion constant for the same model, the asymptotic limit approached only slowly after nonequilibrium excitation","marker":"[59]"}],"fun_headline_variants":["Polaron diffusion onset delayed exponentially with lattice size","Exact simulations show polaron diffusion arrives only exponentially late","Nonequilibrium relaxation exponentially delays polaron diffusion onset","Polarons diffuse only after exponential slow-down in large lattices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption that carries the whole extrapolation is that the memory kernel's cutoffs—$d_U=4r_0$ in space and $\\tau_U=90$ fs in time, fixed from a 10-site, sub-picosecond HEOM reference—stay quantitatively exact when the generator is zero-padded to 1500 sites and stepped for 1.5 ns.","fun_headline_variants_meta":{"raw":{"variants":["Polaron diffusion onset delayed exponentially with lattice size","Exact simulations show polaron diffusion arrives only exponentially late","Nonequilibrium relaxation exponentially delays polaron diffusion onset","Polarons diffuse only after exponential slow-down in large lattices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000468,"raw_usage":{"total_tokens":2389,"prompt_tokens":1060,"completion_tokens":1329,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":1262}},"tokens_in":676,"tokens_out":1329,"duration_ms":11383,"temperature":1.0,"reasoning_tokens":1262,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:38:54.949615+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly compute the memory kernel $U(r,t)$ from a larger exact reference, say 30 sites run well past $\\tau_U=90$ fs, and check whether transitions beyond $r=4r_0$ or times beyond $\\tau_U$ are genuinely negligible; alternatively, compare the STL-GME populations on sites farther than $d_U$ from the injection site against exact dynamics on the same lattice. If non-negligible tails appear or the error grows with propagation time, the exponential $1-\\alpha_{\\max}$ scaling is an artifact of the truncation.","supporting_citations":[{"cited_title":"Space-local memory in generalized master equations: Reaching the thermodynamic limit for the cost of a small lattice simulation","cited_arxiv_id":"2411.08598","evidence_quote":"supplies the space-local memory construction that lets small-lattice dynamics generate thermodynamic-limit dynamics"},{"cited_title":"Sayer and A","cited_arxiv_id":null,"evidence_quote":"defines the time-local GME generator whose finite lifetime enables arbitrary-time propagation"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the HEOM solver used to generate the numerically exact small-lattice reference dynamics"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the n-particle approximation and convergence checks that make the HEOM reference affordable"},{"cited_title":"Bhattacharyya, T","cited_arxiv_id":null,"evidence_quote":"prior observation of long-lived subdiffusion on ~40-site lattices, the finite-size-limited result this paper extends"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"GME treatment of Holstein polaron transport that the paper builds on for long-time dynamics"},{"cited_title":"Takahashi and R","cited_arxiv_id":null,"evidence_quote":"equilibrium linear-response calculation expecting diffusive transport, which the nonequilibrium result contrasts with"},{"cited_title":"Bhattacharyya, T","cited_arxiv_id":null,"evidence_quote":"equilibrium linear-response diffusion constant for the same model, the asymptotic limit approached only slowly after nonequilibrium excitation"}],"review_version":1}