# Relative Entropy, Interaction Energy and the Nature of Dissipation

Bernard Gaveau, Léo Granger, Michel Moreau, Lawrence S. Schulman

► To cite this version:

Bernard Gaveau, Léo Granger, Michel Moreau, Lawrence S. Schulman. Relative Entropy, Interaction Energy and the Nature of Dissipation. Entropy, 2014, 16 (6), pp.3173-3206. 10.3390/e16063173. hal-01316773

HAL Id: hal-01316773

https://hal.sorbonne-universite.fr/hal-01316773

Submitted on 17 May 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Relative Entropy, Interaction Energy and the Nature of Dissipation

Bernard Gaveau<fnref n="1" />, Léo Granger<fnref n="2" />, Michel Moreau<fnref n="3" /> and Lawrence S. Schulman<fnref n="4" />,*

Received: 10 February 2014; In revised form: 20 May 2014 / Accepted: 23 May 2014 /

Published: 6 June 2014

**Abstract:** Many thermodynamic relations involve inequalities, with equality if a process does not involve dissipation. In this article we provide equalities in which the dissipative contribution is shown to involve the relative entropy (a.k.a. Kullback-Leibler divergence). The processes considered are general time evolutions both in classical and quantum mechanics, and the initial state is sometimes thermal, sometimes partially so. By calculating a transport coefficient we show that indeed—at least in this case—the source of dissipation in that coefficient is the relative entropy.

**Keywords:** relative entropy; dissipation; thermodynamic inequalities; interaction energy

## 1. Introduction

The distinction between heat and work, between the uncontrollable flow of energy of molecular processes and the controllable flow of energy usable by an agent, underlies all of thermodynamics, and is implicitly incorporated in the equation $dE = dW + dQ$. This distinction is defined by human subjectivity and by the human technological ability to extract work from the flow of energy of microscopic processes.

On the other hand, assuming that the evolution is given by the fundamental laws of dynamics, classical or quantum, one must represent the evolution of a system at a fundamental microscopic level by the action of a unitary operator on a quantum state (density operator) in quantum situations or by the action

<footnote n="1">¹ Laboratoire analyse et physique mathématique, 14 avenue Félix Faure, 75015 Paris, France</footnote>

<footnote n="2">² Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, D-01187 Dresden, Germany</footnote>

<footnote n="3">³ Université Pierre et Marie Curie, LPTMC, case 121, 4 pl. Jussieu, 75252 Paris cedex 05, France</footnote>

<footnote n="4">⁴ Physics Department, Clarkson University, Potsdam, New York 13699-5820, USA</footnote>

<footnote>* Author to whom correspondence should be addressed; E-Mail: schulman@clarkson.edu.</footnote>
of a symplectic operator on a classical state (probability distribution) in classical situations. The state of the system evolves according to the Heisenberg equation of motion or according to the Liouville equation respectively. An immediate consequence is that *the entropy of the exact microscopic state of a system that is isolated or is coupled only to an external source of work stays constant during the evolution*. In particular, the microscopic state cannot tend to an *equilibrium state*. Thus, the information content of the exact microscopic state stays constant during the evolution. The problem is that in practice it is impossible both to define the microscopic state and to follow its exact evolution. The definition of the state and the representation of the exact evolution using unitary or symplectic dynamics are thus untenable idealizations. Nevertheless, and this constitutes a paradox, these idealizations cannot be ignored or dismissed, because it is precisely the difference between the exact evolution and its standard approximations which explains and can be used to predict dissipative effects, both of energy and information.

In thermodynamics, in kinetic theories or in stochastic dynamics, the exact microscopic state of a system is replaced by an approximate or “coarse-grained” state and the corresponding exact evolution is replaced by an evolution of the corresponding coarse-grained state (or, in standard thermodynamics by a quasi-static or formal evolution). There are two main reasons for using these approximate states and evolutions:

(1) As discussed, it is impossible—even in principle—to specify the exact state of a large system and follow its evolution. An attempt at extremely high precision would modify the system, even in a classical context (related to Maxwell’s demon). And it is even worse for quantum systems. Moreover, this would be useless.

(2) Only slow variables (on the time scales of microscopic processes) can be measured with confidence and stability. As a result, an observer can only describe the system as a state of minimal information (or maximal entropy) compatible with the observed slow variables [1,2].

The coarse-grained state is thus a statistical data structure which summarizes at a given moment the knowledge of the observer. The evolution of this coarse-grained state merely reflects the evolution of the knowledge of the observer about the system. The observer cannot follow the microscopic processes, but only the slow variables which can be measured and used, and as a consequence there is a loss of information about the details of the microscopic processes; in the traditional language of thermodynamics, entropy increases or is produced. The observer, reflecting a particular state of knowledge (or more precisely, a lack of knowledge), describes the state of the system as a state of minimal information (or maximal entropy) compatible with observation. Thus, entropy is not a kind of substance flowing from one part of a system to another part or mysteriously produced internally by the physical system, as is often suggested by many texts of thermodynamics or statistical physics: it is only the observer's partial inability to relate the exact microscopic theory to a reduced macroscopic description in order to use the system as a source of useful work or information. This is what is measured as an increase of entropy or by entropy production. The macroscopic state is the result of a statistical inference (specifically, maximum entropy) for the given, observed, macroscopic variables (which are the slow variables of the system [1,2]). This point of view on the nature of entropy was emphasized by Jaynes, who observed [3,4], “The expression ‘irreversible process’ represents a semantic confusion…”
The difference between the exact evolution of the microscopic ideal state and the evolution of its coarse-grained approximation is what is called “dissipation,” both of information content and of energy or other “useful” variables. Standard thermodynamics uses the maximal coarse graining of equilibrium, and the idealized evolution is not modeled explicitly, so dissipation can be taken into account only by inequalities. For more detailed coarse-graining (as in hydrodynamics, Boltzmann’s equation, kinetic theories or stochastic thermodynamics) one can obtain an estimate for the dissipative effects, for example, by the calculation of transport coefficients.

In this article, our main purpose is to prove that the relative entropy term between the initial and final states measures dissipation. In our approach, “dissipation” is defined as the difference between the maximal work that the physicist thinks could be extracted from a system when using the thermodynamic or quasi-static theory to make predictions, and the work that is actually extracted because the system is evolving according to the exact dynamics, classical or quantum, independently of what the physicist thinks (see also our use of relative entropy in [5] and [6], where the context was more limited [7]). Moreover, in the present context we find that the relative entropy terms are proportional to the square of the interaction energy. In all standard theories, dissipative effects are measured by the transport coefficients of energy or momentum or concentration of chemical species. Thus, we need prove that the relative entropy allows the calculation of transport coefficients. Indeed, we show below that the relative entropy terms provide the calculation of the thermal conductivity between two general quantum systems, initially at thermal equilibrium at different temperatures. This is a kind of Fourier law, except that we do not suppose a linear regime, so that the temperature dependence is more complicated than simple linearity. Moreover, our exact calculation of the transport coefficient shows that it is indeed proportional to the square of the interaction energy, which confirms that for vanishingly small interaction energy no transfer occurs in finite time. In other words, no power or finite rate of information flow can be extracted from a system if one does not have at the same time dissipative effects.

In the following material, we first consider a system comprised of two components, *A* and *B*. We make no specific hypotheses on the size of the systems, and we do not introduce thermal reservoirs. Thus, the identities we derive are in effect exact tautologies. In Sections 3 to 5, we present several identities. We here mention two examples: (1) a derivation of the Brillouin-Landauer estimate of the energy necessary to change the information content of a system; (2) an estimate of the work that can be extracted from a two-part system in interaction with an external source of work in terms of non-equilibrium free energies and relative entropy of the state before and after the evolution. Similar identities were also obtained recently by Esposito *et al.* [8], Reeb and Wolf [9], and Takara *et al.* [10] Continuing, we study the effect of an external agent on an (otherwise) isolated system; again we obtain an identity relating the work to the difference of internal (not the free) energies along with the usual dissipative terms. Then, we derive the relation between the relative entropy and the heat conductivity in a quantum system. Finally, we define a general notion of coarse-graining or reduced description, which includes the usual notions. In some of our examples one or both systems are initially at thermal equilibrium, but only the initial temperatures appear explicitly in the definition of the non-equilibrium free energies. The latter are no longer state functions because they depend explicitly on the initial temperature and not on the actual effective temperature. No coarse graining by an effective final or intermediate thermal state is used, and neither system is a reservoir.
## 2. Notations and Basic Identities

### 2.1. States and Entropy

Many results will be valid both in classical and quantum contexts. We denote by $\rho$ either a probability distribution function over a classical phase space, or a density matrix in the quantum case. We denote by Tr either the integral on the phase space, or the trace operation. Thus $\rho$ is a positive quantity and satisfies Tr $\rho = 1$. The entropy of $\rho$ is

$$S(\rho) = - \mathrm{Tr} \rho \log \rho. \qquad (1)$$

It is defined up to a multiplicative constant. (Classically $\rho$ should be divided by a dimensional constant to render it dimensionless.)

The relative entropy (see [11]) is defined by

$$S(\rho|\rho') = \mathrm{Tr} (\rho (\log \rho - \log \rho')) , \qquad (2)$$

where $\rho$ and $\rho'$ are states.

One has

$$S(\rho|\rho') \ge 0, \qquad (3)$$

and $S(\rho|\rho')$ does not depend on the units in phase space. Moreover $S(\rho|\rho') = 0$ if and only if $\rho = \rho'$.

Writing $S(\rho|\rho')$ as $-\mathrm{Tr} \rho \log \rho' - (-\mathrm{Tr} \rho \log \rho)$, suggests the following interpretation: Suppose the true state is $\rho$, but the observer thinks that the state is $\rho'$. $S(\rho|\rho')$ is then the true average of the missing information minus the estimate of the missing information.

### 2.2. The Basic Identity

If we add and subtract $S(\rho')$ in the second member of Equation (2), we obtain the basic identity

$$S(\rho|\rho') = S(\rho') - S(\rho) - \mathrm{Tr} ((\rho - \rho') \log \rho') . \qquad (4)$$

Most of our results follow from this identity.

When $\rho'$ is a thermal state at (inverse) temperature $\beta$ [12],

$$\rho' = \rho_{\beta} = \frac{e^{-\beta H}}{Z(\beta, H)}, \qquad (5)$$

where

$$Z(\beta, H) = \mathrm{Tr} e^{-\beta H} \qquad (6)$$

is the partition function. With $\rho' = \rho_{\beta}$, the identity (4) reduces to

$$S(\rho|\rho_{\beta}) = S(\rho_{\beta}) - S(\rho) + \beta \mathrm{Tr} ((\rho - \rho_{\beta})H) . \qquad (7)$$

Here $H$ is a given function or operator.

Defining the free energy of state $\rho$ by

$$F(\rho, H) = \mathrm{Tr}(\rho H) - \frac{1}{\beta} S(\rho), \qquad (8)$$
we obtain

$$S(\rho|\rho_{\beta}) = \beta(F(\rho, H) - F(\rho_{\beta}, H)), \qquad (9)$$

and $F(\rho_{\beta}, H)$ is the equilibrium free energy related to the partition function by

$$Z(\beta, H) = \exp(-\beta F(\rho_{\beta}, H)). \qquad (10)$$

Equation (9) is important for applications, because its right hand side can be related to energy dissipation (up to the factor $\beta$), which gives a clear physical meaning to the relative entropy (see Section 3).

## 2.3. Evolution Operators and Entropy

We assume that the system (classical or quantum) evolves under the action of an arbitrary operator $U$ (symplectic or unitary). If $\rho$ is a state, we denote by $\rho^{(U)}$ the new state after the evolution $U$.

Entropy is conserved by the evolution

$$S(\rho^{(U)}) = S(\rho). \qquad (11)$$

For example, in the quantum case, we have $\rho^{(U)} = U\rho U^{\dagger}$, where $U$ is the propagator: $i\frac{dU}{dt} = [H, U]$, $U|_{t=0} = 1$, with $H$ a possibly time-dependent Hamiltonian.

If $\phi(\rho)$ is a functional of $\rho$ which evolves with $U$, and $\phi(\rho^{(U)})$ is the functional after evolution of $\rho$, we denote the variation of $\phi(\rho)$ after the evolution $U$ in the following way

$$\delta^{(U)}(\phi(\rho)) = \phi(\rho^{(U)}) - \phi(\rho). \qquad (12)$$

**Remark 1:** Many of our results are valid for a general evolution $U$ which is not symplectic or unitary, for example stochastic evolution.

# 3. Two Systems in Interaction

A basic procedure in thermodynamics is to consider the evolution and properties of an otherwise isolated two-part system. Although it is often the case that the overall system conserves energy, for the subsystems more general behavior is often seen.

## 3.1. Hypotheses

We assume that the system is formed of two parts, $A$ and $B$, in interaction. At time-0, the state is a product state

$$\rho_0 = \rho_{A,0} \otimes \rho_{B,0}. \qquad (13)$$

After the evolution $U$, the state is $\rho^{(U)}$ and we denote by $\rho_A^{(U)}$ and $\rho_B^{(U)}$ its marginals,

$$\rho_A^{(U)} = \mathrm{Tr}_B \rho^{(U)} \quad \text{and} \quad \rho_B^{(U)} = \mathrm{Tr}_A \rho^{(U)}, \qquad (14)$$

which are then states on $A$ and $B$ respectively. We also assume that there is a quantity $H$ that is conserved by the evolution and $H$ has the form

$$H = H_A + H_B + V_{AB}, \qquad (15)$$
where $H_A$ and $H_B$ are quantities depending only on A and B respectively and $V_{AB}$ is an interaction term.
Then, if we denote

$$E(\rho) = \mathrm{Tr}(\rho H) = E_A(\rho) + E_B(\rho) + E_V(\rho) \qquad (16)$$

$$E_A(\rho) = \mathrm{Tr}(\rho H_A) = \mathrm{Tr}(\rho_A H_A) \qquad (17)$$

$$E_B(\rho) = \mathrm{Tr}(\rho H_B) = \mathrm{Tr}(\rho_B H_B) \qquad (18)$$

$$E_V(\rho) = \mathrm{Tr}(\rho V_{AB}), \qquad (19)$$

our hypothesis is that

$$\delta^{(U(t))} E(\rho) \equiv E(\rho^{(U)}) - E(\rho_0) = 0. \qquad (20)$$

In particular this is the case if $U$ is time-evolution with Hamiltonian $H$.

**Remark 2:** For this situation, certain results are also valid without the assumption that the evolution $U$ preserves the energy $H$.

If $\rho$ is a state corresponding to a system formed of two parts, A and B, and $\rho_A$ and $\rho_B$ are its marginals (as in Equation (14)), then the relative entropy $S(\rho|\rho_A \otimes \rho_B)$ is the same as the mutual information of the associated distributions. It can be interpreted as the amount of information in $\rho$ that comes from the fact that A and B are in interaction (see [11]). This quantity will appear in many of our relations below (e.g., Equation (28)) as part of the dissipation.

## 3.2. Relation between a State and Its Marginals

Assuming Equation (13) (that the initial state is a product state), one has the identity

$$\delta^{(U(t))} S(\rho_A) + \delta^{(U(t))} S(\rho_B) = S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}). \qquad (21)$$

Indeed, using the conservation of the entropy of $\rho$ during the evolution $U$,

$$\begin{aligned} \delta^{(U(t))} S(\rho_A) + \delta^{(U(t))} S(\rho_B) &= S(\rho_A^{(U)}) + S(\rho_B^{(U)}) - S(\rho_0) \\ &= S(\rho_A^{(U)}) + S(\rho_B^{(U)}) - S(\rho^{(U)}) = S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}). \end{aligned}$$

