Anti-Zeno quantum advantage in fast-driven heat machines

Developing quantum machines which can outperform their classical counterparts, thereby achieving quantum supremacy or quantum advantage, is a major aim of the current research on quantum thermodynamics and quantum technologies. Here we show that a fast-modulated cyclic quantum heat machine operating in the non-Markovian regime can lead to significant heat-current and power boosts induced by the anti-Zeno effect. Such boosts signify a quantum advantage over almost all heat-machines proposed thus far that operate in the conventional Markovian regime, where the quantumness of the system-bath interaction plays no role. The present effect owes its origin to the time-energy uncertainty relation in quantum mechanics, which may result in enhanced system-bath energy exchange for modulation periods shorter than the bath correlation-time.


Introduction
The non-equilibrium thermodynamic description of heat machines consisting of quantum systems coupled to heat baths is almost exclusively based on the Markovian approximation [1,2]. This approximation allows for monotonic convergence of the system-state to thermal equilibrium with its environment (bath) and yields a universal bound on entropy change (production) in the system [3]. Yet, the Markovian approximation is not required for the derivation of the Carnot bound on the efficiency of a cyclic two-bath heat engine (HE): this bound follows from the second law of thermodynamics, under the condition of zero entropy change over a cycle by the working fluid (WF), in both classical and quantum scenarios. In general, the question whether non-Markovianity is an asset remains open, although several works have ventured into the non-Markovian domain [4][5][6][7][8][9]. By contrast, it has been suggested that quantum resources, such as a bath consisting of coherently superposed atoms [10], or a squeezed thermal bath [11][12][13], may raise the efficiency bound of the machine. The mechanisms that can cause such a raise include either a conversion of atomic coherence and entanglement in the bath into WF heatup [10,14,15], or the ability of a squeezed bath to exchange ergotropy [11][12][13]16] (alias non-passivity or work-capacity [17][18][19]) with the WF, which is incompatible with a standard HE. However, neither of these mechanisms is exclusively quantum; both may have classical counterparts [20]. Likewise, quantum coherent or squeezed driving of the system acting as a WF or a piston [21] may boost the power output of the machine depending on the ergotropy of the system state, but not on its non-classicality [13].
Finding quantum advantages in machine performance * Electronic address: mukherjeev@iiserbpr.ac.in † Electronic address: kofmana@gmail.com relative to their classical counterparts has been one of the major aims of research in the field of quantum technology in general [22][23][24], and particularly in thermodynamics of quantum systems [25]. Overall, the foregoing research leads to the conclusion that conventional thermodynamic description of cyclic machines based on a (twolevel, multilevel or harmonic oscillator) quantum system in arbitrary two-bath settings may not be the arena for a distinct quantum advantage in machine performance [20]. An exception should be made for multiple identical machines that exhibit collective, quantum-entangled features [26,27]).
Here we show that quantum advantage is in fact achievable in a quantum heat machine (QHM), whether a heat engine or a refrigerator, whose energy-level gap is modulated faster than what is allowed by the Markov approximation. To this end, we invoke methods of quantum system-control via frequent coherent (e.g. phase-flipping or level-modulating) operations [28,29] as well as their incoherent counterparts (e.g. projective measurements or noise-induced dephasing) [30][31][32][33][34]. Such control has previously been shown, both theoretically [30,31,35,36] and experimentally [34,37], to yield non-Markovian dynamics that conforms to one of two universal paradigms: i) quantum Zeno dynamics (QZD) whereby the bath effects on the system are drastically suppressed or slowed down; ii) anti-Zeno dynamics (AZD) that implies the opposite, i.e., enhancement or speedup of the system-bath energy exchange [30,31,38]. It has been previously shown that QZD leads to the heating of both the system and the bath at the expense of the system-bath correlation energy [39], whereas AZD may lead to alternating cooling or heating of the system at the expense of the bath or vice-versa [30,31]. In our present analysis of cyclic heat machines based on quantum systems, we show that analogous effects can drastically modify the power output, without affecting their Carnot efficiency bound. AZD is shown to bring about a drastic power boost, thereby manifesting genuine quantum advantage, as it stems from the timeenergy uncertainty relation of quantum mechanics.

Cold Bath
Hot Bath Figure 1: A two-level system thermal machine: Schematic setup showing a two-level system with periodically modulated level distance ω(t) as the working fluid (WF) in a thermal machine wherein the WF is simultaneously coupled to hot and cold baths with non-overlapping spectra. Possible realizations include a driven-atom WF coupled to filtered heat baths in a cavity or a driven impurity coupled to spectrally distinct phonon baths in a solid structure [40,41].

