Quantum Thermodynamics in the Refined Weak Coupling Limit

TL;DR

This paper develops a thermodynamic framework for quantum systems in the refined weak coupling limit, accounting for non-Markovian dynamics and external drivings, and redefines internal energy accordingly.

Contribution

It introduces a novel thermodynamic approach for the refined weak coupling limit, incorporating non-Markovian effects and external drivings in quantum thermodynamics.

Findings

Derived first and second laws in the refined weak coupling regime.

Extended the framework to include external drivings.

Reconsidered the definition of internal energy in this context.

Abstract

We present a thermodynamic framework for the refined weak coupling limit. In this limit the interaction between system and environment is weak, but not negligible. As a result, the system dynamics becomes non-Markovian breaking divisibility conditions. Nevertheless, we propose a derivation of the first and second law just in terms of the reduced system dynamics. To this end, we extend the refined weak coupling limit for allowing slowly-varying external drivings, and reconsider the definition of internal energy due to the non-negligible interaction.

Equations

$$\frac{d{\tilde{\rho }}_{\mathcal{SE}}(t)}{dt}=-\mathrm{i}[{\tilde{V}}_{\mathcal{SE}}(t),{\tilde{\rho }}_{\mathcal{SE}}(t)],$$

$$\frac{d\tilde{\rho}_{\SS}(\tau)}{d\tau}:=\tilde{\mathcal{L}}[\tilde{\rho}_{\SS}(\tau)]=\mathrm{i}[H_{\rm LS},\tilde{\rho}_{\SS}(\tau)]+\sum_{\omega}\sum_{k,l}\gamma_{kl}(\omega)\left[A_{l}(\omega)\tilde{\rho}_{\SS}(\tau)A_{k}^{\dagger}-\tfrac{1}{2}\big{\{}A_{k}^{\dagger}(\omega)A_{l}(\omega),\tilde{\rho}_{\SS}(\tau)\big{\}}\right].$$

$$\underset{\alpha \to 0}{\lim }\underset{0\le {\alpha }^{2}t\le t_0}{\sup }\Vert {\tilde{\Lambda }}_t({\rho }_{\text{SS}})-{\mathrm{e}}^{{\alpha }^2\tilde{\mathcal{L}}t}({\rho }_{\text{SS}})\Vert =0,t_0<\infty ,$$

$${\Lambda }_t\simeq {\mathrm{e}}^{(-\mathrm{i}\mathcal{H}+{\alpha }^2\mathcal{L})t}:={\mathrm{e}}^{{\mathcal{L}}_{\mathrm{D}}t}$$

$${\Lambda }_t\simeq \mathcal{T}{\mathrm{e}}^{{\int }_0^t{\mathcal{L}}_{\mathrm{D}}(s)ds}$$

$${\tilde{\rho }}_{\mathcal{SE}}(t)=\mathcal{T}{\mathrm{e}}^{{\int }_0^t\tilde{\mathcal{V}}(s)ds}[{\rho }_{\mathcal{S}}(0)\otimes {\rho }_{\mathcal{E}}(0)],$$

$${\tilde{\rho }}_{\text{SS}}(t)={\tilde{\Lambda }}_t[{\rho }_{\text{SS}}(0)]={\mathrm{e}}^{\mathcal{Z}(t)}[{\rho }_{\text{SS}}(0)],$$

$$\mathcal{Z}(t)=\sum _{k=1}^{\infty }{\alpha }^k{\mathcal{Z}}_k(t).$$

$${\mathcal{Z}}_2(t)[{\rho }_{\text{SS}}(0)]=-\frac{1}{2}\mathcal{T}{\int }_0^{t}d{t}_1{\int }_0^{t}d{t}_2{\mathrm{Tr}}_{\mathcal{E}}\left[{\tilde{V}}_{\mathcal{SE}}(t_1),\left[{\tilde{V}}_{\mathcal{SE}}(t_2),{\rho }_{\text{SS}}(0)\otimes {\rho }_{\mathcal{E}}(0)\right]\right].$$

$$\displaystyle\mathcal{Z}_{2}(t)[\rho_{\SS}(0)]=-\mathrm{i}[\Lambda(t),\rho_{\SS}(0)]+\mathrm{Tr}_{E}\big{[}\Upsilon(t)\rho_{\SS}(0)\otimes\rho_{\mathscr{E}}(0)\Upsilon(t)-\tfrac{1}{2}\big{\{}\Upsilon^{2}(t),\rho_{\SS}(0)\otimes\rho_{\mathscr{E}}(0)\big{\}}\big{]},$$

$${\mathcal{Z}}_2(t)[{\rho }_{\text{SS}}(0)]=$$

$$\displaystyle+\sum_{\omega,\omega^{\prime}}\sum_{k,l}\Gamma_{kl}(\omega,\omega^{\prime},t)\big{[}A_{l}(\omega^{\prime})\rho_{\SS}(0)A_{k}^{\dagger}(\omega)-\tfrac{1}{2}\{A_{k}^{\dagger}(\omega)A_{l}(\omega^{\prime}),\rho_{\SS}(0)\}\big{]},$$

$${\Gamma }_{kl}(\omega ,{\omega }^{\prime },t):={\int }_0^{t}d{t}_1{\int }_0^{t}d{t}_2{\mathrm{e}}^{\mathrm{i}(\omega t_1-{\omega }^{\prime }{t}_2)}\mathrm{Tr}[{\tilde{B}}_k(t_1-t_2)B_l{\rho }_{\mathcal{E}}(0)],$$

$$\Lambda (t):=\sum _{\omega ,{\omega }^{\prime }}\sum _{k,l}{\Xi }_{kl}(\omega ,{\omega }^{\prime },t)A_k^{†}(\omega )A_l({\omega }^{\prime }),$$

$${\Xi }_{kl}(\omega ,{\omega }^{\prime },t):=\frac{1}{2\mathrm{i}}{\int }_0^{t}d{t}_1{\int }_0^{t}d{t}_2\mathrm{sgn}(t_1-t_2){\mathrm{e}}^{\mathrm{i}(\omega t_1-{\omega }^{\prime }{t}_2)}\mathrm{Tr}[{\tilde{B}}_k(t_1-t_2)B_l{\rho }_{\mathcal{E}}(0)].$$

$${\Lambda }_t={\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}{\mathrm{e}}^{\mathcal{Z}(t)}\simeq {\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}{\mathrm{e}}^{{\alpha }^2{\mathcal{Z}}_2(t)}.$$

$${\Lambda }_t^{\mathrm{R}}:={\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}{\mathrm{e}}^{{\alpha }^2{\mathcal{Z}}_2(t)}\simeq {\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}{\mathrm{e}}^{{\alpha }^2{\tilde{\mathcal{L}}}_{\mathrm{D}}t}={\mathrm{e}}^{{\mathcal{L}}_{\mathrm{D}}t},t\sim {\alpha }^{-2}.$$

$${\mathcal{L}}_{\mathrm{R}}(t)=\left(\frac{d}{dt}{\Lambda }_t^{\mathrm{R}}\right)({\Lambda }_t^{\mathrm{R}}{)}^{-1}=-\mathrm{i}\mathcal{H}+{\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}\left[\frac{d}{dt}{\mathrm{e}}^{{\alpha }^2{\mathcal{Z}}_2(t)}\right]{\mathrm{e}}^{-{\alpha }^2{\mathcal{Z}}_2(t)}{\mathrm{e}}^{\mathrm{i}\mathcal{H}t}=-\mathrm{i}\mathcal{H}+{\mathrm{e}}^{-\mathrm{i}\mathcal{H}t}{\tilde{\mathcal{L}}}_{\mathrm{R}}(t){\mathrm{e}}^{\mathrm{i}\mathcal{H}t},$$

