Using a quantum work meter to test non-equilibrium fluctuation theorems.

Work is an essential concept in classical thermodynamics, and in the quantum regime, where the notion of a trajectory is not available, its definition is not trivial. For driven (but otherwise isolated) quantum systems, work can be defined as a random variable, associated with the change in the internal energy. The probability for the different values of work captures essential information describing the behaviour of the system, both in and out of thermal equilibrium. In fact, the work probability distribution is at the core of "fluctuation theorems" in quantum thermodynamics. Here we present the design and implementation of a quantum work meter operating on an ensemble of cold atoms, which are controlled by an atom chip. Our device not only directly measures work but also directly samples its probability distribution. We demonstrate the operation of this new tool and use it to verify the validity of the quantum Jarzynksi identity.

Introduction

Classical fluctuation theorems establish surprising relations between non-equilibrium and equilibrium concepts. In particular, the work performed on a system during non-equilibrium processes is connected with key concepts of equilibrium thermodynamics, such as the free-energy[1, 2]. These relations have been verified in various experiments involving microscopic thermodynamic systems[3-5]. Recent advances in quantum technologies enable the control of small quantum systems that can be manipulated far from the regime where the usual thermodynamical laws are obeyed. This triggered the development of the rapidly growing field of non-equilibrium quantum thermodynamics[6-9].

When quantum fluctuations dominate, defining and measuring work and heat, two central concepts in classical thermodynamics, is non-trivial. For driven, but otherwise isolated, quantum systems, work w is a random variable associated with the change in the internal energy[10], as the first law of thermodynamics indicates. Thus, the commonly accepted definition of quantum work requires a two-time measurement strategy, which consists of performing two projective energy measurements, one at the beginning and the other at the end of the process. Then, work is associated with the measured energy difference. However, implementing the two-time measurement is experimentally difficult[11, 12] due to the fact that the two projective measurements are unavoidably disruptive (see ref. [13] for an ion trap implementation). Alternative methods to evaluate the work probability distribution that rely on the direct estimation of its Fourier transform were also proposed in refs. [14, 15] and later implemented in NMR experiments[16].

In this paper we present the design and the experimental implementation of a “quantum work meter” (QWM) operating on an ensemble of cold atoms, combining the idea presented in ref. [17] and the experimental setup used in ref. [18]. Our QWM is conceptually different from previous work-measurement devices. Its main advantage is that the QWM efficiently samples P(w), which is a direct observable in the experiment. Namely, our QWM not only directly measures work but also directly samples its probability distribution P(w) (i.e. the outcome w is obtained with probability P(w)). As the work probability distribution plays a central role in the fluctuation theorems of non-equilibrium quantum thermodynamics, the QWM is an ideal tool to test their validity. In particular, we use it to verify the Jarzynski identity[1, 10, 19–21].

Results

Work measurement and the QWM

A QWM is an apparatus that measures the work performed on a driven quantum system whose Hamiltonian varies from an initial H to a final with eigenvalues E and , respectively. For an isolated system S, with a D-dimensional space of states, the number of different values of work is bounded by D 2. Therefore, the pointer of the QWM has D 2 distinct positions (one for each value of ). The QWM presented here enables us to choose H and (fixing the possible values of w) and to vary the intermediate driving (inducing different evolution operators denoted as ). In this way, we vary the probability P(w), which depends on the intermediate driving .

By sampling P(w), we use the QWM to verify a fundamental result in non-equilibrium quantum thermodynamics: the Jarzynski identity. This identity states that for any initial state with populations identical to the ones associated to a thermal Gibbs state and for any distribution P(w), the linear combination , where β = 1/k T is the inverse temperature of the system, is an equilibrium property (rather than a non-equilibrium one). The Jarzynski identity (see Supplementary Note 1) readswhere ΔF is the free energy difference between the thermal states associated with the Hamiltonians H and . In the absence of degeneracies, this implies that the vector formed by the D 2−1 measured probabilities belongs to a D 2−2 dimensional hyperplane: the ‘Jarzynski manifold’ (as shown in Supplementary Note 1, further constraints restrict this dimensionality to (D−1)2). With the QWM we measure P(w) for different driving fields showing that all probability vectors belong to the same manifold. By characterising this manifold, we not only verify the identity but also independently estimate the free energy difference ΔF [1,3-5].

The work distribution sampled by the QWM[10, 19–21] is:

Thus, P(w) is the probability density of finding the energy difference w after a measurement of H followed by an intermediate driving and a final measurement of . This is indeed the case if p is the probability of obtaining E when measuring H and p is the probability of obtaining when measuring given that E was detected at the beginning. Equation (2) defines a probability density that is independent of the initial coherences in the energy basis. For the discrete D 2 values of w we will use P(w) to denote the probability (not the density) of each w. The concept on which our QWM is based was first discussed in refs. [17, 22] where it was noticed that the work done on S, can be detected by performing a generalised quantum measurement, which enables the number of outcomes to be larger than D. This can be done by entangling with an ancilla that stores a coherent record of w. Then a standard measurement on S can reveal w. Similar strategies have been later studied and extended to other contexts in refs. [22-24].

Design and operation of the QWM

A pictorial representation of the protocol we follow to operate the QWM is shown in Fig. 1a. The QWM is designed to measure the work done on a system whose Hamiltonian changes from H to and which is subjected to a driving in between. We couple to a continuous variable system and use to denote its position (the generator of translations along the momentum p). A coherent record of w is created by an ‘entangling interaction’ between and that must take place before and after the driving . The unitary operators representing these interactions are: and , where λ is a coupling parameter. Thus, U and respectively translate along p by a displacement proportional to (−λH) and . Then, as shown in detail in Supplementary Note 3, the final measurement of p on yields a random result whose distribution is a smeared version of the true work distribution P(w) defined in Eq. (2). In fact, outcome p is obtained with a probability density , where the window function is fixed by , the initial state of (thus, by localising we improve the accuracy in the estimation of P(w)).

Fig. 1

The quantum work meter. a A quantum circuit for the quantum work meter (QWM). and are entangled so that the eigenvalue of the observable of the system is coherently recorded by . Then is driven by . Finally, another entangling operation between and creates a record of w on . In the experiment, is encoded in the motional degree of freedom of the atoms along the vertical direction z, which also evolves while freely falling. is the pseudospin associated with the Zeeman sub-levels of a 87Rb atom. b Physical operations for the QWM on an atom chip: (i) The atoms, prepared in state , are released from the trap, and a RF field generates an initial pseudo-thermal state. (ii) After 2.4 ms, internal and motional degrees of freedom are entangled with a magnetic gradient pulse (U), applied for a duration of τ = 40 μs. (iii) Another RF field drives . (iv) 3.1 ms after the application of U, a second magnetic gradient pulse () is applied for a duration of μs. At this stage, keeps a record of the different work values. (v) After 18.2 ms from the application of , the positions and optical densities of the atomic clouds are measured. The number of atoms in each cloud reveals the work probability in a single experimental realisation. c Image of the four clouds obtained at the end of a single run of the QWM. The four possible values of w fix the position of each cloud

A ‘universal’ QWM is an apparatus which can measure w and sample P(w) for any possible choice of H and . To build it, we need enough control to enforce the entangling operators U and for any choice of H and . Remarkably, this is achieved for a 2-level system by the atom chip implementation we describe below.

Experimental implementation of the QWM

To describe our QWM we should identify the physical systems representing and , the way in which H and can be chosen, and how the associated U and are implemented. In our experiment we represent by the subspace associated with the Zeeman sublevels and of a 87Rb atom that, as in ref. [18] behaves as a two-level system (see below). The motional degree of freedom of the atom plays the role of .

A key element of the QWM presented here is the atom chip[25], which efficiently entangles the internal and motional degrees of freedom of an atom just ~100 μm away from the atom chip surface, through short and strong Stern-Gerlach type magnetic gradient pulses. These pulses are generated using a 3-current-carrying wire setup on the atom chip surface (described in ref. [26] and the Methods section). A gradient pulse along the z direction with amplitude and duration τ, induces a momentum kick m δp on an atom in the m state (, where μ B and g are, respectively, the Bohr magneton and the Landé factor[18]). The evolution of the state of the atom induced by such a pulse is described by the unitary operator , where . This physical operation translates along the momentum p by a displacement (notice that the operator defines the magnetic dipole moment of the atom since ). As described below, we apply two gradient pulses with different amplitudes ( and ) and different durations (τ and ). Thus, defining and , U and implement the required entangling operation U and , respectively. In this implementation λ is consequently replaced by −δp/E and , enforcing . The momentum kicks induced by both pulses are controlled in the experiment, and consequently, by fixing their ratio, we can simulate an arbitrary system with initial and final Hamiltonians H and which are characterised by having the same ratio. Finally, let us note that the two pulses utilise and with opposite signs to ensure that the sequence creates a record of work corresponding to .

