Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2604.13541 · The Journal of Chemical Physics

Beyond the Quantum Regression Theorem in Variational Polaron Master Equations: from Ohmic to Super-Ohmic Environments

Matias Bundgaard-Nielsen, Jake Iles-Smith

quant-phSubmitted Apr 15, 2026• Announced Apr 16, 2026• 0 citations• Revised Aug 21, 2026 (v2)
v1 · Apr 15, 2026v2 · Aug 21, 2026
arXivHTMLPDF

Abstract

While the quantum regression theorem (QRT) is the standard tool for computing multi-time correlation functions in open quantum systems, it relies on system-bath separability and an environment that remains in equilibrium, assumptions that are violated once dynamical correlations develop. Using the projection operator formalism, we derive an extension to the QRT that explicitly incorporates these correlation-induced corrections, enabling us to capture multi-time correlation functions in non-Markovian regimes. We apply this framework to the variational polaron master equation for the spin-boson model in ohmic and super-ohmic regimes, where the polaron transformation mixes system-bath degrees of freedom to produce a non-thermal effective environment. Benchmarking against numerically exact tensor-network simulations, we demonstrate quantitative agreement for single- and two-time observables, including linear-response spectra, even at strong coupling. Our approach broadens the reach of analytic master equations to multi-time observables in strong-coupling and non-Markovian regimes, where environmental memory effects and system-bath correlations are crucial.

Journal reference
J. Chem. Phys. 165, 074111 (2026)
DOI
10.1063/5.0338834
Comments
29 Pages, 9 Figures

Related papers

Capturing non-Markovian polaron dressing with the master equation formalism

Jake Iles-Smith, Owen Diba, Ahsan Nazir

Understanding the dynamics of open quantum systems in strong coupling and non-Markovian regimes remains a formidable theoretical challenge. One popular and well-established method of approximation in these circumstances is provided by the polaron master equation (PME). In this work, we re-evaluate and extend the validity of the PME to capture the impact of non-Markovian polaron dressing, induced by non-equilibrium open system dynamics. By comparing with numerically exact techniques, we confirm that while the standard PME successfully predicts the dynamics of system observables that commute with the polaron transformation (e.g., populations in the Pauli z-basis), it can struggle to fully capture those that do not (e.g., coherences). This limitation stems from the mixing of system and environment degrees of freedom inherent to the polaron transformation, which affects the accuracy of calculated expectation values within the polaron frame. Employing the Nakajima-Zwanzig projection operator formalism, we introduce correction terms that provide an accurate description of observables that do not commute with the transformation. We demonstrate the significance of the correction terms in two cases, the canonical spin-boson model and a dissipative time-dependent Landau-Zener protocol, where they are shown to impact the system dynamics on both short and long timescales.

arXiv:2407.10744•The Journal of Chemical Physics•Oct 3, 2024•13 citations

An efficient variational polaron master equation for non-Markovian spin-boson dynamics: Transformed initial states and general observables

Chien-Chang Chen, Hsi-Sheng Goan

We investigate the effects of initial states on the dynamics of a quantum spin-boson model using an efficient variational polaron master equation approach. The transformation of a factorized system-environment initial state into the variational polaron frame introduces inhomogeneous terms in the master equation, essential for accurately describing non-Markovian dynamics beyond the weak coupling regime. To efficiently account for these effects, we formulate the dynamics in an extended auxiliary Liouville space, converting the resulting integro-differential equations into coupled ordinary differential equations. We further introduce a simplified variational polaron technique for computing the time evolution of expectation values of general observables, including those that do not commute with the polaron transformation, overcoming a major limitation of previous approaches and enabling the unified treatment of coherence and population dynamics. The effects of transformed initial states generated from arbitrary factorized initial states are systematically characterized through expectation values and trace-distance analysis, providing insight into the associated bath-relaxation dynamics over a broad range of tunneling amplitudes, temperatures, and coupling strengths. Comparison with the weak-coupling master equation reveals a delay phenomenon associated with multi-bath-quanta processes, pronounced at intermediate coupling strengths and tunneling amplitudes. We further show that coherent superposition initial states induce stronger nonequilibrium effects than localized states, with inhomogeneous terms playing a crucial role in the short-time dynamics. The proposed formalism applies to both undriven and resonantly driven spin-boson models, providing a computationally efficient framework for open quantum-system dynamics with arbitrary factorized initial states and general physical observables.

quant-pharXiv:2610.08953•Oct 8, 2026•0 citations

Multitime memory beyond the quantum regression theorem in sequential measurement statistics

Paolo Luppi, Claudia Benedetti, Andrea Smirne

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.

quant-pharXiv:2605.06427•May 8, 2026•0 citations

Quantum Regression Theory and Efficient Computation of Response Functions for Non-Markovian Open Systems

Xiantao Li, Chunhao Wang

