Efficiency Fluctuations in a Quantum Battery Charged by a Repeated Interaction Process

Abstract

A repeated interaction process assisted by auxiliary thermal systems charges a quantum battery. The charging energy is supplied by switching on and off the interaction between the battery and the thermal systems. The charged state is an equilibrium state for the repeated interaction process, and the ergotropy characterizes its charge. The working cycle consists in extracting the ergotropy and charging the battery again. We discuss the fluctuating efficiency of the process, among other fluctuating properties. These fluctuations are dominated by the equilibrium distribution and depend weakly on other process properties.

1. Introduction

Repeated interaction schemes, also known as collisional models [1,2,3,4,5,6], have played a vital role in the development of quantum optics [7,8,9,10] and the rapid evolution of quantum thermodynamics [11,12,13,14,15]. The idealized and straightforward formalism has been crucial to designing and understanding quantum devices such as information engines [16,17,18,19], heat engines [12,20,21,22,23], and quantum batteries [24,25,26,27,28,29,30,31,32,33,34]. Recently, it was realized that the framework can be extended to deal with macroscopic reservoirs [23,35], expanding the reach of applications in quantum thermodynamics. For comprehensive reviews of the method and its applications, see [36,37].

In the simplest scenario, many copies of an auxiliary system in the Gibbs equilibrium thermal state interact sequentially with a system of interest. Each interaction step is described by a completely positive trace-preserving (CPTP) map [38]. The repeated interaction process corresponds to concatenations of the map, which eventually will bring the system to a nonequilibrium steady state or an equilibrium state. In equilibrium, heat does not flow to the environment, and entropy is not produced. When the repeated interaction brings the system to an equilibrium state, we say that we iterate a map with equilibrium. In this paper, we apply this framework to study a quantum battery.

Quantum technologies, such as quantum computing, communication, and sensing, are supported by the quantum storage and transfer of energy. Implementing fast and reliable quantum batteries in these technologies may improve their functionality. Different quantum batteries have been proposed to achieve these goals [39,40,41,42]. One paradigmatic setup considers the battery to be composed of noninteracting qubits. Global operations, such as charging or discharging the battery by coupling all qubits to a single optical cavity mode, boost its performance in power [28,29,30,31] and reliability [43].

The most straightforward repeated interaction model for a quantum battery considers nonequilibrium auxiliary systems supplying the energy. However, the process of sustaining the charged state is dissipative. Reference [26] proposed a different kind of quantum battery where the charged state corresponds to the equilibrium state of the process. The work in the recharging stage provides the energy, which is preserved without dissipation in the equilibrium state as long as the battery–environment interaction remains under control. In actual physical implementations, other exchanges can still cause energy leakage. The battery’s charge is characterized by its ergotropy [44], i.e., the maximum amount of energy extracted with a unitary process. Once removed, a repeated interaction process recharges the battery. In this way, we have a thermodynamic cycle.

The recharging energy and the ergotropy delivered by the quantum battery are averaged values that are relevant for several cycles or many batteries working parallel. In a single cycle, one can observe fluctuations when observing these energies. Therefore, their study is relevant for the reliability of the device. The two-point measurement scheme [45] is appropriate for describing these thermal and quantum fluctuations that reveal essential properties of the process [46,47,48,49]. Other sources of randomness in the operation of a battery can arise from changes in the evolution operator [50,51,52], Hamiltonian [53], and initial condition [54]. We do not take them into account. Closer to the spirit of this work are studies of work fluctuations in the charging or discharging process of isolated quantum batteries [55,56,57].

Thus, in this work, we take the dissipative quantum battery [26] and study fluctuations in the thermodynamic quantities such as heat and work during the charging phase and the efficiency fluctuations of the cycle. Efficiency fluctuations are significant in assessing the performance of a machine. They have drawn recent attention in classical [58,59,60,61,62,63,64,65,66,67,68] and quantum [21] engines. Evaluating the fluctuations requires detailed information about the bath and the process [45]. However, a key simplification arises because we deal with maps with equilibrium, allowing us to determine the statistics of the fluctuations. We will illustrate this using two examples.

For completeness, we also consider equilibrium fluctuations. We evaluate the probability of performing work or absorbing heat while keeping the (average) charge in the battery. We compare our findings with the equilibrium fluctuations in a process with a Gibbs equilibrium state.

The remainder of this article is organized as follows. In Section 2, we review the thermodynamics for CPTP maps, emphasizing the results for maps with equilibrium. Then, in Section 3, we introduce our system of study, namely the equilibrium quantum battery proposed in [26]. Section 4 discusses the stochastic versions of the thermodynamic equalities and laws, emphasizing the results for maps with equilibrium again. Subsequently, in Section 5, we evaluate these fluctuations in two illustrative examples. We conclude this article in Section 6.

2. Thermodynamic Description for Completely Positive Trace-Preserving Maps

Consider a system S and a system A that jointly evolve under the unitary $U=e^{-i\frac{\tau }{\hslash }(H_S+H_A+V)}$. The Hamiltonians $H_S$ and $H_A$ of S and A, respectively, are constant in time. The coupling between S and A during the time interval $(0,\tau )$ is given by the interaction energy V and vanishes for $t<0$ and $t>\tau$.

Initially, S and A are uncorrelated; i.e., their density matrix is the tensor product of the respective density matrices ${\rho }_{\mathrm{tot}}={\rho }_S\otimes {\omega }_{\beta }\left(H_A\right)$, where ${\omega }_{\beta }\left(H_A\right)=\frac{e^{-\beta H_A}}{Z_A}$ is the Gibbs thermal state for A with $\beta$ as the inverse temperature, and $Z_A=\mathrm{Tr}{e}^{-\beta H_A}$. After the lapse of time $\tau$, the initial state ${\rho }_{\mathrm{tot}}$ changes to a new state,

$${\rho }_{\mathrm{tot}}^{\prime }=U\left({\rho }_S\otimes {\omega }_{\beta }\left(H_A\right)\right)U^{†}.$$

(1)

In the following, we denote ${\rho }_S^{\prime }={\mathrm{Tr}}_A{\rho }_{\mathrm{tot}}^{\prime }$ and ${\rho }_A^{\prime }={\mathrm{Tr}}_S{\rho }_{\mathrm{tot}}^{\prime }$, where ${\mathrm{Tr}}_X$ is the partial trace over subsystem X. By tracing out A, one obtains a CPTP map $\mathcal{E}$ for the system S evolution

$${\rho }_S^{\prime }=\mathcal{E}\left({\rho }_S\right)={\mathrm{Tr}}_A\left[U\left({\rho }_S\otimes {\omega }_{\beta }\left(H_A\right)\right)U^{†}\right].$$

(2)

The energy change of S

$$\Delta E=\mathrm{Tr}\left[H_S({\rho }_S^{\prime }-{\rho }_S)\right],$$

(3)

can be written as the sum of

$$Q=\mathrm{Tr}\left[H_A({\omega }_{\beta }\left(H_A\right)-{\rho }_A^{\prime })\right],$$

(4)

and

$$W=\mathrm{Tr}\left[(H_S+H_A)({\rho }_{\mathrm{tot}}^{\prime }-{\rho }_{\mathrm{tot}})\right],$$

(5)

