Quantum Signatures in the Quantum Carnot Cycle

TL;DR

This paper explores quantum Carnot-like cycles that operate irreversibly with non-zero power, utilizing shortcuts to equilibrium and coherence effects, revealing quantum signatures and trade-offs between power and efficiency.

Contribution

It introduces and compares three quantum Carnot-analog cycles employing shortcuts and coherence, highlighting quantum signatures and performance trade-offs.

Findings

The Carnot-shortcut cycle shows a maximum power at an optimal cycle time.

The Endo-global cycle exhibits positive power with quantum coherence effects.

Quantum signatures are observed at short cycle times, indicating non-classical behavior.

Abstract

The Carnot cycle combines reversible isothermal and adiabatic strokes to obtain optimal efficiency, at the expense of a vanishing power output. Quantum Carnot-analog cycles are constructed and solved, operating irreversibly at non-vanishing power. Swift thermalization is obtained in the isotherms utilizing a shortcut to equilibrium protocols and the adiabats employ frictionless unitary shortcuts. The working medium in this study is composed of a particle in a driven harmonic trap. For this system, we solve the dynamics employing a generalized canonical state. Such a description incorporates both changes in energy and coherence. This allows comparing three types of Carnot-analog cycles, Carnot-shortcut, Endo-shortcut and Endo-global. The Carnot-shortcut engine demonstrates the trade-off between power and efficiency. It posses a maximum in power, a minimum cycle time where it becomes a dissipator and for a diverging cycle time approaches the ideal Carnot efficiency. The irreversibility of the cycle arises from non-adiabatic driving, which generates coherence. To study the role of coherence we compare the performance of the shortcut cycles, where coherence is limited to the interior of the strokes, with the Endo-global cycle where the coherence never vanishes. The Endo-global engine exhibits a quantum signature at a short cycle-time, manifested by a positive power output while the shortcut cycles become dissipators. If energy is monitored the back action of the measurement causes dephasing and the power terminates.

Figures

Tables

Table 1. Table 1: Cycle parameters

Cycle type	Carnot-shortcut	Endo-Shortcut	Endo-global
Baths temperatures	$T_c=5$ $T_h=8$	$T_c^{\prime }=5.25$ $T_h^{\prime }=7.75$	$T_c=5$ $T_h=8$
Cycle frequencies	${\omega }_1=10$ ${\omega }_2={\omega }_3\frac{T_h}{T_c}=6.25$ ${\omega }_3=5$ ${\omega }_4={\omega }_1\frac{T_c}{T_h}=7.5$	${\omega }_1=10$ ${\omega }_2=6.25$ ${\omega }_3=5$ ${\omega }_4=7.5$	${\omega }_1^g=9.6875$ ${\omega }_2^g=7.75$ ${\omega }_3^g=5.25$ ${\omega }_4^g=6.5625$
Operation type	Branch coherence	Branch coherence	Global coherence

Table 2. Table 2: Cycle parameters

Cycle parameter	Value (atomic units)
Coupling strength, $\frac{{\left|\overrightarrow{d}\right|}^2}{4\pi {ϵ}_0\hslash c}$	0.05
Stroke duration of the adiabatic expansion of the shortcut cycles	5
Stroke duration of the adiabatic compression of the shortcut cycles	5

Equations

$$\hat{H}\left(t\right)=\frac{{\hat{P}}^2}{2m}+\frac{1}{2}m{\omega }^2\left(t\right){\hat{Q}}^2,$$

$$\frac{d\hat{\rho }\left(t\right)}{dt}=\mathcal{L}\left(t\right)\hat{\rho }\left(t\right),$$

$${\mathcal{U}}_{cyc}={\mathcal{U}}_{hc}{\mathcal{U}}_c{\mathcal{U}}_{ch}{\mathcal{U}}_h,$$

$$Coh=\frac{\sqrt{⟨\hat{L}\left(t\right){⟩}^2+⟨\hat{C}\left(t\right){⟩}^2}}{\hslash \omega \left(t\right)}.$$

$$W={\int }_0^t⟨\frac{\partial \hat{H}\left(t^{\prime }\right)}{\partial t^{\prime }}⟩d{t}^{\prime },$$

$$\frac{{\omega }_3}{{\omega }_2}=\frac{{\omega }_4}{{\omega }_1}=\frac{T_c}{T_h}.$$

$$\mathcal{C}>\frac{T_h}{T_c}.$$

$$\displaystyle\begin{array}[]{lll}W_{C}=&\hbar\Delta\omega_{32}\left({n_{2}+1}\right)+\hbar\Delta\omega_{14}\left({n_{1}+1}\right)+k_{B}\left({T_{h}-T_{c}}\right)\rm{ln}\left({\frac{n_{1}}{n_{2}}}\right)~{}~{},\end{array}$$

$$W_C=-k_B{T}_h{\eta }_C\ln \left(\mathcal{C}\right).$$

$$\mathcal{P}-\frac{|W_C|}{{\tau }_{cycle}}\propto -\frac{F}{{\tau }_{cycle}^2},$$

$$\displaystyle\begin{array}[]{ll}\frac{d}{dt}\hat{\rho}\left({t}\right)=&-\frac{i}{\hbar}\left[{\hat{H}\left({t}\right),\hat{\rho}\left({t}\right)}\right]-\gamma_{d}\left[{\hat{H}\left({t}\right),\left[{\hat{H}\left({t}\right),\hat{\rho}\left({t}\right)}\right]}\right]~{}~{},\end{array}$$

$$\displaystyle\begin{array}[]{ll}\frac{d}{dt}\tilde{\rho}_{S}\left({t}\right)={\cal L}\tilde{\rho}=&k_{\downarrow}\left({t}\right)\left({\hat{b}\tilde{\rho}_{S}\left({t}\right)\hat{b}^{\dagger}-\frac{1}{2}\{\hat{b}^{\dagger}\hat{b},\tilde{\rho}_{S}\left({t}\right)\}}\right)\\ &+k_{\uparrow}\left({t}\right)\left({\hat{b}^{\dagger}\tilde{\rho}_{S}\left({t}\right)\hat{b}-\frac{1}{2}\{\hat{b}\hat{b}^{\dagger},\tilde{\rho}_{S}\left({t}\right)\}}\right)~{}~{}~{}.\end{array}$$

$$k_{\downarrow }\left(t\right)=k_{\uparrow }\left(t\right)e^{\alpha \left(t\right)/k_{B}T}=\frac{\alpha \left(t\right)|\vec{d}{|}^2}{4\pi {\epsilon }_0\kappa \hslash c}\left(1+N\left(\alpha \left(t\right)\right)\right),$$

$$\alpha \left(t\right)=\sqrt{1-{\left(\dot{\omega }\left(t\right)/\left(2{\omega }^2\left(t\right)\right)\right)}^2}\omega \left(t\right).$$

$${\tilde{\rho }}_S\left(t\right)=Z^{-1}{e}^{\gamma \left(t\right){\hat{b}}^2}{e}^{\beta \left(t\right){\hat{b}}^{†}\hat{b}}{e}^{{\gamma }^{\ast }\left(t\right){\hat{b}}^{†2}},$$

$${\tilde{\rho }}_S\left(\beta \left(t\right),\mu \left(t\right)\right)=Z^{-1}{e}^{\beta {\hat{b}}^{†}\hat{b}\left(\mu \right)},$$

$$\dot{\beta }=k_{\downarrow }\left(t\right)\left(e^{\beta }-1\right)+k_{\uparrow }\left(t\right)\left(e^{-\beta }-1\right).$$