This is because one evidently has $-\mathrm{Tr}(\rho^{(U)} \log \rho_A^{(U)}) = -\mathrm{Tr}(\rho_A^{(U)} \log \rho_A^{(U)})$. Note that Equation (21) requires that $U$ preserve the entropy. One has also the well-known inequality

$$S(\rho^{(U)}) \le S(\rho_A^{(U)}) + S(\rho_B^{(U)}), \qquad (22)$$

which is a particular case of

$$S(\rho) \le S(\rho_A) + S(\rho_B) \qquad (23)$$

for any state $\rho$.

Note that the stronger result, Equation (21), is obtained by retaining the relative entropy term in this equation. The same remark will apply in most of the following results.
### 3.3. The Case Where A is Initially in a Thermal State

At time 0 we take $\rho_{A,0}$ to be thermal with temperature $\beta_A$,

$$ \rho_{A,0} = \rho_{A,\beta_A} = \frac{e^{-\beta_A H_A}}{Z_A(\beta_A)}, \qquad (24) $$

where $Z_A(\beta_A) = Z_A(\beta_A, H_A)$ is the partition function, (6). From Equation (7) with $\rho \to \rho_A^{(U)}$ and $\rho_\beta \to \rho_{A,\beta_A}$, we deduce (note that this requires that $H_A$ be independent of time)

$$ \delta^{(U(t))} S(\rho_A) - \beta_A \delta^{(U(t))} E_A(\rho_A) = -S(\rho_A^{(U)}|\rho_{A,\beta_A}), \qquad (25) $$

and as a consequence

$$ \delta^{(U(t))} S(\rho_A) - \beta_A \delta^{(U(t))} E_A(\rho_A) \le 0. \qquad (26) $$

The last two equations do not require that $U$ be a unitary evolution conserving the entropy, nor that it conserve the energy.

**Remark 3:** This inequality can be found in [13] as an unnumbered equation. Its consequences were not deduced in that reference.

**Remark 4:** Note that it is the initial temperature that appears in Equations (25) and (26). Moreover, $\rho_A^{(U)}$ is not in general an equilibrium state.

Suppose that $B$ starts in an arbitrary initial state $\rho_{B,0}$, while $A$ begins in the thermal state $\rho_{A,\beta_A}$. Combining Equations (21) and (25), we obtain

$$ \beta_A \delta^{(U(t))} E_A(\rho) + \delta^{(U(t))} S(\rho_B) = S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_A^{(U)}|\rho_{A,\beta_A}). \qquad (27) $$

The last equation requires that $U$ preserve entropy, since that feature is used in the derivation of Equation (21). It also remains valid if the Hamiltonian of $B$, $H_B$, depends on an external parameter varying with time, so that $B$ receives work from an external agent. This is because the entropy-preserving property only depends on $U$ being unitary (or symplectic). On the other hand, $H_A$ should be time independent (see the parenthetical remark before Equation (25)). Then if $U$ conserves energy

$$ \beta_A (\delta^{(U(t))} E_B(\rho) + \delta^{(U(t))} E_V(\rho)) = \delta^{(U(t))} S(\rho_B) - [S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_A^{(U)}|\rho_{A,\beta_A})]. \qquad (28) $$

These relations imply the following inequalities:

(1) If $U$ preserves entropy, even if $H_B$ depends on an external parameter varying with time

$$ \beta_A \delta^{(U(t))} E_A(\rho) \ge -\delta^{(U(t))} S(\rho_B). \qquad (29) $$

(2) If $U$ conserves entropy and the total energy, one has

$$ \beta_A (\delta^{(U(t))} E_B(\rho) + \delta^{(U(t))} E_V(\rho)) \le \delta^{(U(t))} S(\rho_B), \qquad (30) $$
with the following interpretations. Suppose U conserves the entropy; then we couple a system B (initially in an arbitrary state $\rho_{B,0}$) to system A (initially in thermal equilibrium) and that we want to lower the entropy of B so that $\delta^{(U(t))}S(\rho_B) \le 0$. Then, the energy of A must increase by at least

$$ \delta^{(U(t))} E_A(\rho) \geq \frac{1}{\beta_A} |\delta^{(U(t))} S(\rho_B)| \quad (31) $$

even if B receives work from an external source (so that $H_B$ depends on an external parameter). Moreover, if the total energy is conserved, the sum of the energy of B and the coupling energy must decrease by at least:

$$ \delta^{(U(t))} E_B(\rho) + \delta^{(U(t))} E_V(\rho) \leq \frac{1}{\beta_A} \delta^{(U(t))} S(\rho_B) < 0. \quad (32) $$

Thus lowering the entropy of a system B, coupled to a system initially at equilibrium, costs transfers of energy from B to A or to the interaction energy; thus the sum of B's energy and the interaction energy must decrease, but B's energy alone need not decrease. This is a result analogous to those of Brillouin [14] and Landauer [15] (reprinted in Leff and Rex [16]), even if system B receives work from an external source. But note again that only the temperature $\beta_A$ appears. This is the initial temperature at the beginning of the evolution U, so that system A is not necessarily a thermal bath, because its temperature may vary during the evolution U.

### 3.4. The Case of Equality in Equation (31)

It is important to study the case where the previous inequalities are changed into equalities, because this occurs if and only if strong conditions are verified: expressing these conditions is one of the advantages of our approach.

Equation (31) was derived under the hypothesis $\delta^{(U)}S(\rho_B) < 0$. Then by Equation (27) one has $S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) = S(\rho_A^{(U)}|\rho_{A,\beta_A}) = 0$. This implies $\rho_A^{(U)} = \rho_{A,\beta_A}$ and $\rho^{(U)} = \rho_A^{(U)} \otimes \rho_B^{(U)}$. Because $\rho_A^{(U)} = \rho_{A,\beta_A}$, $\delta^{(U)}E_A(\rho_A) = 0$, so by Equation (31), $\delta^{(U)}S(\rho_B) = 0$. But we could also derive $\delta^{(U)}S(\rho_B) = 0$ from Equation (21), because $\delta^{(U(t))}S(\rho_A) = 0$ and $S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) = 0$.

### 3.5. Both Systems A and B are at Equilibrium

Assume that A and B are initially at thermal equilibrium at different temperatures. Then, one has

(1) For a general evolution

$$ \delta^{(U(t))} S(\rho_A) - \beta_A \delta^{(U(t))} E_A(\rho_A) = -S(\rho_A^{(U)}|\rho_{A,\beta_A}), \quad (33) $$

$$ \delta^{(U(t))} S(\rho_B) - \beta_B \delta^{(U(t))} E_B(\rho_B) = -S(\rho_B^{(U)}|\rho_{B,\beta_B}). \quad (34) $$

(2) If U conserves entropy

$$ \delta^{(U(t))} S(\rho_A) + \delta^{(U(t))} S(\rho_B) = S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}). \quad (35) $$

(3) If U conserves energy

$$ \delta^{(U(t))} E_A(\rho_A) + \delta^{(U(t))} E_B(\rho_B) + \delta^{(U(t))} E_V(\rho) = 0. \quad (36) $$
Then, we conclude

(A) If $U$ conserves entropy: Combining Equations (33), (34), and (35) yields

$$ \beta_A \delta^{(U(t))} E_A(\rho_A) + \beta_B \delta^{(U(t))} E_B(\rho_B) = S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}), \quad (37) $$

It is easy to check directly that

$$ S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}) = S(\rho^{(U)} | \rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_A^{(U)} | \rho_{A,\beta_A}) + S(\rho_B^{(U)} | \rho_{B,\beta_B}). $$

Thus

$$ \begin{aligned} S(\rho^{(U)} | \rho_A^{(U)} \otimes \rho_B^{(U)}) &= \delta^{(U(t))} S(\rho_A) + \delta^{(U(t))} S(\rho_B) \\ &= \beta_A \delta^{(U(t))} E_A(\rho_A) + \beta_B \delta^{(U(t))} E_B(\rho_B) - [S(\rho_A^{(U)} | \rho_{A,\beta_A}) + S(\rho_B^{(U)} | \rho_{B,\beta_B})]. \end{aligned} \quad (38) $$

This last identity implies the Clausius-like inequality

$$ 0 \le \delta^{(U(t))} S(\rho_A) + \delta^{(U(t))} S(\rho_B) \le \beta_A \delta^{(U(t))} E_A(\rho_A) + \beta_B \delta^{(U(t))} E_B(\rho_B). \quad (39) $$

(B) For a general evolution $U$: Combining Equations (33) and (34)

$$ \begin{aligned} \delta^{(U(t))} E_A(\rho_A) + \delta^{(U(t))} E_B(\rho_B) = T_A & \left[ \delta^{(U(t))} S(\rho_A) + S(\rho_A^{(U)} | \rho_{A,\beta_A}) \right] \\ & + T_B \left[ \delta^{(U(t))} S(\rho_B) + S(\rho_B^{(U)} | \rho_{B,\beta_B}) \right], \end{aligned} \quad (40) $$

and thus

$$ \delta^{(U(t))} E_A(\rho_A) + \delta^{(U(t))} E_B(\rho_B) \ge T_A \delta^{(U(t))} S(\rho_A) + T_B \delta^{(U(t))} S(\rho_B). \quad (41) $$

This may be viewed as an inequality for free energies of A and B at the respective temperatures $T_A$ and $T_B$. Note again that during the time evolution neither A nor B need remain in thermal states.

### 3.6. Case of Equality in (39) (2nd) and (41)

(A) $U$ conserves entropy. Equality in Equation (39) implies immediately that $\rho_A^{(U)} = \rho_{A,\beta_A}$ and $\rho_B^{(U)} = \rho_{B,\beta_B}$, in which case the energy of A and the energy of B have not changed and $\delta^{(U(t))} S(\rho_A) = \delta^{(U(t))} S(\rho_B) = 0$. From the first equality in Equation (38) one has

$$ \rho^{(U)} = \rho_A^{(U)} \otimes \rho_B^{(U)} = \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}, \quad (42) $$

and one deduces that the state $\rho$ has not changed.

(B) General evolution $U$. If one has equality in Equation (41), it follows from Equation (40) the same results as above: the state $\rho$ has not changed.
### 3.7. Interaction Energy and Relative Entropy

It is often assumed that the interaction energy between the parts of the complete system can be neglected but, obviously, if this were exactly true the subsystems would evolve independently. Of course, an interaction can be small but nevertheless have significant impact when it persists for long times. However, there are cases where even for short times the interaction cannot be neglected. Assume then that $U$ conserves entropy and energy. Divide Equations (33) and (34) by $\beta_B$ and add; then use the conservation of energy Equation (36) to eliminate $\delta^{(U)}E_B(\rho)$ and deduce after some calculations

$$ - \delta^{(U(t))} E_V(\rho) = \left(1 - \frac{\beta_A}{\beta_B}\right) \delta^{(U(t))} E_A(\rho_A) + T_B S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}), \quad (43) $$

so that

$$ - \delta^{(U(t))} E_V(\rho) \geq \left(1 - \frac{\beta_A}{\beta_B}\right) \delta^{(U(t))} E_A(\rho_A). \quad (44) $$

