## Quantum Non-Demolition Measurements and Leggett-Garg inequality

Paolo Solinas 1,2 and Stefano Gherardini 3,4

1 Dipartimento di Fisica, Universit` a di Genova, via Dodecaneso 33, I-16146, Genova, Italy.

2 INFN - Sezione di Genova, via Dodecaneso 33, I-16146, Genova, Italy.

3 Istituto Nazionale di Ottica del Consiglio Nazionale delle Ricerche (CNR-INO), Largo Enrico Fermi 6, I-50125 Firenze, Italy. 4 European Laboratory for Non-linear Spectroscopy, Universit´ a di Firenze, I-50019 Sesto Fiorentino, Italy.

Quantum non-demolition measurements define a non-invasive protocol to extract information from a quantum system that we aim to monitor. They exploit an additional quantum system that is sequentially coupled to the system. Eventually, by measuring the additional system, we can extract information about temporal correlations developed by the quantum system dynamics with respect to a given observable. This protocol leads to a quasi-probability distribution for the measured observable outcomes, which can be negative. We prove that the presence of these negative regions is a necessary and su ffi cient condition for the violation of macrorealism. This is a much stronger condition than the violation of the Leggett-Garg inequalities commonly used for the same task. Indeed, we show that there are situations in which Leggett-Garg inequalities are satisfied even if the macrorealism condition is violated. As a consequence, the quantum non-demolition protocol is a privileged tool to identify with certainty the quantum behavior of a system. As such, it has a vast number of applications to di ff erent fields from the certification of quantumness to the study of the quantum-to-classical transition.

## I. INTRODUCTION

How can we distinguish a quantum system from a classical one? What are the distinctive features of a quantum system? These are two of the questions physicists have asked from the dawn of quantum theory. To give a clear and definite answer is di ffi cult because of the smooth boundary between the classical and quantum worlds, and due to our inability to directly assess quantum features without altering them. For example, when a quantum system interacts with an environment, it quickly loses its quantum features and starts to behave classically. Still, the identification of protocols to spot and quantify a system's quantum features is of paramount importance and it would have a deep impact on both the foundation of quantum mechanics and quantum technologies.

ments cannot a ff ect the present state. All these assumptions are surely satisfied by a classical system, while they might be violated for quantum systems. Therefore, they give a testable way to distinguish between quantum and classical behavior. In this regard, it would be thus desirable to identify a condition that certifies the violation of macrorealism, eventually with certainty.

The milestones in this field are the papers by John Bell in 1964 [1] and by Leggett and Garg [2] in 1985. They both try to answer our initial questions, identifying and giving some criteria to highlight two quantum features. The famous Bell's inequalities [1, 3-5] focus on the non-locality and the quantum correlations related to the entanglement while the Leggett-Garg's inequalities (LGIs), which formally have a similar structure, discuss the violation of the macrorealism (MR) [2].

The concept of MR was introduced to identify the features that a classical system should have and then test if they are present at a quantum level. If not, one would conclude that a quantum system behaves in a non-classical way. In the following, for the sake of continuity with the current literature, we first use the original terminology of MR [2], which is still adopted in several contributions, but then we will give a modern and precise definition [6, 7]. The Leggett and Garg's assumptions to test the MR are 1) Macrorealism per se (MRps), whereby the system is in one of the states available to it at each moment, and 2) Noninvasive measurability (NIM) that identifies the conditions under which one can determine the state of a system without disturbing the subsequent dynamics.

Later, an additional condition (implicit in the original article [2]) was identified [8] as 3) Induction: future measure-

The LGIs are usually stated in terms of quantum correlators of sequential measurements performed on a single quantum system. The violation (resp. validation) of the LGIs is a signature of the quantum (resp. classical) behavior of the system of interest. However, the violation of LGIs is only a su ffi cient condition for the violation of MR [8], which is equivalent to stating that the validity of LGIs is a necessary condition such that even the MR is fulfilled. Practically, this means that there could be situations in which the MR is violated but the LGIs are satisfied.

In this paper, we focus on the MRps condition, comprising the MR, by giving a new general protocol guaranteeing both a necessary and su ffi cient condition to identify its violation. The violation of the MRps implies the violation of the MR's assumption that, in turn, identifies the presence of quantum features. Such a protocol is based on quantum non-demolition (ND) measurements [9-12] and it was proposed in Refs. [1318]. To highlight the practical implications and experimental feasibility of the protocol, we will start by describing its implementation which makes use of a quantum detector to store the desired information. In this way, by measuring the state of the detector, one can construct a quasi-probability distribution (QPD) of the corresponding measurement outcomes. As clarified in Refs. [18, 19] other quasi-probability distributions have been defined in the literature, ranging from KirkwoodDirac quasi-probabilities [19-26] to the full counting statistics [27-30] and Keldysh quasi-probabilities [31, 32]. At the same time, alternative protocols to determine the violation of the LGIs with weak measurements have been proposed [3335]. Nevertheless, in this paper, our focus goes to quasiprobabilities based on ND measurements, as we are going to prove that the corresponding distribution has negative regions if and

only if the MRps condition is violated.

In order to show the e ff ectiveness of our results, we compare the ability of our protocol to reveal a violation of MR with the LGIs, representing so far the proper tool for such a task. For the sake of a fair comparison, we derive, from the full QPD, the multi-time quantum correlators leading to the corresponding LGIs. With an in-depth analysis, we identify the contributions of the QPD that lead to the violation of the LGIs and we show that only in some particular cases the LGIs spot the violation of the MR, thus confirming that such a circumstance is only a su ffi cient condition. To further stress this point and test on a concrete case study, we consider a specific example taken from Ref. [8]. We show how, by using the negative regions of the QPD, we are always able to identify any violation of MRps and thus of MR, while LGIs statistically fail in half of the cases.

To the best of our knowledge, the protocol based on ND presented here is the first allowing for the certain identification of purely quantum behaviors associated with the violation of MR. This feature joined with the fact that it finds application beyond binary observables, for which usually the LGIs are discussed, has implications for the foundation of quantum mechanics and quantum technologies. For example, it can be used (i) to characterize the quantum-to-classical transition [36-39], (ii) to monitor the stability and robustness against noise of quantum devices and computers subject to decoherence process [16, 40, 41], or (iii) to certify quantum random number generator [42]. Such a reliable and quantitative tool would not only enhance our comprehension of quantum phenomena but also inspire innovative methods for preserving the quantum properties essential to technological advancements.

## II. MACROREALISM CONDITIONS

The original definition of MRps given by Leggett and Garg [2] was ambiguous and prone to di ff erent interpretations. During the years, this led to confusion about its meanings and implications. Only recently some authors have clarified that it was intended to understand if a quantum system can be described in terms of classical random variables [6, 7]. More precisely, once the observable to be measured is identified, the MRps assumption is satisfied if quantum superpositions are not allowed, at any time, within the Hilbert space spanned by the eigenbasis of the observable or by a statistical mixture of them [7]. In other words, the MRps assumption is satisfied if there cannot be coherent superpositions and the state of the system is always diagonal in the basis that diagonalizes the observable to be measured [7, 43].

The NIM assumption, entering MR, is equally di ffi cult to formalize. For this reason, it is usually flanked by the nosignaling in time (NSIT) assumption [8, 44], which is regarded as a statistical (and slightly less strong) version of NIM. As extensively discussed by Halliwell [8], this allows us to relate the NIM with the possibility of attaining the marginal probability of single events from the joint probability of the sequence of them.

That is, as an example, if p s ( 1 , s 2 , s 3 ) is a generic joint probability of recording s 1 , s 2 and s 3 , the NSIT condition over, say, the outcome s 3 , reads: P s 1 , s 2 p s ( 1 , s 2 , s 3 ) = P s ( 3 ), where P s ( 3 ) is the probability to measure s 3 in correspondence of the state of the system at the time the measurement is performed, ensuring that no previous measurements in the sequence have a ff ected it. For the sake of clarity, throughout the manuscript, the symbol P will always denote probabilities at single times, given by Born's rule, and probabilities at multi-times respecting the theory of classical probability.

Below, we show that the ND quasi-probability always satisfies the NSIT assumption so that the violation of MR can be reduced only to the violation of the MRps' assumption, for which we give both a necessary and su ffi cient condition.

## III. THREE MEASUREMENT QUASI-PROBABILITY DISTRIBUTION

We consider the following scheme for sequential nondemolition measurements. Suppose we have a quantum system S evolving under a unitary transformation in a time interval 0 ≤ t ≤ T . At times t 0 ≤ t 1 ≤ t 2 = T , we measure a generic observable ˆ A (Hermitian operator).

In general, any projective measurement perturbs a quantum system inducing the collapse of its wave function, thus destroying the quantum coherence between eigenstates of the measured observable. To avoid this drawback in a scheme with sequential measurements, according to the ND approach, we use an auxiliary quantum detector to store the desired information in the phase of the detector's state, which is eventually measured. This scheme allows us to preserve the quantum coherence in the initial density operator and obtain the average value of observables at multi-times, not perturbed by the interaction with the detector [14-17, 45].

Sequential non-demolition measurements rely on a sequence of fast (with respect to the timescale of the system evolution) system-detector couplings. The unitary transformations describing the coupling processes are ˆ uk = exp { i ( λ/ 2) ˆ A ⊗ ˆ p } , which occurs at any time t k with k = 0 1 2 (see Refs. [13-, , 15] for more details). Here, λ is an e ff ective coupling, and ˆ p is an operator acting on the detector's degrees of freedom.

We denote the unitary transformation from time t j -1 to t j acting on the system only as ˆ Uj ⊗ I = ˆ ( U tj , t j -1 ) ⊗ I (with j = 1 2) so that the total unitary transformation acting on , both the system and detector is ˆ Utot = ˆ u U u U u 2 ˆ 2 ˆ 1 ˆ 1 ˆ 0 . In the following, to simplify the notation we will denote ˆ Ui ≡ ˆ Ui ⊗ I . Explicitly, the total evolution operator ˆ Utot reads as

$$\hat { U } _ { t o t } = e ^ { i \frac { 4 } { 2 } \hat { A } \otimes \hat { p } } \hat { U } _ { 2 } \, e ^ { i \frac { 4 } { 2 } \hat { A } \otimes \hat { p } } \hat { U } _ { 1 } \, e ^ { i \frac { 4 } { 2 } \hat { A } \otimes \hat { p } }. \quad \quad ( 1 )$$

Notice that, since in general [ ˆ Uj , ˆ ] A , 0, this scheme can also be applied to describe the sequential measurement of three distinct non-commuting operators ˆ , A B ˆ and ˆ . C Observe also that this ND protocol can be implemented in any quantum platform that allows for a tunable interaction between the quantum system of interest and an auxiliary system acting as the detector. For example, in Ref. [16], the scheme was im-

FIG. 1: A pictorial representation of the evolution in terms of sequential measurements at three times, and Feynman paths built over the measurement outcomes. The measurements occur at time t 0 , t 1 and t 2 . The possible outcomes of the measurements are ai and al at time t 0 , aj and am at time t 1 , and ak at time t 2 . The sequence of outcomes identifies the paths that the monitored system could follow during its evolution. The red dashed curve represents a macrorealistic (or classical) path where, at any time, we have a single outcome (e.g., ai → am → ak ) so that the observable has a determinate outcome. The blue curves present quantum paths in which the system is in a superposition of states with di ff erent outcomes. In this case, the observable ˆ A does not have a unique value and the paths violated the macrorealism requirement.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00088v1-figure-d16dd0b0-c44d-438c-9514-30f8c2449bec.jpg)

plemented using several qubits to describe the system and one qubit for the detector.

We decompose the initial state of the system in the basis {| i ⟩} in which the observable ˆ A is diagonal, i.e. ˆ A i | ⟩ = ai | i ⟩ with ai scalars, so that | ψ 0 ⟩ = P i ψ 0 i | i ⟩ . Throughout the paper, we assume that | ψ 0 ⟩ is not an eigenstate of ˆ . A On the other hand, the state of the detector is taken to be | ϕ 0 ⟩ = GLYPH&lt;16&gt; 1 / √ M GLYPH&lt;17&gt; P p | p ⟩ , where ˆ p p | ⟩ = p p | ⟩ . The calculation of the evolution of the system and detector together is a straightforward extension of the one presented in Refs. [15, 46], and the final total state is

$$| \Psi \rangle = \frac { 1 } { \sqrt { M } } \sum _ { p } \sum _ { i, j, k } e ^ { i \frac { a p } { 2 } ( a _ { k } + a _ { j } + a _ { i } ) } U _ { 2, k j } U _ { 1, j i } \psi _ { i } ^ { 0 } \left | k \right \rangle \left | p \right \rangle \ \ ( 2 )$$

where U kj 2 , = ⟨ k U | ˆ 2 | j ⟩ and U 1 , ji = ⟨ j U | ˆ 1 | i ⟩ .

The information on the evolution of the monitored system is extracted by measuring the phase accumulated in the detector between two generic states |± p ⟩ . This is obtained by taking the o ff -diagonal element of the detector's density operator and tracing out the degrees of freedom of the system. Formally, writing the final total density operator as ˆ R = Ψ Ψ | ⟩ ⟨ | and the density operator of the detector as ˆ r 0 = | ϕ 0 ⟩ ⟨ ϕ 0 | , this is G λ = Tr GLYPH&lt;2&gt; ( I ⊗ ⟨ p | ) R ( I ⊗ | -p ⟩ ) GLYPH&lt;3&gt; / ⟨ p r | 0 |-p ⟩ .

Normalizing p to one (in some opportune unit), G λ explicitly reads as

$$\mathcal { G } _ { \lambda } = \sum _ { i, l, j, m, k } e ^ { i \lambda \left ( a _ { k } + \frac { a _ { j } + a _ { m } + a _ { i } + a _ { l } } { 2 } \right ) } U _ { 2, k j } U _ { 1, j i } \rho _ { i l } ^ { 0 } U _ { 1, l m } ^ { * } U _ { 2, m k } ^ { * } \quad ( 3 ) \quad \pm 2 \left ( \mathfrak { f } _ { i } \right )$$

where ˆ 0 ρ = | ψ 0 ⟩ ⟨ ψ 0 | is the density operator of the system at the initial time t 0 = 0, and ρ 0 il = ⟨ | i ˆ ρ 0 | l ⟩ . The function G λ is also called quasi-characteristic function [14-17, 45].

As usual, the ND quasi-probability distribution is obtained by calculating the inverse Fourier transform of G λ , i.e., P ND ( ∆ = ) (2 π ) -1 R exp {-i λ ∆ }G λ d λ . If we introduce

$$\Delta _ { k, j, m, i, l } = a _ { k } + \frac { a _ { j } + a _ { m } + a _ { i } + a _ { l } } { 2 } \, \quad \ \ ( 4 )$$

and we use the projectors ˆ Π = k | k ⟩ ⟨ k | (with ˆ 2 Π = Π k ˆ k , P k ˆ Π = k I and ˆ Π Π k ˆ j = δ k j ˆ Π k ), P ND ( ∆ ) is

$$\mathcal { P } _ { N D } ( \Delta ) = \sum _ { i, l, j, m, k } P _ { N D } ( k, j, m, i, l ) \delta \left [ \Delta - \Delta _ { k, j, m, i, l } \right ] \quad ( 5 )$$

where PND k ( , j , m i l , , ) = U kjU 2 , 1 , ji ρ 0 il U ∗ 1 , lm U ∗ 2 , mk = Tr h ˆ Π k ˆ U 2 ˆ Π j ˆ U 1 ˆ Π i ˆ ρ 0 ˆ Π l ˆ U † 1 ˆ Π mU ˆ † 2 i , and δ [ ] · denotes the Dirac's delta. Thus, PND k ( , j , m i l , , ) is the probability amplitude to arrive in the state | k ⟩ following the superposition of paths i → j → k and l → m → k (see Fig. 1).

The ND quasi-probability distribution P ND ( ∆ ) comprises both a classical and a quantum contribution, P cl ( ∆ ) and P q ( ∆ ), respectively. The former contribution is the one in which a well-defined value of ˆ A is recorded at any measurement time by measuring the system. That is, for m = j and l = i , so that

$$\mathcal { P } _ { \text{cl} } ( \Delta ) = \sum _ { i, j, k } P _ { N D } ( k, j, j, i, i ) \delta \left [ \Delta - \Delta _ { k, j, j, i, i } \right ]. \quad ( 6 )$$

These can be interpreted as probabilities obeying classical probability theory since PND k ( , j , j , i , i ) = P i (0) P i (1) → j P j (2) → k where P i (0) = ρ 0 ii , P i (1) → j = | U 1 , ji | 2 and P j (2) → k = | U kj 2 , | 2 are the probabilities of finding the system initially in the state | i ⟩ , of having a transition from state i to j and from j to k , respectively.

The quantum contributions in Eq. (5), thus with m , j or l , i , are associated with quantum trajectories where the system is in a coherent superposition of paths defined by the measurement outcomes; see Fig. 1. More specifically, at time t 0 , the measurement outcomes might be ai or al while, at time t 1 , aj or am . Formally, P q ( ∆ ) is

$$\mathcal { P } _ { q } ( \Delta ) = \sum _ { i, l, j, m, k } ^ { \prime } P _ { N D } ( k, j, m, i, l ) \delta \left [ \Delta - \Delta _ { k, j, m, i, l } \right ] \quad ( 7 )$$

where the sum X ′ does not include the terms for which m = j and l = i entering the classical distribution P cl ( ∆ ); this means that in P q ( ∆ ) are allowed contributions with ( l = i , m , j ), ( m = j , l , i ) and ( m , j , l , i ).

We stress that the violation of MRps in the P q ( ∆ ) shall be manifest also in the presence of classically forbidden values of ∆ k j m i l , , , , [14-17, 45], as given in Eq. (4). For example, suppose that ai = ± 1. The possible classical (i.e., m = j and l = i ) values of ∆ are {-3 , -1 1 3 . , , } On the contrary, the quasidistribution P q ( ∆ ) also accounts for terms associated with ∆ = ± 2 (for example, with ak = aj = am = ai = 1 and al = -1) or ∆ = 0 (for example, with ak = aj = 1 and am = ai = al = -1).

The QPD P ND ( ∆ ) in Eq. (5) satisfies the noninvasive measurability condition, in terms of the no-signaling in time assumption [2, 8, 47-49]. This means that marginalizing P ND ( ∆ ) over all the measurement outcomes but the ones at a single time returns the corresponding probability to measure the observable ˆ given by the Born's rule. For example, to obA tain the probability of recording a final eigenvalue ak , we must sum over all the intermediate outcomes ai , al , aj , am . It can be shown (see Methods) that P ak ( ) = P i l , , j m , PND k ( , j , m i l , , ) = Tr h ˆ Π k ˆ U ˆ ρ 0 ˆ U † i . Indeed, P ak ( ) is the probability to measure ˆ A at time T , when the system has evolved under the unitary transformation ˆ U [8]. Similarly, the probability to record the eigenvalues ai and al at time t = 0, is obtained from Eq. (5) by summing over the indices l , j , m k , and i , j , m k , respectively. Thus, we have P ai ( ) = P l , j m k , , PND k ( , j , m i l , , ) = Tr h ˆ Π i ˆ ρ 0 i as expected. With similar calculations, it can be shown that also the intermediate measurements satisfy the NIM condition, such that P aj ( ) = P i l m k , , , PND k ( , j , m i l , , ) = Tr h ˆ Π j ˆ U 1 ˆ ρ 0 ˆ U † 1 i and P am ( ) = P i l , , j k , PND k ( , j , m i l , , ) = Tr h ˆ Π mU ˆ 1 ˆ ρ 0 ˆ U † 1 i . It is also worth observing that summing PND k ( , j , m i l , , ) over three indices returns the Kirkwood-Dirac quasi-probabilities at two-times [18, 19, 25], apart the 3-tuples ( j , m k , ) and ( i , l , k ) that give rise to a single-time probability. By construction, Kirkwood-Dirac quasiprobabilities respect the assumption of no-signaling in time.

## IV. NEGATIVE REGIONS ARE NECESSARY AND SUFFICIENT FOR VIOLATING MACROREALISM

Wenowprove that P ND ( ∆ ) has negative contributions if and only if the macrorealism, i.e., MRps, condition is violated. The proof is composed of three steps. First, we show that P ND ( ∆ ) is real, i.e. given by a distribution of real numbers. Then, we prove that it is normalized to 1 and that this normalization comes uniquely from the classical contribution P cl ( ∆ ). This implies that the terms entering the quantum contribution P q ( ∆ ) cancel out each other and that some of them must be negative. The last part of the proof involves the connection between the macrorelism and the negativity of P ND ( ∆ ).

To prove that P ND ( ∆ ) is real, we group the terms in the distribution that are multiplied by the same δ -function δ [ ∆ -∆ k j m i l , , , , ]. There are indeed identical values of ∆ k j m i l , , , , because of the symmetries between the exchange of the indices ( i , l ) and ( j , m ). Formally, using the properties of the matrix elements, e.g., U km 2 , = ⟨ k U | ˆ 2 | m ⟩ = GLYPH&lt;16&gt; ⟨ mU k | ˆ 2 | ⟩ GLYPH&lt;17&gt; ∗ , we have PND k m ( , , j , i , l ) = ( PND k ( , j , m i l , , )) ∗ so that Eq. (7) can be rewritten as

$$\mathcal { P } _ { q } ( \Delta ) = 2 \sum _ { i, l, j, m, k } ^ { \prime \prime } \text{Re} \left [ P _ { N D } ( k, j, m, i, l ) \right ] \delta \left [ \Delta - \Delta _ { k, j, m, i, l } \right ] \quad ( 8 )$$

where X ′′ sum over half the terms than X ′ as it includes terms for which l &gt; i or m &gt; j . Since the distribution P cl ( ∆ ) contains (positive) classical joint probabilities by construction and P q ( ∆ ) comprises real numbers, we conclude that P ND ( ∆ ) is real.

The proof that P cl ( ∆ ) is normalized comes from noting that its contributions are positive joint probabilities obeying the classical probability theory. In fact, we have that R d ∆ P cl ( ∆ ) = P i , j k , P i (0) P i (1) → j P j (2) → k = P i , j k , Tr h ˆ Π k ˆ U 2 ˆ Π j ˆ U 1 ˆ Π i ˆ ρ 0 ˆ Π i ˆ U † 1 ˆ Π j ˆ U † 2 i = 1.

Regarding the integration of the quantum contribution P q ( ∆ ) of the ND quasi-distribution P ND ( ∆ ), it is worth distinguishing two distinct sums of terms: one with m , j , and the other with ( m = j , l , i ). Then, as above, we integrate both sums of terms in P q ( ∆ ) over ∆ . The integration over ∆ of the sum with m , j leads to

$$\mathfrak { z } _ { \overline { y }. } \quad 2 \sum _ { k, i, l, j, m > j } \text{Re } \text{Tr} \left [ \hat { \Pi } _ { k } \hat { U } _ { 2 } \hat { \Pi } _ { j } \hat { U } _ { 1 } \hat { \Pi } _ { i } \hat { \rho } ^ { 0 } \hat { \Pi } _ { l } \hat { U } _ { 1 } ^ { \dagger } \hat { \Pi } _ { m } \hat { U } _ { 2 } ^ { \dagger } \right ] = 0, \quad ( 9 )$$

where we have used the projectors' properties ˆ Π Π k ˆ j = δ , k j ˆ Π k = 0 and P k ˆ Π k = I . An analogous result is obtained from integrating over ∆ the sum of terms in P q ( ∆ ) with ( m = j , l , i ). In fact, one can get that

$$\mathbf r _ { \mathbf s } \int _ { k, i, j, l > i } \mathbf R e \, \mathbf r r \left [ \hat { \Pi } _ { k } \hat { U } _ { 2 } \hat { \Pi } _ { j } \hat { U } _ { 1 } \hat { \Pi } _ { i } \hat { \rho } ^ { 0 } \hat { \Pi } _ { l } \hat { U } _ { 1 } ^ { \dagger } \hat { \Pi } _ { j } \hat { U } _ { 2 } ^ { \dagger } \right ] = 0 \,. \quad ( 1 0 )$$

Therefore, the normalization of P ND ( ∆ ) uniquely comes from the classical distribution P cl ( ∆ ). Given that the terms in P q ( ∆ ) are real numbers and cancel out upon integration, it follows that at least some of them must be negative. Hence, the negativity of P ND ( ∆ ) derives from contributions in P q ( ∆ ) that are built with the coherent superposition of ˆ 's A eigenstates at times t 0 , t 1 , t 2 . We can thus state that necessary condition for the negativity of P ND ( ∆ ) is the presence of contributions in P q ( ∆ ) coming from the coherent superposition of eigenstates associated to an observable that is measured at multiple times .

Notably, using ND quasi-probabilities, the statement above is also a su ffi cient condition , since the negativity of P ND ( ∆ ) can only derive from P q ( ∆ ). Since the two summations in (9)(10) must always be zero, whenever there is a coherent superposition of the measurement observable's eigenstates, i.e., P q ( ∆ ) , 0, P ND ( ∆ ) must contain negative (real) terms. That is, a su ffi cient condition for the negativity of P ND ( ∆ ) is the presence of the coherent superposition of eigenstates associated with the observable measured at multiple times .

We conclude by noting that if a system is in a coherent superposition of the eigenstates of an observable measured at distinct times, then the MRps assumption is violated. Accordingly, the MRps assumption is violated if and only if the ND quasi-probability distribution P ND ( ∆ ) exhibits negativity .

## V. LEGGETT-GARG INEQUALITIES

Having established that negative regions in the ND quasiprobability distribution are necessary and su ffi cient indicators of the violation of MRps, we now compare our results with the usual tool employed so far to identify such a quantum feature: the Leggett-Garg inequalities.

The LGIs are usually formulated for the so-called binary observables [2, 8, 47, 49], that is, for observables that

can have only two outcomes. Hence, while in the previous discussions, the observable ˆ A had a generic discrete spectrum, here we assume that its eigenvalues can be only ai = ± 1 for any time t i . LGIs can be built using quantum correlators at two times, say t i and t j , which are defined as [8, 47, 49] Cij = Tr hGLYPH&lt;16&gt; ˆ ( A ti ) ˆ ( A t j ) + ˆ ( A t j ) ˆ ( A ti ) GLYPH&lt;17&gt; ˆ ρ 0 i / 2, where ˆ ( A ti ) = ˆ U † ( t i , 0) ˆ ˆ ( AU ti , 0). The expression of Cij can be rewritten as [8, 47, 49]

$$C _ { i j } = \sum _ { i, j } a _ { i } a _ { j } P ( a _ { j }, a _ { i } ) = \sum _ { i, j, k } a _ { i } a _ { j } P ( a _ { k }, a _ { j }, a _ { i } ) \quad ( 1 1 ) \quad \text{ with } \quad \text{ with } \quad \text{ }$$

where, by construction, P aj ( , ai ) = P k P ak ( , aj , ai ) and, as above, P ak ( , aj , ai ) = P i (0) P i (1) → j P j (2) → k is the joint probability to record the measurement outcomes ak , aj and ai at times t i , t j , t k by means of a sequential scheme. The sequential measurements are performed at three times t 0 = 0, t 1 and t 2 .

Following Ref. [49], we can calculate the Leggett-Garg (LG) parameter

$$K \equiv C _ { 0 1 } + C _ { 1 2 } - C _ { 0 2 }. \text{ \quad \ \ } ( 1 2 )$$

Notably, for any classical system fulfilling MR, the LG parameter K reads as

$$K = 1 - 4 \left [ P ( 1, - 1, 1 ) + P ( - 1, 1, - 1 ) \right ] \quad \ \ ( 1 3 )$$

that results in the LGI -3 ≤ K ≤ 1. That is, if the -3 ≤ K ≤ 1 is not satisfied, then at least a condition of MR is violated.

In this context, it would be desirable to find a decomposition of K made of two terms: one equivalent to the right-handside of Eq. (13) (valid under MR), and the other di ff erent from zero whenever the LGI -3 ≤ K ≤ 1 is violated. For this purpose, we can still resort to the ND quasiprobabilities given that, for binary observables (here, the eigenvalues are ± 1), the correlators Cij can be equivalently expressed as [8]

$$C _ { i j } = \sum _ { i, j } a _ { i } a _ { j } P _ { N D } ( j, i ), \quad \ \ \ ( 1 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where PND ( j , i ) = P l m k , , [ PND k ( , j , m i l , , ) + PND k m ( , , j , i , l ) + PND k ( , j , m l i , , ) + PND k m ( , , j , l , i )] with PND ( j , i ) given by Kirkwood-Dirac quasiprobabilities as commented before. Accordingly,

$$C _ { 0 1 } = \frac { 1 } { 4 } & \sum _ { i, j } a _ { i } a _ { j } \sum _ { l, m, k } [ P _ { N D } ( k, j, m, i, l ) + P _ { N D } ( k, m, j, i, l ) \quad \text{$As$} \\ & + P _ { N D } ( k, j, m, l, i ) + P _ { N D } ( k, m, j, l, i ) ]. \quad ( 1 5 ) \quad \text{$current$}$$

Writing the correlator at two-times C 01 as a function of the three-time ND quasiprobabilities brings the advantage of including also the superposition of the observable ˆ 's A eigenstates resulting in multiple di ff erent possible outcomes: ai and aj at times t 0 = 0, and aj and am at time t 1 .

The definition of the correlators C 12 and C 02 involving the (final) measurement at time t 2 are similar to that of C 01 : C 12 = P j k , akajPND ( k , j ) and C 02 = P i k , akaiPND ( k i , ), whose extended expression is in the Methods.

Since we are working with binary variables, each of them can take two values; in our case, ai and -ai . For simplicity, we refer to them with the indices i and ¯ ; i for example, PND k ( ¯ , j , m i i , , ¯ ) shall stand for PND ( -ak , aj , am , ai , -al ). Using the correlators Cij as defined above, the LG parameter K for the ND procedure reads as

$$K \ = \ \frac { 1 } { 4 } \sum _ { i, l, j, m, k } f ( k, j, m, i, l ) P _ { N D } ( k, j, m, i, l ) \quad ( 1 6 )$$

with f ( k , j , m i l , , ) = ( ai + al )( aj + am ) + ak ( aj + am -ai -al ); see Methods for details. It is convenient to separate in K the classical contributions when m = j and l = i , i.e.,

$$K _ { c l } = \frac { 1 } { 2 } \sum _ { k, j, i } \left [ a _ { i } a _ { j } + a _ { k } ( a _ { j } - a _ { i } ) \right ] P _ { N D } ( k, j, j, i, i ), \quad ( 1 7 )$$

from the quantum ones when m = j and l = ¯ i , i

$$K _ { q, 1 } = \frac { 1 } { 2 } \sum _ { i, j, k } a _ { k } a _ { j } P _ { N D } ( k, j, j, i, \bar { i } ) \quad \quad ( 1 8 )$$

and when m = ¯ j , j and l = i

$$K _ { q, 2 } = - \frac { 1 } { 2 } \sum _ { i, j, k } a _ { k } a _ { i } P _ { N D } ( k, j, \bar { j }, i, i ). \quad \quad ( 1 9 )$$

First, we now show that from the classical contributions in Eq. (16) we recover the classical results in the right-handside of Eq. (13). To do this, we recall that PND k ( , j , j , i , i ) = P i (0) P i (1) → j P j (2) → k with P i , j k , P i (0) P i (1) → j P j (2) → k = 1, namely the total probability over a complete set of classical paths must sum to 1. The latter can be written explicitly as P i k , PND k i i i i ( , , , , ) + P k PND k k k k k ( , , , ¯ , ¯ ) + P k PND k k k k k ( , ¯ , ¯ , , ) = 1 where we have split the contributions for j = i (and any k ), j = k and i = ¯ k , k , and j = ¯ k , k and i = k . Using this separation criterion to Eq. (17), the classical contribution Kcl can be written as

$$K _ { c l } = 1 - 4 \sum _ { k } P _ { N D } ( k, \bar { k }, \bar { k }, k, k ). \text{ \quad \ \ } ( 2 0 )$$

Expanding the sum for k = ± 1 leads to the result in Eq. (13) so that Kcl satisfies the LGI, i.e., -3 ≤ Kcl ≤ 1.

As shown in the Methods, the remaining terms in Eq. (16) can be simplified by noting that Kq , 1 = 0. Therefore, when present, the violation of the LGI is due to the second quantum contribution Kq , 2 in Eq. (19). Kq , 2 can be further simplified by observing that P k PND k ( , j , ¯ j , i , i ) = 0 so that PND k ( ¯ , j , ¯ j , i , i ) = -PND k ( , j , ¯ j , i , i ). In this way, setting ak = 1 without loss of generality, Eq. (19) equals to

$$K _ { q, 2 } = - 4 \sum _ { k } \text{Re} \left [ P _ { N D } ( k, j, \bar { j }, k, k ) \right ]. \quad \ \ ( 2 1 )$$

The derivation of Eq. (21) is in the Methods.

Some remarks are due to this point. Using the ND quasiprobability distribution highlights the connection between the

FIG. 2: The top-left panel shows the LG parameter K in Eq. (22) as a function of ωτ . The shaded regions are the ones for which the LGI is violated. The other panels show the quasi-probability distribution P ND ( ∆ ) evaluated at three di ff erent values of the ωτ . The P ND ( ∆ ) always shows negative regions, while only for the green circle point the LGI is violated.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00088v1-figure-fbd57fc4-4981-4227-bd5a-3e5a6642dd1f.jpg)

violation of macrorealism and LGIs. In particular, the contributions with m = ¯ j , j and l = i , denoting the superposition between eigenstates of ˆ , A are the ones responsible for the violation of the LGI. However, there are situations where the MRps condition is violated but the LGI is satisfied. This occurs for any initial condition for which Kq , 2 vanishes. In such a case, the total quantum contribution Kq vanishes as well, Kq = 0, with the result that K = Kcl and the LGI is satisfied. Nevertheless, we can have that l , i , which means that the system is initially in a superposition of eigenstates of ˆ A and this fact naturally violates the MRps. On the contrary, if we use the ND quasi-probability P ND ( ∆ ), the terms PND k ( , j , j , i , ¯ ) i in Eq. (18) always contribute to the negative regions of the quasi-probability distribution. This is a confirmation that the LGIs give only a su ffi cient condition for the violation of the MRps and thus of the MR's assumption, while the ND quasi-probability P ND ( ∆ ) gives both a su ffi cient and necessary condition.

As discussed above, for a classical system -3 ≤ K ≤ 1. Hence, any value of K outside this range implies the violation of such an inequality. If we consider the ND protocol, the whole dynamics of the system and the detector is given by Eq. (1), which leads to the ND quasi-probability distribution in Eq. (5).

## VI. EXAMPLE

Following Ref. [8], we discuss a specific example that clarifies the results obtained. Let us consider a two-level quantum system that is initially in the state | ψ 0 ⟩ = ( |↑⟩ + i |↓⟩ ) / √ 2 where |↑⟩ and |↓⟩ are the eigenstates of the Pauli operator ˆ σ z . The dynamics is generated by the Hamiltonian ˆ H = ωσ / ˆ x 2 and we measure the operator ˆ A = ˆ σ z at times t 0 = 0, t 1 = τ and t 2 = 2 . τ

The correlators Cij [see Eq. (11)] and the corresponding LG parameter K can be calculated analytically:

$$K = C _ { 0 1 } + C _ { 1 2 } - C _ { 0 2 } = 2 \cos ( \omega \tau ) - \cos ( 2 \omega \tau ). \quad ( 2 2 ) \quad \text{the c}$$

The numerical results for the example are shown in Fig. 2. The top-left panel shows the value of K as a function of ωτ with the shaded regions representing the situations in which the LGI -3 ≤ K ≤ 1 is violated. The colored shapes (dot, square, triangle) give the value of K for three specific values of ωτ . The other panels of the figure show the ND quasiprobability distribution [as obtained from Eq. (5)] that are associated with each colored shape. As we can see, the LGI is violated only for some parameter choices, e.g., the green dot in the top-left panel. However, for intermediate values of ωτ (the orange square and purple triangle), the LGI is satisfied but the dynamics exhibit quantum features given by the presence of quantum coherence along the eigenbases of ˆ . A As proved above, this circumstance is identified by determining negative regions of the corresponding distribution P ND ( ∆ ). As expected, the ND quasi-probability distribution is always able to identify the quantum features of the system even when the LGI fails.

From Eq. (22), the violation of the LGI -3 ≤ K ≤ 1 occurs for 0 ≤ ωτ ≤ π/ 2 and 3 π/ 2 ≤ ωτ ≤ 2 π . This means that the LGI correctly identifies the violation of MRps in only half of the cases, while the ND protocol always succeeds.

We have presented a protocol that allows us to identify the presence of genuinely quantum behaviors unambiguously, in terms of the violation of macrorealism per se (MRps). This protocol is based on performing sequential quantum nondemolition measurements of a given observable ˆ , A and this can be attained using a quantum detector. Using the corresponding measurement outcomes recorded at multiple times, we can construct a quasi-probability distribution. We have proven that the presence of negative regions in such a quasiprobability distribution is a necessary and su ffi cient condition for the violation of the MRps (and macrorealism) that outlines the presence of quantum behaviors.

The non-demolition protocol has some additional features that make it interesting for more practical applications. First, it can provide results overcoming some limitations imposed by the LGIs, which are usually employed for binary observables, as the protocol with ND measurements can be implemented for observables with arbitrary discrete spectra. Secondly, being based on an experimental protocol, the ND protocol gives an operational procedure that, in principle, could be implementable on any quantum platform [16].

The advantage of the ND protocol over the LGIs ultimately is a consequence of the more information we extract from the system and its dynamics. While using LGIs we get information only about correlators at two times, with the ND protocol we have access to a full quasi-probability distribution and, therefore, to all their moments [14, 15, 45]. While it is natural to think that the ND protocol gains more information about non-classicality and, thus, needs more resources than other methods as LGIs, the comparison in terms of resources is not straightforward. Regarding the ND quasi-probability distribution, it is obtained by a Fourier transform of the measured data. Therefore, the resources needed to implement the ND protocol depend critically on the requested precision in performing the Fourier transform. To evaluate this a more detailed analysis is needed.

The research fields that can benefit mostly from the implementation of the ND protocol as a tool to identify quantum behavior are the foundation of quantum mechanics and quantum technologies. A few examples are the generation of certified random numbers [42], the study of quantum gravity in mesoscopic systems, and the cooling of optomechanics systems up to a scale where quantum e ff ects play a role [39]. In all these cases, the question of when and how the quantum-to-classical transition occurs is of critical importance, and it might help to improve already-performed experimental results [36, 37] or to better identify causes leading to decoherence in open quantum systems, quantum computers or quantum annealers [16, 40, 41]. This could be possible by showing that a given observed quantum-to-classical transition corresponds to a reduction of negative regions in an ND quasi-probability distribution that is built over a properly chosen observable ˆ . A The results in Refs. [16, 17] seem to suggest that such a circumstance could be in principle attainable.

To compute the probability of measuring an outcome ak (eigenvalue of ˆ ) A from the measurement at time t 2 , we sum over all the intermediate outcomes ai , al , aj , am . In this way, using the properties of the projectors ( ˆ 2 Π = Π k ˆ k , P k ˆ Π = k I and ˆ Π Π k ˆ j = δ k j ˆ Π k ), we get

$$P ( a _ { k } ) \ & = \ \sum _ { i, l, j, m } P _ { N D } ( k, j, m, i, l ) \\ & = \ \sum _ { i, l, j, m } \text{Tr} \left [ \hat { \Pi } _ { k } \hat { U } _ { 2 } \hat { \Pi } _ { j } \hat { U } _ { 1 } \hat { \Pi } _ { i } \hat { \rho } ^ { 0 } \hat { \Pi } _ { l } \hat { U } _ { 1 } ^ { \dagger } \hat { \Pi } _ { m } \hat { U } _ { 2 } ^ { \dagger } \right ] \\ & = \ \text{Tr} \left [ \hat { \Pi } _ { k } \hat { U } \hat { \rho } ^ { 0 } \hat { U } ^ { \dagger } \right ] = \text{Tr} \left [ \hat { \Pi } _ { k } \hat { \rho } ( \mathcal { T } ) \right ]$$

with ˆ U = ˆ U U 2 ˆ 1 and ˆ( ρ T ) = ˆ U ˆ ρ 0 ˆ U † . P ak ( ) is the probability to record ak at time T as given by the Born's rule, consistently with Ref. [8].

Concerning the measurement at time t = 0 of the ND procedure, we have to sum over two possible realizations built with the indices i and l . Thus, renaming the indices, the probability to measure the eigenvalue ai at time t = 0, can be written as

$$\begin{array} { r c l } ; \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cdots \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd. \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd. \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cdk \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cdK \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cdj \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd1 \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd# \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd8 \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd] \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \cd� \\ \$$

In a similar fashion, we obtain the probability P al ( ) = 1 2 P i , j m k , , { PND k ( , j , m i l , , ) + PND k ( , j , m l i , , ) . }

The calculation for the measurement at time t 1 of the ND procedure follows a common reasoning than the one used to get P ai ( ). In this case, indeed, there are two possible realizations over the indices j and m , such that

$$\cdot \quad P ( a _ { j } ) \ & = \ \frac { 1 } { 2 } \sum _ { i, l, m, k } \left \{ P _ { N D } ( k, j, m, i, l ) + P _ { N D } ( k, j, m, l, i ) \right \} \\ \cdot \quad \ & = \ \text{Tr} \left [ \hat { \Pi } _ { j } \hat { U } _ { 1 } \hat { \rho } ^ { 0 } \hat { U } _ { 1 } ^ { \dagger } \right ] = \text{Tr} \left [ \hat { \Pi } _ { j } \hat { \rho } ^ { 0 } ( t _ { 1 } ) \right ]$$

that, as expected, corresponds to the probability of recording the outcome aj after the system has evolved to the state ˆ ρ 0 ( t 1 ). Observing that in general the sum of PND k ( , j , m i l , , ) over three indices returns the Kirkwood-Dirac quasiprobabilities at two times, we can conclude that the quasi-probability distribution P ND ( ∆ ) satisfies the non-invasive measurability condition in terms of the assumption of no-signaling in time.

## Appendix B: Correlators and LG parameter

The two-time quantum correlators C 12 and C 02 , defined in the main text, explicitly read as

$$C _ { 1 2 } \ & = \ \sum _ { j, k } a _ { k } a _ { j } \sum _ { i, l, m } \left [ P _ { N D } ( k, j, m, i, l ) + P _ { N D } ( k, m, j, i, l ) \right ] \quad & \sum _ { i } \\ C _ { 0 2 } \ & = \ \sum _ { i, k } a _ { k } a _ { i } \sum _ { l, j, m } \left [ P _ { N D } ( k, j, m, i, l ) + P _ { N D } ( k, j, m, l, i ) \right ]. \quad & \text{Frc}$$

In the expressions of C 12 and C 02 , we can rearrange the indices to factorize a common PND k ( , j , m i l , , ). For example, we can rewrite C 12 by changing the indices j ↔ m in PND k m ( , , j , i , l ), so that

$$C _ { 1 2 } = \sum _ { i, l, j, m, k } ( a _ { k } a _ { j } + a _ { k } a _ { m } ) P _ { N D } ( k, j, m, i, l ).$$

By repeating this rearrangement of terms also for C 01 and C 02 , we can obtain Eq. (16) in the main text.

Now we are going to show that Kq , 1 = 0. To do this, we have to derive some properties of the ND quasi-probability distribution. Let us consider the case in which m = j and l , i . We want to show that PND k k k i i ( , , , , ¯ ) = -PND k k k i i ( ¯ , ¯ , ¯ , , ¯ ). For a two-level system and binary observables, we have that ˆ Π = k I -ˆ Π ¯ k . Using this relation, by direct calculation, it holds that

$$U _ { 2, k k } U _ { 2, k k } ^ { * } \, & = \, \text{Tr} \left [ \hat { \Pi } _ { k } \hat { U } _ { 2 } \hat { \Pi } _ { k } \hat { U } _ { 2 } ^ { \dagger } \right ] = \text{Tr} \left [ \hat { \Pi } _ { \bar { k } } \hat { U } _ { 2 } \hat { \Pi } _ { \bar { k } } \hat { U } _ { 2 } ^ { \dagger } \right ] \\ & = \, U _ { 2, k k } U _ { 2, \bar { k k } } ^ { * } \,.$$

Analogously, we get that

$$U _ { 1, k i } U _ { 1, \bar { k } } ^ { * } \ = \ \langle \bar { i } | \, \hat { U } _ { 1 } ^ { \dagger } ( \mathbb { I } - \hat { \Pi } _ { \bar { k } } ) \hat { U } _ { 1 } \, | i \rangle = - U _ { 1, \bar { k } } ^ { * } U _ { 1, \bar { k } i }.$$

Hence,

$$P _ { N D } ( k, k, i, \bar { i } ) \ & = \ U _ { 2, k k } U _ { 1, k i } \rho _ { \bar { i i } } ^ { 0 } U _ { 1, \bar { i k } } ^ { * } U _ { 2, k k } ^ { * } \\ & = \ - U _ { 2, \bar { k k } } U _ { 1, \bar { k i } } \rho _ { \bar { i i } } ^ { 0 } U _ { 1, \bar { k } } ^ { * } U _ { 2, \bar { k k } } ^ { * } \\ & = \ - P _ { N D } ( \bar { k }, \bar { k }, \bar { k }, i, \bar { i } ).$$

In a similar way, it can be proven that PND k k k i i ( , ¯ , ¯ , , ¯ ) = -PND k k k i i ( ¯ , , , , ¯ ). At this point, let us take the expression of Kq , 1 [Eq. (18) in the main text]. Separating the contribution for aj = ak and aj , ak and recalling that a 2 k = 1, Kq , 1 can be equivalently written as

$$K _ { q, 1 } = \sum _ { i } \left \{ \sum _ {$$

Using the properties of the ND quasi-probability derived above, we have that both the sums vanish with the result that Kq , 1 = 0.

Let us now derive the simplified expression of Kq , 2 as given by Eq. (21) in the main text. We are going to use the relation P k PND k ( , j , ¯ j , i , i ) = 0, whose proof is in the following 1-line calculation:

$$\sum _ { k } P _ { N D } ( k, j, \bar { j }$$

From this equality, we get PND k ( , j , ¯ j , i , i ) = -PND k ( ¯ , j , ¯ j , i , i ). Moreover, recalling that PND k m ( , , j , i , l ) = ( PND k ( , j , m i l , , )) ∗ , we also have that P j PND k ( , j , ¯ j , i , i ) = 2 Re h PND k ( , j , ¯ j , i , i ) i . As a result, the contribution Kq , 2 in Eq. (19) (with constraints m = ¯ and j l = 1) simplifies to

$$\begin{array} {$$

that corresponds to Eq. (21) in the main text.

## Acknowledgments

The authors acknowledge the fruitful discussion with Andrea Smirne. P.S. acknowledges financial support from INFN. S.G. would like to thank the PRIN project 2022FEXLYB Quantum Reservoir Computing (QuReCo), the PNRR MUR project PE0000023-NQSTI funded by the European UnionNext Generation EU, and the MISTI Global Seed Funds MITFVG Collaboration Grant 'Revealing and exploiting quantumness via quasiprobabilities: from quantum thermodynamics to quantum sensing'.

## Authors' contributions

P.S. designed and performed the research. P.S. and S.G. discussed the results and wrote the paper.

aps.org/doi/10.1103/PhysicsPhysiqueFizika.1.195 .

- [2] A. J. Leggett and Anupam Garg. Quantum mechanics versus macroscopic realism: Is the flux there when nobody looks? Phys. Rev. Lett. , 54:857-860, Mar 1985. doi: 10.1103 PhysRevLett.54.857. /
- [3] Stuart J. Freedman and John F. Clauser. Experimental Test of Local Hidden-Variable Theories. Phys. Rev. Lett. , 28:938-941, Apr 1972. doi:10.1103 PhysRevLett.28.938. / URL https:// link.aps.org/doi/10.1103/PhysRevLett.28.938 .
- [4] Alain Aspect, Philippe Grangier, and G´ erard Roger. Experimental Realization of Einstein-Podolsky-Rosen-Bohm Gedankenexperiment: A New Violation of Bell's Inequalities. Phys. Rev. Lett. , 49:91-94, Jul 1982. doi: 10.1103 PhysRevLett.49.91. / URL https://link.aps.org/ doi/10.1103/PhysRevLett.49.91 .
- [5] Jian-Wei Pan, Dik Bouwmeester, Harald Weinfurter, and Anton Zeilinger. Experimental Entanglement Swapping: Entangling Photons That Never Interacted. Phys. Rev. Lett. , 80:3891-3894, May 1998. doi:10.1103 PhysRevLett.80.3891. / URL https: //link.aps.org/doi/10.1103/PhysRevLett.80.3891 .
- [6] Owen JE Maroney and Christopher G Timpson. Quantum-vs. macro-realism: What does the Leggett-Garg inequality actually test? arXiv preprint arXiv:1412.6139 , 2014.
- [7] David Schmid. A review and reformulation of macroscopic realism: resolving its deficiencies using the framework of generalized probabilistic theories. Quantum , 8:1217, January 2024. ISSN 2521-327X. doi:10.22331 q-2024-01-03-1217. / URL https://doi.org/10.22331/q-2024-01-03-1217 .
- [8] J. J. Halliwell. Leggett-Garg inequalities and no-signaling in time: A quasiprobability approach. Phys. Rev. A , 93:022123, Feb 2016. doi:10.1103 PhysRevA.93.022123. / URL https: //link.aps.org/doi/10.1103/PhysRevA.93.022123 .
- [9] Vladimir B. Braginsky, Yuri I. Vorontsov, and Kip S. Thorne. Quantum Nondemolition Measurements. Science , 209(4456): 547-557, 1980. doi:10.1126 science.209.4456.547. / URL https://www.science.org/doi/abs/10.1126/science. 209.4456.547 .
- [10] V. B. Braginski˘ ı, Farid Ya. Khalili, and Kip S. Thorne. Quantum measurement . Cambridge University Press, Cambridge [England] ;, 1992. ISBN 052141928X.
- [11] Carlton M. Caves. Quantum-Mechanical Radiation-Pressure Fluctuations in an Interferometer. Phys. Rev. Lett. , 45:7579, Jul 1980. doi:10.1103 PhysRevLett.45.75. / URL https: //link.aps.org/doi/10.1103/PhysRevLett.45.75 .
- [12] Carlton M. Caves, Kip S. Thorne, Ronald W. P. Drever, Vernon D. Sandberg, and Mark Zimmermann. On the measurement of a weak classical force coupled to a quantum-mechanical oscillator. I. Issues of principle. Rev. Mod. Phys. , 52:341-392, Apr 1980. doi:10.1103 RevModPhys.52.341. / URL https: //link.aps.org/doi/10.1103/RevModPhys.52.341 .
- [13] Paolo Solinas, Dmitri V. Averin, and Jukka P. Pekola. Work and its fluctuations in a driven quantum system. Phys. Rev. B , 87:060508, Feb 2013. doi: 10.1103 PhysRevB.87.060508. / URL http://link.aps. org/doi/10.1103/PhysRevB.87.060508 .
- [14] P. Solinas and S. Gasparinetti. Full distribution of work done on a quantum system for arbitrary initial states. Phys. Rev. E , 92:042150, Oct 2015. doi: 10.1103 PhysRevE.92.042150. / URL http://link.aps. org/doi/10.1103/PhysRevE.92.042150 .
- [15] P. Solinas and S. Gasparinetti. Probing quantum interference e ff ects in the work distribution. Phys. Rev. A , 94:052103, Nov 2016. doi:10.1103 PhysRevA.94.052103. / URL https: //link.aps.org/doi/10.1103/PhysRevA.94.052103 .
- [16] P. Solinas, M. Amico, and N. Zangh` ı. Measurement of work and heat in the classical and quantum regimes. Phys. Rev. A , 103:L060202, Jun 2021. doi: 10.1103 PhysRevA.103.L060202. URL / https://link.aps. org/doi/10.1103/PhysRevA.103.L060202 .
- [17] P. Solinas, M. Amico, and N. Zangh` ı. Quasiprobabilities of work and heat in an open quantum system. Phys. Rev. A , 105:032606, Mar 2022. doi: 10.1103 PhysRevA.105.032606. / URL https://link.aps. org/doi/10.1103/PhysRevA.105.032606 .
- [18] Stefano Gherardini and Gabriele De Chiara. Quasiprobabilities in quantum thermodynamics and many-body systems: A tutorial. arXiv preprint arXiv:2403.17138 , 2024. doi: 10.48550 arXiv.2403.17138. / URL https://doi.org/10. 48550/arXiv.2403.17138 .
- [19] Matteo Lostaglio, Alessio Belenchia, Amikam Levy, Santiago Hern´ andez-G´mez, o Nicole Fabbri, and Stefano Gherardini. Kirkwood-Dirac quasiprobability approach to the statistics of incompatible observables. Quantum , 7:1128, 2023. doi:10.22331 q-2023-10-09-1128. / URL https:// quantum-journal.org/papers/q-2023-10-09-1128/ .
- [20] N. Yunger Halpern, B. Swingle, and J. Dressel. Quasiprobability behind the out-of-time-ordered correlator. Phys. Rev. A , 97:042105, Apr 2018. doi:10.1103 PhysRevA.97.042105. / URL https://link.aps.org/doi/10.1103/PhysRevA. 97.042105 .
- [21] D. R. M. Arvidsson-Shukur, J. Chevalier Drori, and N. Yunger Halpern. Conditions tighter than noncommutation needed for nonclassicality. J. Phys. A: Math. Theor. , 54:284001, 2021. doi:10.1088 1751-8121 ac0289. / / URL https://iopscience.iop.org/article/10.1088/ 1751-8121/ac0289 .
- [22] S. De Bi` evre. Complete Incompatibility, Support Uncertainty, and Kirkwood-Dirac Nonclassicality. Phys. Rev. Lett. , 127:190404, Nov 2021. doi: 10.1103 PhysRevLett.127.190404. / URL https://link. aps.org/doi/10.1103/PhysRevLett.127.190404 .
- [23] A. Budiyono and H. K. Dipojono. Quantifying quantum coherence via Kirkwood-Dirac quasiprobability. Phys. Rev. A , 107:022408, Feb 2023. doi:10.1103 PhysRevA.107.022408. / URL https://link.aps.org/doi/10.1103/PhysRevA. 107.022408 .
- [24] R. Wagner, Z. Schwartzman-Nowik, I. L. Paiva, A. Te'eni, A. Ruiz-Molero, R. Soares Barbosa, E. Cohen, and E. F. Galv˜ ao. Quantum circuits for measuring weak values, Kirkwood-Dirac quasiprobability distributions, and state spectra. Quantum Sci. Technol. , 9:015030, 2024. doi: 10.1088 2058-9565 ad124c. / / URL https://iopscience. iop.org/article/10.1088/2058-9565/ad124c .
- [25] David R. M. Arvidsson-Shukur, William F. Braasch Jr., Stephan De Bievre, Justin Dressel, Andrew N. Jordan, Christopher Langrenez, Matteo Lostaglio, Je ff S. Lundeen, and Nicole Yunger Halpern. Properties and Applications of the KirkwoodDirac Distribution. arXiv preprint arXiv:2403.18899 , 2024. doi:10.48550 arXiv.2403.18899. / URL https://doi.org/ 10.48550/arXiv.2403.18899 .
- [26] Santiago Hern´ andez-G´mez, o Takuya Isogawa, Alessio Belenchia, Amikam Levy, Nicole Fabbri, Stefano Gherardini, and Paola Cappellaro. Interferometry of quantum correlation functions to access quasiprobability distribution of work. arXiv preprint arXiv:2405.21041 , 2024. URL https://arxiv. org/abs/2405.21041 .
- [27] Leonid S. Levitov, Hyunwoo Lee, and Gordey B. Lesovik. Electron counting statistics and coherent states of