$$\frac{d}{dt}{\mathrm{e}}^{A(t)}={\int }_0^1{\mathrm{e}}^{sA(t)}\left[\frac{d}{dt}A(t)\right]{\mathrm{e}}^{(1-s)A(t)}ds,$$

$${\tilde{\mathcal{L}}}_{\mathrm{R}}(t)=\left[\frac{d}{dt}{\mathrm{e}}^{{\alpha }^2{\mathcal{Z}}_2(t)}\right]{\mathrm{e}}^{-{\alpha }^2{\mathcal{Z}}_2(t)}={\alpha }^2{\int }_0^1{\mathrm{e}}^{s{\alpha }^2{\mathcal{Z}}_2(t)}{\dot{\mathcal{Z}}}_2(t){\mathrm{e}}^{-s{\alpha }^2{\mathcal{Z}}_2(t)}ds.$$

$${\mathrm{e}}^{{\mathcal{L}}_{\mathrm{D}}t}\to \mathcal{T}{\mathrm{e}}^{{\int }_0^t{\mathcal{L}}_{\mathrm{D}}(s)ds}$$

$$\mathcal{T}{\mathrm{e}}^{\overline{{\mathcal{L}}_{\mathrm{D}}(s)t}}=\mathcal{T}{\mathrm{e}}^{{\int }_0^t{\mathcal{L}}_{\mathrm{D}}(s)ds},$$

$${\tilde{\Lambda }}_t^{\mathrm{RA}}:=\mathcal{T}{\mathrm{e}}^{\overline{{\alpha }^2{\mathcal{Z}}_2(t,s)}}=\mathcal{T}{\mathrm{e}}^{\frac{{\alpha }^2}{t}{\int }_0^t{\mathcal{Z}}_2(t,s)ds},$$

$${\Lambda }_t^{\mathrm{RA}}={\mathcal{U}}_t{\tilde{\Lambda }}_t^{\mathrm{RA}},\text{with}{\mathcal{U}}_t:=\mathcal{T}{\mathrm{e}}^{-\mathrm{i}{\int }_0^t\mathcal{H}(s)ds}.$$

$$\mathcal{L}_{\rm D}\big{(}\rho_{\mathscr{S}}^{\beta}\big{)}=0\Rightarrow\mathrm{e}^{\mathcal{L}_{\rm D}t}\big{(}\rho_{\mathscr{S}}^{\beta}\big{)}=\rho_{\mathscr{S}}^{\beta}.$$

$$\lim_{t\to\infty}\mathrm{e}^{\mathcal{L}_{\rm D}t}\big{[}\rho_{\mathscr{S}}(0)\big{]}=\rho_{\mathscr{S}}^{\beta}.$$

$$\mathcal{L}_{\rm D}(t)\Big{[}\rho_{\SS(t)}^{\beta}\Big{]}=0,\quad\text{with}\quad\rho_{\SS(t)}^{\beta}=Z_{\SS(t)}^{-1}\mathrm{e}^{-\beta H_{\mathscr{S}}(t)}.$$

$$\mathcal{T}\mathrm{e}^{\int_{0}^{t}\mathcal{L}_{\rm D}(s)ds}\big{[}\rho_{\mathscr{S}}(0)\big{]}\simeq\rho_{\SS(t)}^{\beta}.$$

$$⟨H(t)⟩=⟨H_{\text{SS}}(t)⟩+⟨H_{\mathcal{E}}⟩+\alpha ⟨V_{SE}⟩=\mathrm{Tr}[{\rho }_{\text{SS}}(t)H_{\text{SS}}(t)]+\mathrm{Tr}[{\rho }_{\mathcal{E}}(t)H_{\mathcal{E}}]+\alpha \mathrm{Tr}[{\rho }_{\mathcal{SE}}(t)V_{\mathcal{SE}}].$$

$$⟨H(t)⟩\simeq \mathrm{Tr}[H_{\mathcal{S}}(t){\rho }_{\mathcal{S}}(t)]+\mathrm{Tr}[H_{\mathcal{E}}{\rho }_{\mathcal{E}}(t)],$$

$$E(t):=\mathrm{Tr}[{\rho }_{\mathcal{S}}(t)H_{\mathcal{S}}(t)].$$

Full text

\pubvolume\issuenum\articlenumber\history\Title

Quantum Thermodynamics in the Refined Weak Coupling Limit \AuthorÁngel Rivas1,2 \AuthorNamesÁngel Rivas

\corresanrivas@ucm.es

1. 1 Introduction

2. 2 Weak Coupling and Refined Weak Coupling Limit

3. 2.1 Weak Coupling Limit

4. 2.2 Refined Weak Coupling Limit

5. 2.3 Refined Weak Coupling Limit under Slowly-Varying Time-Dependent Hamiltonians

6. 3 Standard Thermodynamics of Open Quantum Systems

7. 3.1 The First Law

8. 3.2 The Second Law

9. 3.3 Difficulties beyond the Weak Coupling Limit

10. 4 Thermodynamics in the Refined Weak Coupling Limit

11. 4.1 Time-Independent System Hamiltonian $H_{\SS}$

12. 4.2 Time-Dependent System Hamiltonian $H_{\SS}(t)$

13. 5 Example: Spin-Boson model in the Refined Weak Coupling Limit

14. 6 Conclusions

15. A Internal Energy in the Weak Coupling Limit

16. B Calculation of the Time-Ordered Exponential for the Driven Spin-Boson Model in the Refined Weak Coupling

1 Introduction

The theory of open quantum systems describes the evolution of a quantum system that exchanges energy or information with some environment BrPr02 ; GardinerZoller04 ; Libro . The understanding on whether a thermodynamical framework can be applied to these processes is not only of fundamental, but also of practical interest and constitutes the active field of quantum thermodynamics Qthermo1 ; QthermoKosloff ; Qthermo2 . In analogy with the classical case, one speaks about equilibrium quantum thermodynamics if the open quantum system considered is in equilibrium with its environment, analyzing the change of the system’s thermodynamic properties when there is a change of one equilibrium state into another equilibrium state. On the contrary, nonequilibrium quantum thermodynamics focuses on the change of those properties during the evolution of the open quantum system due to the environmental interaction. Provided that the coupling between system and environment is sufficiently weak, this evolution can be approximated by the celebrated weak coupling limit BrPr02 ; Libro ; Davies , which allows for the formulation of a dynamical equation for the density matrix of the system with the Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) form GKLS1 ; GKLS2 . In absence of external driving, this equation generates a quantum dynamical semigroup AlickiBook with very appealing thermodynamical properties. Namely, the system approaches asymptotically the equilibrium state AlickiBook ; BrPr02 ; Libro , the entropy production is positive SpohnEntropy ; QthermoKosloff , the heat is additive under the presence of several thermal reservoirs in the environment SpohnLebowitz ; QthermoKosloff , etc. Moreover, work can be incorporated into this framework under some conditions, such as slowly-varying or periodically-varying drivings AlickiHeatEngine ; AlickiTutorialPeriodic . However, if the coupling between system and environment is not weak, the evolution of the open system becomes typically non-Markovian, and the formulation of a thermodynamic framework is more complicated Esposito ; Modi ; Tanimura ; Alipour1 ; Alipour2 ; StrasbergPRX ; Pati ; Ghosh ; StrasbergPRE . This is so even for the equilibrium case Gelin ; Seifert ; Hu .

