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.
Adiabatic theorem and non-adiabatic corrections have been widely applied in modern quantum technology. Recently, non-adiabatic linear response theory has been developed to probe the many-body correlations in closed systems. In this work, we generalize the non-adiabatic linear response theory to open quantum many-body systems. We show that, similar to the closed case, the first-order deviation from the instantaneous steady state is memoryless—it depends only on the final parameters and not on the initial state or ramping path. When ramping the Hamiltonian, the linear response of observables is governed by the derivative of the retarded Green’s function, as in closed systems. In contrast, ramping the dissipation gives rise to a different response, characterized by a high-order correlation function defined in the steady state. Our results offer a compact and physically transparent formulation of non-adiabatic response in open systems, and demonstrate that ramping dynamics can serve as a versatile tool for probing many-body correlations beyond equilibrium.
A linear response theory for open quantum systems is formulated by means of the time-local and time-nonlocal quantum master equations, where a relevant quantum system interacts with a thermal reservoir as well as with an external classical field. A linear response function that characterizes how a relaxation process deviates from its intrinsic process by a weak external field is obtained by extracting the linear terms with respect to the external field from the quantum master equation. It consists of four parts. One represents the linear response of a quantum system when system-reservoir correlation at an initial time and correlation between reservoir states at different times are neglected. The others are correction terms due to these effects. The linear response function is compared with the Kubo formula in the usual linear response theory. To investigate the properties of the linear response of an open quantum system, an exactly solvable model for a stochastic dephasing of a two-level system is examined. Furthermore, the method for deriving the linear response function is applied for calculating two-time correlation functions of open quantum systems. It is shown that the quantum regression theorem is not valid for open quantum systems unless their reduced time evolution is Markovian.