Results
Model. We consider a quantum system S that plays the role of a working fluid (WF) in a quantum thermal machine, wherein it is simultaneously coupled to cold and hot thermal baths. The system is periodically driven or perturbed with time period τ S = 2π/∆ S by the timedependent HamiltonianĤ S (t): In order to have frictionless dynamics at all times, we chooseĤ S (t) to be diagonal in the energy basis of S, such that.
The system interacts simultaneously with the independent cold (c) and hot (h) baths viâ where the bath operatorsB c andB h commute: B c ,B h = 0, andŜ is a system operator. For example, for a two level system,Ŝ =σ x , whileŜ =X for a harmonic oscillator, in standard notations. We do not invoke the rotating wave approximation in the system-bath interaction Hamiltonian Eq. (3). As in the minimal continuous quantum heat machine [40], or its multilevel extensions [42], we require the two baths to have non-overlapping spectra, e.g., super-Ohmic spectra with distinct upper cut-off frequencies (see Fig. 1). This requirement allows S to effectively couple intermittently to one or the other bath during the modulation period τ S , without changing the interaction Hamiltonian to either bath. From Markovian to non-Markovian dynamics. In what follows we assume weak system-bath coupling, consistent with the Born (but not necessarily the Markov) approximation. Our goal is to examine the dynamics as we transit from Markovian to non-Markovian time scales, and the ensuing change of the QHM performance as the period duration τ S is decreased. To this end, we have adopted the methodology previously derived in Refs. [28,29,43,44], to account for the periodicity ofĤ S (t), by resorting to a Floquet expansion of the Liouville operator in the harmonics of ∆ S = 2π/τ S [40,45,46]. As explained below, we focus on system-bath coupling durations τ C = nτ S of the order of a few modulation periods, where n > 1 denotes the number of periods. The time-scales of importance are the modulation time period τ S , the system-bath coupling duration τ C , the bath correlation-time τ B and the thermalization time τ th ∼ γ −1 0 , where γ 0 is the system-bath coupling strength. We consider n 1 such that τ C τ S , (ω + q∆ S ) −1 , where ω denotes the transition frequencies of the system S, and q is an integer (see Methods "Floquet Analysis of the non-Markovian Master Equation"). This allows us to implement the secular approximation, thereby averaging over the fast-rotating terms in the dynamics. In the limit of slow modulation, i.e, τ S τ B , we have τ C τ B , which allows us to perform the Born, Markov and secular approximations, and eventually arrive at a time-independent Markovian master equation for τ C τ S , ω −1 , τ B (see Methods "Floquet Analysis of the non-Markovian Master Equation").
On the other hand, in the regime of fast modulation τ S τ B , the Markov approximation becomes inapplicable for coupling durations τ C = nτ S τ B . This gives rise to the fast-modulation form of the master equation (see Methods. "Floquet Analysis of the non-Markovian Master Equation" and "Non-Markovian dynamics of a driven two-level system in a dissipative bath"): For simplicity, unless otherwise stated, we consider = k B = 1. Here, for any modulation period τ S , the generalized Liouville operators L j of the two baths act additively on the reduced density matrix ρ S (t) of S, generated by the ω-spectral components of the Lindblad dissipators D j,ω (see below) for the j = c, h bath acting on ρ S (t). For a two-level system, or an oscillator, D does not depend on ω [1]. For ρ S (t) that is diagonal in the energy basis, which we consider below, the dynamics is dictated by the coefficients I j (ω, t) ≡ Re Ĩ j (ω, t) in Eq. (4), which express the convolution of the j-th bath spectral response function G j (ν) that has spectral width ∼ Γ B ∼ 1/τ B , with the sinc function, imposed by the time-energy uncertainty relation for finite times (see Methods "Non-Markovian dynamics of a driven two-level system in a dissipative bath").
Our main contention is that overlap between the sinc function and G j (ν) at t ∼ τ C τ B may lead to the anti-Zeno effect, i.e., to remarkable enhancement in the convolution I j (ω, t), and, correspondingly, in the heat currents and power. One can stay in this regime of enhanced performance over many cycles, by running the QHM in the following two-stroke non-Markovian cycles: i) Stroke 1: we run the QHM by keeping the WF (system) and the baths coupled over n modulation periods, from time t = 0 to t = nτ S = τ C τ B (n 1, τ S τ B ). The n modulation periods of the WF are equivalent to n cycles of continuous heat machines studied earlier, which have been shown to exploit spectral separation of the hot and cold baths for the extraction of work [40,46], or refrigeration [19,47], in the Markovian regime (see Eq. (8)). By contrast, in the non-Markovian domain a modulation period is not a cycle, since the time-dependent heat currents and the WF state are not necessarily reset to their initial values at the beginning of each modulation period (see below). ii) Stroke 2: In order to reset the WF state and the heat currents to their initial (t = 0) values in the non-Markovian regime, we have to add another stroke: At t = nτ S = τ C , we decouple the WF from the hot and cold baths. One needs to keep the WF and the thermal baths uncoupled (non-interacting) for a time-intervalt τ B , so as to eliminate all the transient memory effects [38]. After this decoupling period, we recouple the WF to the hot and cold thermal baths and continue to drive the WF with the periodically modulated Hamiltonian Eq. (1). Thus the setup is initialized after time τ C +t, provided we choose n to be such that ρ S (τ C +t) = ρ S (0), so as to close the steady-state cycle after n modulation periods, with the WF returning to its state at start of the cycle (see Fig. 2 and Sec. "A minimal quantum thermal machine beyond Markovianity"). The QHM may then run indefinitely in the non-Markovian cyclic regime. By contrast, in the limit of long WF-baths coupling duration τ C = nτ S τ B , the sinc functions take the form of delta functions, and therefore, as expected, the integral Eq. (4) reduces to the standard form obtained  The WF is first connected to the hot and cold baths, whose quasi-Lorentzian spectral functions are given by Eq. (54), at a negative time −tin (tin τ th ), under the initial condition p1(t = −tin) = 0.6, and reaches the steady-state value p1,ss at t + tin τ th . The WF is decoupled from the hot and cold thermal baths at a time −t −τB < 0 after reaching the steady-state, and then recoupled again to the two baths at time t = 0, such that the WF is non-interacting with the hot and cold thermal baths for the time interval −t ≤ t < 0, shown by the red break-line. The quantum heat machine is operated in the anti-Zeno dynamics (AZD) regime for t ≥ 0, wherein it is decoupled from and recoupled to the thermal baths after every AZD cycle, for coupling time duration τC = nτS. The probability p1 remains unchanged at the steadystate value, even after multiple AZD cycles. Inset: Same as the main plot, zoomed in for three consecutive AZD cycles. The WF is non-interacting with the thermal baths for time intervalst τB between two consecutive AZD cycles, shown by the red break lines. Here (see Eqs. in the Markovian regime, given by A minimal quantum thermal machine beyond Markovianity. Here we consider as the QHM a two-level system (TLS) WF with states |0 and |1 , interacting with a hot and a cold thermal bath, described by the Hamiltonian The Pauli matricesσ j (j = x, y, z) act on the TLS, the operatorB c (B h ) acts on the cold (hot) bath, andĤ B denotes the bath Hamiltonian. The resonance frequency ω(t) of the TLS is sinusoidally modulated by the periodiccontrol Hamiltonian where the relative modulation amplitude is small: 0 < λ 1. The periodic modulation Eq. (7) gives rise to Floquet sidebands (denoted by the index q = 0, ±1, ±2, . . .) with frequencies ω q = (ω 0 + q∆ S ) and weights P q , which diminish rapidly with increasing |q| for small λ (see Methods "Non-Markovian dynamics of a driven two-level system in a dissipative bath") [29,40,44].
(c) (c)   A crucial condition of our treatment is the choice of spectral separation of the hot and cold baths, such that the positive sidebands (q > 0) only couple to the hot bath and the negative sidebands (with q < 0) sidebands only couple to the cold bath. This requirement is satisfied, for example, by the following bath spectral functions: which ensures that for small λ, only the q = 1 harmonic exchanges energy with the hot bath at frequencies ±ω 1 = ± (ω 0 + ∆ S ), while the q = −1 harmonic does the same with the cold bath at frequencies We neglect the contribution of the higher order sidebands (|q| > 1) for 0 < λ 1, for which P q → 0 [19,29,40,44,46]. Further, we impose the Kubo-Martin-Schwinger (KMS) detailed-balance condition where β j = 1/T j .
For simplicity, in what follows, G h (ω) and G c (ω) are assumed to be mutually symmetric around ω 0 , i.e. they satisfy where α is a real positive number and 0 ≤ ν < ω 0 (see Methods "Steady states in the anti-Zeno dynamics (AZD) regime").
The WF is first coupled to the thermal baths at an initial time −t in (t in τ th > 0). Irrespective of the value of τ S , at large times t + t in τ th , and under the condition of weak WF-baths coupling, one can arrive at a time-independent non-equilibrium steady state ρ S → ρ ss in the energy-diagonal form (see Methods "Steady states in the anti-Zeno dynamics (AZD) regime"): One can then decouple the WF and the baths, such that they are non-interacting for a time interval exceeding τ B so as to eliminate all memory effects, then recouple them again at t = 0, keeping ρ S = ρ ss , and run the QHM in a cycle (as described in the Section "From Markovian to non-Markovian dynamics").
In general, owing to the finite widths (∼ 1/τ C ) of I h,c (ω q , t) in the frequency domain for short coupling times (τ C τ B ), the WF would be driven away from ρ ss , as follows from Eq. (4), causing ρ S (t) to evolve with time within the time interval 0 < t ≤ τ C . However, in order to generate a cyclic QHM operating in the steady-state, we focus on cycles consisting of n modulation periods that satisfy so that The above conditions Eq. (12) and (13), along with the KMS condition Eq. (9), imply that Equation (14), in turn, guarantees that Eq. (11) yields the steady state even at short times, and thus eliminates any time dependence in ρ S (see Fig. 2). For a QHM operating in the steady state, remains zero even during de-coupling from, and recoupling with the hot and cold baths. This ensures that the system remains in its steady state ρ ss throughout the cycle.
Equations (12) -(14) can be easily satisfied for experimentally achievable parameters; eg., ∆ S ∼ kHz, and n = 10 would imply T c ∆ S /2πnk B ∼ 10 −9 K. The number n = 10 was chosen to be around the minimal number n that allows for the validity of the secular approximation q = q in (23) and hence for a simplification (24) in the master equation. This number should be made as low as possible, since by decreasing n we decrease the cycle duration τ C and hence increase the power boost, as explained above. Since this power boost is then maximized without changing the efficiency, as noted above, the performance is optimized for the chosen n.
From the First Law of thermodynamics, the QHM output powerẆ (t) is given in terms of the hot and cold heat currents J h (t) and J c (t), respectively, by [19] The possible operational regimes of the heat machine, i.e., its being a heat-engine or a refrigerator [19,40], are determined by the signs of the WF-baths coupling duration-averaged J h , J c and W . One can calculate the steady-state efficiency η, average power outputẆ and average heat currents as a function of the modulation speed ∆ S , searching for the extrema of the functions in Eq. (16).
The heat currents J c and J h , flowing out of the cold and hot baths, respectively, are obtained consistently with the Second Law [19,40] in the form where we have used P ±1 = λ 2 /4. In order to study the steady-state QHM performances for different modulation frequencies, we consider the example of two non-overlapping spectral response functions of the two baths displaced by δ with respect to ω q , i.e., G h (ν) (G c (ν)) characterized by a quasi-Lorentzian peak of width Γ B , with the peak at Quasi-Lorentzian bath spectral functions"). Alternatively, we also consider the example of two non-overlapping super-Ohmic spectral response functions G h (ν) and G c (ν) of the two baths, with their origins shifted from ν = 0 by ν h = ω 0 + ∆ S − δ and ν c = ω 0 −∆ S +δ respectively, for 0 < δ ∆ S , ω 0 , ω 0 −∆ S (see Methods "Super-Ohmic bath spectral functions"). The dependence of ν h,c on ∆ S amounts to considering baths with different spectral functions for different modulation frequencies, and ensures that any enhancement in heat currents and power under fast driving results from the broadening (rather than the shift) of the sinc functions, which are centered at ω 0 ± ∆ S .
We plot quasi-Lorentzian bath spectral functions and the sinc functions in Figs. 3a and 3b, and the corresponding time averaged heat currents and power (see Eq. (17)) for the heat engine regime in Fig. 3c. We do the same for super-Ohmic bath spectral functions in Figs. 4a, 4b and 4c. The corresponding heat currents and powers for the refrigerator regimes are shown in Figs. 5a and 5b. The Markovian approximation: sinc(x) ∝ δ(x) in Eq. (4) reproduces the correct heat currents and power only in the limit of slow modulation (τ C τ B ). By contrast, the Markovian approximation reproduces the exact efficiency for both slow and fast modulation rates (see Fig.  6a). Thus, although the efficiency grows as τ S decreases, it is still limited by the Carnot bound.
Anti-Zeno dynamics. The performance of the QHM depends crucially on the relative width of the spectral function and the sinc functions. A slow modulation (τ S τ B ) results in sinc functions which are non-zero only over a narrow frequency range, wherein G j (ν) can be assumed to be approximately constant, which leads to time-independent I j (ω q ) and Markovian dynamics. On the other hand, fast modulation (τ C τ B ) is associated with broad sinc functions, for which G j (ν) is variable over the frequency range ∼ τ −1 C for which the sinc functions are non-zero (see Figs. 3a, 3b, 4a and 4b). This regime is a consequence of the time-energy uncertainty relation of quantum mechanics, and is associated with the anti-Zeno effect [30,31]. This effect results in dynamically enhanced system-bath energy exchange, which we dub anti-Zeno dynamics (AZD). Remarkably, in a QHM, appropriate choices of I h,c (ω q , t) may yield a power and heat-currents boost whenever the sinc functions have sufficient overlap with G h,c (ν) (see Figs. 3c and 4c).
Importantly, we find that spectral functions peaked at frequencies sufficiently detuned from ω 0 ± ∆ S (i.e., δ > Γ B ) may increase the overlap with the sinc functions appreciably under fast modulation in the anti-Zeno regime, for τ −1 C , δ Γ B , thus resulting in substantial output power boost. This regime indicates that finite spectral width of the sinc functions may endow a HE with significant quantum advantage, arising from the timeenergy uncertainty relation, which is absent in the classical regime, be it Markovian or non-Markovian. In the numerical examples shown here, the quantum advantage in the HEs powered by baths with quasi-Lorentzian (super-Ohmic) spectral functions can increase the power by a factor larger than two (seven) (see Figs. 3c and 4c), for the same efficiency (see Fig. 6a and Methods "Efficiency and coefficient of performance").
Quantum refrigeration: AZD can lead to quantum advantage in the refrigerator regime as well, for modulation rates beyond the quantum speed limit [42,48] (see App. A), by enhancing the heat current J c , thus resulting in faster cooling of the cold bath. As for HE, numerical analysis shows that quasi-Lorentzian, as well as super-Ohmic bath spectral functions can lead to significant quantum advantage in the AZD regime (see Fig.  5). On the other hand, as for the efficiency in case of the HE, the coefficient of performance is not significantly affected by the broadening of the sinc function, and on average remains identical to that obtained under slow modulation in the Markovian regime (see Fig. 6b and Methods "Efficiency and coefficient of performance").

Discussion
We have explored the hitherto uncharted domain of quantum heat engines (QHEs) and refrigerator (QRs) based on quantum working fluids (WFs) intermittently coupled and decoupled from heat baths operating on non-Markovian time scales. We have shown that for driving (control) faster than the correlation (memory) time of the bath, one may achieve dramatic output power boost in the anti-Zeno dynamics (AZD) regime.
Let us revisit our findings, using as a benchmark the Markovian regime under periodic driving: In the latter regime, detailed balance of transition rates between the WF levels, as well as the periodic driving (modulation) rate, determine, according to the First and Second Laws of thermodynamics, the heat currents between the (hot and cold) baths, and thereby the power produced or consumed. In our present treatment, the Markovian   (55)). The enhanced overlap resulting from fast modulation (large ∆S) enhances the heat currents Jc to up to a factor larger than 2 for (a) and larger than 9 for (b) in the refrigerator regime (shown by dotted double-arrowed lines), signifying quantum advantage. The green dotted line corresponds to zero power and currents. Here λ = 0.2, ω0 = 20, γ0 = 1, β h = 0.001, βc = 0.002. regime is recovered under slow modulation, such that the WF-baths coupling duration τ C exceeds the bath correlation time τ B . Then, the Markovian approximation is adequate for studying the operation of the QHE or the QR. By contrast, under fast modulations, such that τ C = nτ S τ B , the working fluid interacts with the baths over a broad frequency range of the order of ∼ τ −1 C , according to the time-energy uncertainty relation in quantum mechanics. The frequency-width over which system- bath energy exchange takes place can lead to anti-Zeno dynamics (AZD). The resultant quantum advantage is then especially pronounced for bath spectral functions that are appreciably shifted by δ > Γ B ∼ τ −1 B , from the centers of the sinc functions that govern the system-bath energy exchange rates.
We have explicitly restricted the results to mutually symmetric bath spectral functions (e.g., the experimentally common Lorentzian or Gaussian spectra), in order to ensure time-independent steady-states of the WF. Yet this requirement is not essential, since the WF steadystate may be time-dependent as long as it is periodic so as to allow for cyclic operation. The AZD [28][29][30][31][32][33][34][35] can arise for any bath spectra of finite width ∼ 1/τ B , as long as nτ S τ B . One can therefore operate a thermal machine provided stroke 1 of the cycle is in the AZD regime and achieve a quantum advantage without additional restrictions on the bath spectral functions (see Methods "Thermal machines with arbitrary (asymmetric) spectral functions").
The QHM discussed here is driven by external modulation. As previously shown both theoretically [28][29][30][31][32]35] and experimentally [34,37], periodic perturbations of the TLS state can increase its relaxation rate in the non-Markovian anti-Zeno regime. The reason for the power boost is that at the non-Markovian stage of the evolution which occurs on short time scales, the sinc factors in the convolutions with G(ω), as in (32), are sufficiently broad so as to modify the convolutions and hence the relaxation rates in (34) in comparison with the Markovian case, where these sinc functions are spectrally narrow enough to be approximated by delta-functions. Under the conditions chosen in the paper, this modification leads to an increase in the TLS relaxation rates and hence to a power boost. This boost is of quantum nature, since the broadening of the sinc factors is due to the quantum time-energy uncertainty relation that may lead to the violation of energy conservation at short times. The quantum mechanical time-energy uncertainty relation employed here reflects the fact that the Scrödinger equation for a two-level system coupled to a bath renders the energy transfer probability from the two-level system to the bath and back oscillatory in time. Such oscillation leads at short times (comparable to the required cycle period) to sinc-like deviation from delta-function energy conservation. Classical description of analogous processes, even beyond the Markovian approximation, does not involve discrete energy levels and hence no oscillations of the system-bath transfer rate that deviates from energy conservation. Thus, the effects discussed here are inherently quantum mechanical.
The non-Markovian effect in the present context is quantified by the spectral widths of the sinc functions compared to the bath-response G(ω) spectral width 1/τ B . If the cycle duration is kept fixed, then the non-Markovian effect scales with the spectral width of G(ω). Hence, super-Ohmic bath spectra with their salient cutoff provide realistic examples of the non-Markovian effects described here. Such bath spectra should be contrasted with the flatter and broader Ohmic spectra. Yet, non-Markovian dynamics does not necessarily imply a quantum advantage, as discussed in App. B.
The predicted power boost relies on transient dynamics: the heat fluxes change with time within t = τ C in the non-Markovian AZD regime, even when the WF state hardly changes during that time interval. Yet it is essential that we incorporate this transient dynamics within steady-state cycles by decoupling the WF from the baths, allowing the bath-correlations to vanish within τ B and then recoupling the WF again to the baths when they have all resumed their initial states. These cycles can be repeated without restriction, thereby allowing us to operate the QHM with quantum enhanced performance even for long times, despite the reliance on transient dynamics within the stroke 1 of each cycle.
The quantum advantage of AZD, at zero energetic cost (see App. C), manifests itself in the form of higher output power, for the same efficiency, in the QHE regime (∆ S < ∆ qsl ), as compared to that obtained under Markovian dynamics in the limit of large τ C , all other parameters remaining the same. Alternatively, in the QR regime (∆ S > ∆ qsl ), AZD may lead to quantum advantage over Markovian dynamics in the form of higher heat current J c , or, equivalently, higher cooling rate of the cold bath, for the same coefficient of performance. The latter effect leads to the enticing possibility of quantum-enhanced speed-up of the cooling rate of systems as we approach the absolute zero, and raises questions regarding the validity of the Third Law of Thermodynamics in the quantum non-Markovian regime, if we expect the vanishing of the cooling rate at zero temperature as a manifestation of the Third Law [47,49,50].
The QHE power boost in the anti-Zeno regime results from a corresponding increase in the rates of heatexchange and entropy production, arising from the TLS relaxation by both baths. This is the reason that the efficiency, i.e., the ratio of the work output to the heat input, is unchanged, i.e. is the same as in the standard Markovian regime. Yet, all parameters being equal, the QHM rate of operation (as measured by the power output) speeds up in the anti-Zeno regime, which constitutes a practical quantum advantage.
One can extend the analysis discussed here to Otto cycles [42,51]: Fast periodic modulation during the nonunitary strokes of an Otto cycle can speed up the thermalization through AZD, thereby allowing quantum enhanced performance. Interestingly, fast modulation in the Otto cycle can yield enhanced power or refrigeration rate, even in the Markovian regime [52]. Finally, in the regime of ultrafast modulation with τ −1 C Γ B , δ, quantum Zeno dynamics sets in, leading to vanishing heat currents and power, thus implying that such a regime is incompatible with thermal machine operation (see Fig. 7, App. B and Fig. 8). While Zeno dynamics has commonly been associated with measurements [30,[53][54][55], both the Zeno and the anti-Zeno effects occur under various frequent perturbations, such as phase flips and nonselective (unread) measurements. Generally speaking, to observe the discussed effects, it is sufficient to repeatedly perform cycles of coupling the system (here -the WF) with another system (here -a bath), then destroying or sharply changing the coherence between the two. This sharp change can be effected in different ways, e.g., by a measurement of the system (which can be readout or not) or, as in our case, by abruptly decoupling the WF and the baths, which gives rise to Zeno or anti-Zeno dynamics [28][29][30][31][32][33]. In contrast to previous studies, here the Zeno or anti-Zeno dynamics of the WF arises by an external periodic field, and therefore not around the frequency ω 0 of the unperturbed system, as in previous cases, but at multiple sideband frequencies ω 0 + q∆ S .
The Markovian approximation suffices to find the correct efficiency (for a QHE) or the coefficient of performance (for a QR), even in the non-Markovian regimes, for mutually symmetric bath spectral functions (see Eq. (10)). This guarantees that the efficiency always remains below the Carnot bound, even under fast modulations (see Fig. 6).
Our scenario is conceptually different from that in which work is produced by a QHE on an external quantum system and quantum effects arise from the interaction between the quantum WF and the external quantum system [56]. Such quantum effects are absent in our case, where work input in the QHE is provided by a classical field.
Experimental scenarios where the predicted AZD quantum advantage may be tested are diverse. Since non-Markovianity in general, and AZD in particular, require non-flat bath spectral functions, suitable candidates for the hot and cold baths are microwave cavities and waveguides in which dielectric gratings are embedded, with distinct cut-off and bandgap frequencies [19,57] and the WF is a qubit whose level distance is modulated by fields. The required qubit modulations are then compatible with MHz periodic driving of superconducting transmon qubits [58,59] or NV-center qubits in diamonds [60]. One can effectively decouple the WF from the thermal baths in stroke 2 by abruptly changing the resonance frequency of the two-level WF from ω 0 tõ ω, thus rendering the WF strongly off-resonant with the thermal baths, so that G h (ω) ≈ G c (ω) = 0), thereby precluding any energy flow between the baths and the WF. We can recouple the WF with the thermal baths by reverting this frequency back to ω 0 , and then modulating it periodically, so as to generate either Markovian or non-Markovian anti-Zeno dynamics, as discussed above [34].
The AZD regime was experimentally observed in NMR setups [34]. Micro/nano-scale heat machines have been experimentally realized, for a trapped calcium ion as the WF [61]; nano-mechanical oscillators WF powered by squeezed thermal bath [12]; atomic heat machines assisted by quantum coherence [62]; or a nuclear spin 1/2 as the WF in a quantum Otto cycle [63].
The novel effects and performance trends of QHE and QR in the non-Markovian time domain, particularly the anti-Zeno induced power boost, open new, dynamicallycontrolled pathways in the quest for genuine quantum features in heat machines, which has been a major motivation of quantum thermodynamics in recent years [10-13, 21, 25, 41, 64-68].