It is the aim of this contribution to introduce a valid thermodynamics framework for the evolution obtained in a refined weak coupling limit AlickiRefined ; SchallerBrandes ; Benatti1 ; Refined ; Libro . This technique still considers the system-environment coupling to be weak, but not negligible at short times. As a consequence, the open system dynamics becomes non-Markovian, and it can actually present strong non-Markovian properties like positive indivisibility (P-indivisibility) or “quasieternal” completely positive indivisibility (CP-indivisibility) Refined . We suggest a reformulation of the internal energy in this regime, such that the entropy production is never negative because the obtained dynamical map is, by construction, completely positive (CP). Moreover, in order to allow for the existence of work, we introduce an extension of the refined weak coupling technique for slowly-varying system Hamiltonians, and establish the first and second law in this setting. As a key difference with other approaches, our thermodynamic relations only involve system observables. Namely, they are completely formulated in terms of the system reduced dynamics. This is something natural in the Markovian case, but very challenging to be satisfied for non-Markovian evolution.

The present paper is structured as follows: in Sec. 2 we briefly explain the usual weak coupling limit and the refined one, indicating their principal differences. Sec. 3 is a succinct, but self-consistent review of the standard thermodynamics of weakly-coupled open quantum systems, that emphasizes the approximations taken, which no longer hold for the non-Markovian case. The thermodynamics formulation for the refined coupling limit is presented in Sec. 4, and constitutes the main result of this work. We complete our study by applying these ideas to the case of a two-level system in contact with a thermal bath in Sec. 5. Finally, some discussion and possible future directions of this approach are outlined in the Conclusions section.

2 Weak Coupling and Refined Weak Coupling Limit

Let $\mathscr{S}$ be an open system in contact with some environment $\mathscr{E}$. The total Hamiltonian for system and environment is given by an expression of the form $H=H_{\mathscr{S}}+H_{\mathscr{E}}+V_{\mathscr{SE}}$, where $H_{\mathscr{S}}$, $H_{\mathscr{E}}$, and $V_{\mathscr{SE}}$ denote the system, environment, and interaction Hamiltonian, respectively. In the interaction picture, the von Neumann equation reads:

[TABLE]

where, for any operator $A$, $\tilde{A}(t):=\exp[\mathrm{i}(H_{\mathscr{S}}+H_{\mathscr{E}})t]A\exp[-\mathrm{i}(H_{\mathscr{S}}+H_{\mathscr{E}})t]$, and we have taken units of $\hbar=1$. Unless otherwise stated, we shall consider units of $\hbar=k_{\rm B}=1$ in the following. In addition, we shall take the initial time $t=0$ to be the time where system and environment start interacting, $\rho_{\mathscr{SE}}(0)=\rho_{\mathscr{S}}(0)\otimes\rho_{\mathscr{E}}(0)$, so that the system dynamics can be described in terms of a dynamical map $\Lambda_{t}$, i.e., a linear, trace preserving, and CP map.

2.1 Weak Coupling Limit

In order to obtain $\Lambda_{t}$ from (1), some assumptions are usually required, the weak coupling case being the most common one. In the standard approach to this situation, the strength of the interaction is gauged by some coupling constant $\alpha$, so that $V_{\mathscr{SE}}$ is substituted by $\alpha V_{\mathscr{SE}}$ in (1), and the limit $\alpha\to 0$ is taken on a rescaled time $\tau\to\alpha^{2}t$. This leads to a GKLS equation Davies ; AlickiBook ; Libro :

[TABLE]

Here, $\tilde{\mathcal{L}}$ is the Liouville operator (in the interaction picture), and $A_{k}(\omega)$ are eigenoperators of the system Hamiltonian, $[H_{\mathscr{S}},A_{k}(\omega)]=-\mathrm{i}\omega A_{k}(\omega)$, such that the interaction can always be written in the form of $V_{\mathscr{SE}}=\sum_{k}\sum_{\omega}A_{k}(\omega)\otimes B_{k}$ with $B_{k}$ Hermitian bath operators (see, e.g., Libro ). Moreover, $H_{\rm LS}$ is a Hamiltonian Lamb-shift-type correction $[H_{\mathscr{S}},H_{\rm LS}]=0$, and $\gamma_{kl}(\omega)$ are the elements of a positive-semidefinite matrix (see the details in BrPr02 ; AlickiBook ; Libro ). In addition, $\rho_{\mathscr{E}}(0)$ is assumed to be some Gaussian state of an environment with infinitely many degrees of freedom and $\mathrm{Tr}_{\mathscr{E}}[V_{\mathscr{SE}}\rho_{\mathscr{E}}(0)]=0$. This last condition can always be achieved by a proper redefinition of system and interaction Hamiltonians Libro . Under some relatively mild conditions on the environmental correlation functions, Davies rigorously shows Davies :

[TABLE]

for any initial system state $\rho_{\SS}$. Here, $\tilde{\Lambda}_{t}$ denotes the exact dynamical map (in the interaction picture). Recently, precise bounds to this convergence have been formulated Merkli . The Davies’ result allows for the approximation

[TABLE]

for $t\lesssim\alpha^{-2}$. Here, $\mathcal{L}_{\rm D}$ is the so-called Davies generator, $\mathcal{H}:=[H_{\mathscr{S}},\cdot]$, working again in the Schrödinger picture, and the relation $[\mathcal{H},\mathcal{L}]=0$ is fulfilled. This is the celebrated weak coupling limit.

Furthermore, this technique can be extended for time-dependent Hamiltonians $H_{\SS}(t)$. If $\tau_{H}$ is the typical time for the variation of $H_{\SS}(t)$, in the simplest case, this change is considered to be slow in comparison to the evolution of the open system due to the coupling to the environment, i.e., $\tau_{H}\gtrsim\alpha^{-2}$ DaviesSpohn . Then, we can “adiabatically” deform the generator $\mathcal{L}_{\rm D}\to\mathcal{L}_{\rm D}(t)$, so that $\mathcal{L}_{\rm D}(t)$ is the Davies generator calculated instantaneously for Hamiltonian $H_{\SS}(t)$, and the approximation

[TABLE]

is satisfied, where $\mathcal{T}$ is the time-ordering operator. We should note that the term “adiabatic” is used throughout the text with its standard meaning in quantum mechanics referring to slow transformations. We do not intend to denote lack of heat exchange as in classical thermodynamics.

2.2 Refined Weak Coupling Limit

Despite enjoying very remarkable properties, the Davies’ weak coupling technique has a problem for very short times $t\ll\alpha^{-2}$, as in this case, $\tilde{\Lambda}_{t}\simeq\exp(\alpha^{2}\tilde{\mathcal{L}}t)\simeq\mathds{1}$, and no open system dynamics can be resolved. In order to see this short time scale, a refined weak coupling must be formulated AlickiRefined ; SchallerBrandes ; Benatti1 ; Refined ; Libro . Nonetheless, one has to be careful to prevent violations of complete positivity noCP1 ; noCP2 ; noCP3 . To this end, we can solve (1) as a time-ordered exponential:

[TABLE]