satisfying the first law $\Delta E=W+Q$. Note that $-Q$ is the energy change of A, we call Q the heat, and W is the energy change of the full $S+A$ system, which we call the switching work because it is due to the energy cost of turning on and off the interaction V at the beginning and end of the process, respectively [69,70].

Consider the von Neumann entropy change

$$\Delta S_{\mathrm{vN}}=-\mathrm{Tr}[{\rho }_S^{\prime }ln{\rho }_S^{\prime }]+\mathrm{Tr}[{\rho }_{S}ln{\rho }_S]$$

(6)

of system S and the heat Q given in Equation (4). The entropy production, $\Sigma =\Delta S_{\mathrm{vN}}-\beta Q$, is also given by [71]

$$\Sigma =D\left({\rho }_{\mathrm{tot}}^{\prime }\right||{\rho }_S^{\prime }\otimes {\omega }_{\beta }\left(H_A\right))\ge 0,$$

(7)

with $D\left(a\right|\left|b\right)\equiv Tr[alna]-Tr[alnb].$ The inequality in Equation (7) corresponds to the second law. Note that auxiliary system A does not need to be macroscopic; nevertheless, we will call it the bath.

As in standard thermodynamics, analyzing the process ${\rho }_S\to {\rho }_S^{\prime }=\mathcal{E}\left({\rho }_S\right)$, in terms of $\Delta E=W+Q$ and $\Sigma =\Delta S_{\mathrm{vN}}-\beta Q\ge 0$ with the quantities given in Equations (3)–(7) is very useful. Note that for their evaluation, particularly for the work, Equation (5), and entropy production, Equation (7), we need to know the full state ${\rho }_{\mathrm{tot}}^{\prime }$.

Maps with Thermodynamic Equilibrium

In a repeated interaction process, one concatenates L CPTP maps ${\mathcal{E}}^L\equiv \mathcal{E}\circ \cdots \circ \mathcal{E}(·)$ to describe a sequence of evolutions of a system coupled to an auxiliary thermal system for a given lapse of time $\tau$. With each map $\mathcal{E}$, a new fresh bath is introduced that exchanges heat with the system during the time that the interaction is turned on. The concatenated map ${\mathcal{E}}^L$ is also a CPTP map. The total work performed is the sum of the work performed by switching on and off the interaction energy with each bath. Similarly, the total heat is the sum of the heat exchanged with each bath.

Let us assume that the map $\mathcal{E}$ has an attractive invariant state $\overline{\rho }$, defined as

$$\underset{L\to \infty }{lim}{\mathcal{E}}^L\left({\rho }_S\right)=\overline{\rho },\forall {\rho }_S,$$

and $\overline{\rho }=\mathcal{E}\left(\overline{\rho }\right)$. The process $\overline{\rho }\to \mathcal{E}\left(\overline{\rho }\right)$ is thermodynamically characterized by $\Delta S_{\mathrm{vN}}=0=\Delta E$; see Equations (3) and (6). If the entropy produced by the action of the map $\mathcal{E}$ on $\overline{\rho }$ is $\Sigma >0$, then we say that the invariant state is a nonequilibrium steady state. The invariant state is an equilibrium state if $\Sigma =0$, i.e., if the entropy production, Equation (7), vanishes by the action of $\mathcal{E}$ on $\overline{\rho }$. Maps with these particular states are called maps with equilibrium [72,73].

According to Equation (7), $\Sigma =0$ for the steady state $\overline{\rho }$ if and only if $\overline{\rho }\otimes {\omega }_{\beta }\left(H_A\right)=U\left(\overline{\rho }\otimes {\omega }_{\beta }\left(H_A\right)\right)U^{†}$. Equivalently, if the unitary U in Equation (1) satisfies $[U,H_0+H_A]=0$, where $H_0$ is an operator in the Hilbert space of the system, then the product state ${\omega }_{\beta }\left(H_0\right)\otimes {\omega }_{\beta }\left(H_A\right)$, with ${\omega }_{\beta }\left(H_0\right)=\frac{e^{-\beta H_0}}{Z_0}$, where $Z_0=\mathrm{Tr}\left[e^{-\beta H_0}\right]$, is invariant under the unitary evolution in Equation (1) and $\overline{\rho }={\omega }_{\beta }\left(H_0\right)$ is an equilibrium state for the map in Equation (2).

It follows from $[U,H_0+H_A]=0$ that the heat, Equation (4), and work, Equation (5), simplify to

$$Q=\mathrm{Tr}\left[H_0({\rho }_S^{\prime }-{\rho }_S)\right]$$

(8)

and

$$W={\mathrm{Tr}}_S\left[(H_S-H_0)({\rho }_S^{\prime }-{\rho }_S)\right].$$

(9)

The entropy production also reduces to an expression that does not involve the state of the bath. Indeed, we obtain

$$\Sigma =D({\rho }_S\left|\right|{\omega }_{\beta }\left(H_0\right))-D({\rho }_S^{\prime }\left|\right|{\omega }_{\beta }\left(H_0\right)),$$

(10)

which is positive due to the contracting character of the map [38]. The averaged thermodynamic quantities for a map with equilibrium are only determined by the properties of the system of interest.

If $H_0=H_S$, then the map is called thermal [74,75]. The equilibrium state is the Gibbs state ${\omega }_{\beta }\left(H_S\right)=e^{-\beta H_S}/Z_S$ with $Z_S=\mathrm{Tr}\left[e^{-\beta H_S}\right]$, and the agent is passive because $W=0$ for every initial state ${\rho }_S$; see Equation (9).

When $H_0\ne H_S$, an active external agent has to provide (or extract) work to perform the map on a state ${\rho }_S$. However, once the system reaches the equilibrium state ${\omega }_{\beta }\left(H_0\right)$, the process ${\omega }_{\beta }\left(H_0\right)\to \mathcal{E}\left({\omega }_{\beta }\left(H_0\right)\right)={\omega }_{\beta }\left(H_0\right)$ is performed with $W=0$; see Equation (9), and $\Sigma =0$.

Let us end this section with the following remark. Since the total evolution operator $U=e^{-i\frac{\tau }{\hslash }(H_S+H_A+V)}$ is time-independent, the equilibrium condition is satisfied by finding $H_0$ and V such that $[H_0,H_S]=0$ and $[H_0+H_A,V]=0$ [26]. In this case, $H_S$ and $H_0$ share the same eigenbasis. To simplify the discussion of fluctuations, we consider non-degenerate eigenenergies. We denote the eigensystems as

$$H_S|n\rangle =E_n|n\rangle ,H_0|n\rangle =E_n^0|n\rangle .$$

with increasing order $E_1<E_2<\cdots <E_N$ for the eigenenergies. The eigenvalues $E_n^0$ are not necessarily ordered, but there is always a permutation that we call $\pi$ of $(1,\dots ,N)\to ({\pi }_1,\dots ,{\pi }_N)$ such that $E_{{\pi }_1}^0\le \cdots \le E_{{\pi }_N}^0$.

3. The Battery

As is well known, the Gibbs state ${\omega }_{\beta }\left(H_S\right)$ is passive; i.e., one cannot decrease (extract) its energy with a unitary operation [76,77]. This is not true for the equilibrium state

$${\omega }_{\beta }\left(H_0\right)=\sum _n\frac{e^{-\beta E_n^0}}{Z_0}|n\rangle \langle n|,$$

(11)

