arxiv: v1 [quant-ph] 7 Jan 2019

Transcription

1 Quantum Coherence in a Quantum Heat Engine arxiv: v1 [quant-ph] 7 Jan 2019 Yun-Hao hi, 1, 2, 3 Hai-Long hi, 4, 5 Xiao-Hui Wang, 2, 6, Ming-Liang Hu, 7 i-yuan Liu, 1, 3, 6 Wen-Li Yang, 1, 2, 6 and Heng Fan 3, 1 Institute of Modern Physics, Northwest University, Xi an , China 2 chool of Physics, Northwest University, Xi an , China 3 Institute of Physics, Chinese Academy of ciences, Beijing , China 4 tate Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of ciences, Wuhan , China 5 University of Chinese Academy of ciences, Beijing , China 6 haanxi Key Laboratory for Theoretical Physics Frontiers, Xi an , China 7 chool of cience, Xi an University of Posts and Telecommunications, Xi an , China (Dated: January 8, 2019) We identify that quantum coherence are valuable resource in the quantum heat engine, which is designed in a quantum thermodynamic cycle assisted by a quantum Maxwell s demon. This demon is in a superposed state. The quantum work and heat are redefined as the sum of coherent and incoherent parts in the energy representation. The efficiency of the heat engine potentially can exceed the Carnot efficiency due to the coherence consumption of the demon. In addition, we also discuss a heat engine driven by quantum coherence. This resource-driven viewpoint provides a direct and effective way to clarify the thermodynamic processes where the coherent superposition of states cannot be ignored. Introduction. Quantum information theory [1, 2] can play a crucial role in thermodynamics [3 12]. pecifically, the rigorous frameworks for the quantification of quantum correlations, such as entanglement [13], discord [14] and coherence [15], bring new insights into our understandings about quantum effects in thermodynamics. Those new concepts created from quantum information motivate us to reconsider the physical origin of thermodynamic notions such as work, heat, efficiency of heat engine etc. On the other hand, remarkable progress of quantum thermodynamics began with the fundamental relations in non-equilibrium statistical mechanics known as the Jarzynski equality [19, 20] or fluctuation theorems [21 24], which lead to the second law of thermodynamics as a corollary. Then the Maxwell s demon [28], as an intermediary, bridges information and thermodynamics [30], such as Landauer s principle [31, 32] and generalized second laws [24 27]. All of these support the view that information is a sort of physical existence, contributing to the adaptation of thermodynamics to quantum physics. Existing approaches to quantum thermodynamics are mainly based on operator formalisms [6, 21], which provides a general framework for the thermodynamics of information. Another approach to quantum thermodynamics is path integral formalism [34] based on two-point measurement scheme [20, 35]. This approach provides an effective way to calculate quantum work by utilizing various path integral techniques. However, it is difficult to figure out where information plays an important role and what kind of information is useful in a specific case, since both classical information and quantum information facilitate thermodynamic processes. It is interesting to explore such problems [3, 33]. In this paper, we will give the explicit expressions for some thermodynamic notions by virtue of the resource theory of quantum coherence [15 18]. As is known, coherence marks the departure of quantum mechanics from the classical realm. It is a very direct idea to demystify exactly those thermodynamic notions of quantum system via quantifying coherence. Then, we consider the quantum Maxwell s demon (QMD) [44 46] with off-diagonal elements (see Fig. 1) as a coherent carrier to drive a quantum zilard heat engine [38, 39] (schematic diagram in Fig. 2). The change of internal energy, work and heat after a thermodynamic cycle can be obtained. Then we utilize the resource theory of quantum coherence to interpret why the efficiency of heat engine tend to go beyond the Carnot efficiency, and this interpretation is clear enough so that we can directly obtain the quantum corrections to the classical expression if there is no coherence being consumed. Quantum Coherence. Let us begin with a brief review of quantum coherence. The rigorous framework for quantifying quantum coherence was first introduced by Baumgratz et al. who proposed that the conditions of the coherence measure are split into four [15]. One of the significant measures of quantum coherence based on these conditions is the relative entropy of coherence, which is defined as [15 17] C r (ρ) := min σ I (ρ σ) = (ρ diag) (ρ), (1) where I denotes the set of incoherent states, (ρ) = Tr(ρ ln ρ) is the von Neumann entropy [1] and ρ diag is the state obtained from ρ by deleting all the off diagonal elements. In the reference basis i, ρ diag can be expressed as i i ρ i i i. This entropic measure of coherence has a clear physical interpretation and has many