with $\tilde{\mathcal{V}}(t):=-\mathrm{i}[\tilde{V}(t),\cdot]$, and look for a solution of the reduced dynamics in the form of

[TABLE]

where the exponent $\mathcal{Z}(t)$ is formally defined as the logarithm of $\tilde{\Lambda}_{t}$. Since $\tilde{\Lambda}_{t}$ can be expressed as a power expansion in $\alpha$, we can assume a similar expansion for the exponent:

[TABLE]

The terms $\mathcal{Z}_{k}(t)$ can be computed expanding the exponential and by comparison with the terms of the same exponent of $\alpha$ in the expansion of $\eqref{Tord}$ after taking the partial trace on the environment AlickiRefined ; SchallerBrandes ; Benatti1 ; Refined . The first order vanishes because $\mathrm{Tr}_{\mathscr{E}}[V_{\mathscr{SE}}\rho_{\mathscr{E}}(0)]=0$, so that the first nontrivial contribution is

[TABLE]

Applying $\mathcal{T}$ under the integral signs and reordering terms, we obtain

[TABLE]

with Hermitian operators $\Lambda(t):=\tfrac{1}{2\mathrm{i}}\int_{0}^{t}dt_{1}\int_{0}^{t}dt_{2}\mathrm{sgn}(t_{1}-t_{2})\mathrm{Tr}_{E}[\tilde{V}_{\mathscr{SE}}(t_{1})\tilde{V}_{\mathscr{SE}}(t_{2})\rho_{\mathscr{E}}(0)]$ and $\Upsilon(t):=\int_{0}^{t}\tilde{V}_{\mathscr{SE}}(s)ds$. In terms of eigenoperators of $H_{\SS}$, this is written as AlickiRefined ; Refined

[TABLE]

where

[TABLE]

and

[TABLE]

with

[TABLE]

From (10), we infer that $\mathcal{Z}_{2}(t)$ has the GKLS form GKLS1 ; GKLS2 , so it turns out that the coefficients $\Gamma_{kl}(\omega,\omega^{\prime},t)$ form a positive-semidefinite matrix. This property seems not ensured beyond the second order Schaller2 .

Thus, the refined weak coupling limit is performed by neglecting higher orders on $\alpha$ in the exponent, $\mathcal{Z}(t)\simeq\mathcal{Z}_{2}(t)$:

[TABLE]

By construction, this is a completely positive map [$\mathcal{Z}_{2}(t)$ has the GKLS form for all $t$] which approaches the exact dynamics in the short time limit $t\ll\alpha^{-2}$ (note that, in general, $[\mathcal{Z}_{2}(t),\mathcal{H}]\neq 0$). In the large time scale, $t\sim\alpha^{-2}$, it can be proven AlickiRefined ; SchallerBrandes ; Refined that $\alpha^{2}\mathcal{Z}_{2}(t)\simeq\tilde{\mathcal{L}}_{\rm D}t$, and the refined weak coupling dynamics $\Lambda_{t}^{\rm R}$ approaches the weak coupling limit Davies’ semigroup,

[TABLE]

Hence, thanks to the Davies’ result (3), the dynamical map $\Lambda_{t}^{\rm R}$ provides a nontrivial and consistent approximation for the weak coupling dynamics for the same time scale as the Davies’ semigroup, $0\leq t\lesssim\alpha^{-2}$. However, $\Lambda_{t}^{\rm R}$ resolves the exact dynamics in the small time scale. One may notice that under the assumption that the width of the reservoir correlation functions is negligible in comparison with the time scale of the open system, $\exp(\mathcal{L}_{\rm D}t)$ and $\Lambda_{t}^{\rm R}$ coincide as the integrals in Eqs. (12) and (14) can be effectively extended to infinity for all times.

We can calculate the generator $\mathcal{L}_{\rm R}(t)$ of $\Lambda_{t}^{\rm R}$, i.e., $\frac{d}{dt}\Lambda_{t}^{\rm R}=\mathcal{L}_{\rm R}(t)\Lambda_{t}^{\rm R}$, by applying inverses and derivatives,

[TABLE]

where $\tilde{\mathcal{L}}_{\rm R}(t)$ is the generator in the interaction picture that can be written in terms of the time derivative of $\mathcal{Z}_{2}(t)$ with the help of the Snider-Wilcox formula SWformula ,

[TABLE]

so that

[TABLE]

Here, we have used the “overdot” notation for the time derivative. Since $\alpha^{2}\mathcal{Z}_{2}(t)\simeq\tilde{\mathcal{L}}_{\rm D}t$, from Eq. (17), it is straightforward to check that for large $t$, $\mathcal{L}_{\rm R}(t)\simeq\mathcal{L}_{\rm D}$.

Although we can think of the refined weak coupling as a small correction to the Markovian weak coupling evolution for short times, in general, the dynamical map $\Lambda_{t}^{\rm R}$ turns out to be highly non-Markovian noMarkov1 ; noMarkov2 ; noMarkov3 , breaking conditions such as P-divisibility Refined .

2.3 Refined Weak Coupling Limit under Slowly-Varying Time-Dependent Hamiltonians

As far as we know, the refined weak coupling method for a time-dependent system Hamiltonian $H_{\SS}(t)$ has not been properly studied in the literature. A priori, the simplest situation should be again the case of slowly-varying $H_{\SS}(t)$, i.e., $\tau_{H}\gtrsim\alpha^{-2}$, but the problem becomes now more involved than in the standard weak coupling. We would expect an “adiabatic deformation” of the refined weak coupling generator such that it approaches the aforementioned adiabatically-deformed Davies generator $\mathcal{L}_{\rm D}(t)$ for large times. However, it is not straightforward to obtain this result keeping CP.

In order to do so, notice that, on the one hand, if $\mathcal{Z}_{2}(t,s)$ is the refined weak coupling exponent calculated for the Hamiltonian $H_{\SS}(s)$, for $s$ fixed, we expect $\alpha^{2}\mathcal{Z}_{2}(t,s)\simeq\tilde{\mathcal{L}}_{\rm D}(s)t$ for large $t$ (any $s$). On the other hand, we may observe that, in the weak coupling limit, the substitution

[TABLE]

when passing from $H_{\SS}$ to $H_{\SS}(t)$ in the adiabatic approximation can be seen as a kind of “averaging” approximation. Namely, we construct $\mathcal{L}_{\rm D}(s)$ for $H_{\SS}(s)$, for $s$ fixed, and approximate the evolution as

[TABLE]

where the temporal average of some function of $s$, $f(s)$, has been introduced $\overline{f(s)}:=\frac{1}{t}\int_{0}^{t}f(s)ds$. Thus, the evolution from [math] to $t$ is approximately given, modulo time-ordering, by the exponential of the “time-averaged” value of the exponent $\mathcal{L}_{\rm D}(s)t$ from [math] to $t$.

Similarly, in the refined weak coupling limit, we can perform an adiabatic approximation for slowly-varying $H_{\SS}(t)$ by

[TABLE]

which in the Schrödinger picture yields

[TABLE]

As required, $\Lambda_{t}^{\rm RA}$ approaches (5) for large $t$ and reduces to $\Lambda_{t}^{\rm R}$ for time-independent system Hamiltonians. Furthermore, $\Lambda_{t}^{\rm RA}$ is, by construction, a valid CP map because $\mathcal{Z}_{2}(t,s)$ has the GKLS form for any $t$ and $s$ (and $t>0$).

3 Standard Thermodynamics of Open Quantum Systems