In case of equality in Equation (44), one deduces that $\rho^{(U)} = \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}$ so that the state has not changed and $\delta^{(U)}E_A(\rho_A) = \delta^{(U)}E_V(\rho) = 0$. Moreover if $\delta^{(U)}E_A(\rho_A)$ is positive and $T_A$ is larger than $T_B$, the interaction energy $V$ is necessarily not zero and $\delta^{(U)}E_V(\rho)$ is negative.

Finally, if one could neglect the interaction energy, Equation (44) implies that energy flows from the hot to the cold system.

### 3.8. The Case $\beta_A = \beta_B$

Again assume that $U$ conserves both entropy and energy. From Equation (43) and the conservation of energy, one deduces

$$ - \delta^{(U(t))} E_V(\rho) = \delta^{(U(t))} E_A(\rho_A) + \delta^{(U(t))} E_B(\rho_B) = \frac{1}{\beta} S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}), \quad (45) $$

so that $\delta^{(U)}E_V(\rho) \leq 0$. Thus when $A$ and $B$ are initially at thermal equilibrium at the same temperature, the sum of the energies of $A$ and $B$ can only increase at the expense of the interaction energy [17].

## 4. Two Systems in Interaction With a Work Source

The problem of converting heat into work, first treated by Carnot, was at the origin of classical thermodynamics. Here, to address this issue, we explicitly introduce a work source interacting with two systems A and B, before focusing in Section 5 on the interaction of one system with a work source.

### 4.1. Hypotheses

We consider two systems $A$ and $B$ in interaction, with system $A$ coupled to a work source. We represent the action of the work source by parameters, collectively denoted by $\lambda$, so that $H_A = H_A(\lambda)$. Thus we assume that $H_B$ and $V$ are independent of $\lambda$. The action of the work source is given by an evolution of the parameters $\lambda(t)$ imposed by an external agent. The total system $A+B$ has a unitary or symplectic evolution $U(t)$ depending explicitly on time-$t$. Clearly, $U(t)$ conserves entropy but does not conserve energy, and instead one has the identity

$$ \delta^{(U(t))} E_A(\rho) + \delta^{(U(t))} E_B(\rho) + \delta^{(U(t))} E_V(\rho) + \delta^{(U(t))} W = 0, \quad (46) $$
with the following notation

$$ \delta^{(U(t))} E_B(\rho) = \mathrm{Tr} \left( (\rho^{(U)} - \rho_0) H_B \right) \quad (47) $$

$$ \delta^{(U(t))} E_V(\rho) = \mathrm{Tr} \left( (\rho^{(U)} - \rho_0) H_V \right) \quad (48) $$

$$ \delta^{(U(t))} E_A(\rho) = \mathrm{Tr} \left( \rho^{(U)} H_A(\lambda^{(U)}) - \rho_0 H_A(\lambda_0) \right) . \quad (49) $$

Here, $\lambda_0$ is the initial value of the parameter $\lambda$ and $\lambda^{(U)}$ is its final value at the end of the evolution $U$, this being an abbreviation for $U(t)$, $t$ being the final time. Note that Equation (49) extends the definition given near Equation (12). Such an extension is needed because we now allow changes in the Hamiltonian, represented by the additional variable $\lambda$. Equation (46) defines the work $\delta^{(U)}W$, which is taken to be positive if the source receives work from the system $A+B$. (Note that this is opposite to the usual convention which was implicit in the opening paragraph of this paper.)

We assume that initially A and B are in independent thermal states, but A depends on the work source parameter $\lambda$. The complete initial state at time 0 is thus

$$ \rho_0 = \rho_{A,\beta_A,\lambda_0} \otimes \rho_{B,\beta_B}, \quad (50) $$

with

$$ \rho_{A,\beta_A,\lambda_0} = \frac{e^{-\beta_A H_A(\lambda_0)}}{Z_A(\beta_A, \lambda_0)}, \quad (51) $$

$$ Z_A(\beta_A, \lambda_0) = \mathrm{Tr} \left( e^{-\beta_A H_A(\lambda_0)} \right), \quad (52) $$

and

$$ Z_A(\beta_A, \lambda_0) = \exp(-\beta_A F_A(\beta_A, \lambda_0)). \quad (53) $$

Here, $F_A(\beta_A, \lambda_0)$ denotes the equilibrium free energy for A. For a general state $\rho$ of a system with energy $H$ we define the non equilibrium free energy of the state $\rho$ at temperature $T$ to be

$$ F(\beta, \rho) = \mathrm{Tr}(\rho H) - \frac{1}{\beta} S(\rho). \quad (54) $$

In particular, for subsystem A one can define the non equilibrium free energy of the state $\rho_A^{(U)}$ at temperature $\beta_A$ to be

$$ F_A^{(U)}(\beta_A) = F_A(\beta_A, \rho_A^{(U)}) = \mathrm{Tr} \left( \rho_A^{(U)} H_A(\lambda^{(U)}) \right) - \frac{1}{\beta_A} S(\rho_A^{(U)}). \quad (55) $$

In both of the above formulas temperature is not necessarily related to the state $\rho$.

## 4.2. Identities for the Work

We next establish the following two relations

$$ \begin{aligned} \delta^{(U(t))} W = & -\delta^{(U(t))} E_V(\rho) - \mathrm{Tr} \left( \rho_A^{(U)} [H_A(\lambda^{(U)}) - H_A(\lambda_0)] \right) \\ & - \left(1 - \frac{\beta_B}{\beta_A}\right) \delta^{(U(t))} E_B(\rho) - \frac{1}{\beta_A} S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}) \end{aligned} \quad (56) $$
and

$$ \delta^{(U(t))}W = -\delta^{(U(t))}E_V(\rho) - \left(1 - \frac{\beta_B}{\beta_A}\right) \delta^{(U(t))}E_B(\rho) + \left(F_A(\beta_A, \lambda_0) - F_A^{(U)}\right) \\ -\frac{1}{\beta_A} \left(S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_B^{(U)}|\rho_{B,\beta_B})\right), \qquad (57) $$

with $F_A^{(U)}$ the non equilibrium free energy of $\rho_A^{(U)}$ calculated at the initial temperature $T_A$, namely Equation (55), $F_A^{(U)} = \text{Tr}_A(\rho_A^{(U)} H_A(\lambda^{(U)}) - T_A S(\rho_A^{(U)}))$. We will comment on these relations in Par. 4.3

**Remark 5:** Here the free energy of Equation (55) is not a state function, because it is calculated at the initial temperature of A. Note our notation: When we write $F_A(\beta_A, \lambda_0)$ this is the equilibrium free energy of the thermal state of A at temperature $\beta_A$ and external parameter $\lambda_0$. When we write $F_A^{(U)}$ we mean the non-equilibrium free energy, as defined above.

**Proof of Equation (56):** One again starts from the fundamental identities Equations (7) and (25)

$$ \delta^{(U(t))}S(\rho_B) - \beta_B \delta^{(U(t))}E_B(\rho) = -S(\rho_B^{(U)}|\rho_{B,\beta_B}), \qquad (58) $$

$$ \delta^{(U(t))}S(\rho_A) - \beta_A \text{Tr}_A \left( (\rho_A^{(U)} - \rho_{A,\beta_A,\lambda_0}) H_A(\lambda_0) \right) = -S(\rho_A^{(U)}|\rho_{A,\beta_A,\lambda_0}), \qquad (59) $$

and

$$ \delta^{(U(t))}S(\rho_A) + \delta^{(U(t))}S(\rho_B) = S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}). \qquad (60) $$

Note that Equation (59) is just Equation (7) with the substitutions $\rho \to \rho_A^{(U)}$ and $\rho_B \to \rho_{A,\beta_A,\lambda_0}$. Therefore it contains the initial $H_A(\lambda_0)$ (referring to $\rho_B \to \rho_{A,\beta_A,\lambda_0}$), not the final one. Equation (60) is likewise a rewriting of Equation (35).

Now add Equations (58) and (59) and subtract Equation (60), using the fact that

$$ S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_A^{(U)}|\rho_{A,0}) + S(\rho_B^{(U)}|\rho_{B,0}) = S(\rho^{(U)}|\rho_{A,0} \otimes \rho_{B,0}) \qquad (61) $$

(this is the same as our unnumbered equation between Equations (37) and (38)) we obtain

$$ - \beta_B \delta^{(U(t))} E_B(\rho) - \beta_A \text{Tr}_A \left( (\rho_A^{(U)} - \rho_{A,\beta_A,\lambda_0}) H_A(\lambda_0) \right) = -S(\rho^{(U)}|\rho_{A,0} \otimes \rho_{B,0}). \qquad (62) $$

Conservation of energy Equation (46) gives

$$ \delta^{(U(t))}W = -\delta^{(U(t))}E_V(\rho) - \delta^{(U(t))}E_B(\rho) - \text{Tr}_A \left( \rho_A^{(U)} (H_A(\lambda^{(U)}) - H_A(\lambda_0)) \right) \\ - \text{Tr}_A \left( (\rho_A^{(U)} - \rho_{A,\beta_A,\lambda_0}) H_A(\lambda_0) \right). \qquad (63) $$

We eliminate the second trace in the right hand side of Equation (63) using Equation (62), multiply by $T_A$ to obtain Equation (56).

**Proof of Equation (57):** In Equation (56), we replace the relative entropy term, using

$$ \begin{aligned} S(\rho_A^{(U)}|\rho_{A,\beta_A,\lambda}) &= -S(\rho_A^{(U)}) + \beta_A \text{Tr} \left( \rho_A^{(U)} H_A(\lambda_0) \right) + \log Z_A(\beta_A, \lambda_0) \\ &= -S(\rho_A^{(U)}) + \beta_A \text{Tr} \left( \rho_A^{(U)} H_A(\lambda_0) \right) - \beta_A F_A(\beta_A, \lambda_0), \end{aligned} \qquad (64) $$

and use the definition of $F_A^{(U)}$ of Equation (55)

$$ -\frac{1}{\beta_A} S(\rho_A^{(U)}|\rho_{A,\beta_A,\lambda_0}) - \text{Tr} \left( \rho_A^{(U)} (H_A(\lambda^{(U)}) - H_A(\lambda_0)) \right) = F_A(\beta_A, \lambda_0) - F_A^{(U)}. \qquad (65) $$
### 4.3. Inequalities for the Work

From the identities of Equations (56) and (57), we deduce immediately corresponding inequalities

$$ \delta^{(U(t))}W \leq -\delta^{(U(t))}E_V(\rho) - \mathrm{Tr}_A \left( \rho_A^{(U)} (H_A(\lambda^{(U)}) - H_A(\lambda_0)) \right) - \left(1 - \frac{\beta_B}{\beta_A}\right) \delta^{(U(t))}E_B(\rho) \quad (66) $$

and

$$ \delta^{(U(t))}W \leq -\delta^{(U(t))}E_V(\rho) - \left(1 - \frac{\beta_B}{\beta_A}\right) \delta^{(U(t))}E_B(\rho) + \left(F_A(\beta_A, \lambda_0) - F_A^{(U)}\right). \quad (67) $$

The interpretation of inequality (67) is straightforward. If one can neglect the interaction energy, and if $T_A = T_B$, one gets an analogue of the familiar thermodynamic inequality giving an upper bound between the work received by the work source and the variation of the free energy of A,

$$ \delta^{(U(t))}W \leq F_A(\beta_A, \lambda_0) - F_A^{(U)}. \quad (68) $$

Note that this relation is not restricted to cycles, nor to exchanges with thermal baths (which would stay in their initial thermal states).

**Remark 6:** Equation (57) contains much more information than inequalities Equations (67) and (68), since it expresses the difference between the maximum work that can be delivered by system A and the work effectively extracted from A, which is the energy dissipated in the process. It is expressed in terms of relative entropies, and it will be shown in Section 6 that it can be explicitly estimated, which yields a calculation of transport coefficients from first principles.

### 4.4. The Case of Equalities in Equations (66) and (67)

If one has equality in Equation (66), the relative entropy of Equation (56) must be equal to 0,

$$ S(\rho^{(U)} | \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}) = 0, \quad (69) $$

so $\rho^{(U)} = \rho_{A,\beta_A} \otimes \rho_{B,\beta_B}$ and the final state has come back to its initial value. If we have equality in Equation (67), both relative entropies of Equation (57) are equal to 0. In this case $\rho_B^{(U)}$ has come back to its initial value $\rho_{B,\beta_B}$ and $\delta^{(U)}E_B(\rho) = 0$. Then, one has

$$ \delta^{(U(t))}W = -\delta^{(U(t))}E_V(\rho) + F_A(\beta_A, \lambda_0) - F_A^{(U)}. \quad (70) $$

### 4.5. Case Where A is not Initially in Thermal Equilibrium.

We shall now assume that the initial state is

$$ \rho_0 = \rho_{A,0} \otimes \rho_{B,\beta_0}, \quad (71) $$

$\rho_{A,0}$ being a general state.

The following identity also holds:

$$ -\delta^{(U(t))}W = \delta^{(U(t))}F_A(\beta_B, \rho_A) + \delta^{(U(t))}E_V(\rho) + T_B \left( S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_B^{(U)}|\rho_{B,\beta_B}) \right). \quad (72) $$
This equation, true for any initial state $\rho_A$, can be found in [8,18]. Note the temperature of B appearing in the non equilibrium free energy of A

$$F_A(\beta_B, \rho_A) = E_A(\rho_A) - T_B S(\rho_A). \qquad (73)$$

If no work is performed, Equation (72) reduces to Equation (28) upon exchanging the labels A and B.

Proof: Using $S(\rho^{(U)}) = S(\rho_0)$ and the definition of the thermal state, one has

