A Wigner Quasiprobability Distribution of Work

In this article, we introduce a quasiprobability distribution of work that is based on the Wigner function. This proposal rests on the idea that the work conducted on an isolated system can be coherently measured by coupling the system to a quantum measurement apparatus. In this way, a quasiprobability distribution of work can be defined in terms of the Wigner function of the apparatus. This quasidistribution contains the information of the work statistics and also holds a clear operational definition that can be directly measured in a real experiment. Moreover, it is shown that the presence of quantum coherence in the energy eigenbasis is related with the appearance of features related to non-classicality in the Wigner function such as negativity and interference fringes. On the other hand, from this quasiprobability distribution, it is straightforward to obtain the standard two-point measurement probability distribution of work and also the difference in average energy for initial states with coherences.


I. INTRODUCTION
The notion of work is one of the most basic and fundamental concepts in physics, particularly in thermodynamics.During the last decades, several attempts have been made to obtain the work statistics for nonequilibrium thermodynamic transformations in the quantum regime.These definitions where motivated by the idea of extending classical fluctuation theorems [1][2][3][4][5] to the quantum regime.In order to describe the thermodynamics of general non-equilibrium quantum processes, it is necessary to provide a general definition of work valid for any quantum system and process.However, this task presents serious difficulties.This is due, in the first place, to the fact that many concepts belonging to the classical definition of work cannot be translated to quantum mechanics.For example, the basic definition of the work that a force performs on a particle along a trajectory cannot be used in quantum mechanics because of the lack of meaning of trajectories in the theory, although recently a definition of quantum work was made by considering bohmian trajectories [6].A great advance came in the area with the definition of the two point measurement protocol (TPM) to define work in driven isolated quantum systems [3,4,7,8].This definition is based on the simple observation that, for an isolated system, work is a random variable associated to the difference in energy along the process.Thus, in order to determine this random value one should make an energy measurement at the beginning and another at the end of the process.This definition is not only straightforward in an operational sense, but it also recovers the results of the fluctuation theorems for quantum systems [3][4][5][7][8][9] and was verified experimentally in different platforms [10][11][12][13][14][15].
However, there is a caveat with the TPM when one considers initial states that have coherences in the en- * federico.cerisola@eng.ox.ac.uk ergy basis.This is because the first energy measurement destroys these coherences, and therefore the TPM scheme is insensitive to quantum coherence between different energy subspaces.This leads to undesirable consequences, for instance, related with the fact that the average work done in the process is different from the change in the average energy of the system.Moreover, it has been shown that it is impossible to define a probability distribution of work that satisfies at the same time, the fluctuation theorems and whose mean value of work equals the average energy change for states with coherences [16].This has led the community to consider different approaches to generalize the work distribution [6,[17][18][19] including some proposals for quasiprobability distributions of work [20][21][22][23][24][25][26][27][28].
In this article, we propose a distribution based on the Wigner function.This construction relays on the fact that the work probability distribution can also be coherently measured by coupling the system to a quantum apparatus and making a single measurement over the apparatus, i.e. a single-measurement protocol (SM) [12,29,30].In this way, the final state of the apparatus contains the information about the work distribution and one can define a quasiprobability distribution [31].This approach provides a clear operational definition with an immediate experimental implementation.In addition, the Wigner function is represented using coordinates that have an intuitive interpretation in terms of time and energy associated with the work.Moreover, it can be shown that the presence of quantum coherence is related with the appearance of characteristics related to non-classicality in the Wigner function, such as negativity and interference fringes.On the other hand, for coherence free states this definition agrees with the standard two point measurement probability distribution of work.
The paper is organized as follows.In Sec.II we briefly discuss the two point measurement scheme, and the single measurement protocol.In Sec.III we introduce the quasiprobability distribution of work based on the Wigner function, showing how it works for initial states of the system with and without coherence.In Sec.IV we discuss experimental implementations, and we end with discussions and conclusion in Sec.V.

II. WORK STATISTICS
As mentioned, we are interested in the work distribution for isolated quantum systems that are subjected to an external driving.In this way, the external work can be associated to the energy change of the system.The typical scenario consists of a system S that starts in a given initial state, ρ S , and after that, a driving represented by a unitary evolution U is applied.The driving is such that it changes the Hamiltonian from an initial H to a final one H where Π n are the projectors on each energy subspace.In this case, what we know is that the average change of energy in the system is where U ρ S U † is the final state of the system.Clearly, it would be desirable that the average work obtained from the corresponding probability distribution equals this average energy change.This requisite is equivalent to asking that the first law of thermodynamics for mean values is satisfied for an isolated system.However, it can be shown that if one imposes that the statistics of work is consistent with the standard fluctuation theorems, the distribution of work should be defined by the two-point measurement protocol [3,4,7,8].In this case, although the resulting work average coincides with the mean energy difference for initial equilibrium states (diagonal in the initial energy eigenbasis), it is different for initial states with coherences.

A. The two-point measurement protocol
The two-point measurement protocol allows us to define a probability distribution of work consistent with fluctuation theorems.In order to do that, one should define a work value for each realization.This is done in terms of the difference of two energy values that are obtained by making two energy measurements, one at the beginning and the other one at the end of the driving.In this way, the corresponding probability distribution can be written as where p n is the probability of obtain E n in the first energy measurement, and p m|n is the conditional probability of obtaining Ẽm at the end given that E n was obtained at the beginning.Therefore, if the initial state is already diagonal in the energy eigenbasis, the first measurement does not modify the state, and it can be easily verified that the mean value of work equals the average energy difference.Indeed, from this probability distribution one can evaluate the mean value of work as where ρS = n Π n ρ S Π n is the dephased initial state.This state is obtained by removing all the coherences between different energy subspaces of the initial Hamiltonian, and it is equivalent to the state resulting the following asymptotic temporal average ρS = lim Therefore, unless the initial state is diagonal in the basis of the initial hamiltonian, the average given by the TPM is different from the average energy difference of the system.In fact, if the initial state is diagonal in this basis, then ρ S can be interpreted as a 'classical' probability distribution over the different energies.In that case, the measurement is not invasive, in the sense that it only reveals the value of the energy in each realization of the experiment.On the other hand, for an initial state with coherences, the initial energy is not well defined and this interpretation is not straightforward.

B. The single-measurement protocol
Another method for assessing to the work probability distribution was introduced in [29].The method is based on the idea that the work measurement can be described in terms of a generalized measurement (POVM).That is, by coupling the system to an ancilla that is finally subjected to a 'single measurement' (SM).In this way, it can be shown that one can obtain the same probability distribution provided the ancilla is properly initialized.
Let us now describe briefly the general method that is summarised in the circuit of Fig. 1.Initially, the system is in the state ρ S and there is an auxiliary system (ancilla) A whose state is described terms of a continuous degree of freedom.In the ancilla's space one can consider two canonically conjugated operators, W A and T A , such that [W A , T A ] = i .Thus, the evolution contains two coherent interactions between S and A: one at the beginning e iH⊗T A / and another e −i H⊗T A / at the end of the driving.Thus, each interaction can be viewed either as a coherent translation in the variable w of the ancilla Circuit that describes the single-measurement protocol from where the work probability distribution can be obtained.
in an amount that depends on the energy of the system or, conversely, as a coherent time-translation (free evolution) of the system whose time interval is proportional to the variable τ of the ancilla.Therefore, one can immediately associate the variables w and τ to energy (work) and time, respectively.Following the protocol of the circuit, if the initial state of A is ρ A , after the interaction with the system its reduced state is where w nm = Ẽm −E n are the different work values.The SM protocol finishes by performing a projective measurement of the observable W A .In this case, for highly localized initial pure states of A, the resulting probability distribution is equivalent to the work distribution of the TPM protocol.Notably, within this formulation one can associate work to an observable acting over the ancillary system.Of course, work is not an observable acting on the system's space [4].
It is important to stress at this point that the entangling interaction between system and apparatus establishes a coherent record of the different values of work.Therefore, the reduced state of the ancilla contains information not only about the probability distribution given by the TPM, but also about the initial state of the system.At the end, the type of measurement that is done over the ancilla (or the battery), determines which information is extracted from the protocol.It is also interesting to notice that this type of interaction appears in a very related task: the work extraction from a quantum system.This can be modeled by adding an interaction between the system and an auxiliary system that acts as a battery in which work is stored [32][33][34].In general, the battery can be thought of as a continuous variable system, an ideal weight, with a hamiltonian like the operator W A .The work extraction process consists on some unitary evolution on the joint system (where the driving on the system is included) that can change the system hamiltonian from H to H.The extracted work, in this way, is stored in the battery.There are a few conditions that should be imposed in this framework in order to ensure that the weight does not provide any thermodynamical resource to the work extraction process [33], one of them is of course energy conservation.It has been shown in [33] that the unitary operations that satisfy these conditions are of the form e iH⊗T A / (U ⊗ I A ) e −i H⊗T A / where U is the driving of the system.Therefore, it is straightforward to see that these are the same operations (up to a sign) used in the SM protocol for measuring work.Thus, there is also a clear operational interpretation of the state of the ancilla as the state of a battery where work is stored.

III. THE WIGNER DISTRIBUTION OF WORK
In the following, we will define a generalized work distribution.The general idea is inspired by the SM protocol.As we just mentioned, the state of the ancilla after the interaction not only holds information about work, but also about the coherences present in the initial state.In order to extract such information, we will evaluate their Wigner function [35,36], P W .The Wigner function is a quasi-probability distribution that is used to represent quantum states in phase space.This is a real-valued function that, unlike their classical counterparts, can be negative for generic quantum states.This property has been widely used as an indicator of quantumness in different contexts, for instance in the study of the quantum-classical transition.
In our case, we will define it for the final state of the ancilla and in terms of the conjugate variables w and τ as This expression is valid for a generic initial state of the ancilla.In order to evaluate it, we will assume that the initial state of the ancilla is a coherent Gaussian state.This assumption not only will allow us to easily perform analytical calculations, but is also and appropriate choice for the description of typical experimental situations.Moreover, Gaussian states are classical, in the sense that they have a positive Wigner function.This guarantees that any negativity appearing in the Wigner function of the ancilla comes exclusively from their interaction with the system.Thus, we consider ρ A = |0, σ 0, σ| a coherent Gaussian state with zero mean and variance σ 2 in W A (and hence zero mean and variance 2 2σ 2 in T A ).After replacing this in Eq. ( 7) (see Appendix A) and using that Π n ρ S e iτ En/ = Π n e iτ H/ ρ S we obtain an expression for the quasiditribution of work for a generic process where N (w | µ, σ) is a normal probability density in w with mean value µ and variance σ 2 (analogously for τ ).
From this expression we can appreciate again the operational interpretation of the variables w and τ that characterize the state of the ancilla.
In the following, we will introduce some notation that will be useful to simplify forthcoming expressions.First, let us recall that the distribution P TPM (w) does not take into account any coherence between the different energy subspaces of H in the initial state ρ S .Therefore, we can associate this probability distribution to the dephased state ρS .It would then be convenient to define the probability distribution P N (w|σ) that is the convolution of P TPM (w) with a normal distribution with zero mean and variance σ 2 Notice that P N (w|σ) is simply the TPM distribution, Eq. ( 3), with the Dirac delta replaced by a normal distribution with the corresponding mean values of work and variance σ 2 .Thus, for a highly localized normal distribution, it satisfies P N (w|σ) − −− → σ→0 P TPM (w).
In order to illustrate the effect of initial coherences let us consider Eq. ( 8), and split it in diagonal (n = n ) and non-diagonal (n = n ) contributions The non-diagonal one corresponds to the contribution of the so-called initial coherences and it is easy to see that × N w w nm + w n m 2 , σ N τ 0, √ 2σ .
A. Quasidistribution for initial dephased states When the initial state of the system is diagonal in the energy basis (Π n ρ S Π n = 0 for n = n ), P (c) W (w, τ ) = 0 and the Wigner function is just that is, it's proportional to the convoluted TPM distribution for every value of τ .Moreover, if we calculate the marginal P W (w), we recover the probability distribution that would be obtained if one measures the observable W A .This expression just reflects a characteristic property of the Wigner function: the partial integration provides the probability distribution corresponding to the other variable.Therefore, for initial states without coherences in the initial energy eigenbasis, P W (w) is exactly P N (w|σ).If in addition σ (w nm − w n m ), ∀n, n , m, m then we recover the probability distribution of work given by the TPM protocol.
In Fig. 2(b), upper panel, we show the distribution P W (w, τ ) for a two-level system S without initial coherences.In the lower panel, we show the marginal of the distribution in w, P W (w), and compare it with the discrete probabilities associated to the TPM.Notice that the area under each Gaussian in the marginal is equal to the corresponding probability in the TPM protocol, and we can also see that it effectively reproduces the ideal TPM distribution.On the other hand, in Fig. 2(a) we show the distribution of work obtained for the same system but using an ancilla that has an initial state with a standard deviation five times smaller.One can easily note that this case is much closer to the ideal projective measurement regime.In this case, the Wigner function is invariant under translations in τ , as expected, since the initial state of the system commutes with the initial Hamiltonian.In Fig. 2(c) we show the distribution for a standard deviation even bigger than the one in Fig. 2(b).As we can see, while the position of the peaks matches the correct work values, there is a significant overlap between the different Gaussians.

B. Effects of quantum coherences
Let us now consider a system with initial coherences.From Eq. ( 10) we can notice that in this case the Wigner function also has Gaussian peaks on each work value w nm , just as it happens for the dephased state.However, there are some additional terms centered around the average of two work values with different initial energy, (w nm + w n m )/2.These terms are the ones that hold the non-trivial dependence on the variable τ and,  10) based on the Wigner function.In the lower panel we show the marginal of w, given by Eq. ( 13), along with the discrete probabilities p(w) corresponding to the usual TPM distribution.
as we will see , they can be negative.This can be easily seen from the following argument.If we look at Eq. ( 10), we have that ∞ −∞ dτ dw P W (w, τ ) = 1, and, in addition, also the integral over the phase space of the first term is equal to one, as it is the Wigner function of the initial dephased state.Therefore, the integral of the second term must be zero.In order to do so, the sum should attain some negative values since they are real.In these terms, except for the Gaussian modulation, the variable τ only appears as a time evolution of the state.
The fact that time appears explicitly only for initial states with coherences has a clear interpretation.If the initial state ρ S is diagonal, then it is a steady state of the initial Hamiltonian, and the state is the same for every instant in time before the driving.On the other hand, if ρ S has coherences, the state evolves due to the free evolution with the initial hamiltonian.This time, of course, is irrelevant at the moment of performing the first projecting energy measurement for the TPM distribution.However, it appears in our approach due to the fact we are performing coherent operations between system and ancilla.Notably, one can also observe that the mean energy difference Eq. ( 2) is not invariant under initial time translations for states with coherences.Thus, given a reference state ρ S , the calculated distribution contains, in principle, information about every initial state that is unitarily connected with ρ S by the initial hamiltonian.Nevertheless, the amplitude of the Wigner function decays exponentially to zero when τ → ±∞ due to the Gaussian modulation.
At the same time, given the complementary nature of the variables work w and time τ , when localizing the Gaussian in w we are delocalizing it in τ .We will come back to this issue when we consider the marginals of the distribution.
We have already shown that if the initial state does not have coherence in the energy basis, the resulting quasidistribution of work function is positive, because the diagonal terms in Eq. ( 10) are all positive.Therefore, if the distribution P W (w, τ ) has some negativities, it is a signature of the presence of coherences in the initial state.This can be clearly seen in the upper panel of Fig. 3, where the quasiprobability distribution of work P W (w, τ ) of a two-level system is plotted.The hamiltonians, drivings and ancilla parameters are identical to those of Fig. 2.Moreover, in both cases the initial state of the system has the same probability distribution in the energy basis.The only difference between Figs. 2 and 3 is that in the second one the initial state has coherence between the two energy levels.Comparing both figures, we can notice that effectively we have the same Gaussian distributions over the same work values.The key difference lies in the fact that, for the initial coherent state, the quasidistribution displays additional oscillations that can become negative.This interference fringes indicate the presence of non-classicality in the Wigner function and in the initial state of the system.
In the lower panel of Fig. 3 we show the marginal probability distribution for w.Comparing them with Fig. 2 we can notice that the marginal distributions are equivalent.In Fig. 3 we also show the work distribution for the same system but using an initial state of the ancilla A with different standard deviations.In the case of the smaller standard deviation (corresponding to an ideal Figure 3. Wigner function for work of a two-level system using a Gaussian ancilla.The parameters used are the same as in Fig. 2, except that now the initial density matrix has non-diagonal elements, ρS = (I + σx/2 + σy/2 + σz/4)/2.The upper panel shows the distribution PW (w, τ ) Eq. ( 10) based on the Wigner function.We notice now, because of of the initial coherences, the appearence of negative values in the distribution.The lower panel shows the marginal of w, given by Eq. ( 13), and it is compared to the work values and respective discrete probabilities p(w) that appear in the usual TPM distribution.
projective measurement) we can see that the marginal probability recovers that of the TPM distribution.Notably, in the case of a bigger standard deviation the interference between different Gaussian picks modifies the distribution of the values of w, and there are corrections due to coherences as it is shown in Eq. (A2).In this limit, the marginal distribution may not even coincide with that of the corresponding dephased state.This behaviour is similar to what happen when one makes a weak measurement [37].To understand when it is possible to observe these differences, lets note that when the marginal for w is calculated from Eq. ( 10), the diagonal terms give exactly the convoluted distribution P N (w|σ).For the non-diagonal contributions, we have time-averages of the form This operation is similar to a dephasing map in the energy basis, but there is a significant difference since this average is weighted with a normal distribution with variance 2 /(2σ 2 ) centered in the origin.The bigger the variance of the Gaussian (and therefore the smaller σ), more values of τ enter in the time-average.Therefore, in the limit of small σ we expect the non-diagonal terms to average to zero.Hence, one can show that if σ (E n − E n )/2, ∀n, n , independent of the initial state, Thus, the marginal of the quasiprobability distribution reproduces the TPM distribution.

