Get direct links to references, BibTeX extraction and comments on all arXiv papers with: Librarian
Quasi-probability distributions for observ-
ables in dynamic systems
Patrick P. Hofer
Department of Applied Physics, University of Geneva, 1211 Geneva 4, Switzerland
October 12, 2017
We develop a general framework to investigate ﬂuctuations of non-commuting ob-
servables. To this end, we consider the Keldysh quasi-probability distribution (KQPD).
This distribution provides a measurement-independent description of the observables
of interest and their time-evolution. Nevertheless, positive probability distributions for
measurement outcomes can be obtained from the KQPD by taking into account the
eﬀect of measurement back-action and imprecision. Negativity in the KQPD can be
linked to an interference eﬀect and acts as an indicator for non-classical behavior. No-
table examples of the KQPD are the Wigner function and the full counting statistics,
both of which have been used extensively to describe systems in the absence as well
as in the presence of a measurement apparatus. Here we discuss the KQPD and its
moments in detail and connect it to various time-dependent problems including weak
values, ﬂuctuating work, and Leggett-Garg inequalities. Our results are illustrated
using the simple example of two subsequent, non-commuting spin measurements.
1 Introduction
Quasi-probability distributions such as the Wigner function have been an important tool in quan-
tum mechanics since the early days of the theory [1]. Indeed, the Wigner function lies at the heart
of the phase-space representation of quantum mechanics which is tightly connected to classical
statistical mechanics and provides an alternative route to quantization [2]. In classical statistical
mechanics, averages of observables can be described as averages taken with respect to some un-
derlying probability distribution. This probability distribution fully describes the system under
consideration and reﬂects the lack of knowledge of the observer. If one has perfect knowledge
of the behavior of all degrees of freedom, there would be no need for probability distributions
and measurement outcomes could be predicted in a deterministic manner. This is not possible
in quantum mechanics due to the probabilistic nature of the theory. Similarly to classical sta-
tistical mechanics, the ﬂuctuations in quantum mechanics can be encoded in quasi-probability
distributions. These distributions can take on negative values which is ultimately a consequence
of the non-commutativity of operators, reﬂecting the fact that a measurement inevitably disturbs
all subsequent measurements of observables which do not commute with the ﬁrst one. Negative
quasi-probability distributions therefore reﬂect the impossibility of accessing observables without
disturbing them, indicating a non-classical behavior of the system which requires to take into ac-
count the action of the measurement apparatus to predict its outputs. The concept of negative
quasi-probability distributions, or negativity in short, is therefore closely related to the concept
of contextuality which deals with the dependence of measurement outcomes on the experimental
setup [36].
So far, most quasi-probability distributions describe the quantum state of a system. This holds
true for the Wigner function and its generalizations to ﬁnite dimensional Hilbert spaces [68], as well
as a number of other distributions used in quantum optics [9, 10]. Negativity in these distributions
Patrick P. Hofer: patrick.hofer@unige.ch
Accepted in Quantum 2017-09-19, click title to verify 1
arXiv:1702.00998v4 [quant-ph] 11 Oct 2017
indicates non-classicality of the quantum state. However, a state that exhibits negativity in some
distribution might not do so in another one [11], implying that non-classicality is not necessarily
a general feature but depends on the way a state is interrogated (see also Refs. [35, 12]). A well
studied quasi-probability distribution which does not describe the quantum state of a system but
rather its dynamics is the full counting statistics (FCS) [1319]. Negativity in the FCS has recently
been connected to a peculiar interference eﬀect and it was shown that systems which do not show
any negativity in their Wigner function can still show negativity in the FCS [19]. This further
stresses the fact that systems might exhibit completely classical behavior in some aspects, e.g.
their instantaneous states, while exhibiting non-classicality in others, e.g. their dynamics which
determine the temporal ﬂuctuations within the system. It would therefore be highly desirable
to have access to a quasi-probability distribution which encodes only the aspects of interest of a
system. This is in marked contrast to Refs. [35, 20] which consider a full description of states and
measurements in terms of quasi-probabilities.
The goal of the present paper is to introduce a general quasi-probability distribution that can be
tailored to the problem at hand and thus provides a general framework for investigating ﬂuctuations
of non-commuting observables. Because this distribution is based on the Keldysh formalism, we
call it the Keldysh quasi-probability distribution (KQPD). When the observables of interest are
the position and momentum operators, the KQPD reduces to the Wigner function. Similarly, when
one is interested in the time-integral of an observable, the FCS of the corresponding operator is
obtained. Furthermore, the outcome of von Neumann type measurements [21] can be predicted
by convolving the corresponding KQPD with the Wigner function of the detectors [cf. Sec. 3].
This convolution takes into account the imprecision as well as the back-action of the performed
measurement. An example is provided by the Husimi Q-function which is obtained by convolving
the Wigner function with Gaussians and predicts the outcomes of simultaneous but imprecise
measurements of position and momentum [22].
The rest of the article is structured as follows. In Sec. 2, we introduce the KQPD, discuss
the Wigner function, the FCS, and two subsequent spin measurements as examples, and discuss
negativity of the KQPD as an indicator for non-classical behavior. In Sec. 3, we provide a more
physical motivation for investigating the KQPD by connecting it to measurement outcomes which
requires the introduction of measurement imprecision and back-action. Before concluding in Sec. 5,
we discuss applications of the KQPD in Sec. 4, including the distribution of thermodynamic work,
(anomalous) weak values, and Leggett-Garg inequalities. Throughout this paper, we set ~ = 1.
2 The Keldysh quasi-probability distribution
In this section we brieﬂy introduce the Keldysh formalism, a powerful framework for tackling out-
of-equilibrium problems, where the KQPD arises quite naturally through its moment generating
function. Here we only touch upon the ideas in the Keldysh formalism which are relevant for our
discussion. For a detailed introduction to the subject we refer the reader to Ref. [23]. We will
start by considering the somewhat trivial example of the statistics of a single observable before we
generalize the results to the highly non-trivial scenario of multiple, non-commuting observables at
diﬀerent times.
2.1 Single observable
The main idea that we will be concerned with is the evolution along the closed time contour
illustrated in Fig. 1, an idea that goes back to Schwinger [24]. In this section, we are interested in
the statistics of a single observable
ˆ
A at time τ . Its average value can be written as
hAi
τ
= Tr
n
ˆ
U(0, τ)
ˆ
A
ˆ
U(τ, 0)ˆρ
0
o
=
1
2
Tr
n
ˆ
U(0, T )
ˆ
U(T , τ)
ˆ
A
ˆ
U(τ, 0)ˆρ
0
o
+
1
2
Tr
n
ˆ
U(0, τ)
ˆ
A
ˆ
U(τ, T )
ˆ
U(T , 0)ˆρ
0
o
,
(1)
where
ˆ
U(t
0
, t) =
ˆ
U
(t, t
0
) is the time-evolution operator from time t to time t
0
and ˆρ
0
denotes
the density matrix at time t = 0. Throughout this article, we consider unitary time evolution.
Generalizing the results to non-unitary time evolution is in principle straightforward and discussed
Accepted in Quantum 2017-09-19, click title to verify 2
Figure 1: Illustration of the closed time contour. The average h
ˆ
Ai
τ
can be obtained by evolving the state (or
a mixture of states) up to time τ , where the operator
ˆ
A is evaluated, and back again. Alternatively, one can
extend this closed-time contour to the Keldysh contour (from t = 0 to t = T and back again) and evaluate the
operator either on the forward-in-time or on the backward-in-time branch.
in App. A. Writing the density matrix as a mixture of pure states, the ﬁrst line of the last expression
can be interpreted as time evolving the states up to time τ , evaluating the operator
ˆ
A, and evolving
the states backward to the initial time t = 0 [cf. the left-hand side of Fig. 1]. On the second line,
we inserted identities 1 =
ˆ
U(τ, T )
ˆ
U(T , τ), extending the contour along which the state is time
evolved. We call the resulting contour, from t = 0 to t = T and back again, the Keldysh contour.
We can evaluate the operator
ˆ
A either along the forward-in-time or along the backward-in-time
branch of the contour. As discussed below, the most natural choice is to treat both contours on
an equal footing as illustrated graphically in Fig. 1.
ˆ
A can be obtained by considering its moment generating
function given by
Λ(λ; τ) = he
ˆ
A
i
τ
= Tr
n
ˆ
U(0, τ)e
i
λ
2
ˆ
A
ˆ
U(τ, T )
ˆ
U(T , τ)e
i
λ
2
ˆ
A
ˆ
U(τ, 0)ˆρ
0
o
. (2)
All higher moments can be obtained from this generating function by taking derivatives with
respect to λ
h
ˆ
A
k
i
τ
= i
k
k
λ
Λ(λ; τ)
λ=0
. (3)
In anticipation of the discussion below, we wrote Eq. (2) symmetrically in terms of the Keldysh
contour. Taking the derivative with respect to λ at λ = 0, we recover Eq. (1). Since the moment
generating function is the average of the exponential of
ˆ
A, it corresponds to the Fourier transform
of a probability distribution associated to the ﬂuctuations of
ˆ
A
Λ(λ; τ) =
Z
−∞
dAP(A; τ)e
iλA
, P(A; τ) =
X
A
0
δ(A A
0
)hA
0
|
ˆ
U(τ, 0)ˆρ
0
ˆ
U
(τ, 0)|A
0
i, (4)
where the sum over A
0
goes over all eigenstates
ˆ
A|A
0
i = A
0
|A
0
i and has to be replaced by an
integral if
ˆ
A has a continuous spectrum. In the case of a single observable, the KQPD P(A; τ ) is
a positive probability distribution which predicts the outcome of projectively measuring
ˆ
A at time
τ.
The symmetric treatment with respect to the forward-in-time and the backward-in-time branches
of the Keldysh contour might seem like an unnecessary complication at this point. Below, we show
how such a treatment for multiple observables leads to a time-ordering prescription which is rel-
evant for von Neumann type measurements. Before discussing multiple observables in detail, we
express the moment generating function in terms of the classical and quantum ﬁelds obtained by
a Keldysh rotation [cf. Eq. (8)]. To this end, we assume that the observable
ˆ
A has a discrete
spectrum and that its eigenstates span the whole Hilbert space, i.e. 1 =
P
A
|AihA|. We then
divide the time evolution along the Keldysh contour in Eq. (2) into 2N steps (N for the forward-,
N for the backward-in-time branch) of inﬁnitesimal length δt and insert identities in between the
steps. Writing the trace in terms of the eigenstates of
ˆ
A, we then ﬁnd
Λ(λ; τ) =
X
A
+
,A
δ
A
+
T
,A
T
e
i
λ
2
(A
+
τ
+A
τ
)
A(A
+
)A
(A
)hA
+
0
|ˆρ|A
0
i. (5)
Here, A
+()
t
denotes the state that was inserted at time t on the forward-in-time (backward-in-time)
branch, the vectors A
±
gather all A
±
t
and the amplitudes along the branches read
A(A
±
) = hA
±
T
|
ˆ
U(T , T δt)|A
±
T δt
ihA
±
T δt
|···|A
±
δt
ihA
±
δt
|
ˆ
U(δt, 0)|A
±
0
i. (6)
Accepted in Quantum 2017-09-19, click title to verify 3
The Kronecker delta in Eq. (5) is due to the fact that the point t = T connects the two branches
(cf. Fig. 1). We can write the moment generating function in the form of a path integral
Λ(λ; τ) =
Z
D[A
+
t
, A
t
]e
i
λ
2
(A
+
τ
+A
τ
)
e
iS[A
+
t
,A
t
]
, (7)
where the action S is determined by exp(iS[A
+
t
, A
t
]) = A(A
+
)A
(A
)hA
+
0
|ˆρ|A
0
i and the inte-
R
D[A
+
t
, A
t
] =
P
A
+
,A
δ
A
+
T
,A
T
. Finally, we perform a Keldysh rotation
by introducing the new ﬁelds [23]
A
cl
t
=
1
2
(A
+
t
+ A
t
), A
q
t
= A
+
t
A
t
, (8)
where the superscript stands for classical and quantum respectively. This nomenclature derives
from the fact that the forward- and backward-in-time evolutions are equal for classical processes.
Using these new integration variables, the path integral in Eq. (7) reduces to
Λ(λ; τ) =
Z
D[A
cl
t
, A
q
t
]e
iλA
cl
τ
e
iS[A
cl
t
,A
q
t
]
. (9)
The corresponding probability distribution then reads
P(A; τ) =
1
2π
Z
−∞
e
iλA
Λ(λ; τ) =
Z
D[A
cl
t
, A
q
t
]δ(A A
cl
τ
)e
iS[A
cl
t
,A
q
t
]
. (10)
The moment generating function thus describes the ﬂuctuations of A
cl
τ
and the probability density
P(A; τ) is given by the sum over all paths for which A
cl
τ
= A, weighted by the complex exponential
of the action in complete analogy to Feynman path integrals [25]. In the case of a single observable,
the quantum ﬁeld is zero (i.e. A
+
τ
= A
τ
), making its introduction seem like an unnecessary
complication. However, when multiple observables are involved, the quantum ﬁelds are no longer
bound to be zero. Considering the ﬂuctuations of the classical ﬁelds then corresponds to a given
time-ordering prescription. We now generalize the results of this section to an arbitrary number
of (possibly non-commuting) observables at arbitrary times.
2.2 Multiple observables
The KQPD, denoted by P(A), can be deﬁned via its moment generating function which is obtained
by generalizing Eq. (2) and reads
Λ(λ)
Z
dλe
iλ·A
P(A)
Tr
(
ˆ
T
+
exp
"
i
2
Z
T
0
X
l
f
l
(τ)λ
l
ˆ
A
l
(τ)
#
ˆρ
0
ˆ
T
exp
"
i
2
Z
T
0
X
l
f
l
(τ)λ
l
ˆ
A
l
(τ)
#)
= Tr
(
ˆ
T
K
exp
"
i
2
Z
K
X
l
f
l
(τ)λ
l
ˆ
A
l
(τ)
#
ˆρ
0
)
.
(11)
Here
ˆ
T
+()
is the (anti) time-ordering operator, the variables λ
l
and A
l
are grouped in vectors λ
and A, and the operators are given in the Heisenberg picture, i.e.
ˆ
A
l
(t) =
ˆ
U(0, t)
ˆ
A
l
ˆ
U(t, 0), (12)
where
ˆ
A denotes the operator in the Schrödinger picture (where all considered operators are time-
independent). In the last line of Eq. (11),
ˆ
T
K
denotes time-ordering along the Keldysh contour
(from t = 0 to t = T and back again). The functions f
l
(τ) take into account a possibly time-
dependent coupling to the observables
ˆ
A
l
. When using the Keldysh contour, the functions f
l
take
on the same values on the forward- and the backward-in-time branch. The case of N observables
Accepted in Quantum 2017-09-19, click title to verify 4
probed instantaneously at diﬀerent times τ
l
is obtained by taking f
l
(τ
l
) = δ(t τ
l
) (for l N)
resulting in the moment generating function
Λ(λ) = Tr
n
e
i
λ
N
2
ˆ
A
N
(τ
N
)
···e
i
λ
1
2
ˆ
A
1
(τ
1
)
ˆρ
0
e
i
λ
1
2
ˆ
A
1
(τ
1
)
···e
i
λ
N
2
ˆ
A
N
(τ
N
)
o
. (13)
The apparent symmetry between terms on the left and on the right of the density matrix is equiva-
lent to treating the forward- and the backward-in-time branch of the time evolution symmetrically
and determines the time-ordering of operators when evaluating higher moments (see below).
2.3 Examples
We now introduce two well-known examples of the KQPD: the Wigner function and the Full
counting statistics (FCS). To recover the Wigner function from the moment generating function
in Eq. (11), we set
X
l
f
l
(τ)λ
l
ˆ
A
l
(τ) = δ(τ ) [λ
x
ˆx(τ ) + λ
p
ˆp(τ)] , (14)
which results in
Λ(λ
x
, λ
p
) = Tr
n
e
i
2
(λ
x
ˆx+λ
p
ˆp)
ˆρ
0
e
i
2
(λ
x
ˆx+λ
p
ˆp)
o
, (15)
where we identiﬁed ˆx = ˆx(0) and similarly for ˆp. It is straightforward to verify that the Wigner
function is indeed given by
W(x, p) =
1
(2π)
2
Z
x
p
e
x
x+
p
p
Λ(λ
x
, λ
p
) =
1
2π
Z
p
e
p
p
x +
λ
p
2
ˆρ
0
x
λ
p
2
. (16)
The Wigner function is thus recovered as the KQPD corresponding to a simultaneous observation
of ˆx and ˆp. In the next section, we discuss how the Wigner function can be used to predict a
simultaneous (but imprecise) measurement of position and momentum [26]. The moments obtained
from the Wigner function
hx
n
p
m
i
W
=
Z
dxdp x
n
p
m
W(x, p), (17)
are determined by the Weyl order [27] obtained by complete symmetrization, i.e. every possible
order is given the same weight, e.g.
hxpi
W
=
1
2
Tr {(ˆxˆp + ˆpˆx)ˆρ
0
}, hx
2
pi
W
=
1
3
Tr
(ˆx
2
ˆp + ˆxˆpˆx + ˆpˆx
2
)ˆρ
0
. (18)
This is a general result for the KQPD; when considering multiple operators at equal times, the
moments generated by the KQPD are determined by the quantum averages of the completely
symmetrized operators.
The next example we consider in this section is the FCS which is obtained from the moment
generating function in Eq. (11) upon setting
X
l
f
l
(τ)λ
l
ˆ
A
l
(τ) = λ
ˆ
A(τ), (19)
which results in
Λ(λ) = Tr
ˆ
T
+
e
i
λ
2
R
T
0
ˆ
A(τ)
ˆρ
0
ˆ
T
e
i
λ
2
R
T
0
ˆ
A(τ)
= Tr
ˆ
T
K
e
i
λ
2
R
K
ˆ
A(τ)
ˆρ
0
. (20)
The FCS is given by the Fourier transform of the last expression
F(m)
1
2π
Z
dλe
iλm
Λ(λ), (21)
and gives information on the ﬂuctuations of the time integral of
ˆ
A. The FCS has been investigated
intensively in the context of electron transport [14], where
ˆ
A corresponds to the current operator
Accepted in Quantum 2017-09-19, click title to verify 5
and one is interested in the ﬂuctuations of the charge that is transferred through a conductor
during a time interval [0, T ].
The moments of the FCS can conveniently be expressed using the Keldysh contour
Z
dm m
k
F(m) = i
k
k
λ
Λ(λ)
λ=0
=
Z
K
1
···
k
Tr
n
ˆ
T
K
ˆ
A(τ
1
) ···
ˆ
A(τ
k
)ˆρ
0
o
. (22)
We thus ﬁnd that the time-ordering corresponding to the KQPD is determined by the Keldysh time-
ordering, implying that operators are ordered along the Keldysh contour. The Wigner function as
well as the FCS have been discussed extensively in the literature. Explicit evaluations exhibiting
negativity can be found in, e.g., Refs. [19, 28].
To further illustrate the KQPD, we consider the very simple example of two subsequent spin
measurements. The operators of interest are given by the rotated Pauli matrices ˆσ
1
and ˆσ
2
. For
simplicity, we consider a qubit in a pure state |+
0
i which is an eigenstate of ˆσ
0
with eigenvalue
P
1
, Σ
2
) =
1
(2π)
2
Z
1
2
e
1
Σ
1
+
2
Σ
2
Tr
n
e
i
λ
2
2
ˆσ
2
e
i
λ
1
2
ˆσ
1
|+
0
ih+
0
|e
i
λ
1
2
ˆσ
1
e
i
λ
2
2
ˆσ
2
o
. (23)
Writing the eigenstates of ˆσ
j
in the eigenbasis of ˆσ
1
|+
0
i = α|+
1
i + β|−
1
i |+
2
i = γ|+
1
i + δ|−
1
i, (24)
we can evaluate Eq. (23), which can be written as
P
1
, Σ
2
) =
X
σ
1
=0,±1
X
σ
2
=±1
˜
P(σ
1
, σ
2
)δ
1
σ
1
)δ
2
σ
2
), (25)
with the discrete probability distribution
˜
P(+1, +1) = |α|
2
|γ|
2
,
˜
P(1, +1) = |β|
2
|δ|
2
,
˜
P(0, +1) = 2Re {αβ
γ
δ},
˜
P(+1, 1) = |α|
2
|δ|
2
,
˜
P(1, 1) = |β|
2
|γ|
2
,
˜
P(0, 1) = 2Re {αβ
γ
δ}.
(26)
Note that although the operators ˆσ
j
have eigenvalues ±1, the KQPD is non-zero also for σ
1
= 0.
Furthermore, these terms can become negative. We will discuss the origin of these terms and
their negativity in more detail in the next section. For now, we only mention that the negativity
reﬂects the fact that a measurement of ˆσ
1
disturbs a subsequent measurement of ˆσ
2
. The negative
terms thus vanish if these two operators are simultaneously measurable, [ˆσ
1
, ˆσ
2
] = 0, if the ﬁrst
measurement does not inﬂuence the initial state, [ˆσ
0
, ˆσ
1
] = 0, or if a measurement of ˆσ
2
has a
completely random outcome no matter if ˆσ
1
is measured or not. The latter is the case if ˆσ
0
= ˆσ
x
,
ˆσ
1
= ˆσ
y
, and ˆσ
2
= σ
z
. In all these cases, a measurement of ˆσ
2
yields outcomes which are
independent of the presence or absence of a previous measurement of ˆσ
1
. The KQPD for two
subsequent spin measurements is illustrated in Fig. 2.
2.4 Negativity
As mentioned above, negativity in the KQPD indicates that measuring the considered observables
necessarily inﬂuences the outcomes (i.e. measuring all of the operators results in diﬀerent statistics
than measuring each of them in a diﬀerent realization). Following Ref. [19], we corroborate this
claim by connecting negative values in the KQPD to coherent superpositions of diﬀerent eigenstates
of the observables. In the presence of such superpositions, any measurement (partly) collapses the
state, possibly inﬂuencing subsequent measurement outcomes. For clarity, we consider the case of
N observables at subsequent times as described by the moment generating function in Eq. (13).
Assuming that the eigenstates of each operator
ˆ
A
l
span the whole Hilbert space, we replace
e
i
λ
l
2
ˆ
A
l
(τ
l
)
=
X
A
α
l
e
i
λ
l
2
A
α
l
ˆ
U(0, τ
l
)|A
α
l
ihA
α
l
|
ˆ
U(τ
l
, 0), (27)
where
ˆ
A
l
|A
α
l
i = A
α
l
|A
α
l
i and the superscript α = L, R labels if the replacement happened on the
left (L) or on the right (R) of the density matrix in Eq. (13), i.e. α labels the branch (forward- or
Accepted in Quantum 2017-09-19, click title to verify 6
Figure 2: Keldysh quasi-probability distribution for two subsequent spin measurements on a pure qubit state.
While the second argument is restricted to the eigenvalues of ˆσ
2
, the ﬁrst argument can also be zero. (a)
ˆσ
0
= ˆσ
2
= ˆσ
x
and ˆσ
1
= ˆσ
z
. In this case, a measurement of ˆσ
1
inﬂuences a subsequent measurement of
ˆσ
2
which results in a negative value for
˜
P(0, 1). (b) ˆσ
0
= ˆσ
x
, ˆσ
1
= ˆσ
y
, and ˆσ
2
= σ
z
. In this case, the
measurement outcome of ˆσ
2
is completely random, no matter if the ﬁrst measurement is performed or not.
A measurement of ˆσ
1
does therefore not inﬂuence a subsequent measurement of ˆσ
2
and the KQPD is purely
positive.
backward-in-time). This allows us to take the λ-dependent factors out of the trace and obtain the
KQPD by a Fourier transformation
P(A) =
X
A
L
,A
R
δ
A
L
N
,A
R
N
δ
A
A
L
+ A
R
2
A(A
L
)A
(A
R
)hA
L
1
|ˆρ
τ
1
|A
R
1
i. (28)
Here the vectors A
L/R
have as their entries A
L/R
l
and thus denote a particular manifestation of
eigenvalues that we call a trajectory. The amplitudes associated with such a trajectory read
A(A
α
) = hA
α
N
|
ˆ
U(τ
N
, τ
N1
)|A
α
N1
ihA
α
N1
|···|A
α
2
ihA
α
2
|
ˆ
U(τ
2
, τ
1
)|A
α
1
i, (29)
the Dirac delta for vectors is deﬁned as
δ
A
A
L
+ A
R
2
=
Y
l
δ
A
l
A
R
l
+ A
L
l
2
, (30)
and ˆρ
τ
1
=
ˆ
U(τ
1
, 0)ˆρ
0
ˆ
U(0, τ
1
).
Equation (28) expresses the KQPD as a sum over pairs of trajectories. These pairs, here
denoted left (L) and right (R) trajectory, are equivalent to a single trajectory along the Keldysh
contour. A given pair of trajectories contributes with a weight that is given by the interference
term A(A
L
)A
(A
R
) times a density matrix element hA
L
1
|ˆρ
τ
1
|A
R
1
i corresponding to the “starting
point” of the trajectories. Each pair contributes to the KQPD at the argument which corresponds
to the mean value of the observable on the two trajectories, again underlining the fact that the
KQPD considers the ﬂuctuations of the classical ﬁeld in the Keldysh language (see below). Finally,
the Kronecker delta in Eq. (28) results from the trace and reﬂects the fact that the disturbance of
the last measurement is irrelevant.
From Eq. (28), we see that as long as A
L
= A
R
, the KQPD is strictly greater or equal to
zero. It is therefore the interference terms of trajectories that diﬀer from each other which allow
the KQPD to become negative. This is a genuine quantum eﬀect since these interference terms
are a direct result of a coherent superposition between diﬀerent eigenstates of the observables
ˆ
A
l
.
Although superposition can be encountered in classical (wave) theories, a coherent superposition
between two states which correspond to two diﬀerent values for an observable does not exist in
classical theories (indeed, the presence of such superpositions violates the assumptions made by
Leggett and Garg to obtain their inequalities [29], see Sec. 4.3).
In analogy to Eqs. (8) and (10), we can introduce a Keldysh rotation
A
cl
=
1
2
(A
L
+ A
R
), A
q
= A
L
A
R
, (31)
Accepted in Quantum 2017-09-19, click title to verify 7
and write the KQPD in Eq. (28) as a path integral
P(A) =
Z
D[A
cl
, A
q
]δ(A A
cl
)e
iS[A
cl
,A
q
]
, (32)
with the action
e
iS[A
cl
,A
q
]
= A(A
cl
+ A
q
/2)A
(A
cl
A
q
/2)hA
cl
1
+ A
q
1
/2|ˆρ
τ
1
|A
cl
1
A
q
1
/2i, (33)
and the integration measure
R
D[A
cl
, A
q
] =
P
A
cl
,A
q
δ
A
q
N
,0
. In this language, the KQPD is given
by the sum over all paths for which A = A
cl
, weighted by the complex exponential of the action.
From Eqs. (32) and (33), we see that any negativity in the KQPD necessarily implies the presence
of a non-zero quantum ﬁeld A
q
6= 0. In contrast to Eq. (10), which describes a single observable,
the quantum ﬁeld in Eq. (32) is not necessarily zero implying that the statistics of A
cl
, determined
by the Keldysh time-ordering, are not necessarily equal to the statistics of A
L
or A
R
.
We now illustrate these ideas with the KQPD for two subsequent qubit measurements intro-
duced in the last section. In terms of trajectories, the discrete KQPD in Eq. (26) can be expressed
as
˜
P(σ
1
, σ
2
) =
X
σ
L/R
1
=±
δ
σ
1
,(σ
L
1
+σ
R
1
)/2
hσ
2
|σ
L
1
ihσ
L
1
|+
0
ih+
0
|σ
R
1
ihσ
R
1
|σ
2
i. (34)
We thus see that the terms where σ
L
1
6= σ
R
1
are responsible for the terms
˜
P(0, σ
2
). Note that these
terms are only non-zero if the initial state is in a coherent superposition of the two eigenstates of
ˆσ
1
. One can further show that
X
σ
1
2
=±1
˜
P(σ
1
, σ
2
) = 1, (35)
which directly implies that
˜
P(0, +1) =
˜
P(0, 1). Thus any contribution to Eq. (34) from pairs
of trajectories with σ
L
1
6= σ
R
1
results in negativity in the KQPD as illustrated in Fig. 2.
We close this section with a brief comparison to other non-classicality indicators based on
negative quasi-probabilities. As mentioned above, the Wigner function (at time t) is the KQPD for
the observables ˆx(t) and ˆp(t). A negative Wigner function implies an impossibility of describing
the corresponding state classically in the canonical position-momentum phase space (i.e. as a
mixture of states with well deﬁned positions and momenta). The KQPD generalizes this concept
to arbitrary observables. Negativity in the KQPD implies that it is impossible to describe the
system and its time evolution as a mixture of states with well deﬁned values for the observable
outcomes. We note that negative values in a discrete analogue of the Wigner function have recently
been connected to the computational speed-up of a quantum computer [3032].
In quantum optics, one often considers negativity, or strong singularity, in the Glauber-Sudarshan
P function as an indicator of non-classicality [33, 34]. This criterion captures all states with non-
positive Wigner functions as well as squeezed states and entangled states [35]. A non-classical P
function implies that the state can not be interpreted as a mixture of coherent states which indi-
cates that the state can not be fully described by classical optics [34]. We note that the P function
has recently been extended to multiple times [3638]. In contrast to the KQPD, this distribution
is relevant for photo-detection and its moments correspond to normal ordering and regular time
ordering.
While these QPDs only contain information on the states and possibly their time-evolution,
quasi-probability representations also describe measurements with QPDs [36]. It has been shown
that a positive representation (i.e. a representation that only uses positive probability distributions)
of all states and measurements is impossible. Restricting the states and measurements to a given
experimental procedure, a procedure is then said to require non-classical resources if no positive
representation can be found [4]. It is an interesting question if negativity in the KQPD implies the
absence of a positive representation for a corresponding experimental procedure.
3 von Neumann measurements
In this section, we consider von Neumann type measurements and we show that the corresponding
probability distributions can be obtained by convolving the KQPD with the back-action and impre-
cision associated with the measurement. We start by considering the instantaneous measurement
Accepted in Quantum 2017-09-19, click title to verify 8
of a single observable before moving on to subsequent instantaneous measurements, simultaneous
measurements, and measurements of ﬁnite duration.
3.1 Instantaneous measurements
Following von Neumann [21], we model the measurement of an observable
ˆ
A at time τ by coupling
it to a detector described by the canonically conjugate observables ˆr and ˆπ via the interaction
Hamiltonian
ˆ
H
m
= δ(t τ)χ
ˆ
Aˆπ, (36)
where χ denotes the coupling strength. This interaction induces the time-evolution
ˆ
U
m
= exp(
ˆ
Aˆπ),
which displaces the detector coordinate ˆr by an amount that depends on the state of the system.
Projectively measuring the position of the detector then yields the measurement of
ˆ
A. Measuring
a single observable in this way yields the probability distribution (see App. B)
P (A) =
Z
χdA
0
hχ(A A
0
)|ˆρ
det
|χ(A A
0
)iP(A
0
; τ), (37)
where ˆρ
det
is the density matrix of the detector, |χ(A A
0
)i is an eigenstate of the operator
ˆr, and P(A
0
; τ) is the KQPD for a single observable given in Eq. (4). Note that the position
uncertainty in the initial detector state directly translates into a measurement uncertainty. A
projective measurement can be obtained by choosing a position eigenstate as the initial state for
the detector. In this case, the measured distribution coincides with the KQPD. This is due to the
fact that the back-action of the measurement is irrelevant for a single observable since we are not
interested in the post-measured state.
For multiple subsequent (instantaneous) measurements, we couple each observable to a diﬀerent
detector. The measurement Hamiltonian then reads
ˆ
H
m
=
N
X
l=1
δ(t τ
l
)χ
l
ˆ
A
l
ˆπ
l
. (38)
In contrast to a single observable, each measurement can now inﬂuence all subsequent measure-
ments. To capture this back-action eﬀect, we evaluate the KQPD with the additional term
P
l
δ(t τ
l
)γ
l
ˆ
A
l
in the Hamiltonian, where the γ
j
are random variables which are distributed
according to the detector states (see below). This yields the modiﬁed moment generating function
Λ(λ; γ) = Tr
e
i
λ
N
2
+γ
N
ˆ
A
N
(τ
N
)
···e
i
(
λ
1
2
+γ
1
)
ˆ
A
1
(τ
1
)
ˆρ
0
e
i
(
λ
1
2
γ
1
)
ˆ
A
1
(τ
1
)
···e
i
λ
N
2
γ
N
ˆ
A
N
(τ
N
)
,
(39)
which can be obtained by replacing λ
l
λ
l
± 2γ
l
in Eq. (13), where the + () sign corresponds
to terms on the left (right) side of the density matrix.
The probability distribution which describes subsequent von Neumann measurements can then
be shown to read (for a derivation, see App. B)
P (A) =
Z
dA
0
dγ
"
Y
l
W
l
(χ
l
[A
l
A
0
l
], γ
l
l
)
#
P(A
0
; γ), (40)
where W
l
(r, γ) denotes the Wigner function of detector l and P(A
0
; γ) is the KQPD deﬁned as the
Fourier transform of the moment generating function in Eq. (39). The measurement outcomes are
determined by a convolution of the KQPD and the Wigner functions of the detectors. We can thus
interpret the KQPD as describing the intrinsic ﬂuctuations of the observables of interest which is
being inﬂuenced by the eﬀects of the measurement [17]. The position distributions of the detectors
determine the precision of the measurement while the momentum distributions determine the back-
action via the γ
l
-dependent terms in the KQPD. The Heisenberg uncertainty relation implies a
trade-oﬀ between these two eﬀects, indicating that we can never eliminate both inﬂuences.
As an example, we consider a measurement of position ˆx at τ = 0 followed immediately by a
momentum ˆp measurement. The probability distribution describing such a measurement reads
W
xp
(x, p) =
Z
dx
0
dp
0
x
p
W
p
(χ
p
[p p
0
], γ
p
p
)W
x
(χ
x
[x x
0
], γ
x
x
)W(x
0
, p
0
+ γ
x
). (41)
Accepted in Quantum 2017-09-19, click title to verify 9
Here W
x/p
denotes the Wigner function of the detector measuring x/p and W without a subscript
denotes the KQPD which in this case is given by the Wigner function of the system. The position
measurement disturbs the momentum of the system as reﬂected by the second argument of the
KQPD. Note that it is the momentum distribution of the x-detector which determines the eﬀect
of this disturbance. This reﬂects the trade-oﬀ between the precision of the position measurement
(which requires a peaked position distribution for the x-detector) and the disturbance of the sub-
sequent momentum distribution (which is minimized for a peaked momentum distribution for the
x-detector). Note that the last expression is independent of the momentum distribution of the
p-detector (this can be seen by performing the integral over γ
p
), reﬂecting the fact that the back-
action of the momentum measurement is irrelevant since it is the last measurement. Inverting the
order of the measurements yields a probability distribution which is given by the last expression
upon exchanging W(x
0
, p
0
+ γ
x
) W(x
0
γ
p
, p
0
).
Returning to arbitrary observables and choosing Gaussian Wigner functions for the detectors
W
l
(r, γ) =
1
2πσ
r
l
σ
γ
l
e
r
2
2σ
2
r
l
e
γ
2
2σ
2
γ
l
, (42)
we can write the probability distribution in Eq. (40) in terms of trajectories [cf. Eq. (28)] yielding
P (A) =
X
A
L
,A
R
δ
A
L
N
,A
R
N
Y
l
1
2πσ
imp,l
e
1
2σ
2
imp,l
A
l
A
L
l
+A
R
l
2
2
e
σ
2
BA,l
2
(
A
L
l
A
R
l
)
2
× A(A
L
)A
(A
R
)hA
L
1
|ˆρ
τ
1
|A
R
1
i.
, (43)
where σ
imp,l
= σ
r
l
l
and σ
BA,l
= σ
γ
l
χ
l
. Comparing the last expression with the KQPD in
Eq. (28), we can clearly identify the eﬀects of the measurement imprecision and back-action. The
ﬁrst Gaussian in the last expression replaces the Dirac delta in Eq. (28). It implies that even terms
with A
l
6= (A
L
l
+ A
R
l
)/2 contribute to the probability of measuring A
l
. This eﬀect captures the
imprecision of the measurement which becomes more pronounced as the uncertainty in the position
of detector l grows. The second Gaussian in the last expression reduces the weight of terms which
have A
L
l
6= A
R
l
, i.e. a ﬁnite quantum ﬁeld in the Keldysh language. This eﬀect is due to the
measurement back-action which (partly) collapses the coherent superpositions between eigenstates
of
ˆ
A
l
. The coupling Hamiltonian in Eq. (38) entangles the system with the detector leading to
dephasing in the reduced state of the system. As the uncertainty in the detector momentum grows,
this eﬀect becomes more pronounced. From the Heisenberg uncertainty relation, we ﬁnd
σ
imp,l
σ
BA,l
= σ
r
l
σ
γ
l
1
2
, (44)
ensures that the measured distribution is positive even if the KQPD exhibits negativities. We note
that coupling the system to additional degrees of freedom, such as a thermal reservoir, introduces
dephasing analogously to the measurement back-action [19].
To illustrate the eﬀect of the measurement, we return to our example of two subsequent spin
measurements [cf. Eq. (23)]. Since the back-action of the second measurement is irrelevant, we
consider a measurement of intermediate strength of ˆσ
1
followed by a projective measurement of
ˆσ
2
. In terms of the discrete KQPD given in Eq. (26), we ﬁnd for the distribution describing a von
Neumann measurement (with Gaussian detectors)
P (σ
1
, σ
2
) =
1
2πσ
imp
X
σ
0
1
=0,±1
e
1
2σ
2
imp
(σ
1
σ
0
1
)
2
e
2σ
2
BA
δ
σ
0
1
,0
˜
P(σ
0
1
, σ
2
). (45)
The eﬀect of the measurement is illustrated in Fig. 3. We note that the distribution of measure-
ment outcomes for the ﬁrst measurement depends very strongly on the outcome of the second
measurement. Speciﬁcally, we see from Fig. 3 that only considering events where the second (pro-
jective) measurement yields 1 allows one to distinguish the outcomes of the ﬁrst measurement
much better (for a given measurement strength). This eﬀect, which arises from the negativity of
the KQPD, is responsible for the occurrence of anomalous weak values as discussed in Sec. 4.2.
Note that subsequent von Neumann measurements can also be used to model continuous mea-
surements, a subject which has received a lot of attention recently, see e.g. [4044].
Accepted in Quantum 2017-09-19, click title to verify 10
Figure 3: Eﬀect of measurement imprecision and back-action on the KQPD. Here ˆσ
0
= ˆσ
2
= ˆσ
x
and ˆσ
1
= ˆσ
z
to
maximize the negativity in the KQPD. The measurement imprecision and back-action modify the KQPD resulting
in a positive probability distribution. Here σ
imp
σ
BA
= 1/2, minimizing the inﬂuence of the measurement.
Depending on the measurement outcome of the second (projective) measurement, the distribution for the ﬁrst
measurement results in two clearly distinguishable outcomes or not. This is directly related to the terms
˜
P(0, σ
2
)
which are responsible for the negative values in the KQPD. Note that the peaks are shifted with respect to the
eigenvalues of ˆσ
1
(dashed lines). This eﬀect allows for the occurrence of anomalous weak values as discussed
in Sec. 4.2.
3.2 Simultaneous measurements of position and momentum
We now turn to the simultaneous measurement of position and momentum. In analogy to Eq. (36)
(and following Arthurs and Kelly [26]), we model the measurement by the simultaneous coupling
of two detectors to the position and momentum of the system respectively
ˆ
H
m
= δ(t) (χ
x
ˆxˆπ
x
+ χ
p
ˆpˆπ
p
) . (46)
It can be shown that this measurement is described by the probability distribution [45]
W (x, p) =
Z
dx
0
dp
0
x
p
W
p
(χ
p
[pp
0
], γ
p
p
)W
x
(χ
x
[xx
0
], γ
x
x
)W(x
0
γ
p
/2, p
0
+γ
x
/2). (47)
The arguments of the last Wigner function reﬂect the fact that in a simultaneous measurement,
both position and momentum are being disturbed. Note the close similarity between the last
expression and the probability distribution describing subsequent measurements of position and
momentum [cf. Eq. (41)]. Notably, the KQPD is in both cases given by the Wigner function.
This is a consequence of the commutation relation [ˆx, ˆp] = i and does not hold in general for
non-commuting observables.
For a detailed discussion on the von Neumann type simultaneous measurement of position and
momentum, we refer the reader to Ref. [22]. Here we only discuss the case of Gaussian states of
minimal uncertainty for the detectors (k = x, p)
W
k
(r, γ) =
1
π
e
r
2
2
k
e
γ
2
σ
2
k
, (48)
x
σ
p
/(χ
x
χ
p
) = 1/2. This constrained minimizes the disturbance
of the measurement and yields the Husimi Q-function as the measured probability distribution
[22, 46]
W (x, p) =
1
π
Z
dx
0
dp
0
e
(xx
0
)
2
2σ
2
e
2σ
2
(pp
0
)
2
W(x
0
, p
0
), (49)
where σ = σ
x
x
. In addition to this von Neumann type measurement, the Husimi Q-function
also describes heterodyne measurements of two ﬁeld quadratures [47].
3.3 Measuring the FCS
In this section, we consider a von Neumann type measurement of the time-integral of an observable
[15]. To this end, we couple a detector to the observable of interest for a ﬁnite amount of time
ˆ
H
m
= Θ(t)Θ(T t)χ
ˆ
Aˆπ. (50)
Accepted in Quantum 2017-09-19, click title to verify 11
This measurement can be shown to be described by the probability distribution [15]
F (m) =
Z
W(χ[m m
0
], γ)F(m
0
; γ), (51)
where F(m; γ) is obtained through Eqs. (20) and (21) upon adding the term Θ(t)Θ(T t)γ
ˆ
A
to the Hamiltonian. Again, we ﬁnd a trade-oﬀ between measurement imprecision, depending on
the spread in the detector position, and back-action which depends on the spread in the detector
momentum. As discussed in Ref. [19], the last equation can be expressed in terms of trajectories
yielding (for a Gaussian detector)
F (m) =
1
2πσ
imp
X
A
L
,A
R
δ
A
L
T
,A
R
T
e
(
m
1
2
m
L
1
2
m
R
)
2
2σ
2
imp
e
σ
2
BA
2
(m
L
m
R
)
2
hA
L
0
|ˆρ|A
R
0
iA(A
L
)A
(A
R
),
(52)
where the trajectory amplitudes are obtained by dividing the time evolution into N time-slices of
inﬁnitesimal duration δt and inserting identities resolved in the eigenstates of
ˆ
A after each slice
A(A
α
) = hA
α
T
|
ˆ
U(T , T δt)|A
α
T δt
ihA
α
τδt
|···|A
α
δt
ihA
α
δt
|
ˆ
U(δt, 0)|A
α
0
i. (53)
The discrete time-integral over such a trajectory is denoted by m
α
=
P
l
A
α
l
δt. The Heisenberg
uncertainty relation applied to the detector variables again enforces σ
imp
σ
BA
1/2, reﬂecting
the trade-oﬀ between smearing out the negativities in the FCS by measurement imprecision and
reducing the weight of negative terms (which necessarily have m
L
6= m
R
) by measurement back-
action in complete analogy to Eq. (43).
This concludes our treatment on von Neumann measurements. Expressing the measured prob-
ability distributions in terms of the KQPD suggests that the latter should be interpreted as de-
scribing the intrinsic ﬂuctuations of the observables of interest. In the next section, we will build on
this interpretation and demonstrate the utility of studying the KQPD to understand the quantum
behavior of dynamic systems.
4 Applications of the KQPD
We now turn to situations where it is beneﬁcial to study the KQPD in the absence of any mea-
surement device. We are especially interested in situations where the KQPD fails to be positive.
As discussed in Sec. 2.4, this indicates non-classical behavior. The most prominent KQPD is by
far the Wigner function, which is a representation of the density matrix and is used extensively to
illustrate the state of a quantum system and to capture its non-classicality [28, 48]. An example
for a KQPD which goes beyond instantaneous states is provided by the FCS which describes the
ﬂuctuations of an observable over a certain time-interval. While the FCS has mostly been em-
ployed to describe electronic transport in systems where it is a purely positive function, its utility
for detecting non-classical behavior has been recognized [1619, 39, 49]. Here we discuss three ad-
ditional applications of the KQPD, stressing its versatility and its utility in capturing non-classical
behavior in the dynamics of quantum systems.
4.1 Work ﬂuctuations
Since work is not a state function, it is in general not suﬃcient to interrogate a system at a single
time in order to determine the work performed by the system. As we have seen above, performing
measurements at multiple times generally results in a probability distribution which is inﬂuenced
by the measurement itself. As a consequence, the observed amount of work will in general depend
on the way it is measured [50]. Here we use the KQPD to describe work as a measurement-
independent, ﬂuctuating quantity. We note that the KQPD has been used before to describe work
ﬂuctuations [5153].
In classical statistical mechanics, work is deﬁned as the energy change in the macroscopic,
accessible degrees of freedom, while heat is the energy change in the microscopic, inaccessible
degrees of freedom. In quantum systems, all degrees of freedom are in principle microscopic. We
Accepted in Quantum 2017-09-19, click title to verify 12
can however still deﬁne work as being the energy change mediated by the accessible degrees of
freedom. Here we consider two diﬀerent scenarios: In the ﬁrst, one has access to the energy stored
in a part of the system that we call the work storage device. The energy increase of this device is
then interpreted as the work performed by the system. This is the analog of determining the work
performed on a suspended weight by measuring its height. In the second scenario, one has access
to the power produced by the system and the work is deﬁned as the time-integral of the power.
This captures the situation in thermoelectric heat engines, where one can access the power through
a measurement of the electrical current. In terms of the suspended weight, this corresponds to the
analog of measuring the power by measuring the velocity of the weight.
For the ﬁrst scenario, we write the total Hamiltonian as
ˆ
H =
ˆ
H
w
+
ˆ
H
rest
, where
ˆ
H
w
denotes the
Hamiltonian of the work storage device. Work ﬂuctuations can then be described by the KQPD
corresponding to the diﬀerence of two subsequent energy measurements on the work storage device
P(w; τ) =
1
2π
Z
dλe
iλw
Tr
n
e
i
λ
2
ˆ
H
w
ˆ
U(τ, 0)e
i
λ
2
ˆ
H
w
ˆρ
0
e
i
λ
2
ˆ
H
w
ˆ
U(0, τ)e
i
λ
2
ˆ
H
w
o
=
X
E
τ
,E
L
,E
R
δ
w
E
τ
E
L
+ E
R
2