if a pair $(j,k)$ exists such that $(E_j-E_k)(E_j^0-E_k^0)<0$. In that case, the unitary operator u with matrix elements $u_{ij}=\langle i\left|u\right|j\rangle ={\delta }_{{\pi }_i,j}$ extracts the ergotropy [44]

$$\mathcal{W}\left[{\omega }_{\beta }\left(H_0\right)\right]=\sum _{n=1}^N(E_{{\pi }_n}-E_n)\frac{e^{-\beta E_{{\pi }_n}^0}}{Z_0}>0,$$

(12)

where $\pi$ is the permutation that orders $E_n^0$ increasingly.

Once the ergotropy is extracted, the system is left in the passive state

$${\sigma }_{{\omega }_{\beta }\left(H_0\right)}=u{\omega }_{\beta }\left(H_0\right)u^{†}=\sum _{n=1}^N\frac{e^{-\beta E_{{\pi }_n}^0}}{Z_0}|n\rangle \langle n|.$$

(13)

An equilibrium quantum battery was proposed in [26] based on that observation. The system is driven by a repeated interaction process described by a map $\mathcal{E}$ with equilibrium ${\omega }_{\beta }\left(H_0\right)$. Once the equilibrium is reached, it is kept with no cost ($W=0$), energy does not leak from it, and the battery’s charge, characterized by the ergotropy $\mathcal{W}\left[{\omega }_{\beta }\left(H_0\right)\right]$, is preserved. Equilibrium states with ergotropy are called active.

The thermodynamic cycle is as follows: The battery starts in the active equilibrium state, and then the ergotropy (12) is extracted, leaving the battery in the passive state (13) from which the repeated interaction process ${lim}_{L\to \infty }{\mathcal{E}}^L\left({\sigma }_{{\omega }_{\beta }\left(H_0\right)}\right)$ recharges it. As a consequence of the second law, the recharging work $W_R={\mathrm{Tr}}_S\left[(H_S-H_0)({\omega }_{\beta }\left(H_0\right)-{\sigma }_{{\omega }_{\beta }\left(H_0\right)})\right]$ is never smaller that the extracted ergotropy. In this way, the thermodynamic efficiency

$$\begin{array}{c}\hfill 0\le {\eta }_{\mathrm{th}}\equiv \frac{\mathcal{W}\left[{\omega }_{\beta }\left(H_0\right)\right]}{W_R}\le 1,\end{array}$$

(14)

which is the ratio of the wanted resource over the invested, characterizes the operation of the device.

4. Fluctuations

4.1. Repeated Interaction for a Map with Equilibrium

The thermodynamic quantities in Equations (3)–(7) were obtained as the average over their stochastic versions defined over trajectories using a two-point measurement scheme in [72]. Since all interesting density matrices ${\omega }_{\beta }\left(H_S\right),{\omega }_{\beta }\left(H_0\right)$, and ${\sigma }_{{\omega }_{\beta }\left(H_0\right)}$ are diagonal in the system energy basis, we need only projective energy measurement in this work.

A trajectory $\gamma =\{n;i_1,j_1,\dots ,i_L,j_L;m\}$ for the recharging process is defined by the initial and final, ${\epsilon }_{i_k}$ and ${\epsilon }_{j_k}$, energy results for each auxiliary thermal system and $E_n$ and $E_m$ for the system. According to the two-point measurement scheme [45], its probability is

$$P_{\gamma }^{\left(L\right)}=|\langle j_1\cdots j_{L}m|U_L\cdots U_1|i_1\cdots i_{L}n{\rangle |}^2\frac{e^{-\beta {\sum }_{k=1}^L{\epsilon }_{i_k}}}{Z_A^L}{p}_{\mathrm{ini}}\left(n\right),$$

(15)

where $p_{\mathrm{ini}}\left(n\right)$ is the probability that the initial state of the system is $|n\rangle$; see Appendix A. We now associate the stochastic thermodynamic quantities with these trajectories. The stochastic heat flow to the system $q_{\gamma }$ corresponds to the negative energy change of the bath, i.e., $q_{\gamma }={\sum }_{k=1}^L({\epsilon }_{i_k}-{\epsilon }_{j_k})$. According to the first law of stochastic thermodynamics [47], the stochastic work is given by

$$w_{\gamma }=\Delta e_{\gamma }-q_{\gamma },$$

(16)

where $\Delta e_{\gamma }=E_m-E_n$ is the stochastic energy change. These fluctuating quantities are studied through their distributions

$$p_w^{\left(L\right)}\left(x\right)=\sum _{\gamma }\delta (x-w_{\gamma })P_{\gamma }^{\left(L\right)},p_{\Delta e}^{\left(L\right)}\left(x\right)=\sum _{\gamma }\delta (x-\Delta e_{\gamma })P_{\gamma }^{\left(L\right)},p_q^{\left(L\right)}\left(x\right)=\sum _{\gamma }\delta (x-q_{\gamma })P_{\gamma }^{\left(L\right)},$$

(17)

and, as for the averaged thermodynamic quantities, we need information on the state of the whole system to evaluate them. However, for maps with equilibrium, a stochastic trajectory is determined by the pair $\gamma =\{n,m\}$; see Appendix A. Consequently these formulas simplify and become, $q_{\gamma }=E_m^0-E_n^0,w_{\gamma }=E_m-E_m^0-(E_n-E_n^0)$ with the distributions

$$p_{\Delta e}^{\left(L\right)}\left(x\right)=\sum _{n,m}\delta (x-[E_m-E_n])P_{n\to m}^{\left(L\right)},$$

(18)

$$p_w^{\left(L\right)}\left(x\right)=\sum _{n,m}\delta (x-[(E_m-E_m^0)-(E_n-E_n^0)])P_{n\to m}^{\left(L\right)},$$

(19)

$$p_q^{\left(L\right)}\left(x\right)=\sum _{n,m}\delta (x-[E_m^0-E_n^0])P_{n\to m}^{\left(L\right)},$$

(20)

and the trajectory probability

$$P_{n\to m}^{\left(L\right)}=\langle m|{\mathcal{E}}^L\left(|n\rangle \langle n|\right)|m\rangle p_{\mathrm{ini}}\left(n\right)={\left(T^L\right)}_{m|n}{p}_{\mathrm{ini}}\left(n\right),$$

(21)

in terms of the initial probability $p_{\mathrm{ini}}\left(n\right)$ and of the L power of the stochastic matrix $T_{m|n}=\langle m|\mathcal{E}\left(|n\rangle \langle n|\right)|m\rangle$.

The averages $\int x{p}_{\Delta e}^{\left(L\right)}\left(x\right)dx,\int x{p}_w^{\left(L\right)}\left(x\right)dx,\int x{p}_q^{\left(L\right)}\left(x\right)dx$ reproduce Equations (3), (8) and (9) with ${\rho }_S^{\prime }={\mathcal{E}}^L\left({\rho }_S\right)$ and ${\rho }_S={\sum }_n{p}_{\mathrm{ini}}\left(n\right)|n\rangle \langle n|$.

4.2. Fluctuations in the Equilibrium State

As noted before, all averaged thermodynamic quantities $\Delta E=\Delta S=\Sigma =W=Q=0$ vanish for a process in equilibrium. So, on average, the process ${\omega }_{\beta }\left(H_0\right)\to \mathcal{E}\left({\omega }_{\beta }\left(H_0\right)\right)={\omega }_{\beta }\left(H_0\right)$ has no energy cost. However, if $H_0\ne H_S$, the agent is still active due to non-vanishing work fluctuations. For thermal maps, $H_0=H_S$ and Equation (19) gives $p_w^{\left(L\right)}\left(x\right)=\delta \left(x\right)$. The external agent is truly passive.