$$S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_B^{(U)}|\rho_{B,\beta_B}) = -S(\rho_{A,0}) - S(\rho_{B,\beta_B}) + S(\rho_A^{(U)}) + \beta_B E_B(\rho_B^{(U)}) + \log Z_B(\beta_B). \quad (74)$$

Then,

$$\log Z_B(\beta_B) = -\beta_B E_B(\rho_{B,\beta_B}) + S(\rho_{B,\beta_B}), \qquad (75)$$

and Equation (74) becomes

$$S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_B^{(U)}|\rho_{B,\beta_B}) = \delta^{(U(t))}S(\rho_A) + \beta_B \delta^{(U(t))}E_B(\rho_B). \qquad (76)$$

Using the conservation of energy, Equation (46), one obtains

$$\delta^{(U(t))}W = -\delta^{(U(t))}E_V(\rho) - \delta^{(U(t))}F_A(\beta_B, \rho_A) - T_B \left[ S(\rho^{(U)}|\rho_A^{(U)} \otimes \rho_B^{(U)}) + S(\rho_B^{(U)}|\rho_{B,\beta_B}) \right]. \quad (77)$$

Here

$$\delta^{(U(t))}F_A(\beta_B, \rho_A) = \delta^{(U(t))}E_A(\rho_A) - T_B \delta^{(U(t))}S(\rho_A) \qquad (78)$$

is the variation of the non equilibrium free energy of A calculated at the initial temperature $T_B$ of B and

$$\delta^{(U(t))}E_A(\rho_A) = \mathrm{Tr}_A \left( H_A(\lambda^{(U)}) \rho_A^{(U)} \right) - \mathrm{Tr}_A \left( H_A(\lambda_0) \rho_{A,0} \right). \qquad (79)$$

In particular

$$\delta^{(U(t))}W \leq -\delta^{(U(t))}E_V(\rho) - \delta^{(U(t))}F_A(\beta_B, \rho_A), \qquad (80)$$

which gives a general upper bound for the work production from heat exchanges between an arbitrary system A and a system B initially at equilibrium (not necessarily a heat bath). In this relation, equality is realized if and only if the two relative entropy terms of Equation (77) are zero, which means that

$$\rho_B^{(U)} = \rho_{B,\beta_B} \quad \text{and} \quad \rho^{(U)} = \rho_A^{(U)} \otimes \rho_{B,\beta_B}. \qquad (81)$$

**Remark 7:** By convention, a thermal bath is in a thermal state which is assumed to remain constant during the evolution. Our system B is not a thermal bath in this sense; its state varies during the evolution.

## 5. A System Coupled Only to an External Work Source

While Carnot and many others primarily considered model machines exchanging heat with several reservoirs, new thermodynamic relations have recently been announced [19] concerning exchanges of a single system with a work source. We now focus on this case.
## 5.1. Hypotheses

We consider a system coupled only to an external work source, so that the Hamiltonian of the system is $H(\lambda)$.

At time $t = 0$, the state of the system is supposed to be a thermal state $\rho_{\beta_0}(\lambda_0)$. The external observer imposes an evolution $\lambda(t)$ of the parameter $\lambda$ from $\lambda_0$ to $\lambda^{(U)}$, inducing a unitary or symplectic evolution $U$ of the whole system. The work that the external observer must perform to realize this evolution is obviously the variation of the energy of the system. With the convention of Section 4.1, we denote by $\delta^{(U(t))}W$ the work counted positive if the external source receives it from the system. We are now in a particular case of Section 4.1 when the system is $A$, there is no system $B$ and no $V$. Thus from Equation (46)

$$ \delta^{(U(t))}W = -\delta^{(U(t))}E(\rho) = \mathrm{Tr}\left(\rho_{\beta_0}(\lambda_0)H(\lambda_0) - \rho^{(U)}H(\lambda^{(U)})\right). \quad (82) $$

## 5.2. Identities for the Work

From Equations (56) and (57) we obtain immediately

$$ \delta^{(U(t))}W = -\mathrm{Tr}\left(\rho^{(U)}(H(\lambda^{(U)}) - H(\lambda_0))\right) - \frac{1}{\beta_0}S(\rho^{(U)}|\rho_{\beta_0}(\lambda_0)) \quad (83) $$

and

$$ \delta^{(U(t))}W = F(\beta_0, \lambda_0) - F^{(U)}, \quad (84) $$

with $F^{(U)}$ the non equilibrium free energy at temperature $\beta_0$.

$$ F^{(U)} = \mathrm{Tr}(\rho^{(U)} H(\lambda^{(U)})) - \frac{1}{\beta_0} S(\rho^{(U)}) \quad (85) $$

We now prove the following identity

$$ \delta^{(U(t))}W = F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}) - \frac{1}{\beta_0}S(\rho^{(U)}|\rho_{\beta_0}(\lambda^{(U)})). \quad (86) $$

This is a particular case of the result of [19].

**Proof of Equation (86):** We start from Equation (84) written as

$$ \delta^{(U(t))}W = F(\beta_0, \lambda_0) - F^{(U)} = F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}) + F(\beta_0, \lambda^{(U)}) - F^{(U)}. \quad (87) $$

Now

$$ \beta_0(F^{(U)} - F(\beta_0, \lambda^{(U)})) = -S(\rho^{(U)}) + \beta_0 \mathrm{Tr}(H(\lambda^{(U)})\rho^{(U)}) - \beta_0 F(\beta_0, \lambda^{(U)}). \quad (88) $$

But

$$ S(\rho^{(U)}|\rho_{\beta_0}(\lambda^{(U)})) = -S(\rho^{(U)}) + \beta_0 \mathrm{Tr}(H(\lambda^{(U)})\rho^{(U)}) + \log Z(\beta_0, \lambda^{(U)}), \quad (89) $$

so that comparing Equations (88) and (89), one has

$$ \beta_0(F^{(U)} - F(\beta_0, \lambda^{(U)})) = S(\rho^{(U)}|\rho_{\beta_0}(\lambda^{(U)})), \quad (90) $$

and from Equation (87) we then deduce Equation (86).

**Remark 8:** Since the transition under discussion is adiabatic, free energy is less suitable for inequalities of the form (86) than is internal energy. See Section 5.6.
## 5.3. Inequalities for the Work

### 5.3.1. From Equation (83)

From Equation (83) we deduce

$$ \delta^{(U(t))}W \leq - \mathrm{Tr} \left( \rho^{(U)} (H(\lambda^{(U)}) - H(\lambda_0)) \right), \quad (91) $$

with equality if and only if $\rho^{(U)} = \rho_{\beta_0}(\lambda_0)$, i.e., the final state $\rho^{(U)}$ is the initial state.

### 5.3.2. From Equation (86)

From Equation (86) we deduce

$$ \delta^{(U(t))}W \leq F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}) \quad (92) $$

with equality if and only if

$$ \rho^{(U)} = \rho_{\beta_0}(\lambda^{(U)}). \quad (93) $$

That is, $\rho^{(U)}$ is the thermal state at the initial temperature and final value $\lambda^{(U)}$ of $\lambda$. Note that a necessary condition for this is that the entropy of the final thermal state is the same as the entropy of the initial state.

## 5.4. Relation to the Identity of Jarzynski

Let $z$ denote a point in the phase space of the system. In this section we assume that the dynamics is classical.

We denote by $z(s|z_0)$ the classical trajectory of the phase space point at time $s$ starting from $z_0$ at time $s = 0$, for the classical evolution $U$. The external observer imposes the variation $\lambda(s)$ of $\lambda$ from $\lambda_0$ to $\lambda^{(U)} = \lambda(t)$. The identity of Jarzynski is [20]:

$$ \left\langle e^{-\beta_0 \left(H(z|z_0), \lambda^{(U)}\right) - H(z_0, \lambda)} \right\rangle_{\rho_{\beta_0}(\lambda_0)} = \frac{Z(\beta_0, \lambda^{(U)})}{Z(\beta_0, \lambda_0)} = \exp \left( -\beta_0 \left( F(\beta_0, \lambda^{(U)}) - F(\beta_0, \lambda_0) \right) \right). \quad (94) $$

Because the exponential function is strictly convex, Jensen's inequality implies that

$$ \exp \left( -\beta_0 \left\langle H(z|z_0), \lambda^{(U)} \right\rangle - H(z_0, \lambda_0) \right)_{\rho_{\beta_0}(\lambda_0)} \le \left\langle e^{-\beta_0 \left(H(z|z_0), \lambda^{(U)}\right) - H(z_0, \lambda)} \right\rangle_{\rho_{\beta_0}(\lambda_0)}, \quad (95) $$

so that using Equation (94) and taking the logarithm, one obtains

$$ \delta^{(U(t))}W \leq F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}), \quad (96) $$

which is the inequality (92).

But if the inequality (96) is an equality, we deduce $\rho^{(U)} = \rho_{\beta_0}(\lambda^{(U)})$ as in Equation (93), but we also deduce that the inequality of Jensen (95) is an equality. Because the exponential function is strictly convex, this implies that the differences

$$ H(z|z_0), \lambda^{(U)} - H(z_0, \lambda_0) = C, \quad (97) $$
where *C* is a constant independent of $z_0$ (but obviously dependent on $\lambda_0$, $\lambda^{(U)}$ and *t*); in other words, the “microscopic work” is independent of the microscopic trajectory. Although this equality would seem impossible, it turns out that identity (97) can be realized for certain systems and evolutions of $\lambda$ (see appendix A).

### 5.5. Effective Temperatures

Let $H(\lambda)$ be a Hamiltonian depending on $\lambda$ and $\rho$ a state (classical or quantum) with energy $E(\rho) = \text{Tr}(\rho H(\lambda))$. We can define two temperatures for $\rho$.

(i) The temperature $\beta_e(\rho, \lambda)$ is the temperature such that

$$E(\rho) = E(\beta_e, \lambda), \qquad (98)$$

with $E(\beta_e, \lambda) = \text{Tr}(\rho\beta_e(\lambda)H(\lambda))$. It is known that $\partial E(\beta, \lambda)/\partial\beta < 0$, so that Equation (98) defines $\beta_e$ unambiguously. The basic identity (4) shows that

$$S(\rho\beta_e(\lambda)) - S(\rho) = S(\rho|\rho\beta_e(\lambda)), \qquad (99)$$

so that

$$S(\rho\beta_e(\lambda)) \geq S(\rho), \qquad (100)$$

which is the well known fact that $\rho\beta_e(\lambda)$ maximizes the entropy among all states $\rho$ having a fixed energy. The quantity $\beta_e(\rho, \lambda)$ can be called the *effective temperature*.

(ii) There is a second temperature $\beta_a(\rho, \lambda)$ such that

$$S(\rho) = S(\beta_a, \lambda). \qquad (101)$$

In this definition, $S(\beta, \lambda)$ is the entropy of a thermal state with temperature $\beta$ and external parameter $\lambda$. We call this the *adiabatic temperature*, and by the same arguments as given above it is well-defined. From Equation (100) and Equation (101), one has

$$S(\beta_e, \lambda) \geq S(\beta_a, \lambda). \qquad (102)$$

Because

$$\left. \frac{\partial E}{\partial S(\beta, \lambda)} \right|_{\lambda \text{ fixed}} = \frac{1}{\beta} \qquad (103)$$

we deduce from Equation (102) that

$$E(\beta_a, \lambda) \leq E(\beta_e, \lambda) = E(\rho) \qquad (104)$$

and

$$\beta_a \geq \beta_e. \qquad (105)$$

Because $S$ is a strictly increasing function of $E$ (for $\lambda$ fixed), one sees that in Equation (102) or (104), one has equality if and only if $\beta_a = \beta_e$. Moreover, one has the identity

$$S(\rho|\rho\beta_a(\lambda)) - S(\rho|\rho\beta_e(\lambda)) = S(\rho|\rho\beta_e(\lambda)|\rho\beta_a(\lambda)), \qquad (106)$$

which can immediately be verified.
## 5.6. A More Precise Expression for the Work

In thermodynamics, for an adiabatic evolution, the work is related to the internal energy by $dE = -dW$, rather than to the free energy. Similarly, the work is related to the adiabatic temperature rather than to the effective energy temperature.

Given the state $\rho^{(U)}$ (corresponding to the evolution $U$, the parameter varying from $\lambda_0$ to $\lambda^{(U)}$) we can define the adiabatic temperature $\beta_a^{(U)}$ such that

$$S(\beta_a^{(U)}, \lambda^{(U)}) = S(\rho^{(U)}) = S(\beta_0, \lambda_0). \qquad (107)$$

We prove the following identity

$$\delta^{(U(t))} W = E(\beta_0, \lambda_0) - E(\beta_a^{(U)}, \lambda^{(U)}) - \frac{1}{\beta_a^{(U)}} S(\rho^{(U)} | \rho_{\beta_a^{(U)}}(\lambda^{(U)})). \qquad (108)$$

**Proof of Equation (108):** One has by definition (82)

$$-\delta^{(U(t))} W = E(\rho^{(U)}) - E(\beta_0, \lambda_0) = E(\rho^{(U)}) - E(\beta_a^{(U)}, \lambda^{(U)}) + E(\beta_a^{(U)}, \lambda^{(U)}) - E(\beta_0, \lambda_0). \qquad (109)$$

Then

$$\begin{aligned} S\left(\rho^{(U)} | \rho_{\beta_a^{(U)}}(\lambda^{(U)})\right) &= -S(\rho^{(U)}) + \beta_a^{(U)} E(\rho^{(U)}) + \log Z(\beta_a^{(U)}, \lambda^{(U)}) \\ &= \beta_a^{(U)} (E(\rho^{(U)}) - E(\beta_a^{(U)}, \lambda^{(U)})), \end{aligned} \qquad (110)$$