Linear response functions are a cornerstone concept in physics as they enable efficient estimation of many dynamical properties. In addition to predicting dynamics of observables under perturbations without resimulating the system, these response functions lead to electric conductivity, magnetic susceptibility, dielectric constants, etc. Estimating two-time correlation functions is a key ingredient of measuring linear response functions. However, for open quantum systems, simulating the reduced density operator with a quantum master equation only yields \emph{one-point} observables and is insufficient for this task. In this paper, we develop a memoryless, system-only formulation of two-point correlations for open quantum systems that extends the standard quantum regression theorem (QRT) beyond the Markov limit. We further incorporate the spectral property of the bath and express the time propagators in the response function as the memoryless generators in Lindblad-type forms. The resulting expressions recast the total response function into evolutions generated by time-dependent Hamiltonian and Lindblad primitives together with the more challenging propagation of commutators and anti-commutators. In addition to the derivation of the new QRT, we present quantum algorithms for these primitives and obtain an estimator for two-time correlations whose cost scales poly-logarithmically in the system dimension and $1/ε^{1.25}$ in the target accuracy $ε$. The framework removes the separability (Born-Markov) assumption and offers a pathway to efficient computation of nonequilibrium properties from open quantum systems.

arXiv:2510.05472•Oct 7, 2025•0 citations

Polaron Transformed Canonically Consistent Quantum Master Equation

Juzar Thingna, Xiansong Xu, Daniel Manzano

A central challenge in the theory of open quantum systems is the development of theoretical frameworks capable of accurately describing large, strongly interacting quantum many-body systems in the regime of strong system-bath interactions. In this work, we take a step toward this goal by formulating a polaron-transformed version of the canonically consistent quantum master equation (CCQME) [T. Becker \emph{et~al.}, Phys. Rev. Lett. \textbf{129}, 200403 (2022)]. The CCQME extends beyond standard weak coupling approaches while retaining the same numerical complexity as conventional quantum master equations, thereby enabling the treatment of large quantum systems. The polaron transformation further enhances the accessible system-bath interaction strengths, allowing us to move from moderate to ultra-strong interaction regimes. We present a unified and transparent derivation of these two approaches and combine them to obtain the polar-transformed CCQME (PT-CCQME). Applying our method to the paradigmatic spin boson model, we find excellent agreement with numerically exact time-evolving matrix product operator (TEMPO) simulations. Finally, we predict an initial-state-independent \emph{slowing} down of thermalization in the strong-coupling regime of the spin-boson model.

quant-pharXiv:2604.02731•The Journal of Chemical Physics•Apr 6, 2026•1 citation

Emergence of distinct relaxation behaviour and Quantum Regression Theorem in the Ultra-strong Coupling Limit

Sakil Khan, Bijay Kumar Agarwalla

In the framework of open quantum systems, we derive the dynamical equation governing two-time correlation functions in the ultra-strong coupling (USC) regime between the system and its environment. Unlike the case of the standard weak-coupling regime, in the USC case, we find distinct relaxation behavior for two-time correlators depending on the types of the operators involved in the correlation function. Interestingly, the Quantum Regression Theorem (QRT) emerges after the fastest relaxation time-scale, which is governed by the system-bath coupling strength. We exemplify our findings for the dissipative spin-boson model and further find excellent agreement with the numerically exact hierarchical equations of motion (HEOM) method.