$$y\left(s\right)=y\left(0\right)+c_3{\left(t/t_f\right)}^3+c_4{\left(t/t_f\right)}^4+c_5{\left(t/t_f\right)}^5,$$

$$\hat{\rho }\left(t\right)=Z^{-1}{e}^{-\frac{\hslash \omega \left(0\right)}{k_B{T}_i}},$$

$$\hat{I}\left(t\right)=\frac{1}{2}\left(\alpha \left(t\right){\hat{Q}}^2+\beta \left(t\right){\hat{P}}^2+\gamma \left(\hat{Q}\hat{P}+\hat{P}\hat{Q}\right)\right),$$

$$\frac{d\hat{I}}{dt}=\frac{\partial I}{\partial t}+\frac{1}{i\hslash }\left[\hat{I},\hat{H}\right]=0.$$

$$\left[\hat{H}\left(0\right),\hat{I}\left(0\right)\right]=\left[\hat{H}\left(t_f\right),\hat{I}\left(t_f\right)\right]=0.$$

$$\displaystyle\begin{array}[]{cc}\dot{\alpha}=2m\omega^{2}\gamma\\ \dot{\beta}=-\frac{2}{m}\gamma\\ \dot{\gamma}=-\frac{1}{m}\alpha+m\omega^{2}\beta\end{array}~{}~{}.$$

$$\frac{d}{dt}\left(m^2{\omega }^2\sigma +m^2\ddot{\sigma }\right)\sigma +3\dot{\sigma }\left(m^2{\omega }^2\sigma +m^2\ddot{\sigma }\right)=0.$$

$$\frac{c}{{\sigma }^2}=m^2{\omega }^2{\sigma }^2+m^2\sigma \ddot{\sigma },$$

$$\hat{I}\left(t\right)=\frac{1}{2}\left(\frac{1}{{\rho }^2}{\hat{Q}}^2+\left(m\dot{\rho }\hat{Q}-\rho \hat{P}\right)\right),$$

$$m^2{\omega }^2\left(t\right)\rho +m^2\ddot{\rho }-\frac{1}{{\rho }^3}=0.$$

$$\displaystyle\begin{array}[]{ll}\dot{\rho}\left({0}\right)=\ddot{\rho}\left({0}\right)=\dot{\rho}\left({t_{f}}\right)=\ddot{\rho}\left({t_{f}}\right)=0\\ \rho\left({0}\right)=\frac{1}{\sqrt{m\omega\left({0}\right)}}\\ \rho\left({t_{f}}\right)=\frac{1}{\sqrt{m\omega\left({0}\right)}}\sqrt{\frac{\omega_{0}}{\omega_{f}}}~{}~{}.\end{array}$$

$$\hat{I}\left(t\right)|\lambda \left(t\right)⟩=\lambda |\lambda \left(t\right)⟩,$$

$$\displaystyle\begin{array}[]{ll}\left<\lambda\right|\hat{H}\left|\lambda\right>=\frac{\hbar}{2}\left({m\dot{\rho}^{2}+\frac{1}{m\rho^{2}}+m\omega^{2}\left({t}\right)\rho^{2}}\right)\left({\lambda+\frac{1}{2}}\right)\\ \left<\lambda\right|\hat{L}\left|\lambda\right>=\frac{\hbar}{2}\left({m\dot{\rho}^{2}+\frac{1}{m\rho^{2}}-m\omega^{2}\left({t}\right)\rho^{2}}\right)\left({\lambda+\frac{1}{2}}\right)\\ \left<\lambda\right|\hat{C}\left|\lambda\right>=\omega\left({t}\right)\hbar m\dot{\rho}\rho\left({\lambda+\frac{1}{2}}\right)~{}~{}.\end{array}$$

Full text

• August 2019

1 Introduction

In 1824 Sadi Carnot envisioned a reversible reciprocating heat engine, composed of two adiabatic and two isothermal strokes. He argued that such a reversible engine is universal, and produces optimal work, depending only on the hot and cold bath temperatures, $T_{h}$ and $T_{c}$ [1]. This cycle is termed Carnot cycle and serves as a template for all reversible engines. Generally, a heat engine transforms heat extracted from a hot bath, $Q_{h}$, to work $W$. According to the second law of thermodynamics this work is necessarily accompanied by a heat flow to the cold bath $Q_{c}=|W|-Q_{h}$. The efficiency of the energy transformation is defined by $\eta=-W/Q_{h}$, which is bounded by the Carnot cycle efficiency $\eta<\eta_{C}=1-T_{c}/T_{h}$, for any heat cycle. This universal bound depends only on the bath temperatures, irrespective of the specific working medium. It translates a fundamental limitation, the second law of thermodynamics, to a practical operational limit applying to all engines. Optimal efficiency is achieved only for a reversible operation.

Practically, any finite power engine operates under irreversible conditions, which leads to a trade-off between efficiency and power. This relation has been extensively studied in the framework of finite-time thermodynamics [2, 3]. A prominent emerging result states that under endoreversible conditions the efficiency at maximum power is given by the Curzon-Ahlborn-Novikov efficiency $\eta_{CA}=1-\sqrt{T_{c}/T_{h}}$ [4, 5], which has been generalized for weak dissipation [6, 7].

In the quantum regime energy transfer is constrained as well by the second law of thermodynamics [8, 9, 10, 11]. Implying that quantum heat engines are also bounded by the Carnot efficiency. Reversibility and optimal efficiency is obtained in the quantum adiabatic limit, requiring an infinite cycle-time $\tau_{cycle}$.

In this paper we explore two quantum approaches to obtain non-vanishing power for a ‘Carnot-analog’ engine. The first type is termed Carnot-shortcut, defined by the cycle parameters of the reversible Carnot cycle (bath temperatures and external parameters). At the switching corners between strokes, the working medium is in equilibrium with the baths and out of equilibrium at intermediate times (Fig. 1 Panel (a)). Finite power is achieved by shortcut protocols to each stroke: Shortcuts To Adiabaticity (STA) [12] on the unitaries and Shortcut To Equilibrium (STE) [13] for the thermalization strokes. The cycle is performed by external driving of the system and coupling/decoupling the Working Medium (WM) from the hot and cold baths. Typically, the WM Hamiltonian does not commute with itself at different times $\left[{\hat{H}_{S}\left({t}\right),\hat{H}_{S}\left({t^{\prime}}\right)}\right]\neq 0$, leading to generation of coherence and an accompanied cost in work [14, 15, 16].

An alternative approach to obtain finite power is a Quantum Endoreversible cycle (second type of ’Carnot analog’ engine), for which the WM is in a non-equilibrium state throughout the cycle. In the class of endoreversible cycles, we construct the Endo-shortcut and the Endo-global cycles, which are characterized by Branch and Global coherence operations, respectively. Branch coherence is restricted to the interior of the stroke, while at the cycle’s four corners the coherence vanishes and the WM state is of a Gibbs form, with an internal temperature $T\neq T_{bath}$. In the global coherence operation the WM state exhibits coherence throughout the whole cycle. As a result, coherence generated in one stroke continues to play a role in the subsequent strokes.

Our current aim is to explore the performance and characteristics of the Quantum Carnot-analog engines, i.e., the trade-off between power and efficiency and the role of quantum coherence in the engine operation. For this purpose, we employed an harmonic oscillator WM, allowing an explicit solution of the cycle dynamics. Recently, such a working medium has been employed to experimentally realize quantum heat engines [17, 18].

