Using thermodynamics to identify quantum subsystems

# Identifying quantum subsystems using Clausius' second law of thermodynamics

## Abstract

There are many ways to decompose the Hilbert space of a composite quantum system into tensor product subspaces. Different subsystem decompositions generally imply different interaction Hamiltonians , and therefore different expectation values for subsystem observables. This means that the uniqueness of physical predictions is not guaranteed, despite the uniqueness of the total Hamiltonian and the total Hilbert space . Here we use Clausius’ version of the second law of thermodynamics (CSL) and standard identifications of thermodynamic quantities to identify possible subsystem decompositions. It is shown that agreement with the CSL is obtained, whenever the total Hamiltonian and the subsystem-dependent interaction Hamiltonian commute (i.e. ). Not imposing this constraint can result in the transfer of heat from a cooler to a hotter subsystem, in conflict with thermodynamics. We also investigate the status of the CSL with respect to non-standard definitions of thermodynamic quantities and quantum subsystems.

## I Introduction

More than fifty years ago, it was noted in the context of quantum electrodynamics that there is an inherent relativity present within the definition of a quantum subsystem Power and Zienau (1959); Power and Thirunamachandran (1978). There are many ways of decomposing the total Hilbert space of an atom and its surrounding free radiation field into a tensor product of subspaces. Each decomposition implies a different form of the corresponding atom-field interaction Hamiltonian. The most known forms of this interaction Hamiltonian are the minimal-coupling, the multipolar and the so-called rotating-wave Hamiltonians Drummond (1987); Cohen-Tannoudji et al. (1989); Baxter et al. (1990); Woolley (2000); Stokes et al. (2012); Stokes (2013). Although all of these Hamiltonians have the same spectrum and are unitarily equivalent, they result in different predictions for the time evolution of subsystem expectation values. Overcoming this problem requires identifying physically viable subsystem decompositions, a subject that has now spread out to a much wider quantum physics community Nielsen and Chuang (2000); Zanardi et al. (2004); Barnum et al. (2004); DugiÄ and JekniÄ (2006).

Suppose we decompose the total Hilbert space of two interacting quantum systems and in two different ways,

 H=Ha⊗Hb=H′a⊗H′b. (1)

Moreover, we introduce the orthonormal bases of and the orthonormal bases of with respectively. Both sets of product states, and , form an orthonormal basis in , so they are related by a unitary transformation ;

 |i′a,j′b⟩=U|ia,jb⟩. (2)

Consequently, the operator space , consisting of the linear maps on , possesses bases and , such that

 O′μa⊗O′νb=UOμa⊗OνbU†. (3)

Subsystem relativity arises when the unitary transformation is not of the form . In this case the tensor product decompositions and define physically distinct subsystems and . Adopting the decomposition , observables associated with subsystem have the form

 Oa=∑μcμOμa⊗Ib, (4)

where denotes an identity operator. These operators do not act on but generally have a non-trivial effect in . On the other hand, adopting the decomposition , observables associated with subsystem have the form

 O′a=∑μc′μO′μa⊗I′b. (5)

These operators do not act on but generally have a non-trivial effect in .

In general, the total Hamiltonian of a composite quantum system can be partitioned into free and interaction components in many different ways. Each different partition corresponds to a different subsystem decomposition. For example

 H{Oμa,Oνb} :=Ha{Oμa}⊗Ib+Ia⊗Hb{Oνb}+V{Oμa⊗Oνb} =Ha{O′μa}⊗I′b+I′a⊗Hb{O′νb}+V′{O′μa⊗O′νb} =:H′{O′μa,O′νb} (6)

where , , , , and denote specific functional forms of basis operators. In terms of the distinct operator bases associated with different subsystem decompositions, the total Hamiltonian possesses different interaction components and has different functional forms. Moreover, we have

 ⟨Oa⟩t =⟨ia,jb|eiH{Oμa,Oνb}tOae−iH{Oμa,Oνb}t|ia,jb⟩ ≠⟨i′a,j′b|eiH′{O′μa,O′νb}tO′ae−iH′{O′μa,O′νb}t|i′a,j′b⟩ =⟨O′a⟩t, (7)

which means that the time evolution of subsystem observables depends in general on the chosen subsystem decomposition.

