Unravelling the non-classicality role in Gaussian heat engines

At the heart of quantum thermodynamics lies a fundamental question about what is genuine"quantum"in quantum heat engines and how to seek this quantumness, so that thermodynamical tasks could be performed more efficiently compared with classical protocols. Here, using the concept of $P$-representability, we define a function called classicality, which quantifies the degree of non-classicality of bosonic modes. This function allows us to explore the role of non-classicality in quantum heat engines and design optimal protocols for work extraction. For two specific cycles, a quantum Otto and a generalised one, we show that non-classicality is a fundamental resource for performing thermodynamic tasks more efficiently.

At the heart of quantum thermodynamics lies a fundamental question about what is genuine "quantum" in quantum heat engines and how to seek this quantumness, so that thermodynamical tasks could be performed more efficiently compared with classical protocols. Here, using the concept of P -representability, we define a function called classicality, which quantifies the degree of non-classicality of bosonic modes. This function allows us to explore the role of non-classicality in quantum heat engines and design optimal protocols for work extraction. For two specific cycles, a quantum Otto and a generalised one, we show that non-classicality is a fundamental resource for performing thermodynamic tasks more efficiently. The effort to understand the relationship between quantum mechanics and thermodynamics makes the field of quantum thermodynamics extremely wide. Due to its broad nature, questions underlining the interplay between quantum features, such as coherence and entanglement, have been thoroughly explored using different approaches [1][2][3][4]. One aspect broadly investigated is the production of work in quantum heat engines. These are composed of one or more quantum systems that operate between two different reservoirs with the sole aim of converting heat into work. Although the laws of quantum mechanics rule these engines, this does not necessarily imply a "quantum advantage" [5][6][7][8]. In fact, a remarkable similarity was observed with classical models [9][10][11][12][13][14][15][16][17][18], raising the question of how the intrinsic features of quantum mechanics, such as entanglement and coherence, could be used to enhance the performance of quantum heat engines. The use of non-thermal baths, i.e., engineered reservoirs characterised by their temperatures and additional parameters, has presented a comprehensive scenario to study the relation between quantum effects and thermodynamic efficiency [1,[19][20][21][22][23][24][25][26][27]. In particular, non-thermal baths were previously shown to be used as a resource to exceed the standard Carnot limit [22,28]. In the context of heat engines, squeezed thermal baths have played an essential role, since a proof-of-principle experiment based on a nanobeam heat engine has recently been reported [29]. Subsequently, in [30], the author developed a framework underlying the main properties of this reservoir and showed that a squeezed thermal bath could be recast as a generalised equilibrium reservoir. However, in many of these works, although a non-classical feature was present, its role was not fully investigated, in favour of a previously established concept of non-passivity [31,32]. But how can we quantify non-classicality in heat engines? Can we quantify non-passivity through measurable quantities? An answer may be given from the principles of quantum optics.
Radiation fields with states described by specific features that can only be understood by a quantum mechanical description are known as non-classical states [33][34][35]. The definition of non-classicality is not unique, but for bosonic fields, the most strict way to quantify it is through the concept of P -representability, which states that a given state ρ is said to be P -representable if it can be written as a convex mixture of coherent states with a proper probability distribution function P (α, α * ) over an ensemble of states, i.e., P (α, α * ) is non-negative and is less, or equally, singular than the delta distribution -these states are known as classical states. In contrast, non-classical states corresponds to those that cannot be written as in Eq. (1), because its quasi-probability distribution is negative or highly singular. This definition is particularly relevant in the context of quantum information systems based on bosonic modes since only non-P -representable states can generate entanglement when mixed with vacuum in a beam-splitter [36][37][38]. A prominent example of non-classical states is the set of squeezed states, in which the fluctuation associated with one quadrature component is below the vacuum state [39]. Early theoretical work in the 60s and 80s led to the conclusion that quantum fluctuations can be reduced below the shot noise in many forms of nonlinear optical interactions [40][41][42]. For example, squeezed states are produced in nonlinear processes called degenerate parametric down-conversion, where a "classical" electromagnetic field drives a nonlinear medium and pairs of correlated photons of the same frequency are generated. The non-classical effects of light can be revealed in different ways. In the present case, this manifestation occurs in terms of coherence (off-diagonal elements) since the squeezing operation induces non-diagonal elements in the energy eigenbasis.
We aim to present a discussion by combining tools from quantum optics and continuous variable systems. As a result, we explore the role of non-classicality in quantum thermodynamics by quantifying the degree of non-classicality in terms of the Prepresentability. It is worth stressing that non-passivity is a necessary but not sufficient condition to detect non-classicality, and this distinction is not clear for a single-bosonic mode. However, the converse holds, i.e., all non-classical states are non-passive. For two specific examples, i.e., a quantum Otto and a generalised cycle, we show that the non-classicality is a resource for quantum heat engines, meaning that a given thermodynamic task can be performed more efficiently when this engine operates in the non-classical regime. The first analysed cycle-the Otto cycle-agrees with the results of Refs. [22,30] where they employed a modulation in the frequency of the bosonic mode as a work parameter. Since we are interested in the trade-off between non-classicality and thermodynamic efficiency, the squeezing parameter is used as the work variable.
The paper is organised as follows. In Section II, we start by setting the notation and introduce some preliminaries concepts involving the Gaussian parametrisation and the definition of P -representability. Then, we define the classicality function and briefly recall the concept of passive states. In Section III, the thermodynamic approach is described, and in Section IV we use the defined concepts to analyse two different cycles. Finally, we conclude with an outlook in Section V.

A. Classicality
A multi-mode Gaussian state, with bosonic operators a = (a 1 , a † 1 , ..., a N , a † N ) T satisfying the canonical commutation relation [43,44], is entirely described by its two first statistical moments, the first-order, or displacement, and the second-order moments collected in the covariance matrix (CM). For the problem under consideration, we assume a single bosonic mode with zero first moment and CM given by The diagonal elements of V represent the mean photon number, or occupation, whereas the off-diagonal elements represent correlations between different eigenstates, and these will be referred to as coherence parameter. A system described by a Hamiltonian H, and in thermal equilibrium with a thermal environment at temperature T , is described by a thermal state where k B denotes the Boltzmann constant. This state is an example of a Gaussian state, with null first moment and covariance matrix given by where the thermal average number of photonsn th = (e ω/k B T − 1) −1 is given by the Bose-Einstein distribution. A squeezed thermal state ρ sth is a Gaussian state with zero displacement and CM equal to where now the occupation and coherence parameters are given byn sth =n th cosh 2r + sinh 2 r, with r being the squeezing parameter, and m sth = (n th + 1/2) sinh 2r, respectively. As mentioned before, without loss of generality, the direction in the phase space in which the squeezing occurs is set to zero. Therefore, a squeezed thermal state is completely parameterised by its temperature and the squeezing parameter By definition, a state is P -representable (or classical) if its covariance matrix satisfies the following relation [45] (see Appendix A for more details).
meaning explicitly from Eq. (3) thatn From the condition (8), it is inferred that not all states are classical and their classicality depends on the occupation and coherence parameter. Although the thermal state is always classical (m is always zero), a squeezed thermal state may not be depending on r. Increasing the squeezing parameter, a transition from the classical to the non-classical regime is observed at a critical value r ≡ r c . This threshold, r c , is reached whenn sth = m sth . Consequently, for a fixed temperature, the critical squeezing parameter is given by FIG. 1. (a) A sketch of the parametersnsth and msth as a function of the squeezing parameter. Initially, the system is prepared in a thermal state at temperature T and is classical (blue region). By applying a squeezing operation and modulating the squeezing parameter, the system goes from a classical to a non-classical state (red region). This transition occurs with respect to a critical squeezing parameter rc(T ). In (b) we plotted the classicality as a function of the squeezing parameter for three different temperatures. Observe that after rc the system becomes non-classical, thereby characterising a negative classicality function C < 0 Physically, r c corresponds exactly when one of the quadratures reaches the uncertainty bound, and the state cannot be represented as being merely a mixture of coherent states. In Fig. 1(a), for a fixed T , we plot the behaviour of n sth and m sth as a function of the squeezing parameter. The states between these two curves are classical (blue-shaded area). For any r > r c , the coherence parameter is larger than the mean photon number m sth > n sth , and therefore the state is non-classical (red-shaded area). Following the P -representability criteria (7)- (8), we can define the classicality function as the difference between the mean number of photons and the coherence parameter: C ≡n − |m|. This function not only makes the Eq .(8) operational but also gives a non-classicality quantifier. Hence, C < 0 indicates that the system is described by a non-classical state. For the squeezed thermal states, the classicality function is given by Observe that the classicality of a given state depends exclusively on the temperature and the squeezing parameter. Classical states with different temperatures will typically be characterised by different classicalities. As the temperature increases, the system becomes more classical, and a more considerable value of r is required to achieve its non-classical character (see Fig. 1(b)). As previously mentioned, this is related to the uncertainty in the quadratures of the bosonic mode -a non-classical state is attained when this uncertainty falls below the shot noise.
Note that coherent states also present non-diagonal terms in the covariance matrix, which could also be identified as a source of non-classicality. However, these states do not show any non-classical manifestation, and this is clear from the definition of the classicality function since, for coherence states, its classicality function is zero. Therefore, coherent states can be considered to be on the verge of being classical and non-classical states.
B. Passivity Consider a quantum system described by a Hamiltonian H and prepared in a state ρ. We say that ρ is passive if for all unitaries U that act on the system. States satisfying the above equation cannot have their energy reduced by deterministic unitary transformations. The necessary and sufficient conditions for a state to be passive require that the former commute with the Hamiltonian, i.e., [ρ, H] = 0, and for the eigenvalue decomposition ρ = d−1 n=0 p n |n n|, we find that n ≤ m implies p n ≥ p m for all n and m in {1, ..., d} [46]. Moreover, a state is k-passive if ρ ⊗k , with k ∈ N, is passive. If for all integers k ≥ 1, the state is k-passive, thus it is called completely passive.
A thermal state, given by a diagonal covariance matrix, is the only completely passive state. This means that it does not mind how many copies of a thermal state we have, unitary operations cannot lower its energy. Consequently, any resource state for a cyclic engine must be out of thermal equilibrium [47] and any m = 0 in Eq. (3) defines a non-passive state, but does not necessarily a non-classical state, as we see in Fig. 1(b). For example, a bosonic mode in contact with two reservoirs in thermal states at temperature T 1 = T 2 is passive and, to produce work, it is necessary that the reservoirs are out of equilibrium, T 1 = T 2 . Alternatively, by squeezing the mode, the degree of non-passivity will increase, and also its non-classsicality, which is measured by Eq. (10). However, the state is only non-passive and non-classical for r > r c . This distinction will be crucial during the description of the two thermodynamic protocols.

III. THERMODYNAMIC SETTING
We consider a quantum heat engine based on a single bosonic mode as a working substance. The bosonic mode operates between two reservoirs at different temperatures: a cold thermal bath at temperature T 1 and a hot squeezed thermal bath at T 2 > T 1 with a squeezing parameter r R . The bosonic mode undergoes two different cycles, with the work parameter being the squeezing. More precisely, both processes consist of reversible operations of squeezing and unsqueezing the mode through different protocols designed to explore the role of the classicality function (see Fig. 2(a)).

A. System-bath interaction
We address the problem by assuming a weak coupling regime between the single-mode and heat baths, where the dynamics is modelled adopting a time-independent Lindblad master equationρ = L(ρ). The single-mode, initially isolated and described by a Hamiltonian H = ω 0 a † a, is weakly coupled to a bosonic reservoir, H B = j ω j b † j b j , via a coupling constant k j in the rotating wave approximation according to the interaction The cold reservoir is prepared in a thermal state, and the hot reservoir is prepared in a squeezed thermal state, i.e., a thermal state at temperature T 2 to which the squeezing operator has been applied to each oscillator composing the bath. Under the Born-Markov approximation, the master equation for the single mode when interacting with the squeezed thermal bath is given by [39,48]: where γ is the coupling constant between the single mode and the bath. The parametersn sth =n th cosh 2r + sinh 2 r and m sth = (n th + 1/2) sinh 2r are the same as those in Eq. (6). The stationary solution of Eq. (12) results in the squeezed thermal state [30], where β sth = β cosh(2r) the generalised inverse temperature, and µ = tanh(2r) the chemical-like potential. The operator A has the following form and it is known as the second-order moment's asymmetry since it measures how asymmetric, or compressed, is the state in the phase-space picture. In this work, we will be interested in the regime in which the interaction time between the system and the reservoir goes to infinity, i.e., the asymptotic limit. However, the finite and non-equilibrium regime can also be explored by considering the complete dynamical evolution in Eq. (12).

B. Thermodynamic work
The next step in our description is to define work. Moving from the macroscopic description to the microscopic quantum realm, a clear picture of what is work and heat becomes blurred, as fluctuations and randomness are fundamentally unavoidable. Thus, carrying out these concepts into the quantum realm is not a straightforward task. However, developments in stochastic thermodynamics [49,50], open quantum systems [11,51], and quantum information theory [52,53] have led to significant progress in these efforts.
In a complete analogy with classical thermodynamics, we define work as a controllable (coherent) energy exchange related to a parameter of the system and heat as an (incoherent) flux of energy that cannot be given or subtracted in a controlled manner (or be helpful in some process). As we are interested in a single bosonic mode prepared in a squeezed thermal state with mean energy equal to E(r,n th ) = ω n th + 1 2 cosh 2r , one can differentiate Eq. (14) to obtain two different contributions to the mean energy where we identify the first term as work (an external control parameter) and the second term as heat δW = −2 ω n th + 1 2 sinh 2r dr and δQ = ω cosh 2r dn th .
The minus signal in the work expression is adopted because, during the preparation of the squeezed state, external energy is required.
The connection between the extracted work and non-classicallity (coherence) is made using the relative entropy of coherence [54], which reveals that the amount of coherence present on the Fock basis is proportional to m and inversely proportional to T C(ρ sth ) ≈ β s m .
where C(ρ) = S(ρ diag ) − S(ρ) is the relative entropy of coherence, S(ρ) := Tr(ρ log ρ) is the von-Neumann entropy, and ρ diag denotes the state obtained from ρ by removing all off-diagonal elements. For an isothermal process, the non-equilibrium free energy characterises the optimal amount of work that can be extracted from a given system with the help of a thermal reservoir [55,56]. This result also holds for a squeezed thermal bath [30]. In particular, we can show that the work extracted in this case is directly related to m and is given by Therefore, by squeezing a bosonic mode, coherence between energy eigenstate is created, and this can be quantified by the parameter m. In our description, the presence of non-diagonal elements in the covariance matrix and a negative classicality function will be an indicator of non-classicality. Furthermore, the amount of non-classicality present in a given bosonic system can be quantified by the classicality function (10).

IV. QUANTUM HEAT ENGINES
We now turn our attention to investigating the performance of a heat engine that operates in two different cycles. The first is an Otto cycle with isentropic strokes corresponding to a modulation of the squeezing parameter. In this cycle, the classicality of the system increases during its interaction with the squeezed thermal bath since the work parameter (squeezing) is kept constant and the temperature increases. To understand how the degree of non-classicality affects the performance of the heat engine, we will also consider a second cycle, where instead of keeping a fixed squeezing parameter during the isothermal process, the classicality is maintained constant so that we can analyse the trade-off between the degree of non-classicality and the efficiency of the thermodynamic protocol.  3. (a) Efficiency, ηo, of the Otto heat engine as a function of the squeezing parameter. The efficiency does not depend on the temperature, but only on the squeezing parameter. The difference in colour indicates an increasing degree of non-classicality with r. In (b) we have the classicality function (10) in the Otto cycle as a function of the squeezing parameter for T1 = ω/kB and T2 = 2 ω/kB. In the inset, we present the diagrammatic representation of the Classicality function over the cycle. The blue and red-shaded denotes the classical, and non-classical regime, respectively.

A. Otto Cycle
The Otto cycle consists of four consecutive strokes, two isochoric and two isentropic processes, as shown in Fig. 2(b). During the isentropic stroke, an external agent modulates the squeezing parameter of the harmonic oscillator between 0 and r. Heating and cooling results from the coupling with the two heat baths at temperatures T 1 = ω/k B and T 2 = 2 ω/k B . The hot squeezed thermal bath is prepared with the same squeezing parameter as the system after the first stroke. The state of the working fluid that performs the cycle is denoted by ρ A , ρ B , ρ C , and ρ D , respectively. Furthermore, the mean number of photons will be denoted byn The cycle starts with the working fluid in equilibrium with the cold thermal reservoir. During the first stroke, the working fluid is decoupled from the reservoir and a unitary squeezing operation S is applied to the system, resulting in isentropic compression. The work done in the single mode reads where r is the squeezing parameter of the transformation. In the second stroke, the bosonic mode is put in contact with the squeezed thermal bath, and the squeezing parameter remains constant, resulting in an "isochoric" process. Then, the mode relaxes to the steady state ρ C . In this step, the heat exchanged between the system and the squeezed thermal bath is During the third stroke, the working fluid is again detached from the reservoir, and a unitary unsqueezing operation is applied to the mode S † , bringing its squeezing parameter back to 0. Consequently, the amount of work extracted in this stroke is given by Finally, the cycle is closed by bringing the bosonic mode into contact with the cold thermal bath at temperature T 1 , and consequently the state relaxes to ρ A . During the last isochoric process, the heat transferred from the cold reservoir to the system yields the following results The work over the cycle is given by the two isentropic contributions and the efficiency of the engine, defined as the ratio between the output work and the heat absorbed from the hot reservoir, can be obtained from the previous results. Specifically Eqs. (20) and (23) allow us to compute the efficiency As already expected from Ref. [22], the efficiency increases with the squeezing parameter (see Fig. 3 (a)) and for higher values of squeezing, the efficiency will approach unity, but will never exceed it. Note that, compared to the classical case, the efficiency does not explicitly depend on the temperature of the reservoirs.
The behaviour of the classicality function throughout the cycle is plotted in Fig. 3 (b). Observe that in the first stroke (AB), non-classicality is added to the single mode by the action of the squeezing operator. In other words, a classical state with no coherence was taken to the non-classical regime (red region), now characterised by a non-classical state with coherence in its energy eigenbasis.
In the next stroke (BC), the squeezing parameter stays constant, while the temperature of the bosonic mode increases, meaning an increase in the classicality function. This can be interpreted as if the bath had consumed part of the non-classicality that was added in the first stroke. The last two strokes (CD) and (DA) correspond to the unsqueezing and the lowering of the temperature of the bosonic mode. Thus, in this step, the classicality function increases until its maximum and reaches its initial value. From Fig. 3 (b), we see that non-classicality is an essential resource in thermodynamics as the efficiency of the protocol increases with the degree of squeezing. More precisely, with the degree of non-classicality.
Note that this protocol does not fully explore the non-classicality added in the bosonic mode as the interaction with the hot reservoir causes an increase in the system's classicality. This motivates us to introduce a novel stroke type in which the degree of classicality is kept constant during the system-reservoir interaction.

B. Generalized cycle
We consider the same setup as before: a quantum heat engine operating between a cold thermal bath at temperature T 1 = ω/k B , and a hot squeezed thermal reservoir at T 2 = 2 ω/k B , with squeezing parameter r R (see Fig. 2 (c)). The main difference is that now we want to preserve the same degree of non-classicality added during the interaction of the system with the hot reservoir: this implies that the squeezed thermal bath must be prepared with a different squeezing parameter than the system; essentially, because the system will absorb heat from the hot bath, so its temperature and classicality increases. Consequently, to compensate for the addition of classicality, the system must have its squeezing parameter modulated from r t to r R , where r t is the squeezing parameter related to the first stroke. Consequently, during the interaction between the system and the reservoir, we impose the following condition: The cycle starts with the working fluid in thermal equilibrium with the cold thermal bath, and the procedure followed in the first stroke is the same as in the Otto cycle: the working fluid is disconnected from the cold bath and its squeezing parameter modulated from 0 to r t . In the second stroke, the bosonic mode is put in contact with the squeezed thermal bath. To keep the classicality constant, the stationary state of the system is slowly varied, so that the condition (25) is satisfied. According to the definitions (16), in the second stroke, work, and heat can be obtained from Moreover, to satisfy condition (25), work is performed on the system during its interaction with the squeezed thermal bath. The relation between r t and r R , can be obtained using Eq. (25) and gives In the third stroke BC, the system is again detached from the reservoirs and an unsqueezing operation, S † , is applied to the mode, changing its squeezing parameter adiabatically back to 0. The cycle is closed by putting the bosonic mode in contact with the cold thermal reservoir and relaxing back to ρ A . Once the work and heat are calculated for each stroke, the efficiency of the cycle is determined. Differently from the first protocol, now the efficiency depends on the temperatures and does not have a straightforward expression as before, A comparison between the efficiency of both cycles is shown in Fig. (4)(a), where the generalised cycle is represented by a black line, and the Otto cycle by a dashed grey line. The generalised cycle is always more efficient than the Otto cycle, since the classicality of the system is kept constant along the third stroke. Consequently, work must be done on the system to satisfy FIG. 4. (a) The thermodynamic efficiency of the generalised cycle (black curve), Otto cycle (dashed grey) curve and the Carnot efficiency (dashed white) as a function of the squeezing parameter for T1 = ω/kB and T2 = 2 ω/kB. Also illustrated, a phase diagram with the three regimes of the cycle (i, ii, and iii). The three different colours denote the degree of non-classicality of the bosonic mode. While blue and red indicate the classical, and non-classical behaviour of the mode, the grey colour represents that the mode is non-classical for T1, but classical for T2. In (b) the Classicality function (10) in the generalised cycle as a function of the squeezing parameter for T1 and T2. The blue and red-shaded colours denotes the classical, and non-classical regime, respectively. the condition (25), which means higher extraction of work in the third stroke. Note that in the Otto cycle, the hot bath consumes part of the non-classicality that is given to the bosonic mode. Here, we prevent this by imposing the condition (25).
As can be observed in Fig. (4)(a), the non-classical character of the system depends on its temperature and squeezing parameter. Consequently, for a single-mode operating at two different temperatures, there are two different critical squeezing parameters, denoted by r c1 and r c2 . This implies that this cycle can be investigated in three different regimes. That is, when the state of the single-mode is (as summarised in the phase diagram of Fig. 4(a)): i. Classical for both temperatures.
ii. Non-classical for T 1 , but classical for T 2 .
iii. Non-classical for both temperatures.
Region (i) corresponds to the case where the system is classical, C(r, T ) ≥ 0, but non-passive and only the degree of nonpassivity plays a role. Although there is squeezing, the uncertainty in both quadratures is above the shot-noise limit, and the squeezing can be seen as a modulation of the single-mode frequency from ω 1 to ω 2 , where ω 2 > ω 1 , in a similar fashion to previous works.
Region (ii) corresponds to the case where the state of the system is non-passive and non-classical for T 1 , but classical for T 2 . Compared to Region (i), the non-passivity combined with the non-classicality of the state increases the efficiency of the protocol. Since the non-classical behaviour depends on the external temperature, the single mode undergoes two transitions during the cycle: classical to non-classical (first stroke) and non-classical to classical (third stroke) due to the temperature increasing (see Fig.4 (b)). Observe that Carnot's efficiency (dashed grey line) is surpassed in region (ii), where the working fluid starts the cycle with a non-classical character. Our understanding is that the non-classicality added is used to perform the cycle more efficiently and helps to beat Carnot's limit. However, this comparison is not totally appropriate since the problem in consideration is formulated in a different context from which Carnot's cycle was proposed.
Finally, region (iii) is the most efficient. In this particular case, the state of the system is non-passive and non-classical for both temperatures. As a matter of comparison, the classical Carnot efficiency is also plotted (see Fig.4). As we can see in Fig. (4)(a), with increasing squeezing parameter, the efficiencies in the Otto (grey dashed line) and generalised cycle (black line) become closer to each other. This happens because the absorbed heat increases with the squeezing parameter, while the rejected heat does not. Therefore, for higher values of squeezing (r >> r c2 ), the difference between the two cycles becomes tight. In other words, the degree of non-classicality is so high that the type of cycle does not affect the performance of the quantum heat engine.
It is worth mentioning that in both cycles, we see a monotonic increase in efficiency with the squeezing parameter. There is no visible change in the respective curve when the squeezing parameter reaches the critical value corresponding to the transition from a classical to a non-classical regime. The lack of signature occurs because we are dealing with two different temperatures, and our definition of classicality is formulated for a single temperature. Moreover, the efficiency of the heat engine takes into account all four strokes, and in each, we are interested in the energy exchange. Consequently, the classical or non-classical character of the working fluid contributes only to the average energy, so a discontinuity or a signature discriminating between these two regimes, or their change, is not expected.

V. CONCLUSION
A central aim in quantum thermodynamics is the search for quantum advantages in a given thermodynamical task. Here, we have seen an example of a quantum advantage using a simple measure of classicality for squeezed states. This work presents an analysis of how non-classical features of light will play a role in thermodynamics from a P -representability perspective. Specifically, squeezed states have been shown to be a resource for thermodynamics. This, in turn, allows us to make use of the framework developed in [22,30] to track the trade-off between thermodynamic efficiency and the non-classicality presented in the work substance. The thermodynamic implications of squeezed states on heat engine performance have been studied in numerous previous works, especially in the case of the harmonic quantum Otto engine. However, the external control parameter modified during the isentropic stroke is the frequency, whereas here we is the squeezing parameter. In some sense they are equivalent, as the variation of squeezing can be understood as a variation of frequency, as the energy term E = ω(n sth + 1/2) = ω(n th + 1/2) cosh 2r, can be recast as the energy of a harmonic oscillator, E = ω (r)(n sth + 1/2), with ω (r) = ω cosh 2r, a squeezing dependent frequency. For obvious reasons the same is not true about the coherences m sth , as they are absent for the harmonic oscillator. Our treatment allowed us to identify a novel stroke type in which the degree of classicality is held constant, allowing for the introduction of a new class of generalised heat engine cycles.
In this study, we provided a simple expression, namely the classicality function, which allows us to distinguish classical and non-classical states. Since the squeezing parameter is the route to non-classicality, we treated it as a working parameter. Then, two different cycles where the classicality function is explored were analysed. Our main result shows that the efficiency of a quantum heat engine is enhanced when it operates in the non-classical regime, whereas it is less efficient in the classical one. An important point to be mentioned is the role of passivity/non-passivity. In our description, it can be straightforwardly observed that the anti-squeezing operation S induces maximal work extraction, i.e., ergotropy. However, here we have the interplay between non-passivity and non-classicality. The discussion of the generalised cycle was a way to circumvent this issue and to explicitly show that the degree of non-classicality is indeed a resource for quantum heat engines.
Our work opens many potential avenues for future research. First, in the current work, we have focused exclusively on a single mode as work substance, so a natural question concerns its extension for a multi-mode scenario where entanglement can also be explored. Moreover, since our discussion is only in terms of the equilibrium regime, this manuscript leaves open the possibility of a significant amount of follow-up research. Finite-time analysis (either in the endoreversible or fully nonequilibrium regime) would open up the possibility of studying the consequences of the degree of classicality on a wide range of other engine characterisations, including power output and efficiency at maximum power. As shown in [29], quantum engines employing squeezed states are also experimentally accessible, making experimental implementations of this work a near-term possibility.
derived [45]. For this purpose, let us rewrite the P -function in terms of the parametrisation that we have been discussing, P (a) = da e −a † σza χ(a ) , where σ z is the z-Pauli matrix, and χ N (a) is the normally-ordered characteristic function that can be rewritten as Substituting Eq.(A3) into (A2), and using the fact that −a † σ z a = −a † σ z a [45], the P-function can be written as P (a) = da e a † (V− 1 2 )a = 1 where σ x is the x-Pauli matrix. It can be verified that (σ z a) T = aσ z and aσ z T = −a † σ z , thus where From this result we may conclude that, for a single-mode Gaussian state described by a covariance matrix V is said to be P -representable ifn Note that a P -representable state, P must be non-negative, which means that (V − 1/2) ≥ 0. This condition requires eigenvalues greater than or equal to zero, or alternatively, a non-negative determinant.
It is important to emphasise thatn is not necessarily the mean number of thermal photonsn th , but the mean number of photons corresponds to the distribution in which the system is described.