- electric current. J. Math. Phys. , 37(10):4845-4866, 1996. doi:10.1063 1.531672. / URL https://pubs.aip.org/ aip/jmp/article-abstract/37/10/4845/454669/ Electron-counting-statistics-and-coherent-states? redirectedFrom=fulltext .
- [28] Yu V. Nazarov and M. Kindermann. Full counting statistics of a general quantum mechanical variable. Eur. Phys. J. B , 35:413-420, 2003. doi:10.1140 epjb e2003-00293-/ / 1. URL https://link.springer.com/article/10.1140/ epjb/e2003-00293-1 .
- [29] A. A. Clerk. Full counting statistics of energy fluctuations in a driven quantum resonator. Phys. Rev. A , 84:043824, Oct 2011. doi:10.1103 PhysRevA.84.043824. / URL https: //link.aps.org/doi/10.1103/PhysRevA.84.043824 .
- [30] Patrick P. Hofer and A. A. Clerk. Negative Full Counting Statistics Arise from Interference E ff ects. Phys. Rev. Lett. , 116: 013603, Jan 2016. doi:10.1103 PhysRevLett.116.013603. URL / https://link.aps.org/doi/10.1103/PhysRevLett. 116.013603 .
- [31] Patrick P. Hofer. Quasi-probability distributions for observables in dynamic systems. Quantum , 1:32, October 2017. ISSN 2521-327X. doi:10.22331 q-2017-10-12-32. / URL https: //doi.org/10.22331/q-2017-10-12-32 .
- [32] Patrick P. Potts. Certifying Nonclassical Behavior for Negative Keldysh Quasiprobabilities. Phys. Rev. Lett. , 122:110401, Mar 2019. doi:10.1103 PhysRevLett.122.110401. / URL https:// link.aps.org/doi/10.1103/PhysRevLett.122.110401 .
- [33] Rusko Ruskov, Alexander N. Korotkov, and Ari Mizel. Signatures of quantum behavior in single-qubit weak measurements. Phys. Rev. Lett. , 96:200404, May 2006. doi: 10.1103 PhysRevLett.96.200404. / URL https://link.aps. org/doi/10.1103/PhysRevLett.96.200404 .
- [34] Andrew N. Jordan, Alexander N. Korotkov, and Markus B¨ uttiker. Leggett-garg inequality with a kicked quantum pump. Phys. Rev. Lett. , 97:026805, Jul 2006. doi: 10.1103 PhysRevLett.97.026805. / URL https://link.aps. org/doi/10.1103/PhysRevLett.97.026805 .
- [35] M. E. Goggin, M. P. Almeida, M. Barbieri, B. P. Lanyon, J. L. O'Brien, A. G. White, and G. J. Pryde. Violation of the leggettgarg inequality with weak measurements of photons. Proceedings of the National Academy of Sciences , 108(4):12561261, 2011. doi:10.1073 pnas.1005774108. / URL https: //www.pnas.org/doi/abs/10.1073/pnas.1005774108 .
- [36] Markus Arndt, Olaf Nairz, Julian Vos-Andreae, Claudia Keller, Gerbrand van der Zouw, and Anton Zeilinger. Waveparticle duality of C60 molecules. Nature , 401(6754):680682, 1999. doi:10.1038 44348. URL / https://doi.org/10. 1038/44348 .
- [37] Stefan Gerlich, Sandra Eibenberger, Mathias Tomandl, Stefan Nimmrichter, Klaus Hornberger, Paul J. Fagan, Jens T¨ uxen, Marcel Mayor, and Markus Arndt. Quantum interference of large organic molecules. Nature Communications , 2(1):263, 2011. doi:10.1038 ncomms1263. / URL https://doi.org/ 10.1038/ncomms1263 .
- [38] Markus Aspelmeyer. When Zeh Meets Feynman: How to Avoid the Appearance of a Classical World in Gravity Experiments , pages 85-95. Springer International Publishing, Cham, 2022. ISBN 978-3-030-88781-0. doi: 10.1007 978-3-030-88781-0˙5. / URL https://doi.org/10. 1007/978-3-030-88781-0\_5 .
- [39] Tim M. Fuchs, Dennis G. Uitenbroek, Jaimy Plugge, Noud van Halteren, Jean-Paul van Soest, Andrea Vinante, Hendrik Ul-
- bricht, and Tjerk H. Oosterkamp. Measuring gravity with milligram levitated masses. Science Advances , 10:1-6, 2024. doi: 10.1126 sciadv.adk2949. /
- [40] Youngseok Kim, Andrew Eddins, Sajant Anand, Ken Xuan Wei, Ewout van den Berg, Sami Rosenblatt, Hasan Nayfeh, Yantao Wu, Michael Zaletel, Kristan Temme, and Abhinav Kandala. Evidence for the utility of quantum computing before fault tolerance. Nature , 618(7965):500-505, 2023. doi: 10.1038 s41586-023-06096-3. / URL https://doi.org/10. 1038/s41586-023-06096-3 .
- [41] Sergey Bravyi, Andrew W. Cross, Jay M. Gambetta, Dmitri Maslov, Patrick Rall, and Theodore J. Yoder. High-threshold and low-overhead fault-tolerant quantum memory. Nature , 627(8005):778-782, 2024. doi:10.1038 s41586-024-07107-7. / URL https://doi.org/10.1038/s41586-024-07107-7 .
- [42] Pingal Pratyush Nath, Debashis Saha, Dipankar Home, and Urbasi Sinha. Single-System-Based Generation of Certified Randomness Using Leggett-Garg Inequality. Phys. Rev. Lett. , 133: 020802, Jul 2024. doi:10.1103 PhysRevLett.133.020802. URL / https://link.aps.org/doi/10.1103/PhysRevLett. 133.020802 .
- [43] George C. Knee, Kosuke Kakuyanagi, Mao-Chuang Yeh, Yuichiro Matsuzaki, Hiraku Toida, Hiroshi Yamaguchi, Shiro Saito, Anthony J. Leggett, and William J. Munro. A strict experimental test of macroscopic realism in a superconducting flux qubit. Nature Communications , 7(1):13253, 2016. doi: 10.1038 ncomms13253. / URL https://doi.org/10.1038/ ncomms13253 .
- [44] Johannes Kofler and ˇ aslav C Brukner. Conditions for Quantum Violation of Macroscopic Realism. Phys. Rev. Lett. , 101:090403, Aug 2008. doi: 10.1103 PhysRevLett.101.090403. / URL https://link. aps.org/doi/10.1103/PhysRevLett.101.090403 .
- [45] P. Solinas, H. J. D. Miller, and J. Anders. Measurementdependent corrections to work distributions arising from quantum coherences. Phys. Rev. A , 96:052115, Nov 2017. doi: 10.1103 PhysRevA.96.052115. / URL https://link.aps. org/doi/10.1103/PhysRevA.96.052115 .
- [46] Gabriele De Chiara, Paolo Solinas, Federico Cerisola, and Augusto J. Roncaglia. Ancilla-Assisted Measurement of Quantum Work , pages 337-362. Springer International Publishing, Cham, 2018. ISBN 978-3-319-99046-0. doi: 10.1007 978-3-319-99046-0˙14. / URL https://doi.org/ 10.1007/978-3-319-99046-0\_14 .
- [47] Tobias Fritz. Quantum correlations in the temporal ClauserHorne-Shimony-Holt (CHSH) scenario. New Journal of Physics , 12(8):083055, aug 2010. doi:10.1088 1367-/ 2630 12 8 083055. / / / URL https://dx.doi.org/10.1088/ 1367-2630/12/8/083055 .
- [48] Johannes Kofler and ˇ aslav C Brukner. Condition for macroscopic realism beyond the Leggett-Garg inequalities. Phys. Rev. A , 87:052115, May 2013. doi: 10.1103 PhysRevA.87.052115. / URL https://link.aps. org/doi/10.1103/PhysRevA.87.052115 .
- [49] Clive Emary, Neill Lambert, and Franco Nori. Leggett-Garg inequalities. Reports on Progress in Physics , 77(1):016001, dec 2013. doi:10.1088 0034-4885 77 1 016001. / / / / URL https:// dx.doi.org/10.1088/0034-4885/77/1/016001 .
Paolo Solinas Stefano Gherardini
Quantum Non-Demolition Measurements and Leggett-Garg inequality
Quantum non-demolition measurements define a non-invasive protocol to extract
information from a quantum system that we aim to monitor. They exploit an
additional quantum system that is sequentially coupled to the system.
Eventually, by measuring the additional system, we can extract information
about temporal correlations developed by the quantum system dynamics with
respect to a given observable. This protocol leads to a quasi-probability
distribution for the measured observable outcomes, which can be negative. We
prove that the presence of these negative regions is a necessary and sufficient
condition for the violation of macrorealism. This is a much stronger condition
than the violation of the Leggett-Garg inequalities commonly used for the same
task. Indeed, we show that there are situations in which Leggett-Garg
inequalities are satisfied even if the macrorealism condition is violated. As a
consequence, the quantum non-demolition protocol is a privileged tool to
identify with certainty the quantum behavior of a system. As such, it has a
vast number of applications to different fields from the certification of
quantumness to the study of the quantum-to-classical transition.