arxiv: v3 [quant-ph] 16 Feb 2016

Transcription

1 What is the probability of a thermoynamical transition? arxiv: v3 [quant-ph] 6 Feb 206 Álvaro M. Alhambra,, Jonathan Oppenheim,, 2, an Christopher Perry, Department of Physics an Astronomy, University College Lonon, Gower Street, Lonon WCE 6BT, Unite Kingom 2 Department of Computer Science an Centre for Quantum Technologies, National University of Singapore, Singapore 965 If the secon law of thermoynamics forbis a transition from one state to another, then it is still possible to make the transition happen by using a sufficient amount of work. But if we o not have access to this amount of work, can the transition happen probabilistically? In the thermoynamic limit, this probability tens to zero, but here we fin that for finite-size systems, it can be finite. We compute the maximum probability of a transition or a thermoynamical fluctuation from any initial state to any final state, an show that this maximum can be achieve for any final state which is block-iagonal in the energy eigenbasis. We also fin upper an lower bouns on this transition probability, in terms of the work of transition. As a bi-prouct, we introuce a finite set of thermoynamical monotones relate to the thermo-majorization criteria which governs state transitions, an compute the work of transition in terms of them. The trae-off between the probability of a transition, an any partial work ae to ai in that transition is also consiere. Our results have applications in entanglement theory, an we fin the amount of entanglement require or gaine when transforming one pure entangle state into any other. I. INTRODUCTION Given a quantum system in a state ρ with some Hamiltonian, H, when can it be eterministically transforme into another state σ associate with a potentially ifferent Hamiltonian, H 2? If we can put the system into contact with a heat bath at temperature T, then in the thermoynamical limit, an if interactions are short-range or screene, a transition will occur as long as the free energy of the initial configuration is larger than the free energy of the final configuration. The free energy of the state ρ efine as: F ρ, H = tr [H ρ] T S ρ, were Sρ is the entropy; Sρ = tr ρ log ρ. This is a formulation of the secon law of thermoynamics, if we factor in energy conservation the first law. If we wish to make a forbien transition occur, then we nee to inject an amount of work which is greater than the free energy ifference between initial an final states. However, what if we are intereste in small, finitesize systems? Or in systems with long-range interactions? The thermoynamics of systems in the microregime, where we o not take the thermoynamical limit, has gaine increase importance as we cool an manipulate smaller an smaller systems at the nano scale an beyon [ 5]. Theoretical work has continue a pace, with increase interest in the fiel in recent years [6 32]. If we o not take the thermoynamical limit, then provie σ is block-iagonal in the energy eigenbasis, there is not just one criteria the ecreasing of the free energy, but a family of criteria which etermine whether alvaro.alhambra.4@ucl.ac.uk j.oppenheim@ucl.ac.uk christopher.perry.2@ucl.ac.uk a state transition is possible. A set of such criteria which have been proven to be necessary an sufficient conition for quantum thermoynamical state transformations [20] c.f. [8], are the so-calle thermo-majorization criteria [7, 20]. Thermo-majorization is a set of conitions that are more stringent than the orinary secon laws an ha been conjecture to provie a limitation on the possibility of thermoynamical transformations since 975 [7]. It is relate [8, 33] to a conition known as Gibbs-stochasity[34, 35] a conition which can be extene to inclue fluctuations of work [36]. Once again though, if the iagonal state σ, is not thermo-majorize by ρ, then a transition is still possible, provie sufficient work is use. One can compute the work require or gaine from this transition using thermo-majorization iagrams [20], via a linear program[33], or the relative-mixeness [24]. Suppose however, we want to make a transition from ρ to σ, an it requires work which we cannot, or o not wish to, expen. Can we still nonetheless make the transition with some probability p rather than with certainty? An if so, what is the highest probability, p, that can be achieve? In particular, given ρ an σ, we are intereste in maximizing p in the following process: ρ ρ = pσ + p X, 2 with X being some arbitrary state. Such a transformation can be regare as a fluctuation of a system s state, in the sense that the transformation is only probabilistic. Within the stuy of thermoynamics for small systems, great progress has alreay been mae in analyzing how the work istribution associate with a given transformation of process can fluctuate [37 40] see [4 43] for reviews on both the classical an quantum cases. Fluctuation relations such as the Jarzynski equality [37] an Crooks theorem [38], evelope uner the paraigm of stochastic thermoynamics, have been use

