Probing quantum fluctuation theorems in engineered reservoirs

TL;DR

This paper proposes a realistic platform using an engineered reservoir and a two-level system to experimentally test quantum fluctuation theorems, revealing quantum effects in entropy production and addressing detection inefficiencies.

Contribution

It introduces a practical method to probe quantum fluctuation theorems in open systems, including correction techniques for detection inefficiencies.

Findings

Quantum component in entropy production detected

Photon detection inefficiency correction demonstrated

Platform enables experimental exploration of quantum thermodynamics

Abstract

Fluctuation Theorems are central in stochastic thermodynamics, as they allow for quantifying the irreversibility of single trajectories. Although they have been experimentally checked in the classical regime, a practical demonstration in the framework of quantum open systems is still to come. Here we propose a realistic platform to probe fluctuation theorems in the quantum regime. It is based on an effective two-level system coupled to an engineered reservoir, that enables the detection of the photons emitted and absorbed by the system. When the system is coherently driven, a measurable quantum component in the entropy production is evidenced. We quantify the error due to photon detection inefficiency, and show that the missing information can be efficiently corrected, based solely on the detected events. Our findings provide new insights into how the quantum character of a physical system impacts its thermodynamic evolution.

Figures

Equations

$${\Delta }_{\text{i}}S=(W-\Delta F)/T.$$

$$⟨e^{-W[\gamma ]/k_{\text{B}}T}{⟩}_{\gamma }=e^{-\Delta F/k_{\text{B}}T}$$

$$\dot{\rho }=-\frac{i}{\hslash }\left[H_0+H_{\text{d}}(t),\rho \right]+{\Gamma }_{+}\mathcal{L}[{\sigma }^{†}]\rho +{\Gamma }_{-}\mathcal{L}[\sigma ]\rho ,$$

$$p_{\mathcal{K}(t_n)}=⟨{\psi }_{\gamma }(t_n)|M_{\mathcal{K}(t_n)}^{†}{M}_{\mathcal{K}(t_n)}|{\psi }_{\gamma }(t_n)⟩,\mathcal{K}\in \{0,1,2\}.$$

$$\delta W_{\gamma }(t_n)=⟨{\psi }_{\gamma }(t_n)|d{H}_{\text{d}}(t_n)|{\psi }_{\gamma }(t_n)⟩,\mathcal{K}(t_n)=0.$$

