Measurement Based Quantum Heat Engine with Coupled Working Medium

We consider measurement based single temperature quantum heat engine without feedback control, introduced recently by Yi, Talkner and Kim [Phys. Rev. E 96, 022108 (2017)]. Taking the working medium of the engine to be a one-dimensional Heisenberg model of two spins, we calculate the efficiency of the engine undergoing a cyclic process. Starting with two spin-1/2 particles, we investigate the scenario of higher spins also. We show that, for this model of coupled working medium, efficiency can be higher than that of an uncoupled one. However, the relationship between the coupling constant and the efficiency of the engine is rather involved. We find that in the higher spin scenario efficiency can sometimes be negative (this means work has to be done to run the engine cycle) for certain range of coupling constants, in contrast to the aforesaid work of Yi, Talkner and Kim, where they showed that the extracted work is always positive in the absence of coupling. We provide arguments for this negative efficiency in higher spin scenarios. Interestingly, this happens only in the asymmetric scenarios, where the two spins are different. Given these facts, for judiciously chosen conditions, an engine with coupled working medium gives advantage for the efficiency over the uncoupled one.


Introduction
Unification of seemingly different heat engines in terms of efficiency started in the early 19th century with Sadi Carnot [1]. From then onward, converting heat into useful work with increasing efficiency got a thrust in practical and industrial territories. For a standard heat engine, working cyclically between two heat baths of temperature T 1 and T 2 (T 2 < T 1 ), efficiency of the engine is upper bounded by η = 1 − T 2 /T 1 , the Carnot efficiency [2]. The second law of thermodynamics puts this fundamental limitation on the extent of work that can be converted from heat. The laws of thermodynamics are empirical and were first adopted for classical macroscopic systems. Naturally, the validity of the laws of thermodynamics are questionable and subject to verification in the quantum regime. Moreover, quantum mechanics gives the dynamical viewpoint of thermodynamics [3,4], describing the emergence of thermodynamic laws from quantum mechanics.
The idea of quantum heat engine first appeared in a paper by Scovil and Schulz-DuBois [5], where the authors demonstrated that three level masers can be treated as a working medium for heat engines. Today, the study of heat engines in quantum domain is an active area of research both due to the gradual miniaturization of current technology as well as its theoretical richness.
Within the quantum engines scenario, quantum analog of the classical heat engines [6,7] and many other generalizations [8][9][10][11] have been studied. Analysis of finite power quantum heat engines also have a significant amount of literature (e.g., [12][13][14][15][16][17][18]). With the onset of the quantum effects, many interesting phenomena such as the increase of efficiency beyond Carnot's limit [19][20][21] may occur. However, it is not in contradiction with the second law of thermodynamics. The compensation comes from carefully accounting all the work costs. Gardas and Deffner [22] showed that if one accounts for the work cost to maintain the non-equilibrium reservoir, Carnot's limit cannot be surpassed. However, understanding quantum thermodynamic machines [23] and the role of quantum effects [24][25][26][27][28][29] in quantum thermodynamics is far from fully understood. Where quantum effects set a limit to our ability [30] for practical purposes and where we can actually use the quantum resources are still parts of ongoing research field. The approach to resource theory of quantum thermodynamics [31] tells us about the fundamental corrections to the laws of thermodynamics, setting the limit to the performance of quantum heat engines.
Previously, it was shown that [32,33] a quantum Otto engine with coupled working medium leads to a higher efficiency than that of an uncoupled one. In addition, in information heat engine, e.g., the Szilard engine, where exploiting the information one can extract work from an engine operating at a single temperature [34][35][36][37], it was shown that [38][39][40][41] entanglement can be used to extract work beyond the limit, which is possible using classical correlation only.
Recently, in [42], the authors introduced a new kind of single temperature quantum heat engine without feedback control. The essential part of the engine which replaces feedback is a non-selective quantum measurement on the working medium, changing the average energy of the system, and, thus, enabling one to extract useful work. This engine is similar to a quantum Otto cycle [43][44][45][46] with one thermalization stroke being replaced by a quantum non-selective measurement, whereas, in Maxwell's demon and Szilard engine [34,35], work is extracted from a single heat reservoir using feedback control. Another version of Maxwell's demon engine was introduced in [47,48], where without the presence of any thermal bath, work can be extracted using measurement and feedback control. Thus, quantum measurement plays an important role in quantum thermodynamics. Energetic cost for performing a measurement [49][50][51][52][53] and using the average energy change due to the measurement for extracting useful work are two important facets of quantum thermodynamics. In a subsequent work [54], the authors calculated the detailed fluctuation of work and heat in the above said measurement driven single temperature heat engine without feedback control as well as considered the finite power scenario.
In this paper, we analyze the role of coupled working medium in this single temperature measurement driven quantum heat engine without feedback control [42]. Taking a coupled one-dimensional Heisenberg model as the working medium, we show an advantage for the efficiency over the uncoupled one.
First, we start with a one-dimensional Heisenberg model of two spin-1/2 particles and then generalize that for two spin-d particles, where d can take values of 1/2, 1, and 3/2. We note that, for different choices of non-selective measurements, the efficiency of the heat engine changes. In addition, in the higher-dimensional scenario, another interesting feature is observed: we can either extract work from the engine cycle or have to invest work to run the cycle depending upon the spin configuration we choose. By judiciously choosing all the conditions, such as measurement choices, coupling constant and the dimension of the Hilbert spaces, one can optimize the engine performance, which is better than the uncoupled one in terms of efficiency.
The paper is organized as follows. In the next section, we give a short introduction about the single temperature measurement-based heat engine as introduced in [42]. In Sections 3 and 4, we present our result for the coupled measurement based heat engine taking the working medium to be the Heisenberg model of two spin half particles. In Section 5, we consider the higher-dimensional scenario. In Section 6, we present an analysis of the global and local work. Finally, we conclude in Section 7.