Naturally, the fundamentally important question occurs as to which of the above predictions correctly describes the physical subsystem . Ultimately this question must be answered by comparison with empirical findings. For example, using standard open quantum systems approaches, it has been shown that non-negligible concentrations of bare energy,

 H0 = Ha+Hb, (8)

can occur, when an atom couples to an electromagnetic field Kurcz et al. (2010a); Stokes et al. (2012). Rejecting this bare energy build-up as non-physical (it is not something observed empirically) restricts the form of the atom-field interaction to one which conserves the bare energy. Such a restriction may be especially nontrivial given that it may seem incompatible with the form of interaction obtained a fortiori from elementary first principles.

The laws of thermodynamics are some of the oldest and best-established laws of modern physics and are extensively supported by centuries of empirical evidence. Here we propose applying Clausius’ form of the second law of thermodynamics (CSL) in quantum physics to restrict the possible forms of quantum interactions and subsystem decompositions of composite quantum systems. The CSL Clausius (1854) can be stated as follows:

• No process is possible whose sole result is the transfer of heat from a cooler to a hotter body.

• Since the laws of thermodynamics are usually taken to refer to average quantities pertaining to large numbers of elementary systems, their status within the microscopic domain, which is generally believed to require a quantum-theoretic modelling, is less clear. In the quantum setting, studies of the thermodynamical laws have largely taken place within the context of open-quantum system’s theory, wherein small systems are coupled to thermal reservoirs with infinite degrees of freedom. General quantum versions of basic thermodynamic results have been established in this context Breuer and Petruccione (2007); Alicki and Lendi (2007); Kosloff (2013). Among these are the fluctuation theorem-type results which suggest that the statistical laws of thermodynamics may occasionally be transiently violated for individual systems Evans et al. (1993); Evans and Searles (1994). There are also experimental results which support this possibility Wang et al. (2002); Gieseler et al. (2014).

Situations consisting of only a small number of interacting quantum systems have received widespread attention only relatively recently. Most of the work on quantum thermodynamics focusses on identifying quantum versions of the thermodynamic quantities and their associated laws for a single elementary system in contact with a heat bath Esposito et al. (2010); Esposito and Broeck (2011); Horodecki and Oppenheim (2013); BrandÃ£o et al. (2015). Our focus here is on the Clausius form of the second law for interacting elementary systems initially in contact with separate heat baths so as to be prepared in thermal states. These are then isolated and allowed to evolve. Results relating to this situation are fewer and further between. Tasaki has shown that a general entropic version of the second law of thermodynamics can be proven in this situation Tasaki (2000). Weimer et al. have developed an energy-flux formalism based on the LEMBAS (local measurement basis) principle, which they use to give general definitions of heat and work within the interacting setting Weimer et al. (2008). This formalism is further studied and extended to include open quantum systems in Hossein-Nejad et al. (2015), where it is also shown to give rise to an entropic version of the second law. Esposito et al. have given general identifications of heat and work for a system in contact with several heat baths, and have also identified a useful partition of the change in the system’s von Neumann entropy into an always positive entropy production term and a reversible heat-flux term Esposito et al. (2010); Esposito and Broeck (2011).

In this paper we show that the result of Tasaki can be extended to include Clausius’ version of the second law of thermodynamics, but only provided one assumes that the bare energy is conserved. This condition can be expressed through the vanishing of any one of three equivalent commutators as

 [H0,V]≡[H,V]≡[H0,H]=0 (9)

where . Whether or not condition (9) holds depends on the particular interaction, i.e., subsystem decomposition, selected within a model. Assuming the standard identifications of heat and quantum subsystems are correct, then for the CSL to hold it is sufficient that condition (9) holds. Rejecting condition (9) means accepting at least transient violations of the CSL, or employing non-standard definitions of heat and/or quantum subsystems.

There are five sections in this paper. In section II we recap the standard quantum thermodynamics definitions of heat, work and entropy. We then determine the condition (9) as sufficient in order that the CSL holds for two interacting systems initially prepared in thermal states at different temperatures. We introduce a particularly simple interacting oscillator system which we use to investigate the CSL in the following sections. In section III we demonstrate violations of the CSL using a particular interaction Hamiltonian for the simple two-oscillator system introduced in section II. In section IV we extend our previous analyses by adopting alternative identifications of heat, work and entropy Weimer et al. (2008); Hossein-Nejad et al. (2015); Esposito et al. (2010). Finally we summarise our findings in section V.