Methods
Floquet Analysis of the non-Markovian Master Equation. Let us consider the differential non-Markovian master equation for the system density operator ρ S (t) in the interaction picture [29]: Here ρ B = ρ Bc ⊗ ρ Bh , where ρ Bj is the density operator of bath j. In the derivation of Eq. (19) we have assumed that Tr[B j , ρ Bj ] = 0. We consider commuting bath operators B c (t),B h (t ) = 0, such that the two baths act additively in Eq. (19). Below we focus on only one of the baths and omit the labels c/h for simplicity. We then haveŜ † (t) =Ŝ(t) where q are integers and ω are transition frequencies of the system S. One can use Eq. (20) to write the first term on the r.h.s. of Eq. (19) as In the limit of times of interest, i.e., times larger than the period of driving τ S and the effective periods of the system, t τ S , (ω + q∆ S ) −1 , the terms with the fast oscillating factor before the integral in Eq. (21) become small and can be neglected, i.e., the secular approximation becomes applicable, such that which generally holds only for as long as (q − q) ∆ S is not close to (ω − ω) for any q, q , ω, ω . Condition (23) gives us where µ = t − s, and In the limit of slow modulation, such that t = nτ S τ B , one can perform the Markov approximation, thereby extending the upper limit of the integral in time in Eq. (24) to t → ∞, which finally results in the timeindependent Markovian form [1] T 1 ≈ −π ω,q≥0Ŝ † q,ωŜq,ω ρ S (t)G(ω, q).
On the other hand, in the limit of t ∼ nτ S τ B , the Markovian approximation becomes invalid, and one gets Progressing similarly as above, one can arrive at similar expressions for other terms in Eq. (19) as well.
Non-Markovian dynamics of a driven two-level system in a dissipative bath. The non-Markovian master equation followed by the TLS WF subjected to the Hamiltonian Eq. (6) is (see Eq. (19)) where we have removed the h, c indices for simplicity, and considered the dynamics due to a single bath. Here A ↓ ,Ā ↑ andM are the complex conjugates of A ↓ , A ↑ and M , respectively. The terms corresponding to σ ± ρ S (t)σ ± in Eq. (28) vanish for diagonal steady-state ρ S (t) → ρ ss (see Eq. (11)). Here we focus on times longer than several modulation periods, i.e., t = nτ S τ S , when the fast oscillatory terms corresponding to q = q vanish as well, such that We note that The imaginary part in Eq. (31) acts on terms of the form , which vanish at large times when the off-diagonal elements ρ S (t) approach zero for any initial state. On the other hand, the real part of Eq. (31) gives rise to terms of the form In the limit of slow modulation such that t ∼ nτ S τ B (n ∈ Z, n 1), the function sin ([ν ± (ω 0 + q∆ S )] t) / [ν ± (ω 0 + q∆ S )] assumes a delta function centred at ν = ± (ω 0 + q∆ S ), thus leading to the familiar Markovian form of master equation, with On the other hand, in the anti-Zeno regime of fast modulation: t ∼ nτ S τ B , I j (±ω q , t) is not given by Eq. (33), and one needs to consider the full form Eq. (32).
In particular, for a diagonal state ρ S (t) = p 1 (t)|1 1| + p 1 (t)|0 0|, the dynamics Eq. (28) -Eq. (32) leads us to the rate equationṡ In the Zeno regime of ultra-fast modulation, obtained in the limit of t ∼ nτ S τ B , the integral I(ω q , t) vanishes (see Fig. 8), thus leading to the Zeno effect of no dynamics.
Equations (19), (24) and (28) are of the type known as the differential master equation (DME). An alternative approach is based on the (less convenient) integrodifferential master equation (IME). The two equations are mathematically different and hence require different procedures for reducing them to the Markovian master equation (MME). However, the IME and the DME have the same validity conditions, i.e., generally similar accuracy. The IME and the DME follow from the exact expansions in the totally ordered and partially ordered cumulants, respectively, upon neglecting terms of order higher than 2 in the system-bath coupling, which determines their accuracy [1,35,69] The rates (34) can be negative when the modulation (or measurement) period is short enough to break the rotating wave approximation (RWA). Yet the probabilities p 0 (t), p 1 (t) are never negative, as detailed in Refs. [31,35] and concisely proven in the next section.
Non-Markovian master equation with nonnegative probabilities. The non-Markovian master equations (MEs) for an arbitrarily driven (controlled) two-level system (TLS) presented in (34) have been derived and discussed in Refs. [28][29][30][31][32] and experimentally verified in Refs. [34,37]. These MEs involve the time-dependent relaxation rates R 0 (t) and R 1 (t) which can take negative values, since the quantities I j , (32), are convolutions of a positive spectral response function G j (ν) with a sinc function which takes positive or negative values. As a result, the solutions of the MEs for the populations (probabilities) p 0 (t) and p 1 (t) of the TLS levels are not guaranteed to be non-negative, i.e., to satisfy 0 ≤ p k (t) ≤ 1 (k = 0, 1).
Below we show that the inequalities (35) hold, at least, up to second order in the system-bath coupling strength. This means that for a weak coupling, violations of (35) (if any) are negligibly small. First, we note that at sufficiently long times, t τ B , the MEs become Markovian and coincide with the Lindblad equation. In this case, the rates are constant and positive, R 0 , R 1 ≥ 0, as follows from (33). The inequalities (35) are now known to hold. Generally, the MEs are valid if the couplings of the TLS with the baths are sufficiently weak, so that Consider now the short times, t τ B , where the non-Markovian effects are important. Since p 0 (t) + p 1 (t) = 1, we rewrite where w(t) = p 1 (t) − p 0 (t) is the TLS population inversion. In terms of w(t), inequalities (35) are equivalent to which we now prove. The condition (36) implies that at times τ B , the relaxation can be approximated to first order in the relaxation rates. In this approximation, (34) yields where and From (32) and (34), one can check that From (39) we obtain yielding (38). The second inequality in (43) follows from the assumption |w(0)| ≤ 1 and the relation |J − (t)| ≤ J + (t), resulting from (40) and (42). Steady states in the anti-Zeno dynamics regime. Now we study the regimes which allow us to operate the setup with a time-independent steady state ρ ss even inside the AZD regime. We note that for t → ∞, I j (ω q , t) reduces to the time-independent form πG j (ω q ), thus leading us to the Eq. (11). On the other hand, for t ∼ nτ S τ B , I j (ω q , t) includes contributions from G j (ω q + ν), where |ν| 1/t = 1/ (nτ S ) .
Thermal machines with arbitrary (asymmetric) spectral functions. We consider where, as before, α > 0 andχ(ν) is an arbitrary real function of ν. We then have I h (ω 0 + ∆ S ) = αI c (ω 0 − ∆ S ) + χ(t), anḋ In this case, we get a time-dependent steady-state with Therefore, the rate of change of w(t) with time is given bẏ One can still operate the setup as a cyclic thermal machine for a time t ≤t, as long as χ(t) andχ(t) are small enough so as to ensurė whereẇ max is the maximum value attained by |ẇ(t)| in the time interval 0 ≤ t ≤t.

Data availability
All relevant data are available to any reader upon reasonable request.

Code availability
All relevant codes are available to any reader upon reasonable request.