Initial Correlations in Open Quantum Systems: Constructing Linear Dynamical Maps
and Master Equations
Alessandra Colla,1Niklas Neubrand,1and Heinz-Peter Breuer1, 2
1Institute of Physics, University of Freiburg, Hermann-Herder-Straße 3, D-79104 Freiburg, Germany
2EUCOR Centre for Quantum Science and Quantum Computing,
University of Freiburg, Hermann-Herder-Straße 3, D-79104 Freiburg, Germany
(Dated: October 25, 2022)
We investigate the dynamics of open quantum systems which are initially correlated with their
environment. The strategy of our approach is to analyze how given, fixed initial correlations modify
the evolution of the open system with respect to the corresponding uncorrelated dynamical behavior
with the same fixed initial environmental state, described by a completely positive dynamical map.
We show that, for any predetermined initial correlations, one can introduce a linear dynamical map
on the space of operators of the open system which acts like the proper dynamical map on the set
of physical states and represents its unique linear extension. Furthermore, we demonstrate that this
construction leads to a linear, time-local quantum master equation with generalized Lindblad struc-
ture involving time-dependent, possibly negative transition rates. Thus, the general non-Markovian
dynamics of an open quantum system can be described by means of a time-local master equation
even in the case of arbitrary, fixed initial system-environment correlations. We present some illus-
trative examples and explain the relation of our approach to several other approaches proposed in
the literature.
I. INTRODUCTION
The theory of open quantum systems has since its con-
ception gained more and more traction as a powerful tool
for studying quantum systems [1]. Its formulation for the
case of systems weakly coupled to Markovian baths has
been extended to arbitrary coupling and non-Markovian
behavior with numerous techniques, which brought to
the formulation of popular exact master equations such
as the Nakajima-Zwanzig equation [2, 3] and the time-
convolutionless (TCL) master equation [4, 5]. An often
made assumption in many approaches is that of factor-
izing initial condition between the system and the en-
vironment. Such a requirement has been highly criti-
cized, either on the grounds of it being unphysical [6–8],
or just too restrictive, and the question of initial corre-
lation has become of importance [9–11]. For example,
the treatment of the dynamics of open systems including
correlations in the initial state leads to a general method
for the local detection of correlations between the sys-
tem and an inaccessible environment [12, 13] which has
been realized experimentally both in trapped ion and in
photonic systems [14–16]. Several efforts have thus been
made to extend the theory to allow for correlated initial
states, both in terms of examining conditions for which
the resulting dynamical map is completely positive [17–
23] and of finding alternative methods for the dynami-
cal description of the reduced system [24–28]. The ap-
proaches are numerous and varied, and the wide array of
papers on the subject paints the picture of a complicated
subject matter. Considering initial correlations can mean
very different things, depending on what one assumes of
the initial system-environment state: should the envi-
ronment be in a fixed, specific initial state, or should it
depend on the correlations? Are the system and the en-
vironment initially in a separable, classically correlated,
entangled state? Different choices lead to different re-
duced dynamics, and may induce limitations such as loss
of complete positivity or a restriction of the set of initial
system states that can be studied. Trying to avoid these
choices to allow for any initial total states might instead
result in other losses, like that of a unique dynamical map
[26]. Ultimately, the decision should be made based on
which questions one is trying to answer or which advan-
tage one wants to gain, carefully taking into account the
consequent drawbacks.
A possible limitation of initial correlations between the
system and the environment is that they may impose
non-linearity of the equations describing the evolution of
the reduced system, the dynamical map and the master
equation, depending on the initial conditions chosen – for
example, such non-linearity may appear in the shape of
additional inhomogeneous terms [29, 30]. However, dy-
namical maps in the case of entangled initial state have
been extended to linear maps on matrices by perform-
ing a basis transformation [31, 32], showing that such
a linearized map is, as a consequence, no longer com-
pletely positive. Some groups, primarily Dominy, Sha-
bani and Lidar (DSL) [33, 34], have then more rigorously
approached the subject of initial correlations by develop-
ing general frameworks for which the linearity of the dy-
namical map can be preserved or obtained. The linearity
property, in fact, comes with several advantages, such as
the operator-sum representation due to Choi [35] of linear
and Hermiticity preserving superoperators and the subse-
quent criterion for such a superoperator to be completely
positive, the renowned Kraus representation [36]. More-
over, a dynamical map which is linear – even if not com-
pletely positive (CP) – automatically leads, through its
inverse, to an exact generator of the evolution which can
be put into a time-local generalized Lindblad form with
time-dependent, possibly negative rates [37, 38]. Master
arXiv:2210.13241v1 [quant-ph] 24 Oct 2022