## Ii Theoretical background

In this section we begin by reviewing the standard definitions of the basic thermodynamic quantities. We then define the heat transfer between interacting subsystems. Finally we determine a sufficient condition (9) in order that the CSL holds for the interacting system.

### ii.1 Energy, heat and entropy: the standard definitions

The first law of thermodynamics states that

 ΔH=ΔW+ΔQ (10)

where is a systems total energy, is the change in the external work energy and is the change in heat energy. In quantum theory the expected rate of change of the total energy is given by

 ⟨dHdt⟩=tr(˙Hρ)+tr(H˙ρ) (11)

where is the system’s density matrix. The first term on the right-hand side in (11) can be identified as the work contribution and the second term can be identified as the heat contribution Alicki (1979); Kosloff and Ratner (1984); Kieu (2004); Henrich et al. (2006, 2007). When the total energy is not explicitly time-dependent the work contribution vanishes so that . If we identify the entropy via the Gibbs relation

 ΔS=βΔQ (12)

where is the inverse temperature, then we see that if and only if . In this paper we adopt the above definitions of work, heat and entropy, and apply them on the level of two quantum subsystems and initially prepared in thermal states. The corresponding subsystem quantities are respectively denoted , and with . In the following, we consider the general case in which the subsystems are interacting, but in which the energies , and are not explicitly time-dependent.

For concreteness and simplicity we consider for the most part two interacting quantum harmonic oscillators (QHOs) labelled and , though our main results hold for arbitrary interacting quantum systems. We consider the situation in which the oscillators and are prepared at in the thermal states

 ρth.c=1Zce−βcHc (13)

with . Here is the inverse temperature of oscillator and the ’s are the oscillator energy operators defined by

 Ha=ωaa†a,Hb=ωbb†b. (14)

The operators and in this equation are bosonic annihilation operators satisfying and which act respectively, within the abstract seperable Hilbert spaces and of the corresponding oscillators. The eigenstates of form a basis of orthonormal states in the composite system’s Hilbert space . We denote these states , and refer to them as the bare (as oppossed to dressed) eigenstates. The energy eigenvalue corresponding to the state with is denoted . The functions are included in (13) to ensure the correct normalisation of the thermal states. Starting with the initial state of the composite system

 ρ=ρth.a⊗ρth.b (15)

enables one to form clear questions about heat transfer between the two systems at subsequent times.

### ii.2 Heat transfer according to the standard definitions

The identification of the heat changes and of the two subsystems within a time follows from the standard definitions given above in Section II.1. Since the Hamiltonian is time-independent there are no work contributions. So the heat energies of the subsystems are given by the bare energies and . In what follows the basic quantity we will consider is the heat transfer from oscillator into oscillator within a time . We denote this quantity and define it as follows

 ΔQa→b:=ΔQb−ΔQa, ΔQc=⟨Hc(t)−Hc(0)⟩ρ (16)

with . Here denotes the standard quantum expectation value taken in the state given in (15). In terms of the CSL can be expressed simply as

 sgnΔQa→b=sgn(βb−βa) (17)

which in words, states that the (heat) energy acquired by oscillator from oscillator , within a time , is positive whenever oscillator is initially hotter, and is negative whenever is initially hotter.

### ii.3 Conservation of bare energy and the second law

We will now show that the condition (9) is a sufficient condition in order that the CSL holds. The proof in this section builds upon the work of Tasaki Tasaki (2000), which uses the techniques of Jarzynski Jarzynski (1997). We begin by defining the free entropy change

 ΔS0 :=βa⟨Ha(t)−Ha(0)⟩ρ+βb⟨Hb(t)−Hb(0)⟩ρ =βaΔQa+βbΔQb (18)

where again the initial state is given in (15). Next we define the “classical” average of a given function as

 E[f]t := ∑nmpq e−βaωnae−βbωmbZaZb|Upq;nm(t)|2f(ωna,ωmb,ωpa,ωqb) (19)

