Scaling-up quantum heat engines efficiently via shortcuts to adiabaticity

The finite-time operation of a quantum heat engine that uses a single particle as a working medium generally increases the output power at the expense of inducing friction that lowers the cycle efficiency. We propose to scale up a quantum heat engine utilizing a many-particle working medium in combination with the use of shortcuts to adiabaticity to boost the nonadiabatic performance by eliminating quantum friction and reducing the cycle time. To this end, we first analyze the finite-time thermodynamics of a quantum Otto cycle implemented with a quantum fluid confined in a time-dependent harmonic trap. We show that nonadiabatic effects can be controlled and tailored to match the adiabatic performance using a variety of shortcuts to adiabaticity. As a result, the nonadiabatic dynamics of the scaled-up many-particle quantum heat engine exhibits no friction and the cycle can be run at maximum efficiency with a tunable output power. We demonstrate our results with a working medium consisting of particles with inverse-square pairwise interactions, that includes noninteracting and hard-core bosons as limiting cases.


I. INTRODUCTION
Quantum thermodynamics resembles a fruitful crucible of research fields where the foundations of physics, information science and statistical mechanics merge [1,2]. It is further spurred by the development of quantum technologies that have facilitated the realization and control of thermal machines and related devices exhibiting quantum dynamics. A prominent example is that of quantum heat engines (QHE) that transform thermal energy into mechanical work.
In both classical and quantum domains, the performance of heat engines is characterized by the efficiency and power of the cycle. Studies to date have been limited to the optimization of thermodynamic cycles operating with a single-particle working medium. After the pioneering works [3,4], the quantum Otto cycle has been analyzed in detail [5][6][7][8][9][10][11][12]. It is of relevance to current experimental efforts aimed at realizing a QHE with a single trapped ion [13][14][15]. Nonetheless, it is worth pointing out that a universal behavior emerges among different types of cycles in the limit of small action per cycle [16]. These works show that when a single-particle QHE is operated in a finite amount of time, nonadiabatic excitations act as quantum friction, reducing the efficiency of the engine. As a result, the maximum efficiency is achieved for long cycle times, in the adiabatic limit, when the output power of the QHE vanishes. This state of affairs is already present in classical heat engines and gave rise to the field of finite-time thermodynamics.
Scaling up a QHE arises as a natural strategy to compensate the tradeoff between efficiency and power. Yet, the analysis of QHE with a many-particle working medium remains essentially unexplored. Only recently, a few results have been reported [12,[17][18][19][20][21]. In particular, the efficiency of an adiabatic QHE operating an Otto cycle with non-interacting bosons and fermions has been studied in [12]. It was shown that signatures of quantum statistics are only relevant when the working medium is confined in an anharmonic trap, with an inhomogeneous energy level spacing. More recently, we analyzed a nonadiabatic quantum Otto cycle with an infinite family of quantum fluids as a working medium. Under harmonic confinement, it was shown that nonadiabatic many-particle quantum effects can enhance the efficiency of the cycle operating at maximum output power by up to 50% of the single-particle counterpart [17]. This enhancement was demonstrated in the limit of very fast driving, with a sudden trap frequency change accounting for the compression and expansion strokes. Further efforts in many-particle thermodynamics have shown that the performance of a heat engine can be boosted when the working substance is in the vicinity of a phase transition, an effect that holds in the classical case [21].
An alternative approach to boost the performance of QHE resorts on the use of shortcuts to adiabaticity (STA) that mimic adiabatic dynamics without the requirement of slow driving [22]. STA constitutes a disruptive paradigm in finite-time thermodynamics that avoids the need to sacrifice efficiency for power. This is achieved by tailoring excitations during the nonequilibrium process. STA applicable to the expansion and compression strokes often used in thermodynamic cycles has been extensively studied in recent years, both theoretically [23][24][25][26][27][28][29][30][31][32][33][34][35] and experimentally with ultracold atomic gases [36,37]. The development of STA applicable to many-body systems has been reported in [29,31,[38][39][40][41][42][43] and further demonstrated in the laboratory [44]. At the single-particle level, the application of STA to quantum thermodynamics was pioneered in [9,10]. It was shown that STA can assist the finite-time operation of a QHE with a single trapped particle as a working medium to achieve a tunable output power and zero friction. Similar applications have been discussed in quantum refrigerators [24].
In this manuscript, we combine both strategies to engineer a scaled-up many-particle QHE whose operation is assisted by STA to boost its performance, maximizing its efficiency and enhancing the output power. In particular, we analyze the finitetime thermodynamics of a QHE with a quantum fluid as a working medium. For the latter, we consider an infinite family of scale-invariant many-body systems confined in a time-dependent harmonic trap. We show that the nonequilibrium dynamics of the many-particle thermodynamic cycle can be engineered via STA to run the QHE with zero friction. The enhancement of the output power arises from the many-particle nature of the working medium and the reduced cycle time. The resulting many-particle QHE operates at the maximum possible efficiency, which is shared by classical and quantum heat engines as an upper bound.

A. Trapped Quantum Fluids as Working Media
We introduce a generalization of the single-particle QHE to the many-particle case in which the working medium consists of a quantum fluid that is confined in a time-dependent harmonic trap [17]. In particular, we focus on a broad family of quantum many-body systems described by the Hamiltonian: where N is the total number of particles and ω(t) is the trap frequency of the isotropic harmonic confinement. The only condition imposed on the pairwise interaction is that it is a homogeneous function of degree −2, i.e., V (br) = b −2 V (r). Therefore, Hamiltonian (1) accounts for a wide variety of many-body systems that include free non-interacting gases, quantum fluids with hard-core interactions, inverse-square interactions models and Bose gases with s-wave contact interactions in two spatial dimensions (d = 2), among other prominent examples. We note that under a modulation of the trapping frequency ω(t), the unitary dynamics generated by Equation (1) is governed by scale invariance [45]. As a result, the time-evolution of a manyparticle eigenstate Ψ with energy ε at t = 0 reads [17,38]: where the positive-definite function of time b = b(t) > 0 is known as the scaling factor and fulfills the Ermakov differential To describe the many-particle QHE, we shall consider an initial equilibrium state in the canonical ensemble at (inverse) temperature β described by the density matrix:ρ Tr(e −βĤ ) .
Provided that at t = 0, the state is thermal, the nonadiabatic mean-energy following a variation of the trapping potential is given by the scaling law [17]: where Ĥ (0) β is the mean energy of the initial thermal state, and the nonadiabatic factor Q * ≥ 1 reads: Here, b(t) is the solution of Equation (3) subject to the boundary conditions b(0) = 1 andḃ(0) = 0, to account for the initial equilibrium state. We note that in the adiabatic limitb ≈ 0 and: Since the mean value of the energy following an adiabatic protocol is precisely given by Ĥ (t) β ,ad = Ĥ (0) β /b 2 ad , the nonadiabatic factor Q * (t) equals the ratio between the mean nonadiabatic energy and the corresponding adiabatic value. Therefore, values of Q * (t) > 1 indicate deviations from adiabatic dynamics and can be associated with quantum friction [46], which vanishes whenever Q * (t) = 1.
For the sake of illustration, we shall consider the harmonic Calogero-Sutherland model (CSM), a specific instance of Hamiltonian (1). The CSM describes bosons confined in a time-dependent harmonic trap with inverse-square two-body interactions [47,48] where λ ≥ 0 is the strength of the two-body interactions. While the many-body eigenstates are always symmetric under permutation of the particle coordinates and Pauli's exchange statistics is bosonic for all values of λ , the CSM represents an ideal gas of non-interacting particles obeying generalized (or fractional) exclusion statistics, a concept introduced by Haldane [49] and further extended by Wu [50]. In particular, for λ = 0, the CSM reduces to the trapped ideal Bose gas. Similarly, for λ = 1, the CSM is equivalent to a one-dimensional gas of bosons in the Tonks-Girardeau regime [51,52] that obeys Fermi exclusion statistics as a result of the hard-core interactions. Therefore, when λ = 1, the occupation number of a single particle state (e.g., of the trap) is restricted to the values {0, 1}, and the thermodynamics in this regime is identical to that of non-interacting spin-polarized trapped fermions. As an upshot, the CSM allows one to account for the thermodynamics of non-interacting bosons and fermions on the same footing, given that it is governed by the exclusion statistics as opposed to the exchange statistics. More generally, for values of λ other than zero and one, the CSM describes an ideal gas of Haldane anyons or "geons", particles with fractional exclusion statistics [53]. In the canonical ensemble, an explicit computation of the equilibrium mean-energy of a thermal state of the CSM yields: where E 0 (N, λ ) denotes the ground-state energy: Hence, the effect of the interactions on the thermodynamics is completely absorbed in the renormalization of the ground-state energy of the system [54].

B. Quantum Otto Cycle and Fundamental Limits
In the Otto cycle, isentropic and isochoric strokes alternate with each other, as illustrated in Figure 1. There are two isentropes and two isochores per cycle, where energy in the form of work (W ) and heat (Q) is transferred from and into the heat engine. In order of alternation, these are: isentropic compression (W 1 ), hot isochore (Q 2 ), isentropic expansion (W 3 ) and cold isochore (Q 4 ). The isentropic strokes involve a change in the trapping frequency at constant temperature, while during an isochore, the temperature varies and the frequency kept constant. The efficiency of a heat engine is defined as the total output work per input heat: where . The heat engine is decoupled from the thermal reservoir during the isentropic strokes, when the dynamics is unitary. We use the scaling dynamics of the working medium (Equation (9)) to predict the mean energy according to Equation (5). The efficiency of the many-particle quantum heat engine run in finite time is given by: where the mean energy of the thermal states A and C corresponds to H A = H(0) β c with ω = ω 1 and H C = H(0) β h with ω = ω 2 ; see Equation (10). Equation (13) reduces in the adiabatic limit (Q * AB(CD) → 1) to the Otto efficiency: which is shared as an upper bound by both single-and many-particle quantum Otto cycles.
FIG. 1. Many-particle quantum heat engine. Quantum Otto cycle with a quantum fluid as a working medium confined in a harmonic trap with frequency varying between ω 1 and ω 2 . States at A and C are thermal, while those at B and D are generally nonequilibrium states.
When a compression (expansion) of finite duration τ (τ ) is considered, such that the frequency monotonically varies as a function of time, we have Q * AB(CD) ≤ Q * sq , and Equation (13) further implies that the finite-time efficiency η is bounded from below and above as: where η sq is the efficiency under a sudden quench of the trap frequency for which Q * sq = (ω 2 1 + ω 2 2 )/(2ω 1 ω 2 ) [55]. This applies to driving protocols, such that ω(0) = ω 1 and ω(τ) = ω 2 in the compression stroke and with ω(0) = ω 2 and ω(τ ) = ω 1 in the expansion stroke. For a monotonous frequency, we have that Q * (τ) → Q * sq in the limit τ → 0, as the τ-parameterization of a given protocol ω(t), τ ∈ [0, +∞] approaches then a smooth representation of the step function (e.g., ω 1 θ (−t) + ω 2 θ (t) for the compression). In addition, as proven in [17], the efficiency (Equation (13)) is bounded by a non-adiabatic Otto limit, independent of the number of particles N and interaction potential V .
This formula encodes the "tragedy of finite-time thermodynamics" in the many-particle setting: the maximum efficiency is achieved under slow driving, in the adiabatic limit, when the QHE operates at vanishing output power −( W 1 + W 3 )/τ c as a result of the requirement for a long cycle time τ c . By contrast, realistic engines operated in finite time achieve a finite output power at the cost of introducing nonadiabatic energy excitations that represent quantum friction and lower the efficiency. That nonadiabatic effects generally decrease the engine efficiency follows from the fact that Q * AB(CD) ≥ 1. However, we will show that this tradeoff is not fundamental in nature and can be avoided.

General Method
We first provide a general method to optimize the power of the QHE for any finite-time protocol, which will prove useful in deriving both analytical and numerical results. As shorthand, we introduce the ratio of trap frequencies: as well as the ratio of temperatures of the cold and hot reservoir: For a fixed cycle time τ c , we optimize the output power P = |W |/τ c , where the total output work |W | = −(W 1 +W 3 ) is given by: Assuming the hot bath to be at high temperature, i.e., σ h ≡ Nhβ h ω 2 1, the optimal work is reached when the frequency ratio x is a solution of the following equation: with the factor N F (x), dependent on τ and τ , given by: We note that for an adiabatic protocol, N F = 1, while for a sudden quench driving, we have N F = x 2 . Further, for an adiabatic protocol in the high temperature regime, the mean energy H A in Equation (20) is shifted by the effective zero-point energy term −E 0 (N, λ ) + E 0 (N, 1/2), which results from the expansion of H C in Equation (19) at high temperature.
In [17], it has been shown that for the sudden quench regime, the efficiency at optimal power can be drastically enhanced for a large number of particles, when σ c(h) ≡ Nhβ c(h) ω 1(2) > 1 is achieved at relatively low temperatures, keeping σ c(h) N. In this regime, one can write an explicit approximation for the total work per cycle: where we estimated the Riemann sum in (10) by a Riemann integral usinghβ ω 1 (as the corrections are bounded by N(hβ ω) 2 /2 and N cannot be too large), with: Here, we use the definition of the standard polylogarithm functions Li n (z) = ∑ ∞ j=1 j −n z j , n ≥ 1 [56]. The function µ λ (σ ) monotonically decreases for λ = 0, while it monotonically increases for all λ ≥ 1/2. For 0 < λ < 1/2, µ λ (σ ) decreases until reaching a minimum at σ = σ 0 (λ ) and increases monotonically for larger values of σ . One can interpret the function µ λ (σ c(h) ) as the quantum deviation of the mean energy H A(C) from its classical value N/β c(h) . For any given protocol, Equation (22) allows one to compute the optimal frequency ratio x that maximizes the work done per cycle.

Optimizing Adiabatic Output Power
In this section, we assume that the adiabaticity condition is fulfilled for a fixed value of the duration of the strokes τ. Optimizing the output power is then equivalent to maximizing the output work. We shall provide expressions for the maximum adiabatic output workW (N) ad and the adiabatic efficiency at optimal output power η (N) ad of the scaled-up QHE with N particles as the working medium. Our results provide a generalization of the optimization of the single-particle adiabatic QHE [13,57].
Given that for an adiabatic driving, the nonadiabatic factor reduces to unity Q * = 1, the total output work denoted by |W ad | becomes independent of λ and depends only on the thermal energy, according to Equation (19).
Assuming σ h = Nβ hh ω 2 1, we obtain: where the function E N,λ (ω 1 , β c ) is given by (10) for any λ ≥ 0. To obtain the previous equation, we use a Taylor expansion of the energy H C for small σ h . The optimal frequency ratio is found to be The maximum work and the efficiency at optimal output power read: Using the asymptotic expansion of the thermal energy, one can derive the asymptotic expressions for the optimal work and the corresponding efficiency in different regimes distinguished by the value of σ c . Taking the cold reservoir at high temperature The engine exhibits then a classical performance characterized by:W The latter equations show that many-particle effects are negligible for an adiabatic protocol when the number of particles is not too large, so that σ c(h) 1. The efficiency is then set by the Curzon-Ahlborn result, 1 − √ a, associated with a classical heat engine operated in the adiabatic limit [57].
It has been shown recently that increasing the particle number can lower efficiency drastically; see [17]. Indeed, assuming that σ h 1 and 1 ≤ σ c N, by Equations (22) and (23), we find that: where µ 1/2 (σ c ) is greater than one and increases monotonically with σ c . Consequently, the efficiency is lesser than the Curzon-Ahlborn efficiency, 1 − √ a. Equations (27a) and (27b) are valid only for small values of a = β h /β c consistent with σ h 1. This many-particle applies for a large number of particles N (hβ c ω 1 ) −1 , keeping the temperature of the cold reservoir relatively smallhω 1 β c 1. Now, considering the cold reservoir at very low temperature σ c N (i.e.,hω 1 β c 1), the thermal energy is negligible, and we find x In this quantum limit, the performance of the engine decreases drastically compared to the classical limit, as the term in the square root of Equations (28a) and (28b) is much larger than the corresponding term a = β h /β c in Equations (26a) and (26b).
From the comparison of these limiting cases, it is found that in the adiabatic limit, the efficiency and total work output are maximized in the classical regime when both the cold and hot reservoirs are at high temperature, e.g., σ c(h) 1. It is in this regime that it is most convenient to operate the scaled-up QHE, avoiding quantum corrections that lower the performance. The work per cycle is then scaled by a factor of N with respect to the single-particle QHE, and the efficiency is then set by the Curzon-Ahlborn efficiency, according to Equations (26a) and (26b). Next, we discuss how to achieve this adiabatic performance under nonadiabatic STA protocols that shorten the cycle operation time.

III. SUPERADIABATIC MANY-PARTICLE QUANTUM HEAT ENGINES
As discussed in the Introduction, it has recently been shown that it is possible to boost the performance of a single-particle QHE by reducing the cycle time, replacing the compression and expansion strokes of a quantum Otto cycle by STA protocols [8,10]. In what follows, we shall explore the optimization of the finite-time thermodynamics of a many-particle QHE using the STA for many-body systems [29,31,[38][39][40][41][42][43][44] to engineer a scale-up many-particle QHE operating at maximum efficiency and with high output power. We shall introduce two different approaches depending on the given specifications of the many-particle QHE.
A. First Approach: Finite-Time Optimization and Accidental STA Consider the operation of the many-particle QHE with given cold and hot thermal reservoirs at fixed (inverse) temperatures β c and β h . Assuming that a specific time-dependent protocol ω(t) describes the evolution of the working medium, the conditions for optimal performance can be determined by maximizing the output power for a variable and finite duration τ of the compression and expansion strokes. The duration of the isochores constitutes an overhead that adds to the total cycle time τ c and plays no role in our analysis. The possibility to engineer a superadiabatic many-particle QHE with zero friction relies on the existence of protocols ω(t) during which the nonadiabatic dynamics become adiabatic at a specific set of values for the duration of the strokes {τ n , n = 1, 2, 3, . . . }, which depends on the frequencies ω 1 , ω 2 of the trap. We refer to any such instance as an accidental STA.
To illustrate this approach, we first introduce an accidental STA protocol that interpolates between the adiabatic and the sudden quench limit (τ = 0). In particular, we choose a modulation of the trap frequency for which the adiabaticity coefficientω/ω 2 [58] remains constant during the compression and expansion strokes. We mention that this model has been introduced in the context of quantum refrigerators [24], noisy QHE [11,59] and the dynamics of a single particle in a trap [60]. Here, we introduce this modulation of the trap frequency following a different approach that suits better the many-particle QHE model Equation (9). The condition |ω(t)/ω(t) 2 | = 2γ −1 with γ > 0 is satisfied by the time-dependent frequencies: where ω 1 , ω 2 ,t 1 ,t 2 and τ are fixed by the boundary conditions ω AB (0) = ω 1 , ω AB (τ) = ω 2 and ω CD (0) = ω 2 , ω CD (τ) = ω 1 , leading to t 1 = (1 − ω 1 /ω 2 ) −1 τ and t 2 = (ω 1 /ω 2 )t 1 . Notice that the protocol is symmetric ω CD (t) = ω AB (τ − t) and can be rewritten as ω CD (t) = ω 1 ω 2 τ/(ω 1 τ + t(ω 2 − ω 1 )). The dimensionless parameter γ is given by γ = 2ω 1 t 1 and can be directly related to the duration of the stroke: via the frequency ratio x ≡ ω 1 /ω 2 . The accidental STA with finite γ interpolates between the sudden quench (γ → 0) and adiabatic limit (γ → +∞). This is verified by exact computation of the nonadiabatic factors Q * AB and Q * CD that happen to be equal and given by (see Figure 2): where x = ω 1 /ω 2 . In Figure 2a, we plot the nonadiabatic factor Equation (31) as a function of the time of evolution t for a given duration τ of the compression/expansion stroke and as a function of τ in Figure 2b. For the optimization of the QHE, the value of x is kept fixed, while τ is varied. When τ → 0, we find that Q * → (x 2 + 1)/(2x), equal to the sudden quench factor Q * sq . On the contrary, we obtain the adiabatic factor Q * ad = 1 in the limit τ → +∞. By Equation (31), for γ ≤ 1 (or equivalently τ ≤ (1 − x)/(2ω 1 )), the nonadiabatic factor decreases monotonically from the sudden quench value to a lower one that is still above unity. In Figure 3, one can see that the efficiency behaves as in a sudden quench protocol; see also Figure 2b to see the relation with the nonadiabatic factor. Yet, for γ > 1 (or equivalently τ > (1 − x)/(2ω 1 )), the nonadiabatic factor oscillates above the adiabatic value; see Figure 2b. At the contact points (obtained in [24,60] by a different method), γ n = (1 + 4π 2 n 2 / ln (x) 2 ) 1/2 , with n = 1, 2, 3 . . . , we find an STA given that nonadiabatic factor equals unity, Q * = 1. At these points, the work and the efficiency reproduce the corresponding values found in the slow driving limit. For a given value of x, the accidental STA times τ n are found to be: Using Equation (32) to determine the duration of the compression and adiabatic strokes leads to a many-particle QHE with zero friction, maximum efficiency and reduced cycle time with respect to the adiabatic case. Two fundamental bounds are found on the minimum time τ 1 for the existence of an STA under this driving protocol: where we used that Notice that Equation (33) is independent of N and λ . These bounds are sharp for x ∼ 1, τ n ∼ πn ω 1 and for x ∼ 0, τ n ∼ 1/(2ω 1 ). In Appendix A, it is shown that: as required for an STA protocol. Concerning many-particle effects, we point out that the efficiency at optimal power can be enhanced drastically for τ ≤ 1/ω 1 if λ = 0 (a similar conclusion holds for λ 1, typically λ < 0.1) where 1/ω 1 characterizes the transition between the sudden quench and adiabatic regime. Indeed, taking τ > 1/ω 1 , the nonadiabatic many-particle efficiency at optimal power becomes closer to the adiabatic efficiency at optimal power; see Figure 3.
Surprisingly, the nonadiabatic efficiency at optimal output power can beat the adiabatic value for τ 1 < τ < τ 1 , where τ 1 and τ 1 refer respectively to the first and the second intersection between the nonadiabatic and adiabatic efficiency; see Figure 3. This effective enhancement can be explained by the fact that the optimal value of ω 1 /ω 2 is smaller in the nonadiabatic case, for τ 1 < τ < τ 1 . However, the nonadiabatic efficiency is still bounded by the Otto limit, i.e., the adiabatic efficiency evaluated for the same frequency ratio. As already mentioned, in Figure 3, we find that for τ = τ 1 , as well as for τ = τ 1 , the efficiency at output power of the nonadiabatic protocol equals the value of an adiabatic protocol. For this to happen, τ 1 and τ 1 do not need to correspond to a proper STA with a nonadiabatic factor equal to one. Indeed, τ 1 (first intersection) does not. To find out the value of the minimum adiabatic shortcut time, one has to find the optimal frequency ratio for an adiabatic protocol (following the approach described in Section II C 2) and substitute this value into Equation (32) taking n = 1. We have verified numerically that τ 1 is the first accidental STA. For τ = τ 1 , the nonadiabatic efficiency of the scaled-up QHE matches the adiabatic value; yet, Q * > 1. FIG. 3. Efficiency of a many-particle quantum heat engine (QHE) run in finite-time at optimal power. Top panels (a-c) display the efficiency of a many-particle QHE at maximum power as a function of β h /β c for the accidental protocol (acc) with τ = 1.5/ω 1 , a sudden quench (sq) and adiabatic driving (ad). The interaction strength takes values λ = 0, 1, 2 from left to right and N = 500. A similar representation as a function of τ, in units of In both representations, a transition is observed from the sudden-quench to the adiabatic limits that is governed by the nonadiabatic coefficient computed via Equation (31).

B. Second Approach: Adiabatic Optimization and STA
In this section, we introduce an alternative approach to boost the performance of a many-particle QHE via STA. It exploits the possibility of engineering an STA that mimics an adiabatic expansion and compression stroke in an arbitrary fixed amount of time τ, between any two values of the frequency of the trap, ω 1 and ω 2 . This allows one to opt first for an optimization of the QHE parameters assuming the dynamics to be adiabatic. In particular, for given temperatures of the hot and cold reservoirs, the frequency ratio is chosen to maximize the adiabatic output power. This boils down to maximizing the total work per cycle, as the time of the compression and expansion strokes is fixed; see Section II C 2. The assumption of adiabatic dynamics is then effectively fulfilled by replacing the adiabats by STA protocols. This second approach is more versatile as the duration of the unitary strokes is not restricted to specific values and can be pre-scheduled. We shall consider three different techniques to engineer the required STA.

Reverse Engineering of the Scaling Dynamics
Via reverse engineering, one can formulate a general theory of STA thanks to the consistent description of the dynamics given by Equations (3) and (6). Indeed, if one assume that there exists a time τ, such that the scaling factor b(t) satisfies the boundary conditions [25,38]: it follows that: In particular, this is consistent with the boundary conditions that ω(t) must satisfy in the Ermakov Equation (3), which we rewrite as: and where ω(t) is the solution of the equation above with ω(0) = ω 0 and ω(τ) = ω 0 /b 2 ad . The conditions given by Equations (35a) and (35b) are satisfied if the scaling factor takes the following form: where a 3 = 10(b ad − 1)/τ 3 , a 4 = −15(b ad − 1)/τ 4 , and a 5 = 6(b ad − 1)/τ 5 ; see [25,38]. However, from Equation (37), it is not clear whether ω(t) is real for all t ∈ [0, τ]. There must be a condition on τ (for ω 0 and b ad fixed) to ensure the positivity of ω(t) 2 .
Equation (38) is one instance of STA constructed via reverse engineering of the scaling dynamics. The protocol defined by Equations (29a) and (29b) also satisfies the condition Equations (35a) and (35b); see Appendix A. Following this observation, one can construct many other driving protocols leading to accidental STA, such as: with β = (1 − b ad )/(2nπ) where n ≥ 1 is an integer and where the corresponding adiabatic shortcut times are τ n = 2nπt 0 . For the corresponding protocol obtained via Equation (37), ω(t) need not be real (as ω(t) 2 can be negative). For example, consider Therefore, this protocol can beat the lower bound for the protocol defined by Equations (29a) and (29b), for very small frequency ratio ω 1 ω 2 . It is worth mentioning that after exact numerical computation, we find that τ 1 ≤ 1/ω 1 for all values of x and is also less than the minimal time for STA found in Equation (32) (obtained with n = 1). Actually, if one considers τ and ω 1 (2) fixed, Equation (39) along with Equation (37) define a unique protocol (assuming ω(t) 2 > 0 for all t ∈ [0, τ]), with t 0 being the solution of the equation ω(τ) = ω 2 . In this case, we find another instance for which the minimum time for having an STA is smaller than the one required for the existence of an accidental STA described in Section III A; see Equation (32).

Counterdiabatic Driving
Demirplak and Rice introduced a formalism to "run a fast-motion video" of the adiabatic dynamics of a quantum system by assisting the evolution with counterdiabatic driving fields [61]. An equivalent approach was developed independently by Berry [62]. The technique is known as counterdiabatic driving (CD) and also referred to as transitionless quantum driving. It was extended to many-body systems in [40][41][42][43]. For a single-particle QHE, the use of CD was discussed in [9]. We next consider the performance of a many-particle QHE assisted by STA based on CD.
Let the instantaneous eigenvalue problem of the working-medium Hamiltonian with ω = ω(t) read: In its original formulation, CD ensures that the adiabatic approximation: becomes the exact solution of the time-dependent Schrödinger equation, when the original Hamiltonian is assisted by the additional term:Ĥ That is, given an initial state the dynamics generated by the CD HamiltonianĤ CD =Ĥ +Ĥ 1 is exactly given by the adiabatic approximation toĤ[ω(t)]. CD removes the requirement of slow driving, while the norm of the auxiliary CD term increases with decreasing running time [40,63]. Further, we note that the auxiliary term generally involves multiple-body interactions whose experimental realization is challenging [40,43]. However, in the presence of scaling dynamics, a remarkable simplification occurs; not onlyĤ 1 can be computed in a closed form without the use of the spectral properties of the Hamiltonian as an input, but it takes the form of a one-body counterdiabatic term [41]. For the working medium Equation (9), the driving Hamiltonian is: where p i is the momentum operator canonically conjugated to z i and {z, p} = zp + pz is the anticommutator. The exact solution of the time-dependent many-body Schrödinger equation is then given by: Under CD, the mean energy equals the adiabatic value Ĥ CD (t) = Ĥ (0) /b 2 ad (t) during the whole time-evolution for an arbitrary protocol ω(t), and the nonadiabatic factor reduces to unity, i.e., as shown in Figure 2. Note that the adiabatic evolution Equation (41) associated withĤ(ω) is nonadiabatic with respect to the instantaneous eigenbasis of the CD Hamiltonian (43).

Local Counterdiabatic Driving
The experimental implementation of the squeezing operator ∝ ∑ i {z i , p i } required as a counterdiabatic term discussed in the previous section, understood either as a auxiliary external control or as part of the intrinsic Hamiltonian describing the working medium, is generally a challenging task. As a result, it is natural to look for a unitarily-equivalent Hamiltonian in which the counterdiabatic term is mapped to a local potential. This has been achieved for a large family of many-body systems with scaleinvariant dynamics [41,42], as well as in the single-particle case [28]. The Hamiltonian (9) belongs to such a family. In this case, the driving of the working medium with a local counterdiabatic control (LCD) readŝ with: The time-evolving state is then given by while the corresponding mean energy is Assuming the initial state to be at equilibrium so that {z i , p i }(0) = 0, it follows that Under LCD, the nonadiabatic factor along the process is given by and reduces explicitly to unity at the beginning and end of the LCD protocol, providedω(0) =ω(τ) = 0. We mention that this equation can also be derived using our general equation (6) for Q * with the substitutions b → b ad and ω(t) → Ω(t). This is illustrated in Figure 2 for the driving protocol ω(t) varying from ω 1 -ω 2 and satisfying these boundary conditions together witḧ ω(0) =ω(τ) = 0. The explicit expression for ω(t) resembles that for the scaling factor in Equation (38) and reads where α 3 = 10(ω 2 − ω 1 )/τ 3 , α 4 = −15(ω 2 − ω 1 )/τ 4 and α 5 = 6(ω 2 − ω 1 )/τ 5 . The nonadiabatic nature of the STA designed by LCD becomes apparent in Figure 2a, where transient excitations are generated during the protocol and canceled out upon its completion. At variance with accidental STA, LCD protocols can be engineered for arbitrary small values of τ provided that the modulation of the trap frequency Ω(t) can be implemented; see Figure 2b.

IV. DISCUSSION
Along the compression and expansion strokes of a quantum Otto cycle, the nonadiabatic dynamics of a many-particle working medium can be efficiently characterized by the nonadiabatic factor Q * (t) whenever the dynamics is scale-invariant, as is the case for the family of models described by Hamiltonian (1). The adiabatic performance is matched whenever Q * (τ) = 1 at the end of the compression and expansion strokes governed by unitary dynamics. We have seen that the condition Q * (τ) = 1 can be fulfilled without resorting to adiabatic dynamics. In particular, a variety of STA can be exploited as an alternative. This possibility has prompted us to introduce two alternative schemes to boost the performance of a many-particle QHE. In both approaches, the fixed resources are given by the temperature of the hot and cold reservoirs, together with the particle number and the inter-particle interactions. One can opt for the optimization of the output power of the QHE run in finite-time τ as a function of the frequency ratio of ω 1 /ω 2 . Under the resulting condition for maximum finite-time output power, varying the value of τ for a given specific functional form of ω(t) can lead to accidental STA protocols satisfying Q * (τ) = 1, for a set of discrete values τ = τ n (n = 1, 2, 3, . . . ), as discussed in Section III A. As an alternative, one can directly opt for maximizing the output power of the QHE assuming zero friction (Q * (τ) = 1) and then use STA in many-particle systems to consistently match the adiabatic performance and reach the maximum efficiency of the cycle. Such STA protocols can be engineered by a variety of techniques, including reverse engineering of the scale-invariant dynamics (Section III B 1), counterdiabatic driving (Section III B 2) and local counterdiabatic driving (Section III B 3). Both approaches lead to the operation of the many-particle QHE at maximum efficiency and optimal output power.

V. CONCLUSIONS
Finite-time thermodynamics aims at optimizing the nonadiabatic performance of thermal machines, required for any realistic application. Conditions for optimal performance generally depend on the specific characteristics of the working medium, such as the number of particles and the interaction strength. At the quantum level, the required optimization involves tailoring thermal and quantum fluctuations. The maximum efficiency can be reached in the limit of slow driving, at the expense of a vanishing output power.
We have proposed to scale up quantum heat engines by using a many-particle working medium in combination with shortcuts to adiabaticity, to simultaneously amplify the output power and match the adiabatic performance in a finite cycle time. The resulting heat engine operates with zero friction at the maximum efficiency and exhibits an enhanced output power as a result of the many-particle nature of the working medium and the reduced cycle time. The maximum efficiency is then set by the adiabatic Otto upper bound shared by both single-and many-particle heat engines. This scheme can be implemented realizing a quantum Otto cycle with an interacting trapped Bose or Fermi gas with a suitable modulation of the harmonic trap. Our results provide a new avenue for engineering nonadiabatic friction-free scalable thermal machines that should find applications at the interface of quantum thermodynamics and energy science.

ACKNOWLEDGMENTS
Funding support from UMass Boston (Project P20150000029279) and European Science Foundation (POLATOM-5052) is further acknowledged.

APPENDIX A. NONADIABATICITY OF THE ACCIDENTAL PROTOCOL
Consider the accidental STA protocol defined by Equations (29a) and (29b), with the scaling factor b(t) and nonadiabatic factor Q * AB(CD) discussed in Section III A. In this Appendix, we compute these quantities for the first step of the Otto cycle A → B (compression). One can check that the nonadiabatic factor is symmetric Q * AB = Q * CD and that the scaling factor verifies the consistency conditions for a STA for the expansion C → D, as well as the compression. The scaling factor b(t) is given by the following formula [17,64]: b(t) = G 1 (t) 2 + ω(0) 2 G 2 (t) 2 , for t > 0. (A1) Using Equation (A2), we obtain: , for t = τ, leading to Formula (31).