Multitime memory beyond the quantum regression theorem in sequential measurement statistics
Pith reviewed 2026-05-08 11:09 UTC · model grok-4.3
The pith
For factorized initial states the two-time propagator decomposes into a quantum regression theorem term from the reduced dynamical map plus a memory term from persisting system-environment correlations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For factorized initial states the two-time propagator decomposes exactly into a QRT-like contribution determined entirely by the reduced dynamical map together with a memory term that encodes system-environment correlations persisting across the intervention. In the weak-coupling regime this memory term yields an explicit second-order correction expressed via the reduced map and bath correlation functions. An operational distance between exact and QRT-predicted joint probabilities quantifies the violations, revealing that multitime memory is protocol-dependent and can appear at higher temporal order even when two-time statistics remain compatible with QRT predictions.
What carries the argument
The exact decomposition of the two-time propagator into a QRT-like term fixed by the reduced dynamical map and a memory term that carries system-environment correlations across interventions.
If this is right
- Multitime non-Markovianity witnessed by sequential statistics is related to but inequivalent from one-time non-Markovianity extracted from the reduced state.
- QRT violations can remain invisible at two-time order yet become detectable at higher temporal orders.
- The magnitude of these violations depends on environmental parameters such as spectral density shape and temperature as well as on the chosen measurement protocol.
- Reduced-state non-Markovianity quantifiers therefore provide an incomplete diagnosis of memory effects that appear in multitime measurement records.
Where Pith is reading between the lines
- Full characterization of open-system dynamics for repeated interventions requires tracking multitime correlation functions rather than reduced maps alone.
- Protocols that rely on sequential operations, such as quantum sensing or feedback control, may incur larger errors than single-time estimates suggest once the memory term is included.
- The decomposition offers a route to compute perturbative corrections for any weak-coupling model once the reduced map and bath correlations are known.
- Extension to three-time or higher statistics would likely expose additional memory layers not captured by two-time analysis.
Load-bearing premise
The decomposition assumes factorized initial system-environment states and operates inside the standard reduced dynamical map framework.
What would settle it
A direct numerical or experimental computation of joint probabilities for factorized initial states in the spin-boson model that deviates from both the QRT prediction and the derived second-order memory correction would falsify the claimed decomposition.
Figures
read the original abstract
We investigate the presence of memory in the sequential measurement statistics of an open quantum system, as witnessed by the departure from the quantum regression theorem (QRT), that is, the possibility to predict multitime probabilities from the one-time reduced dynamical map. For factorized initial states, we identify an exact decomposition of the two-time propagator into a QRT-like contribution, fully determined by the reduced dynamical map, and a memory term encoding system--environment correlations across the intervention; in the weak-coupling regime, the memory term yields an explicit second-order correction expressed in terms of the reduced map and bath correlation functions. Furthermore, we introduce an operational quantifier of QRT violations based on the distance between exact and QRT-predicted joint probabilities. Benchmarking the framework on a spin--boson model and using a pseudomode embedding as nonperturbative reference, we comprehensively analyze the impact of spectral-density parameters, environmental temperature, and measurement protocols on the non-Markovianity of the multitime statistics. Comparison with a one-time quantifier shows that reduced-state non-Markovianity and multitime memory are related but inequivalent: the latter, as probed through sequential statistics, is intrinsically protocol dependent and can become visible at higher temporal order even when two-time statistics remain compatible with QRT predictions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a framework for detecting memory effects in the sequential measurement statistics of open quantum systems that go beyond the quantum regression theorem (QRT). For factorized initial states, it derives an exact decomposition of the two-time propagator into a QRT-like term determined solely by the reduced dynamical map and a memory term accounting for system-environment correlations induced by the intervention. In the weak-coupling regime, an explicit second-order perturbative expression for the memory correction is provided in terms of the reduced map and bath correlation functions. An operational quantifier based on the distance between exact and QRT-predicted joint probabilities is introduced. The framework is applied to the spin-boson model using pseudomode embedding as a non-perturbative benchmark to study the dependence on spectral density parameters, temperature, and measurement protocols, revealing that multitime memory is protocol-dependent and inequivalent to one-time reduced-state non-Markovianity.
Significance. If the central decomposition holds, this work offers a valuable tool for quantifying and understanding non-Markovian effects in multitime quantum statistics, which is relevant for quantum control, metrology, and information processing where sequential measurements are common. The separation into QRT and memory contributions clarifies when standard reduced dynamics suffice and when correlations across interventions matter. The perturbative result provides a practical approximation, and the inequivalence finding underscores the limitations of one-time measures. Credit is due for the exact decomposition valid for factorized states and the comprehensive numerical exploration of parameter space.
major comments (1)
- [§4] §4 (Numerical benchmarks and inequivalence claims): The conclusion that multitime memory and reduced-state non-Markovianity are inequivalent relies on numerical comparisons treating the pseudomode embedding as the exact non-perturbative reference. No convergence tests with increasing auxiliary mode number, nor cross-validation against HEOM or TCL methods, are reported. This is load-bearing for the inequivalence and protocol-dependence results, as truncation errors in long-memory or high-temperature regimes could produce spurious QRT departures.
minor comments (2)
- [§2] §2 (Framework and notation): The operational quantifier of QRT violations is introduced via a distance measure on joint probabilities; a short discussion of its invariance properties or relation to other distinguishability measures would improve clarity.
- [§3] §3 (Perturbative expansion): The second-order memory term is expressed using bath correlation functions; ensuring that all intermediate steps explicitly reference the factorized-initial-state assumption would prevent potential misreading.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation of our work and for the constructive comment on the numerical benchmarks. We address the point below and will revise the manuscript to incorporate additional validation.
read point-by-point responses
-
Referee: [§4] §4 (Numerical benchmarks and inequivalence claims): The conclusion that multitime memory and reduced-state non-Markovianity are inequivalent relies on numerical comparisons treating the pseudomode embedding as the exact non-perturbative reference. No convergence tests with increasing auxiliary mode number, nor cross-validation against HEOM or TCL methods, are reported. This is load-bearing for the inequivalence and protocol-dependence results, as truncation errors in long-memory or high-temperature regimes could produce spurious QRT departures.
Authors: We acknowledge that the original manuscript did not include explicit convergence tests with respect to the number of auxiliary modes in the pseudomode embedding or cross-validation against HEOM or TCL methods. This is a valid concern for the robustness of the inequivalence and protocol-dependence claims. In the revised version, we will add a dedicated appendix or subsection presenting convergence data for representative parameter regimes (including long-memory spectral densities and elevated temperatures), demonstrating stability of the multitime statistics with increasing auxiliary modes. We will also include a limited comparison with HEOM results for selected cases to cross-validate the pseudomode benchmarks. These additions will confirm that the observed departures from QRT predictions and the inequivalence to one-time non-Markovianity are not numerical artifacts. The main conclusions remain unchanged. revision: yes
Circularity Check
Derivation chain is self-contained; no load-bearing reductions to inputs or self-citations
full rationale
The central result is an exact decomposition of the two-time propagator for factorized initial states into a term fixed by the reduced dynamical map plus a remainder memory term. This partitioning follows directly from the definitions of the propagators and the assumption of initial factorization; it does not constitute a prediction that is forced by construction or that collapses to a fit. The weak-coupling correction is obtained by standard perturbative expansion in terms of bath correlations, again without circularity. Benchmarking against pseudomode numerics is an external validation step rather than part of the derivation itself. No self-citation is invoked to justify uniqueness or to close the argument. The framework therefore remains independent of its own outputs.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Factorized initial states between system and environment
- standard math Existence of a reduced dynamical map governing one-time evolution
Reference graph
Works this paper leans on
-
[1]
H. J. Carmichael,An Open Systems Approach to Quantum Op- tics, Lecture Notes in Physics Monographs, V ol. 18 (Springer- Verlag, Berlin, 1993). 11 Figure 6. Two- vs three-time QRT error. Left panels: two-time average QRT violation ¯ε QRT(t f ) [Eq. (30)]. Right panels: three-time average QRT violation ¯ε(3) QRT(t f ) [Eq. (61)]. Both quantities are shown o...
1993
-
[2]
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
-
[3]
Breuer and F
H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, Oxford, 2002)
2002
-
[4]
C. W. Gardiner and P. Zoller,Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics, 3rd ed., Springer Series in Synergetics (Springer, Berlin, 2004)
2004
-
[5]
Rivas and S
A. Rivas and S. F. Huelga,Open Quantum Systems: An Intro- duction, SpringerBriefs in Physics (Springer Berlin Heidelberg, 2011)
2011
-
[6]
Vacchini,Open Quantum Systems(Springer Cham, 2024)
B. Vacchini,Open Quantum Systems(Springer Cham, 2024)
2024
-
[7]
Lax, Formal theory of quantum fluctuations from a driven state, Phys
M. Lax, Formal theory of quantum fluctuations from a driven state, Phys. Rev.129, 2342 (1963)
1963
-
[8]
Lax, Quantum noise
M. Lax, Quantum noise. xi. multitime correspondence between quantum and classical stochastic processes, Phys. Rev.172, 350 (1968)
1968
-
[9]
D ¨umcke, The low density limit for an n-level system inter- acting with a free bose gas, J
R. D ¨umcke, The low density limit for an n-level system inter- acting with a free bose gas, J. Math. Phys.24, 311 (1983)
1983
-
[10]
Swain, Master equation derivation of quantum regression theorem, J
S. Swain, Master equation derivation of quantum regression theorem, J. Phys. A: Math. and Gen.14, 2577 (1981)
1981
-
[11]
Alonso and I
D. Alonso and I. de Vega, Multiple-time correlation functions for non-Markovian interaction: Beyond the quantum regression theorem, Phys. Rev. Lett.94, 200403 (2005)
2005
-
[12]
Ivanov and H.-P
A. Ivanov and H.-P. Breuer, Extension of the Nakajima– Zwanzig approach to multitime correlation functions of open systems, Phys. Rev. A92, 032113 (2015)
2015
-
[13]
D. P. S. McCutcheon, Optical signatures of non-Markovian be- havior in open quantum systems, Phys. Rev. A93, 022119 (2016)
2016
-
[14]
M. Ban, S. Kitajima, and F. Shibata, Two-time correlation func- tion of an open quantum system in contact with a Gaussian reservoir, Phys. Rev. A97, 052101 (2018)
2018
-
[15]
Bonifacio and A
M. Bonifacio and A. A. Budini, Perturbation theory for oper- ational quantum non-markovianity, Phys. Rev. A102, 022216 (2020)
2020
-
[16]
Lonigro and D
D. Lonigro and D. Chru ´sci´nski, Quantum regression in dephas- ing phenomena, J. Phys. A: Math. Theor.55, 225308 (2022)
2022
-
[17]
Salamon, O
M. Salamon, O. Dudgeon, A. Nazir, and J. Iles-Smith, Marko- vian approach ton-photon correlations beyond the quantum re- gression theorem, Phys. Rev. Lett.136, 080401 (2026)
2026
-
[18]
Lindblad, Non-Markovian quantum stochastic processes and their entropy, Comm
G. Lindblad, Non-Markovian quantum stochastic processes and their entropy, Comm. Math. Phys.65, 281 (1979)
1979
-
[19]
Accardi, A
L. Accardi, A. Frigerio, and J. T. Lewis, Quantum stochastic processes, Publ. Res. Inst. Math. Sci.18, 97 (1982)
1982
-
[20]
Modi, Operational approach to open dynamics and quantify- ing initial correlations, Sci
K. Modi, Operational approach to open dynamics and quantify- ing initial correlations, Sci. Rep.2, 581 (2012)
2012
-
[21]
F. A. Pollock, C. Rodr ´ıguez-Rosario, T. Frauenheim, M. Pater- nostro, and K. Modi, Operational Markov condition for quan- tum processes, Phys. Rev. Lett.120, 040405 (2018)
2018
-
[22]
Taranto, F
P. Taranto, F. A. Pollock, S. Milz, M. Tomamichel, and K. Modi, Quantum Markov order, Phys. Rev. Lett.122, 140401 (2019)
2019
-
[23]
Milz and K
S. Milz and K. Modi, Quantum stochastic processes and quan- tum non-Markovian phenomena, PRX Quantum2, 030201 (2021)
2021
-
[24]
Giarmatzi and F
C. Giarmatzi and F. Costa, Witnessing quantum memory in non-Markovian processes, Quantum5, 440 (2021)
2021
-
[25]
Vacchini, A
B. Vacchini, A. Smirne, E.-M. Laine, J. Piilo, and H.-P. Breuer, Markovianity and non-Markovianity in quantum and classical systems, New J. Phys.13, 093004 (2011)
2011
-
[26]
Guarnieri, A
G. Guarnieri, A. Smirne, and B. Vacchini, Quantum regression theorem and non-Markovianity of quantum dynamics, Phys. Rev. A90, 022110 (2014). 12
2014
-
[27]
S. Milz, M. S. Kim, F. A. Pollock, and K. Modi, Completely positive divisibility does not mean Markovianity, Phys. Rev. Lett.123, 040401 (2019)
2019
-
[28]
L. Li, M. J. W. Hall, and H. M. Wiseman, Concepts of quantum non-Markovianity: A hierarchy, Phys. Rep.759, 1 (2018)
2018
-
[29]
M. M. Wolf, J. Eisert, T. S. Cubitt, and J. I. Cirac, Assess- ing non-Markovian quantum dynamics, Phys. Rev. Lett.101, 150402 (2008)
2008
-
[30]
Rivas, S
´A. Rivas, S. F. Huelga, and M. B. Plenio, Entanglement and non-Markovianity of quantum evolutions, Phys. Rev. Lett.105, 050403 (2010)
2010
-
[31]
Wißmann, H.-P
S. Wißmann, H.-P. Breuer, and B. Vacchini, Generalized trace-distance measure connecting quantum and classical non- markovianity, Phys. Rev. A92, 042108 (2015)
2015
-
[32]
Breuer, E.-M
H.-P. Breuer, E.-M. Laine, J. Piilo, and B. Vacchini, Collo- quium: Non-Markovian dynamics in open quantum systems, Rev. Mod. Phys.88, 021002 (2016)
2016
-
[33]
Chru ´sci´nski, Dynamical maps beyond Markovian regime, Physics Reports992, 1 (2022)
D. Chru ´sci´nski, Dynamical maps beyond Markovian regime, Physics Reports992, 1 (2022)
2022
-
[34]
B ¨acker, V
C. B ¨acker, V . Link, and W. T. Strunz, Verifying quantum mem- ory in the dynamics of spin-boson models, Phys. Rev. Res.8, 023106 (2026)
2026
- [35]
-
[36]
B. M. Garraway, Nonperturbative decay of an atomic system in a cavity, Phys. Rev. A55, 2290 (1997)
1997
-
[37]
Tamascelli, A
D. Tamascelli, A. Smirne, S. F. Huelga, and M. B. Plenio, Nonperturbative treatment of non-Markovian dynamics of open quantum systems, Phys. Rev. Lett.120, 030402 (2018)
2018
-
[38]
Chru ´sci´nski, A
D. Chru ´sci´nski, A. Kossakowski, and A. Rivas, Measures of non-Markovianity: Divisibility versus backflow of information, Phys. Rev. A83, 052128 (2011)
2011
-
[39]
Kubo, The fluctuation-dissipation theorem, Rep
R. Kubo, The fluctuation-dissipation theorem, Rep. Prog. Phys. 29, 255 (1966)
1966
-
[40]
A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf, Introduction to quantum noise, measurement, and amplification, Rev. Mod. Phys.82, 1155 (2010)
2010
-
[41]
Heinosaari and M
T. Heinosaari and M. Ziman,The Mathematical Language of Quantum Theory: From Uncertainty to Entanglement(Cam- bridge University Press, 2011)
2011
-
[42]
Smirne, D
A. Smirne, D. Egloff, M. Garc´ıa D´ıaz, M. B. Plenio, and S. F. Huelga, Coherence and non-classicality of quantum Markov processes, Quantum Sci. Technol.4, 01LT01 (2019)
2019
-
[43]
S. Milz, D. Egloff, P. Taranto, T. Theurer, M. B. Plenio, A. Smirne, and S. F. Huelga, When is a non-Markovian quan- tum process classical?, Phys. Rev. X10, 041049 (2020)
2020
-
[44]
A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Dynamics of the dissipative two-state system, Rev. Mod. Phys.59, 1 (1987)
1987
-
[45]
Breuer, B
H.-P. Breuer, B. Kappler, and F. Petruccione, Stochastic wave- function method for non-Markovian quantum master equations, Phys. Rev. A59, 1633 (1999)
1999
-
[46]
Vacchini and H.-P
B. Vacchini and H.-P. Breuer, Exact master equations for the non-Markovian decay of a qubit, Phys. Rev. A81, 042103 (2010)
2010
-
[47]
Albarelli, B
F. Albarelli, B. Vacchini, and A. Smirne, Pseudomode treat- ment of strong-coupling quantum thermodynamics, Quantum Sci. Technol.10, 015041 (2024)
2024
-
[48]
Gorini, A
V . Gorini, A. Kossakowski, and E. C. G. Sudarshan, Com- pletely positive dynamical semigroups of n-level systems, J. Math. Phys.17, 821 (1976)
1976
-
[49]
Lindblad, On the Generators of Quantum Dynamical Semi- groups, Commun
G. Lindblad, On the Generators of Quantum Dynamical Semi- groups, Commun. Math. Phys.48, 119 (1976)
1976
-
[50]
F. Chen, E. Arrigoni, and M. Galperin, Markovian treatment of non-Markovian dynamics of open fermionic systems, New J. Phys.21, 123035 (2019)
2019
-
[51]
Smirne, D
A. Smirne, D. Tamascelli, J. Lim, M. B. Plenio, and S. F. Huelga, Non-perturbative treatment of open-system multi-time expectation values in Gaussian bosonic environments, Open Syst. Inf. Dyn.29, 2250019 (2022)
2022
-
[52]
J. R. Johansson, P. D. Nation, and F. Nori, QuTiP: An open- source Python framework for the dynamics of open quantum systems, Comput. Phys. Commun.183, 1760 (2012)
2012
-
[53]
J. R. Johansson, P. D. Nation, and F. Nori, QuTiP 2: A Python framework for the dynamics of open quantum systems, Com- put. Phys. Commun.184, 1234 (2013)
2013
-
[54]
Beth Ruskai, S
M. Beth Ruskai, S. Szarek, and E. Werner, An analysis of completely-positive trace-preserving maps on m2, Linear Al- gebra and its Applications347, 159 (2002)
2002
-
[55]
Swinbank and R
R. Swinbank and R. J. Purser, Fibonacci grids: A novel ap- proach to global modelling, Q. J. R. Meteorol. Soc.132, 1769 (2006)
2006
-
[56]
Gonz´alez, Measurement of areas on a sphere using Fibonacci and latitude–longitude lattices, Math
´A. Gonz´alez, Measurement of areas on a sphere using Fibonacci and latitude–longitude lattices, Math. Geosci.42, 49 (2010)
2010
-
[57]
S. Milz, F. A. Pollock, and K. Modi, An introduction to oper- ational quantum dynamics, Open Syst. Inf. Dyn.24, 1740016 (2017)
2017
-
[58]
Strasberg and M
P. Strasberg and M. G. D´ıaz, Classical quantum stochastic pro- cesses, Phys. Rev. A100, 022120 (2019)
2019
-
[59]
Strasberg, Classicality with(out) decoherence: Concepts, re- lation to Markovianity, and a random matrix theory approach, SciPost Phys.15, 024 (2023)
P. Strasberg, Classicality with(out) decoherence: Concepts, re- lation to Markovianity, and a random matrix theory approach, SciPost Phys.15, 024 (2023)
2023
-
[60]
Luppi, C
P. Luppi, C. Benedetti, and A. Smirne, Temporal nonclassicality in continuous-time quantum walks, Phys. Rev. A113, 042212 (2026)
2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.