where are the transition amplitudes between the bare states.

It is straightforward to verify that the four projection functions defined by , , and satisfy

 E[ωa]t = ∑nmpqe−βaωnae−βbωmbZaZbωna ×⟨na,mb|U−1(t)|pa,qb⟩⟨pa,qb|U(t)|na,mb⟩ = ∑nmωnae−βaωnae−βbωmbZaZb = tr(Ha(0)ρ)=⟨Ha(0)⟩ρ, E[ωb]t = ∑nmpqe−βaωnae−βbωmbZaZbωmb (20) ×⟨na,mb|U−1(t)|pa,qb⟩⟨pa,qb|U(t)|na,mb⟩ = ∑nmωmbe−βaωnae−βbωmbZaZb = tr(Hb(0)ρ)=⟨Hb(0)⟩ρ

where we have used the completeness of the bare states. Similarly

 E[ω′a]t=∑nme−βaωnae−βbωmbZaZb × ⟨na,mb|U−1(t)[∑pqωpa|pa,qb⟩⟨pa,qb|]U(t)|na,mb⟩ =tr(Ha(t)ρ)=⟨Ha(t)⟩ρ, E[ω′b]t=∑nme−βaωnae−βbωmbZaZb × ⟨na,mb|U−1(t)[∑pqωqb|pa,qb⟩⟨pa,qb|]U(t)|na,mb⟩ =tr(Hb(t)ρ)=⟨Hb(t)⟩ρ (21)

where the operators in square brackets are simply the spectral representations of and . Furthermore, it follows from (II.3) that

 E[eβa(ωa−ω′a)+βb(ωb−ω′b)]t=∑pq;nme−βaωpae−βbωqbZaZb ×⟨pa,qb|U(t)|na,mb⟩⟨na,mb|U−1(t)|pa,qb⟩ =∑pqe−βaωpae−βbωqbZaZb=1 (22)

where again we have used the completeness of the bare states. From (II.3), (II.3), (II.3) and the Jensen inequality , it now follows that and hence

 ΔS0≥0. (23)

This result can be viewed as an entropic version of the second law of thermodynamics (ESL) pertaining to the quantum system. The existence of this ESL supports the identifications of heat and entropy initially made in equation (II.3). We remark that these identifications have also been employed elsewhere Kurcz et al. (2010b); Blasone et al. (2011). It is therefore of interest to ask whether or not the CSL also holds when these identifications are made.

In order to go further and make contact with the CSL, we must resort to our assumption (9), which implies

 ΔQa=−ΔQb,ΔQa→b=2ΔQb. (24)

From this it follows using (23) that

 (βb−βa)ΔQa→b≥0, (25)

which is nothing but the CSL (17).

Note that in the present context the CSL is a significantly stronger statement than the ESL, because it requires the additional nontrivial assumption (9). While the ESL (23) holds quite generally for the initial thermal state (15) regardless of the form of , the CSL (17) will not generally hold when condition (9) is not met.

### ii.4 General form of heat transfer

In order to obtain a general expression for the heat transfer we assume that the operators and can be expanded as the following linear combinations of operators at

 a(t)=fa(t)a+ga(t)a†+fb(t)b+gb(t)b† b(t)=pa(t)a+qa(t)a†+pb(t)b+qb(t)b† (26)

where and . Such expansions result whenever the interaction Hamiltonian contains only linear and quadratic terms in the operators . This includes almost all interactions of practical interest. A straightforward calculation now yields

 ΔQa= ωa[(|fa|2+|ga|2−1)Xa +(|fb|2+|gb|2)Xb+|ga|2+|gb|2] ΔQb= ωb[(|pb|2+|qb|2−1)Xb +(|pa|2+|qa|2)Xa+|qa|2+|qb|2]. (27)

Here

 Xc=1eβcωc−1 (28)

with . The calculation of now reduces to the determination of the time-dependent functions in (II.4) for a given choice of interaction Hamiltonian .

## Iii Violations of the CSL

In this section we consider the case of two interacting oscillators, and demonstrate violations of the CSL for a particular choice of interaction Hamiltonian.

### iii.1 Forms of interaction for which the CSL holds