Single Temperature Measurement Driven Quantum Heat Engine without Feedback
In this section, we briefly discuss the recently introduced measurement based single temperature quantum heat engine without feedback control [42]. It is very similar to the Otto cycle except for one thermalization step, which is replaced by a non-selective quantum measurement. If it had been a classical system, in principle, there would be no subsequent effect of the measurement on the system. However, the quantum mechanical system is generally disturbed by measurement and hence average energy of the system changes. Judiciously choosing the measurement operators, as discussed in [42], we can extract work form this type of engine. We now briefly describe the engine cycle.
The working medium of the engine has a Hamiltonian H(λ), which is a function of an external control parameter λ. The system starts from a thermal state of temperature T. This can be achieved with the help of a heat bath of temperature T, which is the only heat bath to be used throughout the action of the engine. The system is brought to the contact with the bath and the system allowed to thermalize. After a long enough time, when the system attains equilibrium thermal state, the heat bath is detached and it gets ready for the first cycle of our heat engine. Thus, the initial state of the system is ρ int = e −βH(λ int ) /Z = ∑ n (e −βE n (λ int ) /Z) |n(λ int ) n(λ int )|, where Z = ∑ n e −βE n (λ int ) , while |n(λ) and E n (λ) are, respectively, the nth eigenstate and eigenvalue of the Hamiltonian H(λ). Now, the engine strokes are as following.
First stroke: The first stroke of the cycle is an adiabatic compression process. The working system is isolated from the heat bath and the Hamiltonian is changed quasi-statically from H(λ int ) to H(λ f in ), with initial occupation probabilities of the state remaining unchanged. For a system defined by a density matrix ρ and Hamiltonian H, its internal energy or average energy is defined as U = Tr [ρH]. Change in internal energy is the sum of two contributions [55], one is heat, defined as dQ = Tr [Hdρ], and the other is work, defined as dW = Tr[ρdH]. However, this identification of heat and work is not always valid, especially in strong system bath coupling [56,57]. We follow here the above definition. Although, during a general adiabatic process, the state of the working medium changes, for the model of engine cycle we consider in this paper, the state of the working medium does not change throughout the adiabatic stroke. Hence, in the first stroke, change in the internal energy of the system is , which can also be written as where p n (λ int ) = (e −βE n (λ int ) /Z). If this is positive, then this is the energy gained by the system. Thus, the average work that is extracted from this stroke is −W 1 .
Second stroke: The next stroke is the most crucial and special one, which involves a non-selective measurement. A measurement [58] corresponding to an observableĜ can be described by a POVM, {G n }, where G n ≥ 0 are the POVM effects, ∑ n G n = 1 and Tr[ρG n ] is the probability of getting nth outcome denoted here as α n . If in addition the POVM elements satisfy G m G n = δ mn G n , then they are orthogonal projectors and the observableĜ can be written asĜ = ∑ n α n G n (spectral value decomposition), where α n s are now the eigenvalues ofĜ. Equivalently, measurement can be completely described by a set of measurement operators {M n }, with M † n M n = G n . If the state of the system before the measurement is |φ , then corresponding to the nth outcome, the state of the system after measurement will be M n |φ √ p n or M n ρ int M † n p n , where ρ int = |φ φ| and p n = Tr[ρ int M † n M n ] = φ| M † n M n |φ is the probability of getting nth outcome. In a non-selective measurement, i.e., if the outcomes of the measurements are not recorded, then the state after measurement is ∑ n M n ρ int M † n . Thus, after the first stroke, we do a non-selective measurement described by the measurement operators {M α } on the system state giving the post-measurement state as ρ M = ∑ α M α ρ int M † α . In this stroke, the Hamiltonian of the system is unchanged at H(λ f in ). Thus, the average energy change of the system is given by , which is reminiscent of heat. This can be written as [42], where T m,n = ∑ α | n(λ f in )| M α |m(λ f in ) | 2 is the transition probability from a eigenstate labeled n before the measurement to an eigenstate labeled m after the measurement, and H M (λ f in ) = ∑ α M α H(λ f in )M α . As shown in [42], Q M is always positive, which follows from the properties of the transition matrix, T m,n = T n,m and ∑ n T m,n = 1. This fact is illustrated more explicitly in the next section. It is also noted that, whenever the Hamiltonian H(λ f in ) of the system does not commute with the measurement operators, we get a nonzero Q M .
Third stroke: This is the second adiabatic process. The Hamiltonian H(λ f in ) is very slowly changed back to the initial Hamiltonian H(λ int ), with unchanged occupation probabilities of the state. As in the previous adiabatic stroke, the average change in energy of the system is W 2 = Tr[ρ M (H(λ int ) − H(λ f in ))], as the state of the system is unchanged throughout the stroke. This is T m,n is the probability of finding the nth eigenstate of H(λ f in ) in ρ M . Thus, the work extracted form this adiabatic stroke is −W 2 , which follows from the similar argument given in the description of the first stroke.
Fourth stroke: In this last stroke of the cycle, the system is brought into contact with the heat bath of temperature T, while keeping the Hamiltonian fixed at H(λ int ) and allowed to thermalize, until it goes back to the initial thermal state ρ int . Thus, heat transfer for this stroke is given by , as the Hamiltonian is fixed in this stroke. This can be written as and is shown to be negative in [42]. This means that heat is going to the heat bath from the working medium at this stage.
The signs of each of these quantities, e.g., W 1 , Q M , W 2 , and Q T , are analyzed explicitly in the next section in the case of coupled working medium. Thus, the whole cycle is, similar to energy Q M , taken by the system, doing a work −(W 1 + W 2 ) and dumping energy Q T to a heat bath. Thus, we have Q M + Q T = −(W 1 + W 2 ), correctly depicting the first law of thermodynamics, i.e., energy conservation. Efficiency of the heat engine is given as the ratio of extracted work −(W 1 + W 2 ) over the average energy change Q M in the measurement stroke: Yi, Talkner and Kim [42] showed that the extracted work −(W 1 + W 2 ) is always positive, relying on the fact that work strokes are either adiabatic expansion or compression of the working medium. However, we show in the next section that, for a coupled working medium, their way of reasoning does not hold well and there can be instances where −(W 1 + W 2 ) is negative.

