Abstract
We report on the experimental quantification of the contribution to non-equilibrium entropy production stemming from the quantum coherence content in the initial state of a qubit exposed to both coherent driving and dissipation. Our experimental demonstration builds on the exquisite experimental control of the spin state of a nitrogen-vacancy defect in diamond and is underpinned, theoretically, by the formulation of a generalized fluctuation theorem designed to track the effects of quantum coherence. Our results provide significant evidence of the possibility to pinpoint the genuinely quantum mechanical contributions to the thermodynamics of non-equilibrium quantum processes in an open quantum systems scenario.
Similar content being viewed by others
Introduction
The irreversible character of most physical processes is, apparently, at odds with the inherent reversibility of the fundamental laws of physics. The way time-reversible quantum laws, which governs the interactions of microscopic systems, gives rise to the irreversible nature of macroscopic phenomena is a very open field of investigation1. In this regard, a breakthrough has been provided by the extension of the second law of thermodynamics into the quantum realm through the so-called fluctuation theorems2,3,4,5. A celebrated instance of this is the integral fluctuation theorem, which stems from Jarzynski’s identity2 and Crooks’ relation3, and connects the non-equilibrium energy fluctuation statistics of unital processes with the corresponding free-energy changes. Operationally, the standard approach to the quantification of energy and entropy fluctuations in non-equilibrium contexts is by using the celebrated two-point measurement (TPM) scheme6, which requires two projective measurements, at the beginning and at the end of the dynamical process under scrutiny7,8.
Despite the clear success of the TPM scheme, evidenced by successful experimental verification in nuclear magnetic resonance9,10, trapped-ion11,12,13, superconducting-qubit14, nitrogen-vacancy (NV) centres15,16 and linear optics settings17,18,19, the scheme has significant limitations when considering the role played by quantum features in the statistics of energy fluctuations. In fact, any quantum coherence in the initial state of the quantum system, and expressed in the measurement basis, is washed away. As a result, the system dynamics following the first measurement is strongly affected. This has motivated recent efforts aimed at modifying the TPM scheme to take into account the presence of quantum features, particularly coherence, in non-equilibrium processes20,21,22,23,24,25,26,27,28,29,30,31,32,33. Ref. 34 introduced an end-point measurement (EPM) approach, where the initial statistics of energy-change fluctuations is inferred from the knowledge of the initial state and the Hamiltonian of the quantum system.
Here, we show the intrinsic operational nature of the EPM approach by considering both the detailed and the integral form of the corresponding fluctuation theorem, and using them to characterise experimentally the entropy production associated to quantum coherence in the state of an open quantum system. In particular, we experimentally make use of a qubit encoded in the spin of an NV centre in diamond, subjected to both a continuous driving and environmental effects15,16,35. In our experiments, we observe a significant increase of the irreversibility of the resulting dissipative map that only originates from the presence of quantum coherence in the initial state of the NV centre. Remarkably, we show that measuring such a quantity provides a tight bound for the average heat exchanged by the system with the environment.
While being valid in principle for arbitrary dynamics, our results establish NV centres as valuable platforms for the exploration of energetics at the quantum level, thus enlarging the already prominent domain of their applications in quantum technologies36,37,38,39.
Results
EPM-based fluctuation theorem
In order to set the stage for the derivation of a generalized fluctuation theorem accounting for the effects of quantum coherence of the initial state, we briefly review the EPM scheme introduced in ref. 34. The results we provide in what follows hold for a generic open quantum system, and we thus consider an arbitrary completely-positive trace preserving (CPTP) map Φ. Let \({\rho }_{0}={{{\mathcal{P}}}}+\chi\) be the initial state of the system, which we have decomposed in its diagonal part \({{{\mathcal{P}}}}\) (expressed in the basis of the initial Hamiltonian \({H}_{{t}_{{\rm{in}}}}\)) and the traceless component χ that accounts for the quantum coherence. The EPM scheme prescribes to perform a single energy measurement, at the end of the process, and to associate to it the stochastic variables \({{\Delta }}{E}_{i,f}\equiv {E}_{f}^{{{{\rm{fin}}}}}-{E}_{i}^{{{{\rm{in}}}}}\) that encodes the energy fluctuations during the open dynamics. Here, \({E}_{k}^{{{{\rm{in}}}}({{{\rm{fin}}}})}\) denotes the eigenvalues of the initial (final) Hamiltonian. We thus introduce the probability distribution PEPM(ΔE) = ∑i,f PΓ(i, f)δ(ΔE − ΔEi,f) for the EPM scheme, where δ( ⋅ ) denotes the Dirac delta function and PΓ(i, f) are the joint probabilities of the energy-change in the forward trajectory Γ, defined as
Here \({{{\Pi }}}_{{{{\rm{j}}}}}^{{{{\rm{in}}}}({{{\rm{fin}}}})}\) is the projector on the j-th initial (final) energy eigenstate. Computing the characteristic function of the probability distribution PEPM(ΔE) leads immediately to an integral fluctuation relation. In fact, let us take the diagonal part \({{{\mathcal{P}}}}\) of the initial state ρ0 as the thermal state of the initial Hamiltonian \({H}_{{t}_{{\rm{in}}}}\) with inverse temperature β. One can then show that
where \({\rho }_{{{{\rm{th}}}}}^{{{{\rm{fin/in}}}}}\equiv {{{{\rm{e}}}}}^{-\beta {H}_{{t}_{{\rm{fin/in}}}}}/{Z}_{{\rm{fin/in}}}\) with \({Z}_{{\rm{fin/in}}}\equiv {{{\rm{tr}}}}({{{{\rm{e}}}}}^{-\beta {H}_{{t}_{{\rm{fin/in}}}}})\), and d is the dimension of the system’s Hilbert space. The second term in the right-hand-side of Eq. (2) showcases the contribution coming from the initial quantum coherence, while the first term represents a classical deviation from the Jarzynski’s equality. Such deviation would be present also in the absence of initial coherence due to the non-linearity of the EPM’s probability distribution for convex combination of states34. Notice that the non-linearity of the probability distribution stems from enforcing the positivity and reality of the joint distribution PΓ(i, f) in the case the initial state and the quantum system Hamiltonian are not commuting31,32.
Consider now an open quantum map Φ that admits a non-singular fixed point ρ* = Φ(ρ*). Using the results in40, we can define the corresponding backwards dynamics (cf. Methods – Derivation of the detailed fluctuation theorem). In this way, one can compare the joint probabilities \({P}_{{{\Gamma }}}(i,f)=p({E}_{i}^{{{{\rm{in}}}}})p({E}_{f}^{{{{\rm{fin}}}}})\equiv {p}_{i}^{{{{\rm{in}}}}}{p}_{f}^{{{{\rm{fin}}}}}\) and \({P}_{\widetilde{{{\Gamma }}}}(f,i)={{{\rm{tr}}}}({{{\Pi }}}_{f}^{{{{\rm{fin}}}}}{\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}){{{\rm{tr}}}}({{{\Pi }}}_{i}^{{{{\rm{in}}}}}\widetilde{{{\Phi }}}({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}))\equiv {\widetilde{p}}_{f}^{{{{\rm{in}}}}}{\widetilde{p}}_{i}^{{{{\rm{fin}}}}}\) for measuring the energy of the quantum system in the forward and backward trajectories, Γ and \(\widetilde{{{\Gamma }}}\) respectively. Notice that, in \({P}_{\widetilde{{{\Gamma }}}}(f,i)\), \({\rho }_{{{{\rm{in}}}}}^{{{{\rm{B}}}}}\) denotes the initial state of \(\widetilde{{{\Gamma }}}\), and \(\widetilde{{{\Phi }}}\) is the corresponding time-reversed map.
By extending the derivation of the Jarzynski equality2 by means of the Crooks’ formalism3,40, we consider the case in which (i) the initial quantum state of the forward dynamics is \({\rho }_{0}={\rho }_{{{{\rm{th}}}}}^{{{{\rm{in}}}}}+\chi\), i.e., ρ0 is written as the sum of a thermal state at inverse temperature β and the traceless component χ encoding the initial coherence in the energy basis, and (ii) the initial state of the backward quantum dynamics is the state that is thermal in the final Hamiltonian at the same inverse temperature β of the forward process, i.e., \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}={Z}_{{{{\rm{fin}}}}}^{-1}\exp [-\beta {H}_{{t}_{{{{\rm{fin}}}}}}]={\rho }_{{{{\rm{th}}}}}^{{{{\rm{fin}}}}}\). These assumptions allow to write the balance equation
where
and \({p}_{f}^{{{{\rm{fin}}}}}(\chi )\equiv {{{\rm{tr}}}}({{{\Pi }}}_{f}^{{{{\rm{fin}}}}}{{\Phi }}(\chi ))\) (cf. Methods – Derivation of the detailed fluctuation theorem). Accordingly, averaging Eq. (3) over the forward probability distribution, one obtains
Eq. (3) and Eq. (5) are the detailed and the integral forms of the EPM’s fluctuation theorem, respectively.
It is worth noting that, by construction, \({p}_{f}^{{{{\rm{fin}}}}}(\chi )\) can be negative as it operates on the quantum coherence operator χ, which is Hermitian but traceless. However, the quantity ΔΣi,f is well-defined, as the term inside the logarithm is always non-negative. In fact, \({p}_{f}^{{{{\rm{fin}}}}}({\rho }_{{{{\rm{th}}}}}^{{{{\rm{in}}}}})+{p}_{f}^{{{{\rm{fin}}}}}(\chi )={p}_{f}^{{{{\rm{fin}}}}}({\rho }_{0})\ge 0\). In addition, thanks to the assumption (ii), ΔΣi,f depends only on the forward dynamics and it satisfies a fluctuation theorem by its own (cf. Methods – Derivation of the integral fluctuation theorem), i.e.,
Resorting to the Jensen’s inequality, one has \(\left\langle {{\Delta }}\Sigma \right\rangle \ge 0\). ΔΣ encodes the entropic contribution of the initial quantum coherence of the system, and we thus identify it as the coherence-affected irreversible entropy production for a non-equilibrium dynamical process. At the same time, the quantity Δσ represents a completely classical contribution to the entropy production that comes from adopting the EPM formalism, namely from the extra uncertainty implied by the factorization condition in Eq. (1) (see also discussion in ref. 34).
The experimental system
We consider the spin qubit associated to a negatively charged NV centre –a localized impurity in a diamond lattice based on a nitrogen substitutional atom next to a vacancy– which forms an electronic spin S = 1 in its orbital ground state41,42,43,44. A magnetic bias field aligned with Sz removes the degeneracy of the spin eigenstates, so as to allow for the selective coherent manipulation of the transition \(\left\vert {S}_{z}=0\right\rangle \leftrightarrow \left\vert {S}_{z}=+1\right\rangle\). The Hamiltonian H of this effective two-level system is determined by a continuous nearly resonant microwave field and, in the frame rotating at the microwave frequency, the Hamiltonian is \(H=\hslash \omega (\cos \alpha \,{\tilde{\sigma }}_{z}-\sin \alpha \,{\tilde{\sigma }}_{x})/2\), where \(\omega =\sqrt{{{{\Omega }}}^{2}+{\delta }^{2}}\) and \(\tan \alpha =-{{\Omega }}/\delta\); Ω denotes the bare Rabi frequency and δ ∈ [0, Ω] is the microwave detuning. We have used the tilde for the Pauli matrices in view of the change of basis to the Hamiltonian eigenstates, i.e., \(\{\left\vert 0\right\rangle ,\left\vert 1\right\rangle \}\equiv \{\cos \frac{\alpha }{2}\,\left\vert {S}_{z}=0\right\rangle -\sin \frac{\alpha }{2}\,\left\vert {S}_{z}=+1\right\rangle ,\sin \frac{\alpha }{2}\,\left\vert {S}_{z}=0\right\rangle +\cos \frac{\alpha }{2}\,\left\vert {S}_{z}=+1\right\rangle \}\) with eigenvalues ± ℏω/2. In this new basis, the Hamiltonian becomes H = ωσz/2, where ℏ is set to 1 from here on.
The qubit is governed by an alternated sequence of unitary and non-unitary (controlled-dissipative) dynamics, as follows. The system is repeatedly subjected to a sequence of pulses, occurring regularly at time intervals τ. Among two consecutive pulses, the evolution of the NV centre is unitary and described by the operator \(U\equiv \exp [-{{{\rm{i}}}}H\tau ]\). As depicted in Fig. 1, the NV spin also undergoes open dynamics due to its interaction with a train of short laser pulses with a duration tL that is much shorter than the characteristic time-scale of the unitary dynamics (tL ≪ 2π/ω). The short laser pulses trigger cycles of spin preserving and non-preserving transitions between different orbital levels15,45,46. This entails non-unitary dynamics that project the state of the system into the eigenstates of \({\tilde{\sigma }}_{z}\) and partially transfer the spin population to \(\left\vert {S}_{z}=0\right\rangle\). Such spin amplitude damping along the \({\tilde{\sigma }}_{z}\) axis, also known as optical pumping, is caused by the spin non-preserving transitions, and can be modelled as a controlled dissipative channel toward \(\left\vert {S}_{z}=0\right\rangle\)15,35. The overall dynamics takes the NV centre into an asymptotic fixed point ρ* that we can use to define the backward dynamics [see Supplementary Note 1 for more details on the modelling of the NV centre dynamics in the experiment]. In the experiments shown throughout this work, we set α = π/4 (i.e., δ = −Ω), τω ≃ (2π)0.9, and τ = 190 ns.
Coherence-affected entropy production
We now show the results obtained with the NV qubit subjected to the dissipative dynamics introduced above. Here, our aim is to characterize the thermodynamic role of initial quantum coherence in terms of the coherence-affected entropy production. We have thus performed a series of experiments to determine both the EPM probability distribution of energy-change fluctuations as well as the usual TPM one.
At the beginning of each experimental realization, the electronic spin is initialized into the \(\left\vert {S}_{z}=0\right\rangle\) eigenstate of the spin operator Sz via optical spin pumping under long laser excitation. Then, the system is brought in each of the four different pure states that correspond to the eigenvectors of σz and σy, by applying rotation gates (on-resonant microwave pulse) to the state \(\left\vert {S}_{z}=0\right\rangle\). After n short laser pulses, with n ∈ [0, N], we measure the energy of the system in the Hamiltonian basis σz. To achieve this, we optically readout the intensity of the NV photo-luminescence (PL). The PL intensity (averaged over ~ 106 repetitions of the experiment) determines the probability for the system to be in the energy eigenstate \(\left\vert 1\right\rangle\) (see Methods for a more detailed description of the readout process). This approach allows us to achieve an experimental estimate of \({p}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}(\rho )={{{\rm{tr}}}}({{{\Pi }}}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}{{\Phi }}(\rho ))\), where ρ is one of the eigenstates of either σz or σy. From this, we obtain the probability for f = 0 as \(1-{p}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}(\rho )\). As we collect data from several experiments (in which the system is initialized in one of the four pure qubit-states \(\{\left\vert 0\right\rangle ,\left\vert 1\right\rangle ,{\left\vert +\right\rangle }_{y},{\left\vert -\right\rangle }_{y}\}\)), it is a matter of data processing to compute the EPM and TPM statistics for every classical mixture of such states (for more details see Supplementary Note 2). In particular, looking at a Bloch sphere representation for the qubit, this implies that we are able to obtain the statistics of the quantum process for any initial state included in the y − z equatorial plane of the Bloch sphere, i.e. any initial state expressed as a convex combination of the states \(\left\vert 0\right\rangle ,\left\vert 1\right\rangle ,{\left\vert +\right\rangle }_{y},{\left\vert -\right\rangle }_{y}\). Specifically, in the following we show experimental results corresponding to the initial states (expressed in the Hamiltonian basis)
with p ∈ [0,1], such that ρ0 is the convex mixture of \(\left\vert 1\right\rangle \,\left\langle 1\right\vert\) and \({\left\vert +\right\rangle }_{y}\,\left\langle +\right\vert\) with probability p and 1 − p, respectively. We also need to consider data associated to the state \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}=\frac{1}{2}({\mathbb{1}}+p{\sigma }_{z})\). Notice that in ref. 15 we reported similar experimental results, but only considering diagonal states \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}\). While the working point defined by the choice of experimental parameters (δ, Ω, τ, tL) is unique of the experiments reported here, our investigation goes significantly beyond the context set in ref. 15 in light of the fact that we have addressed the case of coherent (non-diagonal) states.
The first quantity we are interested in characterizing experimentally is the coherence-affected entropy production encoded in the average of ΔΣ as given in Eq. (4). Notice that this average is defined solely in terms of the forward trajectory probability; thus, we can fully characterize it by resorting to the experimental data acquired during the forward dynamics. The results of our experiments are shown in Fig. 2: in panel (a) we present the experimental values of the stochastic quantum entropy production ΔΣi,f, and in panel (b) we show the experimental verification of the fluctuation theorem in Eq. (6). Instead, in panel (c), the behaviour of 〈ΔΣ〉 is plotted as a function of the number of laser pulses. While the error bars are quite large, it can be observed how the corresponding experimental points nicely follow the theoretical predictions and how the experimental data show a positive coherence-affected entropy production.
Another quantity that can be investigated directly from the available data on the forward dynamics is the average of the energy-change fluctuations ΔE. In this regard, it is worth noting that this quantity is identified with the average work when considering time-dependent unitary processes like in the Jarzynski’s original work2. In our case, as the Hamiltonian is time-independent, we can unambiguously interpret this quantity as the average heat that the system exchanges with its environment, by means of the open dynamics to which the NV centre is subjected. The average of the stochastic variable ΔEi,f = Ef − Ei in the EPM approach is related to the TPM scheme via the following relation:
The second term on the right-hand side of Eq. (8) represents a contribution to the average heat ascribable to the quantum coherence of the initial state, which is not deleted by a first energy measurement over ρ0. Figure 3 displays the three quantities entering in Eq. (8), by comparing theoretical expectations with the experimental results. It is shown that the experimental data are able to discern the coherence contribution to the heat exchanged between the system and the environment due to the pulsed dynamics. Among the protocols that allows to account for quantum features in energy fluctuations20,21,22,23,27,28,29,31,34, the EPM scheme requires only a final energy measurement at the time t and does not require any knowledge of the quantum map Φ that models the forward dynamics. This makes the method particularly suitable for application in open quantum systems.
Finally, we have also verified the validity of the full EPM fluctuation theorem in Eq. (5) as well as its consequences for the expectation values of the involved thermodynamic quantities. In principle, this requires having access to the backward trajectories of the system, so as to characterize Δσi,f in Eq. (4). This can be problematic when the system dynamics is non-unitary: while the backward trajectories can be easily simulated numerically, implementing them at the experimental level is not currently possible with our set-up. However, for the range of experimental parameters, the choice of the initial state of the backward process and the form \(\widetilde{{{\Phi }}}\) obtained from40 thanks to the existence of a non-singular fixed point of the dynamics, the backward dynamics are such that \({\widetilde{p}}_{{{{\rm{j}}}}}^{{{{\rm{fin}}}}}({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}})={p}_{{{{\rm{j}}}}}^{{{{\rm{fin}}}}}({\rho }_{{{{\rm{th}}}}}^{{{{\rm{in}}}}})\). Therefore, by assuming that this property holds and thus that we can estimate 〈Δσ〉 from the sole data of the forward trajectories, in Fig. 3b we show the experimental verification of the integral form of the EPM fluctuation theorem in Eq. (5).
The application of the Jensen inequality to Eqs. (2) and (5) leads to
with \({{{{\mathcal{G}}}}}_{{{{\rm{EPM}}}}}=2{{{\rm{tr}}}}\left({\rho }_{{{{\rm{th}}}}}^{{{{\rm{fin}}}}}{{\Phi }}({\rho }_{{{{\rm{th}}}}}^{{{{\rm{in}}}}}+\chi )\right)\) denoting the EPM characteristic function34. In Fig. 4, the theoretical expectations of these quantities are compared with the corresponding experimental results. For the initial state ρ0 considered for this figure (i.e., the state in Eq. (7) with p = 0.38), the results obtained by using the characteristic function of the EPM approach are quite distinct from the ones coming from the integral fluctuation theorem (5). Specifically, for N ∈ [0, 10] and for the initial states accessible from the experimental data, the bound on 〈ΔE〉 derived from the integral fluctuation theorem is tighter than the one resulting from \({{{{\mathcal{G}}}}}_{{{{\rm{EPM}}}}}\). In addition, we can also see in Fig. 4 that the inequality in Eq. (9) is almost saturated, meaning that the average heat exchange from the EPM method offers even a good estimate of the sum 〈ΔΣ + Δσ〉 in the regime experimentally analyzed.
Discussion
We have used the end-point measurement (EPM) approach for the characterization of energy-change fluctuations arising from a non-equilibrium process, with the aim to quantify the contribution to the entropy production that is originated by the presence of quantum coherence in the initial density matrix of the quantum system under scrutiny.
The results obtained from the EPM scheme have to be compared with the ones arising from the two-point measurement (TPM) scheme. While in the TPM an initial projective measurement cancels out any effect of the quantum coherence of the initial state in the basis of the measured observable, the EPM preserves such coherence. In order to compare these two schemes, in this work we followed the standard set-up of the TPM scheme while allowing for initial coherence in the energy basis applying the EPM approach. As a general remark, notice that the coherence of an initial quantum state crucially influences the system evolution, and consequently the physical quantities associated with the system dynamics, such as the quantum entropy production.
Our formalism gives a physical characterization of the entropy production associated with an arbitrary preparation of the system. Besides the quantum entropy originated during the dynamics, our results clearly single-out the presence of a change of entropy production due to the absence of the first measurement, and of an entropy production solely depending on the quantum coherence of the initial state. By focusing the attention on the latter contribution, we have shown that the EPM formalism provides –quite naturally– the fluctuation theorem for the quantum entropy production linked to the quantum coherence of the initial state.
The operational nature of our approach has enabled a successful experimental assessment of such coherence-affected entropy production in a solid state platform, at room-temperature and inherently dissipative, where an NV spin qubit is addressed by a sequence of laser pulses so that it undergoes a controlled dissipative dynamics.
Our study grounds the EPM approach as a powerful framework to assess the role of quantum coherence in the energetics of quantum systems and devices47. Specifically, our findings add a crucial ingredient for the analysis of the thermodynamic role of quantum coherence and allow the characterization of the introduced entropy production due to quantum coherence. It seems natural to associate the entropy production heralded by quantum coherence in the initial state, to the system dissipation. While our work has provided robust empirical and experimental evidence for this connection, further investigations are required to formalise and establish it definitively.
Methods
Here we report a detailed derivation of the detailed and integral fluctuation theorems analyzed in the section Results.
Derivation of the detailed fluctuation theorem
Let us consider a CPTP map, written in Kraus representation as \({{\Phi }}(\bullet )={\sum }_{\alpha }{K}_{\alpha }\bullet {K}_{\alpha }^{{\dagger} }\) with {Kα} the set of Kraus operators, that allows a non-singular fixed point ρ* such that ρ* = Φ(ρ*). As ρ* is non-singular, following Crooks’ formalism40 we can define the time-reversed map as
where \({\widetilde{K}}_{\alpha }\equiv {({\rho }^{* })}^{1/2}{K}_{\alpha }^{{\dagger} }{({\rho }^{* })}^{-1/2}\).
Then, in accordance with the EPM formalism34, we can introduce the expression of the joint probabilities \({P}_{{{\Gamma }}}({E}_{i}^{{{{\rm{in}}}}},{E}_{f}^{{{{\rm{fin}}}}})\equiv {P}_{{{\Gamma }}}(i,f)\) and \({P}_{\widetilde{{{\Gamma }}}}({E}_{\ell }^{{{{\rm{fin}}}}},{E}_{k}^{{{{\rm{in}}}}})\equiv {P}_{\widetilde{{{\Gamma }}}}(\ell ,k)\), respectively, for the quantum trajectories Γ and \(\widetilde{{{\Gamma }}}\) of the forward and backward process, respectively. Let us observe that the latter is obtained by reversing the arrow of time by means of a transformation that obeys the time-reversal symmetry. Specifically, one has
and
where, we recall, ρfin ≡ Φ[ρ0] and \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{fin}}}}}\equiv \widetilde{{{\Phi }}}[{\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}]\) with \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}\) denoting the initial state of the backward process and, without loss of generality, we have assumed that each projector on the energy eigenstates is invariant to the application of the time-reversal transformation. As a result, the expression of the detailed balance equation for the trajectories Γ and \(\widetilde{{{\Gamma }}}\) is formally equal to
Before proceeding, it is worth observing that if (i) we apply the EPM scheme for characterising energy-change fluctuations and (ii) we choose the initial state of the backward process as the time-reversal of the final state of the forward process, then \({\widetilde{p}}_{\ell }^{{{{\rm{in}}}}}={p}_{f}^{{{{\rm{fin}}}}}\) for ℓ = f. As a matter of fact,
where \(\widetilde{(\bullet )}\equiv {{\Theta }}(\bullet ){{{\Theta }}}^{{\dagger} }\) with Θ denoting the time-reversal operator. The latter, by construction, is anti-unitary, i.e., it is an anti-linear operator and satisfies the relations \({{{\Theta }}}^{{\dagger} }{{\Theta }}={{\Theta }}{\Theta }^{{\dagger} }={\mathbb{I}}\). In this context, we obtain the fluctuation relation
where we have identified
This relation encodes information only on the initial stochastic quantum entropy production due to the open system dynamics. In fact, for the special case in which the dynamics is unitary, and under the assumption of micro-reversibility, i.e., \({{\Theta }}\,{{{\mathcal{U}}}}(\lambda ){{{\Theta }}}^{{\dagger} }={{{{\mathcal{U}}}}}^{{\dagger} }(\widetilde{\lambda })\) with λ(t) generic time-dependent transformation such that the system Hamiltonian is invariant under time-reversal, one finds that \({{{{\rm{e}}}}}^{-{{\Delta }}{\overline{\sigma }}_{i}^{{{{\rm{in}}}}}}=1\,\,\forall i\), i.e.,
Notice also that, in determining the integral form of Eq. (15), one gets \({\langle {{{{\rm{e}}}}}^{{{\Delta }}{\overline{\sigma }}^{{{{\rm{in}}}}}}\rangle }_{{{\Gamma }}}=1\) with
where, we recall, \({p}_{i}^{{{{\rm{in}}}}}\) and \({\widetilde{p}}_{i}^{{{{\rm{fin}}}}}\) are the probabilities to measure the i-th energy value of the system, respectively, at the initial and final time instants of the forward and backward process. Here, S(q∣∣p) is the classical relative entropy between the two probability distributions q and p, and thus it naturally corresponds to a measure of how far is the final state of the inverse quantum dynamics from the initial quantum state.
Coming back to the derivation of Eq. (3) in the main text, let us assume that (i) the initial quantum state ρ0 of the forward dynamics has thermal populations but also non-zero coherence terms (in the energy basis of the system), and (ii) the initial quantum state \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}\) of the backward quantum dynamics has thermal populations (with respect to the final Hamiltonian) and once again non-vanishing off-diagonal elements. In other terms,
where \({{{\rm{tr}}}}({\chi }^{{{{\rm{in}}}}})={{{\rm{tr}}}}({\chi }^{{{{\rm{fin}}}}})=0\), \({Z}_{{\rm{fin/in}}}\equiv {{{\rm{tr}}}}(\exp [-\beta {H}_{{t}_{{\rm{fin/in}}}}])\) and the system Hamiltonian H is not necessarily assumed as a time-independent operator. Thus, by substituting Eqs. (19) and (20) into (13) with k = i and ℓ = f, one finds
where
\({p}_{f}^{{{{\rm{fin}}}}}({{{\mathcal{A}}}})\equiv {{{\rm{tr}}}}\left({{{\Pi }}}_{f}^{{{{\rm{fin}}}}}{{\Phi }}({{{\mathcal{A}}}})\right)\), and \({\widetilde{p}}_{i}^{{{{\rm{fin}}}}}({{{\mathcal{A}}}})\equiv {{{\rm{tr}}}}\left({{{\Pi }}}_{i}^{{{{\rm{fin}}}}}\widetilde{{{\Phi }}}({{{\mathcal{A}}}})\right)\) with \({{{\mathcal{A}}}}\) a generic linear operator. It is worth noting that \({p}_{f}^{{{{\rm{fin}}}}}({{{\mathcal{A}}}})\) and \({\widetilde{p}}_{i}^{{{{\rm{fin}}}}}({{{\mathcal{A}}}})\) denote the probability to measure the f-th and i-th final energy values of the quantum system in the forward and backward process, respectively, conditioned to have evolved the thermal contribution of the initial state (without coherence terms in the energy eigenbasis). Finally, Eq. (21) can be simplified as
by introducing the quantities
Let us observe that, if no quantum coherence is present neither in ρ0 nor in \({\rho }_{{{{\rm{B}}}}}^{{{{\rm{in}}}}}\) (i.e., χin = χfin = 0), then ΔΣi,f = 0. Thus, ΔΣi,f can be considered as a correction, due to initial coherence in the energy basis of the system, to the entropy difference Δσi,f obtained by propagating initial thermal states in the forward and backward process, respectively.
The form of the detailed fluctuation theorem used in the main text is a particular case of Eq. (26) where χfin is assumed to be vanishing. This choice is motivated by our aim to consider the minimal modification to the “Jarzynski set-up”, respect to which only coherence in the initial state of the forward dynamics is added.
Derivation of the integral fluctuation theorem for ΔΣ
From Eq. (26) the integral fluctuation theorem of Eq. (5) can be easily obtained.
Instead, for what concerns the integral fluctuation theorem involving the sole coherence induced entropy production, let us consider the case in which χfin = 0. Upon substitution, one has
Thus, taking the average over the EPM probability distribution of the forward process Γ, we can conclude that
where we have used the fact that \({p}_{f}^{{{{\rm{fin}}}}}({\rho }_{{{{\rm{th}}}}}^{{{{\rm{in}}}}})\) is itself a probability distribution normalized to 1.
Description of the readout process
As described in the main text, at the end of the protocol we want to measure the probability \({p}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}(\rho )\equiv {{{\rm{tr}}}}({{{\Pi }}}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}{{\Phi }}(\rho ))\), where ρ is one of the four different pure states corresponding to the eigenvectors of σz and σy, and Φ is the map associated to a train of n short laser pulses (n ∈ [0, N]) with a unitary evolution U between consecutive pulses. First, the electronic spin is prepared in the desired state ρ by means of a rotation gate (on-resonant microwave pulse) to the state \(\left\vert {S}_{z}=0\right\rangle\). After n short laser pulses, we apply another rotation gate (on-resonant microwave pulse) such that \({\tilde{\sigma }}_{z}\to {\sigma }_{z}\). Finally, we measure the photo-luminescence (PL), where the photon counts are averaged over ~ 106 repetitions of the experiment. Then, this PL intensity is normalized with respect to the reference PL intensity of the eigenstates \(\left\vert {S}_{z}=0\right\rangle\) and \(\left\vert {S}_{z}=+1\right\rangle\) (also measured after each repetition of the experiment).
Let us name PL0 and PL1 the reference PL intensity values for the states \(\left\vert {S}_{z}=0\right\rangle\) and \(\left\vert {S}_{z}=+1\right\rangle\), respectively. Thus, the normalized signal s we obtain from our experiments corresponds to s = (PL − PL0)/(PL1 − PL0), which is the probability to find the system in the state \(\left\vert {S}_{z}=+1\right\rangle\). Hence, \(s=\frac{1}{2}(1+\langle {\tilde{\sigma }}_{z}\rangle )=\left\langle {S}_{z}=+1\right\vert {\tilde{\rho }}^{{{{\rm{fin}}}}}\left\vert {S}_{z}=+1\right\rangle\), where \({\tilde{\rho }}^{{{{\rm{fin}}}}}\) is the state of the spin qubit at the very end of the protocol. Notice though that the final rotation gate needs the change of basis \({\tilde{\sigma }}_{z}\to {\sigma }_{z}\). This rotation is described by the unitary operator \({R}_{y}(\alpha )\equiv \exp [-{{{\rm{i}}}}{\sigma }_{y}\alpha /2]\). Therefore, \(s=\left\langle {S}_{z}=+1\right\vert {R}_{y}(\alpha ){\rho }^{{{{\rm{fin}}}}}{R}_{y}(-\alpha )\left\vert {S}_{z}=+1\right\rangle =\left\langle 1\right\vert {\rho }^{{{{\rm{fin}}}}}\left\vert 1\right\rangle\), where we have used that \({R}_{y}(-\alpha )\left\vert {S}_{z}=+1\right\rangle =\left\vert 1\right\rangle\) and ρfin = Φ(ρ) is the state of the system before the final rotation. As a result, the probability that we want to measure is provided by the normalized signal \(s={{{\rm{Tr}}}}({{\Phi }}(\rho )\left\vert 1\right\rangle \left\langle 1\right\vert )={p}_{{{{\rm{f = 1}}}}}^{{{{\rm{fin}}}}}(\rho )\).
Data availability
Data are available upon request from the authors.
Code availability
Codes are available from the authors upon reasonable requests.
References
Landi, G. T. & Paternostro, M. Irreversible entropy production: from classical to quantum. Rev. Mod. Phys. 93, 035008 (2021).
Jarzynski, C. Nonequilibrium equality for free energy differences. Phys. Rev. Lett. 78, 2690–2693 (1997).
Crooks, G. E. Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. Phys. Rev. E 60, 2721–2726 (1999).
Campisi, M., Hänggi, P. & Talkner, P. Colloquium: Quantum fluctuation relations: Foundations and applications. Rev. Mod. Phys. 83, 771 (2011).
Esposito, M., Harbola, U. & Mukamel, S. Nonequilibrium fluctuations, fluctuation theorems, and counting statistics in quantum systems. Rev. Mod. Phys. 81, 1665–1702 (2009).
Talkner, P., Lutz, E. & Hänggi, P. Fluctuation theorems: Work is not an observable. Phys. Rev. E 75, 050102 (2007).
Gherardini, S., Müller, M., Trombettoni, A., Ruffo, S. & Caruso, F. Reconstructing quantum entropy production to probe irreversibility and correlations. Quantum Sci. Technol. 3, 035013 (2018).
Manzano, G., Horowitz, J. M. & Parrondo, J. M. R. Quantum fluctuation theorems for arbitrary environments: adiabatic and nonadiabatic entropy production. Phys. Rev. X 8, 031037 (2018).
Batalhão, T. B. et al. Experimental reconstruction of work distribution and study of fluctuation relations in a closed quantum system. Phys. Rev. Lett. 113, 140601 (2014).
Batalhão, T. B. et al. Irreversibility and the arrow of time in a quenched quantum system. Phys. Rev. Lett. 115, 190601 (2015).
An, S. et al. Experimental test of the quantum Jarzynski equality with a trapped-ion system. Nat. Phys. 11, 193–199 (2015).
Smith, A. et al. Verification of the quantum nonequilibrium work relation in the presence of decoherence. New J. Phys. 20, 013008 (2018).
Xiong, T. P. et al. Experimental verification of a jarzynski-related information-theoretic equality by a single trapped ion. Phys. Rev. Lett. 120, 010601 (2018).
Zhang, Z. et al. Experimental demonstration of work fluctuations along a shortcut to adiabaticity with a superconducting Xmon qubit. New J. Phys. 20, 085001–13 (2018).
Hernández-Gómez, S. et al. Experimental test of exchange fluctuation relations in an open quantum system. Phys. Rev. Research 2, 023327 (2020).
Hernández-Gómez, S., Staudenmaier, N., Campisi, M. & Fabbri, N. Experimental test of fluctuation relations for driven open quantum systems with an NV center. New J. Phys. 23, 065004 (2021).
Cimini, V. et al. Experimental characterization of the energetics of quantum logic gates. npj Quantum Inf. 6, 1–8 (2020).
Ribeiro, P. H. S. et al. Experimental study of the generalized jarzynski fluctuation relation using entangled photons. Phys. Rev. A 101, 052113 (2020).
Aguilar, G. H. et al. Two-point measurement of entropy production from the outcomes of a single experiment with correlated photon pairs. Phys. Rev. A 106, L020201 (2022).
Allahverdyan, A. E. Nonequilibrium quantum fluctuations of work. Phys. Rev. E 90, 032137 (2014).
Deffner, S., Paz, J. P. & Zurek, W. H. Quantum work and the thermodynamic cost of quantum measurements. Phys. Rev. E 94, 010103 (2016).
Lostaglio, M. Quantum fluctuation theorems, contextuality, and work quasiprobabilities. Phys. Rev. Lett. 120, 040602 (2018).
Santos, J., Celeri, L., Landi, G. & Paternostro, M. The role of quantum coherence in non-equilibrium entropy production. npj Quantum Inf. 5, 23 (2019).
Kwon, H. & Kim, M. S. Fluctuation theorems for a quantum channel. Phys. Rev. X 9, 031029 (2019).
Wu, K.-D. et al. Experimentally reducing the quantum measurement back action in work distributions by a collective measurement. Sci. Adv. 5, eaav4944 (2019).
Wu, K.-D. et al. Minimizing backaction through entangled measurements. Phys. Rev. Lett. 125, 210401 (2020).
Sone, A., Liu, Y.-X. & Cappellaro, P. Quantum jarzynski equality in open quantum systems from the one-time measurement scheme. Phys. Rev. Lett. 125, 060602 (2020).
Micadei, K., Landi, G. T. & Lutz, E. Quantum fluctuation theorems beyond two-point measurements. Phys. Rev. Lett. 124, 090602 (2020).
Micadei, K. et al. Experimental validation of fully quantum fluctuation theorems using dynamic Bayesian networks. Phys. Rev. Lett. 127, 180603 (2021).
Yada, T., Yoshioka, N. & Sagawa, T. Quantum fluctuation theorem under quantum jumps with continuous measurement and feedback. Phys. Rev. Lett. 128, 170601 (2022).
Lostaglio, M. et al. Kirkwood-dirac quasiprobability approach to quantum fluctuations: theoretical and experimental perspectives. arXiv https://doi.org/10.48550/arXiv.2206.11783 (2022).
Hernández-Gómez, S. et al. Projective measurements can probe non-classical work extraction and time-correlations. arXiv https://doi.org/10.48550/arXiv.2207.12960 (2022).
Bellini, M. et al. Demonstrating quantum microscopic reversibility using coherent states of light. Phys. Rev. Lett. 129, 170604 (2022).
Gherardini, S., Belenchia, A., Paternostro, M. & Trombettoni, A. End-point measurement approach to assess quantum coherence in energy fluctuations. Phys. Rev. A 104, L050203 (2021).
Hernández-Gómez, S. et al. Autonomous dissipative Maxwell’s demon in a diamond spin qutrit. PRX Quantum 3, 020329 (2022).
Maze, J. R. et al. Nanoscale magnetic sensing with an individual electronic spin in diamond. Nature 455, 644 (2008).
Rondin, L. et al. Magnetometry with nitrogen-vacancy defects in diamond. Rep. Prog. Phys. 77, 056503 (2014).
Degen, C. L., Reinhard, F. & Cappellaro, P. Quantum sensing. Rev. Mod. Phys. 89, 035002 (2017).
Hernández-Gómez, S. & Fabbri, N. Quantum control for nanoscale spectroscopy with diamond nitrogen-vacancy centers: a short review. Front. Phys. https://doi.org/10.3389/fphy.2020.610868 (2021).
Crooks, G. E. Quantum operation time reversal. Phys. Rev. A 77, 034101 (2008).
Gruber, A. et al. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science 276, 2012–2014 (1997).
Jelezko, F., Gaebel, T., Popa, I., Gruber, A. & Wrachtrup, J. Observation of coherent oscillations in a single electron spin. Phys. Rev. Lett. 92, 076401 (2004).
Doherty, M. W. et al. The nitrogen-vacancy colour centre in diamond. Phys. Rep. 528, 1–45 (2013).
Aharonovich, I. & Neu, E. Diamond nanophotonics. Adv. Opt. Mater. 2, 911–928 (2014).
Wolters, J., Strauß, M., Schoenfeld, R. S. & Benson, O. Quantum zeno phenomenon on a single solid-state spin. Phys. Rev. A 88, 020101 (2013).
Manson, N. B., Harrison, J. P. & Sellars, M. J. Nitrogen-vacancy center in diamond: Model of the electronic structure and associated dynamics. Phys. Rev. B 74, 104303 (2006).
Gianani, I. et al. Diagnostics of quantum-gate coherences via end-point-measurement statistics. arXiv https://doi.org/10.48550/arXiv.2209.02049 (2022).
Acknowledgements
We acknowledge support from the European Union’s Horizon 2020 FET-Open project TEQ (766900), the Horizon Europe EIC Pathfinder project QuCoM (Grant Agreement No. 101046973), the Leverhulme Trust Research Project Grant UltraQuTe (grant RGP-2018-266), the Royal Society Wolfson Fellowship (RSWF/R3/183013), the UK EPSRC (grant EP/T028424/1), the Department for the Economy Northern Ireland under the US-Ireland R&D Partnership Programme, the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) project number BR 5221/4-1, the MISTI Global Seed Funds MIT-FVG Collaboration Grant “Non-Equilibrium Thermodynamics of Dissipative Quantum Systems (NETDQS)”, the Blanceflor Foundation through the project “The theRmodynamics behInd thE meaSuremenT postulate of quantum mEchanics (TRIESTE)”, the CNR-FOE-LENS-2020, the Horizon Europe RIA project MUQUABIS GA n. 101070546, and the PNRR MUR project PE0000023-NQSTI. We also acknowledge support from the European Union’s Next Generation EU Programme with the I- PHOQS Infrastructure IR0000016 ID D2B8D520 CUP B53C22001750006 “Integrated infrastructure initiative in Photonic and Quantum Sciences”, and from Fondazione CR Firenze through the SALUS project.
Author information
Authors and Affiliations
Contributions
N.F. and S.H.-G. designed and engineered the experimental protocol; S.H.-G. carried out the experiment and collected the data; S.G., A.B., A.T., and M.P. developed the theory and carried out the theoretical analysis; S.H.-G., S.G., and A.B. analyzed the data; All authors have contributed equally to the writing of the manuscript. S.H.-G., S.G., and A.B. are co-first authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hernández-Gómez, S., Gherardini, S., Belenchia, A. et al. Experimental signature of initial quantum coherence on entropy production. npj Quantum Inf 9, 86 (2023). https://doi.org/10.1038/s41534-023-00738-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41534-023-00738-0