We have established that whenever condition (9) is met the CSL will hold. In many areas of quantum theory, such as quantum optics, interactions satisfying this condition are usually obtained through approximations of Hamiltonians derived from first principles. The most obvious example is the famous Jaynes-Cummings model, which has become the paradigm of solvable light-matter systems. The Jaynes-Cummings interaction results as a so-called rotating-wave approximation (RWA) of a more fundamental Hamiltonian for which condition (9) does not hold. The RWA eliminates terms in the interaction which do not conserve the bare energy .

In the context of the interacting oscillator system we are considering the RWA could for example take the form

 V=ig(a†+a)⊗(b†−b)⟶ig(a⊗b†−a†⊗b) (29)

where is a real coupling parameter. Note that before the RWA , while afterwards . Within this RWA the dynamics of the interacting system are especially simple to solve. Assuming for further simplicity the resonance condition we obtain the solution

 a(t) =e−iωt[acosgt−bsingt], b(t) =e−iωt[bcosgt+asingt]. (30)

Comparison with (II.4) shows that for this choice of interaction

 fa=e−iωtcosgt=pb,fb=−e−iωtsingt=−pa, (31)

while . Substitution of these expressions into (II.4) then yields

 ΔQa→b=2ω(Xa−Xb)sin2gt. (32)

Since according to (28) we have that we have verified explicitly that the CSL holds for the rotating-wave approximated interaction. Figure 1 illustrates the behavior of , and clearly shows that it is quite in-line with one’s intuition regarding the CSL.

### iii.2 Violations of the CSL

There are numerous interactions one might consider in order to exhibit a violation of the CSL. Indeed, any of the more fundamental interactions to which the RWA is usually applied do not have the property . They therefore open the door to violations of the CSL. Examples of such interactions include minimal coupling of oscillators and linear-coupling of oscillators. The former of these is the more physical of the two, as the minimal coupling prescription is a major ingredient in elementary theories that are built from first principles. However, because it contains additional self energy terms (terms quadratic in for example) its associated dynamics are more difficult to solve (though still quite possible).

For simplicity we consider the straightforward approach of avoiding the RWA as prescribed in (29) thereby adopting the linear interaction

 V=ig(a†+a)⊗(b†−b). (33)

Again, for simplicity of solution, we retain the resonance condition . This is despite the fact that considering off-resonant oscillators actually allows one to exhibit additional violations of the CSL.

When on resonance the introduction of new modes and decouples the equations of motion, which can then be solved simply giving

 c±(t)= c±[cosμ±t−i(ω±g)μ±sinμ±t]+c†±igμ±sinμ±t (34)

where . From these solutions the solutions for and are easily found and the time-dependent functions in (II.4) can then be read-off. Subsequent substitution of these functions into (II.4) then yields expressions for and .

Figure 2 shows how the behavior in can differ in an essential way when . Specifically, it is clear that , and that actually the hotter (as well as the cooler) oscillator begins to absorb heat (bare energy). This behavior by itself does not necessarily imply a violation of the CSL, because when is not conserved one cannot necessarily interpret the heat absorbed by the hotter oscillator as having been transferred from the cooler oscillator. Put differently, we have proved that is a sufficient condition in order that the CSL holds, but we have not proved that it is also a necessary condition.

To determine whether the CSL is truly violated the relevant quantity is . Figure 3 shows how oscillates with and becomes negative even when . In fact also oscillates when is fixed but varies. A critical value occurs at , at which point making the quantities of interest singular.

The amplitudes of the oscillations in begin to rapidly diverge as the value of surpasses , a regime which can be classed as the strong coupling regime. However, for sufficiently weak coupling, , the non-divergent oscillatory nature of suggests an interpretation of the violations of the CSL as transient phenomena. This interpretation is bolstered by the behavior of the time-averaged heat transfer

 ¯¯¯¯¯¯¯¯¯ΔQa→b(τ)=1τ∫τ0dtΔQa→b(t) (35)

plotted in figure 4, which shows that the CSL is only violated for very small values of .

The time-averaged heat transfer is ill-defined when , and for large values of , i.e., in the strong coupling regime such that , behaves quite differently. The interpretation of the violations as transient is no longer available in this regime as the behavior of indicates extensive violations of the CSL for large (figure 5).