To achieve universality we only need to be able to fix the energy splitting E and of H and , as well as their eigenbasis. The traces of H and (the sum of their eigenvalues) do not affect P(w) but only add a constant to all values of w. As arbitrary E and can be simulated and any change of basis can be absorbed into , we conclude that our atom chip QWM can sample P(w) for an arbitrary 2-level system and is thus universal.

The 87Rb atoms are magnetically trapped in state and evaporatively cooled to a Bose–Einstein condensation (BEC). The BEC is released from the trap and a radio-frequency (RF) pulse is used to prepare a superposition of and . A strong homogeneous magnetic field (created by external coils) suppresses the transitions taking into the state (due to the non-linear Zeeman effect[18]). The initial populations (p 1 and p 2) are chosen so that the temperature is determined via βE = ln(p 1/p 2). The initial motional state is a wave-packet localised in position and momentum.

It should be noted that the initial internal state of the atom, while having the same populations as defined by the temperature of a thermal state, is still a pure state. However, the quantum coherences of this initial state do not affect the results of the QWM. As explained in Supplementary Note 3, the contribution of the initial coherences to the final probability is multiplied by the overlap between the motional states of the atom associated with the different values of work. Thus, when the atomic clouds associated with the different work values are well separated, the effect of initial coherences is negligible. In this regime, our experiment gives the same result as the one we would obtain by preparing an initial thermal state (with no coherences). The study of the importance of the initial coherences in the definition of work is an interesting topic in itself, which is beyond the scope of our paper (see, for example, refs. [24, 27–30]).

The experimental sequence, presented in Fig. 1b, is: (i) prepare the initial state and release the cloud (which then freely falls along z, the direction of gravity), (ii) apply the magnetic gradient U along z, (iii) apply the driving by exposing the atoms to a RF field resonant with the Zeeman splitting induced by the homogeneous bias field, (iv) apply the gradient , (v) obtain an image of the four clouds after a time-of-flight and count the number of atoms in each cloud. More details of the experiment can be found in the Methods section and Supplementary Note 3. For the experimental demonstration presented here we set the ratio between the measured momentum kicks induced by the two pulses to . Hence, our realisation of the QWM samples the work distribution of a simulated system in which the energy splitting is reduced to 56% of its original value, from E to , while driven by .

Figure 1c shows a typical image obtained by the QWM. Four clouds are visible. From the positions of the centre of each cloud, , we infer the total momentum shift, , induced by the pulses on that cloud (we take into account both the free fall and the kicks induced by the pulses, see Supplementary Note 3). Then, we obtain the corresponding value of work as (w is proportional to E, whose value, together with the measurement of the remaining quantities, determine the work w). Furthermore, the probability P(w) for each w is directly measured by the number of atoms in each cloud. Notably, this experiment determines the entire P(w) distribution in a single shot.

Testing the Jarzynski identity

We repeat the experiment fixing the timing, duration and pulse strength. We consider three initial β’s and vary the intermediate driving by changing the duration of the RF field. In this way, we obtain many sets of probability distributions, each of which defines a 3D-vector (as there are three independent probabilities). When we represent all these vectors in the same 3D-plot, we find that they all belong to the same β-dependent manifold. Figure 2a shows that this manifold is a β-dependent line (the dimensionality of this ‘Jarzynski manifold’ is (D−1)2, which in this case equals 1).

Fig. 2

The Jarzynski identity. a Each point defines a probability vector (with its experimental error) measured for a certain driving. The three lines correspond to three temperatures: (blue circle), 1.11 ± 0.02 (red square) and 1.75 ± 0.04 (grey triangle). For each temperature all points lie in the same Jarzynski manifold (which in this case is a line). Reported errors are the SEM of three independent experiments with the same initial parameters and driving. The projections onto the three different axes of the probabilities are shown in detail in Supplementary Fig. 2. b becomes independent of the duration of the intermediate driving (for three temperatures). The dots are the calculated values using the measured work distribution in the Jarzynski identity, and the solid line is the theoretical estimate of βΔF (with an uncertainty due to the uncertainties in the temperature and energy splitting). Error bars are the SEM

Using the measured work probabilities we calculate the exponential average of the work for each driving field. Figure 2b displays the value of as a function of the duration of the intermediate RF field, that parametrises . As established by the Jarzynski identity, G is independent of the driving field and only depends on β. The horizontal lines in Fig. 2b are the theoretically predicted values of βΔF, obtained from a direct calculation (with its own theoretical uncertainty, due to the error in the estimation of βE). This calculation simply involves computing the initial and final partition functions, respectively, denoted as Z and , and using the identity . We find, as the Jarzynski identity establishes, G = βΔF. From Fig. 2b one can notice that the largest errors in the estimation of βΔF appear for βE = 1.75. In this case, for two values of w and, therefore, the relative error in the atom number estimation is large, inducing a larger error in the estimation of βΔF.