2 2 applications [17]. For instance, quantum coherence can enhance the success probability in Grover search quantum algorithm [40]. In Ref. [41], physical situations are considered where the resource theories of coherence and thermodynamics play competing roles, especially the creation of quantum coherence comes at a cost of energy. Therefore, coherence can be viewed as a potential quantum resource in quantum physical processes. Probing more applications of quantum coherence is not just about basic interests. Ideas. In general, we consider a quantum system described as H ψ n = E n ψ n, where H, E n and ψ n are the Hamiltonian, n-th eigenenergy and eigenstate, respectively. In the energy representation, a general density matrix of this system can be given by ρ = n P n ψ n ψ n + n m P nm ψ n ψ m, (2) then the internal energy of the system can be expressed as E = Tr[ρH] = n P n E n (3) where P n is the probability distribution of n-th energy eigenstate. In equilibrium P n obeys the canonical distribution. Note that although the density matrix contains off-diagonal part, the internal energy of the system is connected with the diagonal part due to the energy representation. In the following discussion, we will always choose the energy representation (basis). From the derivative of E, one obtains d E = n (E ndp n + P n de n ). Here, the previous viewpoint is analogizing it to the classical thermodynamic first law, i.e., dw c = de c + dq c = de c + T d c, where E c, W c, Q c and c are the classical internal energy, work, heat and entropy, respectively. Then, the quantum work can be identified as d W = n P nde n [36 38], and quantum heat is d Q = n E ndp n associated with T d c since the classical entropy c is defined as c = k B n P n ln P n. Nevertheless, one must have noticed that the classical entropy is actually the von Neumann entropy of diagonal density matrix in energy representation, namely c = k B (ρ diag ), (4) where k B is the Boltzmanns constant and ρ diag = n ψ n ρ ψ n ψ n ψ n denotes the diagonal part of ρ. Therefore, with the analogy above, this quantum system must have lost part of the information concerning its nonequilibrium (off-diagonal) state. To describe the whole quantum state, the entropy must be the whole k B (ρ) instead of k B (ρ diag ) which contains only the information of equilibrium (diagonal). Hence, combined with Eq. (1), it is natural to define the quantum heat as Q = E n dp n + k B T C r, (5) ρ σ n (a) Figure 1. (color online) Classical Maxwell s demon (CMD) vs quantum Maxwell s demon (QMD). (a) CMD can be considered as a classical mixture of two definite eigenstates g and e, while QMD with coherence should be preferred to a qubit described in the Bloch sphere, as shown in (b). Here, p g = (1 + z)/2, p e = (1 z)/2, F = (x iy)/2 and F = (x + iy)/2 are the classical probability distribution of ground state and excited state, and its off-diagonal elements, respectively. Note that x 2 + y 2 + z 2 1, thus F 2 p gp e, where the equality is achieved if it is a pure quantum state. where C r = C r (ρ) C r (σ) is the coherence consumption of the system from the initial state ρ to the final state σ. Then the first term of the right-hand side of Eq. (5) only determined by the diagonal (incoherent) part P n can be considered as the incoherent heat, i.e., Q incoh = ρ σ n E ndp n = T c = k B T (ρ diag ), and the second term can be understood as the coherent heat, i.e., Q coh = k B T C r. Thus, (b) Q = Q incoh + Q coh, (6) which gives the physical origin of quantum heat: the change of diagonal (incoherent) distribution and the coherence consumption. By virtue of the definition of quantum heat in Eq. (5), we can define the quantum work as W = P n de n + k B T C r, (7) ρ σ n since the first law E = W + Q. Likewise, the quantum work can be written as W = W incoh + W coh, (8) with the incoherent work W incoh = ρ σ n P nde n and the coherent work W coh = k B T C r. Thus, the quantum work can be viewed as the contribution of both the change of energy level under the invariable diagonal (incoherent) distribution and the coherence consumption. Quantum Maxwell s Demon assisted Quantum Thermodynamic Cycle. Now return to a basic problem of statistical mechanics, namely the Maxwell s demon [28] whose role of the information was firstly studied by zilard [29]. A classical Maxwell s demon (CMD) can be