### iii.3 Discussion

There are at least two ways in which the violations of the CSL exhibited in section III.2 can be interpreted. If one assumes the bare energy operators and represent the observable (heat) energies of the subsystems, then the results in section III.2 show that the CSL is genuinely violated, albeit transiently, within certain parameter regimes and for certain interaction Hamiltonians. In quantum optics such violations might be explained in terms of the energy-time uncertainty principle, which allows for the violation of bare energy conservation, and therefore the spontaneous production of heat, over very short timescales. There is experimental evidence to support such violations Wang et al. (2002); Gieseler et al. (2014). If however, such violations are viewed as untenable, even transiently, then either the physically available interactions between quantum subsystems become limited by the CSL, or one must reject the physicality of the free energy operators as representing observable energies. With regard to the former of these interpretive options, it is worth remarking that the RWA interaction in (29) can be obtained as a unitary transformation of the minimal coupling interaction Drummond (1987); Baxter et al. (1990). Thus, forms of interaction that necessarily obey the CSL need not be the result of an approximation.

## Iv Non-standard definitions of heat, entropy and subsystems

In this section we extend the analysis of the previous sections by considering non-standard definitions of the basic thermodynamic quantities and quantum subsystems.

### iv.1 True subsystem energies via non-standard quantum subsystems

The first law of thermodynamics requires the conservation of the total energy of an isolated system. The term isolated is preferred in classical thermodynamics, whereas in quantum theory the same kind of system is referred to as closed. A closed quantum system is characterised by its unitary dynamics, and the composite -oscillator system that we are considering is of precisely this type.

Since the entropy in (II.3) is defined in terms of the free energies and rather than , the first law that corresponds to the ESL (23) could be taken as the requirement , which is nothing but condition (9). One can therefore understand the violations of the CSL as a consequence of violations of this first law. However, the condition , which expresses the conservation of the true total energy, and which necessarily holds for a closed quantum system with time-independent Hamiltonian, could also be interpreted as the first law of thermodynamics in the present context.

Our results show that the CSL will hold when the first and second laws are assumed to refer to the same type of energy , and the first law is assumed to hold in the form . Alternatively, if one takes as expressing the first law, we see that the first and second laws refer to two distinct types of energy whenever there are interactions. In this latter case, it is possible to violate the CSL even though the first law necessarily holds, because the sum of subsystem energies does not coincide with the true total energy . In order that , we must have either that , or we must redefine the subsystem energies and .

To investigate the option of redefining the subsystem energies we consider the case of minimally coupled oscillators. Introducing position and momentum variables as

 xc=√12mω(c†+c),pc=i√mω2(c†−c) (36)

where and , and minimally coupling oscillator to through a coupling parameter gives the Hamiltonian

 H=12m(pa−qxb)2+12mω2x2a+p2b2m+12mω2x2b (37)

where for notational simplicity tensor products with the identity operators for and have been omitted. Alternatively, one could minimally couple oscillator to so as to give the Hamiltonian

 H′=p2a2m+12mω2x2a+12m(pb+qxa)2+12mω2x2b. (38)

The two Hamiltonians and are unitarily related by the transformation . As a result they produce identical equations of motion written in terms of the derivatives of and 1. However the non-local (in the sense of tensor product structure) transformation , which does not commute with and results in these operators taking on different physical meanings with respect to the distinct Hilbert space representations connected by . Because of this the free energies and are also physically different within the two different Hilbert space representations. We see that conventional subsystem observables are highly non-unique and therefore physically ambiguous as noted in Section I.

One can define unique subsystem observables in terms of , and their velocities. In particular both and can be written in the invariant form

 H=12m(˙x2a+ω2x2a)+12m(˙x2b+ω2x2b)=H′ (39)

showing that they are in fact identical Hamiltonians, with alternative forms in (37) and (38) that result from their superficial expression in terms of physically distinct canonical momenta. It is clear from (39) that one can define unique and physically unambiguous subsystem energies and such that

 Htruec=12m(˙x2c+ω2x2c) (40)

with . These definitions are especially natural in that they express the total energy as the sum of subsystem energies. The subsystems are moreover defined in terms of precisely the same physical variables as in the free (non-interacting) case, namely the positions and , and the mechanical momenta and .