$$\displaystyle\delta Q_{\text{cl},\gamma}(t_{n})=\left\{\begin{array}[]{c}-\hbar\omega_{1}(t_{n}),\quad{\cal K}(t_{n})=1\\ \hbar\omega_{1}(t_{n}),\quad{\cal K}(t_{n})=2.\end{array}\right.$$

$$\displaystyle\delta Q_{\text{q},\gamma}(t_{n})=\left\{\begin{array}[]{c}\hbar\omega_{1}(t_{n})|\langle g|\psi_{\gamma}(t_{n})\rangle|^{2},\quad{\cal K}_{\gamma}(t_{n})=1\\ -\hbar\omega_{1}(t_{n})|\langle e|\psi_{\gamma}(t_{n})\rangle|^{2},\quad{\cal K}_{\gamma}(t_{n})=2,\end{array}\right.$$

$$\delta Q_{\text{q},\gamma }(t_n)$$

$$\Delta U[\gamma ]=U_{\gamma }(t_f)-U_{\gamma }(t_i)=W[\gamma ]+Q_{\text{cl}}[\gamma ]+Q_{\text{q}}[\gamma ],$$

$${\Delta }_{\text{i}}s[\gamma ]=\frac{1}{T_{\text{eff}}}(W[\gamma ]-\Delta F)+\frac{Q_{\text{q}}[\gamma ]}{T_{\text{eff}}}.$$

$${\mathcal{E}}_0^{(\eta )}[{\rho }_{\gamma }]$$

$$\displaystyle\delta Q_{\text{cl}}^{\eta}(t_{n})=\left\{\begin{array}[]{c}-\hbar\omega_{1}(t_{n}),\quad{\cal K}(t_{n})=1\\ \hbar\omega_{1}(t_{n}),\quad{\cal K}(t_{n})=2.\end{array}\right..$$

$$⟨e^{-{\Delta }_{\text{i}}{s}^{\eta }[\gamma ]/k_{\text{B}}-{\sigma }_{\eta }[\gamma ]}{⟩}_{\gamma }=1.$$

$$\displaystyle P_{\text{d}}[\gamma]=p_{\text{i}}\text{Tr}\left\{|f\rangle\langle f|\left(\overrightarrow{\prod_{n=1}^{N}}{\cal E}_{{\cal K}_{\gamma}(t_{n})}^{(\eta)}\right)\big{[}|i\rangle\langle i|\big{]}\right\},$$

$$\displaystyle P_{\text{r}}[\gamma^{\text{r}}]=p_{\text{f}}\text{Tr}\left\{|i\rangle\langle i|\left(\overleftarrow{\prod_{n=1}^{N}}{\cal E}_{{\cal K}_{\gamma}(t_{n})}^{(\eta),\text{r}}\right)\big{[}|f\rangle\langle f|\big{]}\right\},$$

$${\mathcal{E}}_1^{(1),\text{r}}[{\rho }_{\gamma }](t_n)$$

$${\mathcal{E}}_2^{(1),\text{r}}[{\rho }_{\gamma }](t_n)$$

$${\mathcal{E}}_0^{(1),\text{r}}[{\rho }_{\gamma }](t_n)$$

$${\mathcal{E}}_0^{(\eta )}={\mathcal{E}}_{00}^{(\eta )}+{\mathcal{E}}_{01}^{(\eta )}+{\mathcal{E}}_{02}^{(\eta )},$$

$${\mathcal{E}}_{01}^{(\eta )}[{\rho }_{\gamma }]$$

$${\mathcal{E}}_{02}^{(\eta )}[{\rho }_{\gamma }]$$

$${\mathcal{E}}_{00}^{(\eta )}[{\rho }_{\gamma }]$$

$$\sum _{\gamma }{P}_{\text{r}}[{\gamma }^r]$$

$${\sigma }_{\eta }[\gamma ]=-\log \sum _{\mathcal{F}[\gamma ]}P({\gamma }_F|\gamma )e^{(Q_{\text{cl}}[{\gamma }_F]-Q_{\text{cl}}^{\eta }[\gamma ])/k_{\text{B}}{T}_{\text{eff}}}.$$

$${\dot{P}}_e=-{\Gamma }_0(2\bar{n}[{\omega }_1(t)]+1)\left(P_e(t)-\frac{\bar{n}[{\omega }_1(t)]}{2\bar{n}[{\omega }_1(t)]+1}\right).$$

$${\dot{P}}_e=-({\Gamma }_{+}+{\Gamma }_{-})\left(P_e(t)-\frac{{\Gamma }_{+}}{{\Gamma }_{+}+{\Gamma }_{-}}\right).$$

$$\Omega (t)$$

$${\Gamma }_{-}(t)$$

$$\delta Q_{\text{q}}$$

$$\delta Q_{\text{q}}(t_n)$$

Full text

I Introduction

The existence of some preferred direction of time is a fundamental concept of physics, captured by the second Law of thermodynamics Eddington (1928); Lebowitz (1993). The irreversibility of physical phenomena is manifested by a strictly positive entropy production, classically identified with the change of entropy of the closed system considered. In the textbook situation of a system ${\cal S}$ exchanging work $W$ with some external operator and heat $Q_{\text{cl}}$ with a thermal reservoir ${\cal R}$ of temperature $T$, entropy production equals the entropy change of the system $\Delta S_{\text{sys}}$ and the reservoir $\Delta S_{\text{res}}=-Q_{\text{cl}}/T$, and reads

[TABLE]

$\Delta F$ is the change of the system’s free energy, satisfying $\Delta S_{\text{sys}}=(W+Q_{\text{cl}}-\Delta F)/T$.

Focusing on the transformations of microscopic systems, stochastic thermodynamics has extended the concepts of thermodynamics to single stochastic realizations or “trajectories” of the system in its phase space Seifert (2005); Sekimoto (2010); Jarzynski (2011). In this framework, heat and work exchanges become stochastic quantities defined for single trajectories $\gamma$, which can give rise to a negative entropy production $\Delta_{\text{i}}s[\gamma]<0$. The fluctuations of $\Delta_{\text{i}}s[\gamma]$ verify the central Fluctuation Theorem (FT) Seifert (2008) $\langle e^{-\Delta_{\text{i}}s[\gamma]/k_{\text{B}}}\rangle_{\gamma}=1$, while the Second Law of thermodynamics remains valid on average. A famous example of such central FT is provided by Jarzynski’s equality (JE) Jarzynski (1997, 2011) characterizing isothermal transformations of systems driven out of equilibrium. Eq.(1) remains valid at the trajectory level, such that entropy produced along a single trajectory $\gamma$ verifies $\Delta_{\text{i}}s[\gamma]=(W[\gamma]-\Delta F)/T$. $W[\gamma]$ is the stochastic work exchange for the trajectory $\gamma$, yielding JE:

[TABLE]

In the classical regime, the stochastic work exchanges $W[\gamma]$ are inferred from the knowledge of the trajectory, which can be recorded without perturbing the system in principle with an arbitrary precision (Fig.1a). Using this protocol, JE has been experimentally probed in various setups, e.g. with Brownian particles Toyabe et al. (2010); Berut et al. (2012) or electronic systems Saira et al. (2012).

The recent developments of quantum and nano-technologies have urged to explore the validity of FTs for out-of-equilibrium quantum systems. In the particular case of Jarzynski’s protocol, the system coupled to a thermal reservoir can now be coherently driven. However in the quantum world, monitoring a system perturbs its evolution (Fig.1b), posing severe difficulties to measure and even to define work exchanges. In a series of pioneering papers Talkner P. and Hanggi (2007, 2009); Esposito et al. (2009); Campisi et al. (2011) it was shown that JE is still valid and can be probed, provided that the system’s (resp. reservoir’s) internal energy change $\Delta U[\gamma]$ (resp. $-Q_{\text{cl}}[\gamma]$) can be known from projective energy measurements performed in the beginning and at the end of the transformation. The work is then inferred from the First Law: $\Delta U[\gamma]=W[\gamma]+Q_{\text{cl}}[\gamma]$. Known as the two-points measurement scheme, this seminal extension of JE in the quantum world has triggered numerous investigations, motivated by two major challenges. The first challenge is practical: How to measure energy changes induced by a microscopic system on a large reservoir? One possible strategy relies on finite size reservoirs whose temperature is sensitive to heat exchanges with the system Pekola et al. (2013); Hekking and Pekola (2013); Campisi et al. (2015), such that both emission and absorption events can be detected using fine calorimetry. But such reservoirs are by essence non Markovian, therefore affecting the dynamics of the system and the resulting work distribution Suomela et al. (2016).

Another strategy involves the convenient platforms provided by superconducting circuits and semi-conducting quantum photonics. In both cases, the radiation produced by quantum emitters (superconducting qubits or quantum dots) is efficiently funneled into well-designed waveguides, and is thus recorded with high efficiency. This recent experimental ability has lead to the development of bright single photon sources Claudon et al. (2010); Somaschi et al. (2016) and to the monitoring of quantum trajectories of superconducting quantum bits Ristè et al. (2013); Murch et al. (2013); Weber et al. (2014). However, standard schemes based on photo-counters do not allow for the recording of photons absorbed by the system, a major drawback for quantum thermodynamics purposes.

The second challenge is of fundamental nature, and consists in identifying quantum signatures in quantum FTs. For instance, the coupling to a thermal bath erases quantum coherences during the system’s evolution. This is an irreversible process, expected to induce some genuinely quantum entropy production Auffèves (2015); Elouard et al. (2017). Furthermore, the incoherent energy exchange between the system and the bath (heat) is carried by quantized excitations. Thus any information depending on the recording of these excitations may be significantly hindered by the eventual inefficiency of the detection scheme. So far, an experimental observation of such genuinely quantum effects has remained elusive.

Here we propose a strategy that addresses both the practical and fundamental issues raised above (Fig.1c and Fig.2). Our proposal is based on an effective two-level system (TLS) coupled to an engineered thermal bath already introduced by some of us Carvalho and Santos (2011); Santos and Carvalho (2011) to perform incoherent quantum computing Santos et al. (2012). A great interest of this scheme is that it allows for the detection of the emission and the absorption of photons in the reservoir, while still keeping its Markovian character. Eventually, both the heat exchanges and the quantum trajectory followed by the system can be fully reconstructed. Note that a quite similar strategy is proposed in Horowitz (2012) to access quantum trajectories of a forced Harmonic oscillator.

We first introduce and validate the platform by checking JE in the simple case where the external drive does not create any coherences in the system, and heat exchanges are perfectly detected. We then turn to the case of a coherence-inducing drive for which we provide a full thermodynamic analysis. Genuinely quantum signatures in thermodynamic quantities are evidenced, as well as how to measure them. Finally, we investigate the influence of the detectors finite efficiency on the expression of the measured FTs. The non-detected heat results in a modified expression of JE that includes an extra term in the entropy production. We show that this term is related to the measurement record and can be computed from the available data.

II Engineered thermal bath

We consider a system whose three levels are denoted by $|m\rangle$, $|e\rangle$ and $|g\rangle$ of respective energy $E_{m}>E_{e}>E_{g}$ (see Fig.2a). This three-level system is coupled to an electro-magnetic reservoir at zero temperature, such that $|m\rangle$ decays to state $|e\rangle$ (resp. $|e\rangle$ decays to $|g\rangle$) with a spontaneous emission rate $\Gamma$ (resp. $\Gamma_{-}$). We assume that $|m\rangle$ is a metastable level, such that $\Gamma\gg\Gamma_{-}$. We denote $\omega_{1}=(E_{e}-E_{g})/\hbar$ and $\omega_{2}=(E_{m}-E_{e})/\hbar$. In addition, the atom is weakly driven by a laser resonant with transition $\omega_{1}+\omega_{2}$ and of Rabi frequency $\Omega$, satisfying $\Omega\ll\Gamma$. Due to its short lifetime, level $|m\rangle$ can be adiabatically eliminated, resulting in an effective incoherent transition rate $\Gamma_{+}=4\Omega^{2}/\Gamma$ from state $|g\rangle$ to state $|e\rangle$ Carvalho and Santos (2011). The states $|e\rangle$ and $|g\rangle$ define our effective qubit of interest (Fig.2b). Without external drive, this qubit relaxes towards some effective thermal equilibrium characterized by the temperature $T_{\text{eff}}$, satisfying $e^{-\hbar\omega_{1}/k_{\text{B}}T_{\text{eff}}}=\frac{\Gamma_{+}}{\Gamma_{-}}$.

In what follows we drive the qubit transition with a Hamiltonian $H_{\text{d}}(t)$, such that $\omega_{1}$ may depend on time: the rates $\Gamma_{\pm}$ can be adjusted accordingly to keep the effective temperature constant. These adjustments are discussed in Ref.Carvalho and Santos (2011) and the Methods Section at the end of this manuscript. The dynamics of the effective qubit density matrix $\rho$ is ruled by the master equation:

[TABLE]

where $H_{0}=(\hbar\omega_{1}/2)\sigma_{z}$ with $\sigma_{z}=|e\rangle\langle e|-|g\rangle\langle g|$. We have introduced the dissipation super-operator $\mathcal{L}[X]\rho=X\rho X^{\dagger}-\tfrac{1}{2}\{X^{\dagger}X,\rho\}$, with $\sigma=|g\rangle\langle e|$. Assuming that a photo-counter ${\cal D}_{1}$ (resp. ${\cal D}_{2}$) detects the photons emitted at frequency $\omega_{1}$ (resp. $\omega_{2}$), one can formulate the evolution of the qubit state conditioned to the measurement records of the detectors in terms of quantum jumps Wiseman and Milburn (2010). Assuming an initial known pure state $|\Psi_{0}\rangle$ of the qubit, and discretizing the time between $t_{\text{i}}$ and $t_{\text{f}}$ such as $t_{n}=t_{\text{i}}+ndt$ (with $n\in\llbracket 0,N\rrbracket$ and $t_{N}=t_{\text{f}}$), the evolution of the system features a stochastic trajectory $\gamma$ of pure states $|\psi_{\gamma}(t_{n})\rangle$. The trajectory is generated by applying a sequence of operators $M_{{\cal K}(t_{n})}$, where ${\cal K}(t_{n})$ stands for the stochastic measurement outcome at time $t_{n}$.

Namely, detecting a photon on detector ${\cal D}_{1}$ (resp. ${\cal D}_{2}$) at time $t_{n}$ corresponds to applying the operator $M_{1}=\sqrt{\Gamma_{-}dt}\sigma$ (resp. $M_{2}=\sqrt{\Gamma_{+}dt}\sigma^{\dagger}$) on the qubit state $|\psi_{\gamma}(t_{n})\rangle$: remarkably in this scheme, the absorption of a photon from the effective heat bath (Fig.2b) is detectable and actually corresponds to the emission of a photon of frequency $\omega_{2}$ (Fig.2a). If no photon is detected, the no-jump operator $M_{0}=\mathds{1}-idtH(t)-\tfrac{1}{2}M_{1}^{\dagger}M_{1}-\tfrac{1}{2}M_{2}^{\dagger}M_{2}$ is applied. The qubit state is then renormalized. Each of these possible evolutions occurs at time $t_{n}$ with probabilities

[TABLE]

III Stochastic thermodynamic quantities

We now define the thermodynamic quantities associated with the system’s quantum trajectory $|\psi_{\gamma}(t_{n})\rangle$. Following Elouard et al. (2017), the internal energy of the qubit at time $t_{n}$ is defined as $U_{\gamma}(t_{n})=\langle\psi_{\gamma}(t_{n})|H(t_{n})|\psi_{\gamma}(t_{n})\rangle$, with $H(t)=H_{0}+H_{\text{d}}(t)$ the total Hamiltonian of the system. The variation of this quantity along trajectory $\gamma$ splits into three terms: there is a work increment $\delta W_{\gamma}(t_{n})$ due to the Hamiltonian time-dependence and given by:

[TABLE]

There is also a classical heat $Q_{\text{cl}}[\gamma]$ corresponding to the energy exchanged with the effective thermal bath under the form of emitted or absorbed excitations of energy $\hbar\omega_{1}(t_{n})$, where $\omega_{1}(t_{n})$ is the qubit effective transition frequency defined as $\omega_{1}(t)=(\langle e|H(t)|e\rangle-\langle g|H(t)|g\rangle)/\hbar$. This increment reads

[TABLE]

Finally, a third contribution must be introduced whenever the system’s dynamics creates coherences in the $\{|e\rangle;|g\rangle\}$ basis. Defined as quantum heat in Elouard et al. (2017), its increment when a jump takes place reads:

[TABLE]

while if no jump takes place, it reads

[TABLE]

The first Law expressed for the trajectory $\gamma$ is then given by

[TABLE]

where $O[\gamma]$ denotes the integrated quantity over the whole trajectory.

IV Jarzynski equality in the quantum regime

Jarzynski’s protocol generally consists in preparing the qubit in the thermal equilibrium state $\rho_{\text{i}}=Z_{\text{i}}^{-1}e^{H(t_{\text{i}})/k_{\text{B}}T_{\text{eff}}}$, then performing a strong measurement of the qubit energy state to project it on an initial pure state $|i\rangle$ of internal energy $U_{\text{i}}$. The qubit is driven out of equilibrium until time $t_{\text{f}}$ when a final strong energy measurement is performed, projecting the qubit onto the final state $|f\rangle$ of internal energy $U_{\text{f}}$. Note that the final density matrix is in general different from the final thermal equilibrium state $Z_{\text{f}}^{-1}e^{H(t_{\text{f}})/k_{\text{B}}T_{\text{eff}}}$. We have introduced the initial and final partition functions $Z_{\text{i}}$ and $Z_{\text{f}}$. The two strong measurements allow inferring the qubit internal energy change $\Delta U[\gamma]=U_{\text{f}}-U_{\text{i}}$, while the bath monitoring allow recording the classical heat exchange $Q_{\text{cl}}$. As stated in the introduction (the demonstration is provided in Section V), the entropy produced in a single trajectory $\gamma$ equals $\Delta_{\text{i}}s[\gamma]=(\Delta U[\gamma]-Q_{\text{cl}}[\gamma]-\Delta F)/T_{\text{eff}}$ Horowitz and Parrondo (2013); Elouard et al. (2017). Inserting Eq.(13), it becomes

[TABLE]

While fully compatible with former expressions for entropy production, Eq.(14) reveals two components respectively involving the amounts of work and quantum heat exchanged. The second component represents a genuinely quantum contribution due to the presence of quantum coherences in the qubit bare energy basis, that quantitatively relates entropy production and energetic fluctuations. In the following, we show that this quantum contribution can be measured in our setup.

We first validate our platform by considering as driving Hamiltonian $H_{\text{d}}(t)=(\hbar\omega_{1}/2)\epsilon\sigma_{z}(t-t_{\text{i}})$ for $t_{\text{i}}\leq t\leq t_{\text{f}}$. The work and heat exchanged in this case take particularly simple forms, the work being due to time-varying energy levels, and the heat to stochastic population changes induced by the bath. More precisely, we have $\delta W_{\gamma}(t_{n})=\langle\sigma_{z}(t_{n})\rangle\hbar\epsilon dt$ and $\delta Q_{\text{cl},\gamma}(t_{n})=(\langle\sigma_{z}(t_{n+1})\rangle-\langle\sigma_{z}(t_{n})\rangle)\hbar\omega_{1}(t_{n})$, while the quantum heat increment and the quantum component of entropy production are zero. Note that this transformation boils down to the famous Landauer’s protocol Landauer (1961); Berut et al. (2012). In Fig.3a and b, we have plotted the quantity $\langle e^{-\Delta_{\text{i}}s[\gamma]/k_{\text{B}}}\rangle_{\gamma}$ and mean entropy production $\langle\Delta_{\text{i}}s[\gamma]\rangle_{\gamma}$ for different values of the coupling constant $\epsilon$ and temperatures. First note that, as expected, JE is verified for the thermodynamic quantities previously defined. Also note that the entropy production vanishes for $\epsilon\ll\Gamma_{-}$. This is also expected since, in this limit, the drive is quasi-static and the qubit is always in equilibrium with the heat bath. This situation corresponds to a reversible transformation and, in our example, amounts to a large exchange of photons between the system and the bath before the energy of the system significantly changes. Finally, entropy production drastically increases when $\epsilon\gtrsim\gamma$. In this case, very few photons are exchanged and the system is driven far from equilibrium.

We now turn to a case with no classical counterpart by considering a drive $H_{\text{d}}(t)=(\hbar g/2)(\sigma e^{i\omega_{1}t}+\sigma^{\dagger}e^{-i\omega_{1}t})$ for $t_{\text{i}}<t<t_{\text{f}}$ and $H_{\text{d}}(t_{\text{i}})=H_{\text{d}}(t_{\text{f}})=0$. This can be implemented by coupling the atom with a classical light field resonant with the transition $\omega_{1}$. It gives rise to Rabi oscillations, i.e. to the reversible exchange of energy (work) between the qubit and the field under the form of the coherent and periodic emission and absorption of photons at the Rabi frequency $g$. Note that this scenario deviates from the previous one as, here, work exchanges induce coherences and population changes in the qubit bare energy basis. According to our definitions, the coupling of this coherently driven qubit to the bath will now give rise to both classical and quantum heat exchanges, as well as a quantum component in the entropy production.

We first check that JE is recovered from the detection events in this quantum regime, as it appears in Fig3c. This shows that our generalized thermodynamic quantities are properly defined. A more thorough analysis can be drawn from Fig.3d, where we plot the average entropy production as a function of the ratio $g/\Gamma_{-}$. There, one can identify three different situations, two of which corresponding to reversible transformations ($\langle\Delta_{\text{i}}s\rangle_{\gamma}\rightarrow 0$). For a very weak drive $g\ll\Gamma_{-}$, the qubit is always in a thermal equilibrium state. This regime corresponds to reversible quasi-static isothermal transformations. In the other extreme, i.e. when $g\gg\Gamma_{-}$, the drive is strong enough to induce almost unperturbed Rabi oscillations. In this case, the transformation is too fast to allow for a stochastic event to take place and heat to be exchanged, i.e. it is adiabatic in the thermodynamic sense and once again reversible.

The intermediate scenario where $g\sim\Gamma_{-}$ gives rise to a maximal average entropy production. Note that, in this situation, the thermodynamic time arrow has a completely different nature from the classical case studied before. Here, entropy is mostly produced by the frequent interruptions of the qubit Rabi oscillations caused by the stochastic exchanges of photons with the bath (a.k.a. heat) (Inset of Fig3d). During such quantum jumps, the quantum coherences induced by the driving process are erased, giving rise to the exchange of quantum heat, and to entropy production of quantum nature. The deep connection between quantum heat and entropy production is captured in Eq.(14) where they are quantitatively related, and in Fig3(e to h), where both histograms concentrate around zero in the reversible regime, and take finite values in the irreversible regime. Such a connection between quantum irreversibility and energy fluctuations has already been evidenced in Elouard et al. (2017) in a different context, where stochasticity is caused by quantum measurement instead of a thermal bath. In both cases, quantum heat exchanges and entropy production appear as two thermodynamic signatures of the same phenomenon, i.e. coherence erasure.

V Finite detection efficiency.

Up to now, we have assumed that the detection of heat is $100\%$ efficient, i.e. that all the photons emitted by the three level atom can be detected. In practice, state-of-the art devices only allow for a fraction $\eta$ of these photons to be indeed detected Natarajan et al. (2012); Gol’tsman et al. (2001). A more realistic scenario is now investigated by considering that the photo-detectors have a finite efficiency $\eta$. As a consequence, between two detected photons, there exist several different possible (or fictitious) trajectories of pure states, which cannot be distinguished by the measurement record: Each of these trajectories corresponds to a particular sequence of undetected emissions and absorptions. This effect induces a decrease of the purity of the qubit state which has to be described by a stochastic density matrix $\rho_{\gamma}$, rather than a wavevector. In the limit where $\eta=0$ (no detector), the evolution of this density matrix is captured by Eq.(3), such that $\rho_{\gamma}(t)=\rho(t)$.

For $0<\eta<1$, the evolution of $\rho_{\gamma}(t)$ is still conditioned on the stochastic measurement outcome of the detector ${\cal K}_{\gamma}(t_{n})$, that now corresponds to applying a set of super-operators $\{{\cal E}^{(\eta)}_{\cal K}[\rho_{\gamma}]\}$ (${\cal K}\in\{0,1,2\}$). Detecting one photon in detector ${\cal D}_{1}$ (resp. ${\cal D}_{2}$) between time $t_{n}$ and $t_{n}+dt$ corresponds to applying the super-operator ${\cal E}^{(\eta)}_{1}[\rho_{\gamma}]=\eta M_{1}\rho_{\gamma}M_{1}^{\dagger}$ (resp. ${\cal E}^{(\eta)}_{2}[\rho_{\gamma}]=\eta M_{2}\rho_{\gamma}M_{2}^{\dagger}$), which occurs with probability $p_{1}^{(\eta)}(t)=\eta\text{Tr}\{M_{1}^{\dagger}M_{1}\rho_{\gamma}(t)\}$ (resp. $p_{2}^{(\eta)}(t)=\eta\text{Tr}\{M_{2}^{\dagger}M_{2}\rho_{\gamma}(t)\}$). When no photon is detected, which occurs with probability $p_{0}^{(\eta)}=1-p_{1}^{(\eta)}-p_{2}^{(\eta)}$, a super-operator decreasing the state purity is applied:

[TABLE]

Note that ${\cal E}_{0}^{(\eta)}$ is a linear interpolation between ${\cal E}_{0}^{(1)}[\rho]=M_{0}\rho_{\gamma}M_{0}^{\dagger}$ applied in the perfect efficiency quantum jump formalism and ${\cal E}_{0}^{(0)}[\rho]$ which is the Lindbladian (right-hand terms in Eq.(3)). After applying the super-operator ${\cal E}^{(\eta)}_{\cal K},$ the density matrix has to be divided by $p_{\cal K}^{(\eta)}$ in order to be renormalized.

Generalized Jarzynski equality.

We now extend the definitions of thermodynamic quantities into the finite efficiency regime. For the sake of simplicity, we restrict the study to a driving Hamiltonian of the form $H_{\text{d}}(t)=(\hbar\omega_{1}/2)\epsilon(t-t_{i})\sigma_{z}$, such that no quantum heat is exchanged. With the definitions proposed above, the increment of measured classical heat reads:

[TABLE]

Despite its apparent similarity with the definition in the perfect efficiency regime, this energy flow does not capture the entire classical heat flow dissipated in the heat bath, as it does not take into account the energy carried by the undetected photons. Therefore, the measured entropy production $\Delta_{\text{i}}s^{\eta}[\gamma]=(\Delta U-Q_{\text{cl}}^{\eta}[\gamma]-\Delta F)/T_{\text{eff}}$ does not check Jarzynski equality as evidenced in Fig.4a) and b). The violation of JE is all the larger as the efficiency $\eta$ is smaller, and JE is recovered in the limit $\eta\to 1$, showing that our definitions for the finite-efficiency case still hold in the perfect efficiency limit. Remarkably, however, it is possible to formulate another equality taking into account the finite efficiency, converging towards the standard JE when $\eta\rightarrow 1$:

[TABLE]

This equality is reminiscent of JE, but it includes a trajectory-dependent correction term $\sigma_{\eta}[\gamma]$. We now show that $\sigma_{\eta}[\gamma]$ can be computed from the measurement record only.

The entropy creation along one trajectory $\gamma$ is usually defined as $\Delta_{\text{i}}s[\gamma]=k_{\text{B}}\log(P_{\text{d}}[\gamma]/P_{\text{r}}[\gamma^{\text{r}}])$, where $P_{\text{d}}[\gamma]$ is the probability of a trajectory $\gamma$ in the direct protocol defined by $H_{\text{d}}(t)$:

[TABLE]

We have introduced $p_{\text{i}}=Z^{-1}_{\text{i}}e^{-U_{\text{i}}/k_{\text{B}}T_{\text{eff}}}$ the probability of drawing the initial pure state $|i\rangle\langle i|$ of internal energy $U_{\text{i}}$ in the initial thermal distribution. The arrow indicates the order in which the sequence of super-operators is applied. $P_{\text{r}}[\gamma^{\text{r}}]$ is the probability of the time-reversed trajectory $\gamma^{\text{r}}$ corresponding to $\gamma$. It is generated by the time-reversed drive $H_{\text{d}}(t_{\text{f}}-t)$ and the reversed sequence of super-operators ${\cal E}^{(\eta),r}_{{\cal K}_{\gamma}(t_{n})}$ Manzano et al. (2015); Crooks (2008); Elouard et al. (2017), such that:

[TABLE]

where $p_{\text{f}}=Z^{-1}_{\text{f}}e^{-U_{\text{f}}/k_{\text{B}}T_{\text{eff}}}$ the probability of drawing the final pure state $|f\rangle\langle f|$ of internal energy $U_{\text{f}}$ in the final thermal distribution.

When $\eta=1$, the time-reversed operators at time $t_{n}$ reduce to

[TABLE]

By inserting these expressions in Eq.(21), we express the ratio $P_{\text{r}}[\gamma^{\text{r}}]/P_{\text{d}}[\gamma]=e^{-(\Delta U[\gamma]-\Delta F)/k_{\text{B}}T_{\text{eff}}}e^{Q_{\text{cl}}[\gamma]/k_{\text{B}}T_{\text{eff}}}$, which is JE.

When $\eta<1$, the ratio $P_{\text{r}}[\gamma^{\text{r}}]/P_{\text{d}}[\gamma]$ is modified. It is useful to decompose ${\cal E}^{(\eta)}_{0}$ into three super-operators conserving purity:

[TABLE]

with

[TABLE]

From this decomposition, we generate a set ${\cal F}[\gamma]$ of fictitious trajectories $\gamma_{F}$ of pure states, compatible with $\gamma$: a trajectory $\gamma_{F}\in{\cal F}[\gamma]$ is built by applying a sequence of super-operators $\{{\cal E}_{{\cal K}_{\gamma_{F}}(t_{n})}\}_{n}$, with ${\cal K}_{\gamma_{F}}(t_{n})={\cal K}_{\gamma}(t_{n})$ if ${\cal K}_{\gamma}(t_{n})$ belongs to $\{1,2\}$ (i.e. when a photon has been detected), and ${\cal K}_{\gamma_{F}}(t_{n})\in\{00,01,02\}$ when ${\cal K}_{\gamma}(t_{n})$ is zero (no photon detected). ${\cal F}[\gamma]$ thus contains $3^{{\cal N}_{0}[\gamma]}$ fictitious trajectories, where ${\cal N}_{0}[\gamma]$ is the number of time-steps during which no photon is detected.

The perfect-efficiency fictitious trajectories are interesting because they can be time-reversed using the rules Eq.(22)-(24), while their probability distribution checks $P_{\text{d}}[\gamma]=\sum_{\gamma_{F}\in\mathcal{F}[\gamma]}P_{\text{d}}[\gamma_{F}]$. Eventually, the following equality is derived:

[TABLE]

We have introduced the conditional probability $P_{\text{d}}[\gamma_{F}|\gamma]=P_{\text{d}}[\gamma_{F}]/P_{\text{d}}[\gamma]$ of the fictitious trajectory $\gamma_{F}\in{\cal F}[\gamma]$, given the detected trajectory $\gamma$. In order to find Eq.(19), we now define:

[TABLE]

$\sigma_{\eta}[\gamma]$ can be numerically computed, based solely on the measurement outcomes. The computation involves the simulation of a sample of the trajectories $\gamma_{F}\in{\cal F}[\gamma]$ for each trajectory $\gamma$. The corrected expression $\langle e^{-\Delta_{\text{i}}s[\gamma]/k_{\text{B}}-\sigma_{\eta}[\gamma]}\rangle_{\gamma}$ is plotted in Fig.4a): JE is verified, showing that the experimental demonstration of a FT in a realistic setup is within reach.