3 3 Figure 2. chematic diagram of the quantum thermodynamic cycle assisted by quantum Maxwell s demon (QMD). (i) A system is prepared in a closed box with size L while QMD with coherence is isolated from the box. (ii iii) A wall, depicted by a vertical blue bar, is isothermally inserted at location l. And the measurement is performed later to register the information of the system into QMD s memory. (iii iv) According to the state of the quantum memory, the wall controlled by QMD undergoes a reversible and isothermal expansion until it reaches equilibrium. (v) To complete this thermodynamic cycle, the wall is removed isothermally. At the moment, the state of QMD is overwhelmingly different from the initial state and its coherence has been consumed. described as an incoherent state, whereas a quantum Maxwell s demon (QMD) with off-diagonal elements is in a coherent state, namely ρ D = p g g g + p e e e + F g e + F e g, (9) where p g and p e = 1 p g are the probability distributions in the ground state g and excited state e, respectively. It can be arbitrary one qubit system in the Bloch sphere (see Fig. 1). Here, we introduce the off-diagonal elements F and F to embody the coherence of QMD. The QMD and its related equalities (inequalities) have already been experimentally verified by using various systems in recent years [45 51]. As a carrier of coherence, QMD is typically designed to drive a zilard engine, which consists of a closed box with definite size L and one system (such as a single molecule). The whole thermodynamic cycle is briefly shown in Fig. 2, which is split into five stages: (i) Initial state, (ii) Insertion, (iii) Measurement, (iv) Expansion, and (v) Removal. At the stage of insertion, a wall is isothermally inserted at location l. We denote P L as the quantum probability of the system on the left after insertion, then P R = 1 P L is the probability of on the right. The things that QMD should do are performing a global measurement both on the system and itself (the controlled-not operation [42 44]), and then controlling the expansion of. When the wall is removed and this thermodynamic cycle is completed, returns to its initial state of equilibrium. In the presence of quantum coherence, it is insufficient to know its diagonal part. Instead, we describe the evolution of the whole system in terms of the full density matrix. In upplemental Material, we analyze the whole thermodynamic cycle and demonstrate that the quantum work done by the total system to the outside is given by W tot = T ( c + k B C r ) P R (p g p e ), (10) where T represents the temperature of the system, the separate sign denotes the energy difference between the two levels of QMD, the change of classical entropy is c, and C r = C r (ρ i D ) C r(ρ f D ) is the coherence consumption from the initial state ρ i D to the final ρf D of QMD. Here, the first term of Eq. (10) is the quantum heat absorbed from the outside, i.e., Q tot = T ( c + k B C r ), and the second term is the total change of the internal energy, i.e., E tot = P R (p g p e ). The explicit expressions for c and C r are given in the upplemental Material, but for our present purposes the crucial point is that the total work can be increased just because QMD consumes its coherence in the cycle. As we mentioned in ideas, the internal energy has nothing to do with quantum coherence in energy representation, while the work involves a coherent part that can be described by the relative entropy of coherence mathematically. Dividing the total work by the heat, one can obtain the efficiency of the quantum thermodynamic cycle, namely η = 1 P R(p g p e ) T ( c + k B C r ). (11) Note that the coherence consumption appears in the denominator of the minuend, which implies the efficiency can be improved due to the coherence consumption of QMD. Indeed, we strictly prove that the upper bound of the efficiency is greater than the Carnot efficiency,

4 F = p g p e F = 0.7 p g p e F = 0 Carnot efficiency replace the two memories with one memory M in which quantum coherence exists only and obtain the work extractable from the engine (see upplemental Material): W ext F + T c + k B T C r, (12) η P R Figure 3. (color online) The efficiency η with respect to the quantum probability P R of finding the system on the right side for different off-diagonal (coherent) element F. In the classical case, the efficiency must be limited by Carnot efficiency (the red line). However, with the coherence-assisted QMD, this limit may be allowed for a breakthrough (the case of pure sate with F = p gp e is depicted by black line and the blue line represents one of the mixed states with F = 0.7 p gp e). Here we choose T = 1, T D = 0.5, = 0.5, and F = F. i.e., η Carnot = 1 T D /T, where T D = k B ln(p g/p e) is the temperature of QMD (see upplemental Material and Fig. 3). In particular, there exists a critical probability PR cri, below which the efficiency goes beyond the Carnot efficiency. According to Eq. (11), PR cri satisfies T D ( c +k B C r ) = PR cri(p g p e ), which means QMD s absorption of heat maintains balance with the energy flow of the total composite system. Furthermore, when P R is larger than a certain value PR 0, the efficiency is zero. However, with the increase of coherence, the value of PR 0 will become higher. This quirk is completely the quantum effect as it depends on how much resource of quantum coherence we expend, in the light of the quantitative resource theory of coherence. Heat Engine Driven by Quantum Coherence. Before concluding this paper, we discuss a heat engine driven by quantum coherence, which is inspired by Ref. [53]. The authors introduced an information heat engine (IHE) with a system and a reservoir R, which is controlled by a demon consisting of two memories A and B, then one can extract work from this IHE by using quantum mutual information (or split into the classical correlation and quantum discord one of the nonclassical resources) between these two memories (the case of multiple reservoirs is studied in Ref. [52]). A natural question is that whether work can be extracted from the inherent property of one quantum system coherence, rather than the quantum correlation between two or more systems. We in which F = F (f) F (i) is the difference of free energy of the system, c = k B diag = [ρ (f) Mdiag ] [ρ (i) Mdiag ] is the diagonal (classical) entropy change, C r = C r [ρ (i) M ] C r[ρ (f) M ] is the coherence consumption of the memory. From Eq. (12), one can notice that if the contribution of F is ignored, the bound of the work extractable is given by the total entropy change of the memory, which is split into the change of incoherent part and coherence consumption. Even if no classical entropy changes, i.e. diag = 0, we still find the maximum work extractable is given by coherence consumption. And this is also in support of the resource-driven viewpoint, specifically, one can extract work from quantum coherence, a potential quantum physical resource, to drive a heat engine. ummary. We propose and study quantum thermodynamics by utilizing the resource theory of coherence based on its mathematically rigorous fashion. With the analogy in definition for related concepts, the subtle connection between quantum coherence and thermodynamic notions (work and heat) is probed quantitatively. Then a quantum thermodynamic cycle assisted by a quantum Maxwell s demon with off-diagonal elements is analyzed, where coherence consumption of the demon enhances the quantum work and the efficiency. By virtue of the resource-driven viewpoint, we explain that the efficiency of heat engine tend to go beyond the Carnot efficiency due to the decoherence of the demon. In addition, we also discuss a heat engine driven by quantum coherence based on the information heat engine. The maximum work extractable from which is given by coherence consumption, no doubt to the extraction of work via vanishing quantum coherence provides a strong support [56]. Our results are relevant for the quantum information processing in thermodynamic systems where the coherent superposition of states cannot be ignored. Yun-Hao hi thanks Z. C. Tu for insightful discussions. This work was supported by NFC (Grants Nos and ), the Key Innovative Research Team of Quantum Many-body theory and Quantum Control in haanxi Province (Grant No. 2017KCT-12) and the Major Basic Research Program of Natural cience of haanxi Province (Grant No. 2017ZDJC-32). Hu was supported by NFC (Grant No ), New tar Project of cience and Technology of haanxi Province (Grant No. 2016KJXX-27) and New tar Team of XUPT.