Previous studies of the Carnot cycle lacked coherence. These studies consider a Hamiltonian that commutes with itself at all times [19, 20], or cycles performed in the adiabatic [21, 22] or stochastic [6] limits. The relatively few studies on the quantum Carnot-analog cycle is in contrast to the popularity of the quantum Otto cycle. Analysis of the Otto cycle has been a major source of insight on quantum reciprocating engines [14, 23, 24, 25, 17, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35]. The vast difference in popularity between the quantum Otto and Carnot cycles arises from the difficulty to describe the open-system dynamics of non-adiabatically driven systems. Recent development of the Non-Adiabatic Master Equation (NAME) [36] and the inertial theorem [37] enable this study.

The analysis of the quantum Carnot-analog engine, demonstrates two roles of coherence in the energy representation. First, the control of coherence allows optimizing the efficiency for a finite cycle-time, minimizing quantum friction [14]. Moreover, for very short cycle-times global coherence becomes crucial to obtain non-vanishing power. This phenomena is a quantum signature [38, 39, 40], which can be unraveled from thermodynamics observables. We consider a thermodynamic quantum signature as a measurable thermodynamic quantity, such as work or efficiency, which indicates a unique quantum behaviour. Such behaviour, by definition, does not have a classical counterpart. Specifically, we find that under fast driving the coherent work extraction mechanism outperforms the stochastic work extraction ̵ֿ[38].

The paper is organized as follows: We begin by addressing the general construction and features of the various Carnot analog cycles, Sec 2. Section 3 then focuses on the Carnot-shortcut cycle, describing its construction and thermodynamics behaviour. In this framework, we discuss the role of coherence in the efficiency-power tradeoff, which naturally emerges in the performance of the Carnot-shortcut cycle. In Sec. 4 we introduce two additional Carnot analog cycles, the endoreversible shortcut and global Carnot cycles. Following the description, we compare between branch and global coherence operation modes, associated with these two cycles. Such comparison, highlights the role of coherence in the fast driving regime, and illuminates the emergence of quantum thermodynamics signatures. We conclude in Sec. 5 with a summary and compare the quantum Carnot and Otto engines.

2 Dynamics of the Carnot-analog cycle

We choose a particle of mass $m$ confined by a varying Harmonic potential as the engine’s working medium. The associated Hamiltonian reads

[TABLE]

where $\hat{Q}$ and $\hat{P}$ are the position and momentum operators, and $\omega\left({t}\right)$ is the oscillator frequency and the control parameter. The explicit time-dependence of $\omega\left({t}\right)$ defines the cycle’s protocol.

The dynamics during the cycle strokes is associated with a time-dependent completely positive map [41, 42], generated by

[TABLE]

For the adiabatic strokes ${\cal{L}}\left({t}\right)\hat{\rho}\left({t}\right)=-i\left[{\hat{H}\left({t}\right),\hat{\rho}\left({t}\right)}\right]$, generating the unitary maps of the adiabatic expansion and compression ${\cal{U}}_{ch}$ and ${\cal{U}}_{hc}$.

In the isothermal-type strokes ${\cal{L}}\left({t}\right)\hat{\rho}\left({t}\right)=-i\left[{\hat{H}\left({t}\right),\hat{\rho}\left({t}\right)}\right]+{\cal{L}}_{D}\left({t}\right)\hat{\rho}\left({t}\right)$, where ${\cal{L}}_{D}\left({t}\right)$ is an explicitly time-dependent Lindblidian [36], generating the isotherms dynamical maps, ${\cal{U}}_{h}$ and ${\cal{U}}_{c}$. Theses strokes incorporate non-adiabatic driving protocols.

The generators of the isotherms can be derived from first principles under the following restrictions: (i) Weak system-bath coupling; (ii) Fast bath dynamics relative to the system and driving [36]; (iii) The protocol complies with the Inertial Theorem [37]. Specifically, for the parametric harmonic oscillator, Eq. (1), the inertial theorem is satisfied in the limit of a slow change in the adiabatic parameter, $d\mu/dt\rightarrow 0$, where $\mu=\dot{\omega}/\omega^{2}$ . The adiabatic parameter can be large, indicating fast driving relative to the Bohr frequencies, while varying slowly. The resulting master equation is consistent with the thermodynamic laws [36, 43].

The Carnot analog cycle is constructed by concatenating four strokes:

• •

An expansion stroke (Fig. 1 Panel (a) corners $1\rightarrow 2$), decreasing the oscillator frequency. while the system is coupled to a hot bath of a temperature $T_{h}$. This stroke is termed open-expansion, referring to the fact that during this segment the working medium is an open quantum system, exchanging energy and entropy with the bath. The evolution of the WM is governed by the propagator ${\cal{U}}_{h}$.

• •

Adiabatic expansion stroke ($2\rightarrow 3$) in which the system is isolated from the baths and the oscillator frequency is decreased. The stroke is associated with the propagator ${\cal{U}}_{ch}$.

• •

Open-compression ($3\rightarrow 4$), the particle is brought in contact with the cold bath of temperature $T_{c}$ and ‘compressed’ towards a higher frequency, associated with the propagator ${\cal{U}}_{c}$.

• •

Adiabatic compression ($4\rightarrow 1$), the final stroke restores the system to its initial state while keeping the particle isolated from the baths. The stroke propagator is ${\cal{U}}_{hc}$.

Formally, the complete cycle propagator can be decomposed to stroke propagators:

[TABLE]

Each propagator ${\cal{U}}_{i}$ is a completely positive map, associated with the cycle stroke. For a closed Lie algebra, the system dynamics can be described within a restricted operator vector space. Within this framework, we can represent the state as a generalized Gibbs form, which is the maximum entropy state under the constraint imposed by the observables of the Lie algebra [44, 45, 46]. The generalized Gibbs structure is preserved under the full dynamics of the cycle, this property is termed canonical invariance [47, 13, 36]. For the analysis, we choose operators motivated by thermodynamical insight: The instantaneous Hamiltonian and two additional operators to represent coherence the Lagrangian $\hat{L}\left({t}\right)$ and ${\hat{C}}\left({t}\right)$ (Cf. D). These operators do not commute with $\hat{H}\left({t}\right)$, together they define the coherence measure [24]:

[TABLE]

Repeated operation of the cycle propagator, Eq. (3), converges to a unique limit cycle [48]. The thermodynamic analysis is then carried out on the limit cycle in terms of work $W$ and heat $Q$. External work preformed on the WM is defined in terms of the instantaneous power, $\rm{tr}\left({\dot{\hat{H}}\hat{\rho}}\right)$, [8]

[TABLE]

and heat is obtained from the first law of thermodynamics $E=W+Q$.

3 Carnot-shortcut cycle

The Carnot-shortcut cycle is characterized by the four corners of the ideal Carnot cycle. Utilizing shortcut protocols, the WM transitions between the thermal states associated with the cycle corners in finite time. Starting from a thermal state the non-adiabatic protocols induce coherence generation and excitations leading to a non-equilibrium state. Generated coherence is then transferred to energy at the end of the stroke leading to the desired thermal state. In the long cycle time limit, the state follows the adiabatic change in the Hamiltonian and the cycle converges to the reversible Carnot cycle.

The switching points between strokes are defined in terms of the bath temperatures $T_{c},\,T_{h}$ and oscillator frequencies $\omega_{1},\omega_{2},\omega_{3}$ and $\omega_{4}$. Defining the compression ratio ${\cal{C}}=\omega_{max}/\omega_{min}=\omega_{1}/\omega_{3}$ , the cycle is completely determined by $\omega_{min}=\omega_{3}$, $\cal C$ and bath temperatures.