To analyze equilibrium fluctuations, we use Equations (18)–(20) with $p_{\mathrm{ini}}\left(n\right)=\frac{e^{-\beta E_n^0}}{Z_0}$.

4.3. Recharging Process

Since the recharging process starts from ${\sigma }_{{\omega }_{\beta }\left(H_0\right)}$, we take $p_{\mathrm{ini}}\left(n\right)=e^{-\beta E_{{\pi }_n}^0}/Z_0$—see Equation (13)—in the distribution Equations (18)–(20).

Since the charged state ${\omega }_{\beta }\left(H_0\right)$ is reached asymptotically, we take $L\to \infty$ to charge the battery fully.

Moreover, since $\mathcal{E}$ has a unique equilibrium state, we will find that T is a regular stochastic matrix [78], implying that ${lim}_{L\to \infty }{\left(T^L\right)}_{m|n}=e^{-\beta E_m^0}/Z_0,\forall n$. Therefore, the limit in Equation (21)

$$P_{n\to m}^{(\infty )}=p_{\mathrm{ini}}\left(n\right)e^{-\beta E_m^0}/Z_0=e^{-\beta (E_{{\pi }_n}^0+E_m^0)}/Z_0^2,$$

(22)

is independent of the map’s details. Interestingly, the rate of convergence of $T^L$ to the equilibrium distribution depends on the map $\mathcal{E}$ parameters. We later discuss the fluctuations of a concatenated process ${\mathcal{E}}^L$ with finite L.

The average of the stochastic energy change in the recharging process

$${\langle \Delta e_{\gamma }\rangle }^{(\infty )}\equiv \sum _{n,m}(E_m-E_n)P_{n\to m}^{(\infty )}=Tr\left[H_S({\omega }_{\beta }\left(H_0\right)-{\sigma }_{{\omega }_{\beta }\left(H_0\right)})\right]=\mathcal{W}\left({\omega }_{\beta }\left(H_0\right)\right)$$

(23)

is the ergotropy. The average stochastic work

$${\langle w_{\gamma }\rangle }^{(\infty )}\equiv \sum _{n,m}((E_m-E_m^0)-(E_n-E_n^0))P_{n\to m}^{(\infty )}=Tr\left[(H_S-H_0)({\omega }_{\beta }\left(H_0\right)-{\sigma }_{{\omega }_{\beta }\left(H_0\right)})\right]=W_R$$

(24)

is the recharging work.

4.4. Extracting Process

The extracting process also fluctuates when we measure the battery’s energy in the charged state and the discharged state. We call $\kappa$ the stochastic trajectory in the ergotropy extracting process and ${\varpi }_{\kappa }$ the stochastic extracted energy. The probability $p_{\kappa }$ of $\kappa =(m^{\prime },n)$ is the product of the transition probability from $|m^{\prime }\rangle$ to $|n\rangle$ under the permutation u, $P_{m^{\prime }\to n}^{\mathrm{ext}}={\left|\langle n\left|u\right|m^{\prime }\rangle \right|}^2={\delta }_{{\pi }_n,m^{\prime }}$, with the initial probability $e^{-\beta E_{m^{\prime }}^0}/Z_0$; see Equation (11). The averaged extracted energy,

$$\langle {\varpi }_{\kappa }\rangle =\sum _{\kappa }{\varpi }_{\kappa }{p}_{\kappa }=\sum _{m^{\prime },n}(E_{m^{\prime }}-E_n)P_{m^{\prime }\to n}^{\mathrm{ext}}\frac{e^{-\beta E_{m^{\prime }}^0}}{Z_0}=\sum _n(E_{{\pi }_n}-E_n)\frac{e^{-\beta E_{{\pi }_n}^0}}{Z_0}=\mathcal{W}\left({\omega }_{\beta }\left(H_0\right)\right)$$

(25)

is the ergotropy Equation (12).

Equations (23) and (25) show the cycle’s consistency, where two processes, recharging ($\gamma$) and extracting ($\kappa$), connect the same states, ${\omega }_{\beta }\left(H_0\right)$ and ${\sigma }_{{\omega }_{\beta }\left(H_0\right)}$.

4.5. Fluctuating Efficiency for the Cycle

In terms of Equations (24) and (25), we have the thermodynamic efficiency ${\eta }_{\mathrm{th}}=\frac{\mathcal{W}}{W_R}=\frac{\langle {\varpi }_{\kappa }\rangle }{{\langle w_{\gamma }\rangle }^{(\infty )}}$.

As the thermodynamic efficiency is the ratio of the ergotropy over the recharging work, the fluctuating efficiency [21] should be the ratio of their fluctuating equivalents. The fluctuating extracted energy is ${\varpi }_{\kappa }=E_{m^{\prime }}-E_n$, and the fluctuating work is $w_{\gamma }=E_m-E_m^0-(E_n-E_n^0)$. Therefore, we define the fluctuating efficiency as

$${\eta }_{\gamma \kappa }=\frac{{\varpi }_{\kappa }}{w_{\gamma }}=\frac{E_{m^{\prime }}-E_n}{E_m-E_m^0-(E_n-E_n^0)}.$$

(26)

Given the extracting trajectory $\kappa$, the probability of the recharging trajectory $\gamma$ is $P_{m^{\prime }\to n}^{\mathrm{ext}}{P}_{n\to m}^{\infty }$. Thus, the joint probability for the processes $\kappa$ and $\gamma$ is

$$p_{\gamma \kappa }=\frac{e^{-\beta E_{m^{\prime }}^0}}{Z_0}{P}_{m^{\prime }\to n}^{\mathrm{ext}}{P}_{n\to m}^{\infty }=\frac{e^{-\beta E_{m^{\prime }}^0}}{Z_0}{\delta }_{{\pi }_n,m^{\prime }}\frac{e^{-\beta E_m^0}}{Z_0},$$

and the distribution of the fluctuating efficiency is

$$p_{\eta }\left(x\right)=\sum _{\gamma ,\kappa }\delta (x-{\eta }_{\gamma \kappa })p_{\gamma \kappa }=\sum _{n,m}\delta \left(x-\frac{E_{{\pi }_n}-E_n}{E_m-E_m^0-(E_n-E_n^0)}\right)\frac{e^{-\beta (E_m^0+E_{{\pi }_n}^0)}}{Z_0^2}.$$

(27)

To simplify the notation, we write this as

$$p_{\eta }\left(x\right)=\sum _{n,m}\delta \left(x-{\eta }_{nm}\right)P_{n\to m},$$

(28)

with

$${\eta }_{nm}=\frac{E_{{\pi }_n}-E_n}{E_m-E_m^0-(E_n-E_n^0)},\mathrm{and}{P}_{n\to m}=\frac{e^{-\beta (E_m^0+E_{{\pi }_n}^0)}}{Z_0^2}.$$

(29)

The probability $P_{n\to m}$ corresponds Equation (22), and we omit the superscript.