5 5 [1] M. A. Nielsen, and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, 2000). [2] V. Vedral, Rev. Mod. Phys. 74, 1 (2002). [3] J.M.R. Parrondo, J. M. Horowitz and T. agawa, Nat. Phys. 11, (2015). [4] J. Goold, M. Huber, A. Riera, L. del Rio, and P. krzypczyk, J. Phys. A 49, (2016). [5] F. Binder,. Vinjanampathy, K. Modi, and J. Goold, Phys. Rev. E 91, (2015). [6] P. trasberg, G. challer, and T. Brandes, Phys. Rev. X 7, (2017). [7] M. Horodecki, J. Oppenheim, Nat. Commun. 4, 2059 (2013). [8] F. Brandaõ, M. Horodecki, N. Ng, J. Oppenheim, and. Wehner, Proc. Natl Acad. ci. 112, (2015). [9] P. Ćwikliński, M. tudziński, M. Horodecki, and J. Oppenheim Phys. Rev. Lett 115, (2015). [10] Á. M. Alhambra,. Wehner, M. M. Wilde, and M. P. Woods Phys. Rev. A 97, (2018). [11] G. Francica, J. Goold, and F. Plastina, arxiv: [12] J. P. antos, L. C. Céleri, G. T. Landi, and M. Paternostro, arxiv: [13] V. Vedral, M. B. Plenio, M. A. Rippin, and P. L. Knight, Phys. Rev. Lett 78, 2275 (1997). [14] H. Ollivier and W. H. Zurek, Phys. Rev. Lett 88, (2002). [15] T. Baumgratz, M. Cramer, and M. B. Plenio, Phys. Rev. Lett 113, (2014). [16] A. treltsov, G. Adesso, and M. B. Plenio, Rev. Mod. Phys. 89, (2017). [17] M.-L. Hu, X. Hu, J. Wang, Y. Peng, Y.-R. Zhang, and H. Fan, Phys. Rep ,1 (2018). [18] A. Winter, D. Yang, Phys. Rev. Lett 116, (2016). [19] C. Jarzynski, Phys. Rev. E 78, 14 (1997). [20] H. Tasaki, arxiv: cond-mat/ [21] M. Campisi, P. Hänggi, and P. Talkner, Rev. Mod. Phys. 83, 771 (2011). [22] P. P. Potts and P. amuelsson, Phys. Rev. Lett 121, (2018). [23] E. Iyoda, K. Kaneko, and T. agawa, Phys. Rev. Lett 119, (2017). [24] J. J. Park,. W. Kim, and V. Vedral, arxiv: [25] T. agawa, Thermodynamics of Information Processing in mall ystems (pringer, 2013). [26] D. Mandal, H. T. Quan, and C. Jarzynski, Phys. Rev. Lett 111, (2013). [27] L. Mancino, V. Cavina, A. D. Pasquale, M. broscia, R. I. Booth, Phys. Rev. Lett 121, (2018). [28] J. C. Maxwell, Theory of Heat (Longmans, London, 1871). [29] L. zilard, Z. Phys. 53, 840 (1929). [30] K. Maruyama, F. Nori, and V. Vedral, Rev. Mod. Phys. 81, 1 (2009). [31] R. Landauer, IBM J. Res. Dev. 5, 183 (1961). [32] A. Bérut, A. Arakelyan, A. Petrosyan,. Ciliberto, R. Dillenschneider, and E. Lutz, Nature (London) 483, 187 (2012). [33] V. Vedral, arxiv: [34] K. Funo, and H. T. Quan, Phys. Rev. Lett 121, (2018). [35] J. Kurchan, arxiv:cond-mat/ [36] T. D. Kieu, Phys. Rev. Lett 93, (2004). [37] M. Esposito and. Mukamel, Phys. Rev. E 73, (2006). [38]. W. Kim, T. agawa,. D. Liberato, and M. Ueda, Phys. Rev. Lett 106, (2011). [39] M. Plesch, O. Dahlsten, J. Goold, and V. Vedral, Phys. Rev. Lett 111, (2013). [40] H.-L. hi,.-y. Liu, X.-H. Wang, W.-L. Yang, Z.-Y Yang, and H. Fan, Phys. Rev. A 95, (2017). [41] A. Misra, U. ingh,. Bhattacharya, and A. K. Pati, Phys. Rev. A 93, (2016). [42] W. H. Zurek, Phys. Rev. A 67, (2003). [43] H. T. Quan, Y. D. Wang, Y.-X. Liu, and C. P. un, Phys. Rev. Lett 97, (2006). [44] H. Dong, D. Z. Xu, C. Y. Cai, and C. P. un, Phys. Rev. E 83, (2011). [45] M. Naghiloo, J. J. Alonso, A. Romito, E. Lutz and K. W. Murch, Phys. Rev. Lett 121, (2018). [46] C. Elouard, D. Herrera-Mart a, B. Huard, and A. Auffèves, Phys. Rev. Lett 118, (2017). [47] P. A. Camati, J. P.. Peterson, T. B. Batalhão, K. Micadei, A. M. ouza, R.. arthour, I.. Oliveira, and R. M. erra, Phys. Rev. Lett 117, (2016). [48] M. D. Vidrighin, O. Dahlsten, M. Barbieri, M.. Kim, V. Vedral, and I. A. Walmsley, Phys. Rev. Lett 116, (2016). [49] J. V. Koski, A. Kutvonen, I. M. Khaymovich, T. Ala- Nissila, and J. P. Pekola, Phys. Rev. Lett 115, (2015). [50] J. V. Koski, V. F. Maisi, T. agawa, and J. P. Pekola, Phys. Rev. Lett 113, (2014). [51]. Toyabe, T. agawa, M. Ueda, E. Muneyuki and M. ano, Nat. Phys. 6, (2010). [52] L.-H. Ren, and H. Fan, Phys. Rev. A 96, (2017). [53] J. J. Park, K.-H. Kim, T. agawa, and. W. Kim, Phys. Rev. Lett 111, (2013). [54] T. agawa, and M. Ueda, Phys. Rev. Lett 100, (2008). [55] T. agawa, and M. Ueda, Phys. Rev. Lett 102, (2009). [56] M. O. cully, M.. Zubairy, G.. Agarwal, H. Walther, cience 299, 862 (2003).