2 2 to calculate the work fluctuations of non-equilibrium processes. Investigating fluctuation in a system s state provies a natural, complementary stran of research which we are able to formulate an analyze in this paper by applying techniques from quantum information theory evelope in [20]. In relate work [36], we shall aress the problem of fluctuating work within this information theoretic framework. This shall serve to unify the two approaches to thermoynamics for small systems an exten an provie insight into previous work base on the stochastic thermoynamics perspective. Here, we will upper boun the maximum probability of a fluctuation between any given ρ an σ. When σ is block-iagonal in the energy eigenbasis, we will show that this boun can be achieve an furthermore, that there exists a two outcome measurement that can be performe on ρ such that we obtain σ with the maximum probability p. Of course, measurements o not come for free in thermoynamics - it costs work to erase the recor of the measurement outcome [44]. That this measurement can be performe is note for completeness however, we take Eq. 2 as our primary goal, efining what we mean by a thermoynamical transition. We will iscuss measurements in Section III B as they only provie a small correction of kt log 2 to the work cost of a probabilistic transformation. Our main result will be Theorem 5, which upper bouns the probability p in terms of a minimization over a finite set of ratios between thermoynamical monotones, which are quantities that can only ecrease uner the set of allowe operations. When the final state is block-iagonal, this boun is achievable, but this may not be the case if the final state has coherences in energy. These monotones, which we will show are given by Eq. 4, can be thought of as analogous to free energies.this is proven in Theorem 4 an is equivalent to the thermo-majorization criteria of [7, 20]. The set of ratios that we use to boun p thus gives an alternative way of verifying if the thermo-majorization criteria are satisfie. Rather than consiering the thermo-majorization curves [20] or consiering a continuous set of monotones [24] we provie a finite set of conitions to check. Inee this set provies a strengthening of results from the theory of relative majorization [45, 4.B.4c] by reucing the number of constraints that nee to be consiere. Before proving Theorem 5, we will consier in Section II the simpler case where the Hamiltonian of the system is trivial, i.e. H I. Solving the problem in this regime, referre to as Noisy Operations [5, 46], will provie us with insight into the solution for non-trivial Hamiltonians. In this simplifie situation, p is given by Theorem. The result is similar in form to [47] which consiers the analogous problem of probabilistic pure state entanglement manipulation using Local Operations an Classical Communication LOCC. However, care must be taken the class of operations allowe uner LOCC is very ifferent to what is allowe in thermoynamics. For example, uner LOCC one can bring in pure states for free which can be a source of work in thermoynamics an one is allowe to make measurements for free which costs work. Perhaps more importantly, many of the LOCC monotones are concave, which is not the case in Noisy Operations, thus we will require some ifferent techniques. It shoul also be note that in entanglement manipulation, the maximum probability achievable will be zero if the target state has a larger Schmit rank than the starting state. Uner Noisy Operations, we will see that p is always non-zero though it can be arbitrarily small. In Section III we consier the general case of arbitrary initial an final Hamiltonians an states. We will prove our results using the paraigm of Thermal Operations TO [20, 34, 35]. There are a number of ifferent paraigms one can use to stuy thermoynamics e.g. allowing interaction Hamiltonians or changing energy levels, however, these other paraigms are equivalent to Thermal Operations [20, 48], an thus Thermal Operations are the appropriate paraigm for stuying funamental limitations. We introuce Thermal Operations at the beginning of Section III. In the case of a trivial Hamiltonian, Thermal Operations reuce to Noisy Operations, the regime consiere in Section II. Our expression for the cost of a transition between any two states using only a finite number of monotones is given in Lemma 2 for Noisy Operations an Lemma 6 for Thermal Operations. The Noisy Operations result can be aapte to give an expression for the amount of entanglement require or gaine when transforming any pure bipartite state into another uner LOCC. This is given in Appenix A an generalizes existing expressions for the istillable entanglement [49, 50] an cost of entanglement formation [5]. We also show how p can be upper an lower boune using the work of transitions from ρ to σ an σ to ρ. This is one in Lemma 3 for the case of a trivial Hamiltonian, an in Lemma 7 for the general case. Finally, we conclue in Section IV with a iscussion on other goals, relate to Eq. 2, which one coul attempt when making a probabilistic transition. One such goal, the optimization of the herale probability, is iscusse in etail in Appenix B where we obtain bouns on it, even in the presence of coherence or catalysts. The heraling probability can be thought of as a generalization of the case where one achieves Eq. 2 with a measurement i.e. ρ 0 0 T O ˆρ = pσ p X. an the transition is conclusive. This allows us to analyze state fluctuations in the presence of measurements, coherence an catalysis. We also pose some open questions. One of these regars how p varies if we supply aitional work to rive the transition from ρ to σ or eman that aitional work be extracte. The solution for qubit systems with trivial Hamiltonian is given in Appenix C.

3 3 II. PROBABILITY OF TRANSITION UNDER NOISY OPERATIONS Before investigating Eq. 2 in the context of Thermal Operations, we will first consier a simpler, special case - Noisy Operations. In this particular instance of thermoynamics, the Hamiltonian of the system uner consieration is trivial. Noisy Operations were first efine in [5] where the problem of whether a transition between two given states uner a particular set of operations was consiere. Within Noisy Operations, the following actions are allowe: i a system of any imension in the maximally mixe state can be ae, ii any subsystem can be iscare through tracing out an iii any unitary can be applie to the global system. Throughout this paper, we shall use η i to enote the eigenvalues of ρ an ζ i to enote those of σ. For a comprehensive review of Noisy Operations, see [46]. Given two states, ρ an σ, it was shown in [5] that transition from ρ to σ is possible uner Noisy Operations if an only if ρ majorizes σ written ρ σ. That is, if we list the eigenvalues of ρ an σ [52] in ecreasing orer an enote these orere lists by η = {η,..., η n } an ζ = {ζ,..., ζ n } respectively, the transition is possible if an only if: where: V l ρ V l σ, l {,..., n}, 3 V l ρ = l η i. 4 Lorenz curves are a useful tool for visualizing these criteria Figure. For a given state ρ, its Lorenz curve is forme by plotting the points: { k k n, η i}n, 5 k= an connecting them piecewise linearly together with the point 0, 0 to form a concave curve. If ρ majorizes σ, the Lorenz curve for ρ is never below that of σ. The functions efine in Eq. 4, an their analogue in Thermal Operations, will be crucial for the rest of the paper. They are monotones of the theory, only ecreasing uner Noisy Operations. Excellent reviews regaring the theory of majorization an Lorenz curves can be foun in [45, 46]. A. Non-eterministic transitions We will now consier transitions when the conitions given in Eq. 3 are not necessarily fulfille. Here, rather than transforming ρ to σ with certainty, we shall o so with some probability as formulate in Eq. 2. In particular, we are intereste in the maximum probability, 0 Elbow 2 I Σ Σ Ρ s I Σ s log Maximally mixe FIG.. Lorenz Curves. { a The Lorenz curve for ρ is efine k by plotting the points:, } n k n ηi b The transition k=. from σ to ρ is possible uner NO as the curve for σ is never below that of ρ. c The Lorenz curve for a maximally mixe state is given by the ashe line from 0, 0 to,. All other states majorize it. s log 5 is an example of a sharp state. e 4 s I σ is the least sharp state that majorizes σ. p, that can be achieve. A similar problem is consiere in [47] for entanglement manipulation an aapting its techniques the following theorem can be shown: Theorem. Suppose we wish to transform the state ρ to the state σ uner Noisy Operations. The maximum value of p that can be achieve in the transition: is given by: ρ NO ρ = pσ + p X, 6 p = V l ρ min l {,...,n} V l σ. 7 Proof. The proof is split into two parts: first we apply Weyl s inequality an the efinition of majorization to erive a contraiction if it were possible to achieve a value of p large than p. Next, we aapt the techniques of [47] to provie a protocol achieving p = p. To achieve our first goal we begin by showing that given Eq. 6: V l ρ pv l σ, l. 8 To prove this, we make use of Weyl s inequality [53, 54]. Given n n Hermitian matrices, A, B an C such that A = B + C, let {a i } n, {b i} n an {c i} n be their respective eigenvalues arrange in escening orer. Weyl s inequality then states that: b i + c n a i b i + c, 9 for all i. Applying this to ρ, σ an X, we obtain: η i pζ i + p x n, i, 0