Trajectories with $w_{\gamma }=0$ and ${\varpi }_{\kappa }\ne 0$ have $|{\eta }_{\gamma \kappa }|=\infty$. Therefore, the average $\langle {\eta }_{\gamma \kappa }\rangle$ does not always exist, and if it does, ${\eta }_{\mathrm{th}}\ne \langle {\eta }_{\gamma \kappa }\rangle$, unless the stochastic work and efficiency are uncorrelated. In fact, $\langle {\eta }_{\gamma \kappa }{w}_{\gamma }\rangle =\langle {\varpi }_{\kappa }\rangle =\mathcal{W}$. So only if $\langle {\eta }_{\gamma \kappa }{w}_{\gamma }\rangle =\langle {\eta }_{\gamma \kappa }\rangle W_R$ do we have $\langle {\eta }_{\gamma \kappa }\rangle ={\eta }_{\mathrm{th}}$. The thermodynamic and fluctuating efficiency can be very different.

The following section discusses efficiency fluctuations for the cycle, heat and work fluctuations for the recharging process and equilibrium fluctuations in two examples.

5. Examples

We illustrate our results in two simple examples. The first example is a single-qubit battery that we use to discuss equilibrium fluctuations (Section 4.2). The second example is a two-qubit battery where we compute heath and work distributions in a partial recharging process (Section 4.3). In both, we compute the fluctuating efficiency distribution (Section 4.5).

5.1. Single-Qubit Battery

An interesting protocol, with $H_0=-H_S$, was discussed in [26] for a system S interacting with systems A, which are copies of S. The corresponding process $\mathcal{E}$ has the remarkable equilibrium state

$${\omega }_{\beta }(-H_S)=\sum _{n=1}^N\frac{e^{\beta E_n}}{Z_{+}}|n\rangle \langle n|,$$

with $Z_{+}=Tr\left[e^{+\beta H_S}\right]$ between a system in the state ${\omega }_{\beta }(-H_S)$ with copies of itself in the Gibbs state ${\omega }_{\beta }\left(H_S\right)$.

In this subsection, we consider the battery S and auxiliary systems A identical qubits; i.e., the battery Hamiltonian is $H_S=(h/2){\sigma }_S^z$, and the baths Hamiltonians are $H_A=(h/2){\sigma }_A^z$, with $h>0$. Hereafter, ${\sigma }^x,{\sigma }^y$ and ${\sigma }^z$ are Pauli matrices.

The coupling between the system and the bath qubit is

$$V=a({\sigma }_S^{+}{\sigma }_A^{+}+{\sigma }_S^{-}{\sigma }_A^{-}),$$

with ${\sigma }^{\pm }=({\sigma }^x\pm {\sigma }^y)/2$, and is such that $[{\sigma }_A^z-{\sigma }_S^z,V]=0$, i.e., $H_0=-H_S$.

In the basis defined by ${\sigma }^z|\uparrow \rangle =|\uparrow \rangle$ and ${\sigma }^z|\downarrow \rangle =-|\downarrow \rangle$, the eigenvalues and eigenvectors of $H_S$ and $H_0$ are

$$\begin{array}{cccccc}\hfill E_2& =h/2,\hfill & \hfill E_2^0& =-h/2,\hfill & \hfill |2\rangle & =|\uparrow \rangle \hfill \end{array}$$

(30)

$$\begin{array}{cccccc}\hfill E_1& =-h/2,\hfill & \hfill E_1^0& =h/2,\hfill & \hfill |1\rangle & =|\downarrow \rangle \hfill \end{array}$$

(31)

and the ordering permutation is $({\pi }_1,{\pi }_2)=(2,1)$. Thus, on the above basis, the equilibrium state is

$${\omega }_{\beta }\left(H_0\right)={\omega }_{\beta }(-H_S)=\frac{e^{\beta \frac{h}{2}}}{Z}|2\rangle \langle 2|+\frac{e^{-\beta \frac{h}{2}}}{Z}|1\rangle \langle 1|,$$

and the passive state for the system is

$${\sigma }_{{\omega }_{\beta }\left(H_0\right)}={\omega }_{\beta }\left(H_S\right)=\frac{e^{-\beta \frac{h}{2}}}{Z}|2\rangle \langle 2|+\frac{e^{\beta \frac{h}{2}}}{Z}|1\rangle \langle 1|,$$

where $Z=Z_{+}=2cosh(\beta h/2).$ With Equations (30) and (31), and the permutation $\pi$, we can evaluate the transition probabilities in Equation (29). The ergotropy of the battery in the equilibrium state ${\omega }_{\beta }(-H_S)$ is $\mathcal{W}=htanh\beta h/2$. From Equations (23) and (24), we see that the thermodynamic efficiency of the process is ${\eta }_{\mathrm{th}}=1/2,$ independent of the inverse temperature $\beta$.

The recharging process in this single-qubit battery (1Q) is determined by the stochastic matrix (see Equation (21))

$$T_{1Q}=\left(\begin{array}{cc}1-\frac{e^{\beta \frac{h}{2}}}{Z}g(a,h)& \frac{e^{-\beta \frac{h}{2}}}{Z}g(a,h)\\ \frac{e^{\beta \frac{h}{2}}}{Z}g(a,h)& 1-\frac{e^{-\beta \frac{h}{2}}}{Z}g(a,h)\end{array}\right)$$

(32)

where $g(a,h)=\frac{a^2{sin}^2(\tau \sqrt{h^2+a^2}/\hslash )}{h^2+a^2}$ and $Z=e^{\beta \frac{h}{2}}+e^{-\beta \frac{h}{2}}$. It is a regular stochastic matrix if $g(a,h)\ne 0$.

5.1.1. Fluctuating Efficiency

The fluctuating efficiency (see Equation (29)) takes the values

$${\eta }_{11}=-{\eta }_{22}=\infty ,{\eta }_{12}={\eta }_{21}=\frac{1}{2}$$

Its distribution Equation (28) is

$$p_{\eta }\left(x\right)=\delta (x-\infty ){\mathrm{P}}_{\infty }+\delta (x+\infty ){\mathrm{P}}_{-\infty }+\delta \left(x-\frac{1}{2}\right){\mathrm{P}}_{\frac{1}{2}}$$

with

$${\mathrm{P}}_{\infty }=P_{1\to 1}={\mathrm{P}}_{-\infty }=P_{2\to 2}=\frac{1}{Z^2},{\mathrm{P}}_{\frac{1}{2}}=P_{1\to 2}+P_{2\to 1}=\frac{e^{\beta h}+e^{-\beta h}}{Z^2}$$

(33)

The explicit formulas at the right follow from Equation (29), which is valid if $g(a,h)\ne 0$ in $T_{1Q}$.

In Figure 1a, we depict the probabilities ${\mathrm{P}}_{\eta }$ as functions of $\beta h$ and see that for $\beta h\gg 1$ with probability 1; the fluctuating efficiency equals the thermodynamic efficiency $1/2$, because, as we see in Figure 1b, $P_{1\to 2}\to 1$, reflecting the charging character of the process.

The diagrams in Figure 2 depict $\kappa$ transitions (left, up to down), followed by $\gamma$ transitions (right, down to up). The values of all variables and their probability are given underneath.

The numbers correspond to the energy levels 1 and 2. In the limit of large temperature, $\beta \to 0$, all these processes have the same probability $1/4$, while at low temperature $\beta \to \infty$, the probability of the second process goes to one and the others to zero. Only the third diagram has a transition assisted by heat, $q_{\gamma }=h$. We extract energy $w_{\gamma }=-2h<0$ in the $\gamma$ process and invest ${\varpi }_{\kappa }=-h$ in the $\kappa$ process. This cycle is the least likely. Its probability is $e^{-\beta h}/Z^2$ and decreases quickly as $\beta h$ increases.