6 6 upplemental Material for: Quantum Coherence in a Quantum Heat Engine QUANTUM MAXWELL DEMON AITED QUANTUM THERMODYNAMIC CYCLE tage (i): Initial state Initially, the system, prepared in a closed box with size L, contacts with the reservoir at temperature T = (k B β) 1 so that it is in equilibrium. With the initial Hamiltonian H = n E n(l) ψ n (L) ψ n (L), the density matrix of the initial state of the total system can be written as a product state ρ i = ρ i ρ i D = n e βen(l) ψ n (L) ψ n (L) ρ i Z(L) D, (13) where Z(L) = n e βen(l), E n (L) and ψ n (L) are the partition function, n-th eigenenergy and eigenstate of the system in this box with size L. Remark 1 : The system will always be in the isothermal process with the reservoir, leading to its thermalization so that all the coherence among energy levels vanishes [38]. Even so, it is necessary to describe the dynamics of the total system in terms of the full density matrix, because our QMD with off-diagonal elements is still coherent at another temperature T D = (k B β D ) 1. Remark 2 : Here we emphasize that the box doesn t have to be an infinite potential well model, i.e., the eigenenergies and eigenfunctions of the system are E n (L) = 2 π 2 n 2 /(2mL 2 ) and x ψ n (L) = 2/L sin(πnx/l), where the quantum number n ranges from 1 to. But, it is convenient to utilize such a simple model to further expound the results by numerical methods. tage (ii): Insertion After preparing initial state, a wall is then isothermally inserted at a certain position l. In the classical situation, the position of the system is definite, so either it is on the left side or stay on the right side, while in the quantum case the system is simultaneously on both sides of the wall, until its location is determined by the measurement of the observer. Therefore, when the wall insertion is completed, the total density matrix can be expressed as a mixture, namely ρ ins = [P L ρ L (l) + P R ρ R (L l)] ρ i D, (14) where ρ L (l) = n e βen(l) Z(l) ψn L (l) ψn L (l) and ρ R (L l) = n e βen(l l) Z(L l) ψn R (L l) ψn R (L l) are the independent densities of on the left and right sides, P L = Z(l)/[Z(l) + Z(L l)] and P R = 1 P L are the corresponding quantum probabilities, which differ from the classical probabilities PL classical = l/l and PR classical = (L l)/l [44]. Remark: In this stage, QMD is still out of work, which is the same as CMD. And because the system has no coherence, thus the results (work, heat, etc.) only relate with diagonal part of, as in the case of no coherence in Ref. [38]. However, it will make a real difference in the following stages. tage (iii): Measurement Without impacting the original performance of the demon in ZE [29], the measurement is done by using the controlled-not operation [43, 44], namely U = n [ ψ L n (l) ψ L n (l) ( g g + e e ) + ψ R n (L l) ψ R n (L l) ( e g + g e )]. (15) It is easy to proof that U is a unitary operator, which is different from the positive operator valued measure (POVM) in previous discussions [53 55]. Here, the physical process corresponding to U is that the demon will keep the original state if the system with n levels is on the left side, while the demon flips its state when is on the right. After this operation, the total density matrix is given by