Coupled Single Temperature Measurement Engine
In this section, we present an analysis of a coupled measurement-based single temperature heat engine. We consider the working medium of the system to be a one-dimensional Heisenberg model of two particles with the following Hamiltonian: For two spin half particles, S A = S B = 1 2 σ, where σ = (σ x , σ y , σ z ) are the Pauli matrices. Thus, in this case, we can write the Hamiltonian as, where J is the coupling constant and B is the external magnetic field. The entanglement between two qubits for this model has been studied in [59]. J > 0 and J < 0 cases correspond to the anti-ferromagnetic and ferromagnetic interactions, respectively. In this paper, we restrict ourselves to the anti-ferromagnetic case only. Eigenvalues and eigenstates of this Hamiltonian are listed in Table 1, are the eigenstates of σ z . As described above, the engine cycle of the measurement based heat engine has four steps. The first stroke of the cycle is an adiabatic compression, where the Hamiltonian of the working system, as described above, is quasi-statically changed from an initial parameter value to a final parameter value. External magnetic field B is the parameter of the Hamiltonian here. It is changed quasi-statically from the initial value B 1 to the final value B 2 . As this process is done adiabatically, the state of the system has same occupation probabilities throughout the stroke. At the beginning of the cycle, we take the working medium of the heat engine to be in a thermal equilibrium state of temperature T. Thus, here, ρ int = ∑ 4 n=1 P n |ψ n (B) ψ n (B)|, where P n = exp(−E n /k B T)/Z, Z = ∑ n exp(−E n /k B T), E n s and corresponding |ψ n (B) are given in Table 1. Then, in the second stroke of the cycle, the Hamiltonian of the system is kept unaltered but a measurement of an observable is performed on the system. As discussed above, the observable has to be non-commuting with the Hamiltonian to get a positive work output. In this case, we have a distributed system and we show that the efficiency of the heat engine depends on the local measurements we are performing. The third stroke is again an adiabatic process changing the external magnetic field B 2 back to B 1 . The final stage of the cycle is a thermalization step and in this stage the system is brought to contact with a heat bath of the starting temperature T and the system allowed to thermalize for a sufficiently long time, after which it again goes back to the initial thermal equilibrium state. Now, as mentioned above, the initial state of the system is a thermal state of temperature T, where P n = exp(E n /k B T)/Z, Z = ∑ n exp(E n /k B T), E n and |ψ n (B) s are the energy eigenvalues and eigenstates, respectively, as listed in Table 1. From now on, we take k B T = 1 throughout the paper. The energy eigenvalues and hence the probabilities P n depend on the changing parameter, which is the external magnetic field B, but the eigenstates are independent of the parameter. Thus, we omit the dependence of B from the notations of the eigenstates. The average work in the first adiabatic stroke is, as the state remains in its instantaneous eigenstate with same probability. For the system we consider and the initial thermal state of the system, we have