Tr
n
hE
τ
|
ˆ
U(τ, 0)|E
L
ihE
L
|ˆρ
0
|E
R
ihE
R
|
ˆ
U(0, τ)|E
τ
i
o
.
(54)
Here
ˆ
H
w
|E
j
i = E
j
|E
j
i and τ denotes the time interval between the two measurements. For initial
states which commute with
ˆ
H
w
, we ﬁnd E
L
= E
R
and we recover the probability for obtaining w as
the diﬀerence of two projective energy measurements [54]. Any negativity in the above distribution
is thus due to the coherence in the initial state. We note that the probability distribution obtained
by performing two Gaussian energy measurements, as discussed in Ref. [55], is obtained from
Eq. (54) by including the imprecision and back-action of the measurement as discussed in Sec. 3.
The moments generated by the KQPD are given by
hw
k
i =
k
X
l=0
1
2
k
k
l
Tr
ˆ
T
ˆ
H
w
(τ)
ˆ
H
w
l
ˆ
T
+
ˆ
H
w
(τ)
ˆ
H
w
(kl)
ˆρ
0
, (55)
where
ˆ
H
w
(τ) =
ˆ
U(0, τ)
ˆ
H
w
ˆ
U(τ, 0). For the ﬁrst two moments, i.e. k = 1, 2 in the last expression,
this reduces to
hw
k
i = Tr
ˆ
H
w
(τ)
ˆ
H
w
k
ˆρ
0
, for k = 1, 2. (56)
The mean work obtained through the KQPD is thus just given by the diﬀerence of the average
energy in the work storage device at t = τ and t = 0. This is in contrast to the probability
distribution obtained by performing two projective energy measurements, where the mean value
only coincides with Eq. (56) for initial states that commute with
ˆ
H
w
. For an explicit evaluation
of Eq. (54), we refer to Ref. [53]. See Ref. [56] for an experimental reconstruction of the work
distribution by measuring its moment generating function (for an initial state that commutes with
the Hamiltonian).
ˆ
P , the KQPD describing work ﬂuctuations is given
by the FCS of the power operator. The KQPD and its moments are then given by Eqs. (20-22)
upon replacing
ˆ
A with
ˆ
P . This scenario is relevant for thermoelectric devices, where the power
operator is given by the electrical current operator times the voltage across the device
ˆ
P =
ˆ
IV .
In this case, the work ﬂuctuations reduce to the problem of charge FCS [14]. We now return to
systems described by
ˆ
H =
ˆ
H
w
+
ˆ
H
rest
, where we take the Hamiltonians to be time-independent for
simplicity (see Ref. [57] for a discussion on the power operator for a time-dependent Hamiltonian).
In this case, the power operator is given by the energy change of the work storage device per
unit time, i.e.
ˆ
P = i[
ˆ
H,
ˆ
H
w
]. If the power operator commutes with the total Hamiltonian, i.e.
[
ˆ
P ,
ˆ
H] = 0, we ﬁnd
e
i
λ
2
ˆ
H
w
ˆ
U(τ, 0)e
i
λ
2
ˆ
H
w
=
ˆ
T
+
e
i
λ
2
R
τ
0
dt
ˆ
P (t)
. (57)
This implies that the FCS of the power operator is equal to the KQPD for two energy measurements
given in Eq. (54) [51]. The diﬀerence of the two scenarios considered here is then only relevant when
taking into account the disturbance from an actual measurement. For a thermoelectric device, the
Accepted in Quantum 2017-09-19, click title to verify 13
energy in the work storage device corresponds to the electric potential energy stored in an electronic
contact, i.e.
ˆ
H
w
= V
ˆ
Q, where
ˆ
Q denotes the charge operator in the contact with
ˆ
I = i[
ˆ
H,
ˆ
Q].
The work provided by a thermoelectric device can thus either be determined by measuring the
electrical current (determining the power) or by measuring the number of electrons in an electric
contact (determining the energy in the work storage device). We note that both approaches have
been considered in the theory of charge FCS [14].
4.2 Weak values
The weak value of an observable
ˆ
A is obtained by preparing the system in some state |Ii, performing
a weak measurement of
ˆ
A, and post-selecting the measurement outcomes on the ﬁnal state |F i
[58]. Describing the weak measurement as a von Neumann measurement, the post-measured state
of the detector is determined by the weak value
A
w
=
hF |
ˆ
A|Ii
hF |Ii
, (58)
and the average of the measurement outcome is given by the real part of the last expression [58, 59].
Any possible time-evolution between pre-selection, measurement and post-selection is absorbed in
the deﬁnition of the initial and ﬁnal states. In the KQPD framework, we can describe the weak
value measurement using the distribution describing two subsequent measurements, one for the
weak measurement and one for the post-selection. For simplicity, we do not treat the pre-selection
as an additional measurement but impose it by setting ˆρ
0
= |IihI|. The relevant KQPD then reads
P(A, B) =
1
(2π)
2
Z
A
B
e
A
A
e
B
B
Tr
n
e
i
λ
B
2
ˆ
B
e
i
λ
A
2
ˆ
A
|IihI|e
i
λ
A
2
ˆ
A
e
i
λ
B
2
ˆ
B
o
=
X
B
F
,A
L
,A
R
δ (B B
F
) δ
A
A
L
+ A
R
2
hB
F
|A
L
ihA
L
|IihI|A
R
ihA
R
|B
F
i,
(59)
where |B
F
i are eigenstates of
ˆ
B. The KQPD corresponding to the measurement of the weak value
is then obtained by taking into account the post-selection, i.e. ﬁxing the argument corresponding
to the second measurement
P
w
(A) =
P(A, F )
R
dAP(A, F )
=
1
|hF |Ii|
2
X
A
L
,A
R
δ
A
A
L
+ A
R
2
hF |A
L
ihA
L
|IihI|A
R
ihA
R
|F i. (60)
The statistics of weak value measurements (with post-selection performed by a projective mea-
surement) can be described using the last expression by including measurement imprecision and
back-action as discussed in Sec. 3. In the weak measurement limit, the disturbance of the mea-
surement does not alter the average value of the KQPD, which is equal to the real part of Eq. (58)
hAi
w
=
Z
dA AP
w
(A) = Re
(
hF |
ˆ
A|Ii
hF |Ii
)
. (61)
From Eq. (60), we see that the KQPD vanishes for arguments that lie outside the spectrum
of
ˆ
A. However Eq. (61) can take on values which lie outside the spectrum of
ˆ
A. In such cases,
the weak value is called anomalous. We now show that an anomalous weak value implies that the
distribution P
w
(A) takes on negative values. To see this, assume that P
w
(A) 0. In that case,
we ﬁnd
hAi
w
=
Z
A
max
A
min
dA AP
w
(A) A
max
Z
A
max
A
min
dAP
w
(A) = A
max
, (62)
where A
min (max)
is the smallest (largest) eigenvalue of
ˆ
A and a similar inequality restricts hAi
w
A
min
. Note that the inequality in Eq. (62) is only valid for a strictly positive KQPD. An anomalous
weak value thus necessarily implies that the KQPD P
w
(A) [and also P(A, B), cf. Eq. (60)] exhibits
negativity indicating non-classical behavior. This is an agreement with Ref. [60], where anomalous
weak values have been shown to imply contextuality. We note that weak values have been connected
Accepted in Quantum 2017-09-19, click title to verify 14
to quasi-probability distributions before [59]. In particular, the average momentum at a given
position that is obtained from the Wigner function has been shown to be determined by the weak
value of the momentum
R
dp pW(x, p)
R
dpW(x, p)
= Re
hx|ˆp|Ii
hx|Ii
, (63)
which follows from our approach when taking
ˆ
A = ˆp and
ˆ
B = ˆx.
Finally, we note that the distribution in Eq. (60) is relevant for a measurement of
ˆ
A with an
arbitrary strength. It has the variance
h
ˆ
A
2
i
w
h
ˆ
Ai
w
2
=
1
2
Re
(
hF |
ˆ
A
2
|Ii
hF |Ii
)
1
2
Re
(
hF |
ˆ
A|Ii
hF |Ii
)!
2
. (64)
For
ˆ
A = ˆp and |F i = |xi, the last expression is proportional to the quantum potential energy [59], a
central object in Bohmian mechanics [61]. See also Ref. [62] for a connection between weak values
and von Neumann measurements of arbitrary strength and post-selection.
For our example of two subsequent spin measurements, we ﬁnd [using Eqs. (23), (26), and (60)]
P
w
(σ
1
) =
|α|
2
|δ|
2
δ(σ
1
1) + |β|
2
|γ|
2
δ(σ
1
+ 1) 2Re{αβ
γ
δ}δ(σ
1
)
|αδ βγ|
2
, (65)
where we chose |−
2
i as the ﬁnal state |F i. We see that for a suﬃciently small denominator, the
probability density associated with σ
1
= ±1 can become very large as long as the presence of
the negative term ensures the normalization of P
w
(σ
1
). In this way, the weak value, determined
by P
w
(1) P
w
(1) can be ampliﬁed resulting in values far outside the range [1, +1]. This can
however only occur for a suﬃciently large negative term in the KQPD as can be seen from Eq. (65)
which implies
hσ
1
i
w
> 1 Re{αβ
γ
δ} > |β|
2
|γ|
2
. (66)
4.3 Leggett-Garg inequalities
In Ref. [29], Leggett and Garg derived inequalities which hold under the assumption of macroscopic
realism [i.e., a (macroscopic) system will at all times be in one (and only one) of the states
available to it] and non-invasive measurability. Clearly quantum systems do not fulﬁll either of
these assumptions and can thus violate Leggett-Garg inequalities. In contrast to Bell inequalities,
where multiple spatially separated parties apply local operations, Leggett-Garg inequalities are
obtained by probing the same system at multiple times. We note that even if a Leggett-Garg
inequality is violated, one can always argue that this happened because the measurement did
inﬂuence the system, and therefore all subsequent measurements, in an unexpected (but in principle
avoidable) way. This is known as the clumsiness loophole. Since it is not possible to certify the
non-invasiveness of a measurement, the clumsiness loophole can never be closed [63].
We will now show that the violation of a particular Leggett-Garg inequality directly implies
the occurence of negative values in the corresponding KQPD. The considered inequality reads
K = C
21
+ C
32
C
31
1, (67)
where the correlators are given by
C
ij
=
X
Q
i
,Q
j
=±1
Q
i
Q
j
P
ij
(Q
i
, Q
j
). (68)
Here P
ij
(Q
i
, Q
j
) is the probability of obtaining the outcome Q
i
at time t
i
and Q
j
at t
j
if the
system is measured at these two times only. The measurement outcomes are restricted to Q
i
= ±1.
Violation of the inequality in Eq. (67) implies that the correlators C
ij
can not be obtained from a
positive probability distribution which describes all three measurements and is independent of the
choice of measurements that are performed [29, 63], i.e.
K > 1 C
ij
6=
X
Q
1
,Q
2
,Q
3
=±1
Q
i
Q
j
P (Q
3
, Q
2
, Q
1
), (69)
Accepted in Quantum 2017-09-19, click title to verify 15
with P (Q
3
, Q
2
, Q
1
) 0.
For projective measurements of a dichotomic observable
ˆ
Q, the correlator to be used in Eq. (67)
is given by [64]
C
ij
= Tr
1
2
n
ˆ
Q(t
i
),
ˆ
Q(t
j
)
o
ˆρ
0
. (70)
We now introduce the KQPD for three subsequent measurements of
ˆ
Q
P(Q
3
, Q
2
, Q
1
) =
1
(2π)
3
Z
3
2
λ
1
e
3
Q
3
+
2
Q
2
+
1
Q
1
× Tr
n
e
i
λ
3
2
ˆ
Q(t
3
)
e
i
λ
2
2
ˆ
Q(t
2
)
e
i
λ
1
2
ˆ
Q(t
1
)
ˆρ
0
e
i
λ
1
2
ˆ
Q(t
1
)
e
i
λ
2
2
ˆ
Q(t
2
)
e
i
λ
3
2
ˆ
Q(t
3
)
o
.
(71)
It is straightforward to show that the correlators in Eq. (70) can directly be obtained from the
KQPD
C
ij
=
Z
dQ
1
dQ
2
dQ
3
Q
i
Q
j
P(Q
3
, Q
2
, Q
1
). (72)
If the KQPD was strictly positive, it would thus provide a positive probability distribution which
describes all three measurements and is independent of the choice of measurements that are per-
formed. The existence of such a probability distribution prevents the Leggett-Garg inequality to
be violated [63] (see also App. C). Inverting the argument, we conclude that a violation of the
Leggett-Garg inequality in Eq. (67) implies negative values for the KQPD in Eq. (71) and thus
non-classical behavior in the time-evolution of
ˆ
Q.
While we discussed projective measurements, the KQPD could also ﬁnd applications when
considering Leggett-Garg inequalities obtained from weak continuous measurements [65, 66].
5 Conclusions
The KQPD provides a generalization of the Wigner function and the FCS to describe arbitrary
observables within a quasi-probability formalism. This provides a uniﬁed framework for the de-
scription of ﬂuctuations in dynamic quantum systems. Negative values of the KQPD can be directly
linked to an interference between paths through Hilbert space which visit diﬀerent eigenstates of
the observables of interest. This requires coherent superpositions of states that correspond to
diﬀerent measurement outcomes, implying that negativity in the KQPD indicates non-classical
behavior. As we have shown, this negativity is witnessed by anomalous weak values and violated
Leggett-Garg inequalities (see also Refs. [67, 68] for a connection between anomalous weak values
and Leggett-Garg inequalities). The versatility of the KQPD allows to describe the problem at
hand with a tailor-made quasi-probability distribution which focuses on the observables of inter-
est. When measuring these observables in <