because $S(\rho^{(U)}) = S(\beta_a^{(U)}, \lambda^{(U)})$ by the definition (107). From this result and Equation (109) we deduce Equation (108).

As a consequence of Equation (108), we deduce the inequality

$$\delta^{(U(t))} W \le E(\beta_0, \lambda_0) - E(\beta_a^{(U)}, \lambda^{(U)}). \qquad (111)$$

In standard thermodynamics, for system thermally isolated and coupled to a work source, one has $dE = -dW$, because $\delta^{(U(t))} S = 0$ for an adiabatic (thermally isolated) process and we recover equality in Equation (111). In this situation, the inequality (92) comparing the work to the difference of free energies is not relevant, because the temperature does not remain constant.

Note that the work upper bound (111), given in terms of energy and the adiabatic temperature, is sharper than the bound given by (92), which is in terms of free energy. This is proved in the next subsection.

## 5.7. Upper Bounds on the Work Delivered by a System. Comparison of Equations (92) and (111)

We next show that using internal energy for the work inequality gives a sharper result than using the free energy. Specifically,

$$\delta^{(U(t))} W \le E(\beta_0, \lambda_0) - E(\beta_a^{(U)}, \lambda^{(U)}) \le F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}). \qquad (112)$$

**Proof of Equation (112):** We need only prove that

$$\Delta \equiv F(\beta_0, \lambda_0) - E(\beta_0, \lambda_0) - (F(\beta_0, \lambda^{(U)}) - E(\beta_a^{(U)}, \lambda^{(U)})) \ge 0. \qquad (113)$$
Using the definition of the equilibrium free energy and Equation (107) we have

$$ \Delta = -\frac{1}{\beta_0} [S(\beta_a^{(U)}, \lambda^{(U)}) - S(\beta_0, \lambda^{(U)})] + [E(\beta_a^{(U)}, \lambda^{(U)}) - E(\beta_0, \lambda^{(U)})]. \quad (114) $$

Note that in Equation (114) all terms involving $\lambda$ are evaluated at $\lambda^{(U)}$. Therefore

$$ \Delta = \int_{\beta_0}^{\beta_a^{(U)}} \left[ \frac{\partial E(\beta, \lambda^{(U)})}{\partial \beta} - \frac{1}{\beta_0} \frac{\partial S(\beta, \lambda^{(U)})}{\partial \beta} \right] d\beta. \quad (115) $$

But

$$ \frac{\partial S(\beta, \lambda)}{\partial \beta} = \frac{\partial S}{\partial E} \frac{\partial E(\beta, \lambda)}{\partial \beta} = \beta \frac{\partial E(\beta, \lambda)}{\partial \beta}. \quad (116) $$

Using Equation (116) in Equation (115), we obtain

$$ \Delta = \int_{\beta_0}^{\beta_a^{(U)}} \frac{\partial E(\beta, \lambda^{(U)})}{\partial \beta} \left[ 1 - \frac{\beta}{\beta_0} \right] d\beta. \quad (117) $$

But $\frac{\partial E(\beta, \lambda^{(U)})}{\partial \beta} < 0$, so that $\Delta \ge 0$. Note that this does not depend on which of $\beta_0$ and $\beta_a^{(U)}$ is larger.

## 5.8. The Case of Equalities in Equations (111) and (92)

### 5.8.1. Equality in Equation (111)

In this case, one has $S(\rho^{(U)}|\rho_{\beta_a^{(U)}}(\lambda^{(U)})) = 0$ in Equation (108) so

$$ \rho^{(U)} = \rho_{\beta_a^{(U)}}(\lambda^{(U)}). \quad (118) $$

In particular, $\rho^{(U)}$ is a thermal state so that

$$ \beta_e^{(U)} = \beta_a^{(U)}. \quad (119) $$

However, if one has equality in Equation (111), this does not improve the upper bound of Equation (92) for the free energy,

$$ \delta^{(U(t))}W \le F(\beta_0, \lambda_0) - F(\beta_0, \lambda^{(U)}). \quad (120) $$

In other words, the optimal bound for $\delta^{(U(t))}W$ is given by the internal energy and not the free energy (so that the internal energy in general yields a better bound than that given in [20]).

### 5.8.2. Equality in Equation (92)

From Equation (93) we deduce that

$$ \rho^{(U)} = \rho_{\beta_0}(\lambda^{(U)}), \quad (121) $$

so that $\rho^{(U)}$ is a thermal state and thus

$$ \beta_e^{(U)} = \beta_0 = \beta_a^{(U)}. \quad (122) $$

This implies that we also have equality in Equation (111)

$$ \delta^{(U(t))}W = E(\beta_0, \lambda_0) - E(\beta_0, \lambda^{(U)}). \quad (123) $$
5.9. The Case $\lambda^{(U)} = \lambda_0$

If one assumes that the final value $\lambda^{(U)}$ of $\lambda$ is equal to its initial value, we see immediately that
$\beta_a^{(U)} = \beta_0$. Indeed

$$
S(\beta_a^{(U)}, \lambda_0) = S(\rho^{(U)}) = S(\beta_0, \lambda_0), \qquad (124)
$$

so that the temperatures are equal $\beta_a^{(U)} = \beta_0$. In this case, one has from Equation (108)

$$
\delta^{(U(t))} W = - \frac{1}{\beta_0} S(\rho^{(U)} | \rho_{\beta_0}(\lambda_0)) \le 0, \quad (125)
$$

with equality if and only if

$$
\rho^{(U)} = \rho_{\beta_0}(\lambda_0), \tag{126}
$$

so that the state has returned to its initial value.

Remark 9: If the external observer imposes a variation $\lambda(t)$ of the control parameter with $\lambda(0) = \lambda_0$, $\lambda(t_1) = \lambda_1$, $\lambda(t_f) = \lambda_0$, inequality (125) says that at the end of the cycle, the observer has always lost work. In particular, the work that the external observer has put in the system in the time interval $[0, t_1]$ cannot be entirely recovered in the time interval $[t_1, t_f]$ whatever one does, except if the final state $\rho^{(U)}$ is the initial state.

Remark 10: When $\lambda^{(U)} = \lambda_0$, one can also recover Equation (125) from the identity (86). This identity reduces to

$$
\delta^{(U(t))} W = - \frac{1}{\beta_0} S(\rho^{(U)} | \rho_{\beta_0}(\lambda_0)). \qquad (127)
$$

**6. Relative Entropy, Energy Dissipation and Fourier’s Law**

In this Section we derive dissipation in the quantum context and show it to be intimately related to
the relative entropy.

6.1. *The Born Approximation*

A quantum system has a Hamiltonian

$$
H = H_0 + V. \tag{128}
$$

Let $\psi_k^{(0)}$, $E_k^{(0)}$ be the eigenstates and eigenvalues of $H_0$. In the Born approximation, the state $|\psi_n^{(0)}\rangle$
becomes at time $t$ a state $|\psi_n(t)\rangle$ with

$$
|\psi_n(t)\rangle = \sum_k a_k^{(n)}(t) e^{-iE_k^{(0)}t/\hbar} |\psi_k^{(0)}\rangle . \quad (129)
$$

The quantities $a_k^{(n)}(t) = \delta_{k,n} + \tilde{a}_k^{(n)}(t)$ satisfy

$$
i\hbar \frac{d\tilde{a}_k^{(n)}}{dt} = \sum_l V_{k,l}(t) (\delta_{l,n} + \tilde{a}_l^{(n)}), \quad (130)
$$

where

$$
V_{k,l}(t) = V_{k,l} \exp \left( \frac{i}{\hbar} (E_k^{(0)} - E_l^{(0)}) t \right). \quad (131)
$$
We assume here $V_{n,n} = 0$ for all $n$. One readily deduces that in the Born approximation

$$a_k^{(n)}(t) = \frac{V_{k,n}}{E_k^{(0)} - E_n^{(0)}} \left(1 - e^{i(E_k^{(0)} - E_n^{(0)})\frac{t}{\hbar}}\right) \quad (k \neq n) \qquad (132)$$

and by unitarity $\sum_k |a_k^{(n)}(t)|^2 = 1$, so to second order in $\tilde{a}_k^{(n)}$

$$2 \operatorname{Re} \tilde{a}_n^{(n)}(t) = - \sum_{k \neq n} |\tilde{a}_k^{(n)}(t)|^2. \qquad (133)$$

Let $\rho_0$ be an initial state diagonal in the basis $\psi_n^{(0)}$

$$\rho_0 = \sum p_{0,n} |\psi_n^{(0)}\rangle\langle\psi_n^{(0)}|. \qquad (134)$$

Then, at time $t$, the state becomes

$$\begin{aligned} \rho^{(U(t))} = & \rho_0 + \left( \sum_n p_{0,n} \sum_{l \neq n} e^{-i(E_n^{(0)} - E_l^{(0)}) \frac{t}{\hbar}} \tilde{a}_l^{(n)*} |\psi_n^{(0)}\rangle\langle\psi_l^{(0)}| + \text{c.c.} \right) \\ & + \sum_n p_{0,n} (\tilde{a}_n^{(n)}(t) + \tilde{a}_n^{(n)*}(t)) |\psi_n^{(0)}\rangle\langle\psi_n^{(0)}| \\ & + \sum_n p_{0,n} \sum_{k,l \neq n} \tilde{a}_k^{(n)}(t) \tilde{a}_l^{(n)}(t)^* e^{-i(E_k^{(0)} - E_l^{(0)}) \frac{t}{\hbar}} |\psi_k^{(0)}\rangle\langle\psi_l^{(0)}| + \dots \end{aligned} \qquad (135)$$

If $L$ is a Hermitian operator diagonal in the basis $\psi_n^{(0)}$ with eigenvalues $\lambda_n$, using Equation (135) one obtains in the Born approximation

$$\operatorname{Tr} (L(\rho^{(U(t))} - \rho_0)) = \frac{1}{2} \sum_{k \neq n} |\tilde{a}_k^{(n)}(t)|^2 (\lambda_k - \lambda_n) (p_{0,n} - p_{0,k}). \qquad (136)$$

## 6.2. Two Interacting Systems

We consider two quantum systems A, B with Hamiltonians $H_A, H_B$ respectively, interacting. Denote by $V = V_{A,B}$ the interaction energy and

$$H = H_A + H_B + V. \qquad (137)$$

We call $|\psi_{A,k}^{(0)}\rangle, E_{A,k}^{(0)}$ (resp. $|\psi_{B,l}^{(0)}\rangle, E_{B,l}^{(0)}$) the eigenstates and eigenvalues of $H_A$ (resp. $H_B$), and we apply the Born approximation to $H$, with $H_0 = H_A + H_B$. The non perturbed Hamiltonian $H_0$ has eigenstates $|\psi_{A,k}^{(0)}\rangle|\psi_{B,l}^{(0)}\rangle$ with eigenvalues $E_{A,k}^{(0)} + E_{B,l}^{(0)}$.

We assume that at time $t = 0$, the state of the system $A+B$ is $\rho_0 = \rho_A \otimes \rho_B$ with

$$\begin{aligned} \rho_A &= \sum p_{A,k} |\psi_{A,k}^{(0)}\rangle\langle\psi_{A,k}^{(0)}| \\ \rho_B &= \sum p_{B,l} |\psi_{B,l}^{(0)}\rangle\langle\psi_{B,l}^{(0)}|, \end{aligned} \qquad (138)$$

so that they are diagonal in the eigenbasis of $H_A$ and $H_B$ and therefore commute with $H_A + H_B$. At time $t$, the initial state $\rho_0 = \rho_A \otimes \rho_B$ evolves to $\rho(t)$. Then

$$S(\rho(t)|\rho_A \otimes \rho_B) = \operatorname{Tr}(\rho(t) \log \rho(t)) - \operatorname{Tr}(\rho(t) \log(\rho_A \otimes \rho_B)). \qquad (139)$$
But $S(\rho(t)) = S(\rho_0)$ by unitarity of the evolution, so that

$$S(\rho(t)|\rho_A \otimes \rho_B) = - \mathrm{Tr} \left[ (\rho(t) - \rho_A \otimes \rho_B) (\log \rho_A + \log \rho_B) \right]. \quad (140)$$

This is of the form of Equation (136) with

$$L = -(\log \rho_A + \log \rho_B). \quad (141)$$

$L$ has eigenvectors $|\psi_{A,k}^{(0)}\rangle|\psi_{B,l}^{(0)}\rangle$ with eigenvalues $\log p_{A,k} + \log p_{B,l}$; $\rho_A \otimes \rho_B$ has the same eigenvectors with eigenvalues $p_{A,k} + p_{B,l}$. Applying Equation (136), one obtains in the Born approximation

$$S(\rho(t)|\rho_A \otimes \rho_B) = \frac{1}{2} \sum_{(k,l) \neq (n,m)} |\tilde{a}_{(k,l)}^{(n,m)}(t)|^2 (p_{A,n}p_{B,m} - p_{A,k}p_{B,l}) \log \frac{p_{A,n}p_{B,m}}{p_{A,k}p_{B,l}}. \quad (142)$$

Notice that the quantity in the right hand side is automatically non-negative. Here, we have