At the four corners the WM is in a thermal state: $\hat{\rho}_{th}\left({\omega,T}\right)\equiv Z^{-1}\exp\left({-\hat{H}\left({\omega}\right)/k_{B}T}\right)$, where $Z$ is the partition function (see Fig. 1 Panel (a)). The adiabatic strokes conserve the working medium’s entropy. Therefore, the populations satisfy $n_{2}=1/\left({\exp\left({\frac{\hbar\omega_{2}}{k_{B}T_{h}}}\right)-1}\right)=n_{3}=1/\left({\exp\left({\frac{\hbar\omega_{3}}{k_{B}T_{c}}}\right)-1}\right)$ and $n_{1}=n_{4}$, where $n_{i}$ is the population at the $i$’th corner. This leads to the condition

[TABLE]

From the cycle definition $\omega_{1}>\omega_{2}$ and using Eq. (6) we obtain a lower bound for the compression ratio

[TABLE]

The Carnot-shortcut obtains non-vanishing power by utilizing Shortcut To Equillibration (STE) protocols [13] during the open-expansion and open-compression strokes (see A), and Shortcuts To Adiabaticity (STA) protocols for the adiabatic expansion and compression (see B).

Shortcut to equilibrium protocols are designed to manipulate an open quantum system between thermal states. The construction of the these protocols is based on the inertial theorem [37] and the NAME [36]. Starting from a stationary Gibbs state the driving is first accelerated (complying with the inertial theorem) and then decelerated, leading to the target Gibbs state at the end of the stroke. In the present open-expansion process we modify the frequency from $\omega_{1}$ to $\omega_{2}$ while the system interacts with a bath at temperature $T_{h}$ (or in the open-compression from $\omega_{3}$ to $\omega_{4}$ at bath temperature $T_{c}$). This protocol balances coherence generation and dissipation to achieve a Gibbs state, diagonal in the energy representation.

The protocol duration can be varied within the framework of the inertial approximation with negligible deviations from optimal fidelity [13]. Once the protocol duration is cut short, the STE requires a rapid change of the Hamiltonian, generating more intermediate coherence. This leads to larger dissipation, resulting in an increase of the entropy production and work cost.

On the adiabats, frictionless protocols are constructed employing the Lewis-Riesenfeld invariant [49, 12, 50]. These protocols vary the oscillator frequency non-adiabatically, generating coherence at intermediate times and storing energy in the WM. At the end of the stroke, the coherence is fully extracted, leading to a total zero work cost. Since all stored energy is retrieved when the protocol is completed, we consider these processes as an analog to catalysis [24]. In principle, these protocols can be achieved almost instantaneously. However, this implies a temporary storage of an infinite amount of energy in the WM [12, 50, 51]. To comply with practical physical considerations, we choose a constant stroke duration for the adiabats, which is consistent with the inertial theorem. Moreover, the cycle-time is dominated by the open-strokes, therefore the time allocated to the adiabats does not alter the performance qualitatively. A more complete description accounts for the power dissipated in the controller [52, 53], this contribution is out of the scope of the present analysis.

Control of the cycle time allows a transition between two extreme operation modes. Short cycle times produces finite power ${\cal{P}}=-W/\tau_{cycle}$, on account of a reduction in efficiency, due to quantum friction associated with fast driving. The non-adiabatic driving generates coherence, which is dissipated to the bath. As a result, the useful work decreases, leading to a decrease in efficiency. In contrast, the long cycle time regime optimizes the efficiency while leading to a vanishing power output. Such cycle performance demonstrates the trade-off between power and efficiency.

Correspondence between the Carnot-shortcut cycle and the ideal classical result is obtained in the quantum adiabatic limit (diverging stroke times). In this limit, the WM state remains on the energy shell along the whole cycle and the STE and STA strokes converge to reversible isothermals and adiabats. The optimal work extraction becomes

[TABLE]

where $\Delta\omega_{ij}=\omega_{i}-\omega_{j}$. In the high temperature limit Eq. (9) obtains the elegant form

[TABLE]

3.1 Performance of the Carnot-shortcut cycle

The performance of the Carnot-shortcut cycle is analysed by varying the duration of the open strokes, while keeping the time allocation of the adiabatic strokes fixed. We find two characteristic operational modes: Engine and dissipator. These are associated with different driving speed regimes. For ‘slow’ to ‘moderate’ driving speeds (‘long’ or ’medium’ cycle-times) the cycle operates as an engine (negative work output and positive efficiency, using the convention that outgoing energy is negative). For decreasing cycle-times the work output of the open-expansion decreases (in absolute value) and the work required to perform the open-compression stroke increases. This leads to a reduced efficiency, see Fig. 2 Panel (a). Once the cycle-time is reduced below $\tau_{trans}=23.87\,\left({2\pi/\omega_{min}}\right)$, where $\omega_{min}=5\,\rm{a.u}$, the cycle operates as a dissipator, converting net positive work to heat, which is in turn dissipated to the cold bath. This leads to negative efficiency for $\tau_{cycle}<\tau_{trans}$, as shown in Fig. 2 Panel (a).

The transition between dissipator and engine operational modes is caused by the large energy dissipation, which accompanies fast driving. Fig. 1 displays a map of the cycle in the energy frequency plane. For slow driving, Panel (a), the open-strokes follow the isotherms, leading to a Carnot-analog engine that obtains close to optimal efficiency $\eta_{C}=0.375$. When the cycle-time shortens, the energy deviates significantly from the isotherms, Panel (b), and the cycle operates as a dissipator. Decrease in efficiency is attributed to fast driving, that requires large generation of coherence. Consequently, dissipation is increased, eventually canceling the useful extracted work. Geometrically, the reduction in work can be witnessed as a reduction in the area confined by the cycle, in the $\langle{\hat{H}}\rangle,\omega$ plane.

Since the adiabats are preformed by frictionless protocols, the cycle’s operational mode is dictated by the open strokes. In these strokes, decreasing the stroke times costs work, eventually damaging the performance. Such an effect is naturally described in the simplified asymptotic large cycle time limit. Asymptotically, the work cost scales as $1/t$ [13] and the power behaves as, (Cf. C)

[TABLE]

where $W_{C}$ is the ideal work, Eq. (9), and $F$ is the friction action (see C). The three terms of Eq. (11) are positive for an engine operation, inferring that there exists a maximum in power for a finite cycle-time $\tau_{cycle}$ 111Eq. (11) can be written as a second order polynomial in the parameter $x=1/\tau_{cycle}$, with $d^{2}{\cal{P}}/dx^{2}<0$. . Furthermore, for sufficiently small $\tau_{cycle}$ Eq. (11) implies that the power output is negative and the cycle becomes a dissipator. This leads to the approximate relation between the cycle-time at maximum power $\tau^{*}$ and the transition time $\tau_{trans}$: $\tau^{*}=2\tau_{trans}$. Fig. 2 demonstrates this approximate relation. Overall, reducing the cycle-time increases the invested work during the open-compression stroke, and reduces the work output in the open-expansion, in accordance with the thermodynamic principles.