Our results explicit a new quantum-classical border in the thermodynamic framework, which is illustrated in Fig.4(b to f) where the deviation from JE is studied as a function of the reservoir temperature. In the quasi-static limit characterized by $\epsilon\ll\Gamma_{-}$ (high temperature case), many photons are exchanged and the system is allowed to thermalize before any significant work is done: the transformation is reversible. We see from Fig.4c) and d) that in this limit, the effect of detection inefficiency in the entropy production is merely to widen the spread of the distribution but conserving its Gaussian shape around zero. The statistical properties of the distribution are not affected, even if many photons are missed, and JE is verified (Fig.4b). This corresponds to a semi-classical situation where, even though the exchange of heat is quantized, the overall thermodynamic behavior of the system is equivalent to a classical ensemble. In this limit, there is no significant information loss due to the undetected photons.

On the other hand, as soon as work and heat are exchanged at similar rates $\epsilon\sim\Gamma_{-}$, the system is brought out of equilibrium and the transformation is irreversible. In this limit, few photons are exchanged before a non negligible amount of work is done on the system. Each missed photon represents, then, a significant information loss, such that the measured distribution of entropy production is severely affected (see Fig.4e) and f)). Such information loss is quantifiable by the violation of the non-corrected JE, and is a direct consequence of the quantization of heat exchanges.