5.1.2. Equilibrium Fluctuation

Let us analyze the fluctuations when maintaining the charged state, i.e., those of the process ${\omega }_{\beta }\left(H_0\right)\to {\mathcal{E}}^L\left({\omega }_{\beta }\left(H_0\right)\right)={\omega }_{\beta }\left(H_0\right)$; see Section 4.2. As we can verify in the examples above, and as shown in [72], the transition matrices T for maps with equilibrium satisfy the detailed balance condition $T_{m|n}{e}^{-\beta E_n^0}=T_{n|m}{e}^{-\beta E_m^0}$. From this fact, it is simple to show that $P_{n\to m}^{\left(L\right)}=P_{m\to n}^{\left(L\right)}$ with $p_{\mathrm{ini}}\left(n\right)=e^{-\beta E_n^0}/Z_0$ in Equation (21).

We are interested in distinguishing fluctuations in an active equilibrium state from fluctuations in a Gibbs equilibrium state. The main difference is that the probability distribution of equilibrium work fluctuation is $p_w\left(x\right)\ne \delta \left(x\right)$ for the former, reflecting an active agent, and $p_w\left(x\right)=\delta \left(x\right)$ for the latter, reflecting a passive agent.

To investigate other differences, we consider our charging map $\mathcal{E}$ and a thermal map ${\mathcal{E}}^{\mathrm{Thm}}$ for a qubit. The map ${\mathcal{E}}^{\mathrm{Thm}}$ is obtained by coupling the qubit to an auxiliary thermal qubit with $V=a({\sigma }_S^{+}{\sigma }_A^{-}+{\sigma }_S^{-}{\sigma }_A^{+})$ and tracing out the auxiliary system. The resulting map is thermal (i.e., a map with the Gibbs equilibrium state), and the transition matrix for this process is

$$T^{\mathrm{Thm}}=\left(\begin{array}{cc}1-\frac{e^{-\beta \frac{h}{2}}}{Z}g(a,0)& \frac{e^{\beta \frac{h}{2}}}{Z}g(a,0)\\ \frac{e^{-\beta \frac{h}{2}}}{Z}g(a,0)& 1-\frac{e^{\beta \frac{h}{2}}}{Z}g(a,0)\end{array}\right)$$

where $g(a,0)={sin}^2(\tau a/\hslash )$ and $Z=e^{\beta \frac{h}{2}}+e^{-\beta \frac{h}{2}}.$ $T^{\mathrm{Thm}}$ is a regular stochastic matrix if $g(a,0)\ne 0$. The most crucial difference between $T^{\mathrm{Thm}}$ and $T_{1Q}$ in Equation (32) is the position of the factors $e^{\pm \beta h/2}$.

For the charging map, one can show $P_{2\to 2}^{\left(L\right)}>P_{1\to 1}^{\left(L\right)}$, reflecting the higher population of the excited state in the active equilibrium. Instead, for the thermal map, $P_{1\to 1}^{\left(L\right)\mathrm{Thm}}>P_{2\to 2}^{\left(L\right)\mathrm{Thm}}$, reflecting the higher population of the ground state in Gibbs equilibrium. On the other hand, energy fluctuations due to $1\leftrightarrow 2$ transitions are qualitatively similar if $g(a,h)\approx g(a,0)$ for processes with finite L but are indistinguishable for $L\to \infty$. Indeed, for $L\to \infty$, we have

$$P_{1\to 2}^{(\infty )\mathrm{Thm}}=P_{2\to 1}^{(\infty )\mathrm{Thm}}=\frac{1}{Z^2},P_{1\to 1}^{(\infty )\mathrm{Thm}}=\frac{e^{\beta h}}{Z^2},P_{2\to 2}^{(\infty )\mathrm{Thm}}=\frac{e^{-\beta h}}{Z^2}$$

and for the charging map,

$$P_{2\to 1}^{(\infty )}=P_{1\to 2}^{(\infty )}=\frac{1}{Z^2},P_{2\to 2}^{(\infty )}=\frac{e^{\beta h}}{Z^2},P_{1\to 1}^{(\infty )}=\frac{e^{-\beta h}}{Z^2}.$$

Thus, these processes are very similar at the level of energy fluctuations.

5.2. Two-Qubit Battery

We consider a two-qubit battery with Hamiltonian [26]

$$H_S=\frac{h}{2}\left({\sigma }_1^z+{\sigma }_2^z\right)+J\left({\sigma }_1^x{\sigma }_2^x+{\sigma }_1^y{\sigma }_2^y\right),$$

coupled with

$$V=J^{\prime }({\sigma }_A^x{\sigma }_1^x+{\sigma }_A^y{\sigma }_1^y),$$

to auxiliary systems with Hamiltonian $H_A=\frac{h}{2}{\sigma }_A^z$ in the thermal state. The corresponding map $\mathcal{E}$ has the equilibrium state ${\omega }_{\beta }\left(H_0\right)$ with $H_0=\frac{h}{2}\left({\sigma }_1^z+{\sigma }_2^z\right)$.

The eigenvalues and eigenvectors of $H_S$ and $H_0$ in the basis defined by ${\sigma }^z|\uparrow \rangle =|\uparrow \rangle$ and ${\sigma }^z|\downarrow \rangle =-|\downarrow \rangle$ are

$$\begin{array}{cccccc}\hfill E_3& =h,\hfill & \hfill E_3^0& =h,\hfill & \hfill |3\rangle & =|\uparrow \uparrow \rangle ,\hfill \end{array}$$

(34)

$$\begin{array}{cccccc}\hfill E_4& =2J,\hfill & \hfill E_4^0& =0,\hfill & \hfill |4\rangle & =(|\uparrow \downarrow \rangle +|\downarrow \uparrow \rangle )/\sqrt{2},\hfill \end{array}$$

(35)

$$\begin{array}{cccccc}\hfill E_1& =-2J,\hfill & \hfill E_1^0& =0,\hfill & \hfill |1\rangle & =(|\uparrow \downarrow \rangle -|\downarrow \uparrow \rangle )/\sqrt{2},\hfill \end{array}$$

(36)

$$\begin{array}{cccccc}\hfill E_2& =-h,\hfill & \hfill E_2^0& =-h,\hfill & \hfill |2\rangle & =|\downarrow \downarrow \rangle .\hfill \end{array}$$

(37)

We take $2J>h>0$ such that $E_{i+1}>E_i$. The permutation that orders $E_{{\pi }_{i+1}}^0\ge E_{{\pi }_i}^0$ is $({\pi }_1,{\pi }_2,{\pi }_3,{\pi }_4)=(2,1,4,3)$. Thus, on the above basis, the equilibrium state is

$${\omega }_{\beta }\left(H_0\right)=\frac{e^{-\beta h}}{Z_0}|3\rangle \langle 3|+\frac{1}{Z_0}(|1\rangle \langle 1|+|4\rangle \langle 4|)+\frac{e^{\beta h}}{Z_0}|2\rangle \langle 2|,$$

and the passive state for the system is