$$|\tilde{a}_{(k,l)}^{(n,m)}(t)|^2 = |V_{(k,l)}^{(n,m)}|^2 \frac{\sin^2\left(\frac{E_{A,k}^{(0)} + E_{B,l}^{(0)} - E_{A,n}^{(0)} - E_{B,m}^{(0)}}{2\hbar} t\right)}{\left(\frac{E_{A,k}^{(0)} + E_{B,l}^{(0)} - E_{A,n}^{(0)} - E_{B,m}^{(0)}}{2}\right)^2}, \quad (143)$$

with $V_{(k,l)}^{(n,m)} = \langle \psi_{A,n}^{(0)} \otimes \psi_{B,m}^{(0)} | V | \psi_{A,k}^{(0)} \otimes \psi_{B,l}^{(0)} \rangle$.

We also deduce from this result that in this approximation $S(\rho(t)|\rho_A \otimes \rho_B) = 0$ if and only if $V = 0$ (recall that the diagonal elements of $V$ are 0).

### 6.3. The Case Where Both Initial States are Thermal

Assume that at time $t = 0$, $\rho_A = \rho_{A,\beta_A}$ and $\rho_B = \rho_{B,\beta_B}$ are the thermal states of A and B respectively. From Equation (37) one has

$$\beta_A \delta^{(U(t))} E_A(\rho_A) + \beta_B \delta^{(U(t))} E_B(\rho_B) = S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B}). \quad (144)$$

Moreover, from conservation of energy

$$\delta^{(U(t))} E_A(\rho_A) + \delta^{(U(t))} E_B(\rho_B) + \delta^{(U(t))} E_V(\rho) = 0, \quad (145)$$

so that eliminating $\delta^{(U(t))} E_B(\rho_B)$, one obtains

$$(\beta_A - \beta_B) \delta^{(U(t))} E_A(\rho_A) = \beta_B \delta^{(U(t))} E_V(\rho) + S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B}). \quad (146)$$

We now estimate both terms on the right hand side of Equation (146).

#### 6.3.1. Estimate of the Relative Entropy

From Equation (142) we deduce

$$S(\rho(t)|\rho_A \otimes \rho_B) = \frac{1}{2Z_A Z_B} \sum_{(k,l) \neq (n,m)} |\tilde{a}_{(k,l)}^{(n,m)}(t)|^2 e^{-(\beta_A E_{A,n} + \beta_B E_{B,m})} \\ \times (1 - e^{-(\beta_A(E_{A,k} - E_{A,n}) + \beta_B(E_{B,l} - E_{B,m})))} \\ \times (\beta_A(E_{A,k} - E_{A,n}) + \beta_B(E_{B,l} - E_{B,m})). \quad (147)$$
Moreover, when $t \to \infty$, as in the usual Born approximation, Equation (143) shows that

$$|\tilde{a}_{(k,l)}^{(n,m)}(t)|^2 \approx \frac{2\pi}{\hbar} |V_{(k,l)}^{(n,m)}|^2 \delta(E_{A,k}^{(0)} + E_{B,l}^{(0)} - E_{A,n}^{(0)} - E_{B,m}^{(0)}) t. \quad (148)$$

Thus if $f_A$ and $f_B$ denote the density of states for A and B, we obtain from Equation (147)