Eigenvalues Eigenstates
The next stroke is the non-selective measurement part. We can choose any arbitrary observable for measurement with only constraint that the observable must be non-commuting with the Hamiltonian. For the time being, we restrict ourselves for projective measurements. In our case, for coupled measurement based heat engine, we take most general measurement operators as, where |+ a + a | and |− a − a | are the eigenstate projectors for the observable σ.â for one party and |+ b + b | and |− b − b | are the eigenstate projectors for the observable σ.b for the other. Now, if the initial state of the working medium is given by Equation (4), then, for a non-selective measurement given by the above measurement operators, the post measurement state will be, The average energy change of the system during this measurement stroke is given by, where T m,n is the transition probability, written as For the most general form of the measurement operators, the expressions are quite complicated. Thus, we first write the expressions for some special choices of measurements. The first case is when For σ.â = σ y and σ.b = σ z , the expression for the average change in energy remains same. Another case is when σ.â = σ x and σ.b = σ y . In this scenario, There are many other different choices of measurement operators, such as both measuring σ x or σ y or σ z , etc. In each case, the expression for Q M will change accordingly. This average energy change is similar to heat in the conventional quantum heat engine. Now, the third step of the cycle is again an adiabatic process, where the magnetic field B 2 is changed back into B 1 very slowly. For this part of the cycle, the work done is, where P n = ψ n | ρ PM |ψ n = ∑ m T m,n P m (B 1 ) is the probability of getting nth eigenstate in the post measurement state. When σ.â = σ x and σ.b = σ z , the expression of this work for our system is, The last step of the cycle is to bring the system in contact with the heat bath of temperature T and let it thermalize back to the initial thermal state ρ int . Heat exchanged in this step is given by, Again, for σ.â = σ x and σ.b = σ z , heat dumped into the heat bath in this last step is, Thus, as discussed in the previous section, the total work that can be obtained from the cycle is given by the sum of the works that can be extracted in the first and third strokes. For σ.â = σ x and σ.b = σ z , total extracted work is given by, The quantities we have calculated so far are global (from the perspective of two systems together), i.e., global heat, global work or global energy change. In the next section, we discuss the global efficiency of the heat engine. In addition, the measurements we choose are all projective measurements. We investigated some cases of POVM, namely the SIC POVM [60], and some other examples. However, in all cases, we found that projective measurements are more effective so far as the efficiency is concerned. Thus, we be restricting ourselves with projective measurements only.