C. Calculation of mean values
Given the formalism associated to the Wigner function [36], one can easily obtain average values from this quasidistribution.In fact, using the Wigner-Weyl representation [36] of an operator A acting on the ancillary space, their mean value is just (17) For instance, the mean value of work is just the mean value of the operator W A , and it is obtained by integrating the function w over the phase space The other typical average that is calculated in the context of fluctuation theorems, where the system is initially in thermal equilibrium at inverse temperature β, is e −βw .This is easily done by integration of the function e −βw .In all cases, the calculated mean values depend on the initial state of the ancilla.As it can be easily proven, for any observable of the type f (W A ), in the limit of σ → 0 their averages converge to the values associated with the TPM distribution.

D. Energy difference in the presence of coherences
Finally, we will show another interesting property that is contained in this quasidistribution.We have seen that, unless the initial state of the system is diagonal in the energy eigenbasis, the difference in mean energy and the mean value of work (Eqs.( 2) and ( 4)) do not coincide.Thus, the TPM distribution does not provide any information about the initial coherences.Notably, as we will show, this information is also contained in this quasiprobability distribution.
In order to do so, let us consider the average in phase space of the function g τ0 (w, τ ) = w δ(τ − τ 0 ) (see Appendix A).Using the Wigner-Weyl transform [36], it corresponds to the expectation value of the operator Ĝp0 = (W A |τ 0 τ 0 | + |τ 0 τ 0 | W A ) /2 measured over the ancilla.It can be easily shown that this average, which is equivalent to the integral of the function w weighed by the Wigner function along an horizontal line at where ∆E τ0 = tr HUρ S (−τ 0 ) U † − tr[Hρ S (−τ 0 )] is the mean energy difference for a situation where the driving U is turned on at time −τ 0 , and the proportionality constant is just equal to the Gaussian modulation at τ 0 , N τ 0 0, √ 2σ (see Appendix A).Therefore, when τ 0 = 0, this is just proportional to the 'initial' energy difference ∆E in Eq. ( 2).As we have shown, from this quasiprobability distribution we can calculate, not only the energy difference corresponding to the actual initial state, but also for the set of states ρ S (τ ), τ ∈ R.This set can be viewed as different 'initial times' at which the driving is turned on starting from a reference state ρ S at time zero.This is so because this set of initial states is connected with ρ S by a free hamiltonian evolution.
Interestingly, for a Gaussian initial state of the ancilla, one obtains the correct value ∆E independent of their initial variance σ.However, since there is a Gaussian modulation centered around τ 0 = 0 (the proportionality constant), the error in its determination increases as one localizes the initial state of the ancilla in the variable w.But, if one reduces the value of σ, the estimation of P TPM (w) gets worse.Therefore, one can also appreciate in this case the complementary nature of the variables w and τ .