$$S(\rho(t)|\rho_A \otimes \rho_B) = \frac{2\pi}{2\hbar Z_A Z_B} (\beta_A - \beta_B) \int dE_A dE'_A dE_B dE'_B f_A(E_A) f_A(E'_A) f_B(E_B) f_B(E'_B) \\ \times e^{-(\beta_A E'_A + \beta_B E'_B)} |V_{(E_A, E'_B)}^{(E'_A, E'_B)}|^2 \delta(E_A + E_B - E'_A - E'_B) (E_A - E'_A) \\ \times (1 - e^{-(\beta_A - \beta_B)(E_A - E'_A)}) t. \quad (149)$$

### 6.3.2. Estimate of the Interaction Energy

Because $\delta^{(U(t))}V(\rho) = -\delta^{(U(t))}E_A(\rho) - \delta^{(U(t))}E_B(\rho)$, one has

$$\delta^{(U(t))}V_{A,B}(\rho) = \mathrm{Tr}(-(H_A + H_B)(\rho(t) - \rho_{A,\beta_A} \otimes \rho_{B,\beta_A})). \quad (150)$$

This is of the form of Equation (136) with $L = -H_A - H_B$, and so

$$\beta_B \delta^{(U(t))} V_{A,B}(\rho) = -\frac{\beta_B}{2Z_A Z_B} \sum_{(k,l) \neq (n,m)} |\tilde{a}_{(k,l)}^{(n,m)}(t)|^2 (E_{A,k} + E_{B,l} - E_{A,n} - E_{B,m}) \\ \times e^{-(\beta_A E'_A + \beta_B E'_B)} (1 - e^{-(\beta_A(E_{A,k} - E_{A,n}) + \beta_B(E_{B,l} - E_{B,m})))}) t. \quad (151)$$

Up to a sign, this expression is formally identical to the expression Equation (142), except that the difference of energies ($E_{A,k} + E_{B,l} - E_{A,n} - E_{B,m}$) replaces the quantity $\beta_A(E_{A,k} - E_{A,n}) + \beta_B(E_{B,l} - E_{B,m})$. As a consequence $E_{A,k} + E_{B,l} - E_{A,n} - E_{B,m}$ partially cancels the denominator of $|\tilde{a}_{(k,l)}^{(n,m)}(t)|^2$ and one sees that $\beta_B \delta^{(U(t))} E_V(\rho)$ is negligible when $t \to \infty$.

Then from Equations (146) and (149), one sees that

$$\delta^{(U(t))} E_A(\rho_A) \approx \frac{S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B})}{\beta_A - \beta_B} \approx (\beta_A - \beta_B)Kt. \quad (152)$$

In Equation (152) K is the positive constant

$$K = \frac{2\pi}{2\hbar Z_A Z_B} \int dE_A dE'_A dE_B dE'_B \varphi(E_A, E'_A, E_B, E'_B), \quad (153)$$

with

$$\varphi = f_A(E_A) f_A(E'_A) f_B(E_B) f_B(E'_B) e^{-(\beta_A E'_A + \beta_B E'_B)} |V_{(E_A, E'_B)}^{(E'_A, E'_B)}|^2 \delta(E_A + E_B - E'_A - E'_B) \\ \times \frac{E_A - E'_A}{\beta_A - \beta_B} (1 - e^{-(\beta_A - \beta_B)(E_A - E'_A)}). \quad (154)$$

It is obvious that $\varphi \ge 0$. Note that K does not vanish for $\beta_A$ close to $\beta_B$.

The expression (152) is a form of Fourier's law for heat transport from B to A, $(\beta_A - \beta_B)K$ being the rate of dissipation. In this case, one sees that the significance of the relative entropy $S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B})$ is that of a transport coefficient, here the transport of energy from one part of a system to another part.
## 7. Coarse Grained States

Coarse-graining is omnipresent in macroscopic and mesoscopic physics, since microscopic variables are often not what is observed. In general coarse-graining represents a loss of information, hence an increase in entropy. In this section we consider a variety of coarse-graining procedures, and consistent with our work in this article, relative entropy plays a significant role both in the definition of coarse-graining and in the measures of entropy increase.

### 7.1. Definition

Let $\rho$ and $\rho'$ be two states of the same system (classical or quantum). We say that $\rho'$ is obtained from $\rho$ by a coarse graining operation if

$$ \mathrm{Tr}((\rho - \rho') \log \rho') = 0. \qquad (155) $$

The idea is that the information associated with $\rho'$ (namely $\log \rho'$) is the same whether one averages with $\rho'$ or with the more detailed distribution $\rho$. Using the basic identity, Equation (4), we can say that $\rho'$ is obtained from $\rho$ by a coarse graining operation, if and only if

$$ S(\rho') - S(\rho) = S(\rho|\rho'). \qquad (156) $$

In particular, $S(\rho') \ge S(\rho)$, so that the entropy increases by coarse-graining. (See the comment after Equation (3).)

A coarse-graining mapping is a mapping $\Gamma$ which associates to any state $\rho$ (or to some states of a given class), a coarse grained state $\rho' = \Gamma(\rho)$.

### 7.2. Examples of Coarse-graining Mappings

#### Example 1: Maximum entropy.

Let $A_1, \dots, A_n$ be observables of the system, so they are either functions in the phase space or hermitian operators on the Hilbert space of the system. We consider the class of states $\rho$ such that

$$ \mathrm{Tr}(A_i \rho) < \infty, \quad i = 1, \dots, n. \qquad (157) $$

One can then consider the state $\rho'$ such that $\rho'$ has maximal entropy given the relation

$$ \mathrm{Tr}(A_i \rho') = \mathrm{Tr}(A_i \rho). \qquad (158) $$

It is immediately seen that

$$ \rho' = C \exp \left( \sum_{i=1}^{n} \alpha_i A_i \right), \qquad (159) $$

where $C$ is a normalization constant and $\alpha_i$ are the “conjugate parameters”, (provided $\rho'$ is normalizable). The mapping $\Gamma: \rho \to \rho'$ is indeed a coarse grain mapping in the sense of the previous definition, because by Equation (159)

$$ \mathrm{Tr}((\rho - \rho') \log \rho') = \sum_{i=1}^{n} \mathrm{Tr}((\rho - \rho') A_i) = 0. \qquad (160) $$
In particular, one has Equation (156).

The case of the thermal state is the best known, where one takes $A = H$, the Hamiltonian of the system.

*Example 2:* Naive coarse graining; the observables as characteristic functions.

(i) Classical case: Let $Z$ be the phase space of the system and $\{Z_i\}$ a finite partition of $Z$ ($Z = \bigcup_{i=1}^n Z_i$ and $Z_i \cap Z_j = \emptyset$ for $i \neq j$). We choose $A_i = \chi_{Z_i}$ (i.e. the characteristic function of $Z_i$). This is a particular case of example 1 and if $\rho$ is a state

$$ \rho' = C \exp \left( \sum_{i=1}^{n} \alpha_i \chi_{Z_i} \right). \qquad (161) $$

Using the condition (156), namely,

$$ \int_{Z_i} \rho' dz = \int_{Z_i} \rho dz, \qquad (162) $$

one can deduce from Equations (161) and (162)

$$ \rho'|_{Z_i} = \frac{1}{\mathrm{Vol}(Z_i)} \int_{Z_i} \rho \, dz \quad \text{or} \quad \rho' = \sum_{i=1}^{n} (\rho'|_{Z_i}) \chi_{Z_i}. \qquad (163) $$

This equation implies that $\rho'$ is normalized $\int \rho' dz = 1$. We recover the usual coarse graining.

(ii) Quantum case: Let $\mathcal{H}$ be the Hilbert space of the system and $P_i$ a resolution of the identity by orthogonal projectors

$$ \mathrm{Id} = \sum P_i \quad \text{and} \quad P_i P_j = P_i \delta_{i,j}. \qquad (164) $$

Then the analogue of Equation (163) is

$$ \rho' = \sum_{i=1}^{n} \frac{\mathrm{Tr}(\rho P_i)}{\mathrm{dim} P_i(\mathcal{H})} P_i. \qquad (165) $$

*Example 3:* Coarse graining by marginals.

(i) Classical case: We assume that the system consists of several parts, and that its phase space is a Cartesian product, $Z = \prod_{i=1}^n Z_i$, corresponding to various subsystems with phase space $Z_i$. If $\rho$ is a state on $Z$, we denote by $\rho_i$ its marginal probability distribution on $Z_i$, so

$$ \rho_i(z_i) = \int \dots \int \rho(\zeta_1, \dots, \zeta_{i-1}, z_i, \zeta_{i+1}, \dots, \zeta_n) \prod_{j \neq i} d\zeta_j. \qquad (166) $$

Let $\Gamma$ be the mapping that associates the product of its marginals to $\rho(z)$

$$ \Gamma(\rho)(z_1, \dots, z_n) = \prod_{i=1}^{n} \rho_i(z_i). \qquad (167) $$

Then the condition (155) is satisfied. It is easy to see that $\prod_{i=1}^n \rho_i(z_i)$ is the state $\rho'$ that maximizes the entropy among all the states $\rho''$ such that $\rho''_i = \rho_i$ for any $i$.
(ii) Quantum case. The Hilbert space of the system is

$$ \mathcal{H} = \bigotimes_{i=1}^{n} \mathcal{H}_i, \qquad (168) $$

where the $\mathcal{H}_i$ are the Hilbert spaces of the subsystems. If $\rho$ is a state, then its marginal state on $\mathcal{H}_i$ is the partial trace on the Hilbert space $\mathcal{K}_i$, which is the tensor product of the Hilbert spaces $\mathcal{H}_j$ for $j$ different from $i$

$$ \rho_i = \mathrm{Tr}_{\mathcal{K}_i} \rho \qquad (169) $$

and the mapping $\Gamma$,

$$ \Gamma(\rho) = \bigotimes_{i=1}^{n} \rho_i, \qquad (170) $$

is a coarse grained mapping. $\Gamma(\rho)$ is again the state $\rho'$ which maximizes the entropy among all states $\rho''$ such that $\rho''_i = \rho_i$ for all $i$.

**Example 4:** Decomposition of $Z$.

If $Z = \bigcup_{i=1}^n Z_i$, but the $Z_i$ do not form a partition of $Z$ (they can have intersections of non-zero measure), one can still apply Example 1 to $A_i = \chi_{Z_i}$ and obtain

$$ \rho' = C \exp \left( \sum \alpha_i \chi_{Z_i} \right). \qquad (171) $$

But now Equation (163) is no longer valid because, for given $z$, there will be in general several $i$ with $z \in Z_i$.

## 7.3. Coarse Graining and Relative Entropy

(i) The case of the naive coarse-graining is distinguished among all types of coarse-graining by the following property. Let $Z = \bigcup_{i=1}^n Z_i$ a partition of the phase space and $p, q$ two probability distributions on $Z$. Let $\bar{p} = \Gamma p$ and $\bar{q} = \Gamma q$ be the coarse grained states of $p$ and $q$ associated to this partition. Then one has

$$ S(\bar{p}|\bar{q}) \le S(p|q). \qquad (172) $$

*Proof:* call $p_i = \int_{Z_i} p \, dz$ and $q_i = \int_{Z_i} q \, dz$. We have, using the definition of $\bar{p}$ and $\bar{q}$

$$ S(\bar{p}|\bar{q}) = \sum_{i=1}^{n} p_i \log \frac{p_i}{q_i}. \qquad (173) $$

Now

$$ \begin{aligned} \sum_i p_i \log \frac{p_i}{q_i} &= \sum_i \left( \int_{Z_i} p(z) \, dz \right) \log \frac{\int_{Z_i} p(z) \, dz}{\int_{Z_i} q(z') \, dz'} \\ &= \sum_i \left( \int_{Z_i} q(z') \, dz' \right) \left( \int_{Z_i} \frac{p(z)}{q(z)} \frac{q(z)}{\int_{Z_i} q(z') \, dz'} \, dz \right) \\ &\quad \times \log \left( \int_{Z_i} \frac{p(z)}{q(z)} \frac{q(z)}{\int_{Z_i} q(z') \, dz'} \, dz \right). \end{aligned} \qquad (174) $$
But $\int_{Z_i} \frac{q(z)}{\int_{Z_i} q(z') dz'} dz = 1$. We use the fact that the function $x \log x$ is convex, so that for each $i$

$$ \left( \int_{Z_i} \frac{p(z)}{q(z)} \frac{q(z)}{\int_{Z_i} q(z') dz'} dz \right) \log \left( \int_{Z_i} \frac{p(z)}{q(z)} \frac{q(z)}{\int_{Z_i} q(z') dz'} dz \right) \\ \le \int_{Z_i} \frac{q(z)}{\int_{Z_i} q(z') dz'} \left( \frac{p(z)}{q(z)} \log \frac{p(z)}{q(z)} \right) dz. \quad (175) $$

Therefore from Equation (175)

$$ \sum_i p_i \log \frac{p_i}{q_i} \le \sum_i \int_{Z_i} p(z) \log \frac{p(z)}{q(z)} dz = S(p|q). \quad (176) $$

(ii) For the coarse-graining associated to subsystems one has $Z = \prod_{i=1}^n Z_i$ and if $p, q$ are states on $Z$, the coarse grained states are $\tilde{p} = \bigotimes_{i=1}^n p_i$, $\tilde{q} = \bigotimes_{i=1}^n q_i$ and we deduce immediately that

$$ S(\tilde{p}|\tilde{q}) = \sum_{i=1}^{n} S(p_i|q_i). \quad (177) $$

Consider the case $i=1$, and call $z = (z_1, z')$ with $z' = (z_2, \dots, z_n)$ and call $Z' = Z_2 \times \dots \times Z_n$. Then

$$ \begin{aligned} S(p_1|q_1) &= \int_{Z_1} dz_1 \left( \int_{Z'} p(z_1, z') dz' \right) \log \frac{\int_{Z'} p(z_1, z') dz'}{\int_{Z'} q(z_1, z') dz'} \\ &= \int_{Z_1} dz_1 \left( \int_{Z'} q(z_1, z'') dz'' \right) \left( \int_{Z'} dz' \frac{p(z_1, z')}{q(z_1, z')} \frac{q(z_1, z')}{\int_{Z'} q(z_1, z'') dz''} \right) \\ &\qquad \times \log \int_{Z'} dz' \frac{p(z_1, z')}{q(z_1, z')} \frac{q(z_1, z')}{\int_{Z'} q(z_1, z'') dz''} . \end{aligned} \quad (178) $$

Now $\int_{Z'} dz' \frac{q(z_1, z')}{\int_{Z'} q(z_1, z'') dz''} = 1$. As in Equation (175), we use the convexity of $x \log x$ and deduce that

$$ S(p_1|q_1) \le \int_{Z_1} dz_1 \int_{Z'} dz' p(z_1, z') \log \frac{p(z_1, z')}{q(z_1, z')} = S(p|q). \quad (179) $$

From Equation (177) we deduce that for the coarse graining mapping associated to the division of $Z = \prod_{i=1}^n Z_i$ in $n$ subsystems, one has

$$ S(\tilde{p}|\tilde{q}) \le nS(p|q). \quad (180) $$

**Remark 11:** The upper bound of Equation (180) cannot be improved. Indeed consider the case where: $p(z_1, \dots, z_n) = p_1(z_1)\delta(z_1 - z_2) \dots \delta(z_{n-1} - z_n)$, $q(z_1, \dots, z_n) = q_1(z_1)\delta(z_1 - z_2) \dots \delta(z_{n-1} - z_n)$. Then $p_i = p_1$ and $q_i = q_1$, but $S(p|q) = S(p_1|q_1)$ and $S(\tilde{p}|\tilde{q}) = nS(p_1|q_1)$.

(iii) Thermal coarse graining.

Let $Z$ be a phase space, and $p$ and $q$ two probability distributions on $Z$, $H(z)$ a function of $z \in Z$.
Let $\tilde{p}$ and $\tilde{q}$ be the thermal coarse grained probability distributions of $p$ and $q$, respectively, with respect to $H$. So

$$ \tilde{p}(z) = \frac{1}{Z(\beta(p))} \exp(-\beta(p)H(z)). \qquad (181) $$

where $\beta(p)$ is the effective temperature of $p$, i.e., $\langle H \rangle_p = \langle H \rangle_{\tilde{p}}$.

Assuming that $p - q$ is small, an obvious bound, after straightforward calculations (expanding to second order in $p - q$), is

$$ S(\tilde{p}|\tilde{q}) \le \frac{\langle H^2 \rangle_p}{\langle H^2 \rangle_{\tilde{p}} - (\langle H \rangle_p)^2} S(p|q). \qquad (182) $$

This bound is surely not optimal, because if $p$ and $q$ are already thermal states, $S(\tilde{p}|\tilde{q}) = S(p|q)$. Note though that even without the hypotheses on $p$ and $q$, $S(\tilde{p}|\tilde{q}) \le S(p|q)$.

# 8. Conclusions

The results in this article are used to obtain upper bounds for entropy production or energy variation in various situations of thermodynamic interest, with many such results either new or sharper than similar known bounds. Furthermore, the energy dissipated in these processes is expressed in terms of relative entropies, which not only gives a general microscopic interpretation of dissipation, but also, in relevant examples, leads to an explicit, first principles, evaluation of dissipation terms, analogous to the Fourier law.

Although relative entropy has made appearances in many contexts, especially with respect to information theory, our results on a generalized Fourier heat law relates it in a direct way to the notion of dissipation as understood in physics.

# Acknowledgments

We (B.G. and L.S.) are grateful to the visitor program of the Max Planck Institute for the Physics of Complex Systems, Dresden, for its hospitality during the time that much of this work was performed. We also thank one of the referees for an extremely careful reading of the manuscript, leading we believe to considerable improvement in presentation and content.

# Author Contributions

The authors contributed equally to the presented mathematical framework and the writing of the paper.

# Appendixes

## A. An Example of Trajectory-Independent Microscopic Work.

We exhibit a Hamiltonian $H(z, \lambda)$ and an evolution $\lambda(t)$ of the external parameter such that

$$ H(z(t_f|z_0), \lambda(t_f)) - H(z_0, \lambda_0) = C, \qquad (A1) $$

with $C$ independent of $z_0$.
Take the harmonic oscillator

$$H(x, p, \lambda) = \frac{p^2}{2} + \frac{\omega^2 x^2}{2} - \omega^2 \lambda x. \qquad (A2)$$

Call

$$\begin{aligned}
\bar{x}(t|x_0, p_0) &= x_0 \cos \omega t + p_0 \frac{\sin \omega t}{\omega} \\
\bar{p}(t|x_0, p_0) &= -\omega x_0 \sin \omega t + p_0 \cos \omega t
\end{aligned} \qquad (A3)$$

the solution with $\lambda = 0$.

For $\lambda(s)$ a function of time $s$, the solutions of the Hamiltonian equations starting from $(x_0, p_0)$ at $s = 0$ are

$$x(t|x_0, p_0) = \bar{x}(t) + \omega \int_0^t \lambda(s) \sin(\omega(t-s)) \, ds \qquad (A4)$$

$$p(t|x_0, p_0) = \bar{p}(t) + \omega^2 \int_0^t \lambda(s) \cos(\omega(t-s)) \, ds \qquad (A5)$$

with $\frac{dx}{dt} = p$ and $\frac{dp}{dt} = -\omega^2 x + \omega^2 \lambda(t)$. Assume that $\lambda_0 = 0$. Define $\Lambda_c(t) \equiv \int_0^t \lambda(s) \cos \omega(t-s) \, ds$ and $\Lambda_s(t) \equiv \int_0^t \lambda(s) \sin \omega(t-s) \, ds$. Then

$$\begin{aligned}
H(x(t), p(t), \lambda(t)) \quad &-H(x_0, p_0, 0) = -\omega^2 \lambda(t) \bar{x}(t) + \omega^2 \bar{p}(t) \Lambda_c(t) \\
&+\omega^3 \bar{x}(t) \Lambda_s(t) + \frac{1}{2} \omega^4 (\Lambda_c(t))^2 \\
&+\frac{1}{2} \omega^4 (\Lambda_s(t))^2 - \omega^3 \lambda(t) \Lambda_s(t).
\end{aligned} \qquad (A6)$$

We can impose a condition on $t$ such that this quantity does not depend on $x_0$ and $p_0$, namely

$$\begin{aligned}
\Lambda_c(t) &= \int_0^t \lambda(s) \cos \omega(t-s) \, ds = 0 \\
-\lambda(t) + \omega \Lambda_s(t) &= -\lambda(t) + \omega \int_0^t \lambda(s) \sin \omega(t-s) \, ds = 0.
\end{aligned} \qquad (A7)$$

Then using these two equalities, one has

$$H(x(t), p(t), \lambda(t)) - H(x_0, p_0, 0) = -\frac{\omega^2}{2} (\lambda(t))^2. \qquad (A8)$$

Thus if $\lambda(t) \neq 0$, we can arrange that the microscopic work is independent of the initial condition and is non zero.

## B. An Exactly Solvable Model

The system $A + B$ is formed of two two-levels atoms. The Hamiltonians of $A$ and $B$ are

$$H_A = \begin{pmatrix} 0 & 0 \\ 0 & E_A \end{pmatrix} \quad \text{and} \quad H_B = \begin{pmatrix} 0 & 0 \\ 0 & E_B \end{pmatrix}, \qquad (A9)$$
with eigenstates $|0_A\rangle, |+_A\rangle, |0_B\rangle, |+_B\rangle$, so that the total Hamiltonian is in the basis $|0_A, 0_B\rangle, |+_A, 0_B\rangle, |0_A, +_B\rangle, |+_A, +_B\rangle$:

$$H = \begin{pmatrix} 0 & 0 & 0 & w \\ 0 & E_A & 0 & 0 \\ 0 & 0 & E_B & 0 \\ w^* & 0 & 0 & E_A + E_B \end{pmatrix} \qquad (A10)$$

where $w$ is the interaction energy.

Calling $E_0 = E_A + E_B$, the eigenvalues of $H$ are

$$\lambda_{\pm} = \frac{E_0 \pm \sqrt{E_0^2 + 4|w|^2}}{2}. \qquad (A11)$$

as well as $E_A$ and $E_B$. The eigenstates of $E_A$ and $E_B$ are $|+_A, 0_B\rangle, |0_A, +_B\rangle$, and the eigenstates of $\lambda_{\pm}$ are

$$|\varphi_{\pm}\rangle = \frac{1}{N_{\pm}} (w|0_A, 0_B\rangle + \lambda_{\pm}|+_A, +_B\rangle), \qquad (A12)$$

so that

$$|0_A, 0_B\rangle = \frac{N_+ N_-}{w(\lambda_- - \lambda_+)} \left( \frac{\lambda_-}{N_-} |\varphi_+\rangle - \frac{\lambda_+}{N_+} |\varphi_-\rangle \right) \qquad (A13)$$

$$|+_A +_B\rangle = \frac{N_+ N_-}{\lambda_+ - \lambda_-} \left( \frac{1}{N_-} |\varphi_+\rangle - \frac{1}{N_+} |\varphi_-\rangle \right). \qquad (A14)$$

Here $N_{\pm} = \sqrt{|w|^2 + |\lambda_{\pm}|^2}$ is the normalization factor.

The initial state is $\rho_{A,\beta_A} \otimes \rho_{B,\beta_B}$:

$$\rho_{A,\beta_A} \otimes \rho_{B,\beta_B} = \frac{1}{Z_A Z_B} \left( |0_A, 0_B\rangle \langle 0_A, 0_B| + e^{-\beta_A E_A} |+_A, 0_B\rangle \langle+_A, 0_B| \right. \\ \left. + e^{-\beta_B E_B} |0_A, +_B\rangle \langle 0_A, +_B| + e^{-\beta_A E_A - \beta_B E_B} |+_A +_B\rangle \langle+_A +_B| \right) \qquad (A15)$$

Using these formulas one can compute

$$\rho(t) = e^{-iHt} \rho_{A,\beta_A} \otimes \rho_{B,\beta_B} e^{iHt} = U(t) \rho_{A,\beta_A} \otimes \rho_{B,\beta_B} U(t)^+ \qquad (A16)$$

and verify that

$$\mathrm{Tr} (H_A \rho(t)) = \frac{E_A}{Z_A Z_B} \left( e^{-\beta_A E_A} + e^{-\beta_A E_A - \beta_B E_B} \frac{|\lambda_+ e^{-i\lambda_+ t} - \lambda_- e^{-i\lambda_- t}|^2}{(\lambda_- - \lambda_+)^2} |w|^2 \frac{|e^{-i\lambda_+ t} - e^{-i\lambda_- t}|^2}{(\lambda_- - \lambda_+)^2} \right). \qquad (A17)$$

Then

$$S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B}) = \mathrm{Tr}((\beta_A H_A + \beta_B H_B)(\rho(t) - \rho_{A,\beta_A} \otimes \rho_{B,\beta_B})) \qquad (A18)$$

$$\delta E_V(\rho) = - \mathrm{Tr}((H_A + H_B)(\rho(t) - \rho_{A,\beta_A} \otimes \rho_{B,\beta_B})) \qquad (A19)$$

and

$$S(\rho(t)|\rho_{A,\beta_A} \otimes \rho_{B,\beta_B}) = \frac{\beta_A E_A + \beta_B E_B}{Z_A Z_B} \left(1 - e^{-\beta_A E_A - \beta_B E_B}\right) \frac{|w|^2 \sin^2 \left(\frac{(\lambda_+ - \lambda_-)t}{2}\right)}{\left(\frac{\lambda_+ - \lambda_-}{2}\right)^2} \qquad (A20)$$
$$ \delta^{(U(t))} E_V(\rho) = -\frac{E_A + E_B}{Z_A Z_B} (1 - e^{-\beta_A E_A - \beta_A E_B}) |w|^2 \frac{\sin^2((\lambda_+ - \lambda_-)t)}{(\lambda_- - \lambda_+)^2}. \quad (A21) $$

Using Equation (146), one obtains

$$ \delta^{(U(t))} E_A(\rho) = \frac{E_A}{Z_A Z_B} (1 - e^{-\beta_A E_A - \beta_A E_B}) |w|^2 \frac{\sin^2((\lambda_+ - \lambda_-)t)}{(\lambda_- - \lambda_+)^2}. \quad (A22) $$

Here these quantities are periodic functions of period $\frac{2\pi}{\lambda_+-\lambda_-} = \frac{2\pi}{\sqrt{E_0^2+4|w|^2}}$. Near resonance, where $\lambda_+ \approx \lambda_-$, $w \approx 0$, $E_0 = E_A + E_B \approx 0$ and we recover that $\delta^{(U(t))} E_A(\rho) \approx K(\beta_A - \beta_B)t$ from Equation (A22).

### C. Example: Forced Harmonic Oscillator

We take the Hamiltonian

$$ H = -\frac{1}{2} \frac{d^2}{dx^2} + \frac{\omega^2 x^2}{2} + \lambda(t)x = H_0 + \lambda(t)x, \quad (A23) $$

with the condition $\lambda(0) = 0$. The classical action is

$$ S(x, t|x') = \frac{\omega}{2\sin\omega t} \left( (x^2 + x'^2) \cos\omega t - 2xx' - \frac{2x}{\omega} \int_0^t \lambda(s) \sin\omega s \, ds - \frac{2x'}{\omega} \int_0^t \lambda(s) \sin\omega(t-s) \, ds + C(t) \right), \quad (A24) $$

where $C(t)$ does not depend on $x$ or $x'$. The quantum propagator is

$$ G(x, t|x', 0) \approx \exp(iS(x, t|x', 0)), \quad (A25) $$

where “$\simeq$” indicates that we have not written the normalization factor. This factor does not depend on $x$ or $x'$ and is at the moment unimportant. The thermal state for $\lambda = 0$ is

$$ \rho_{\beta} \approx \exp \left( - \frac{i\omega}{2\sin(i\omega\beta)} \left( (y^2 + y'^2) \cos i\omega\beta - 2yy' \right) \right). \quad (A26) $$

The time-evolved state at time-$t$ is

$$ \rho(t, x|x') = \iint G(x, t|y) \rho_{\beta}(y|y') G(y't|x')^* dy dy'. \quad (A27) $$

The energy at time-$t$, using $\lambda(t) = 0$, is

$$ E(t) = \int dx H_{0,x} \rho(t, x|x')|_{x'=x}, \quad (A28) $$

with $H_{0,x} = -\frac{1}{2} \frac{d^2}{dx^2} + \frac{\omega^2 x^2}{2}$. Define

$$ I_1 = \int_0^y \lambda(s) \frac{\sin \omega s}{\sin \omega t} ds \quad (A29) $$

$$ I_2 = \int_0^t \lambda(s) \frac{\sin(\omega(t-s))}{\sin \omega t} ds \quad (A30) $$

$$ A = -iI_1 + iI_2 \left( \frac{\sin \omega t}{\sin i\omega\beta} - \frac{\sin(\omega t + i\omega\beta)}{\sin i\omega\beta} \right) \quad (A31) $$

$$ A' = iI_1 + iI_2 \left( \frac{\sin \omega t}{\sin i\omega\beta} + \frac{\sin(i\omega\beta - \omega t)}{\sin i\omega\beta} \right). \quad (A32) $$
The calculation of the double Gaussian integral in Equation (A27) gives

$$ \rho(t, x|x') = \frac{1}{N(t)} \exp \left( -\frac{\omega}{2}(x^2 + x'^2) \coth \omega\beta + \frac{\omega x x'}{\sinh \omega\beta} + Ax + A'x' \right), \quad (A33) $$

where $N(t)$ is the normalization factor

$$ N(t) = \exp \left( \frac{1}{4} \frac{(A + A')^2}{\omega \coth \omega\beta - \frac{\omega}{\sinh \omega\beta}} \right) \frac{\sqrt{2\pi}}{\sqrt{2(\omega \coth \omega\beta - \frac{\omega}{\sinh \omega\beta})}}. \quad (A34) $$

The action of the Hamiltonian on the propagated state is

$$ H_{0,x} \rho(t, x|x') = \left( \frac{1}{2} \omega \coth \omega\beta - \frac{1}{2} \left( \left( -\omega x \coth \omega\beta + \frac{\omega x'}{\sinh \omega\beta} \right) + A \right)^2 + \frac{\omega^2 x^2}{2} \right) \rho(t, x|x'). \quad (A35) $$

We define the variable $X$ as

$$ X = x - \frac{A + A'}{2 \left( -\omega \coth \omega\beta + \frac{\omega}{\sinh \omega\beta} \right)}. \quad (A36) $$

Then the energy of the propagated state at time $t$ is

$$ E(t) = \frac{\sqrt{2 \left( \omega \coth \omega\beta - \frac{\omega}{\sinh \omega\beta} \right)}}{\sqrt{2\pi}} \int dX \exp \left( - \left( \omega \coth \omega\beta - \frac{\omega}{\sinh \omega\beta} \right) X^2 \right) \times \left[ \frac{1}{2}\omega \coth \omega\beta - \frac{1}{2} \left( X \left( -\omega \coth \omega\beta + \frac{\omega}{\sinh \omega\beta} \right) + \frac{1}{2}(A - A')^2 \right) + \frac{\omega^2}{2} \left( X + \frac{1}{2} \frac{A + A'}{\omega \coth \omega\beta - \frac{\omega}{\sinh \omega\beta}} \right)^2 \right], \quad (A37) $$

and $E(0)$ is the value of $E(t)$ at $t = 0$, so that

$$ E(t) - E(0) = \frac{1}{8} \left( \frac{(A + A')^2}{\left( \coth \omega\beta - \frac{1}{\sinh \omega\beta} \right)^2} - (A - A')^2 \right). \quad (A38) $$

Finally using the values of $A$ and $A'$ in terms of $I_1$ and $I_2$, we obtain

$$ E(t) - E(0) = \frac{1}{2} \left( I_1^2 + I_2^2 + 2I_1 I_2 \cos \omega t \right). \quad (A39) $$

This is independent of $\beta$ and is positive. As a corollary, this result is valid if one propagates any eigenstate of the Hamiltonian $H_0$. One can also derive the classical energy

$$ E(t) - E(0) = \langle H(x(t|x_0, p_0), p(t|x_0, p_0), \lambda = 0) - H(x_0, p_0, \lambda = 0) \rangle_{\rho_{\beta}(\lambda=0)}, \quad (A40) $$

where $\rho_{\beta}(\lambda = 0)$ is the classical thermal state. One uses the equations of motion

$$ x(t|x_0, p_0) = x_0 \cos \omega t + p_0 \frac{\sin \omega t}{\omega} + \omega \int_0^t \lambda(s) \sin \omega(t-s) \, ds \quad (A41) $$

$$ p(t|x_0, p_0) = -\omega x_0 \sin \omega t + p_0 \cos \omega t + \omega^2 \int_0^t \lambda(s) \cos \omega(t-s) \, ds \quad (A42) $$
$$ \rho_{\beta}(\lambda = 0) = \frac{1}{N_{\beta}} \exp \left( -\beta \left( \frac{p^2}{2} + \frac{\omega^2 x^2}{2} \right) \right), \qquad (A43) $$

and then

$$ E(t) - E(0) = \frac{1}{2} \omega^4 \left( \left( \int_0^t \lambda(s) \cos \omega(t-s) \, ds \right)^2 + \left( \int_0^t \lambda(s) \sin \omega(t-s) \, ds \right)^2 \right). \qquad (A44) $$

If $\lambda(0) = 0$ but $\lambda(t) \neq 0$, one gets

$$ E(t) - E(0) = \frac{1}{2} \omega^4 \left( \left( \int_0^t \lambda(s) \cos \omega(t-s) \, ds \right)^2 + \left( \int_0^t \lambda(s) \sin \omega(t-s) \, ds \right)^2 - \frac{2\lambda(t)}{\omega} \int_0^t \lambda(s) \sin \omega(t-s) \, ds \right). \qquad (A45) $$

This can be negative, for example if $\lambda(t) = t$:

$$ E(t) - E(0) = 1 - \frac{t^2 \omega^2}{2} - \cos \omega t < 0. \qquad (A46) $$

## Conflicts of Interest

The authors declare no conflict of interest.

## References and Notes

1. Landau, L.D.; Lifshitz, E.M. *Statistical Physics*; Pergamon Press: Oxford, UK, 1980.

2. Schulman, L.S.; Gaveau, B. Coarse grains: The emergence of space and order. *Found. Phys.* **2001**, *31*, 713–731.

3. Jaynes, E.T. Information theory and statistical mechanics. *Phys. Rev.* **1957**, *106*, 620–630.

4. Jaynes, E.T. Information theory and statistical mechanics II. *Phys. Rev.* **1957**, *108*, 171–190.

5. Gaveau, B.; Schulman, L.S. Master equation based formulation of non-equilibrium statistical mechanics. *J. Math. Phys.* **1996**, *37*, 3897–3932.

6. Gaveau, B.; Schulman, L.S. A general framework for non-equilibrium phenomena: The master equation and its formal consequences. *Phys. Lett. A* **1997**, *229*, 347–353.

7. In [5] and [6] we found the relative entropy to be related to dissipation, but we there dealt with stochastic dynamics and our conclusions were limited to short times.

8. Esposito, M.; Lindenberg, K.; van den Broeck, C. Entropy production as correlation between system and reservoir. *New J. Phys.* **2010**, *12*, 013013.

9. Reeb, D.; Wolf, M.M. (Im-)Proving Landauer’s Principle. **2013**, arXiv:1306.4352v2.

10. Takara, K.; Hasegawa, H.H.; Driebe, D.J. Generalization of the second law for a transition between nonequilibrium states. *Phys. Lett. A* **2010**, *375*, 88–92.

11. Cover, T.M.; Thomas, J.A. *Elements of Information Theory*; Wiley: New York, NY, USA, 1991.

12. Henceforth we suppress the qualifier “inverse” in referring to various $\beta$'s as temperatures.

13. Partovi, M.H. Quantum thermodynamics. *Phys. Lett. A* **1989**, *137*, 440–444.
14. Brillouin, L. *Science and Information Theory*, 2nd ed.; Academic Press: New York, NY, USA, 1962.

15. Landauer, R. Irreversibility and heat generation in the computing process. *IBM J. Res. Dev.* **1961**, *5*, 183–191.

16. Leff, H.S.; Rex, A.F. *Maxwell's Demon: Entropy, Information, Computing*; Princeton University Press: Princeton, NJ, USA, 1990.

17. Schulman, L.S.; Gaveau, B. Ratcheting Up Energy by Means of Measurement. *Phys. Rev. Lett.* **2006**, *97*, 240405.

18. Esposito, M.; van den Broeck, C. Second law and Landauer principle far from equilibrium. *Europhys. Lett.* **2011**, *95*, 40004.

19. Kawai, R.; Parrondo, J.M.R.; Van den Broeck, C. Dissipation: The Phase-Space Perspective. *Phys. Rev. Lett.* **2007**, *98*, 080602.

20. Jarzynski, C. Nonequilibrium Equality for Free Energy Differences. *Phys. Rev. Lett.* **1997**, *78*, 2690–2693.

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).