We shall assume from now on the environment is a thermal reservoir or bath. Namely, $\mathscr{E}$ is a system with infinitely many degrees of freedoms, which remains initially in the canonical Gibbs state $\rho_{\mathscr{E}}=\rho_{\mathscr{E}}^{\beta}=Z_{\mathscr{E}}^{-1}\exp(-\beta H_{\mathscr{E}})$. If the system Hamiltonian $H_{\SS}$ is time independent, the system Gibbs state $\rho_{\mathscr{S}}^{\beta}=Z_{\SS}^{-1}\exp(-\beta H_{\mathscr{S}})$ is a steady state of the Davies semigroup AlickiBook ; BrPr02 ; Libro ; Davies ,

[TABLE]

In fact, provided that system-reservoir coupling $V_{\mathscr{SE}}$ only allows $\rho_{\mathscr{S}}^{\beta}$ as a steady state, the evolution of any system eventually approaches this unique steady state AlickiBook ; Libro ; BrPr02 ; Spohn2 ; SpohnLebowitz ; Wolf ,

[TABLE]

In the case of a slowly-varying time-dependent system Hamiltonian $H_{\SS}(t)$, the adiabatically-deformed Davies $\mathcal{L}_{\rm D}(t)$ generator fulfills

[TABLE]

Provided that $\rho_{\SS(t)}^{\beta}$ is the only state satisfying this, and $H_{\SS}(t)$ remains almost constant during the relaxation time –very slow variation limit–, for large enough $t$, any initial state approaches $\rho_{\SS(t)}^{\beta}$,

[TABLE]

According to the first law of thermodynamics, the change of the internal energy of $\mathscr{S}$ can be divided into work and heat. The former is usually identified with a controllable and measurable change, for instance by means of variations in the parameters of $H_{\mathscr{S}}$ (which must be then time-dependent), whereas the latter is a change generally out of our capability to control and observe, and it is often identified with reservoir properties. However, the specific definition of work and heat is in general difficult because the identification of internal energy for a given system-reservoir Hamiltonian $H(t)=H_{\SS}(t)+H_{\mathscr{E}}+\alpha V_{\mathscr{SE}}$ is also difficult. It is natural to include $H_{\SS}(t)$ inside the system internal energy, but it is not clear which part of the interaction term $\alpha V_{SE}$ should be considered as “internal” to the system. In general, we have

[TABLE]

Nevertheless, in the weak coupling limit $\alpha^{2}\to 0$,

[TABLE]

hence the internal energy can be defined via

[TABLE]

Actually, under the condition $\mathrm{Tr}_{\mathscr{E}}(V_{\mathscr{SE}}\rho_{\mathscr{E}}^{\beta})=0$, it can be proven that the neglected term in (29) is order $\alpha^{2}$, instead of order $\alpha$ (see Appendix A).

3.1 The First Law

Taking the time derivative in the internal energy (30),

[TABLE]

and heat and work become defined via its time derivative:

[TABLE]

so that

[TABLE]

with the heat flow $\dot{Q}(t)$ into the system and applied power $\dot{W}(t)$ at time $t$ positive if they increase the system energy. Consequently, the integrated form of the first law reads

[TABLE]

with $\Delta E(t):=E(t)-E(0)$, and

[TABLE]

3.2 The Second Law

In order to derive the second law, we define the thermodynamic entropy of the system by

[TABLE]

where $S(\rho)=-\mathrm{Tr}(\rho\log\rho)$ is the von Neumann entropy, so that

[TABLE]

in units of $k_{B}=1$. The quantum relative entropy between two states $\rho_{1}$ and $\rho_{2}$ is defined by

[TABLE]

It can be proven Monotonic that $S(\rho_{1}\|\rho_{2})$ is monotonic under any CP and trace preserving map $\Phi$,

[TABLE]

Recently, the proof has been extended to any positive and trace preserving map $\Phi$ PMonotonic . Suppose $\mathcal{L}(t)$ is any time-dependent GKLS generator with a steady state $\rho_{\rm ss}(t)$, $\mathcal{L}(t)[\rho_{\rm ss}(t)]=0$. Considering the CP map $\exp[\mathcal{L}(t)r]$ with $r\geq 0$, we have

[TABLE]

which implies

[TABLE]

and, particularly for $r=0$, the inequality SpohnEntropy

[TABLE]

Here, for the time derivative of the thermodynamic entropy, one has the general result Wilde :

[TABLE]

The bound (44) is sometimes often called Spohn’s inequality QthermoKosloff ; Qthermo2 ; AlickiTutorialPeriodic ; StrasbergPRX . In the (adiabatic) weak coupling limit (5), Eq. (45) becomes

[TABLE]

Since the canonical Gibbs state $\rho_{\SS(t)}^{\beta}$ is a steady state of $\mathcal{L}_{\rm D}(t)$, Eq. (26), the inequality (44) combined with (46) and (32) leads to

[TABLE]

This is the differential form of the second law of nonequilibrium thermodynamics in the weak coupling limit. It implies the integrated form:

[TABLE]

with $\Delta S(t)=S(t)-S(0)$. Namely, the entropy production [here also its rate (47)] due to the interaction with the reservoir cannot be negative for any final time $t$.

3.3 Difficulties beyond the Weak Coupling Limit

When the interaction term $\alpha V_{\mathscr{SE}}$ cannot be neglected, we cannot expect inequality (47) to be satisfied. In general, the open system evolution becomes non-Markovian noMarkov1 ; noMarkov2 ; noMarkov3 , and the generator of such a dynamics, $\mathcal{L}(t)$, provided it is well defined, does not have the GKLS form. Spohn’s inequality (44) cannot be applied.

The problems go a step back because, as aforementioned, it is not clear which part of the total energy $\langle H(t)\rangle=\langle H_{\SS}(t)\rangle+\langle H_{\mathscr{E}}(t)\rangle+\alpha\langle V_{\mathscr{SE}}(t)\rangle$ must be considered as internal energy of the system. A possible “extreme-type” of splitting is $\dot{E}(t)=\mathrm{Tr}[(H_{\SS}+\alpha V_{\mathscr{SE}})\dot{\rho}_{\mathscr{SE}}(t)]$, and so $\dot{Q}(t)=\mathrm{Tr}[H_{\mathscr{E}}\dot{\rho}_{\mathscr{SE}}(t)]$; then, it can be proven that the inequality (48) is always satisfied StrasbergPRX ; Esposito ; Tanimura . However, it is satisfied for a reservoir of any size, and this might generate some criticism reading its strict equivalence to the phenomenological second law of nonequilibrium thermodynamics.

Another strategy for non-negligible $\alpha$ in the case of time-independent $H_{\SS}$ is based on the fact that

[TABLE]

is the global equilibrium state of system and reservoir. Therefore, the reduced state

[TABLE]

is a fixed point of the system dynamics. Formally, one can write this reduced state as a Gibbs state

[TABLE]

a Hamiltonian of “mean” force Gelin ; Seifert ; Hu . This suggests a possible choice of internal energy as

[TABLE]

which approaches the weak coupling internal energy as $\alpha\to 0$, Eq. (30). Then, in absence of work, $\tilde{Q}=\mathrm{Tr}[H_{\SS}^{*}\rho_{\SS}(t)]$. However, this choice does not allow for the use of the monotonicity of the relative entropy (41) to obtain the second law. If the system is initially in the state $\rho_{\SS}^{\rm ss}$, it remains invariant provided that the initial system-reservoir state is in the total canonical ensemble (49). Since this is not a product state, the reduced dynamics would not be given by a general (CP) dynamical map.