For increased cycle-times, the efficiency of the cycle improves, obtaining the Carnot bound asymptotically (see Fig. 2 Panel (a)). Optimal power is achieved for relative short cycle-times $\tau_{max{\cal{P}}}=43\left({2\pi/\omega_{min}}\right)$, (see Fig. 2 Panel (b)). We have also analyzed the performance at elevated temperatures, maintaining a constant temperature ratio. Asymptotically, the power increased linearly, Eq. (10), and the efficiency at maximum power converged to the value $\eta_{max{\cal P}}^{*}$ slightly larger than $\eta_{CA}$. This implies that this cycle does not operate in the weak dissipation limit [7].

4 Quantum signature: Endoreversible Shortcut and Global Carnot cycles

Endoreversible cycles are defined by four corners for which the WM is in a Gibbs state with temperatures $T\neq T_{bath}^{\prime}$. As a result, finite heat flow between the engine and bath occurs and the WM is never in equilibrium with the heat bath. In the following, primes designate the bath temperatures of the endoreversible cycle.

The Endo-shortcut cycle is constructed in a similar fashion to the Carnot-shortcut cycle Sec. 3, including two STE and two STA protocols. For such a cycle, the STE protocols are modified to transform the state between non-equilibrium Gibbs states, Cf. A.1. As a result, the four corners are maintained in a non-equilibrium Gibbs state $\hat{\rho}\left({\omega,T}\right)=Z^{-1}\exp\left({-\hat{H}\left({\omega}\right)/k_{B}T}\right)$, keeping the same frequencies ($\omega_{1}-\omega_{4}$), where $T$ is the working medium’s internal temperature which differs from the bath temperatures $T_{c}^{\prime}$ and $T_{h}^{\prime}$. The bath temperatures of the Endo-shortcut cycle satisfy $T_{c}^{\prime}<T_{c}$ and $T_{h}^{\prime}>T_{h}$, relative to the Carnot-shortcut bath temperatures $T_{c}$ and $T_{h}$, see Table 1 for a summary of the various cycle parameters. During the open-expansion the STE transfers the oscillator from $\hat{\rho}\left({\omega_{1},T_{h}}\right)$ to $\hat{\rho}\left({\omega_{2},T_{h}}\right)$, and similarly in the open-compression ($\omega_{4}\rightarrow\omega_{1}$ with $T=T_{c}$), see Fig. 1 Panel (d).

Alternatively, the Endo-global cycle operates as a Carnot-analog cycle where the strokes are performed with a pre-defined constant adiabatic speed $|\mu|=|\dot{\omega}/\omega^{2}|=\rm{const}$. The strokes are defined by four frequencies $\omega_{1}^{g},\omega_{2}^{g},\omega_{3}^{g},\omega_{4}^{g}$, and $\mu$. These determine the stroke duration $t_{f}$ and protocols $\omega\left({t}\right)=\omega_{i}/\left({1-\mu\omega_{i}t}\right)$ for a stroke starting at $\omega_{i}$ and ending at $\omega\left({t_{f}}\right)$, Cf. D.

Propagators of the free dynamics, Eq. (1), can be obtained explicitly in terms of an operator basis in Liouville space Cf. D. These propagators construct the WM unitary transformations of the adiabats. Next, the same free dynamics solution is used to derive the NAME, employed to generate the propagator of the open-expansion, ${\cal{U}}_{h}$, and open-compression, ${\cal{U}}_{c}$, strokes [36]. Combining the four strokes forms the Global cycle.

In the following, we analyze an Endo-global Carnot cycle operating between baths of temperatures $T_{c}$ and $T_{h}$, with frequencies $\omega_{1}^{g}=T_{h}^{g}\left({\omega_{1}/T_{h}}\right)$, $\omega_{3}^{g}=T_{c}^{g}\left({\omega_{3}/T_{c}}\right)$, $\omega_{2}^{g}=\omega_{3}^{g}\left({T_{h}^{g}/T_{c}^{g}}\right)$ and $\omega_{4}=\omega_{1}^{g}(T_{c}^{g}/T_{h}^{g})$, where $T_{c}^{g}=5.25$ and $T_{h}^{g}=7.75$. The frequencies are chosen so as to comply with an ideal Carnot cycle (or Carnot-shortcut cycle) operating between temperatures $T_{c}^{g}$ and $T_{h}^{g}$, in the long cycle-time limit. In this limit, the state at the four corners becomes isoentropic with the states of the Carnot-shortcut cycle. As a result, the power of the two cycles coalesce at long cycle times, see Fig. 2 Panel (b). Thermodynamic analysis is carried out on the limit-cycle [48], defined by the cycle parameters above.

4.1 Cycle performance comparison

Performance of the three cycles is compared with an emphasis on the role of coherence. The cycle output is determined by the operational mode and the cycle parameters: Cycle frequencies, bath temperatures and cycle-time, Cf. Table 1. The Shortcut cycles (Carnot and endoreversible) are characterized by a branch coherence operation, where coherence is created ‘locally’ during each stroke, with initial and final diagonal Gibbs states. This can be seen in Fig. 5, where the dashed lines represent the coherence of the Carnot shortcut cycle which vanish after each stroke. Conversely, in the limit-cycle of the Endo-global engine, coherence is maintained throughout the cycle (globally) (non-vanishing continuous line in Fig. 5). As a result, coherence generated in one stroke can be utilized along the subsequent strokes.

This comparison is unbiased, since the same bath temperatures are considered for the Carnot-shortcut and Endo-global cycles. In the Endo-shortcut cycle, despite of the larger temperature difference there is no foreseen advantage.

The Carnot-shortcut cycle shows a superior performance at moderate and long cycle-times, Figs. 2 and 3. In this operational regime, the efficiency, power, maximum power, and efficiency at maximum power exceed both the Endo-global and Endo-shortcut cycles. However, for short cycle-times, both shortcut cycles become dissipators, producing negative power, while the Endo-global cycle continues to operate as an engine.

This behaviour compels us to reexamine the question: When is a quantum heat engine quantum? An intuitive approach defining a quantum signature, is that by observing the engine operation it will cease to produce useful output. To be more precise we can imagine continuously monitoring the energy of the WM on the adiabatic strokes. As a result, the state will relax to a diagonal state in the energy representation and coherence vanishes. The back action of the weak measurement is therefore equivalent to pure dephasing [54, 55], described by the master equation:

[TABLE]

where the parameter $\gamma_{d}$ scales the dephasing rate. The explicit derivation is found in F. The influence of the dephasing rate on the engine efficiency is shown in Fig. 6. For long cycle-times where the nonadiabatic effects are small the influence of dephasing is minor. Conversely, for short cycle-times the engine operation requires coherence and the strong dephasing nulls the power output. Similar behaviour has been observed in the sudden Otto cycle and the two-stroke NV engines [24, 56, 39].

Another viewpoint on quantum signatures arises from a comparison with the classical counterpart. Introducing strong dephasing effectively leads to a pure stochastic mode of operation, where the system is solely characterized by its population in the energy representation and the energy levels. The improved efficiency, in the small cycle time regime, demonstrates that for fast driving the coherent operation mode outperforms the stochastic work extraction mechanism. This highlights the quantum advantage in this regime. The relative increase of efficiency under coherent operation also constitutes a thermodynamic measurable quantity that indicates the existence of a pure quantum phenomena. Therefore, it is considered as a quantum signature. In essence, coherence emerges when the WM state is in a superposition of eigenstates of the instantaneous Hamiltonian, implying that the superposition nature is the origin of the increased efficiency and quantum signature.

4.2 Geometric representation of coherence