IV. POSSIBLE EXPERIMENTAL IMPLEMENTATIONS
The measurement of the quasi-probability distribution that we propose requires two fundamental ingredients: (i) coherent control of two degrees of freedom of system and ancilla in order to implement the interactions of the SM protocol; (ii) being able to measure the Wigner function of the ancilla.In particular, implementing the SM requires the ability of performing translations of the ancilla conditioned on the energy degree of freedom on which the work is performed.There is a great variety of systems where this sort of interactions can be implemented, and an experimental realization of the SM protocol has been realised using cold atoms [12].However, it is not clear how one can implement the measurement of the Wigner function in such platform.Nevertheless, there are systems where both requirements are in principle satisfied and in what follows we will briefly describe two of them.
The first example is given by superconducting qubits coupled to a cavity, i.e. circuit quantum electrodynamics.In these systems, one has a charge qubit formed by a superconducting island coupled to a Josephson junction, and the two states correspond to the presence or absence of excess Cooper pairs in the island [38].Furthermore, the qubit circuit can be coupled to a wave-guide that acts as a microwave cavity where coherent states or states with a well defined number of photons can be stored [39].For instance, in Ref. [40] they generate coherent displacements of the state of the cavity depending on the state of the qubit.This interaction is exactly what is needed for implementing the protocol where the qubit acts as the system and the cavity as the ancilla.On the other hand, in a different coupling regime between qubit and cavity, this same scheme has been used to measure the Wigner function of the state of the field in the cavity [41].
The second possible platform are trapped ions.In this case, ions are trapped in an electric potential such that the motional degrees of freedom of the ion are subjected to an effective harmonic oscillator potential [42].At the same time, using the interactions between the electronic degree of freedom and the position of the ion, it is possible to generate coherent, squeezed and Fock states of the oscillator [42].In particular, in different experiments [43,44] it has been shown that one can apply forces on the ion depending on its electronic state, and in this way displacements in phase space depending on the qubit state can be coherently implemented.Again, this is the interaction needed to perform the protocol.At last, the Wigner function of the motional degree of freedom of trapped ions has been successfully measured [45].

