1. Introduction
Much like the celebrated Bell inequalities, which shed light on the deeply non-classical properties of spatial correlations encountered in entangled systems, quantum mechanics posseses a rich temporal structure that distinguishes it from classical physics. In 1985 Leggett and Garg explored this structure by introducing the concept of macrorealism [
1]. In essence, macrorealism can be condensed into two main assumptions about the temporal properties of physical observables within any
classical description of physics [
1,
2,
3,
4,
5,
6,
7,
8];
Macrorealism per se: physical observables take on well-defined values at all times independent of the act of observation.
Non-invasive measurability: in principle it is possible to measure the value of an observable without changing the subsequent evolution of the system.
Assumptions (i) and (ii) can be used to derive mathematical inequalities, the so-called Leggett-Garg inequalities, that serve to test the macrorealism of physical observables. Violations of these inequalities subsequently rule out what would be expected in a classical system, and this quantum behaviour has now been confirmed experimentally in a variety of settings [
9,
10,
11,
12,
13,
14].
In the quantum regime, thermodynamic quantities such as fluctuating work and heat cannot be represented by hermitian observables, but are conventionally defined via multi-time projective measurements performed on the system [
15,
16,
17]. For a closed quantum system, one way of defining the fluctuating work done on the system driven out of equilibrium is by the difference in energy eigenvalues observed at the start and end of its evolution. This framework is commonly referred to as the two-projective-measurement scheme, and serves as a route to many of the known fluctuation theorems such as the Jarzynski equality [
15,
16] and Tasaki-Crooks relation [
18]. Given that these results mirror the corresponding classical fluctuation relations [
19,
20], it is often assumed that work is simply a classical stochastic variable even within the quantum regime. However, the influence of non-classical temporal correlations that arise from two-time quantum measurements on the statistics of fluctuating work have yet to be fully understood. One aspect of this is the fact that work measurements remove coherences in the energy basis and can affect the future evolution of the quantum system, modifying the average work done during the process [
21,
22]. Alternative definitions of quantum work related to weak measurement have been proposed in order to circumvent this effect of measurement disturbance on the statistics of work [
22,
23,
24]. However, it has been shown that the resulting quantum work distributions are not generally positive-definite [
25]. The emergence of negative quasi-probabilities is a signature of quantum behaviour, and hints at a link to violations of the Leggett-Garg inequalities [
6,
26,
27]. In a similar vein, violations of macrorealism have also been related to the presence of anomalous weak values in quantum systems [
8]. Recent work by Blattmann and Mølmer [
28] has successfully linked violations of macrorealism to quantum work in the standard TPM approach by utilising the entropic Leggett-Garg inequalities. In their approach one compares the Shannon entropy of the work distribution over different intervals of time. However, the Shannon entropy is not well-defined if the work distribution fails to be positive, and so the entropic Leggett-Garg inequalities cannot be applied to situations in which the work distribution becomes a quasi-probability.
In this paper we will utilise the assumptions of Leggett and Garg to demonstrate that there exists no general macrorealistic description of work for quantum systems driven out of equilibrium. In particular, we show that quantum temporal correlations between energy measurements performed at different times influence the statistical moments of the fluctuating work done on the system during a non-equilibrium process. This result is shown to hold for three different definitions of quantum work: the two-projective measurement (TPM) scheme [
15], the full-counting statistics (FCS) [
24] and the Margenau-Hill (MH) work distribution [
23]. Crucially the inequalities that we derive can be used to test for violations of macrorealism in both strong and weak measurement schemes, regardless of whether or not the work distribution is positive or not.
The paper is organised as follows: we first introduce a set of Leggett-Garg inequalities for the moments of fluctuating work, and then consider a driven two-level system and show that the inequalities can be violated. Following that we introduce an alternative set of Leggett-Garg inequalities for the moment-generating function, and apply these inequalities to alternative definitions of quantum work that are based on weak measurement, namely the FCS and MH definitions, subsequently showing that violations of macrorealism can also occur. Finally we conclude with a discussion of our results.
2. Inequalities for Moments of Work
We first recall the setup for the original Leggett-Garg inequalities [
1]. First consider performing three protocols in which the spin
of a qubit is projectively measured at two times within a set of three times
. For each of the three protocols one can obtain the temporal correlation function for the values of the spin at times
and
, denoted
. The macrorealism assumptions (i) and (ii) imply that there exists a three-time probability distribution
such that the distributions describing the statistics of each individual protocol can be obtained as marginals of this three-time distribution, eg.
and so on [
4]. Note that while assumption (i) implies the existence of a three-time probability with the correct marginals, assumption (ii) guarantees that this distribution is the same for all three separate experiments [
4]. Finally, using the marginal properties of
yields the following Leggett-Garg inequality relating the correlation functions for the three protocols [
1];
This inequality holds for any dichotomic observable. For a simple qubit, the spin at time
can be represented by a combination of Pauli matrices;
. By performing successive projective measurements of the spin, the correlation functions can be obtained for each of the three protocols. This leads to a violation the RHS of Equation (
1) which can take a maximum value of
[
4]. This example illustrates the failure of macrorealism for quantum systems.
It is also possible to derive a set of Leggett-Garg inequalities reminiscent of Equation (
1) for the moments of fluctuating work in a closed quantum system driven in time. For simplicity we consider a system that can occupy one of two fixed energy states, which we denote by
and
, at three points in time
during the driving process. For the statistics of work measured during a particular time interval
, the
k’th moment of fluctuating work is defined as
Here the energies occupied by the system at time
are denoted by
and we assume that the possible work values are given by the energy changes
. The probability
governs the statistics of energy at times
and
, and we make no assumptions about the exact definition of
aside from assuming it is normalised and non-negative. As with the standard Leggett-Garg experiment described above, the aim is to compare the work statistics observed within different time intervals along the driving process, as shown in
Figure 1. Thus in analogy with Equation (
1) we will consider the following quantity;
This quantity can be measured over many runs of the driving process during each of the three time intervals, where the system is prepared in same state at time
for each experiment. For example
is obtained by measuring the average work done on the system sequentially during intervals
and
, and then subtracting the average work done during the total time interval
. We now seek to bound Equation (
3) through the assumption that the fluctuating work is a macrorealistic variable, as defined by (i) and (ii).
These assumptions imply the existence of a global probability distribution
describing the energy statistics at all points in time along the driving process, where each two-time distribution can be obtained as a marginal;
We can show that this condition immediately leads to the following Leggett-Garg inequality for the moments of work (see
Appendix A);
For
this implies
. This makes intuitive sense; in classical thermodynamics one would not expect to observe any difference between the sum of each intermediate average amount work done and the total average work done between the initial and final points in time. However, we will subsequently show that this does not generally hold for quantum systems, as the bounds in Equation (
5) can be violated for certain driving processes. It should also be noted that while we have assumed a discrete energy spectrum for the time-dependent Hamiltonian, this is not crucial to the derivation of Equation (
5). Indeed, so long as one assumes that the energy moments are always finite then inequalities of the form Equation (
5) can be derived. However, we will restrict our attention to a two-dimensional quantum system throughout the paper for simplicity, as this is sufficient to demonstrate violations of macrorealism in the statistics of fluctuating work.
3. Violations of the Leggett-Garg Inequalities for Work Moments
We will now utilise the inequalities Equation (
5) to show that fluctuating work can not generally be described by a macrorealist theory for quantum systems. The relevant situation that we consider is a standard setup for the thermodynamics of work extraction; an isolated system is initially thermalised and then driven out of equilibrium via changing its Hamiltonian in time, extracting work in the process [
15,
16,
17]. Suppose that we have a two-level system described by a time-dependent Hamiltonian
in the Heisenberg picture with initial state
such that;
where
and
are vectors. Note that the time-dependence of the Hamiltonian is attributed only to each vector
. In an experimental setup Equation (
6) describes a spin-1/2 particle coupled to an external classical magnetic field, with the direction of the applied field adjusted by the experimenter in time. Without loss of generality we will set the initial Hamiltonian along the z-axis of the Bloch sphere, i.e.,
, and choose an initial thermal state with respect to
at inverse temperature
;
. This in turn implies that
. To obtain the moments of work in Equation (
3) two projective energy measurements are performed at the start and end of the driving process within the fixed time intervals shown in
Figure 1. This method is commonly referred to as the two-projective measurement scheme (TPM) [
15,
16,
29]. The joint probability to observe energy
at
and then
at
is given by
Here we have denoted
as the projector onto the relevant energy state of the Hamiltonian at time
. Substituting Equation (
7) into Equation (
2) gives the moments of work from successive projective energy measurements (see
Appendix B):
To identify the conditions under which fluctuations in work violate Equation (
5), we substitute the above expression into Equation (
3) and parameterise the driving process by introducing
. For even
k we have the following quantum bound for
after minimising over all normalised vectors
for
(see
Appendix B);
which is saturated by choosing
. Secondly, for odd
k we find that maximising over
gives
where the bound is saturated by choosing
. From Equations (
9) and (
10) it is now apparent that there exists unitary protocols that violate the bounds in Equation (
5). Given that the bounds Equation (
5) necessarily follow from assumptions (i) and (ii), we conclude that quantum fluctuating work generally lacks a macrorealistic description. Let us note that Equation (
5) cannot be violated for odd
k in the high temperature limit. However, Equation (
8) shows that the even work moments are independent of
, and so Equation (
5) can indeed be violated for even
k regardless of temperature. This is not surprising, as the standard Leggett-Garg inequality can be violated in the case of an initially maximally mixed qubit [
4]. Violations of the work Leggett-Garg inequality, Equation (
5), can still occur at all temperatures since the system may acquire coherences in energy at intermediate times due to the unitary driving.
4. Inequalities for the Characteristic Function of Work
While the inequalities Equation (
5) provide a simple identification of non-classicality for the moments of work, we will show in this section that it is possible to condense this information into two inequalities related to the characteristic function for work rather than the moments themselves. A similar method has previously been applied to investigate the non-classical properties of electron-transport through conductors [
30]. The benefit of this approach is two-fold. While the projective energy measurements used to obtain the moments in Equation (
2) may be difficult to implement in practice, measurements of the characteristic function for work can be performed via ancilla-assisted measurement as shown in [
31,
32,
33,
34]. Secondly, inequalities for the characteristic function allow us to consider alternative non-invasive measurement schemes such as the full-counting statistics approach to the quantum work distribution proposed in [
24], as we later show in the next section.
The characteristic function uniquely defines the probability distribution for work in a unitarily driven system, and is obtained through the Fourier transform of the work distribution;
with work values
. Assuming the same protocol given by Equation (
6), we can consider a linear combination of characteristic functions for the three intervals of time shown in
Figure 1;
The assumptions (i) and (ii) for macrorealism imply the following upper bound on the real part of Equation (
12) for the qubit system (see
Appendix C);
whilst the imaginary part of Equation (
12) becomes an equality;
The bounds Equations (
13) and (
14) now constitute a pair of Leggett-Garg-type inequalities for any classical characteristic function for fluctuating work, assuming that the energy is given by either
or
at all times along the driving protocol.
Within the TPM scheme the characteristic function for work is given by [
15,
16];
where
is the initial state
decohered in the basis of
. The quantum upper bound for the real part of Equation (
12) is as follows (see
Appendix D);
Consequently the upper bound exceeds the classical inequality Equation (
13) for all
. Secondly, the upper bound for the imaginary part of Equation (
12) is
with a maximum violation of
for
. As one would expect, the quantum bounds Equations (
16) and (
17) are obtained by choosing the same protocol-dependent parameters used to obtain Equations (
9) and (
10) respectively.
To summarise this section we have presented a Leggett-Garg inequality for the characteristic function of fluctuating work, and shown that the non-classicality observed in the moments of work is also exhibited in the characteristic function itself for the TPM scheme.
5. Generalisation to Weak Measurements of Work
In the previous section we simply recast the original violations of the Leggett-Garg inequality Equation (
5) into the form relevant to the characteristic function for the work statistics of a qubit. While this was applied to the TPM protocol, the inequalities Equations (
13) and (
14) apply to any qubit with a fixed energy spectrum, and this ultimately allows us to investigate alternative measurement schemes and their resulting non-classical violations. Due to the invasive nature of projective measurements, some have argued that the definition of fluctuating work is thermodynamically inconsistent when applied to states with initial coherences in energy [
21,
23,
24,
25]. In turn this has inspired formulations of non-invasive work statistics that remain consistent with energy conservation in closed systems [
22,
23,
24]. In particular, we will consider the full-counting statistics for fluctuating work [
24,
26,
30,
35,
36,
37]. To obtain the full-counting statistics for work, one couples the system’s Hamiltonian to the momentum of an external detector, and subsequently measures the phase change acquired by the detector’s momentum during the driving process given by the system’s time-dependent Hamiltonian
[
36]. In turn this allows one to reconstruct a characteristic function of the following form [
24,
38];
For example, as shown in [
24] the first moment of Equation (
18) yields
; in other words the average work done on the system is simply the difference in average energy evaluated in the Heisenberg picture. Consequently one can easily see that the first moment in this measurement scheme still obeys the Leggett-Garg inequality Equation (
5). However, we demonstrate that the same is not true for the higher order moments contained in the characteristic function Equation (
18). Considering again the quantum two-level system described by Equation (
6), it can be shown that the real part of Equation (
18) is equal to the real part of Equation (15) obtained via the TPM scheme (see
Appendix E);
This is a surprising result, as it suggests that the same violations of Equation (
13), i.e., the upper bound Equation (
16), can be obtained non-invasively. It is only the imaginary part of
that differs from the TPM scheme, in which we find the following upper bound for Equation (
12);
which is again obtained by choosing
.
The full-counting statistics are not the only way to characterise non-invasive measurements of work. As proposed by Allahverdyan [
23], an alternative characteristic function describing the statistics of work derived from the Margenau-Hill distribution for successive energy measurements [
39,
40] is as follows (see
Appendix F):
where
denotes the symmetric Jordan product. Notably the corresponding probability distribution can be obtained via sequential weak measurement [
41]. While the first and second moments,
and
, are the same as those obtained from the full-counting statistics, in general higher order moments differ. However, for the isolated driven qubit we again find precisely the same violations of Equation (
13) because the real part of
is also equivalent to the real part of
. On the other hand the inequality for the imaginary term, Equation (
14), cannot be violated in our setup (
Appendix F).