4 4 where x n is the smallest eigenvalue of X. As X is a positive semiefinite matrix, x n 0 an: Hence: V l ρ V l ρ = η i pζ i, i. l η i p l ζ i = pv l σ, 2 where the first inequality uses Eq. 3 an the secon follows from Eq.. Now suppose it was possible to achieve a value of p greater than p in Eq. 7. Then there woul exist an l such that V l ρ < pv l σ, contraicting Eq. 8. To show that p is obtainable, we efine the following quantities. First, efine l by: { l = max l : V } l ρ V l σ = p r. 3 Then we procee iteratively an, provie l i < n, efine: so we have: r i r i V l ρ V li ρ = min l>l i V l σ V li σ, 4 l j=l i + ζ j l j=l i + η j, l > l i. 5 Define l i by: { l i = max l : l > l i, V } l ρ V li ρ V l σ V li σ = ri. 6 Note that we have r i > r i. To see this, first observe that for a, b, c, > 0: Setting: a b < a + c b + a b < c. 7 a = V li ρ V li 2 ρ, b = V li σ V li 2 σ, c = V li ρ V li ρ, = V li σ V li σ, so a b = ri an c = ri, then: a + c b + = V l i ρ V li 2 ρ V li σ V li 2 σ > ri = a b, where the inequality follows from the efinition of r i. Using Eq. 7, the claim that r i > r i now follows. Overall, this protocol generates a set of l i such that 0 = l 0 < l <... < l k = n an a set of r i such that p = r <... < r k. Now we split ρ an σ into blocks an efine: ρ i = iag η li +,..., η li, 8 σ i = iag ζ li +,..., ζ li. 9 Then from Eq. 5 an the fact that equality occurs when l = l i, ρ i majorizes r i σ i an we can perform: NO ρ i r i σ i = p σ i + r i p σ i, i. 20 With a bit of massaging an recombining the blocks, this is the same form as Eq. 6 with p = p an the blocks of X being efine by: X i = ri p p σ i. 2 Note that as the enpoints of the Lorenz curves coincie at, an η > 0, we are guarantee that 0 < p. If we want to obtain σ from ρ with probability p rather than have it as part of a probabilistic mixture as per Eq. 6, we can o so by performing a two outcome measurement, with measurement operators { M, I M}, where the blocks of M are given by: p p M i = iag,...,. 22 r i r i To see that M is a vali measurement, we note that in general 0 < p. Hence both { M, I M} are r well efine, an i their squares trivially a up to the ientity. After applying this measurement to ρ an reaing the result, we will have either: or M ρ M = p σ, 23 I M ρ I M = p X. 24 However, performing this measurement is outsie of the class of Noisy Operations an hence costs work. As such, if a general two outcome measurement is allowe without taking its cost into account, it can be possible to transform ρ into σ with probability greater than p. For example, if ρ an σ are qubits, we can convert ρ into σ with certainty using this extra resource. Firstly we a an aitional qubit in the maximally mixe state an measure it in the computational basis. This results in a pure state, either 0 or. As these majorize all other qubit states we can use it to obtain any σ with certainty.