It is important to recognise that this definition of a quantum subsystem cannot coincide with the conventional one, because and do not commute. This means that there are no operators and such that and simultaneously . As such the physical subsystems are not determined through a (physically non-unique) mathematical decomposition of the Hilbert space based on tensor product structure, but rather in terms of their assumed to be fundamental dynamical variables.

With this definition of a quantum subsystem the interaction between and is described by the non-commutativity of the velocity operators, rather than through an interaction Hamiltonian , and it is not possible to simultaneously prepare the subsystems in states of well-defined energy, or well-defined temperature as in (15). One cannot then pose sharp questions about heat transfer between the two systems and the apparent violations of the CSL in section III.2 lose any physical significance.

It is natural to ask whether there are, from the new point of view, any conditions under which questions about heat transfer can be meaningfully posed. By way of an answer, note that according to (37) while , and yet according to (38) while . Since performing the RWA and neglecting the self-energy terms in either or produces the same resulting Hamiltonian (when the oscillators are resonant), consistency requires that within the regime in which these approximations are justified we must have that and . Thus, the rotating-wave form of the Hamiltonian is the unique form in which both of the bare energy operators approximately coincide with the corresponding true energies.

From the current perspective, the initial state in (15) cannot truly specify the individual states of the subsystems, because the true subsystems do not specify a tensor product decomposition of the composite Hilbert space such that the factors represent the physical subsystem state spaces. At best the initial state in (15) can be considered an approximation, and the analysis above demonstrates that the conditions under which this approximation is valid are precisely those conditions under which the CSL cannot be violated. We can therefore explain the apparent violations of the CSL as artifacts of an unwarranted attribution of sharp temperatures to the subsystems, when such attributions are not strictly permissible, and can hold only approximately. This explanation is similar in nature to that offered by the energy-time uncertainty relation, which suggests that there is always some uncertainty in the (heat) energies of the subsystems, so that on short enough timescales apparent violations of the CSL may occur.

### iv.2 Entropy production and heat exchange

In Esposito et al. (2010); Esposito and Broeck (2011) Esposito et al. investigate the production of entropy within a system coupled to several thermal reservoirs in thermal equilibrium. Since no constraints are placed on what actually constitutes a reservoir, the treatment there is easily adapted to the situation considered here such that each oscillator is considered as a reservoir for the other. The authors in Esposito et al. (2010); Esposito and Broeck (2011) show that the total change in the von Neumann entropy of system can be written assuming an initially separable state

 ΔSa(t) =−tr[ρa(t)lnρa(t)]+tr[ρalnρa] =ΔiSa(t)+ΔeSa(t) (41)

where denotes the irreversible entropy production within system and denotes the reversible entropy-flux due to the reservoir system . Explicitly these contributions are given by

 ΔiSa(t) =S(ρ(t)||ρa(t)⊗ρb) ΔeSa(t) =βb(⟨Hb(0)⟩−⟨Hb(t)⟩)=−βbΔQb(t) (42)

where denotes the relative entropy between states and , and where denotes the contribution from system to the free entropy change given in (II.3).

One possibility at this point is to use the Gibbs relation (12) and the identifications made in (IV.2) to define the heat flux’s into systems and by

 ~ΔQa(t) =1βaΔeSa(t)=−βbβaΔQb(t), ~ΔQb(t) =1βbΔeSb(t)=−βaβbΔQa(t). (43)

These definitions are essentially the reverse of those given in section II. As such the proof of the CSL given these definitions follows almost exactly as in II. Assuming we obtain the CSL in the form

 (βa−βb)~ΔQa→b=−[βaβb+βbβa]ΔS0≤0. (44)

Here, as in section II, when the CSL does not follow from the ESL , and the results in section III.2 remain valid.