$${\sigma }_{{\omega }_{\beta }\left(H_0\right)}=\frac{e^{\beta h}}{Z_0}|1\rangle \langle 1|+\frac{1}{Z_0}(|2\rangle \langle 2|+|3\rangle \langle 3|)+\frac{e^{-\beta h}}{Z_0}|4\rangle \langle 4|,$$

where $Z_0=2+2cosh\left(\beta h\right).$ The ergotropy of the equilibrium state $\mathcal{W}=Tr\left[H_S({\omega }_{\beta }\left(H_0\right)-{\sigma }_{{\omega }_{\beta }\left(H_0\right)})\right]$ is

$$\mathcal{W}=(2J-h)\frac{sinh\beta h}{1+cosh\beta h}.$$

The work performed in the charging process ${\sigma }_{{\omega }_{\beta }\left(H_0\right)}\to {\omega }_{\beta }\left(H_0\right)$ is

$$W_R=2J\frac{sinh\beta h}{1+cosh\beta h}.$$

We see that the thermodynamic efficiency is ${\eta }_{\mathrm{th}}=\mathcal{W}/W_R=1-\frac{h}{2J}$ independently of the inverse temperature $\beta$.

The recharging process in this two-qubit battery (2Q) is determined by the stochastic matrix (see Equation (21))

$$T_{2Q}=\frac{1}{{(J^2+J^{\prime 2})}^2}\left(\begin{array}{cccc}{\Phi }^2& \frac{2}{(1+e^{\beta h})}\Phi \Psi & \frac{2e^{\beta h}}{(1+e^{\beta h})}\Phi \Psi & {\Psi }^2\\ \frac{2e^{\beta h}}{(1+e^{\beta h})}\Phi \Psi & \frac{e^{\beta h}{(J^2+J^{\prime 2})}^2+\Delta }{(1+e^{\beta h})}& 0& \frac{2e^{\beta h}}{(1+e^{\beta h})}\Phi \Psi \\ \frac{2}{(1+e^{\beta h})}\Phi \Psi & 0& \frac{{(J^2+J^{\prime 2})}^2+e^{\beta h}\Delta }{(1+e^{\beta h})}& \frac{2}{(1+e^{\beta h})}\Phi \Psi \\ {\Psi }^2& \frac{2}{(1+e^{\beta h})}\Phi \Psi & \frac{2e^{\beta h}}{(1+e^{\beta h})}\Phi \Psi & {\Phi }^2\end{array}\right),$$

(38)

with

$$\Phi =J^2+J^{\prime 2}{cos}^2\left(\frac{\tau }{\hslash }\sqrt{J^2+J^{\prime 2}}\right),\Psi =J^{\prime 2}{sin}^2\left(\frac{\tau }{\hslash }\sqrt{J^2+J^{\prime 2}}\right),\Delta ={(\Phi -\Psi )}^2,$$

which is a regular stochastic matrix excepts at points with $\Psi =0$ or $\Phi =0$, as one can check by computing $T^2$.

5.2.1. Fluctuating Efficiency

For the fluctuating efficiency Equation (29), we have

$$\begin{array}{cc}\hfill {\eta }_{12}& ={\eta }_{13}={\eta }_{21}={\eta }_{34}={\eta }_{42}={\eta }_{43}=1-\frac{h}{2J}\hfill \end{array}$$

(39)

$$\begin{array}{cc}\hfill {\eta }_{14}& ={\eta }_{41}=\frac{1}{2}(1-\frac{h}{2J})\hfill \end{array}$$

(40)

$$\begin{array}{cc}\hfill {\eta }_{32}& =-{\eta }_{23}=\infty \hfill \end{array}$$

(41)

$$\begin{array}{cc}\hfill {\eta }_{24}& ={\eta }_{31}=-(1-\frac{h}{2J})\hfill \end{array}$$

(42)

and ${\eta }_{11}={\eta }_{33}=-{\eta }_{22}=-{\eta }_{44}=\infty$. Its distribution follows from Equation (28), and it is

$$\begin{array}{c}\hfill p_{\eta }\left(x\right)=\delta (x-\infty ){\mathrm{P}}_{\infty }+\delta (x+\infty ){\mathrm{P}}_{-\infty }+\delta \left(x-1+\frac{h}{2J}\right){\mathrm{P}}_{(1-\frac{h}{2J})}+\\ \hfill \delta \left(x+1-\frac{h}{2J}\right){\mathrm{P}}_{-(1-\frac{h}{2J})}+\delta \left(x-\frac{1}{2}+\frac{h}{4J}\right){\mathrm{P}}_{(1/2)(1-\frac{h}{2J})}\end{array}$$

with

$$\begin{array}{cc}\hfill & {\mathrm{P}}_{\infty }=P_{3\to 2}+P_{1\to 1}+P_{3\to 3}=\frac{2e^{\beta h}+e^{-\beta h}}{Z_0^2},\hfill \end{array}$$

(43)

$$\begin{array}{cc}\hfill & {\mathrm{P}}_{-\infty }=P_{2\to 3}+P_{2\to 2}+P_{4\to 4}=\frac{e^{\beta h}+2e^{-\beta h}}{Z_0^2},\hfill \end{array}$$

(44)

$$\begin{array}{cc}\hfill & {\mathrm{P}}_{(1-\frac{h}{2J})}=P_{1\to 2}+P_{1\to 3}+P_{2\to 1}+P_{3\to 4}+P_{4\to 2}+P_{4\to 3}=\frac{2+{(e^{\beta h}+e^{-\beta h})}^2}{Z_0^2},\hfill \end{array}$$

(45)

$$\begin{array}{cc}\hfill & {\mathrm{P}}_{-(1-\frac{h}{2J})}=P_{3\to 1}+P_{2\to 4}=\frac{2}{Z_0^2},\hfill \end{array}$$

(46)

$$\begin{array}{cc}\hfill & {\mathrm{P}}_{(1/2)(1-\frac{h}{2J})}=P_{1\to 4}+P_{4\to 1}=\frac{(e^{-\beta h}+e^{\beta h})}{Z_0^2}.\hfill \end{array}$$

(47)

The explicit formulas on the right follow from Equation (29) and are valid for parameters $\tau ,J$ and $J^{\prime }$ in which $T_{2Q}$ is regular.

In Figure 3a, we plot the probabilities ${\mathrm{P}}_{\eta }$ in Equations (44)–(47) as a function of $\beta h$. We see that for small $\beta h$, the average efficiency does not exist. On the other hand, when $\beta h\gg 1$, the efficiency goes to the thermodynamic efficiency with a probability of one because the work becomes deterministic.

The diagrams in Figure 4, summarize all possible extracting–recharging cycles. The numbers correspond to the energy levels $1,2,3$, and 4. Since extracting the ergotropy only allows transitions $\kappa :m^{\prime }\to n$ with $m^{\prime }={\pi }_n$, we have four possible processes $\kappa$. From Figure 3b, we see that the only process with a high probability for large $\beta h$ is the sequence $2\stackrel{\kappa }{\to }1\stackrel{\gamma }{\to }2$ contained in the first diagram. Its efficiency equals thermodynamic efficiency. Green arrows are processes assisted by heat ($q_{\gamma }>0$). These have very low probabilities, as depicted in Figure 3b. In Figure 3b, we see that $P_{1\to 2}$ goes to one in that limit. Second in importance are $P_{1\to 4}$, associated with the largest charge, but in the extracting $\kappa$ process, one has $4\stackrel{\kappa }{\to }3$, and the $\gamma$ process starting in 3 reaches $1,2$, or 4 with similar probabilities and the battery is noisy.

