# Extraction of Work via a Thermalization Protocol

Nicolò Piccione, Benedetto Militello, Anna Napoli, Bruno Bellomo

► To cite this version:

Nicolò Piccione, Benedetto Militello, Anna Napoli, Bruno Bellomo. Extraction of Work via a Thermalization Protocol. 11th Italian Quantum Information Science Conference, Sep 2018, Catania, Italy. pp.22, 10.3390/proceedings2019012022. [https://doi.org/10.3390/proceedings2019012022](https://doi.org/10.3390/proceedings2019012022). [hal-02418405](https://HAL.inistre.fr/hal-02418405)

HAL Id: hal-02418405

https://hal.science/hal-02418405v1

Submitted on 22 Jul 2025

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
<figure><img src="image_2.png" /></figure>

# Extraction of Work via a Thermalization Protocol †

Nicolò Piccione 1,*, Benedetto Militello 2,3, Anna Napoli 2,3 and Bruno Bellomo 1,*

1 Institut UTINAM - UMR 6213, CNRS, Université Bourgogne Franche-Comté, Observatoire des Sciences de l'Univers THETA, 41 bis Avenue de l'Observatoire, F-25010 Besançon, France; nicolo.piccione@univ-fcomte.fr

2 Dipartimento di Fisica e Chimica - Emilio Segrè, Università degli Studi di Palermo, via Archirafi 36, I-90123 Palermo, Italy; benedetto.militello@unipa.it (B.M.); anna.napoli@unipa.it (A.N.)

3 INFN Sezione di Catania, via Santa Sofia 64, I-95123 Catania, Italy

* Correspondence: bruno.bellomo@univ-fcomte.fr

† Presented at the 11th Italian Quantum Information Science conference (IQIS2018), Catania, Italy, 17–20 September 2018.

Published: 4 July 2019

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

**Abstract:** This extended abstract contains an outline of the work reported at the conference IQIS2018. We show that it is possible to exploit a thermalization process to extract work from a resource system $R$ to a bipartite system $S$. To do this, we propose a simple protocol in a general setting in the presence of a single bath at temperature $T$ and then examine it when $S$ is described by the quantum Rabi model at $T = 0$. We find the theoretical bounds of the protocol in the general case and we show that when applied to the Rabi model it gives rise to a satisfactory extraction of work and efficiency.

**Keywords:** quantum thermodynamics; extraction of work; quantum Rabi model

## 1. Introduction

In this extended abstract, based on Ref. [1], we present a work-extraction protocol which exploits a simple thermalization process as the main ingredient. To this purpose, it is necessary to define the quantity work and the meaning of “work extraction”. Defining work is a difficult task in quantum mechanics and a general definition has not yet been adopted [2]. We base our protocol on a definition [3] given inside the theoretical framework of thermodynamic resource theory. This deals with operations performed on a system and a thermal bath, by imposing that the sum of the energies of the system and the bath is conserved through the process (not only on average but also in the entire distribution) [4]. This recently developed theoretical framework provided interesting results such as generalizations of the second law of thermodynamics [5,6] and a way to exploit coherences in the energy basis [7].

In the definition of Ref. [3], we need to fix the temperature of the thermal bath and divide the system of interest in two systems. We name one of them $R$, the resource from which we want to extract the work, and the other one $S$, the system on which we want to do work. The work done or extracted is quantified by the difference of a quantity, which is inherent to the state of the system $S$, between the start and the end of the process. The work quantifier can be chosen among different quantities and its choice depends on the physical constraints that we consider for the process at hand. Furthermore, whatever the chosen quantity is, the amount of work extracted from system $S$ has to be always equal or lower than the amount of work lost by system $R$. In Ref. [3] there are a lot of possible quantifiers of work done on (or extracted from) a system, among which we choose to adopt the following one:

$$W = \Delta F(\rho', H') - \Delta F(\rho, H), \quad (1)$$
<figure><img src="image_4.png" /><figcaption>Figure 1. This figure illustrates the phases of the thermalization protocol, as described in the text.</figcaption></figure>

where

$$ \Delta F(\rho, H) = F(\rho, H) - F(\rho^{\text{th}}, H), \quad F(\rho, H) = \text{Tr}(H\rho) - k_B T S(\rho) \quad \text{and} \quad S(\rho) = -\text{Tr}\{\rho \ln(\rho)\}. \qquad (2) $$

Here, $F(\rho, H)$ is the free energy of the state $\rho$ when the system is governed by the Hamiltonian $H$, $\rho^{\text{th}}$ is the thermal state of the system at temperature $T$ equal to the temperature of the thermal bath which is used in the process, $k_B$ is the Boltzmann's constant and $S(\rho)$ is the Von Neumann entropy of the state $\rho$. The quantities which have the symbol ′ are connected to the end of the process, while those not marked are connected to the start of the process. If $H' = H$, Equation (1) reduces to the simpler form

$$ W = F(\rho', H) - F(\rho, H). \qquad (3) $$

Contrarily to other definitions of work, this one is not based on measurements.

## 2. Results and Discussion

### 2.1. The Thermalization Protocol

The main idea of our protocol is to exploit the difference between the thermal states of a bipartite system in two different cases: when the interaction between the two subsystems is turned off and when it is not. When there is no interaction, the global thermal state is the tensor product of the two local thermal states. On the other hand, when the interaction is turned on, the global thermal state is more complex. In particular, the local state of each subsystem, named reduced thermal state, is different from the original local thermal state and has a free energy higher than the thermal one [3,4] when both are computed with respect to the same free Hamiltonian of the subsystem. This translates into formulas as follows:

$$ \rho_A^{\text{th}} \otimes \rho_B^{\text{th}} \rightarrow \rho_S^{\text{th}} \neq \rho_A^{\text{th}} \otimes \rho_B^{\text{th}} \quad \text{and} \quad F(\rho_{A(B)}^{\text{th}}, H_{A(B)}) > F(\rho_{A(B)}^{\text{th}}, H_{A(B)}), \qquad (4) $$

where we named $A$ and $B$ the two subsystems of $S$, $H_{A(B)}$ is the free Hamiltonian of subsystem $A(B)$, $\rho_{A(B)}^{\text{th}}$ is the corresponding local thermal state and $\rho_{A(B)}^{\text{rth}} = \text{Tr}_{B(A)}\{\rho_S^{\text{th}}\}$ is its reduced thermal state and $\rho_S^{\text{th}}$ is computed with respect to the Hamiltonian of system $S, H_A + H_B + H_I$, being $H_I$ the interaction term between $A$ and $B$. The protocol is composed of the following main steps (see Figure 1 and find more information in Ref. [1]).

* At the start of the protocol the two subsystems of $S$ are in their local thermal states and they are non-interacting ($t_1$).

* Then, the resource $R$ turns on the interaction between them while not changing their states (in a small time interval preceding $t_2$, represented by the small box before $t_2$).

* The two subsystems of $S$, in contact with the thermal bath at temperature $T$, thermalize together to the new thermal state of $S$ (from $t_2$ to $t_3$).

* Lastly, the resource $R$ turns off the interaction (in a small time interval following $t_3$) while not changing the states of the subsystems of $S$.
If one follows the evolution of the free energy of the subsystems of S, he finds that at the end of the protocol (once the interaction is off) the extracted work, W, satisfies [1]

$$W = \Delta F(\rho_A^{\text{rth}}, H_A) + \Delta F(\rho_B^{\text{rth}}, H_B) + k_B TS(A:B) \leq \langle H_I \rangle_{t_2} - \langle H_I \rangle_{t_3}, \quad (5)$$

where $S(A:B) = S(\rho_A^{\text{rth}}) + S(\rho_B^{\text{rth}}) - S(\rho_S^{\text{th}})$ is the mutual information of the A and B subsystems. Moreover, the average of $H_I$ is calculated as soon as the interaction is on ($t_2$) and right before it is turned off ($t_3$). We note that the term which contains the mutual information is a non-local term which, in some cases, could then be neglected because of the difficulty of its exploitation.

It is also possible to define an ideal efficiency of this process by taking into account the minimum amount of work that the system R has lost. Then, the efficiency is given by the ratio of the extracted work and the work lost by R:

$$\eta = \frac{W}{\langle H_I \rangle_{t_2} - \langle H_I \rangle_{t_3}} \leq 1. \quad (6)$$

Once the thermalization protocol is completed, one can use the resulting states to perform other tasks of various kind. For example, one could want to charge an external battery. In this case one needs a transfer protocol to move the work extracted by S to the battery so that S can be reinitialized and can begin another cycle of the thermalization protocol. We studied this issue in the specific case of the application of the thermalization protocol to a system S governed by the Rabi Hamiltonian at zero temperature. We will include this study in the revised version of Ref. [1].

We now apply the thermalization protocol to the Rabi model at zero temperature.

## 2.2. Application to the Rabi Model at Zero Temperature

The Rabi model describes the interaction between a two-level system and a quantum harmonic oscillator. Its Hamiltonian has the following form [8]:

$$H_{\text{Rb}} = \frac{\hbar\omega}{2}\sigma_z + \hbar\omega\hat{n} + \hbar g\sigma_x(a^\dagger + a), \quad (7)$$

where $\sigma_z$ and $\sigma_x$ are the Pauli operators of the two-level system, $a(a^\dagger)$ is the annihilation (creation) operator of the harmonic oscillator, $\hat{n} = a^\dagger a$, $\omega$ is the frequency of both the two-level system and the harmonic oscillator (they are resonant) and $g$ is the coupling strength between the two subsystems.

In this model, the extracted work and the ideal efficiency take the following form [1]:

$$W = \frac{\hbar\omega}{2} (1 + \langle\sigma_z\rangle_{t_3}) + \hbar\omega \langle\hat{n}\rangle_{t_3} = \hbar v_0 + \frac{\hbar\omega}{2} - \hbar g \langle\sigma_x (a^\dagger + a)\rangle_{t_3} \quad \text{and} \quad \eta = -\frac{W}{\hbar g \langle\sigma_x (a^\dagger + a)\rangle_{t_3}}, \quad (8)$$

where $v_0$ is the frequency of the ground state of the system when the interaction is on. Both $\hbar v_0$ and $\hbar g \langle\sigma_x (a^\dagger + a)\rangle_{t_3}$ are negative quantities. Also, it is worth noting that the non-local term $S(A:B)$ does not contribute to the work extraction at zero temperature.

The Rabi Hamiltonian has been recently exactly diagonalized [8,9]. This let us write its ground state in the bare basis and, thus, calculate the extracted work and the efficiency of the protocol when applied to the Rabi model at zero temperature. The graph of Figure 2 shows both the amount of extracted work and the ideal efficiency as a function of the coupling parameter in units of $\omega$. As one can see, the extracted work is comparable with the typical energies of the two subsystems when $g/\omega$ is sufficiently high. We remark that these high values of the coupling constant became recently experimentally attainable [10]. Moreover, the efficiency of the thermalization protocol is always higher than one half, with a peak of roughly 0.60.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** This graph shows the extracted work (solid black line), $W/(\hbar\omega)$, and the ideal efficiency (dashed blue line), $\eta$, as a function of the coupling parameter $g/\omega$, for the Rabi model at $T = 0$.</figcaption></figure>

These results show that the protocol we propose can achieve good results when applied to realistic physical models such as the Rabi model.

**Author Contributions:** N.P. performed the main part of the investigation and prepared a first version of the manuscript, mainly under the supervision of B.B. All authors contributed to the conceptualization and the analysis.

**Acknowledgments:** N.P. acknowledges the financial support of the Erasmus+ programme of the European Union and of the UTINAM Institute for the development of this programme. B.B. acknowledges the financial support of the Observatoire des Sciences de l'Univers THETA Franche-Comté / Bourgogne for his participation to the conference IQIS2018, including publication charges.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Piccione, N.; Militello, B.; Napoli, A.; Bellomo, B. Work extraction exploiting thermalization with a single bath. *arXiv* 2018, arXiv:1806.11384.

2. Perarnau-Llobet, M.; Bäumer, E.; Hovhannisyan, K.V.; Huber, M.; Acin, A. No-Go Theorem for the Characterization of Work Fluctuations in Coherent Quantum Systems. *Phys. Rev. Lett.* **2017**, *118*, 070601, doi:10.1103/PhysRevLett.118.070601.

3. Gallego, R.; Eisert, J.; Wilming, H. Thermodynamic work from operational principles. *New J. Phys.* **2016**, *18*, 103017, doi:10.1088/1367-2630/18/10/103017.

4. Horodecki, M.; Oppenheim, J. Fundamental limitations for quantum and nanoscale thermodynamics. *Nat. Commun.* **2013**, *4*, 2059, doi:10.1038/ncomms3059.

5. Brandão, F.; Horodecki, M.; Ng, N.H.Y.; Oppenheim, J.; Wehner, S. The second laws of quantum thermodynamics. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 3275-3279, doi:10.1073/pnas.1411728112.

6. Lostaglio, M.; Jennings, D.; Rudolph, T. Description of quantum coherence in thermodynamic processes requires constraints beyond free energy. *Nat. Commun.* **2015**, *6*, 6383, doi:10.1038/ncomms7383.

7. Åberg, J. Catalytic coherence. *Phys. Rev. Lett.* **2014**, *15*, 150402, doi:10.1103/PhysRevLett.113.150402.

8. Braak, D. Integrability of the Rabi model. *Phys. Rev. Lett.* **2011**, *107*, 100401, doi:10.1103/PhysRevLett.107.100401.

9. Chen, Q.H.; Wang, C.; He, S.; Liu, T.; Wang, K.L. Exact solvability of the quantum Rabi model using Bogoliubov operators. *Phys. Rev. A* **2012**, *86*, 023822, doi:10.1103/PhysRevA.86.023822.

10. Yoshihara, F.; Fuse, T.; Ashhab, S.; Kakuyanagi, K.; Saito, S.; Semba, K. Superconducting qubit-oscillator circuit beyond the ultrastrong-coupling regime. *Nat. Phys.* **2017**, *13*, 44-47, doi:10.1038/nphys3906.

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