VI Conclusion

We presented a realistic setup based on a driven three-level atom allowing to simulate a qubit in equilibrium with an engineered reservoir, giving full access to the distribution of the heat dissipated by the qubit. We exploited this scheme to characterize irreversibility in the case of a coherently driven qubit, coupled to a thermal bath: we evidenced genuinely quantum contributions to the entropy production, and showed these contributions do have some energetic counterpart. In a second part, we took into account finite detection efficiency and evidenced deviations from Jarzynski equality that can be tested and corrected. We derived a modified equality, involving only quantities computed from the finite-efficiency measurement record. This work opens avenues for the experimental verification of fluctuation theorems in quantum open systems. Owing to the versatility of the scheme, various reservoirs could be simulated, including non-thermal ones. Furthermore, quantitative relations between entropy production and energy fluctuations are the basis for energetic bounds for classical computation Landauer (1961); Bennett (1982); Faist et al. (2015). In this work we study situations where equivalent relations can be derived in the quantum regime, providing new tools to explore energetic bounds for quantum information processing.

Acknowledgments

MFS acknowledges CNPq (Project 305384/2015-5). NKB thanks the support from the Brazilian agency CAPES. This work was supported by the Fondation Nanosciences of Grenoble under the Chair of Excellence ”EPOCA”, by the ANR under the grant 13-JCJC-INCAL and by the COST network MP1209 ”Thermodynamics in the quantum regime”.