V. CONCLUSIONS
In this work we introduced a generalization of the probability distribution of work based on the Wigner function.The starting point is the SM protocol proposed in [29], where an ancilla is coupled to the system whose work one wants to measure in order to keep a coherent record of all possible work values.Following this idea, we define the Wigner function of the final state of the ancilla.This quasi-probability distribution contains all the information regarding both work and coherence in the initial state of the system.In fact, initial quan-tum coherence in the system results in negativities in the quasi-probability distribution of work, a clear signature of non classicality.In this case, we can also recover the mean value of energy, which is different from the average work for states with coherences.Moreover, we show that from this quasi-probability distribution one can easily recover the standard TPM distribution simply by integrating over the time variable.In addition, we show that given that the average work and other quantities of interest can be obtained as the mean value of an operator acting on the ancillary space, it is easy to calculate mean values using the formalism of the Wigner function.The quasi-probability distribution here defined has certain similarities with the one proposed in [21,22,46].The way in which the distribution is defined there, is also inspired on the SM scheme [47], and requires the preparation of a coherent superposition of the ancilla between two momentum eigenstates, |p + |−p , together with the implementation of an interaction analogous to that of the SM.At the end of the protocol, the relative phase between these states is measured and a quasiprobability distribution that contains information about work and coherence is obtained [47].In contrast, our proposal has a clear operational interpretation and direct experimental application, as it is simply the Wigner function of the final state of the measurement apparatus.Moreover, our protocol not only contains all the information of Ref. [21,22,46], but also for coherent initial states it has additional information on the dependence of the time variable, τ .From a practical point of view, our protocol does not need ideal (non-physical) states and it is easy to adapt to any initial state of the ancilla.Here, we have just developed the case of Gaussian states given that they are easy to treat analytically and are typically appropriate to model experimental conditions.However, this whole analysis can be repeated for any initial state.We hope that this approach to the work distribution can shed some light to elucidate the effects of quantum coherences in thermodynamic transformations.

Figure 2 .
Figure 2. Wigner function of work for a two-level system using a Gaussian ancilla.The initial hamiltonian is H = Eσ+σ−, with σ± the Pauli creation and annihilation operators.The unitary driving is given by U = ( √ 2I + iσx + iσz)/2 and the final hamiltonian is H = 2Eσ+σ−.The initial state of the system is ρS = (I + σz/4)/2 and the variance of the initial Gaussian packets of the ancilla are (a) σ = 0.02E, (b) σ = 0.1E and (c) σ = 0.35E.The upper panel shows the distribution PW (w, τ ) of Eq. (10) based on the Wigner function.In the lower panel we show the marginal of w, given by Eq. (13), along with the discrete probabilities p(w) corresponding to the usual TPM distribution.