5 5 B. Nonuniformity of transition uner Noisy Operations If it is not possible to eterministically convert ρ into σ using Noisy Operations, to perform the transformation with certainty will cost some resource, in the form of nonuniformity. For instance, if we a some pure states of sufficiently high imension, a previously impossible transition will become possible. Aing these aitional pure states can be thought of as the analogue to aing work. Similarly, if ρ can be converte into σ using Noisy Operations, it may be possible to extract some nonuniformity e.g. by transforming some maximally mixe states into pure states. This is the analogue of extracting work. More generally, we shall extract or expen the equivalent of work using sharp states. These sharp states, as iscusse in the next subsection, will serve as a natural unit for the nonuniformity resource. We will compute the nonuniformity of transition in terms of a finite set of ratios of monotones. This is one in a similar manner to [46], although we show that the minimization can be one over fewer points.. Sharp States Quantifying the optimal amount of work of transition for the more general Thermal Operations was consiere in [20, 24]. We shall enote the Noisy Operations equivalent of work, the nonuniformity of transition, by I ρ σ. If nonuniformity must be ae, the quantity is negative, while if we can extract nonuniformity, it will be positive. For I ρ σ = log j, we efine an associate sharp state [46] by: s Iρ σ = iag j,...,, 0,..., j }{{}}{{} j j Appening a sharp state I log to the system is equivalent j to introucing log j units of nonuniformity. See Figure for an example of a sharp state s Lorenz curve. The state s Iρ σ is such that: NO ρ s Iρ σ σ, if I ρ σ 0, ρ NO σ s Iρ σ, if I ρ σ > In terms of Lorenz curves, tensoring a state ρ with a sharp state s I has the effect of compressing the Lorenz curve of ρ by a factor of 2 I with respect to the x-axis [20]. 2. Monotones for Noisy Operations an the nonuniformity of transition The function V l ρ is equal to the height of the Lorenz curve of ρ at x = l n. An alternative set of monotones, L y ρ where 0 y, can be efine as the shortest horizontal istance between the Lorenz curve of ρ an the y-axis at y. Note that these functions never ecrease uner Noisy Operations. In particular: L yk ρ = k n, for y k = L ρ = rank ρ. n k η i, k < rank ρ, If we efine the set D σ by: 27 { k } rankσ D σ = ζ i, 28 k= then a transition from ρ to σ is achievable with certainty uner Noisy Operations if an only if: L y ρ L y σ, y D σ. 29 That it is sufficient to consier only y D σ will be justifie below. The horizontal monotones, L y, also allow us to quantify the optimal work of transition that is require or extracte in going from ρ to σ: Lemma 2. Given two states ρ an σ, uner Noisy Operations: L 2 Iρ σ y ρ = max y Dσ L y σ. 30 Proof. To prove this, we make use of the geometrical structure of Lorenz curves an the properties of I ρ σ. Note that we have: L 2 Iρ σ y ρ = max y [0,] L y σ, 3 as this follows from the fact that to obtain the optimal value of I ρ σ, we wish to rescale the Lorenz curve of ρ with respect to the x-axis in such a way that it just majorizes that of σ - the curves shoul touch but not cross. The amount that we nee to rescale by is given by Eq. 3. We now show that it is sufficient to maximize over y D σ. Let s 0 = 0 an s k = k ζ i for k rank σ. Then, for j rank σ, as the Lorenz curve of σ is a straight line on the interval [s j, s j ] an the Lorenz curve of ρ is concave: L y ρ max y [s j,s j] L y σ max r [0,] rl sj ρ + r L sj ρ r j n + r j. n 32 It is straightforwar to check that the maximum value occurs at either r = 0 or r =. We can thus replace the inequality in Eq. 32 with an equality an it follows that it suffices to maximize over y D σ.

6 6 As ρ NO σ is possible if an only if I ρ σ 0, the finite set in Eq. 29 is justifie. Note that in [46] it was shown that it is possible to calculate I ρ σ by performing an optimization over the ratios calculate at the elbows see Figure for a efinition of both ρ an σ. In Lemma 2 we have shown that it suffices to consier just the elbows of σ. 3. Bouns on the transition probability The quantities I ρ σ an I σ ρ can be use to boun p as follows: Lemma 3. Given two states ρ an σ, uner Noisy Operations: 2 Iρ σ p 2 Iσ ρ, 33 where as p, we assume I ρ σ 0. If I ρ σ 0, p = an the transformation from ρ to σ can be one eterministically, potentially extracting a finite amount of nonuniformity. Proof. We start proving with the lower boun, giving a protocol which achieves p = 2 Wρ σ. The upper boun is erive by consiering properties of the purity of the least sharp state that majorizes ρ. Assuming W ρ σ = log j for simplicity, an efining I to be the maximally mixe state of a level system: ρ ρ NO I, 34 = j ρ s log + j j ρ s log, j NO j σ I + j Y, NO j σ + j Tr BY, where Y is the state obtaine by applying the secon Noisy Operation to ρ s log. Using this protocol, we j obtain something of the form Eq. 6 with p = 2 Iρ σ an X = Tr B Y. As p is the maximum value of p obtainable in Eq. 6, we erive the lower boun. We now consier the upper boun an to obtain a useful boun, assume I σ ρ > 0. We efine I ρ as the nonuniformity of formation of ρ uner NO[8], given by I ρ = log η n, an hence let s I ρ be the least sharp state that majorizes ρ see Figure. Note that I ecreases uner Noisy Operations an is aitive across tensor proucts [46]. In terms of the eigenvalues of ρ an σ: By efinition, as I σ ρ > 0: s I ρ = s logηn, s I σ = s logζn. 35 σ NO ρ s Iσ ρ. 36 Now, using first the monotonicity of I an then the aitivity: I σ I ρ siσ ρ, monotonicity = I ρ + I σ ρ. aitivity I σ ρ I σ I ρ, = log ζ n log η n, ζ = log. 2 Iσ ρ η ζ, as require. η = V ρ V σ, p, by efinition From Eq. 33 we can see that when I ρ σ = I σ ρ I that is, in a reversible transition then p = 2 I. This occurs when either σ NO = ρ s I or ρ NO = σ s I epening on whether I is positive or negative when I 0 the transition is eterministic. In terms of Lorenz curves, this means that the curves of ρ an σ have the same shape up to re-scaling by a factor 2 I. In particular, this is the case when both ρ an σ are sharp states, where both Lorenz curves are straight lines. This result can be applie in the thermoynamic regime of many inepenent copies. If we want to perform a transition such as: ρ N σ N, 37 we nee an amount of work given by NI ρ σ. Hence, the probability of success in such a case is boun by: 2 NIρ σ p 2 NIσ ρ, 38 which tens to 0 for large N. This can be seen as a way in which in the thermoynamic limit statistical fluctuations are suppresse. 4. Lorenz curve interpretation In terms of Lorenz curves, aing I ρ σ nonuniformity to ρ to make the transition possible is equivalent to compressing the Lorenz curve with respect to the x-axis by a ratio 2 Iρ σ, such that the curve of ρ lies just above an touches that of σ. Hence, a compression by p 2 Iρ σ must mean that there is at least a point of the compresse curve just below or touching σ. A proof of this is given in Figure 2. Extracting I σ ρ nonuniformity from σ before performing NO into ρ is equivalent to compressing the curve of ρ by a ratio of 2 Iσ ρ such that the curve of σ lies just above an touches that of ρ. Hence, to prove the upper boun in Eq. 33, it suffices to show that in compressing

7 7 A is lower an upper boune respectively by the minimum an the maximum horizontal ratio of the two. 0 6 C 3 B Σ Ρ Ρs log p FIG. 2. We plot the curves of ρ, σ an ρ compresse by p with respect to the x-axis. The points A an B at which the vertical ratio between the curves of ρ an σ is maximum which sets l an p, an the sharp states that pass through those points are also shown as ashe lines. After compressing the Lorenz curve of ρ by a ratio of p, the point B will be taken to C, which will always either be below the curve of σ or just touching it. This proves the lower boun in Eq. 33. C' A' B' A B Σ Ρ 5 6 Ρs log p 0 FIG. 3. We plot the curves of ρ, σ an ρ compresse by p with respect to the x-axis. The points A an B at which the vertical ratio between the curves of ρ an σ is maximum which sets l an p an the sharp states I ρ an I σ are also shown as ashe lines. Given that for sharp states all bouns are saturate, the appropriate maximum vertical an horizontal ratios coincie, an are η /ζ, the ratio of the heights of B an A. But this ratio is, by efinition, bigger than or equal to p, the ratio between A an B. This means that if the curve of ρ is compresse by p, the point B is mappe to C just above or touching the curve of σ, proving the upper boun of Eq. 33. the curve of ρ by p at least one point of the new curve must lie above or touch that of σ. In Figure 3 we show a iagrammatic version of the proof given in Section II B. It shoul be note that with Lemma 3 we are proving a general statement about convex Lorenz curves. This is, that the minimum vertical ratio of two given curves p III. PROBABILITY OF TRANSITION UNDER THERMAL OPERATIONS Noisy Operations can be generalize to inclue systems with arbitrary, finite Hamiltonians. This is the resource theory of Thermal Operations [20, 34, 35, 48]. Within this scheme, the allowe operations are: i a system with any Hamiltonian in the Gibbs state of that Hamiltonian can be ae, ii any subsystem can be iscare through tracing out an iii any energy-conserving unitary, i.e. those unitaries that commute with the total Hamiltonian, can be applie to the global system. These operations moel the thermoynamics of a system in the presence of an ieal heat bath [20, 48]. Note that while the heat bath the system is in contact with is assume to be large, thermal operations inclue processes that only interact with a small part of the bath. As such, limitations erive with respect to such an iealize bath can be regare as truly funamental. Even though the bath size can be large, the system of interest is fixe, an can for example, be only a single system. They thus escribe processes beyon the thermoynamic limit. In general, the initial an final systems may have ifferent Hamiltonians but, by making use of the switching qubit construction in [20], we can w.l.o.g. assume that the initial an final Hamiltonians are the same. As such, the results in this section will assume this but in Section III D we will iscuss how a changing Hamiltonian affects them. In Appenix H of [48] it was shown that other mainstream thermoynamical paraigms such as time epenent Hamiltonians, the insertion of interaction terms between system, bath an work systems an various master equations are all inclue within the scope of Thermal Operations. In the absence of catalysts, an provie the final state is block-iagonal in the energy eigenbasis, it was establishe in [20] that a transition from ρ to σ is possible uner Thermal Operations if an only if ρ thermomajorizes σ. This is similar in form to the majorization criteria of Noisy Operations an can be visualize in terms of thermo-majorization iagrams which are similar to Lorenz curves but with two crucial ifferences. Suppose ρ is also block-iagonal in the energy eigenbasis with eigenvalue η i associate with energy level E i, for i n. Firstly, rather than orering accoring to the magnitue of η i, we instea β-orer them, listing η i e βei in escening orer. The secon ifference is that we no longer plot the β- orere η i at evenly space intervals. Instea we plot the points: { k e βeρ i, k η ρ i } n k=, 39

8 8 Σ Ρ Pure State Thermal State Β E Β E Β E 2... Β E k FIG. 4. We show the β-orere thermo-majorization iagrams for various states of the system. Note that ifferent states may have ifferent β-orerings an the markings on the x-axis correspon to one particular β-orering. The curves always en at Z,. The thermo-majorization criteria states that we can take a state to another uner Thermal Operations if an only if the curve of the initial state is above that of the final state. Hence, in this case provie ρ is block-iagonal in the energy eigenbasis there is a set of operations such that σ T O ρ, but not for the reverse process. where the superscript ρ on E i an η i inicates that they have been β-orere an this orering epens on ρ. Thermo-majorization states that ρ can be eterministically converte into a block-iagonal σ if an only if its thermo-majorization curve never lies below that of σ, as is shown in Figure 4. This is analogous to the case of Noisy Operations. In what follows, we assume that the η i have been β-orere unless otherwise state. If ρ is not block-iagonal in the energy eigenbasis, to etermine if a transition is possible we consier the thermo-majorization curve associate with the state forme by ecohering ρ in the energy eigenbasis. This state, ρ D, is given by: ρ D = n E i E i ρ E i E i, 40 where E i is the eigenvector of the system s Hamiltonian associate with energy level E i. The operation of ecohering ρ to give ρ D is a Thermal Operation an commutes with all other Thermal Operations [48]. A transition from ρ to σ, where σ is block-iagonal in the energy eigenbasis, can be mae eterministically if an only if the thermo-majorization curve of ρ D is never below that of σ. Finally, if σ is not block-iagonal, a transition from ρ to σ is possible only if ρ D thermo-majorizes σ D an fining a set of sufficient conitions is an open question. In what follows, the thermo-majorization curve of a state with coherences is efine to be the thermomajorization curve of that state ecohere in the energy eigenbasis as per Eq. 40. Similarly to how Eq. 4 efines monotones for the Noisy Operations resource theory, the height of the β- orere thermo-majorization curves provies monotones for Thermal Operations. If we enote the height of the thermo-majorization curve of ρ at x by Ṽx ρ, for 0 x Z where Z is the partition function, then by the thermo-majorization criteria, this function is nonincreasing uner Thermal Operations. In particular, for block-iagonal ρ, we have: Ṽ xk ρ = k η ρ i, where x k = k e βeρ i. 4 These monotones also give us an alternative way of stating the thermo-majorization criteria: Theorem 4. Suppose σ is block-iagonal in the energy eigenbasis. Let L σ = { k } n e βeσ i k=. Then ρ can be eterministically converte into σ uner Thermal Operations if an only if: Ṽ x ρ Ṽx σ, x L σ. 42 Proof. To prove this theorem, we make use of the concavity properties of thermo-majorization curves. Suppose ρ T O σ. Then by thermo-majorization, Ṽx ρ Ṽx σ, for 0 x Z an in particular Eq. 42 hols. Conversely, suppose Eq. 42 hols an, setting t 0 = 0, label the elements of L σ arrange in increasing orer by t i for i = to n. Then on the interval [t i, t i ], for i n, the thermo-majorization curve of σ is given by a straight line. From ρ, efine the block-iagonal state ρ σ by the thermo-majorization curve: { } n t i, Ṽt i ρ, 43 an note that ue to the concavity of thermomajorization curves, ρ thermo-majorizes ρ σ. On the interval [t i, t i ], i n, the thermo-majorization curve of ρ σ is also given by a straight line. The construction of this state ρ σ is shown in Figure 5. As Ṽt i ρ σ = Ṽt i ρ, i by construction, Eq. 42 implies that Ṽt i ρ σ Ṽt i σ, i. Hence on the interval [t i, t i ], i n, the thermo-majorization curves for ρ σ an σ, an therefore ρ an σ, o not cross. As this hols for all i an the intervals cover [0, Z] the thermomajorization curve of ρ is never below that of σ an we can perform ρ T O σ eterministically. If we efine the number of elbows in the thermomajorization curve of σ to be j, this reuces thermomajorization to checking j criteria an generalizes Lemma 7 of [46] to Thermal Operations. Note also that if σ is not block-iagonal in the energy eigenbasis, Eq. 42 gives a necessary but not sufficient conition for the transition from ρ to σ to be possible.

9 9 Β E Β E Β E 2... Β E k Ρ Σ Ρ Σ Thermal state FIG. 5. Here we illustrate the construction of the state ρ σ use in the proof of Theorem 4. The points of the curve ρ that are at the same horizontal position as the elbows of σ are joine, an by concavity the resultant curve is always below ρ. A. Non-eterministic transformations Having efine the appropriate monotones for Thermal Operations, we are now in a position to investigate non-eterministic transformations an prove a theorem analogous to Theorem. Theorem 5. Suppose we wish to transform the state ρ to the state σ uner Thermal Operations. The maximum value of p, p, that can be achieve in the transition: is such that: ρ TO ρ = pσ + p X, 44 p Ṽ x ρ min x Lσ Ṽ x σ. 45 Furthermore, if σ is block-iagonal in the energy eigenbasis, there exists a protocol that achieves the boun. Proof. Proving this result is more complicate than proving Theorem ue to the fact that ρ an σ may have ifferent β-orerings. We procee as before, first showing the boun in Eq. 45 an then giving a protocol that achieves the boun when σ is block-iagonal. We prove the boun in Eq. 45 by constructing useful intermeiate curves between those of ρ an pσ to eal with iffering β-orers. With these in place, the result will follow in a similar manner to Theorem. We begin by showing that given Eq. 44: Ṽ x ρ pṽx σ, x [0, Z]. 46 First consier for general σ the maximum value of p that can be achieve in attempting to convert ρ into σ. As ecohering is a Thermal Operation, this value of p can also be achieve when attempting to convert ρ into σ D : ρ TO ρ = pσ + p X, ecohere ρ D = pσ D + p X D. Thus, to upper boun p, it suffices to show that Eq. 45 hols for block-iagonal σ. Furthermore, w.l.o.g. we can assume that ρ an X are also block-iagonal. Using Weyl s inequality as per Theorem to eal with egenerate energy levels, for block-iagonal ρ, σ an X, we have: η i pζ i, i. 47 Now consier the sub-normalize thermo-majorization curve of pσ given by the points: { k e βeσ i, p k ζ σ i } n k=, 48 an the possibly non-concave curve forme by plotting the eigenvalues of ρ accoring to the β-orering of σ. This is given by the points: { k e βeσ i, k η i σ } n k=. 49 By Eq. 47, the curve efine in Eq. 49 is never below that efine in Eq. 48. Finally, the thermo-majorization curve of ρ is given by: { k e βeρ i, k } n η i ρ k=. 50 Note that attempting to construct a thermo-majorization curve for ρ with respect to the β-orering of another state, as we o in Eq. 49, has the effect of rearranging the piecewise linear segments of the true thermomajorization curve. This means that they may no longer be joine from left to right in orer of ecreasing graient. Such a curve will always be below the true thermomajorization curve. To see this, imagine constructing a curve from the piecewise linear elements an in particular, trying to construct a curve that woul lie above all other possible constructions. Starting at the origin, we are force to choose the element with the steepest graient - all other choices woul lie below this by virtue of having a shallower graient. We then procee iteratively, starting from the enpoint of the previous section ae an choosing the element with the largest graient from the remaining linear segments. The construction that we obtain is the true thermo-majorization curve. A graphical escription of this proof is shown in Figure 6.

10 0 Β E Β E Β E 2... Β E k FIG. 6. Here we show graphically the steps of the proof of the first part of Theorem 5. In the ecomposition of Eq. 44 the curve pσ must always be below that of ρ an hence also ρ. This sets the maximum probability p as efine in Eq. 45. Both pσ an the isorere ρ have the same β-orering. As such, the curve in Eq. 50 is never below that in Eq. 49. This gives us: Ṽ x ρ Ṽx ρ pṽx σ, 5 where the first inequality hols as, by efinition, ρ thermo-majorizes ρ. In particular we have: p Ṽ x ρ min x Lσ Ṽ x σ. 52 When σ is block-iagonal in the energy eigenbasis, a protocol that saturates the boun is: ρ TO ρ σ, TO ρ σ = p σ + p X, where ρ σ was efine in Eq. 43 an is thermomajorize by ρ. As ρ σ an σ have the same β-orering an: Ṽ x ρ Ṽ x σ = Ṽx ρ σ, x L σ, 53 Ṽ x σ applying the same construction use in Theorem gives a strategy to prouce ρ σ that achieves: p = Ṽ x ρ min x Lσ Ṽ x σ. 54 B. Measuring whether the transition occurre uner Thermal Operations For block-iagonal σ, after obtaining ρ through Thermal Operations we may apply the measurement efine by Eq. 22 to extract our target state with probability p. This can be one through a process that uses an ancilla qubit system, Q, that starts an ens in the state 0 an has associate Hamiltonian, H Q = I 2, a unitary that correlates the system with the ancilla an a projective measurement on the ancilla qubit. As the measurement operators are iagonal in the energy eigenbasis, we will fin that the unitary is energy conserving an within the set of Thermal Operations. Furthermore, the ancilla that is use to perform the POVM can be returne back into it s original state. Hence the only cost we have to pay is to erase the recor of the measurement outcome itself. As is well known [55], the cost of erasing the recor is kt log 2, although if one is repeating the process many times, then it is kt h p with h p the binary entropy h p = p log p p log p [23]. The unitary that we shall use is given by: M I M U SQ =, 55 I M M where M is efine as per Eq. 22. Note that U SQ = U SQ. Its effect on the initial joint state is: U SQ ρ 0 0 U SQ, M I M = ρ 0 M I M, I M M 0 0 I M M Mρ M Mρ I M = I Mρ M I Mρ, I M p = σ Mρ I M I Mρ M p. X If we now measure the ancilla in the computational basis, the joint state will collapse to σ 0 0 when the 0 outcome is observe. This happens with probability p. If the outcome is observe, the joint state collapses to X an this happens with probability p. In aition, if the outcome is observe, we can then apply a Pauli Z to the ancilla qubit to return it to its initial state. To see that U SQ commutes with the total Hamiltonian an belongs to the class of Thermal Operations, first note that the total Hamiltonian is given by: H SQ = H S I 2 + I n I The unitary trivially commutes with the secon term so focusing on the first term, an noting that M an H S are both iagonal matrices so commute, it is easy to check that: M I M [U SQ, H S I 2 ] = HS 0 I M M 0 H S HS 0 M I M, 0 H S I M M = 0.

11 Hence [U SQ, H SQ ] = 0. Observe that this reasoning can be generalize to measurements with s outcomes [56]. Provie the measurement operators commute with H S, the measurement can be performe using a s-level ancilla system with trivial Hamiltonian an a joint energy-conserving unitary. Such a measurement can be performe for free up to having to spen work to erase the recor of the measurement outcome at a cost of kt ln s. On the other han, channels that are not compose of Thermal Operations incluing some measurements characterize by non-iagonal operators can be seen as a resource [57]. Ρ WW Ρ 00 Thermal state C. Work of transition uner Thermal Operations. Work systems In general, if we want a transition ρ σ to be possible, work may have to be supplie. Alternatively, if a transition can be achieve with certainty, it may be possible to extract work. For the thermoynamics of small systems, the concept of eterministic work also referre to in the literature as single-shot or worst-case work has been introuce [8, 20, 2]. Within the Thermal Operation paraigm, the optimal amount of work that must be ae or gaine can be quantifie using the energy gap, W, of a 2-level system with groun state 0 an excite state W with energy W. The associate Hamiltonian is: H = W W W. 57 The work of transition, W ρ σ, is such that: if W ρ σ 0, ρ W ρ σ W ρ σ T O σ 0 0, if W ρ σ > 0, ρ 0 0 T O σ W ρ σ W ρ σ. 58 Defining work in such a way enables the quantification of the worst-case work of a process. When W ρ σ is negative, it can be interprete as the smallest amount of work that must be supplie to guarantee the transition takes place. If it is positive, it is the largest amount of work we are guarantee to extract in the process. As the work system is both initially an finally in a pure state, no entropy is containe within it an its energy change must be completely ue to work being exchange with the system. Given the energy-conservation law that Thermal Operations follow equivalent to the first law, this iea of work automatically yiels a efinition of what heat is. In a given operation, the change in energy of work bit, system an heat bath must be zero, an hence we can straightforwarly ientify heat as the change in energy of the heat bath, or minus the change in energy on system an work bit. FIG. 7. We show the thermo-majorization curves of a state to which a work qubit in one of two pure states has been tensore. Aing this work system takes Z Z + e βw, extening the x-axis. When we tensor with the groun state to form ρ 0 0, the curve is the same as for ρ alone, but when the excite state is tensore, there is a change in the energy levels of the β-orering, an as a result the curve of ρ is compresse by a ratio of e βw. As we illustrate in Figure 7, the effect of appening a pure state of work to ρ is equivalent to stretching the thermo-majorization curve by a factor of e βw, an tensoring by the corresponing groun state to σ oes not change the curve [20]. In both cases the β-orer is preserve, an the new curves will have a lengthene x-axis [ 0, Z + e βw ]. These ifferent stretchings can serve to place the curve of ρ just above that of σ, in which case W will be the work of transition, in a similar way to the case of nonuniformity within Noisy Operations. 2. Monotones uner Thermal Operations, an the work of transition In Thermal Operations, the horizontal istance between a state s thermo-majorization curve an the y-axis is again a monotone for each value of y [0, ]. We enote these by L y an, as before, they never ecrease uner Thermal Operations. In particular, for blockiagonal ρ, we have: L yk ρ = L ρ = k rankρ e βeρ i, for y k = e βeρ i, k where all sums have been properly β-orere. Similarly to Lemma 2 we have: η ρ i, k < rank ρ, 59 Lemma 6. Given two states ρ an σ, where σ is blockiagonal in the energy eigenbasis, uner Thermal Oper-

12 ations: e βwρ σ L y ρ = max y Dσ L y σ Now, if we appen an iealize weight with Hamiltonian H W = w w w w as a work storage system initially R in the state 0, by efinition there exists a ifferent set of thermal operations that extracts work W σ τ from σ The proof is near ientical to that given in Lemma 2 for Noisy Operations an so we omit it here. If σ is not block-iagonal, the right han sie of Eq. 60 lower bouns e βwρ σ. To see this, recall that ecohering commutes with Thermal Operations, an hence if the transition ρ 0 0 σ W ρ σ W ρ σ is possible, so is ρ 0 0 σ D W ρ σ W ρ σ, an hence W ρ σ W ρ σd. 3. Bouns on the transition probability We can prove a result analogous to Eq. 33 for the thermal case: Lemma 7. Given two states ρ an σ, where σ is blockiagonal in the energy eigenbasis, uner Thermal Operations: e βwρ σ p e βwσ ρ D, 6 where as p, we assume W ρ σ 0. If W ρ σ 0, p = an the transformation from ρ to σ can be one eterministically, potentially extracting a finite amount of work. Proof. The previous Lemma 3 can be seen as a general statement about pairs of concave Lorenz-like curves: the minimum vertical ratio is lower an upper boune by the minimum an maximum horizontal ratios of the two. Given our previous efinitions of the work of transition, an the fact that p is the minimum vertical ratio of the two Lorenz curves as shown in Theorem 5, the result follows. The upper boun of Lemma 7 can be relate to the Jarzynski equality, which is foun to hol for general thermal operations applie to the system in an initial thermal state, see [36] for further etails. The equality states that, for a given thermal operation that extracts work w with some probability p w, we have that: e βw = w e βw pw =. 62 The Jarzynski equation is vali if the initial state is thermal, so let us take the special case of Lemma 7, of a process where we start with a thermal state τ an probabilistically go to some σ iagonal in energy, with optimal probability p. Because τ is the fixe point, the effect of that operation is trivial: τ T O ρ = τ = p σ + p X. 63 σ 0 0 T O τ W σ τ W σ τ. 64 By linearity, applying this set of TO to τ = p σ + p X yiels: p τ W σ τ W σ τ + p X sw, 65 where X SW is some joint system-weight state, with the weight in some work istribution p X w. Note that this operation is applie on both system an weight, an oes not nee to conserve the thermal state of the system alone. The Jarzynski equality for this operation reas: p e βwσ τ + p p X we βw =. 66 w The secon term in this sum is positive, an hence we have: p e βwσ τ, 67 which is the upper boun of Lemma 7. Note that in situations where the upper boun is saturate such as reversible processes with W σ τ = W τ σ, when the thermomajorization curve of σ is also a straight line the operation in Eq. 64 costs a ivergent amount of work in the case of failure i.e. from the state X in Eq. 63. D. Changing Hamiltonian Our results so far have assume that ρ an σ are associate with the same Hamiltonian. Suppose the initial system has Hamiltonian H an the final system Hamiltonian H 2. Following [20], this scenario can be mappe to one with ientical initial an final Hamiltonian, H, if we instea consier the transition between ρ 0 0 an σ where: H = H H Note that the partition function associate with H is Z = Z + Z 2. The height of the thermo-majorization curve of ρ 0 0 with respect to H, is ientical to that of ρ with respect to H on [0, Z ] an equal to on [Z, Z]. Similarly, the height of the thermo-majorization curve of σ is ientical to that of σ on [0, Z 2 ] an equal to on [Z 2, Z]. Hence by extening the efinition of Ṽx ρ so that Ṽx ρ = for x Z, we can reaily apply Theorems 4 an 5 to the case of changing Hamiltonians. Note that as L y ρ = L y ρ 0 0 for 0 y an similarly for σ, changing Hamiltonians oes not affect the results of Section III C.

13 3 IV. CONCLUSION Here, we have introuce a finite set of functions which, like the free energy, can only go own in the resource theory of Thermal Operations. We use these to compute the work of transition, an the maximum probability of making a transition between two states. Finally, we saw that the work of transition between the two states, an vice-versa, can be use to boun the maximum probability of making the transition. In maximizing the value of p in Eq. 2 to obtain p, we have attempte to maximize the fraction of σ present in a state obtainable from ρ. With access to a single two outcome measurement, σ can also be obtaine from ρ with probability at least p. There are other measures that one coul quantify in attempting to obtain a state that behaves like σ. For example, one coul consier the fielity between σ an a state reachable from ρ: F TO ρ, σ max ρ { F ρ, σ : ρ TO ρ }, 69 ] where F ρ, σ = tr [ ρσ ρ is the fielity between the two states. Investigating this problem is an open question, but note that for iagonal σ we have F TO ρ, σ F ρ, σ p. Another alternative woul be to consier herale probabilistic transformations. Here a 2-level flag system with trivial Hamiltonian an starting in the state 0 is provie with the initial state ρ. The goal is to transform both system an flag so that a measurement on the final flag state woul reveal that the system is in state σ with probability p an some other state with probability p. More concretely, one woul be intereste in maximizing the value of p in the transformation: ρ 0 0 T O ˆρ = pσ p X. 70 Due to the results in Section III B it is clear that the maximum value of p achievable in the herale case Eq. 70, is at least at large as p in the unherale case Eq. 2 for block-iagonal σ. In follow-up work to the initial version of this manuscript the converse was proven [58] an thus the two maximum probabilities are equal. In Appenix B we exten this analysis to consier the achievable herale probability when σ contains coherences or when one may use a catalyst to assist in the transformation. At the moment, although our results regaring maximum extractable work are general, little is known about transitions when the final state is not block-iagonal in the energy eigenbasis. In such a situation, our results provie necessary conitions but are not sufficient. Fining sufficient conitions is expecte to be ifficult, as we o not know such conitions even for non-probabilistic transformations. For recent results on the role of coherences in quantum thermoynamics, see for example [26 28, 32]. Nonetheless, we are able to utilize some of these results to provie bouns on the achievable herale probability when the target state contains coherences in energy. This is one in Appenix B. Our analysis has focuse on Noisy an Thermal Operations in the absence of a catalyst, i.e. an ancilla which is use to ai in a transition but returne in the same state. In Catalytic Thermal Operations, CTO, given ρ an σ, we are intereste in whether there exists a state ω such that: ρ ω TO σ ω. 7 If such an ω exists, we say ρ CT O σ. There exist instances where ρ TO σ an yet ρ CTO σ. Investigating when such catalytic transitions exists has le to a family of secon laws of thermoynamics that apply in the single-shot regime [25]. Having access to catalysts has the potential to achieve higher values of p than that efine by p an it woul be interesting to fin an expression or boun for the maximum value of p in the process: ρ CT O ρ = pσ + p X. 72 Note that a boun can be obtaine from any nonincreasing monotone of CTO, M say, that satisfies M pσ + p X pm σ. Bouning the maximum transition probability uner Catalytic Thermal Operations is mae more ifficult by the fact that the generalize free energies foun in [25] are not concave. However, for the case of herale probability, the situation is somewhat easier an in Appenix B, we completely characterize what is achievable uner CTO when the target state is block-iagonal in energy. Another avenue of research is to generalize our result to the case where one is intereste in not only maximizing the probability of obtaining a single state, but rather, fining the probability simplex of going to an ensemble of many states. Again, the fact that the monotones use in thermoynamics are not in general concave, means that straight application of the techniques use in entanglement theory [59] cannot be immeiately applie. Finally, by supplying more work or emaning that extra work is extracte, the value of p achieve can be raise or lowere. For W 0, one coul calculate p as a function of W for the states ρ W W an σ 0 0. For W > 0 the states to consier woul be ρ 0 0 an σ W W. What is the traeoff between p an W? As an example, the solution for qubit systems in the Noisy Operations framework is given in Appenix C. This work has focuse on the probability with which a given state can fluctuate into another uner a thermoynamical process. The term fluctuation is usually applie within thermoynamics to the concept of fluctuating work, a notion most famously capture by the Jarzynski equality an Crooks theorem. These were erive uner the framework of stochastic thermoynamics while our research was base on applying ieas from quantum information theory. Fining common groun between the two paraigms is likely to be beneficial to