A slightly different choice, originally due to Seifert Seifert , defines internal energy at equilibrium by using a classical thermodynamic relation:

[TABLE]

which suggests the out of equilibrium definition StrasbergPRE :

[TABLE]

with $H_{\SS}^{*}(t)$ defined as $H_{\SS}^{*}$ in Eq. (51) for the canonical Gibbs state with a time-dependent $H_{\SS}(t)$. In such a case, the classical definition of free energy $F:=-\beta^{-1}\log Z_{\SS}$ jointly with the relation $F=E-TS$, motivates the following redefinition of free energy and thermodynamic entropy:

[TABLE]

so that $F^{*}=E^{*}-TS^{*}$. It is possible to obtain an equation with the form of (48) for the dynamics with initial system-reservoir states to be either the total canonical ensemble (49) or belonging to a class of zero discord states StrasbergPRE .

Other approaches beyond the weak coupling framework that are based on e.g., reservoir full counting statistics or coordinate mappings can be found in Chaps. 11 and 15, and 22-25 of Qthermo2 and the references therein, respectively.

4 Thermodynamics in the Refined Weak Coupling Limit

The analysis of thermodynamic properties in the refined weak coupling limit is challenging because of several reasons. Since the $\Lambda_{t}^{\rm R}$ is in general non-Markovian, the breaking of P-divisibility precludes the use of the Spohn’s inequality (44) to formulate the second law Pati . For non-Markovian dynamics, where there is a back-and-forth of information between system and environment, it is reasonable to consider that global system-reservoir information measures must be included in the second law Esposito ; Alipour2 ; StrasbergPRX . However, we will take a different approach here.

Moreover, in the refined weak coupling limit, we cannot neglect the interaction term $\alpha V_{\mathscr{SE}}$ in the total energy (28) for short times. This can be seen with the following argument. If the initial system-reservoir state is taken to be

[TABLE]

under the hypothesis that the $\alpha V_{\mathscr{SE}}$ is negligible, we have for the total Gibbs state (49)

[TABLE]

Since (49) is a stationary state of the total dynamics, Eq. (58) would imply $\rho^{\beta}_{\SS}$ is a steady state of the reduced system dynamics in the refined weak coupling limit. However, $\rho^{\beta}_{\SS}$ is not a steady state in the refined weak coupling for finite $t$,

[TABLE]

but

[TABLE]

Hence, the identifications (30) and (32) for internal energy and heat in absence of work, respectively, cannot be assumed to be true in the refined weak coupling limit for finite times.

4.1 Time-Independent System Hamiltonian $H_{\SS}$

Let us consider first the case of a time-independent system Hamiltonian $H_{\SS}$, so that no work is applied to/performed by the system. We could adopt the definition $E^{*}(t)$ (54) for internal energy. This has the drawback that, since (30) does not coincide with (53) unless the weak interaction term $\alpha V_{\mathscr{SE}}$ is neglected, the equilibrium internal energy in the refined weak coupling would not coincide with the weak coupling one (30) for large times, despite the fulfillment of (60). In addition, the refined weak coupling limit assumes an initial system-reservoir product state, which implies $\langle H(0)\rangle=\langle H_{\SS}(0)\rangle+\langle H_{\mathscr{E}}(0)\rangle$, as assumed $\mathrm{Tr}_{\mathscr{E}}[V_{\mathscr{SE}}\rho_{\mathscr{E}}(0)]=0$. Then, the internal energy at $t=0$ must be $\langle H_{\SS}(0)\rangle$, and this is not obtained with the choice (53). Therefore, by construction, the internal energy in the refined weak coupling regime should be defined such that it fits the weak coupling internal energy $E(t)=\mathrm{Tr}[\rho_{\SS}(t)H_{\SS}(t)]$ at $t=0$ and for asymptotic times, but differs from $E(t)$ at finite times.

Such a definition is nevertheless possible by following a similar argument to the Hamiltonian of a mean force (51). That is, we can write the time-evolution of the system Gibbs state $\rho_{\SS}^{\beta}=Z_{\SS}^{-1}\exp(-\beta H_{\SS})$ in a canonical form:

[TABLE]

with $Z_{\SS}^{\rm R}(t)=Z_{\SS}$ ($\Lambda_{t}^{\rm R}$ is trace preserving) and

[TABLE]

Then, we introduce the refined weak coupling internal energy by analogy to (30) as

[TABLE]

This definition differs from the standard weak coupling energy $E(t)$ (30) for finite times, but it presents the required limiting behavior:

$E^{\rm R}(0)=E(0)$. Thus, the deviation of $E^{\rm R}(t)$ from $E(t)$ at finite times is unambiguously caused by the interaction term $V_{\mathscr{SE}}$ in the Hamiltonian. 2. 2.

$E^{\rm R}(t)$ approaches $E(t)$ for $t$ large. In such a case $\Lambda_{t}^{\rm R}$ approaches the Davies’ semigroup (16) and then, $H_{\SS}^{\rm R}(t)$ approaches $H_{\SS}$.

Given that no work is considered, the first law reads

[TABLE]

for the refined weak coupling heat, or alternatively

[TABLE]

with

[TABLE]

the refined weak coupling heat flow.

The second law can be derived in the integrated form. Since $\Lambda_{t}^{\rm R}$ is trace preserving and CP for any $t$, the monotonicity of the relative entropy (41) gives

[TABLE]

which can be straightforwardly recast in the form

[TABLE]

The equality is reached for $\rho_{\SS}(0)=\rho^{\beta}_{\SS}$. This is a property shared with the standard weak coupling limit, but, as a key difference, $\rho_{\SS}^{\beta}$ is not invariant for the refined weak coupling [see Eqs. (59) and (60)].

Therefore, the present choice of internal energy is thermodynamically consistent for a system coupled with a thermal reservoir in the refined weak coupling limit (in absence of work). Note that, as aforementioned, since $\Lambda^{\rm R}_{t}$ is in general not P-divisible, a differential form for the second law as in (47) cannot be expected to be true in the refined case for finite times.

4.2 Time-Dependent System Hamiltonian $H_{\SS}(t)$

The situation becomes considerably more intricate if we allow for a time-dependent system Hamiltonian $H_{\SS}(t)$. This is so even for slowly-varying $H_{\SS}(t)$ such that the dynamics is well approximated by the adiabatically-deformed refined weak coupling method, $\Lambda_{t}^{\rm RA}$, Eqs. (22) and (23).

In analogy to (30) and (63), we can define internal energy by an equation of the form

[TABLE]

with some appropriate choice of Hermitian operator $H_{\SS}^{\rm RA}(t)$ satisfying $H_{\SS}^{\rm RA}(0)=H_{\SS}(0)$, and $H_{\SS}^{\rm RA}(t)=H_{\SS}^{\rm R}(t)$ in (62) for a time-independent $H_{\SS}$. However, the straightforward generalization of (62) with the changes $\Lambda_{t}^{\rm R}\to\Lambda_{t}^{\rm RA}$ and $H_{\SS}\to H_{\SS}(t)$ is problematic. We should notice that the refined heat flow in (66) has a term with the time derivative $\dot{H}_{\SS}^{\rm R}(t)$, which is not zero because the dynamics is not given by a semigroup in the refined weak coupling limit. Hence, the choice of $H_{\SS}^{\rm RA}(t)$ is a subtle point, because the term $\mathrm{Tr}[\rho_{\SS}(t)\dot{H}_{\SS}^{\rm RA}(t)]$ in the time derivative of the internal energy $E^{\rm RA}(t)$ cannot be unambiguously identified with work.