To illuminate the interplay between coherence and energy throughout the Carnot analog cycles, we present a geometrical representation of the WM state. Generally, a Gaussian state of the Harmonic oscillator is fully characterized by a basis of four operators forming a closed Lie algebra. Moreover, assuming an initial Gaussian form, the WM state is canonical invariant under the equations of motion of the various strokes [47, 24]. We chose a time-dependent operators basis $\{\hat{H}\left({t}\right),\hat{L}\left({t}\right),\hat{C}\left({t}\right)\}$, see Sec. 2 and D. This leads to the geometrical representation in terms of the vector $\vec{v}\left({t}\right)=\{\langle{\hat{H}}\rangle,\langle{\hat{L}}\rangle,\langle{\hat{C}}\rangle\}^{T}$, shown in Figure 7.

Coherence in the energy basis is related to the operators $\hat{L}\left({t}\right)$ and $\hat{C}\left({t}\right)$ which do not commute with $\hat{H}\left({t}\right)$, D Eq. (51). The shortcut cycles, shown in Panels (a) and (b) display significant coherence generation, with an overshoot in energy during the open-strokes. The Lewis-Riesenfeld protocols of the adiabatic strokes lead to $\langle{\hat{L}\left({t}\right)}\rangle\approx 0$ [50]. In this cycle, all the connecting corners between strokes are on the line $\langle{\hat{L}}\rangle=\langle{\hat{C}}\rangle=0$. The trajectory of the Endo-global cycle encircles the zero coherence line $\langle{\hat{L}}\rangle=\langle{\hat{C}}\rangle=0$, while never coinciding with it. As expected, short cycle-times lead to greater coherence, Panel (c).

5 Discussion

A family of finite power quantum Carnot-analog engines is investigated to determine the role of coherence in the engine operation. In the construction of the models, we have overcome the difficulty of describing thermalization in the presence of non-adiabatic driving [36], including an explicit description of the working medium dynamics. This enables a comparison to other finite-time quantum cycle as the Otto cycle [24].

Carnot and Otto engines have been the keystones in the study of quantum heat devices. These engines differ by their thermalization strokes and share their unitary adiabatic strokes. For a fair comparison, we assume the same working medium and bath temperatures. The Otto and Carnot cycles have different permissible compression ratio range. The Otto cycle is constrained by $\frac{T_{h}}{T_{c}}\geq{\cal{C}}>1$, while the compression ratio of the Carnot cycle obeys ${\cal{C}}>\frac{T_{h}}{T_{c}}$. Optimal work for the Otto cycle is achieved when ${\cal{C}}=\sqrt{\frac{T_{h}}{T_{c}}}$ and vanishes in limits ${\cal C}=1$ and ${\cal C}=\frac{T_{h}}{T_{c}}$. On the other hand, the power of the Carnot-analog has a local optimum, and the work output diverges with ${\cal C}$, Eq. (10).

A major performance measure of engines is the efficiency. Otto and the Carnot cycles have different efficiency characteristics, in particular, their dependence on the compression ratio and bath temperatures. The Otto efficiency $\eta_{O}=1-{\cal C}^{-1}$ converges to $\eta_{C}$ in the limit of zero work. Moreover, optimal efficiency increases conjointly with the work, achieving $\eta_{O}=\eta_{CA}$ at the optimum power point in the high temperature limit [47]. Conversely, for the Carnot-shortcut cycle, efficiency $\eta_{C}$ is constant as function of compression ratio, while the work monotonically increases with $\cal C$. For this cycle the power obtains its optimum at a finite compression ratio. In addition, the efficiency at maximum power is found to be greater than the Curzon-Ahlborn efficiency $\eta_{max{\cal P}}>\eta_{CA}$.

Control of coherence is the key to obtain finite-power in quantum engines. Shortcut protocols utilize coherence temporally to accelerate the dynamics and enhance the power. On the adiabatic unitary strokes, shortcuts (STA) can be obtained with no direct cost, Sec. 2 [28, 29]. For the thermalization strokes, STE protocols enable swift equilibration. This is achieved by generating significant coherence which increases the dissipation of energy and information to the bath. As a result, such a speedup is accompanied by an additional work cost. These protocol are incorporated in the operation of the Carnot-shortcut and Endo-shortcut. The dissipated potential work in the cycle explains the mechanism leading to the power-efficiency tradeoff, and the existence of a minimum cycle time [47, 57, 58, 59]. Similar STE protocols can also be used to accelerate the thermalization in the Otto and Stirling cycle. In addition, STE protocols can be experimentally realized in engines with harmonic working fluid as well as a spin system [17, 18, 60].

The Endo-global cycle is characterized by non-vanishing coherence throughout the cycle operation, Fig. 5. This is in contrast to the shortcut cycles, which are diagonal in energy representation at the four corners of the cycle. Global coherence enables engine operation at short cycle-times, where the Endo-Shortcut and Carnot-Shorcut become dissipators. This indicates a quantum signature[39]. In the presence of pure-dephasing, the coherence vanishes, reducing power and efficiency. Eventually, strong dephasing transforms the cycle from an engine to a dissipator, a signature of a pure quantum operational mode, see Fig. 6.

To conclude, the motto of this study follows Carnot’s tradition of learning by example. The analysis of a solvable cycle has unraveled the role of coherence. On one hand, coherence is responsible for quantum friction, and on the other hand it allows rapid thermalization and a coherent mode of operation.

We thank Peter Salamon for inspiring this study. We also gratefully acknowledge KITP for their hospitality, and the support of the Adams Fellowship Program of the Israel Academy of Sciences and Humanities. Partial funding was supplied by the National Science Foundation under Grant No. NSF PHY-1748958 and the Israel Science Foundation Grant No. 2244/14.

Appendix A Shortcut to equilibration protocol

The STE protocol rapidly transfers the working medium between two Gibbs states with different entropies. Construction of the protocol relies on the Non-Adiabatic Master Equation (NAME) [36] and the Inertial Theorem [37]. Derived from first principles, the NAME describes the reduced dynamics of a non-adiabatically driven open quantum system in the weak coupling Markovian regime. The NAME incorporates, as a limit, both the adiabatic and periodic Floquet master equations.

The derivation of the NAME requires an explicit solution of the closed system dynamics. Utilizing the inertial theorem, we obtain the free propagator $\hat{U}_{S}\left({t}\right)$ for slowly ‘accelerated’ external control. For this study, we considered a harmonic oscillator working medium, where the control is achieved by varying the potential. In this case, the inertial condition is associated with a slow change in the adiabatic parameter $\mu=\dot{\omega}/\omega^{2}$. The explicit solution of free dynamics propagator allows to transform the system-bath coupling to the interaction representation. Using the Born-Markov approximation we derive the generator of the open system dynamics. For this case, we consider an Ohmic Boson bath with Markovian properties. The derivation leads to the reduced system dynamics in the interaction representation

[TABLE]

Here, the density operator in the interaction picture reads $\tilde{\rho}_{S}\left({t}\right)=\hat{U}_{S}\left({t,0}\right)\hat{\rho}_{S}\left({t}\right)\hat{U}_{S}^{\dagger}\left({t,0}\right)$.

[TABLE]

where $N$ is the occupation number of the Bose-Einstein distribution, $\kappa=\sqrt{4-\mu^{2}}$, and $\alpha$ is a modified frequency, determined by the non-adiabatic driving protocol [36]. In terms of the oscillator frequency, the modified frequency is given by

[TABLE]