5.2.2. Heat and Work Fluctuations in the Partial Recharging Process

Here, we consider the process ${\mathcal{E}}^L$ starting in the state ${\sigma }_{{\omega }_{\beta }\left(H_0\right)}$ and evaluate the heat and work distributions. Hence, we consider Equations (19) and (20) with $P_{n\to m}^{\left(L\right)}={\left(T^L\right)}_{m|n}\frac{e^{-\beta E_{{\pi }_n}^0}}{Z0}$, with the permutation $\pi$ ordering the eigenvalues of $H_0$ by increasing values.

For the two-qubit battery, we obtain

$$\begin{array}{cc}\hfill p_w^{\left(L\right)}\left(x\right)& =\delta \left(x\right)A_0^{\left(L\right)}+\delta (x-4J)A_1^{\left(L\right)}+\delta (x-2J)A_3^{\left(L\right)}+\delta (x+4J)A_2^{\left(L\right)}+\delta (x+2J)A_4^{\left(L\right)},\hfill \end{array}$$

(48)

$$\begin{array}{cc}\hfill p_q^{\left(L\right)}\left(x\right)& =\delta \left(x\right)B_0^{\left(L\right)}+\delta (x-h)B_4^{\left(L\right)}+\delta (x-2h)B_2^{\left(L\right)}+\delta (x+h)B_3^{\left(L\right)}+\delta (x+2h)B_1^{\left(L\right)}\hfill \end{array}$$

(49)

with

$$\begin{array}{cccc}\hfill A_0^{\left(L\right)}& =P_{2\to 3}^{\left(L\right)}+P_{3\to 2}^{\left(L\right)},\hfill & \hfill B_0^{\left(L\right)}& =P_{1\to 4}^{\left(L\right)}+P_{4\to 1}^{\left(L\right)},\hfill \end{array}$$

(50)

$$\begin{array}{cccc}\hfill A_1^{\left(L\right)}& =P_{1\to 4}^{\left(L\right)},\hfill & \hfill B_1^{\left(L\right)}& =P_{3\to 2}^{\left(L\right)},\hfill \end{array}$$

(51)

$$\begin{array}{cccc}\hfill A_2^{\left(L\right)}& =P_{4\to 1}^{\left(L\right)},\hfill & \hfill B_2^{\left(L\right)}& =P_{2\to 3}^{\left(L\right)},\hfill \end{array}$$

(52)

$$\begin{array}{cccc}\hfill A_3^{\left(L\right)}& =P_{1\to 2}^{\left(L\right)}+P_{1\to 3}^{\left(L\right)}+P_{2\to 4}^{\left(L\right)}+P_{3\to 4}^{\left(L\right)},\hfill & \hfill B_3^{\left(L\right)}& =P_{1\to 2}^{\left(L\right)}+P_{3\to 1}^{\left(L\right)}+P_{4\to 2}^{\left(L\right)}+P_{3\to 4}^{\left(L\right)},\hfill \end{array}$$

(53)

$$\begin{array}{cccc}\hfill A_4^{\left(L\right)}& =P_{2\to 1}^{\left(L\right)}+P_{3\to 1}^{\left(L\right)}+P_{4\to 2}^{\left(L\right)}+P_{4\to 3}^{\left(L\right)},\hfill & \hfill B_4^{\left(L\right)}& =P_{2\to 1}^{\left(L\right)}+P_{1\to 3}^{\left(L\right)}+P_{2\to 4}^{\left(L\right)}+P_{4\to 3}^{\left(L\right)},\hfill \end{array}$$

(54)

where $A_i^{\left(L\right)}\ne B_i^{\left(L\right)}$ for finite L but $A_i^{(\infty )}=B_i^{(\infty )}$ with

$$A_0^{(\infty )}=\frac{6cosh\beta h}{Z_0^2},A_1^{(\infty )}=\frac{e^{\beta h}}{Z_0^2},A_2^{(\infty )}=\frac{e^{-\beta h}}{Z_0^2},A_3^{(\infty )}=\frac{e^{2\beta h}+3}{Z_0^2},A_4^{(\infty )}=\frac{3+e^{-2\beta h}}{Z_0^2}.$$

This means that the average work $W^{\left(L\right)}$ and average heat $Q^{\left(L\right)}$

$$W^{\left(L\right)}=2J(A_3^{\left(L\right)}-A_4^{\left(L\right)})+4J(A_1^{\left(L\right)}-A_2^{\left(L\right)})\underset{L\to \infty }{\to }\frac{2Jsinh\beta h}{1+cosh\beta h}$$

$$Q^{\left(L\right)}=h(B_4^{\left(L\right)}-B_3^{\left(L\right)})+2h(B_2^{\left(L\right)}-B_1^{\left(L\right)})\underset{L\to \infty }{\to }\frac{-hsinh\beta h}{1+cosh\beta h}$$

become proportional when $L\to \infty$.

Since Markov chains converge exponentially quickly to the stationary state, it is unnecessary to consider a large L to observe the asymptotic distribution. However, since the convergence rate depends on the map’s parameters, we see deviations from it near the points where $\Phi =0$ or $\Psi =0$ in Equation (38). To illustrate this point, we plot in Figure 5 the probabilities $A_0^{\left(L\right)},B_0^{\left(L\right)},A_2^{\left(L\right)}$, and $B_2^{\left(L\right)}$ for various values of L and varying map parameters.

Figure 5 shows that for $L=20$, convergence is achieved with a duration of each iteration $\tau /\hslash =1$. Note that the dependence on $\tau$ in $T_{2Q}$ is periodic; see Equation (38). In the limit $L\to \infty ,\tau \to 0$ and $J^{\prime }=j/\sqrt{\tau }$, the dynamics of the battery has the Lindblad form [69] and converges exponentially quickly to the equilibrium distribution.

6. Discussion

We have studied stochastic fluctuations in repeated interaction processes subjected to the two-point energy-measurement scheme. Because map $\mathcal{E}$ has an equilibrium state, all quantities are expressed in terms of system properties simplifying their study because one does not require measuring the environment. We have shown that the equilibrium distribution of the map dominates the distributions, except at particular points in the parameter space of the map, where its details become essential. Near these zones, the convergence rate towards the asymptotic value is low, requiring larger values of L to reach it. The quantum aspect of the system is relevant near these zones since the Planck constant appears in the parameters that set the convergence rate to the stationary state. We have applied these results to study active equilibrium fluctuations, fluctuations in the charging process of a quantum battery, and efficiency fluctuations of the cycle charging and extracting energy for the battery in two examples. The fluctuating efficiency converges to the thermodynamic efficiency of these examples in the low-temperature limit, where the batteries operate in the cycle $2\stackrel{\kappa }{\to }1\stackrel{\gamma }{\to }2$ and are reliable. On the other hand, at large temperatures, where heat assists some transitions, all cycles are probable, and the battery is unreliable.

For future research, it would be interesting to extend the results obtained here for single-cycle efficiency to the case of an arbitrary number of cycles. As this number increases, universal statistical behaviors have been shown to appear in other machines [58,59,68]. Likewise, considering the collective boost in power for dissipative quantum batteries [79] and the result in [43], studying fluctuations as the number of batteries increases is of similar interest.