In order to find a suitable definition work, we may argue that the composition of the system and reservoir forms a closed system, so any energy change in this global system can be unambiguously identified with work. Thus, for $H(t)=H_{\SS}(t)+H_{\mathscr{E}}+V_{\mathscr{SE}}$, using the von Neumann equation for $\dot{\rho}_{\mathscr{SE}}(t)$, we obtain

[TABLE]

for the power applied to/performed by the global system. However, since it acts only on $\SS$ via the variation of $H_{\SS}(t)$, we adopt it as the definition of power for $\SS$ Esposito ; Modi ; Tanimura ; StrasbergPRX ; StrasbergPRE also in the refined weak coupling limit:

[TABLE]

Having in mind this definition of work, a convenient definition of $H_{\SS}^{\rm RA}(t)$ turns out to be:

[TABLE]

where $\Lambda^{\star}$ denotes the Heisenberg adjoint of $\Lambda$, $\mathrm{Tr}[\Lambda(A)B]=\mathrm{Tr}[A\Lambda^{\star}(B)]$. This is a Hermitian operator which satisfies the requirements $H_{\SS}^{\rm RA}(0)=H_{\SS}(0)$, and for time-independent $H_{\SS}$, $H_{\SS}^{\rm RA}(t)=H_{\SS}^{\rm R}(t)$, as we desired.

Thus, we state the first law in the form

[TABLE]

with refined weak coupling heat flow given by

[TABLE]

This way, for time-independent $H_{\SS}$, $\dot{Q}^{\rm RA}(t)=\dot{Q}^{\rm R}(t)$ in (66), and the first law (73) reduces to the case of the previous section (64). For time-dependent $H_{\SS}(t)$, the last term in (74) aims to subtract the direct dependence of the first two terms with the time derivative $\dot{H}_{\SS}(t)$ which we have considered to define work, Eq. (71). Hence, the integrated heat is

[TABLE]

In order to derive the second law, we define the auxiliary object

[TABLE]

which satisfies $\Omega(t,t)=H_{\SS}^{\rm RA}(t)$. A straightforward computation gives the heat written in terms of $\Omega(t,r)$ as

[TABLE]

For the state

[TABLE]

the monotonicity of the relative entropy (41) gives

[TABLE]

which can be recast in the form

[TABLE]

Since this is fulfilled for all $r$, and particularly for $r=t$, according to (77) we obtain the second law:

[TABLE]

This completes the thermodynamic formulation of the refined weak coupling limit.

On the other hand, it is worth noticing that the refined internal energy $\dot{E}^{\rm RA}(t)$ approaches $\dot{E}(t)$ at large times for time-dependent $H_{\SS}(t)$ in the very slow variation limit where Eq. (27) holds true. However, in this limit, the work is actually approaching zero in the time scale $t\lesssim\alpha^{-2}$, where the refined and the standard weak coupling limit are good descriptions of the exact dynamics. This manifests that the asymptotic thermodynamic behavior of the refined weak coupling limit for slowly-varying $H_{\SS}(t)$ is a nontrivial point.

5 Example: Spin-Boson model in the Refined Weak Coupling Limit

As an example of the previous results, we can consider the (transverse) spin-boson problem studied in Refined in the refined weak coupling limit. Consider first the case of a time-independent system Hamiltonian,

[TABLE]

with $H_{\SS}=\frac{\omega_{0}}{2}\sigma_{z}$, $H_{\mathscr{E}}=\sum_{k}\omega_{k}a_{k}^{\dagger}a_{k}$, and $V_{\mathscr{SE}}=\sum_{k}\sigma_{x}g_{k}(a_{k}+a_{k}^{\dagger})$. The exponent $\mathcal{Z}_{2}(t)$ is obtained in the form Refined :

[TABLE]

Here, $\sigma_{\pm}=(\sigma_{x}\pm i\sigma_{y})/2$ are the lowering and raising Pauli matrices and the coefficients are given by

[TABLE]

where $J(\omega)$ is the spectral density of the bath, $\bar{n}_{\beta}(\omega)=[\exp(\beta\omega)-1]^{-1}$ is the mean number of bosons in the bath with frequency $\omega$, and $\mathrm{sinc}(\omega):=\tfrac{\sin\omega}{\omega}$.

We can then compute the refined weak coupling internal energy. In this case, the map $\Lambda^{\rm R}_{t}=\mathrm{e}^{-\mathrm{i}\mathcal{H}_{z}t}\mathrm{e}^{\mathcal{Z}_{2}(t)}$ is an incoherent map in the $\sigma_{z}$ eigenbasis ($\mathcal{H}_{z}:=\tfrac{\omega_{0}}{2}[\sigma_{z},\cdot]$), so that the operator $H_{\SS}^{\rm R}$ in (62) remains diagonal in the same basis as $H_{\SS}$:

[TABLE]

A straightforward computation yields

[TABLE]

[TABLE]

The non-Markovianity of this model has been studied in Refined , showing that it breaks P-divisibility and it is “quasieternal” CP-indivisible.

We have examined the entropy production (68) and the refined internal energy (63) for three different initial states, the excited $|e\rangle$ and ground $|g\rangle$ states, and the superposition state $|+\rangle_{y}$, $\sigma_{y}|+\rangle_{y}=|+\rangle_{y}$. The results are shown in Fig. 1 for three different temperatures, jointly with the results obtained for the Davies quantum dynamical semigroup in the standard weak coupling limit. We find the expected convergence towards the internal energy and the entropy production of the weak coupling. Despite the model is not P-divisible, we do not see oscillations in the entropy production. This is probably caused by the diagonal form of the operator $H_{\SS}^{\rm R}(t)$. These results seem to indicate that the larger difference between the standard weak coupling limit and the refined one arises at low bath temperatures. This is in agreement of what was obtained in Refined , and somehow expected because the width of the reservoir correlation functions increases when $T$ decreases MME .

As a further example including a system time-dependent Hamiltonian, we can consider the same model with a diagonal driving given by the frequency modulation:

[TABLE]

For $\nu$ small enough, the dynamics can be approximated by the adiabatically-deformed refined weak coupling limit $\Lambda^{\rm RA}_{t}$, Eqs. (22) and (23). Thus, the structure of the dynamics is the same as before with the generator $\mathcal{Z}_{2}(t,s)$ as in Eqs. (83)–(87), but $\omega_{0}(s)$ in the place of the formerly constant $\omega_{0}$.

Despite involving an intricate time-ordering operation, if the system is initially in a diagonal state in the $\sigma_{z}$-eigenbasis, e.g., the ground state $|g\rangle\langle g|$, the calculation of $\Lambda^{\rm RA}_{t}(|g\rangle\langle g|)$ can be simplified by a series of algebraic manipulations explained in detail in Appendix B. The same methods can be applied in the computation of the operator $H^{\rm RA}_{\SS}(t)$, Eq. (72), which also remains diagonal for all $t$. This was done analytically up to the numerical computation of two (standard) integrals. The results are depicted in Fig. 2 for several temperatures and modulation frequencies.