Appendix

VI.1 Simulating a heat bath with an engineered environment

The rates $\Gamma_{\pm}$ induce the same dynamics as a heat bath, if and only if

[TABLE]

We have introduced the population of the excited qubit’s state $P_{e}=(\text{Tr}[\rho(t)|e\rangle\langle e|]$, the mean number of thermal photons in the effective bath $\bar{n}[\omega]=(e^{\hbar\omega/k_{\text{B}}T_{\text{eff}}}-1)^{-1}$ and $\Gamma_{0}$ the spontaneous emission rate of the qubit. The dynamics induced by the engineered environment is:

[TABLE]

By identification, we find the conditions:

[TABLE]

Eq.(29) can be fulfilled by tuning accordingly the intensity of the laser drive. Eq.(30) requires to tune the incoherent rate of transition from $|e\rangle$ to $|g\rangle$. This task can be performed by embedding the three level atom in a quasi-resonant optical cavity and tuning the frequency between state $|e\rangle$ and $|g\rangle$ e.g. with a non-resonant laser (Stark effect) according to a protocol $\omega_{1}^{\prime}(t)$ designed to fulfill Eq.(30). Note that the implemented protocol $\omega_{1}^{\prime}(t)$ is in general different from the target protocol $\omega_{1}(t)$ which defines the constraints Eq.(29)-(30) on the qubit dynamics.

VI.2 Expression of the quantum heat

For a general driving Hamiltonian $H_{\text{d}}(t)=\delta(t)\sigma_{z}+\mu(t)\sigma+\mu^{*}(t)\sigma^{\dagger}$, the quantum heat $\delta Q_{\text{q}}(t_{n})$ exchanged when a jump takes place reads (for $\eta=1$):

[TABLE]

Between two photon detections, the quantum heat increment reads:

[TABLE]

Note that $\delta Q_{\text{q}}(t_{n})$ is zero when the qubit state has no coherence in $|e\rangle-|g\rangle$ basis, such that $\langle\sigma\rangle=\langle\sigma^{\dagger}\rangle=0$, $(P_{e}(t_{n}),P_{g}(t_{n}))=(1,0)$ (resp. $(0,1)$) when ${\cal K}_{\gamma}(t_{n})=1$ (resp. ${\cal K}_{\gamma}(t_{n})=2$): the qubit has zero probability to absorb a photon at time $t$ when $P_{e}(t_{n})=1$, and zero probability to emit a photon when $P_{g}(t_{n})=1$.