Efficiency of the Heat Engine, Global Analysis
In this section, we evaluate the efficiency of the measurement based coupled heat engine and compare it with the uncoupled one. Before that, it is necessary to determine the signs of the quantities Q M (Equation (12)), Q T (Equation (18)), W 1 (Equation (5)) and W 2 (Equation (16)) for this coupled working medium. The average energy change Q M of the system during the measurement stroke, given in Equation (12), can also be written as [42], This is obtained by employing the properties of the transition matrix, T m,n = T n,m and ∑ m T m,n = 1. From the fact that T m,n ≥ 0 and the equilibrium occupation probability P n for an energy level E n decreases with the increase of energy E n , it turns out that Q M ≥ 0. That means, in the measurement step, heat enters into the working medium. Similarly, the expression for Q T given in Equation (18) can be written as, By similar arguments made for Q M , it is evident that Q T ≤ 0, which means that heat goes from the working medium to the heat bath of temperature T at the last step of the cycle. Now, to determine the signs of W 1 and W 2 , let us first write down the alternative expression for the work as done in [42].
where ∆ α m,n denotes the difference between the mth and nth energy eigenvalues and is given by, In [42], the authors designed the adiabatic strokes to be compression and expansion, in the sense that spacing between the energy levels of the Hamiltonian either increases or decreases. Then, they argued that for the compression stroke, ∆ f m,n ≥ ∆ i m,n , and as the canonical probability decreases monotonically with the increase of energy, for ∆ i m,n > 0, P m (B 1 ) − P n (B 1 ) is negative. Together with the fact that T m,n is positive, every term in the above expression of W is non-positive and hence the total work extracted W t = −W is always positive. Similar arguments hold for the expansion stroke. However, for a coupled working medium, this argument does not hold good. In the presence of coupling J, the uniform increase or decrease of spacing between the energy levels does not happen. For example, for the two spin-1/2 scenario, the energy eigenvalues are −6J, 2J − 2B, 2J, 2J + 2B, which are ordered from low to high for J > 0 (and for 4J ≥ B). Throughout the paper, we consider the anti-ferromagnetic case, i.e., J > 0, as stated above. Now energy difference between lowest two energy levels is 8J − 2B and highest two energy levels is 2B. Thus, for a fixed J, with increasing B, spacing between lowest two energy levels decreases, where spacing between highest two energy level increases. Thus, as a whole, we cannot say that the adiabatic stroke is compression or expansion for the working medium and consequently the line of reasoning in [42] does not hold well here. We simply say that the work strokes are the first and second adiabatic work strokes.
Thus, whether the extracted work is positive or negative would depend upon the structure of the energy levels of the Hamiltonian. Let us illustrate this by considering the working medium to be a two spin-1/2 system. In Equation (23), some terms are negative and some terms are positive. Sign of the total work depends upon the positive and negative contributions of these terms. We denote each term of the expression as W mn = 1 and note that W mn = W nm . We plot different terms with the coupling constant J in Figure 1 for B 2 > B 1 . We notice that −W 12 (−W 21 ) gives negative contributions, −W 14 (−W 41 ) and −W 23 (−W 32 ) give positive contributions and the other terms are zero. The balances between these terms decide the sign of the work output. In this case, it is positive. We also plot the contributions for first and second adiabatic work strokes in Figure 2 for B 2 > B 1 , and notice that the contribution of the first work stroke is positive, whereas the contribution of the second work stroke is negative but their sum is positive and hence we extract positive work. It is important to note that the arguments behind Q M ≥ 0 and Q T ≤ 0 are independent of the choice of initial and final magnetic field B 1 and B 2 , respectively. Indeed, for both B 2 > B 1 and B 2 < B 1 , Q M is always positive and Q T is always negative. It is the extracted work W t , whose sign depends upon B 1 and B 2 . From Equation (20), it is evident that, when B 2 > B 1 , we extract positive work, i.e., W t ≥ 0 and, when B 2 < B 1 , extracted work is negative for all J. This can be seen in Figure 3. This plot shows that when the values of B 1 and B 2 are interchanged, i.e., B 1 = 4 and B 2 = 3, the nature of the plots are just opposite to each other. From the definition of efficiency as the ratio of total work extracted (W t ) over Q M , η = W t /Q M , it is evident that the efficiency can be negative only if the extracted work W t is negative, which means that we cannot extract work from the engine but have to do work to run it. We have seen that, for B 2 > B 1 , the extracted work is positive. For the rest of the analysis, we stick to this scenario, i.e., B 2 > B 1 . From the expressions derived above, we have different efficiencies depending upon the measurement choices. When σ.â is σ z and σ.b is σ x , we have, Thus, the efficiency depends upon the coupling constant J, where J = 0 corresponds to the uncoupled scenario with efficiency (1 − B 1 B 2 ). As J increases, for a certain range of J, efficiency is also increased over the uncoupled value. Expression for the efficiency remains same as in Equation (24) for σ.â = σ y and σ.b = σ z , For, σ.â=σ x and σ.b=σ y ,  Now, let us examine those cases where same observables are being measured on both sides, e.g. σ.â = σ.b = σ z , σ x or σ y . When both observables are σ x , we have, Exactly the same expression is obtained when both observables are σ y . When both observables are σ z , the work contributions from two adiabatic branches are equal and opposite to each other. Consequently, the total work done is zero and hence the efficiency is zero. After calculating the efficiency for different measurement choices, we plot them in Figure 4 together with B 1 = 3 and B 2 = 4. Next, we plot the efficiency for a fixed observable σ z on one side and varying the parameters for the observable on the other side. We can write σ.m = sin θ cos φσ x + sin θ sin φσ y + cos θσ z . In this case, we calculate the efficiency and it turns out to be independent of the parameter φ but depends on θ. From the 3D plot in Figure 5, it is clear that, when θ is π/2, the efficiency is optimum and it is exactly equal to the case where σ x is measured for one spin and σ z is measured for the other. Our results show that using non-zero coupling J, we can actually get an advantage over the no-coupling scenario. We argue that, for B 2 > B 1 , the extracted work is positive and hence we get a positive efficiency. Now, we show that, for any B 1 and B 2 (B 1 , B 2 > 0) with B 2 > B 1 , the efficiency of the coupled engine is greater than that of an uncoupled one, for a certain range of J. To show this, we rewrite Equation (24) as,  As already mentioned, 1 − B 1 B 2 is the efficiency of an uncoupled (J = 0) engine. From the above expression, it turns out that the efficiency will be greater than that of the uncoupled one if From this inequality, it is evident that for any B 1 > 0 ( B 2 can have any value greater than B 1 ), we always get a positive value for J, below which we get an advantage for the efficiency over the uncoupled one. To give an example, let us consider a small value of B 1 , e.g., B 1 = 0.1. When B 1 = 0.1, we have to find the value of J for which the above inequality holds. Solving the Equation (28) for B 1 = 0.1, one can show that when J < 0.00166 approximately, the efficiency will be greater than that for the uncoupled one. For B 1 = 0.1, and B 2 = 4, the uncoupled efficiency is 0.975. Let us take J = 0.0014. For that, we have η = 0.975011 for the same values of B 1 and B 2 , i.e., B 1 = 0.1 and B 2 = 4. Thus, the coupled engine is still more efficient but it is in such a small region of that nature of J that cannot be seen from the plot unless the plot has a very fine scaling. With the increase in the value of B 1 , the cutoff value of J, above which the coupled engine is more efficient increases and the lower is the value of the ratio of B 1 over B 2 , the greater is the efficiency of the engine. For the uncoupled efficiency closer and closer to one, the region in which there is an advantage of the coupled engine will be narrower and narrower. However, in principle, we always have the coupled engine as the efficient one, compared to the uncoupled one. In addition, in Figure 4, it is evident that different measurement choices give different efficiencies for the heat engine. Here, choosing σ x and σ z or σ y and σ z as measurement operators, we get the maximum efficiency. Thus, judiciously choosing measurement operators is important for optimum performance of the heat engine.

