Calorimetric measurement of work in a quantum system

We propose a calorimetric measurement of work in a quantum system. As a physical realization, we consider a superconducting two-level system, a Cooper-pair box, driven by a gate voltage past an avoided level crossing at charge degeneracy. We demonstrate that, with realistic experimental parameters, the temperature measurement of a resistor (environment) can detect single microwave photons emitted or absorbed by the two-level system. This method would thus be a way to measure the full distribution of work in repeated measurements, and to assess the quantum fluctuation relations.

demonstrate that a calorimetric measurement gives both a theoretical and experimental tool to test the Jarzynski equality (JE) and other fluctuation relations in an open quantum system, and to analyze the distribution of dissipation in them, based on the very principle of conservation of energy. We focus on an experimentally feasible two-level system, a superconducting Cooperpair box (CPB) [19][20][21] subject to Landau-Zener (LZ) interband transitions [22][23][24]. Because of the small heat capacity and weak relaxation to the phonon bath, the calorimetric measurement on the electron gas (resistor) turns out to be a very feasible experimental method.
Consider the basic setting for non-equilibrium fluctuation relations [14,15]: the system is coupled to a bath and is initially in thermal equilibrium with it. Thereafter it is driven by a control field in a non-equilibrium fashion. We are interested in the distribution of work performed on the system in repeated experiments with identical driving protocol among the realizations. In order to secure thermal equilibrium in the next realization, the interval between the repetitions is ideally infinitely long. We mainly consider the protocols where the drive signal returns to its initial position or to an equivalent one with equal free-energy at the end. Our main statement in this respect is then that, due to the energy conservation, the work W done on the system under these conditions is equal to the heat Q dissipated to the environment in the time interval from the beginning of the driving till the end of the equilibration period. With this premise we can avoid the formidable task of analyzing the work itself during the quantum evolution, characterized by difficult-to-analyze time-orderings of the operators involved. The calorimetric principle can, at least in principle, be extended to more general protocols by bringing the system back to the initial position adiabatically at the end of the protocol.
Before going into the detailed analysis of the calorimetric method, let us first define the quantum system and present results in certain limiting cases for the work distribution, which can be obtained by known methods. We focus on a quantum two-level system, shown in figure 1 and characterized by an avoided level crossing under the operation of the control parameter q. For simplicity, we consider a symmetric set-up, where the maximum energy spacing between the ground (g) and the excited (e) level is E max in the beginning and the end of the drive, and E min at the minimum in the middle. We normalize q such that it obtains the values ∓1/2 at the beginning and the end of the drive, respectively, and assumes a value q = 0 at the avoided crossing. Suppose for the basic illustration that the system evolves in a unitary way from the initial state: this represents very weak coupling to the environment. In other words, the system is initially thermalized by coupling it to a bath with inverse temperature β, and decoupled from the bath for the duration of the driven evolution. This corresponds to a common situation for considering the fluctuation relations [10,11,14]. The value of work in each realization is then determined by the internal energy stored in the system during the operation. There are two possibilities: the system makes an LZ transition around q = 0 between the two states with probability P LZ , or it does not make the transition with probability 1 − P LZ . If the system makes the transition, work ±E max is done on the system (+ for g → e transition, and − for the e → g transition) by the end of the driving, otherwise the performed work vanishes. The distribution of work W can then be written as in agreement with the results in [2,9]. Here δ(E) refers to the Dirac delta function, and ρ 0 gg = (1 + e −β E max ) −1 is the initial thermal population in the ground state. We may also state that the same distribution as for work applies for heat in the sense discussed above, since after A two-level system with avoided crossing between the ground (g) and the excited (e) state. In (a) the control parameter q is swept from left to right and the system is initially in the ground state and evolves either in the ground state (probability 1 − P LZ ) or makes an LZ transition (probability P LZ ). (b) Examples of the distribution of work in repeated sweeps under the conditions that the system starts in thermal equilibrium at various temperatures, and it is not coupled to the environment during the evolution. The distributions correspond to (1) with artificial broadening of the peaks, and for P LZ = 0.3, and β E max = 0.1, 0.3 and 1.0 (red, blue and black curves, respectively). The peak at W = −E max arises from e to g LZ transitions, the peak at W = 0 is for adiabatic evolution, and that for W = +E max arises from g to e transitions, as indicated by the diagrams within the (b) panel.
an infinitely long time, for even weak coupling, the system has relaxed back to the thermal equilibrium state by exchanging energy with the environment. Examples of work distributions (with broadening) are shown in figure 1(b) at a few bath temperatures. As for quantum systems evolving unitarily from a thermal state, the JE holds naturally for the two-level system. Indeed we see immediately that for the distribution of (1). Here the brackets · refer to averaging over the experimental realizations. Equations (1) and (2) simply confirm the earlier understanding of JE in a closed quantum system. The true challenge, which can be addressed by the here-proposed calorimetric detection, is, however, the measurement of work in an open quantum system coupled to the environment also during the driving period.  Although our proposal and arguments are quite general, as will be discussed below, we next specify a particular physical system, a CPB, where concrete results about the distribution of work, and a description and analysis of the calorimetric measurement are immediately possible and can be implemented experimentally. The CPB [19][20][21] consists of a superconducting island connected to a superconducting lead by a Josephson tunnel junction. The system is described by the circuit scheme in the inset of figure 2 and it is characterized by a voltage source V g , coupling gate capacitance C g , a Josephson junction with energy E J and capacitance C J . We denote C ≡ C g + C J . Resistor R forms the dissipative environment of the box. In the regime ≡ E J /(2E C ) 1, where E C = 2 e 2 /C is the charging energy of the box, we can treat the CPB as a two-level quantum system [21]. Denoting with |0 and |1 the states with zero and one excess Cooper-pairs on the island, respectively, the Hamiltonian of the system readsĤ where q = C g V g /(2e) is the normalized gate voltage. We assume a linear gate ramp q(t) = −1/2 + t/τ over a period τ = 1/q starting at t = 0, when the state |0 is the ground state of the CPB. The energy gap separating the ground and the excited state of the system is given by E = 2E C q 2 + 2 . Thus the system passes a minimum energy gap E min = E J at q = 0, see figure 1(a), where it makes a transition with the probability P LZ = exp(− π E 2 J τ 4h E C ) according to the standard LZ model [22][23][24]. The value of the energy gap at the ends of the trajectory is We proceed further by applying the calorimetric principle to an open driven two-level system, the main example being a CPB. Full description of the statistics of the work done and the heat generated by the CPB under the driven evolution with LZ transitions in the dissipative regime is beyond the focus of the present work. We restrict the discussion to the regime when energies on the order of the minimum gap E J or smaller can be neglected. In this approximation, dynamics of the two-level system can be described by the master equation for probabilities p 0 and p 1 for the system to be in the charge states |0 and |1 : with the rates of transitions between these states given by where E = 2E C q is the energy dissipated to or absorbed from the environment (resistor) in the is the corresponding probability density of the system to emit energy E to the environment [25].
The standard correlation function in (6) is given by Here, R Q ≡h/e 2 and the relevant impedance of the CPB reads where R eff = λ 2 R, C eff = C J C /C g and λ ≡ C g /C is the coupling parameter. This masterequation-based approach to the description of the system dynamics can be justified most directly in the 'incoherent' regime, when the broadening of the energy levels due to the coupling to environment is larger than E J . However, as discussed below, it also quantitatively reproduces main results for the coherent LZ dynamics.
As an illustration of this, we consider first the case of vanishing temperature T . In this case, the system dynamics reduces to the transitions out of the state |0 only, and equation (4) can be solved for the probability p 0 as follows: and p 0 = 1 for t < τ/2. Using the normalization condition satisfied by P(E) (6), dEP(E) = 1, and taking into account that at vanishing temperature, P(E) = 0 for E < 0, one can immediately see from equation (9) that in the situation of indefinitely increasing energy difference E between the states |0 and |1 , the LZ probability, i.e. the probability p 0 (t)| t→∞ to stay in the state |0 at large times, is independent of the coupling to environment: Such independence of the zero-temperature LZ probability of the dissipation strength is one of the main known properties of the dissipative LZ dynamics [26][27][28][29], and is reproduced here from very elementary considerations.
In the evolution protocol relevant for the CPB, when the energy difference is limited by the maximum E C , the final probability P 0 = p 0 (t = τ ) to be in the state |0 at the end of each ramp follows from equation (9) as Equations (9) and (10) make it possible to find immediately the zero-temperature distributions of work W done on the CPB by the changing gate voltage and the heat Q released into the resistor in each cycle. Since at T = 0 the evolution of the CPB starts from the definite charge state |0 and, after complete relaxation, comes back to the same state, the work W and the heat Q coincide not only on average but for each cycle. This means that the distributions of heat and work coincide, and the work distribution p(W ) is given by the probability that the transition from state |0 to state |1 happens when the energy difference between these states is W : Some of the qualitative features of this expression, e.g. existence of the δ-peak at W = E C with the intensity dependent on the dissipation strength, and a broader peak in the distribution at W = 0 associated with the known behavior of the energy-loss distribution P(E) at small energies, P(E)∝E α , with α = −1 + 4R eff /π R Q (see e.g. [25]), agree with the numerical results presented in figure 2 even with non-vanishing temperatures. The δ-peak at W = 0 that can be seen in this figure is the trace of the δ-peak at W = 0 in the dissipation-free distribution (1). Another case where our master-equation approach to the LZ transitions allows a simple evaluation of the distribution of work is the limit of weak dissipation, R eff R Q , and relatively large temperatures T E J . For small R eff (and not too small coupling parameter λ) the cut-off frequency ω C = 1/R eff C eff of the impedance (8) is large, ω C E C /h. In this case, the tunneling rates ± in the relevant energy range up to E C can be calculated from equations (5)-(8) as with − (E) related to + (E) by the detailed balance condition, − (E) = e −β E + (E), that can be seen in general from equations (5)- (7). Qualitatively, the rate (12) consists of the resonant peak at E = 0 broadened to width γ E C by the thermal noise in the resistor, and the 1/E tail on the 'emission' side (E > 0) of the peak. Strong variation of the tunneling rate with energy, from roughly max ∼ E 2 J /hγ at resonance to min ∼ E 2 J γ /hT E C away from the resonance, when E ∼ E C , creates different regimes in the system dynamics. If the rate of change of the gate voltage is very small, i.e.q min , the dynamics is adiabatic and the system maintains local equilibrium throughout the gate voltage ramp. In such an adiabatic evolution, the generated heat vanishes proportionally to drive rateq, not only on average but for individual evolution trajectories-see e.g. [18]. Also, in this quasistatic regime, the total evolution time is much larger than the correlation time of the occupation probabilities of the charge states, i.e. the total number of the back-and-forth transitions is large, and therefore there are no correlations between the initial and final occupations of the charge states. This means that if one neglects completely small adiabatically generated heat, the work distribution in this regime is given by the expression similar to equation (1) but with different intensities of the δ-peaks: For a more rapid drive,q ∼ max min , the total ramp time is sufficiently small to neglect completely transitions that could happen away from the resonance. In this regime, one can keep only the resonant peak in the rate (12), so that both rates ± coincide: With these rates, the solution of the master equation (4) with the initial condition p 0 (0) = 1 is Taking the integral of the rate (14) in this expression we obtain the LZ transition probability p 0 (t)| t→∞ ≡ P T in this 'high-temperature' regime: This is another well-known result for the dissipative LZ transitions [26,[30][31][32] reproduced here by a much more elementary means. In the resonant approximation (14) for the transition rates, used to obtain this result, all transitions happen at energies E ∼ γ T , i.e. the heat exchanged with the reservoir is small regardless of the drive rate. Neglecting this heat, one can see that the work distribution in this regime is again given by equation (1) but with modified LZ transition probability, P LZ → P T .
For parameter values away from the limiting cases, we have analyzed the distribution of work numerically by stochastic simulations tracing the paths, assuming rates of (5). Results of such simulations are shown in figure 2. For all these distributions JE is numerically valid within 0.1%, and it should naturally be valid identically since P(E), and thus the ± rates in (5) obey detailed balance [33]. These simulations correspond to the pumping trajectory described above, with parameters given in the figure 2 caption. The distributions have several features that can be explained by simple arguments as follows. First of all, the peaks in (a) and the steps in (b) at W = ±E C , and at 0 correspond to the peaks arising from excitation/relaxation, and from the adiabatic dynamics of the unitary evolution in (1), respectively, and they get gradually weaker toward faster relaxation (i.e. toward increasing R). In general, the mean of the work is positive, i.e. the evolution is dissipative. Furthermore, we have checked numerically that the Crooks equality [15] is valid in the form p(−W ) = e −βW p(W ) for these distributions.
Next we discuss the calorimetric measurement of the energy deposited on the resistor. The electronic system of the resistor forms the environment that is then coupled to a superbath formed, e.g. by phonons to which these electrons couple by standard electron-phonon coupling. This leads to the following sequence once an energy quantum E is emitted/absorbed by the two-level system due to ∓ transitions, see figure 3. We assume the standard hotelectron regime [34,35]. (i) The temperature T of the resistor (electrons) changes abruptly by T = ± E/C at the moment of the photon exchange; here C is the heat capacity of the electron gas in the resistor. (ii) The electronic temperature of the resistor starts to relax back toward the equilibrium T 0 ≡ (k B β) −1 at the rate governed by the thermal conductance G that is determined, e.g. by electron-phonon coupling. (iii) New energy pulses may occur during the evolution. In this evolution the deposited energy (heat) to the bath can be obtained generally fromQ = CṪ + G T , where T = T − T 0 . Note that the pulses can have either + or − sign, i.e. the resistor may cool as well. In the considered example, where E i is the magnitude of heat released in the ith absorption/emission event, and t i are the time instants of these events. Note that unlike in the analysis leading to (1) and (2), we allow here different values for E i for transitions occurring along the driven trajectory. Integrating the equation forQ for the situation where the system is in thermal equilibrium in the beginning and at the end, i.e. T (0) = T (∞) = T 0 , and assuming small temperature rise in the pulses, we have for the heat in each realization To make full use of (17), the following practical issues need to be taken care of [35]. The thermometer is first calibrated against the bath temperature under the equilibrium conditions (without applying the driving field), or one uses a primary thermometer that serves without calibration. The only remaining calibration is then the determination of G: this can be obtained by applying a known constant power P to the calorimeter, and measuring the temperature rise δT of the electrons in the thermometer, since G = P/δT . Before turning to the detailed practical realization and numbers, one further consideration is in order. Namely, biasing the gate capacitor in the ramp that we consider leads to a current that dissipates Joule power E d in the resistor R. This dissipation can be estimated by elementary circuit analysis with the result E d = e 2 R/τ . For typical values of τ and R, for instance those in figure 2, this energy is at most of order 0.01 E C . Since E C is the scale of work performed, the Joule heating produces only a small error in the calorimetric measurement, even if no precautions are taken.
We analyze next the feasibility of measuring calorimetrically the distribution of quantum work in a superconducting two-level system considered. The most straightforward realization would involve a small lithographic metal resistor, ∼(0.1 µm) 3 , on the chip [36]. The electronic heat capacity at T 0 = 100 mK temperature is C ∼ 10 −20 J K −1 , yielding a jump in temperature of E C /C ∼ 3 mK upon absorbing the relaxation heat at the end of the ramp, assuming a realistic E C /k B = 3 K. This is a sizable change of temperature (3%) which can be detected over a relatively long time of order 10 −4 s that is determined by the equilibration of the resistor back to its initial temperature via electron-phonon coupling. The large sensitivity of the calorimeter is naturally due to the standard T 0 /T F ∼ 10 −6 suppression of the heat capacity of the electron gas with Fermi temperature T F ∼ 10 5 K. If needed, longer relaxation times can be obtained by, for instance, suspending the resistor on a silicon nitride membrane [37], which would make the detection even easier. A tunnel probe, in form of a superconducting lead connected to the resistor, i.e. a superconductor-insulator-normal metal (SIN) junction, measured by RF reflection or transmission techniques [38], is fully compatible with the requirements of the presented calorimetric detection of heat.
In the current work we considered only one quantum system (CPB) in a specific (incoherent) regime as an illustration of the calorimetric method. Yet the same detection method can be applied to a large class of systems, where the validity of the fluctuation relations is not necessarily obvious unlike in the current example. First of all, the CPB in a weak dissipation regime is a qubit that obeys nearly coherent evolution and to be precise, an analysis respecting energies down to E J E C is needed. The calorimetric measurement applies, with equal conditions, also to this regime, but here the distribution of dissipation would have more structure at low energies and the theoretical assessment of the fluctuation relations would be less obvious. Secondly, the calorimetry can be applied to any system, where (i) the coupling between the quantum system and the calorimeter is the dominant energy relaxation channel, and where (ii) the sensitivity and the bandwidth of the thermometer are sufficiently high for observing the exchange of photons. These conditions are valid for mesoscopic electronic circuits quite generally. Instead of a CPB, one may envision other kinds of superconducting qubits [39], such as flux qubits (SQUIDs), where the coupling to the calorimeter can be provided inductively. Moreover, (semiconducting) quantum dot circuits [40] can be coupled capacitively to the dissipative calorimeter, similarly to what was presented for CPB here. This method can measure dissipation in many-body systems as well, as long as the systems can be coupled to a calorimeter. This way the number of possible applications of the method would be further increased [41][42][43].
In conclusion, the calorimetric measurement yields a direct determination of quantum work and quantum fluctuation relations. It can be used to analyze theoretically and experimentally the work and its distribution in a quantum system.