In Table 1 we compare the measured and estimated values of βΔF. The uncertainty in the estimation of βΔF and ΔF is close to 10%, which is enough to distinguish the three values of βΔF. On the other hand, in the case of ΔF, there is a significant overlap in the measured values which does not allow to properly distinguish between the three different cases due to the error in the estimation of βE.

Table 1

Estimates of βΔF and ΔF for three different temperatures

βE	βΔF (JI)	βΔF (PF)	ΔF/E (JI)	ΔF/E (PF)
0.58 ± 0.02	−0.36 ± 0.04	−0.35 ± 0.03	−0.62 ± 0.07	−0.60 ± 0.06
1.11 ± 0.02	−0.63 ± 0.05	−0.63 ± 0.04	−0.57 ± 0.05	−0.57 ± 0.04
1.75 ± 0.04	−0.92 ± 0.09	−0.93 ± 0.06	−0.53 ± 0.05	−0.53 ± 0.04

We show the estimation obtained using the Jarzynski identity (JI) and from a direct calculation of the partition function (PF)

Discussion

We presented and implemented a QWM, a new device directly sampling the work distribution on an ensemble of cold atoms. Our QWM can be used to simulate the behaviour of an arbitrary 2-level system. We implemented it with an atom chip and verified the Jarzynski identity over a wide range of non-equilibrium processes. This is the first experiment, and so far the only one, directly sampling P(w) offering advantages and different perspectives over previous work-measurement schemes. Remarkably, in this cold atom experiment, the QWM extracts full statistical information about the work distribution in a single shot.

Methods

Initial state preparation

After preparing the BEC, a homogeneous magnetic field of 36.7G (25 hMHz/μ B, where h is Planck’s constant) is used to push the transition to out of resonance by ~180 kHz due to the non-linear Zeeman effect, which is larger than the power broadened driving RF field of . This ensures that the atoms behave as 2-level systems. The BEC is released from the trap and a RF pulse is used to prepare a superposition of and . By varying the relative populations we consider three different pseudo-thermal states. The initial motional state is a wave packet , well localised at z 0 = 91 ± 1.2 μm from the chip with momentum ~0.

Entangling operations and measurement

An inhomogeneous magnetic field is used to couple spin and motional degrees of freedom. This is generated by a current I = 0.85A in the 3-wire setup during a time τ. The three parallel gold wires lie on the x direction of the chip surface (Fig. 1b). They are 10 mm long, 40 μm wide and 2 μm thick. Their centres are at y = −100, 0, 100 μm and the same current runs through them in alternating directions (−I, I, −I, respectively), creating a 2D quadrupole field at z = 98 μm below the chip. After a time of flight of 2.4 ms the atoms are at z ~ 119 μm. At this point the first gradient pulse implements U: τ = 40 μs with an amplitude of G/mm, such that the momentum kick is along +z. Then, after 3.1 ms the atoms are at mm and the second gradient pulse implements : μs, G/mm, such that the momentum kick is along −z (since the gradient direction is inverted for z > 200 μm). The relative strengths of the spin-dependent forces sets the energy splitting of the Hamiltonians which in this case is on average (this is the measured value, where the error takes into account fluctuations in the initial position of the cloud and in the gradient pulses). In between the entangling operation is applied with a RF pulse. Finally, an image of the atomic clouds is obtained after a time-of-flight of 18.2 ms after the second gradient (the clouds are centred around z ~ 3 mm). The position and number of atoms of each cloud are determined. The momentum shifts of each cloud (that codifies the value of w) are obtained from the difference in positions between the clouds, that follow approximately classical trajectories (Supplementary Note 3).

Uncertainties

The main source of position error is the initial distance of the cloud from the atom chip, whose uncertainty is  ~ 1%. This error is later translated to momentum uncertainty, since the field gradients are position dependent. The field gradients have a fractional uncertainty of 10−3 due to current fluctuations[18]. The central position of each cloud is estimated by fitting a Gaussian profile. Each work probability is estimated as a normalised sum of the measured optical density in a relevant region around the cloud, introducing a probability uncertainty (due to atom numbers uncertainty). Our ~5 μm optical resolution also induces an error in the determination of the position for each cloud. We perform three different runs for each combination of initial state population ratios and intermediate driving and use the average values of position and probability. This gives us a position uncertainty of ~0.015 mm and a probability uncertainty of ~0.015 (standard error).

Data availability

The data that supports the findings of this study are available from the corresponding author upon request.