Since $\Lambda^{\rm RA}_{t}$ approaches (5) for large times, convergence to the adiabatically-deformed weak coupling limit was found for the power (71). For the internal energy, there is a $\nu$-decreasing small mismatch between the internal energy given by both techniques even at large time. As aforementioned, the full equivalence is found for very small $\nu$. The entropy production is positive and again its rate of change is also positive. This was something expected after the results for the time-independent case and the fact that the modulation keeps $H_{\SS}(t)$ diagonal in the same basis. Finally as $\nu$ approaches zero we obtain consistency with the internal energy for the refined weak coupling with a time-independent system Hamiltonian.

6 Conclusions

In this contribution, we have formulated a thermodynamic framework for the refined weak coupling limit, which represents a non-Markovian approach to the evolution of open quantum systems. We have included the case of time-independent system Hamiltonians, as well as slowly-varying time-dependent ones. To this end, we have extended the refined weak coupling techniques to deal with slowly-varying time-dependent Hamiltonians and redefine the internal energy because of the non-negligible effect of the system-environment interaction. In the appropriate limiting situation, this refined internal energy approaches the usual internal energy in the weak coupling. Finally, we have illustrated the results by analyzing the refined internal energy and the entropy production of a two-level system in contact with a thermal bath. In this example, both operators $H^{\rm R}_{\SS}(t)$ and $H^{\rm RA}_{\SS}(t)$ are diagonal in the same basis as $H_{\SS}$, which is probably a reason way we see not only positive entropy production but also a positive entropy production rate for the initial states analyzed.

As a future work, several extensions of this approach may be considered. For instance: the cases of periodic time-dependent $H_{\SS}(t)$ under the Floquet formalism AlickiTutorialPeriodic , and systems coupled with several heat baths. Furthermore, since some of the arguments employed in Sec. 4 only rely on the fact that $\Lambda_{t}^{\rm R}$ and $\Lambda_{t}^{\rm RA}$ are CP dynamical maps, one may wonder if the same proposal can be used to formulate a thermodynamic framework for any CP dynamical map. This seems an interesting possibility to be analyzed in the near future.

Acknowledgements.

The author acknowledges the Spanish MINECO grants FIS2015-67411, FIS2017- 91460-EXP, the CAM research consortium QUITEMAD S2018/TCS-4342, and US Army Research Office through grant W911NF-14-1-0103 for partial financial support. He is grateful to Aurelia Chenu and Adolfo del Campo for sharing their insight on this subject during his visit to DIPC. He also thanks them for their kind invitation and hospitality. \conflictsofinterestThe author declares no conflict of interest. \appendixtitlesyes

Appendix A Internal Energy in the Weak Coupling Limit

Consider the projection operators $\mathcal{P}\rho_{\mathscr{SE}}=\mathrm{Tr}_{\mathscr{E}}(\rho_{\SS})\otimes\rho_{\mathscr{E}}^{\beta}$, and $\mathcal{Q}=\mathds{1}-\mathcal{P}$. By applying $\mathcal{Q}$ on the the von Neumann equation (1):

[TABLE]

where we have inserted the identity between $\tilde{\mathcal{V}}(t)$ and $\tilde{\rho}_{\mathscr{SE}}(t)$. The formal integration of this equation gives BrPr02 ; Libro :

[TABLE]

with

[TABLE]

For $\rho_{\mathscr{SE}}(0)=\rho_{\SS}(0)\otimes\rho_{\mathscr{E}}^{\beta}$, $\mathcal{Q}\rho_{\mathscr{SE}}(0)=0$ and the first term of the above equation vanishes:

[TABLE]

On the other hand, in the weak coupling limit:

[TABLE]

By introducing the projection operators in the mean value of $H_{\SS}+\alpha V_{\mathscr{SE}}$ (in the interaction picture), we obtain

[TABLE]

The last term in this equation is simplified because $\mathrm{Tr}[A\mathcal{Q}\tilde{\rho}_{\mathscr{SE}}(t)]=0$ for any system operator $A$ and using (95):

[TABLE]

where we have taken the coupling constant $\alpha$ outside of $\tilde{\mathcal{V}}(s)$, and after (96),

[TABLE]

Furthermore, since \mathrm{Tr}\big{[}\tilde{V}_{\mathscr{SE}}\rho_{\mathscr{SE}}^{\beta}\big{]}=0, \mathrm{Tr}\big{[}\tilde{V}_{\mathscr{SE}}(t)\mathcal{P}\tilde{\rho}_{\mathscr{SE}}(t)\big{]}=0, and Eq. (A) yields

[TABLE]

hence

[TABLE]

Appendix B Calculation of the Time-Ordered Exponential for the Driven Spin-Boson Model in the Refined Weak Coupling

For the spin-boson model with diagonal driving $H_{\SS}(t)=\tfrac{\omega_{0}(t)}{2}\sigma_{z}$, the operator $\mathcal{Z}_{2}(t,s)$ has the same algebraic structure as (83). It is a linear combination of the operators $\mathcal{Z}_{z}:=[\sigma_{z},\cdot]$, $\mathcal{D}_{+-}:=\sigma_{+}(\cdot)\sigma_{-}-\{\sigma_{-}\sigma_{+},\cdot\}/2$, $\mathcal{D}_{-+}:=\sigma_{-}(\cdot)\sigma_{+}-\{\sigma_{+}\sigma_{-},\cdot\}/2$, $\mathcal{D}_{--}:=\sigma_{-}(\cdot)\sigma_{-}$, and $\mathcal{D}_{++}:=\sigma_{+}(\cdot)\sigma_{+}$. As commented in Refined , these operators close a Lie algebra:

[TABLE]

with the zero value for the rest of the cases. Actually, it is clear that this algebra is the direct sum of the subalgebras generated by $\{\mathcal{D}_{+-},\mathcal{D}_{-+}\}$ and $\{\mathcal{Z}_{z},\mathcal{D}_{++},\mathcal{D}_{--}\}$, so that the temporally-ordered exponential in (22) splits into

[TABLE]

The second one leaves invariant any diagonal operator $D$ in the $\sigma_{z}$ eigenbasis, so that from (23)

[TABLE]

Therefore, the application of $\Lambda_{t}^{\rm RA}$ on diagonal states only requires the calculation of the first time-ordered exponential between curly brackets of (B). It is straightforward to prove that the same is true for the Heisenberg adjoint $(\Lambda_{t}^{\rm RA})^{\star}$ acting on diagonal operators.

In addition, the fact that $\{\mathcal{D}_{+-},\mathcal{D}_{-+}\}$ closes a Lie algebra allows for writing the required temporally ordered exponential as a product of exponentials of these operators Puri ,

[TABLE]

In order to compute the numbers $C_{j}(t)$, we introduce a new ones $C_{j}(t,r)$ for the decomposition

[TABLE]

such that $C_{j}(t,t)=C_{j}(t)$. Differentiating with respect to $t$ and multiplying by the inverse (B5) at both sides yields

[TABLE]

The computation of the exponential products on the right-hand side, and a comparison between the coefficients of $\mathcal{D}_{+-}$ and $\mathcal{D}_{-+}$ on both sides leads to a pair of first order non-linear differential equations,

[TABLE]

These equations are easy to solve. Under the initial conditions $C_{j}(0,r)=0$, we obtain the solution

[TABLE]

Finally, the coefficients $C_{j}(t)$ are

[TABLE]

which allows us to compute the time-ordered exponential at (B4) by composing ordinary matrix exponentials.

\reftitle

References