arXiv:2508.08378•Aug 11, 2025•0 citations
While the quantum regression theorem (QRT) is the standard tool for computing multi-time correlation functions in open quantum systems, it relies on system-bath separability and an environment that remains in equilibrium, assumptions that are violated once dynamical correlations develop. Using the projection operator formalism, we derive an extension to the QRT that explicitly incorporates these correlation-induced corrections, enabling us to capture multi-time correlation functions in non-Markovian regimes. We apply this framework to the variational polaron master equation for the spin-boson model in ohmic and super-ohmic regimes, where the polaron transformation mixes system-bath degrees of freedom to produce a non-thermal effective environment. Benchmarking against numerically exact tensor-network simulations, we demonstrate quantitative agreement for single- and two-time observables, including linear-response spectra, even at strong coupling. Our approach broadens the reach of analytic master equations to multi-time observables in strong-coupling and non-Markovian regimes, where environmental memory effects and system-bath correlations are crucial.
Understanding the dynamics of open quantum systems in strong coupling and non-Markovian regimes remains a formidable theoretical challenge. One popular and well-established method of approximation in these circumstances is provided by the polaron master equation (PME). In this work, we re-evaluate and extend the validity of the PME to capture the impact of non-Markovian polaron dressing, induced by non-equilibrium open system dynamics. By comparing with numerically exact techniques, we confirm that while the standard PME successfully predicts the dynamics of system observables that commute with the polaron transformation (e.g., populations in the Pauli z-basis), it can struggle to fully capture those that do not (e.g., coherences). This limitation stems from the mixing of system and environment degrees of freedom inherent to the polaron transformation, which affects the accuracy of calculated expectation values within the polaron frame. Employing the Nakajima-Zwanzig projection operator formalism, we introduce correction terms that provide an accurate description of observables that do not commute with the transformation. We demonstrate the significance of the correction terms in two cases, the canonical spin-boson model and a dissipative time-dependent Landau-Zener protocol, where they are shown to impact the system dynamics on both short and long timescales.
We investigate the effects of initial states on the dynamics of a quantum spin-boson model using an efficient variational polaron master equation approach. The transformation of a factorized system-environment initial state into the variational polaron frame introduces inhomogeneous terms in the master equation, essential for accurately describing non-Markovian dynamics beyond the weak coupling regime. To efficiently account for these effects, we formulate the dynamics in an extended auxiliary Liouville space, converting the resulting integro-differential equations into coupled ordinary differential equations. We further introduce a simplified variational polaron technique for computing the time evolution of expectation values of general observables, including those that do not commute with the polaron transformation, overcoming a major limitation of previous approaches and enabling the unified treatment of coherence and population dynamics. The effects of transformed initial states generated from arbitrary factorized initial states are systematically characterized through expectation values and trace-distance analysis, providing insight into the associated bath-relaxation dynamics over a broad range of tunneling amplitudes, temperatures, and coupling strengths. Comparison with the weak-coupling master equation reveals a delay phenomenon associated with multi-bath-quanta processes, pronounced at intermediate coupling strengths and tunneling amplitudes. We further show that coherent superposition initial states induce stronger nonequilibrium effects than localized states, with inhomogeneous terms playing a crucial role in the short-time dynamics. The proposed formalism applies to both undriven and resonantly driven spin-boson models, providing a computationally efficient framework for open quantum-system dynamics with arbitrary factorized initial states and general physical observables.
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.
Linear response functions are a cornerstone concept in physics as they enable efficient estimation of many dynamical properties. In addition to predicting dynamics of observables under perturbations without resimulating the system, these response functions lead to electric conductivity, magnetic susceptibility, dielectric constants, etc. Estimating two-time correlation functions is a key ingredient of measuring linear response functions. However, for open quantum systems, simulating the reduced density operator with a quantum master equation only yields \emph{one-point} observables and is insufficient for this task. In this paper, we develop a memoryless, system-only formulation of two-point correlations for open quantum systems that extends the standard quantum regression theorem (QRT) beyond the Markov limit. We further incorporate the spectral property of the bath and express the time propagators in the response function as the memoryless generators in Lindblad-type forms. The resulting expressions recast the total response function into evolutions generated by time-dependent Hamiltonian and Lindblad primitives together with the more challenging propagation of commutators and anti-commutators. In addition to the derivation of the new QRT, we present quantum algorithms for these primitives and obtain an estimator for two-time correlations whose cost scales poly-logarithmically in the system dimension and 1/ε1.251/ε^{1.25}1/ε1.25 in the target accuracy εεε. The framework removes the separability (Born-Markov) assumption and offers a pathway to efficient computation of nonequilibrium properties from open quantum systems.
A central challenge in the theory of open quantum systems is the development of theoretical frameworks capable of accurately describing large, strongly interacting quantum many-body systems in the regime of strong system-bath interactions. In this work, we take a step toward this goal by formulating a polaron-transformed version of the canonically consistent quantum master equation (CCQME) [T. Becker \emph{et~al.}, Phys. Rev. Lett. \textbf{129}, 200403 (2022)]. The CCQME extends beyond standard weak coupling approaches while retaining the same numerical complexity as conventional quantum master equations, thereby enabling the treatment of large quantum systems. The polaron transformation further enhances the accessible system-bath interaction strengths, allowing us to move from moderate to ultra-strong interaction regimes. We present a unified and transparent derivation of these two approaches and combine them to obtain the polar-transformed CCQME (PT-CCQME). Applying our method to the paradigmatic spin boson model, we find excellent agreement with numerically exact time-evolving matrix product operator (TEMPO) simulations. Finally, we predict an initial-state-independent \emph{slowing} down of thermalization in the strong-coupling regime of the spin-boson model.
In the framework of open quantum systems, we derive the dynamical equation governing two-time correlation functions in the ultra-strong coupling (USC) regime between the system and its environment. Unlike the case of the standard weak-coupling regime, in the USC case, we find distinct relaxation behavior for two-time correlators depending on the types of the operators involved in the correlation function. Interestingly, the Quantum Regression Theorem (QRT) emerges after the fastest relaxation time-scale, which is governed by the system-bath coupling strength. We exemplify our findings for the dissipative spin-boson model and further find excellent agreement with the numerically exact hierarchical equations of motion (HEOM) method.