Higher-Dimensional Case
In this section, we are interested in the higher-dimensional Heisenberg model as considered in [33], where one spin half particle is coupled to a spin s particle. It would be interesting to observe the effect of higher spin as an additional parameter along with the coupling constant. We have the system Hamiltonian, where S A = (S x A , S y A , S z A ) and S B = (S x B , S y B , S z B ) are two spin operators, J is the system-bath coupling constant and B is the external magnetic field. Two spin-1/2 cases have already been discussed. Now, the detailed calculation is carried out for S A/B = 1 and S A/B = 3/2. One can obviously go on to calculate the cases for 2 and 5/2 and so on, but the essential points can be observed by studying the following cases: S A/B = 1 and S A/B = 3/2. We first start with a spin-1/2 and a spin-1 operator, which is an asymmetric case in the sense that the spins on two sides are different. We call it symmetric when two spins at two sides are the same. We deal with the different cases one by one. We start with the asymmetric cases.

Asymmetric Case
First, we take S A = 1/2 and S B = 1. Spin operators for spin-1/2 particle are 1 2 σ, where σ = (σ x , σ y , σ z ) are the Pauli matrices. These spin operators for spin-1/2 particle form the fundamental irreducible representation of SU (2). Spin operators for the spin 1 particle are the three-dimensional irreducible representation of SU (2). With standard representation in canonical bases [61], the eigenvalues and the eigenstates of the Hamiltonian given in Equation (29) are listed in Table 2. As in the previous cases, for the measurement step of the engine cycle, we have a number of choices for the measurement operators and we explore those options in the previous sections for the case of two spin-1/2 particles. Now, for brevity, we consider one particular measurement setup and observe the effect of higher spin, such that this spin can also be a controlling parameter for the efficiency. We choose the measurement operators as, are the eigenstates of the operator S x for the spin-half particle A. In other words, we are doing measurement of the operator S x on side A for the spin-1/2 and measurement of the operator S z on the side B for spin-1. In the plots (see Figures 6-8), we compared different scenarios for the same measurement settings on the two sides. By the same measurement settings, we mean that on the spin half side the measurement operators will be the projectors constructed from the eigenstates of the operator S x and, on the higher spin side, it will be the projectors of the eigenstates of S z . For the above measurement operators, we calculate the quantities as work and heat (during measurement process and thermalization step) and evaluate the efficiency of the heat engine. The expression for the total work extracted for this case is,

Eigenvalues Eigenstates
From the expression above, one can see that only B 2 > B 1 will not guarantee the positivity of W t for all J, unlike the previous case of two spin-1/2 particles. Specifically, the work extracted is negative when 4 + 3e 2B1 − e 12J < 0. This is not the case for two spin-1/2 particles, where expression of the extracted work is such that B 2 > B 1 always gives positive extracted work and B 1 > B 2 gives negative extracted work. However, in the present case, for B 2 > B 1 , extracted work can be negative and, for B 1 > B 2 , extracted work can be positive when 4 + 3e 2B1 − e 12J < 0. As we show below, this is special for the asymmetric cases only. In Figure 6, we plot the efficiency and compare it with the spin half scenario for B 2 > B 1 , with same values of magnetic fields previously considered, i.e., B 1 = 3 and B 2 = 4. We note that ( Figure 6) the efficiency for the spin-1 scenario can be higher than that of spin half scenario for some range of non-zero values of J. For the uncoupled case, both engines give same efficiencies, which is 1 − B 1 B 2 . Another point to note is that the efficiency can go to negative for the spin 1 scenario. As discussed above, the negative efficiency comes entirely from the negative extracted work, because Q M is always positive. The total extracted work W t can now be negative starting from a certain value of J, as shown in Figure 7, for B 2 > B 1 , with B 2 = 4 and B 1 = 3, whereas Q M is always positive, as shown in Figure 8. Now, one can further analyze the work strokes and see which work stroke is contributing more for negative work. From the Figure 9, one can note that, as before, the first adiabatic work stroke gives the positive work output, whereas, for second adiabatic work stroke, we get negative work. If we take B 1 > B 2 , the situation is reversed. We have shown the case of B 1 > B 2 , with B 1 = 3 and B 2 = 4, in Figure 10. The plot is exactly opposite to the previous one with B 2 = 4 and B 1 = 3 (see Figure 9). Thus, in the range of J where the efficiency is negative, the situation appears as follows: the average energy Q M is entering to the working medium, W t work is being done and thereby heat Q T goes into the heat bath of temperature T. In some sense, one can associate a refrigerator action for this negative work scenario. In a conventional quantum Otto refrigerator, the system is first prepared in the thermal state with temperature corresponding to the cold bath. Then, in the first and third adiabatic strokes, a total work W is added to the working medium. In the second and fourth steps, Q 2 heat is taken from the cold bath and Q 1 heat is added to the hot heat bath and eventually cooling the cold bath more. The co-efficient of performance (COP) for the refrigerator is given as η COP = Q 2 /W. In our case, Q M is always positive and Q T is always negative. This indicates that we can consider a cold bath of effective temperature T 2 , such that, after the second stroke, we have, where ρ eq (T , B 2 ) is the thermal state at temperature T with the Hamiltonian H(B 2 ), given in Equation (29), where T = (T, T 2 ). Solving the above equation we can associate an effective temperature T 2 with Q M . Thus, now one can read the cycle as transferring heat from the cold bath of temperature T 2 to hot bath of temperature T and, for this, W t work has to be done. This means that the COP of this refrigerator action is Q M /(−W t ). Thus, to get a positive work output, we have to judiciously choose the value of the coupling constant J, such that efficiency is not negative. Then, we get an advantage for higher efficiency over the uncoupled one. Let us now consider the next asymmetric scenarios and see whether a similar trend, i.e, increase in efficiency and occurrence of negative efficiency is present or not. We start with the case where the spins on two sides are S A = 1/2 and S B = 3/2. With standard representation of spin-3/2 in canonical bases [61], the eigenvalues and eigenstates of the Hamiltonian in Equation (29) are given in Table A1 of Appendix A. Now, we choose the same kind of measurement operators as in the previous case, i.e., on the spin-1/2 side we measure S x and on the spin-3/2 side, S z : , and |0 B , |1 B , |2 B , and |3 B are the eigenstates of S z B corresponding to the eigenvalues 3, 1, −1, and −3, respectively. We also consider the scenario of S A = 1 and S B = 3/2. For this case, the eigenvalues and eigenstates of the Hamiltonian in Equation (29) is given in Table A3 of Appendix A. Again, we take the same measurement choices as before, i.e., measurement of S x spin operator on the side A and S z spin operator on the side B. We calculate W, Q M and the engine efficiency for each case. As in the previous asymmetric case, the condition B 2 > B 1 does not guarantee the positivity of the extracted work W t in both present cases. Starting from a certain value of J, W t can be negative for B 2 > B 1 and positive for B 1 > B 2 .
Thus, for asymmetric cases, B 1 > B 2 and B 2 > B 1 , both situations give rise to the extracted work to be negative staring from certain ranges of J. We plot the efficiency of the heat engine for the aforesaid three asymmetric cases together in Figure 11 for B 2 > B 1 with B 2 = 4 and B 1 = 3. B 1 > B 2 case can be calculated in a similar way. We observe that, for asymmetric situation, the efficiency becomes negative after a certain value of J. In addition, in Figure 11, we can observe that, as the difference of spins increases between the two sides, the efficiency goes to more negative value. Thus, from these observations, it is clear that, if the two spins on both sides are not the same, then efficiency can be negative. Another interesting feature to notice from the plot is that, within the range of J where efficiency is positive, higher differences of the spin values give larger gain in efficiency over the uncoupled one. We have to take correct coupling strength J, to have a higher but positive work output from these measurement-based coupled higher spin coupled heat engines. As in the case of S A = 1/2 and S B = 1, we also plot the work done in two adiabatic strokes for these two asymmetric cases in Figures 12 and 13 and note that the negative contribution in the extracted work is due to the second work stroke.