The Lindblad jump operators become $\hat{b}\equiv\hat{b}\equiv\hat{b}\left({0}\right)=\sqrt{\frac{m\omega\left({0}\right)}{\kappa\hbar}}\frac{\kappa+i\mu}{2}\left({\hat{Q}\left({0}\right)+\frac{\mu+i\kappa}{2m\omega\left({0}\right)}\hat{P}\left({0}\right)}\right)$.

In the adiabatic driving limit, $\mu\rightarrow 0$, the Lindblad jump operators converge to the adiabatic creation and annihilation operators $\hat{b}^{\dagger},\hat{b}\rightarrow\hat{a}^{\dagger},\hat{a}$. Therefore, in this limit, Eq. (16) reduces to the adiabatic master equation.

The completely positive map generated by Eq. (16) preserves the Gaussian form. This property is termed canonical invariance [44, 61, 47, 24]. Formally, the Gaussian state can also be expressed in a product form

[TABLE]

where $Z={\tr}\left({e^{\gamma\left({t}\right)\hat{b}^{2}}e^{\beta\left({t}\right)\hat{b}^{\dagger}\hat{b}}e^{\gamma^{*}\left({t}\right)\hat{b}^{\dagger 2}}}\right)$, the operators $\hat{b}^{\dagger}\hat{b},\,\hat{b}^{2},\hat{b}^{\dagger 2}$ vary with $\mu$. Parameters $\gamma$ and $\beta$ are time-dependent functions.

For an initial thermal state, the construction is simplified and the reduced system remains in the following form throughout the entire evolution

[TABLE]

with initial conditions $\beta\left({0}\right)=-\frac{\hbar\omega\left({0}\right)}{k_{B}T}$ and $\mu\left({0}\right)=0$. Substituting Eq. (20) into Eq. (16) leads to a non-linear differential equation for $\beta\left({t}\right)$

[TABLE]

Equation (21) constitutes the basis for the STE control scheme. We define $y\left({t}\right)\equiv e^{\beta\left({t}\right)}$ and guess a polynomial solution for $y$. This solution should transfer an initial thermal state of frequency $\omega\left({0}\right)$ to a final thermal state with a frequency $\omega\left({t_{f}}\right)$, which leads to the boundary conditions $\beta\left({0}\right)=-\frac{\hbar\omega\left({0}\right)}{k_{B}T}$, $\beta\left({t_{f}}\right)=-\frac{\hbar\omega\left({t_{f}}\right)}{k_{B}T}$ and $\mu\left({0}\right)=\mu\left({t_{f}}\right)=0$. Furthermore, the condition on $\mu$ implies that the state and protocol are stationary at initial and final times, leading to $\dot{\beta}\left({0}\right)=\dot{\beta}\left({t_{f}}\right)=\ddot{\beta}\left({0}\right)=\ddot{\beta}\left({t_{f}}\right)$.

A fifth-degree polynomial is sufficient to comply with all the boundary conditions. The solution reads

[TABLE]

where $c_{3}-c_{5}$ are determined from the boundary conditions $y\left({0}\right)=e^{\beta\left({0}\right)}$, $y\left({t_{f}}\right)=e^{\beta\left({t_{f}}\right)}$, and $\dot{y}\left({0}\right)=\dot{y}\left({t_{f}}\right)=\ddot{y}\left({0}\right)=\ddot{y}\left({t_{f}}\right)=0$. Next, we insert the solution for $y$, Eq. (22), into Eq. (21) and solve for $\alpha\left({t}\right)$. The control $\omega\left({t}\right)$ is obtained by solving Eq. (18) by numerical means.

A.1 STE for non-thermal initial and final states

The Endo-shortcut cycle, Sec. 4, includes open-compression and open-expansion strokes between non-equilibrium states. During these strokes, the working medium state is of the Gibbs form

[TABLE]

where the internal temperatures $T_{i}$ and $T_{f}$ differs from the bath temperatures, i.e., $T_{i}=T_{f}\equiv T_{0}\neq T$.

The control protocol $\omega\left({t}\right)$ for the open-strokes is obtained by a similar reverse-engineering method as in the case of an initial and final equilibrium states ($T_{i},T_{f}=T$). However, since the system-bath interaction leads to non-vanishing decay rates, the initial and final states are non-stationary, which implies $\dot{\beta}\left({0}\right),\dot{\beta}\left({t_{f}}\right)\neq 0$. As in the previous construction, we require a continuous change in the control protocol, associated with the restriction $\dot{\omega}\left({0}\right)=\dot{\omega}\left({t_{f}}\right)=0$. Substituting the condition $\dot{\omega}\left({0}\right)=0$ into the decay rates of Eq. (21) we obtain the initial value for $\dot{\beta}$ (this is in accordance with the dynamics of the time-independent master equation [62]). This leads to the boundary conditions for $y$: $y\left({0}\right)=y\left({0}\right)=\exp\left({-\frac{\hbar\omega\left({0}\right)}{k_{B}T_{0}}}\right)$, $\dot{y}\left({0}\right)=\dot{\beta}\left({0}\right)\exp\left({-\frac{\hbar\omega\left({0}\right)}{k_{B}T_{i}}}\right)$, $y\left({t_{f}}\right)=\exp\left({-\frac{\hbar\omega\left({t_{f}}\right)}{k_{B}T}}\right)$, $\dot{y}\left({t_{f}}\right)=\dot{\beta}\left({t_{f}}\right)\exp\left({-\frac{\hbar\omega\left({t_{f}}\right)}{k_{B}T_{0}}}\right)$ and $\ddot{y}\left({0}\right)=\ddot{y}\left({t_{f}}\right)=0$. Following a similar derivation as the previous section, we introduce a fifth order polynomial that satisfies the boundary conditions to obtain $\omega\left({t}\right)$.

Appendix B Adiabatic strokes - Shortcut to adiabaticity protocols utilizing the Lewis Riesenfeld invariant

The (thermodynamic) adiabatic strokes are achieved utilizing shortcut to adiabaticity (STA) protocols. These transform a diagonal state in the energy basis, of a frequency $\omega_{i}$, to a state with the same population, with a final frequency $\omega_{f}$. These protocols are engineered utilizing the Lewis-Riesenfeld invariant [49]. We follow a similar procedure as presented in Refs. [49, 12] to construct the STA protocols.

We introduce an Hermitian invariant $\hat{I}$ for the harmonic oscillator algebra, SU(1,1). Generally, such an invariant can be expressed as sum of the algebra operators

[TABLE]

and must satisfy the condition

[TABLE]

To obtain transition-less driving we desire a protocol for which the invariant commutes with the Hamiltonian at initial and final times:

[TABLE]

Since the operators share a common eigenstate basis at these times and the eigenvalues of $\hat{I}$ are stationary, the engineered protocol induces transition between two diagonal states in the energy representation.

Substituting Eq. (24) into Eq. (25) leads to three coupled differential equations for the time-dependent coefficients

[TABLE]

By defining $\beta\equiv\sigma^{2}$ and conducting some algebraic manipulations (see Ref [49]) Eqs. (30) can be represented in terms of a single differential equation

[TABLE]

Solving for $m^{2}\omega^{2}\sigma+m^{2}\ddot{\sigma}$ and substituting the solution into Eq. (31) gives

[TABLE]

where, $c$ is an arbitrary real integration constant. Next, we define $\sigma=c^{1/4}\rho$, and substitute the definition into Eq. (32), to obtain

[TABLE]