Since the above identification (IV.2) uses the Gibbs relation as in section II it is not suprising that the proof of the CSL requires the same additional assumption (9). However, in Esposito et al. (2010); Esposito and Broeck (2011) the connection between the heat flux defined in terms of and the first law of thermodynamics is discussed without ever involving the Gibbs relation. If we consider system as a reservoir for then the heat flux from the reservoir into can be identified as . The authors in Esposito et al. (2010); Esposito and Broeck (2011) use this definition and the fact that to identify the heat and work contributions to the energy change within system as follows

 ⟨ΔHtruea⟩ =⟨ΔQa⟩+⟨ΔWa⟩, ⟨ΔQtruea⟩ =∫t0dt′tr[(Ha(t)+V(t))˙ρa(t)], ⟨ΔWtruea⟩ =∫t0dt′tr[(˙Ha(t)+˙V(t))ρa(t)]. (45)

In our situation is not explicitly time-dependent in the Schrödinger picture so the work contribution vanishes and

 ⟨ΔHtruea⟩ =⟨ΔQtruea⟩ =tr[(Ha+V)ρa(t)]−tr[(Ha+V)ρa]. (46)

According to (IV.2) the physical energy of system is an observable as in section IV.1. This true energy is defined in terms of operators that according to the standard definition of a quantum subsystem are associated with system . Therefore, according to this definition the true subsystems and are not standard quantum subsystems. Furthermore, since we are considering a situation in which the labels “system” and “reservoir” can be interchanged, we see that naively applying the same method as above in order to obtain the energy-flux into system from system gives

 ⟨ΔQtrueb⟩=tr[(Hb+V)ρb(t)]−tr[(Hb+V)ρb]. (47)

Now, if the energy lost by system () must equal that gained by system () then (47) cannot generally be reconciled with (IV.2) and (IV.2) using only a single tensor-product decomposition of the Hilbert space into subsystem state spaces. This means one must either identify subsystems in a non-standard way using more than one tensor product decomposition, or allow the non-conservation of the sum of subsystem energies. The latter option means allowing , provided of course, that we use the standard definition of subsystem energies. If, instead, one uses within one and the same tensor-product decomposition, the non-standard definition where

 Htruea=H−Hb,   Htrueb=H−Ha (48)

as supplied in Esposito et al. (2010); Esposito and Broeck (2011) and used in (IV.2) and (47), then one obtains , which is not generally a conserved quantity. Moreover, adopting (48) the heat-flux between the systems actually remains unchanged;

 ΔQtruea→b :=⟨Htruea(t)⟩−⟨Htrueb⟩−⟨Htrueb(t)⟩+⟨Htrueb⟩ ≡ΔQa→b, (49)

and we know that the CSL for the quantity will not generally hold unless condition (9) is met.

As with the energies , the sum of the alternative energies in (48) will only coincide with the true total energy when . Thus, as explained in section IV.1, one might argue that within a single tensor-product decomposition the two equalities in (48) can only approximately hold simultaneously when the coupling is sufficiently weak; . The CSL can then be recovered within the regime for which using (48) in conjunction with a single tensor-product decomposition constitutes a valid approximation.

In summary then, if we adopt as in Esposito et al. (2010); Esposito and Broeck (2011), the alternative identifications of entropy and heat given by (IV.2) and (IV.2) we reach essentially the same conclusions as before; either we accept that the CSL is violated, or we accept the restrictions it places upon the form of the interaction, or we adopt a non-standard definition of quantum subsystem. It is this last option that appears to have been implicitly assumed in Esposito et al. (2010).

### iv.3 The CSL with non-standard subsystem energies

The identification of the true heat flux made in (IV.2) Esposito et al. (2010); Esposito and Broeck (2011), is consistent with the analysis made in section IV.1 whereby it was argued that the initial thermal state in (13) cannot be justified when (9) does not hold. It is obviously justified however, if at the interaction is adiabatically switched on. This gives a time-dependent Hamiltonian such that . Since the Hamiltonian is time-dependent the work contribution to the expected energy no longer vanishes.

We now adopt the non-standard definition of subsystem energies used in section IV.1, and assume without loss of generality that work is performed on system . We can visualise this situation through some mechanism that gradually introduces system to system , which remains stationary. Such a situation could be modelled in the case of minimally coupled oscillators by using the Hamiltonian in (37) with the coupling constant replaced by a time-dependent function . The energy operators for the oscillator subsystems then become

 Htruea(t)=12m(˙xa(t)2+x2a),Hb=12m(˙x2b+x2b) (50)

with