Symmetric Case
In this section, we consider the symmetric case. One particular characteristic to look at is whether the efficiency gets negative for symmetric situation as well. We show above one symmetric situation, namely the case of two spin-1/2 particles, where we see positive efficiency for all natures of J (see Figure 4). Let us investigate the case for higher spin symmetric situations. Three cases can arise for the symmetric scenario, if we restrict ourselves up to spin-3/2. Among these, two spin-1/2 case is discussed at the very beginning of the paper. The remaining two cases are the cases of two spin-1 and two spin-3/2 particles. Eigenvalues and the eigenstates of the Hamiltonian Equation (29) for these two cases are given in Tables A2 and A4, respectively, in Appendix A. We take the same measurement settings as before, i.e., measurement of S x spin operator on the side A and S z spin operator on the side B. We calculate W, Q M , and the engine efficiency for all these symmetric cases. Interestingly, for the symmetric cases, for B 2 > B 1 , we get the extracted work W t = −W to be positive for all J, and, for B 1 > B 2 , W t is negative for all J. This is the same as in the case of two spin-1/2 particles (which is the simplest example of a symmetric case). For two spin-1 particles, the expression for W t is given by, From the above expression for W t , it is evident that, whenever B 2 > B 1 , we have the extracted work is positive and, for B 1 > B 2 , the extracted work is negative. The situation for the other symmetric case is similar, i.e., for two spin-3/2 particles. One can also analyze the contributions of two adiabatic work strokes and find that sum of these two strokes gives rise to the extracted work a positive quantity for B 2 > B 1 . We plot the efficiency of the engine for all symmetric cases together in Figure 14 for B 2 > B 1 , with B 2 = 4 and B 1 = 3. Given the fact that Q M is always positive, for symmetric case, we always get positive efficiency for B 2 > B 1 . We also see that efficiency gets higher with the increase of spin value, although the efficiency decreases more quickly for higher spins with respect to the coupling strength J. Thus, from these observations, we can conclude that, for the asymmetric scenario, the work output is not always positive. However, for the symmetric scenario, this is not the case. For this latter scenario, we always get positive work output and hence positive efficiency for the heat engine with B 2 > B 1 . In addition, for higher spin scenario, the efficiency is always greater (as long as it is positive for the asymmetric case) than that of the two spin half case. Thus, along with the coupling J, spin also plays an important role for the increase in efficiency for the measurement driven single temperature coupled heat engine.