with the subsidiary condition (Eq. (32)

[TABLE]

This equation introduces constraints on the protocol $\omega\left({t}\right)$, which comply with the boundary conditions of $\rho\left({t}\right)$. The strategy to engineer a transition-less control protocol is to choose $\rho\left({t}\right)$ such that relations Eq. (26) are satisfied. This leads to the following boundary conditions

[TABLE]

Equation (34) is now solved by introducing a polynomial solution. A fifth order polynomial in $t$ is sufficient to satisfy the conditions of Eq. (38). We substitute the polynomial solution into Eq. (34) and solve for $\omega\left({t}\right)$.

To obtain the expectation values of $\hat{H},\,\hat{L}$ and $\hat{C}$, we introduce the eigenstates of $\hat{I}\left({t}\right)$. These obey the eigenvalue equation

[TABLE]

with time-dependent eigenstates and time-independent eigenvalues $\lambda$. Next, we define the annihilation operator $\hat{c}={2\hbar}^{-1}\left({\hat{Q}/\rho-i\left({m\dot{\rho}\hat{Q}-\rho\hat{P}}\right)}\right)$ and matching annihilation operator. These operators satisfy: $\hat{c}\left|\lambda\right>=\sqrt{\lambda}\left|\lambda-1\right>$ and $\hat{c}^{\dagger}\left|\lambda\right>=\sqrt{1+\lambda}\left|\lambda+1\right>$. By expressing $\hat{H},\,\hat{L}$ and $\hat{C}$ in terms of $\hat{c}$ and $\hat{c}^{\dagger}$ one obtains

[TABLE]

The final step is to sum over the contribution of each state $\left|\lambda\right>$. For an initial Gibbs state the energy reads

[TABLE]

The derivation for $\langle{\hat{L}}\rangle$ and $\langle{\hat{C}}\rangle$ follows a similar procedure.

Appendix C Friction action

Shortcut to equilibrium processes rely on non-adiabatic driving of an open quantum system. The driving incorporates the dissipative and unitary dynamics to lead the system towards a target thermal state. These protocols accelerate the system thermalization rate by generating coherence at intermediate times (transforming energy to coherence), and terminating them at the end of the protocol. The protocol duration $\tau_{stroke}$ of the STE can be varied within the inertial approximation, with a negligible influence on the final fidelity. For increasing protocol duration, the generation of coherence reduces, converging to the quantum-adiabatic result in the limit $\tau_{stroke}\rightarrow\infty$. As a result, as the protocol duration increases less coherence dissipates to the bath and the work output improves (reduced in the convention of $W<0$ for work extraction). Asymptotically, the work output scales as $\tau_{stroke}^{-1}$ [13], therefore, one can introduce the ‘friction action’ $F$ and express the total cycle work output $W$ as a function of the ideal work $W_{C}$ and $F$

[TABLE]

In the quantum adiabatic limit the process is optimal and the work output converges to $W_{C}$, Eq. (9). Substituting Eq. (49) into the expression for the power output ${\cal{P}}=-W/\tau_{cycle}$, leads to the asymptotic relation

[TABLE]

Equation (50) is consistent with the general argument that the frictional forces should be independent of the sign of $\dot{\omega}$. Hence, to lowest order the power against friction is proportional to $\mu^{2}$. Since $\mu\propto 1/\tau_{cycle}$ we expect that the assymptotic power against friction scales as Eq. (50) [15].

Appendix D Endo-global cycle construction

The Endo-global cycle is constructed by combining two open strokes and two adiabats, that operate at a constant adiabatic speed. During the open strokes, the working medium dynamics is generated by the non-adiabatic master equation Eq. (16). This equation incorporates both the dissipative and unitary effects. It is exact for Markovian dynamics in the weak coupling regime and when the external driving is slow relative to the bath dynamics (see Ref. [36]). The condition of constant adiabatic speed, $|\mu|=\rm{const}$, leads to protocols of the following form $\omega\left({t}\right)=\omega_{i}/\left({1-\omega_{i}\mu t}\right)$, with initial frequency $\omega_{i}$.

The unitary dynamics are given in terms of an operator basis in Liouville space, $\vec{v}\left({t}\right)=\{\hat{H}\left({t}\right)$, $\hat{L}\left({t}\right),\hat{C}\left({t}\right),\hat{I}\}^{T}$. Here, $\hat{H}\left({t}\right)$ is the Hamiltonian, Eq. (1),

[TABLE]

is the Lagrangian, and

[TABLE]

is the position-momentum correlation operator and $\hat{I}$ is the identity operator. These operators completely determine the dynamics of the harmonic oscillator Gaussian state (19), and give a clear physical interpretation of energy and coherence,

[TABLE]

The solution for $\vec{v}\left({t}\right)$ is given by [24]:

[TABLE]

with

[TABLE]

where $\kappa=\sqrt{4-\mu^{2}}$ and $c=\cos\left({\kappa\theta\left({t}\right)}\right)$, $s=\sin\left({\kappa\theta\left({t}\right)}\right)$. The free propagation, governed by Eq. (53), mixes coherence and populations due to non-adiabatic driving [24].

In order to judge the performance for varying cycle-times, we construct different cycles with the same cycle frequencies, and vary $|\mu|$. The adiabatic parameter determines the stroke duration according to $\tau_{stroke}=\mu^{-1}\left({\omega_{f}-\omega_{i}}\right)/\left({\omega_{f}\omega_{i}}\right)$, where $\omega_{f}$ is the final frequency. The cycles are then propagated until convergence to the limit-cycle [48, 63], where the performance is evaluated.

Appendix E Coherence measures

Coherence is associated with non-diagonal elements in the energy representation. This means that a state possessing coherence is non-stationary under the free Hamiltonian dynamics. There have been many proposals to quantify coherence [64]. In this context, two measures have been employed to analyze quantum heat engines: (i) Divergence, which becomes the difference between the energy and the von Neumann entropies [65, 66, 56, 35], and (ii) The algebraic definition of coherence, utilized here, Eq. (51).

Appendix F Addition of pure dephasing

Pure dephasing is introduced during the adiabats of the Endo-global cycle to evaluate the importance of coherence in the engine operation. Such dephasing can be caused by noise in the driving field and weak measurement [67, 68, 55, 69]. A Lindbladian describing dephasing is added to the free dynamics. In the Heisenberg picture the master equation becomes

[TABLE]

Here, the last term induces pure dephasing in the instantaneous energy basis, with a dephasing strength $\gamma_{d}$. We solve the dynamics, by representing the system in terms of the operator basis $\vec{v}\left({t}\right)=\{\hat{H}\left({t}\right)$, $\hat{L}\left({t}\right),\hat{C}\left({t}\right),\hat{I}\}^{T}$, Cf. D. Substituting the basis operators into Eq. (55) we express the dynamics of $\vec{v}\left({t}\right)$ in a matrix-vector notation

[TABLE]

where $\mu$ is the adiabatic parameter, $\cal{I}$ is the identity matrix and $\cal M$ is given by

[TABLE]

Equation (56) is solved with a standard Runge-Kutta-Dormand-Prince propagator, leading to the system dynamics in the presence of pure-dephasing.

The reconstruction of the density operator, Eq. (19), from the $\{\hat{H}\left({t}\right),\hat{L}\left({t}\right),\hat{C}\left({t}\right),\hat{I}\}$ operator basis assumes that the working medium is described by a generalized canonical state. Strictly, the dephasing dynamics, Eq. (55), does not conserve the Gaussian structure. Nevertheless, for the cases studied, where the coherence is relatively small, the Gaussian state is a valid representation.

Appendix G Model parameters

The cycle parameters are described in Tables 1 and 2.