7 7 ρ mea = Uρ ins U = P L ρ L (l) (p g g g + p e e e + F g e + F e g ) + P R ρ R (L l) (p e g g + p g e e + F g e + F e g ) = [p g P L ρ L (l) + p e P R ρ R (L l)] g g + [p e P L ρ L (l) + p g P R ρ R (L l)] e e + [F P L ρ L (l) + F P R ρ R (L l)] g e + [F P L ρ L (l) + F P R ρ R (L l)] e g, (16) which means QMD is now correlated with the system, contributing to the work to the outside. Remark 1 : Compared with CMD, QMD seem to be blurry-eyed, so that it can t exactly distinguish which side the system stays on. In the limits of F 0 and T D 0, ρ mea P L ρ L (l) g g + P R ρ R (L l) e e, which implies that QMD is cured, i.e., its ground state g corresponds to the system on the left and vice versa, as the case of CMD. Remark 2 : For registering the information of the system into QMD s memory, the outside agent have to do work on the total system in this process, leading to the change of the total internal energy. In Appendix, we calculate the mean value of this work, namely W mea = P R (p g p e ). Note that this work has nothing to do with the off-diagonal elements, because the flip of off-diagonal part makes no contribution to the change of internal energy. However, the off-diagonal part has affected the total density matrix, giving rise to a more complex form about entropy and heat as the following discussions. In this stage, the expansion of the system is slowly enough to ensure the process reversible and isothermal, which is controlled by QMD. The diagonal bases of QMD, i.e., g g and e e, represent the definite ground state and excited state, respectively. The evolution of the states of controlled by these two bases follow the rules that the wall at position l is allowed to move to the left side at final position l g if QMD is on g g and move to the right side at final position l e if QMD is on e e. Then, how about the coherent (offdiagonal) part? From Eq. (16), one can see that the off-diagonal bases of QMD, i.e., g e and e g, corresponding to the states [F P L ρ L (l) + F P R ρ R (L l)] and [F P L ρ L (l) + F P R ρ R (L l)], are unable to change the density distribution [P L ρ L (l) + P R ρ R (L l)], at the moment when the measurement is completed but the expansion is not performed yet, in spite of the coefficients F and F. And that implies the off-diagonal bases have no role in controlling these two states to make the wall move in a definite direction during expansion, but the post-expansion densities must be related to l g and l e. tage (iv): Expansion Mathematically, the total density matrix after expansion is expressed as ρ exp = [p g P L ρ L (l g ) + p e P R ρ R (L l g )] g g + [p e P L ρ L (l e ) + p g P R ρ R (L l e )] e e + [F P L ρ L l g,l e + F P R ρ R l g,l e ] g e + [F P L ρ L l e,l g + F P R ρ R l e,l g ] e g, (17) where ρ L l g,l e and ρ R l g,l e are the post-expansion densities controlled by g e, while ρ L l e,l g and ρ R l e,l g are the postexpansion densities controlled by e g. Remark: In classical physics, the system must be on one side of the wall with another side empty after insertion, thus the wall is doomed to be moved to an end boundary of the box due to expansion. In contrast, quantum system can be simultaneously on both sides of the wall, thus the equilibrium position is somewhere in the box rather than the boundary. The condition of equilibrium, in mechanics, is that the wall has equal and opposite forces on the two sides, i.e. F L eq = F R eq. tage (v): Removal For purpose of the thermodynamic cycle, the system must be reset to its own initial state, i.e. ρ i =