Local vs. Global Work
In this section, we briefly touch upon the status of "local" and "global" works and how they are related. Until now, we have been discussing the global aspect of the heat engine, i.e., global work output, global efficiency, etc. Here, global means that we consider the coupled working medium as a whole and do not take into account the subsystems' contribution. In the scenario of coupled quantum Otto cycle, the idea of local and global works was introduced by Thomas and Johal [32] and further analyzed by Altintas and Müstecaplıoglu [33]. During the two thermalization strokes, states of the subsystems were defined by the reduced density matrices of the global thermal state of the coupled working medium. Local heat exchange is defined as the average energy transferred between the subsystem and the hot or cold heat bath with respect to the local Hamiltonian of the subsystem. According to the first law of thermodynamics, local work done by a subsystem is defined as the sum of the local heat exchange with the hot and cold heat bath. This definition is valid in the weak coupling approximation between the working medium and the heat bath. For the measurement driven heat engine, we define the local work in a similar way. Let us assume that, before the measurement, total state of the system is ρ int and after the measurement it becomes ρ M . We denote ρ A int = Tr B (ρ int ) to be the reduced density matrix for the subsystem A and ρ B int = Tr A (ρ int ) to be the reduced density matrix of subsystem B, and similarly for the reduced states after the measurement. In the same way, we can calculate the reduced density matrix for the subsystems before and after the thermalization step. The local work outputs for the subsystems are defined as [32,33], q i 1 represents the average energy exchange for the subsystem A or B (i = A, B) for the measurement step after the first adiabatic expansion and q i 2 is the conventional heat exchange with the heat bath in the last step, i.e., the thermalization step. H i (B 1 ) and H i (B 2 ) are the local Hamiltonians for the subsystems (i = A, B) for external magnetic field B 1 and B 2 , respectively. After the first adiabatic expansion, the parameter of the Hamiltonian is changed from B 1 to B 2 . In the next adiabatic stroke, the magnetic field is changed back to the initial value B 1 . Total local work done by the two subsystems is w = w A + w B . Nevertheless, the validity of these definitions in the scenario of measurement driven engine is under question as, during the measurement stroke, the coupling between the working medium and the apparatus is in general not weak. Moreover, after the non-selective measurement, the working medium is driven out of equilibrium, such that the state of the whole working medium is no longer a thermal state. In Figures 15-17, we plot the local works and global works for three different spin combinations with same measurement settings considered above (S x on side A and S z on side B). The plots show very interesting behavior. Unlike the coupled quantum Otto cycle [32,33], the sum of the local works for the subsystems is not always equal to the global work. The nature of the plots also change with the change of spin values. For the two spin half case, the sum of the local works start from a negative value and goes up to zero, whereas, for the two other scenarios, it starts with a positive value. Interestingly, local work for one subsystem is zero for first two scenarios (Figures 15 and 16). Moreover, except for the two spin half scenario, local work for one subsystem exactly matches with the global work (at least for the cases considered here). As a result, global work matches with the total local work for the second scenario ( Figure 16) (black, blue and green curve have merged together). In contrast to the conventional quantum Otto cycle, the relationship between local and global works (assuming the above definition of local work) is complex. This also opens up the avenue for suitable definition of local work, taking into account the work cost for measurement, which demands a more detailed description of the engine cycle including the measuring apparatus. This may be a potential future area to explore.

Conclusions
In this paper, we investigate the effect of coupled working medium in the measurement based single temperature quantum heat engine without feedback. We consider z one-dimensional Heisenberg model of two spins and calculate the efficiency of the heat engine. We show that, when the coupling constant J is nonzero, the efficiency is increased over the uncoupled scenario, i.e., J = 0 case for certain range of non-zero J. Thus, interaction enhances the efficiency of this type of heat engine, similar to the coupled Quantum Otto cycle [32,33]. We also consider the higher-dimensional scenario, where the two spins are not only spin-1/2 but also 1 or 3/2. In these cases, we observe a very interesting situation that is absent in the conventional coupled quantum Otto engine as well as in the case of uncoupled measurement driven heat engine. When the two spins of the two subsystems are the same, we always get useful work output from the engine, and hence the efficiency is positive. However, this is not true for the asymmetric situation. If the spins for the two subsystems are different, we get negative efficiency after a certain nonzero value of the coupling constant, implying that we cannot extract work, but we have to invest work to run the cycle. It is very similar to the case of a refrigerator in the absence of cold reservoir, which is replaced here by a measurement protocol. However, as long as the efficiency is positive, for higher spin systems, efficiency is higher than the lower-dimensional system. Thus, both coupling and dimension of the Hilbert space decide the efficiency of the measurement-based heat engine. Next, we consider the local work and global work for the engine cycle. Using the same definition of local work used for quantum otto cycle, we show that global work is never the same as total local work. This is unlike the results obtained in [32,33], which may be due to inappropriate definition of local work in the measurement driven heat engine. Throughout the paper, we focuss on the quasistatic regime for the engine. In [54], the authors considered the scenario of imperfect thermalization stroke and analyzed the power of a single temperature measurement driven engine. Thus, the effect of a coupled working medium on the power of this engine might be a good candidate for subsequent study. Table A1. S A = 1/2, S B = 3/2. Table A2. S A = 1, S B = 1.