{"id":"25de8565-5b0a-498a-992f-41adf01ac68e","arxiv_id":"2505.21017","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A time-local extrapolation scheme based on stationary one-step maps converges at least as fast as the transfer tensor method in spin-boson model tests.","lead":"Researchers compared two ways to predict the long-time behavior of quantum systems with memory, and found that a simpler time-local method matches or beats the standard time-nonlocal method. This matters because long-time simulations of such systems are expensive, and a cheaper route could speed up the design of quantum devices.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The abstract's 'invariably' claim is contradicted by the paper's own sub-ohmic data: at τc ≈ 12 time-local extrapolation is erratic while TTM gives a finite error; the claim holds only for cutoffs beyond the last singularity of E_{t,0}.","rationale":"The reader's verdict of CONDITIONAL is appropriate. My analysis identifies a closely related but distinct load-bearing concern: the abstract's 'invariably' is not merely an overstatement—it is contradicted by the paper's own Fig. 1(b), where time-local extrapolation fails for a range of cutoff times (τc ≈ 12) while TTM remains finite. The reader's weakest assumption focused on the theoretical justification of stationarity (Markovian embedding, diagonalizability, linear independence); my concern is more directly empirical and internal to the paper's evidence. The two are connected: the failures at τc ≈ 12 are caused by singularities of E_{t,0}, which also threaten the invertibility in Eq. (4) and thus the stationarity construction. The paper's own Discussion acknowledges that time-local extrapolation 'can be recommended unreservedly' only for cutoffs beyond which the dynamical maps are nearly stationary, but the abstract does not carry this qualification. Therefore, the paper needs a revised, qualified claim; the core numerical demonstration that time-local extrapolation can outperform TTM after the memory time remains credible. A conditional acceptance with mandatory abstract revision is the right outcome—exactly what the reader recommended. No change to the reader's verdict is needed, hence UNCHANGED.","tokens_in":9958,"tokens_out":3281,"duration_ms":35582,"concrete_test":"For the sub-ohmic model of Sec. III.A, scan τc in fine increments (e.g., 0.2) over the interval [8, 30] and compute both the time-local and TTM extrapolation errors at t = 80. Then determine the set (or count) of τc values for which |⟨σz⟩_TL(τc) − exact| > |⟨σz⟩_TTM(τc) − exact|. If this set is non-empty (as suggested by Fig. 1(b) near τc ≈ 12), the 'invariably' statement in the abstract is false as written; the claim must be restricted to τc beyond the last singularity of E_{t,0}, or explicitly redefined as an asymptotic statement about convergence rates.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim, as stated in the abstract, is that 'our simple time-local extrapolation invariably converges at least as fast as time-nonlocal extrapolation.' This is contradicted by the paper's own results in Sec. III.A and Fig. 1(b). For the sub-ohmic spin-boson model, at a cutoff time around τc ≈ 12, the time-local extrapolation is described as 'erratic and unreliable,' and the extrapolation error diverges or becomes very large, whereas the TTM still produces a finite approximation with a moderate error. Thus, for finite cutoff times, it is not true that time-local extrapolation invariably converges at least as fast as TTM; there is a finite interval of cutoff times where the time-local scheme is substantially worse. The only way to rescue the 'invariably' wording is to reinterpret it asymptotically as τc → ∞, or to restrict the claim to cutoff times beyond the last singularity of the dynamical map E_{t,0}. Neither qualification appears in the abstract or in Sec. III.A. The stationarity derivation in Sec. II.B also relies on unverified assumptions (existence of a Markovian embedding, diagonalizable Liouvillian, linear independence of M ≤ D² projected dual vectors), and the singularities observed in E_{t,0} directly violate the invertibility required in Eq. (4). Hence the paper's own evidence shows that the practical claim is conditional, and the unconditional 'invariably' is unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper compares two schemes for extrapolating short-time non-Markovian quantum dynamics to long times: the time-nonlocal transfer tensor method (TTM) and a time-local scheme based on the observation that time-dependent time-local dynamical maps become stationary well before the system reaches steady state. The authors derive the time-local extrapolation from an assumed Markovian embedding of the system, and then test both methods on driven spin-boson models with sub-ohmic, ohmic, and super-ohmic spectral densities, using the ACE process-tensor code as a numerically exact benchmark. The central claim, stated in the abstract and repeated in the introduction, is that the time-local extrapolation 'invariably converges at least as fast as' TTM, and that time-nonlocality is not a prerequisite for accurate long-time extrapolation. The numerical comparisons use identical short-time data for both methods and validate against an independent solver, so the comparisons are fair; however, the extent to which the stated claim is supported by the results is limited by the sub-ohmic example, where time-local extrapolation is erratic for cutoff times near a singularity of the dynamical map.","tokens_in":10226,"tokens_out":2699,"duration_ms":31799,"significance":"If the findings are robust, they have clear practical value: practitioners using TTM-style extrapolation could instead apply a stationary time-local map once the map has settled, avoiding the need to construct and store a sequence of transfer tensors. The paper is commendable for its clean numerical methodology: both extrapolations start from the same short-time dynamical maps, no free parameters are fitted, and the long-time benchmark is an independent numerically exact calculation. The claim that time-local extrapolation can outperform TTM in certain regimes is well supported by the examples. However, the manuscript's headline claim of 'invariably converges at least as fast' is not supported by the presented evidence, and the theoretical justification for stationarity rests on a number of unverified assumptions. The paper therefore contains a useful observation and a practical proposal, but the central claim requires substantial qualification before publication.","major_comments":[{"comment":"The abstract's claim that time-local extrapolation 'invariably converges at least as fast as time-nonlocal extrapolation' is contradicted by the paper's own sub-ohmic results. In Fig. 1(b), for cutoff times around τc ≈ 12, the time-local extrapolation is 'erratic and unreliable' and its error is orders of magnitude larger than that of TTM; TTM gives a finite, moderate error while the time-local error becomes very large or divergent. The unconditional claim is therefore false for finite cutoff times. The paper should either replace 'invariably' with a qualified statement, e.g. for cutoff times beyond the last singularity of E_{t,0}, or explicitly state the asymptotic interpretation τc → ∞. This is a load-bearing issue because the abstract and Sec. IV present the unqualified statement as the main conclusion.","section":"Abstract and Sec. III.A, Fig. 1(b)"},{"comment":"The stationarity argument relies on several assumptions that are not established for the systems treated: the existence of a time-independent Markovian embedding, diagonalizability of the extended Liouvillian, and linear independence of the M projected dual vectors with M ≤ D². Footnote 34 correctly allows for exceptional points via generalized eigenvectors, but the more serious problem is that Eq. (4) requires E_{t_n,t_0} to be non-singular, and Fig. 1(c) shows that this condition is violated at isolated intermediate times in the sub-ohmic example. The paper notes this and advises checking for singularities, but the theoretical derivation in Sec. II.B does not incorporate this limitation. The manuscript should state explicitly that Eq. (4) and the stationarity argument apply only when E_{t_n,t_0} is invertible, and that the derivation does not prove stationarity before the last singularity.","section":"Sec. II.B, Eqs. (4)-(6)"},{"comment":"The paper's own recommended operational criterion is that time-local extrapolation should be applied only when the dynamical map is nearly stationary, monitored by |E_{t+Δt,t} - E_{t,t-Δt}|, and that cutoff times near singularities should be avoided. This is a sensible practical guideline, but it is not the same as the abstract's 'invariably' claim. The manuscript should reconcile the practical recommendation with the abstract wording. In particular, the statement in Sec. IV that time-local extrapolation 'can be recommended unreservedly for cutoff times beyond which the dynamical maps are nearly stationary' is conditional and should be reflected in the abstract and introduction.","section":"Sec. III.A, Fig. 1(d) and Sec. IV"}],"minor_comments":[{"comment":"There are several typos: 'prerequiste' should be 'prerequisite', 'simulation open qunatum systems dynamics' should be 'simulate open quantum system dynamics', and 'indepedent' should be 'independent'. The introduction also uses 'typically' where the abstract uses 'invariably'; the wording should be made consistent throughout.","section":"Abstract and Sec. I"},{"comment":"The caption for Fig. 1(b) states 'Difference between extrapolated (from cutoff time τc) and exact value of ⟨σz⟩ at time t = 80.' This is clear, but the vertical axis label in the figure is given as '|extr.−exact|' without stating the observable; please clarify in the caption that the error is for ⟨σz⟩.","section":"Sec. III.A, Fig. 1 caption"},{"comment":"The caption says 'The extrapolation error in (b) is evaluated with respect to reference time t = 100 ps^{-1}.' Since the horizontal axis of panel (b) is the cutoff time and the reference time is a time, the unit should be ps, not ps^{-1}. Please correct the unit.","section":"Sec. III.B, Fig. 3 caption"},{"comment":"The expression for J_QD(ω) has an unbalanced parenthesis: it reads ω^3 ( c_ee^{-ω^2/ω_e^2} - c_he^{-ω^2/ω_h^2} ) but the opening parenthesis appears before c_e and the closing after the second exponential, which is consistent, but the notation is ambiguous because the ω^3 is outside the parenthesis yet no multiplication sign is shown. Please add parentheses or multiplication dots for clarity.","section":"Sec. III.B, Eq. (13)"},{"comment":"The definition of the projected stationary map in Eq. (6) omits an explicit projection onto the system subspace on the right-hand side: the right side uses Pv_j and ṽ_j^† P, but the left side E_{t+Δt,t} acts on system density matrices. It would be clearer to write E_{t+Δt,t} = P [ sum_j e^{(-λ_j+iω_j)Δt} v_j ṽ_j^† ] P and to state that the maps are understood as restricted to the system subspace.","section":"Sec. II.B, Eq. (6)"},{"comment":"The discussion proposes constructing transfer tensors from dynamical maps E_{τ_c, τ_c - nΔt} around the cutoff time, but does not give a concrete test or estimate of when this would outperform time-local extrapolation. This is a reasonable outlook, but it would strengthen the paper to include at least one numerical example or a clear argument for why such a scheme could be advantageous.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the quantum-dynamics community, and the numerical results are clean and reproducible in principle, since the ACE code is cited and the parameters are given. The main issue is that the abstract overstates the findings: the sub-ohmic case itself provides a counterexample to the 'invariably' claim, and the stationary-map derivation depends on assumptions that are not verified. These are fixable by rewording the central claim and adding a clear statement of the validity conditions. I do not see evidence of circularity or inappropriate citation behavior; the self-citations to the ACE code are appropriate because the code is used in the numerical experiments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read on 2505.21017: the useful result is the systematic numerical comparison. They take the known TTM and the known time-local map reconstruction, and show on three spin-boson models that time-local extrapolation, once the local map has settled, converges at least as fast as TTM and often orders of magnitude better. The benchmarking is fair: both methods start from the same short-time dynamical maps, the long-time reference comes from an independent exact solver, and no free parameters are fit. That is a genuinely useful message for practitioners, and it deserves a serious referee. The time-local route is simpler than TTM, and this paper gives concrete evidence that it is not worse.\n\nThe soft spots are real but concentrated. The abstract's 'invariably' is too strong, and their own Fig. 1(b) shows it. In the sub-ohmic case at tau_c approx 12, the time-local extrapolation is erratic and divergent while TTM gives a finite error. The paper later explains this via singularities of E_{t,0} and recommends checking for them, but then the claim should be phrased as 'for cutoff times beyond the last singularity' or 'once the local map is nearly stationary.' As written, the unconditional wording in the abstract is not supported by the evidence.\n\nThe theoretical justification for stationarity in Sec. II B is heuristic: it assumes a time-independent Markovian embedding, a diagonalizable Liouvillian, and linear independence of the projected dual vectors. None of these are verified for the models, and the singularities directly violate the invertibility in Eq. (4). That is not fatal to the practical claim, because the numerical evidence stands on its own, but it means the paper's contribution is empirical, not a proof of when the method works.\n\nMinor: the sub-ohmic example is the only strongly non-Markovian case, and the erratic region is described as 'typically' avoidable. The discussion does contain the right caveats, so the overclaim is mainly in the abstract and the word 'invariably'—an easy fix. Self-citation of the ACE code is not a problem; the solver is published and independently usable.\n\nRecommendation: send it to peer review. Ask the authors to soften the abstract and to state explicitly that reliable time-local extrapolation requires checking for singularities and stationarity of the local map. With that revision, this becomes a solid reference for practitioners choosing between TTM and time-local schemes.","headline":"Useful numerical comparison with an overstrong abstract: time-local extrapolation often beats TTM, but the paper's own sub-ohmic data contradict 'invariably.'","tokens_in":10750,"tokens_out":1657,"would_cite":true,"duration_ms":17807,"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":"This paper argues that long-time extrapolation of non-Markovian open quantum systems can be done with a simple, stationary time-local map instead of the memory-laden transfer tensor method, and that the simpler scheme converges at least…","keywords":["open quantum systems","non-Markovian dynamics","transfer tensor method","time-local dynamical maps","long-time extrapolation","spin-boson model","time-convolutionless master equations","PT-MPO"],"falsifier":"Compute the Frobenius norm $\\|\\mathcal{E}_{t+\\Delta t,t} - \\mathcal{E}_{t,t-\\Delta t}\\|/\\Delta t$ for a strong-coupling sub-ohmic spin-boson model with a small cutoff frequency $\\omega_c$; if this norm has not decayed well below its short-time values by the time the observable $\\langle \\sigma_z \\rangle$ stops changing, the claimed stationarity-before-equilibrium window does not exist.","tokens_in":1669,"feed_emoji":"⚛️","tokens_out":5958,"duration_ms":96770,"temperature":0.7,"pith_summary":"Simulating a small quantum system coupled to a structured environment is expensive because the environment remembers the past, so researchers often learn the short-time dynamics and then extrapolate to long times. The transfer tensor method is the standard tool for this, building a time-nonlocal propagator from short-time dynamical maps. This paper claims that a much simpler time-local recipe works at least as well: invert the short-time dynamical maps to obtain a step-to-step map, wait until that map stops changing, and then apply the fixed map repeatedly. The authors test this on the spin-boson model with sub-ohmic, ohmic, and super-ohmic spectral densities and find that the time-local extrapolation converges at least as fast as the transfer tensor method, sometimes orders of magnitude better. If this is right, practitioners can skip the memory-kernel machinery entirely and use one stationary local map.","feed_headline":"A fixed local map predicts long-time non-Markovian quantum dynamics","feed_subtitle":"Short-time data gives a step map that goes constant long before equilibrium, so long-time behavior costs almost nothing extra.","key_machinery":"The central object is the time-local dynamical map $\\mathcal{E}_{t+\\Delta t,t} = \\mathcal{E}_{t+\\Delta t,t_0} \\mathcal{E}_{t,t_0}^{-1}$, obtained by inverting the short-time propagator from the initial time. The argument is that, for a time-independent Markovian embedding of the open system, all but $M \\le D^2$ eigenmode contributions decay or interfere away after a transient, so this map becomes time-independent as $\\mathcal{E}_{t+\\Delta t,t} \\to \\mathcal{E}_s$ while the reduced density matrix is still evolving. The extrapolation then applies this fixed map indefinitely, with the cutoff $\\tau_c$ chosen after the map becomes nearly stationary and after any singularities of $\\mathcal{E}_{t,t_0}$ have passed.","core_discovery":"The central claim is that time-dependent time-local dynamical maps become stationary long before the open quantum system itself reaches its steady state, so long-time dynamics can be extrapolated by repeated application of a constant local map. In the spin-boson model, the time-local extrapolation defined by $\\rho_{t_{n+1}} = \\mathcal{E}_s \\rho_{t_n}$ with $\\mathcal{E}_s$ taken from the short-time propagators converges at least as quickly as the time-nonlocal transfer tensor method. The paper presents this as evidence that time-nonlocality is not a prerequisite for accurate and efficient long-time extrapolation of non-Markovian quantum dynamics.","pith_inferences":["A testable extension would be to apply the same stationarity check to time-dependent driving: if the map settles to a periodic rather than constant form, a Floquet-like extension of the time-local extrapolation could handle driven systems.","The stationarity-before-equilibrium property suggests a general diagnostic for any open-system solver: compute the difference between consecutive time-local maps and use its decay as a convergence criterion, independent of the specific extrapolation method.","If the stationarity window is generic, then any numerically exact short-time method, not just process-tensor solvers, could be upgraded to a long-time solver by a trivial postprocessing step that requires no memory-kernel fitting.","The observed failure near singularities of $\\mathcal{E}_{t,t_0}$ suggests a practical rule of thumb: extrapolation should begin only after the last singularity in the singular-value spectrum, since the time-local map is otherwise ill-defined."],"forward_implications":["For common spin-boson models, short-time propagation up to the memory time suffices to extrapolate accurately to times far beyond the cutoff, with errors comparable to or smaller than those of the transfer tensor method.","The recipe for safe use is concrete: examine the singular values of the short-time dynamical maps, avoid cutoffs near singularities, and monitor the difference $\\|\\mathcal{E}_{t+\\Delta t,t} - \\mathcal{E}_{t,t-\\Delta t}\\|$ to verify that the map has become stationary.","Time-local dynamical maps give direct access to canonical Lindblad rates, and persistent negative stationary rates indicate that the system remains non-Markovian even at long times.","The transfer tensor method remains useful for other tasks such as tomography, classification of open systems, and reconstruction of memory kernels, even if it is not needed for simple long-time extrapolation.","If a time-local map becomes stationary but the system has not equilibrated, the same data can be used to build time-local master equations and to diagnose the flow of information between system and environment."],"supporting_citations":[{"why":"Introduces the transfer tensor method, the time-nonlocal extrapolation baseline that the paper compares against.","marker":"[13]"},{"why":"Shows how to obtain time-local dynamical maps by inverting dynamical maps from the initial time, the basis of Eq. (4).","marker":"[22]"},{"why":"Supplies the prior observation that time-local dynamical maps become stationary on the scale of the memory time.","marker":"[24]"},{"why":"Establishes the process tensor matrix operator formalism used to obtain numerically exact short-time propagators.","marker":"[25]"},{"why":"Provides the spin-boson PT-MPO algorithm used to generate the benchmark dynamical maps.","marker":"[26]"},{"why":"Describes the implementation of the PT-MPO solver used as the black-box exact short-time propagation method.","marker":"[29]"},{"why":"Defines the canonical Lindblad form and its connection to non-Markovianity, used to interpret the extracted rates.","marker":"[23]"}],"fun_headline_variants":["Time-local beats time-nonlocal for non-Markovian long-time","Constant local map predicts long-time non-Markovian","Short-time data yields stationary map for long-time prediction","Non-Markovian dynamics predicted without time-nonlocal memory","A fixed local map suffices for long-time non-Markovian"],"cache_read_input_tokens":12928,"weakest_assumption_plain":"The whole scheme rests on the assumption that, after a short transient, the step-to-step map of the open system stops changing while the system itself is still far from equilibrium, and that the map inversion does not run into a singularity before that happens.","fun_headline_variants_meta":{"raw":{"variants":["Time-local beats time-nonlocal for non-Markovian long-time","Constant local map predicts long-time non-Markovian","Short-time data yields stationary map for long-time prediction","Non-Markovian dynamics predicted without time-nonlocal memory","A fixed local map suffices for long-time non-Markovian"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001292,"raw_usage":{"total_tokens":5252,"prompt_tokens":896,"completion_tokens":4356,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":4281}},"tokens_in":512,"tokens_out":4356,"duration_ms":29731,"temperature":1.0,"reasoning_tokens":4281,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:42:00.819913+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Frobenius norm $\\|\\mathcal{E}_{t+\\Delta t,t} - \\mathcal{E}_{t,t-\\Delta t}\\|/\\Delta t$ for a strong-coupling sub-ohmic spin-boson model with a small cutoff frequency $\\omega_c$; if this norm has not decayed well below its short-time values by the time the observable $\\langle \\sigma_z \\rangle$ stops changing, the claimed stationarity-before-equilibrium window does not exist.","supporting_citations":[{"cited_title":"Cerrillo and J","cited_arxiv_id":null,"evidence_quote":"Introduces the transfer tensor method, the time-nonlocal extrapolation baseline that the paper compares against."},{"cited_title":"Andersson, J","cited_arxiv_id":null,"evidence_quote":"Shows how to obtain time-local dynamical maps by inverting dynamical maps from the initial time, the basis of Eq. (4)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the prior observation that time-local dynamical maps become stationary on the scale of the memory time."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the process tensor matrix operator formalism used to obtain numerically exact short-time propagators."},{"cited_title":"Cygorek and E","cited_arxiv_id":null,"evidence_quote":"Describes the implementation of the PT-MPO solver used as the black-box exact short-time propagation method."}],"review_version":1}