8 8 n e βen(l) Z(L) ψ n (L) ψ n (L). The corresponding physical process is the wall inserted on the stage (ii) will be removed and all the ensembles of the system will evolve into ρ i, namely where ρ rem D ρ rem = n e βen(l) ψ n (L) ψ n (L) ρ rem D. (18) Z(L) = [(p g P L + p e P R ) g g + (p e P L + p g P R ) e e + (F P L + F P R ) g e + (F P L + F P R ) e g ]. (19) Here we emphasize that the system has been reset to thermal equilibrium state with no coherence, but it doesn t mean to say this requirement will also apply to QMD. The state of QMD remains coherent, but obviously it s a lot less coherence than it was at the beginning. And that s probably a linkage to work and heat in the whole thermodynamic cycle, as discussed in the next section. Beyond Carnot efficiency via Quantum Maxwell s Demon with Coherence What we have talked above is the evolution of state of the total system. Because of the coherence, it is insufficient to know its diagonal part. Instead, we describe the evolution of our system in terms of the full density matrix, which is split into five, namely tage (i-v). As seen in the follows, this approach is in support of studying what we are more concerned, i.e., whether this quantum thermodynamics cycle assisted by QMD with coherence could make difference to the efficiency (compared with the the classic case, which is limited by Carnot cycle) or not. First of all, let us calculate the change of total internal energy E tot by taking ρ rem as the final density ρ f at the end of this thermodynamics cycle. Focusing on the initial and final states, we can obtain E tot = E f E i = Tr[(ρ f ρ i )(H + H D )] = P R (p g p e ), (20) where = E e E g is the gap of the two levels of QMD. Note that E tot = W mea, i.e., the change of total internal energy merely results from the work done by the outside agent during measurement, as what mentioned before. Remark: From Eq. (20), one can notice that the change of total internal energy E tot does not involve the offdiagonal part at all. In fact, it is the same to the internal energy E of an arbitrary quantum system with the Hamiltonian H and the density ρ, namely Another important thermodynamic quantity is the total heat Q tot absorbed from the outside, which is associated with the total entropy change due to the reversibility, namely Q tot = k B T [(ρ f ) (ρ i )] = k B T [( e βen(l) ψ n (L) ψ n (L) ρ f D Z(L) ) n ( e βen(l) ψ n (L) ψ n (L) ρ i Z(L) D)] n = k B T [(ρ f D ) (ρi D)] = T ( c + k B C r ), (21) where the second to the last line follows from the subadditivity equality for von Neumann entropy [1] and the last line follows from the definition of the coherence consumption C r := C r (ρ i ) C r (ρ f ) = C r (ρ i D ) C r(ρ f D ), and c := k B diag = k B [(ρ f diag ) (ρi diag )] = k B [(ρ f Ddiag ) (ρi Ddiag )]. According to what we discussed above, the initial density matrix and the final density matrix are expressed as ( ) ρ i pg F D = F (22) p e and ( ρ f D = pg P L + p e P R F P L + F P R F P L + F P R p g P R + p e P R ), (23) respectively. Hence, one can obtain the change of classical entropy c = k B [(ρ f Ddiag ) (ρi Ddiag )] = k B [p g ln p g + p e ln p e (p g P L + p e P R ) ln (p g P L + p e P R ) (p e P L + p g P R ) ln (p e P L + p g P R )] and the coherence consumption C r = C r (ρ i D) C r (ρ f D ) = (ρ f D ) (ρi D) c /k B. (24) (25) Note that in the classical process, the heat without coherence can be denoted as Q incoh = T diag, which just depends on the diagonal part. In general, we specify Q tot = Q incoh + Q coh (26) as the total heat absorbed from the outside, where the coherent item Q coh = k B T C r is proportional to the coherence depletion C r.

9 9 Then if we still believe that the first law holds, i.e., E tot = Q tot W tot, the total work W tot done by the system to the outside will be expressed as W tot = T ( c + k B C r ) P R (p g p e ). (27) An alternative derivation of this result can be constructed by using W tot = F with the standard free energy F (ρ) = Tr(Hρ) T (ρ) [3, 9], which is a special case of the general free energy F α (ρ) = F eq + T D α (ρ ρ eq ) based on quantum Renyi entropies D α (ρ ρ eq ) when α 1 [10 12]. In analogy to Eq. (26), the total work can be split into W tot = W incoh + W coh, (28) where the incoherent work is W incoh = Q incoh + E tot = T c P R (p g p e ) and W coh = Q coh = k B T C r represents the coherent part. Now turn to the crucial efficiency, which is defined as dividing the total work done by the system to the outside by total heat absorbed from the outside, namely η = W tot Q tot = 1 P R(p g p e ) T ( c + k B C r ). (29) Here, the efficiency η is a monotonous increasing function of the coherence depletion C r. The more coherence depleted, the higher efficiency we can obtain. If there is no coherence being consumed, i.e. C r = 0, it will be limited by Carnot efficiency { η max max 1 P } R(p g p e ) k B T c { } P R (p g p e ) = lim 1 P R 0 k B T [p g ln p g + p e ln p e (p g P L + p e P R ) ln (p g P L + p e P R ) (p e P L + p g P R ) ln (p e P L + p g P R )] { } P R (p g p e ) = lim 1 P R 0 k B T [p g ln p g + p e ln p e (p g P L ) ln (p g P L ) (p e P L ) ln (p e P L ) + P R (p g p e ) ln (p g /p e ))] = 1 k B T ln (p g /p e ) = 1 T D T, (30) where the third line follows from the Taylor expansion (c 1 + c 2 x) ln (c 1 + c 2 x) = c 1 ln c 1 + c 2 (ln c 1 + 1)x + o(x 2 ) with x 0 (here c 1 and c 2 are constants), and the last line follows from the definition of T D := /[k B ln (p g /p e )] due to the initial probability distribution p g = 1/(1 + e β D ) and p e = 1 p g. This result speaks for itself, i.e., one may go beyond the Carnot efficiency η Carnot = 1 T D /T by consuming a certain amount of quantum coherence, see Fig. 3. mutual information (or split into the classical correlation and quantum discord one of the nonclassical resources) between these two memories (the case of multiple reservoirs is studied in Ref. [52]). A natural question is that whether work can be extracted from the inherent property of one quantum system coherence, rather than the quantum correlation between two or more systems. We give the answer by replacing the two memories (as shown in Fig. 1 in Ref [53]) with one memory M in which quantum coherence exists only. HEAT ENGINE DRIVEN BY QUANTUM COHERENCE In this section, we will discuss a heat engine driven by quantum coherence, which is inspired by Ref. [53]. It has been demonstrated in Ref. [53] that if an information heat engine (IHE) with a reservoir is controlled by a demon consisting of two memories A and B, then one can extract work from this IHE by using quantum tage (i): Initial state The demon consisting of one memory M is off-line and the system contacts with the reservoir R at temperature T = (k B β) 1 in thermodynamic equilibrium, i.e. the total density matrix of the initial state is expressed

10 10 as ρ (i) = ρ (i) M ρ(i) R = ρ (i) (i) M e βh Z (i) e βh R Z R, (31) where Z (i) = Tr[e βh(i) ] and ZR = Tr[e βh R ] are the initial partition functions of the system and the reservoir R, respectively. In general, there is no restriction on ρ (i) M. tage (ii): Unitary evolution The system begins to interact with R (a unitary evolution), and M is still off-line. Then the density matrix is given by with U (1) = I M U (1) R. ρ (1) = U (1) ρ (i) U (1), (32) Note that the final state is not necessarily the canonical distribution, i.e. ρ (f)can R = exp ( βh (f) )/Z(f) exp ( βh R )/Z R. Maximum work extractable and general second law ince the POVMs increases the entropy while the unitary evolution ahead of the measurement keep it invariant, i.e., [ρ (2) ] [ρ (1) ] = [U (1) ρ (i) U (1) ] = [ρ (i) ], one can obtain [ρ (i) R ] + [ρ(i) M ] [ρ(2) M ] + k p k [ρ (2)k R ] (36) [ρ (2) M ] + [ρ(2) R ], where the second line follows from the subadditivity of von Neumann entropy [1]. Considering the final state of R, i.e., ρ (f) R = Tr M [ρ (f) ] = k U R k ρ(2)k R U k R, we can obtain the following inequality by virtue of the concavity of the von Neumann entropy [1] tage (iii): POVM [ρ (f) R ] k p k [ρ (2)k R ]. (37) This stage is where the demon consisting of M started working. The measurement is performed by M with POVMs (positive operator valued measures) [1, 25]. Explicitly, the measurement process is implemented by performing a unitary transformation U (2) on the whole system followed by a projection measurement {Π k M : k k } (a rank-1 projector) only on M, namley ρ (2) = k = k Π k M U (2) ρ (1) U (2) Π k M p k k M k ρ (2)k R, (33) where p k = Tr[Π k M U (2) ρ (1) U (2) Π k M ] is the measurement outcome registered by the memory and the post-measurement state of R is ρ (2)k R = Tr A [Π k M U (2) ρ (1) U (2) Π k M /p k]. tage (iv): Feedback control After measurement, the demon will control R according to the outcome p k. Mathematically, this feedback control is performed by a unitary operator U (3) = k then the final state becomes k M k U k R, (34) ρ (f) = U (3) ρ (2) U (3). (35) Likewise, we have ρ (f) M = Tr R[ρ (f) ] = k k M k = Tr R [ρ (2) ] = ρ (2) M, thus [ρ (f) M ] = [ρ(2) M ]. (38) By substituting Eq. (37) and Eq. (38) into Eq. (36), one obtains [ρ (i) R ] [ρ(f) R ] [ρ(f) M ] [ρ(i) M ], (39) which implies the entropy decrease of the heat engine and the reservoir cannot exceed the entropy increase of the memory, or should say, the whole system satisfies a general principle of entropy increase, i.e., R + 0, where R = [ρ (f) R ] [ρ(i) R [ρ (i) M ]. ] and = [ρ(f) M ] According to the quantum version of Klein s inequality, i.e. Tr[ρ (f) R ln ρ(f)can R ] Tr[ρ (f) R ln ρ(f) R ] [25], Eq. (39) becomes [ρ (i) R ] + Tr[ρ(f) R ln ρ(f)can R ]. (40) To investigate the work extractable from the heat engine, we turn the inequality above into the following form E (i) E(f) + E(i) R E(f) R with E (i) = Tr[ρ(i) E (i) R = Tr[ρ(i) R H(i) R F (i) H(i) ], E(f) ], E(f) R F (f) + k BT (41) = Tr[ρ(f) H(f) ], (i) ], F = Tr[ρ(f) R H(f) R =

11 11 k B T ln Z (i) (f), and F = k B T ln Z (f). Thus, by the definition of the work extractable W ext := E + Q, where E = E (f) E(i) is the change of the internal energy of and Q = E (f) R E(i) R is the heat exchange between and R, we obtain W ext F + T c + k B T C r, (42) with F = F (f) F (i) and = diag + C r, where diag = [ρ (f) Mdiag ] [ρ(i) Mdiag ] is the diagonal (classical and incoherent) entropy change, C r = C r [ρ (i) M ] C r [ρ (f) M ] is the coherence depletion of the memory. From Eq. (42), one can notice that if the contribution of F is ignored, the bound of the work extractable is given by the total entropy change of the memory, which is split into the change of incoherent part and coherence depletion. Even if no classical entropy changes, i.e. diag = 0, we still find the maximum work extractable is given by coherence depletion. And this gives the answer to the question above: one can extract work from quantum coherence, a potential quantum physical resource, to drive a heat engine.