Quantum Otto heat engine based on a multiferroic chain working substance

We study a quantum Otto engine operating on the basis of a helical spin- 1/2 multiferroic chain with strongly coupled magnetic and ferroelectric order parameters. The presence of a finite spin chirality in the working substance enables steering of the cycle by an external electric field that couples to the electric polarization. We observe a direct connection between the chirality, the entanglement and the efficiency of the engine. An electric-field dependent threshold temperature is identified above which the pair correlations in the system, as quantified by the thermal entanglement, diminish. In contrast to the pair correlations, the collective many-body thermal entanglement is less sensitive to the electric field, and in the high temperature limit converges to a constant value. We also discuss the correlations between the threshold temperature of the pair entanglement, the spin chirality and the minimum of the fidelities in relation to the electric and magnetic fields. The efficiency of the quantum Otto cycle shows a saturation plateau with increasing electric field amplitude.


Introduction
With the advances in nanotechnology enabling a controlled miniaturization and a functionalization of nanostructured materials, questions related to the thermodynamical properties are gaining an increased attention.Several theoretical proposals were put forward for nanoscale Brownian motors [1], refrigerators [2] and quantum heat engines [3,4,5,6,7,8,9,10,11,12,13].On the other hand, for finite systems, the application of the laws of thermodynamics is the subject of an ongoing debate [14].One of the fundamental questions concerns the size limit to which the working substance might be scaled down.Recent studies point out that the quantum nature of a size-quantized working substance, e.g. a quantum heat engine, may lead to a close connection between the efficiency of the cycle and quantum correlations [15], which can be quantified in terms of the entanglement [16,17,18], a behavior that is atypical for classical engines.According to the fundamental laws of thermodynamics, the efficiency of a classical engine is independent of its detail and is solely determined by the character of the cycle itself and the temperatures of the heat baths.The quantum nature of the working substance, however, has key consequences for the engine output power as well.Recently it was shown that purely quantum phenomena, such as noise-induced coherence, yields greater engine output power [19,20].In general, physical phenomena at the crossover of quantum mechanics and thermodynamics are the subjects of the emergent field of quantum thermodynamics where, among other topics, questions are addressed as to what extent standard classical thermodynamic cycles, such as Carnot or Otto cycles, can be reformulated for quantum systems [4].A key issue thereby is the difference between thermodynamic and quantum adiabatic processes.For example, a thermodynamical adiabatic process does not necessarily mean that the occupation probabilities are invariant during an adiabatic transition.As usual, thermodynamical adiabatic processes are identified in terms of the conservation of the entropy and the isolation of the system from the heat exchange with the thermal bath.An essential requirement for the quantum adiabatic process is that the population distributions remain unchanged.Thus, quantum adiabaticity is a stricter requirement than the thermodynamic one.The adiabatic quantum process is also adiabatic in the thermodynamic sense, the opposite is however not true in general.Therefore, quantum adiabaticity entails a relatively low power output from a slowly operating quantum engine, unless the energy spectrum of the working substance has nodal crossing points, however, this is not a generic feature of realistic physical systems.Landau Zener transitions are avoided during an adiabatic segment of the cycle by slowly varying the control parameters [21].As mentioned above, the subtlety of quantum engines is related to the internal connection between essentially quantum phenomena such as entanglement and the thermodynamic characteristics of the cycle.In this regard, the choice of the working substance for the operating quantum engine is an important issue [22].Recently, there has been a great interest in composite multiferroic (MF) materials that possess coupled ferromagnetic (FM) and ferroelectric (FE) properties [23,24,25,26,27,28,29,30,31,32,33,34,35,36,37] (for a review we refer to [38]).These materials allow for a multitude of novel applications based on the control of magnetism (ferroelectricity) with electric (magnetic) fields.They offer new opportunities for the design and control of new circuits for quantum information processing.For an interesting class of magnetoelectrics, the ME coupling is rooted in a chiral magnetic ordering that is coupled to an electric polarization P such that [39] Above r i,i+1 , is the relative spatial vector between the effective spins σ i and σ i+1 localized at neighboring sites.Though Eq.( 1) was derived initially phenomenologically, a fully microscopic theory based on the electronic states was developed shortly thereafter [40].The emergence of an electric polarization, when coupling the spatial degrees of freedom to the spin chirality, renders possible an efficient manipulation and control of the spin order parameter via an applied external electric field.
In the present project, we will study a model for a quantum Otto engine operating on the basis of a one-dimensional (1D) finite size MF chiral spin chain that acts as a working substance.The possibility of controlling the efficiency of the quantum heat engine via an external electric field motivates our choice of the working substance.For small working substance consisting of four spins with periodic boundary conditions, we provide an analytical solution to the problem.For a larger size of the working substance an exact numerical diagonalization reveals the connection between the thermal entanglement and the cycle efficiency.
Our theoretical model is experimentally feasible.Recently discovered materials such as the quantum S = 1/2 spin chain magnets LiCu 2 O 2 , CoCr 2 O 4 , LiV CuO 4 , possess simultaneously ferroelectric and ferromagnetic properties.1D patterns of the chain magnets can be manufactured using CuO 2 powders and a P t stove [41].Then, a spin chain magnet doped on the Ir(001) or nonmagnetic Zn 2+ substrate could serve as a working substance, while a P t stove could be implemented as a thermal bath for controlling the temperature of the working substance.The cycle we are going to study consists of two thermodynamic adiabatic and two isochoric strokes.During the two isochoric strokes, the multiferroic spin chain interacts with the heat baths.During the two thermodynamic adiabatic strokes the amplitude of the electric field is changed (cf.Fig. 1 for a schematic illustration).Since the energy levels of the system depend on the electric field, a change in the electric field modifies the energy levels.In this way work is done by the engine during the thermodynamic adiabatic strokes.In what follows, we will study the connection between the cycle efficiency, the entanglement and the electric field.
The paper is organized as follows.In section II we introduce our model.In section III we present analytical results obtained for the multiferroic chain working substance consisting of four spins.In particular we study: the dependence of the pair and the nonlocal many-body entanglements on the temperature and the amplitude of the electric field and the temperature dependence of the quantum chirality and the electric and magnetic susceptibilities.Each of these quantities has a particular meaning: Quantum chirality is a measure for the spin frustration and allows to drive the cycle and to control the cycle efficiency by an external electric field.The electric and magnetic susceptibilities can be used to detect the thermal phase transitions in the working substance, and the local and many-body entanglements are useful to observe the connection between the cycle efficiency and the quantum correlations.In section IV we discuss details of the thermodynamic cycle and in section V we evaluate the scaling of the cycle efficiency with the size of the working substance.In section VI we study the quantum Otto cycle in the semiclassical limit and wrap up in section VII.We envisage the application of a MF helical chain in one spatial dimension as an electric-field controlled heat engine.For definiteness we take the x axis as the chain direction.An effective model that captures the physics of the MF chain [40,42] is based on the Hamiltonian

Model
We assumed here that the chain is subjected to an electric field (℘ = (0, ℘, 0) applied along the y axis and to a magnetic B along the z axis.The exchange interaction constant between the nearest neighbor spins is chosen ferromagnetic J 1 > 0 while the next-nearest interaction is antiferromagnetic J 2 < 0. Pauli matrices are used in standard notations σ i , and γ e is the gyromagnetic ratio for electron spin.h is Planck's constant.With the help of Eq.( 1) the coupling of the electric field to the MF chain can be written as as discussed in the literature [43,44,45,46].For LiCu 2 O 2 the values of the model parameters are J 1 ≈ 81K, J 2 ≈ 44K see [45].
The problem of an electric-field control of the magnetic chirality of a ferroaxial MF system was addressed in a recent paper [36].Information transfer by the vector spin chirality in magnetic chains was discussed in [37].The effect of the electric field, or the DM anisotropy variation on the quantum information processing as well as on many-body quantum ground states and quantum-phase transitions of MF helical chain we addressed recently in [42].
In what follows, we suppose that J 1 = −J 2 = J and go over to dimensionless units such that B → γ e hB/J, ℘ → g M E | − → ℘ |/J, i.e., we measure the Zeeman energy and the interaction energy with electric field in units of the exchange constant.As was mentioned in the introduction, the purpose of the present project is to investigate a possible control of the cycle efficiency and thermal entanglement via an external electric driving field.For clarity we combine analytical and full numerical approaches and start with a solvable model consisting of a four spins with periodic boundary conditions as a working substance.

Four spins case: Analytical treatment
In the case of four spins, the Hamiltonian (2) can be diagonalized analytically.All technical details are provided as supplementary materials to this paper.Here we present the main results.To study the thermal entanglement and the cycle efficiency we construct the density matrix ρ corresponding to an equilibrium Gibbs distribution, where |ψ n and E n are the eigenfunctions and eigenvalues of (2) given explicitly in the supplementary materials (See (A1), (A2)).Using the density matrix (3) we calculate the mean value of the z component of the VC.The only nonzero component for the considered configuration of the system and the chosen direction of the electric field is [17] we calculate the pair concurrence between two arbitrary spins of the working substance C nm = max(0, R Here R (α)  nm , α = 1, 2, 3, 4 are the eigenvalues of the matrix 2 ) and ρ R nm is the reduced density matrix of the system of two spins obtained from the density matrix of the system ρ (3) after tracing out two remaining spins ρ R nm = tr sp (ρ), where s, p = m, n.After some rather straightforward calculations we obtain Similarly we calculate for the other components In view of Eqs. ( 4)-( 6) for the different pair concurrences we infer Explicit expressions of the parameters a 1,2 , b 1,2 , c 1,2 , d 1,2 are quite involved and therefore are presented in the supplementary materials (see Eqs. (A3), (A4)).The pair concurrences C nm depend on the chosen spins n and m (more precisely, due to the translational symmetry of periodic chain, on the distance between the two spins) and as one can see from Eq.( 7), C nm are quite different from each other.Therefore, the more informative and universal quantity seems to be the two tangle [47] τ 2 which contains information on the total pair correlations in the spin chain τ 2 = 2C 2  12 + C 2 13 .Taking into account Eq.( 7) for the two tangle we deduce The degree of the pair correlations depends thus on several inequalities between the parameters a 1,2 , b 1,2 , c 1 , d 1,2 which are functions of the temperature T , and the amplitudes of the driving electric and magnetic fields ℘, B. The explicit expressions of the parameters entering in Eqs.( 8) can be found in the supplementary materials as Eqs.(A3) and (A4).In effect, depending on the values of the three parameters ℘, B, T the system can be entangled or disentangled.The threshold temperature T c (℘, B) for the given amplitudes of the driving fields defines the regimes of entangled and disentangled states.Our principle interest now is to see how the threshold temperature of the system scales with the electric field T c (℘).
Another interesting object is the one tangle [47] τ 1 = 4detρ 1 as it quantifies the nonlocal many-body correlations in the spin chain.Here ρ1 = tr 2,3,4 (ρ) is the reduced density matrix of the first spin after tracing out the states of all other spins.Explicit expressions for the one tangle can be found in the analytical form τ 1 = 4 Z 2 Q, where the explicit form of Q is presented in the supplementary materials (see (A5)).One tangle and the collective nonlocal entanglement is related to the complex spiral spin structure of the frustrated MF chain and therefore is generated by the external electric field.An indication of the existence of the spiral spin structure in the system is the nonzero chirality.Therefore, the values of the z component of the vector chirality tr(ρ 4 i=1 [ê x × (σ i × σ i+1 )] y ) and the values of the one tangle τ 1 = 4 Z 2 Q are in correlation with each other.Both of them depend on the amplitude of the electric field.Based on the definition of the parameter Q (supplementary materials Eq.(A5)) it is easy to see that in the high temperature limit T → ∞ for one tangle we obtain τ 1 = 1.Considering the eigenvalues and eigenfunctions of the Hamiltonian (2) (see supplementary materials Eqs.(A1) and (A2)) for the chirality of the MF chain we obtain: The magnetic field dependence of both types of the entanglement, namely the shortrange pair correlations as quantified by τ 2 and the many-body collective entanglement described by τ 1 , is transparent (this is due to the fact that in our model the magnetic field couples to the magnetization which is a conserved quantity in our model).With increasing B the entanglements decrease.The dependence on the amplitude of the electric field is less obvious and deserves a detailed consideration.First we focus on the thermal pair entanglement τ 2 .As we see in Fig. 2, τ 2 is finite only for a very large amplitude of the electric field.The pair thermal concurrence τ 2 is practically zero until the electric field amplitude becomes quite substantial as can be seen in Fig. 2 (recall we are operating in scaled units).We also determine the threshold temperature below which τ 2 is finite and above which τ 2 = 0, see Fig. 3.In contrast to the pair correlations and the entanglement, the collective entanglement τ 1 is different from zero for an arbitrary electric field.(see Fig. 4).Another remarkable difference is that the collective entanglement τ 1 is very robust and is practically not effected by the temperature (see Fig. 5).Therefore, the amount of the thermal entanglement stored in the nonlocal correlations τ 1 is always larger than the thermal entanglement of the pair correlations τ 1 > τ 2 , as shown in Fig. 6.With increasing the size of the working substance the ratio between τ 2 /τ 1 becomes smaller.This means that the many-body entanglement τ 1 is increasing with the size of the system N faster than the total two pair correlations τ 2 .The situation with respect to the thermal chirality is different.In particular, we observe that with the increase of the temperature the thermal chirality undergoes a strong change and above the threshold temperature T c of the two tangle τ 2 , the thermal chirality is almost zero, as depicted in Fig. 7. Thus, we see that the thermal chirality is correlated with the pair correlations in the system.As for the temperature, the magnetic field also has a negative influence on the chirality.The dependence of the thermal chirality on B is plotted in Fig. 8, where one can see that with increasing B, the thermal chirality decreases.Another quantity of interest that quantifies the sensitivity to perturbations of quantum systems near a critical region is the fidelity [48].A zero temperature fidelity is a measure for the overlap between two ground states corresponding to slightly different values of the controlling parameters.A dip in fidelity reflects changes in the structure of the ground state at a quantum critical point [48].The finite temperature thermal state extension of the quantum fidelity was considered in [49].The fidelity of a mixed state at a finite temperature characterizes a second-order thermal phase transition [5] and is defined in the following way where ρ0 (β, ζ 0 ), ρ1 (β, ζ 1 ) are the density matrixes of the system corresponding to slightly different control parameters The expression for the fidelity related to the electric and magnetic fields can be simplified to the following form where χ(ζ) = − ∂ 2 F ∂ζ 2 is the susceptibility to the corresponding external field ζ = ℘, B at constant temperature.Analytical expressions of the susceptibilities are presented in the supplementary materials (Eq.(A6)-(A7)).
Finally, the dependence of the electric χ(℘) and the magnetic χ(B) susceptibilities on the temperature are depicted in Fig. 9, Fig. 10, Fig. 11 and Fig. 12.We choose B = 1.The peak in the susceptibility shifts towards lower temperatures as we increase the system size.A similar peak in the magnetic susceptibility at finite temperatures was observed for B = ℘ = 0 and N = 24 in [50] and was interpreted as the result of a competition between antiferromagnetic and ferromagnetic correlations in the system.According to the definition Eq.( 11), the maxima of the susceptibilities correspond to the minima in the fidelities that are related to the electric and magnetic fields.Comparing Fig. 9 -Fig.12 with Fig. 3, we see a direct correlation between the threshold temperature of the pair entanglement τ 2 and the minima of the fidelities related to the electric and magnetic fields.The maxima of the electric and magnetic susceptibilities are related to threshold temperatures of the pair correlations.For larger threshold temperatures of the pair correlations, the maxima of the electric and magnetic susceptibilities are shifted towards higher temperatures.Interestingly for the electric susceptibility, the correlation between the threshold temperature of the pair entanglement τ 2 and the minimum of the fidelity is not only qualitative but also quantitative as well.As we see for large enough electric fields the maximum of the electric susceptibility is observed almost on the threshold temperatures T c ≈ 24, T c ≈ 45 of the pair entanglement τ 2 .
In Figs. 10 and 12 we present the system size dependence of the electric and the magnetic susceptibilities.The heights of the peaks of the electric and the magnetic susceptibilities increase with the system size.One can as well observe that for N = 8 the location of the peak of the magnetic (electric) susceptibility shifts towards lower (higher) temperatures.

Efficiency of the multiferroic heat Otto engine
In analogy to the classical Otto cycle, the quantum Otto cycle also consists of two quantum isochoric and two adiabatic processes [4], as sketched in Fig. 1.The quantum isochoric process corresponds to a heat exchange between the working substance and cold and hot heat baths.During the quantum isochoric process only level populations P n are reshuffled, while during the adiabatic process the working substance produces work and in this case the energy levels are changed.Therefore, the work produced by the engine depends on the amplitude of electric field which causes changes in energy levels.An adiabatic process can be thermodynamic adiabatic or quantum adiabatic.A process is thermodynamic adiabatic if the working substance is thermally isolated from the heat exchange with the heat bath.However, this does not exclude inter-level transitions of a purely quantum nature, while in the case of a quantum adiabatic process the level populations are fixed.
As was mentioned in introduction, in our case the working substance is a MF spin frustrated chain with a discrete energy spectrum of 16 levels.The first law of thermodynamics for a system with discrete energy spectrum reads Here dU is the change of the system energy U(E n , T ) = tr (ρ Ĥ) = 16 n=1 E n P n (E n , T ).The first term on the right hand side of Eq. ( 12), δQ = E n dP n can be viewed as the heat exchange and is related to the change of the level populations P n (E n , T ) occurring due to a change of the temperature for E n = const, while the second and the third terms correspond to the produced work.If the adiabatic strokes of the cycle are quantum adiabatic then ∂Pn ∂En T =const = 0 and Eq.( 12) reduces to the form given in [4].The work produced during the quantum adiabatic process reads δW = P n dE n .The working substance produces work due to the change of the amplitude of electric field ℘.This leads to a modification in the energy levels with ∆E n = E n (℘) − E n (℘ 1 ).Our goal is so to study the dependence of the cycle efficiency on the modulation of the control parameter, i.e., the electric field amplitude ℘.
To this end we considered two slightly different quantum Otto cycles.As shown in Fig. 1, the first cycle consists of four strokes [22]: Step A → B and C → D are two isochoric processes.During the step A → B the system couples to the hot bath at temperature T H and the energy levels are unchanged.After absorbing energy from the hot bath, the system reaches a thermodynamic equilibrium state, which can be described by the level populations Step C → D is the reverse of the step A → B. Namely, the system is brought to couple to a sink at the temperature T L .After energy exchange with the heat bath a thermodynamic equilibrium state is established, which can be described by the level populations P D n (E n (℘ 1 ), T L ).Steps B → C and D → A are quantum adiabatic processes in which the level populations are unchanged, i.e.P A n = P D n , P B n = P C n .During the process B → C amplitude of the electric field is changed ∆E n = E n (℘) − E n (℘ 1 ).The working substance performs a positive work.Therefore, the heat absorbed by the working substance and the heat released read [22] In the second scenario for the cycle, the quantum adiabatic strokes of the cycle are replaced by thermodynamic adiabatic strokes.The heat absorbed by the working substance Q in and the heat released in the quantum isochoric cooling process Q out in the case of the thermodynamic adiabatic cycle are defined in the following form: while the expression for the cycle efficiency reads The dependence of the quantum Otto cycle efficiency on the modulation of the electric field amplitude is presented in Fig. 13.For both types of the cycles (quantum adiabatic and thermodynamic adiabatic) we observed qualitatively similar dependencies on the electric field.In both cases the maximal efficiency is reached for certain optimal values of the modulation of the electric field amplitude.However, the maximal efficiency obtained for thermodynamic adiabatic cycle is higher compared to the efficiency corresponding to the quantum adiabatic case.Let us concentrate on the thermodynamic adiabatic cycle.As one can see reasonably high efficiency of around 75% is achievable already for ℘/℘ 1 ≈ 5. We also observe a saturation of the cycle efficiency with a further increase of the electric field amplitude.Depending on the amplitude of electric field, the efficiency of the quantum Otto engine might be higher or lower as compared to the maximal efficiency of the Carnot cycle η C = 1 − T L /T H .The reason why the efficiency of the quantum cycle exceeds the maximal efficiency of the Carnot cycle is of an entirely quantum origin, as it can be traced back to the entanglement of the working substance [15].To illustrate this, we plotted the efficiency of the cycle as a function of the entanglement.Fig. 14 evidences that an increase of the entanglement in the system result in an enhanced cycle efficiency.The efficiency of the quantum Otto cycle as a function of the pair entanglement, for the following parameters: ℘ 1 = 3.5, B = 1.We find that a strong entanglement in the system is related to an enhanced cycle efficiency.

Scaling of the cycle efficiency with the size of working substance
In order to investigate the scaling of the Otto cycle efficiency with the size of the working substance for different values of the electric field we plot the dependence of the cycle efficiency on the length of the spin frustrated MF chain N, as shown in Fig. 15.The efficiency for ℘ 1 = 3.5 and ℘ = 4, ℘ = 5, ℘ = 10 is plotted as a function of system size N .We find that the efficiency jumps abruptly for N = 3.For N > 4 no significant variation in the efficiency is observed when increasing N .This is we observed for all cases.
The case with a smaller size of the working substance shows a different behavior in comparison to a large system size.An abrupt increase in the efficiency is found for N = 3 .For N > 4 almost the efficiency hardly change with a larger system size.For N > 4 and a very high electric field ℘ = 10, we encounter the efficiency corresponding to the saturation value.

Four spins semi-classical limit
To conclude our analytical considerations we inspect the semi-classical limit utilizing the canonical thermodynamic perturbation theory [51].We proceed by writing for the Hamiltonian (2) Here again ℘ is assumed to be the small parameter.The eigenvalues and the eigenfunction of Ĥ0 are denoted by E 0 n , |φ n .We will utilize the normalization condition is the free energy and the left side of the normalization condition is a function of temperature.Taking the derivative of the normalization condition and after straightforward algebraic manipulations we obtain On the other hand, using the canonical thermodynamic perturbation theory [51] in the semi-classical, high temperature limit E 0 n − E 0 m < T we deduce Here V nn (℘) = n P n V nn (℘) is the mean value of the matrix element of the perturbation V evaluated in the basis of the unperturbed Hamiltonian Ĥ0 and The level populations are described in terms of the Gibbs distribution function P n .For four spins according to (17) we have Where For other energy values of the unperturbed Hamiltonian Ĥ0 we refer to the supplementary materials.
Combining (16), (17) one can infer that the change in the entropy during the heating and cooling processes depends not only on the initial and the final temperature but also on the electric field.In detail one finds S(T, ℘) = ln We remark that the entropy is defined by the partial derivative of the free energy with respect to the temperature at constant values of the electric field (16).Therefore, the dependence of the entropy on the electric field is parametrical.If the temperature is constant the entropy is constant, as well.However, the change in the entropy due to a change of the temperature depends on the values of the electric field.Eq. (19) tells that this dependence is quadratic in the field S(T, ℘) − S(T, ℘ = 0) ≈ ℘ 2 (• • •).The entropy as a function of the temperature for different values of the parameters of electric field is plotted at the Fig. 16.The maximum of the entropy is observed for small values of the electric field and in the high temperature limit.We note that Eq.( 19) is obtained via a thermodynamic perturbation theory and negative values of the entropy correspond to values of the parameters beyond the range where perturbation theory is viable.Taking into account Eq.( 19) we can express the semi-classical efficiency in terms of the electric field as The semiclassical efficiency as a function of the electric field ℘ and the temperature difference between the hot and the cold baths ∆T = T H − T L for T L = 100 are presented in Fig. 17  The results exhibit a larger sensitivity of the semiclassical efficiency to electric field ℘ than to temperature variations ∆T .

Conclusions
In the present project we studied a quantum Otto engine operating with a working substance consisting of electrically controlled multiferroic spin chain.We have shown that due to the existence of a nonzero spin chirality coupled to an emergent electric polarization in the working substance, the efficiency of the cycle is sensitive to the applied external electric field.We analyzed the dependence of the cycle efficiency on the size of the working substance.In particular, for a small working substance consisting of N = 3 spins, the efficiency reaches considerably high values, i.e., slightly below 100%.With increasing the size of the working substance, the efficiency of the quantum Otto cycle, shows a saturation plateau.Another interesting finding is the robustness of the many-body collective entanglement to temperature variations.Thereby, the many-body entanglement is quantified in terms of one tangle τ 1 and is always larger than the total pair concurrence, as described by the two tangle τ 2 .This indicates that a major amount of the entanglement of the MF working substance is stored in the long-range, multi-spin correlations (see Fig. 6).In contrast to the one tangle, the pair correlation is sensitive to the increase in the temperature.In particular, we observe the existence of a threshold temperature for the two tangle τ 2 .The stronger the electric field amplitude, the higher is the threshold temperature, as illustrated in Fig. 3.The same behavior we encounter for the chirality, as well (cf.Fig. 7).In particular, with increasing the temperature the thermal chirality undergoes strong changes and for the threshold temperature T C of the two tangle τ 2 , the thermal chirality reaches its minima.Thus, the thermal chirality is related rather to the pair correlations in the system, not to the nonlocal entanglement.We also studied the relations of the magnetic and the electric susceptibilities to the temperature.According to the definition Eq.( 11), the maximum of the susceptibilities corresponds to a minimum of the fidelity.Comparing Fig. 9 -Fig.12 with Fig. 3, we conclude on a direct correlation between the threshold temperature of the pair entanglement τ 2 and the minimum of the fidelities.The maximum of the electric and the magnetic susceptibilities are related to the corresponding threshold temperatures of the pair correlations.The maximum of the electric susceptibility is observed almost at the threshold temperatures T c ≈ 24, T c ≈ 45 of the pair entanglement τ 2 .Indeed, the quantum engine with a MF working substance is much more sensitive to the electric field than to the magnetic field which we think is favorable from an experimental point of view.
Supplementary Materials to: Quantum Otto heat engine based on a multiferroic chain working substance

A. Eigenfunctions and eigenvalues
Eigenfunctions and corresponding eigenvalues of the Hamiltonian (2) for four spins are,  Thereby the following notations are used .
Note, since the magnetic field couples to a conserved quantity, the eigenfunctions do not depend on the magnetic field amplitude B.
The expressions for the parameters entering in the (2)-( 7) read,

Figure 1 .
Figure 1.(Color online) A schematic of the considered quantum Otto cycles based on a chiral multiferroic chain.The cycle has four strokes: Step A → B and C → D are two isochoric processes.During the step A → B the system is attached to a hot bath with a temperature T H . Step C → D is inverse to the step A → B. After releasing (absorbing) energy to the cold (hot) heat bath with temperature T L (T H ) the system reaches a thermodynamic equilibrium state associated with the level populations P D n E n (℘ 1 ), T L (P B n E n (℘), T H . Steps B → C and D → A are two thermodynamic adiabatic processes.During the process B → C the amplitude of the electric field is changed so that ∆E n = E n (℘) − E n (℘ 1 ) and the working substance performs a positive work.

Figure 2 .
Figure 2. (Color online) Dependence of the two tangle τ 2 (℘) on the electric field for B = 1, and for three different Temperatures: red cross line T = 10, blue solid line T = 20, black triangular line T = 30.

Figure 3 .
Figure 3. (Color online) Dependence of the two tangle τ 2 (T ) on the temperature for B = 1, and for fixed values of the electric field: blue solid line ℘ = 1, black triangular line ℘ = 10, red cross line ℘ = 20 The threshold temperatures are T c = 7.37, T c = 22.31, T c = 44.45,respectively.

Figure 4 .
Figure 4. (Color online) Dependence of the one tangle τ 1 (℘) on the electric field for B = 1, and for three different temperatures: red cross line T = 10, blue solid line T = 20, black triangular line T = 30.

Figure 5 .
Figure 5. (Color online) Dependence of the one tangle on the temperature for B = 1, and fixed values of the electric field: black solid line ℘ = 1, blue triangular line ℘ = 5, Red cross line ℘ = 10.

Figure 6 .
Figure 6.(Color online) Ratio between two and one tangle τ 2 /τ 1 as a function of external electric field for the following values of the parameters B = 1, T = 7.37.Result for N = 4 spins is plotted using obtained analytical solutions.Result for N = 8 spins is plotted using numerical solutions.

Figure 7 .
Figure 7. (Color online) Chirality as a function of the temperature, for B = 1, and for the fixed values of the electric field: cyan solid line ℘ = 1, red cross line ℘ = 1.5.

Figure 8 .
Figure 8. (Color online) Chirality as a function of the magnetic field for the following values of the parameters ℘ = 1, T = 20.

Figure 9 .
Figure 9. (Color online) Magnetic susceptibility as a function of the temperature for B = 1 and for the fixed values of the electric field: blue solid line ℘ = 1, black triangular line ℘ = 10, red cross line ℘ = 20.We see that with the increasing of the electric field amplitude the maximum of the magnetic susceptibility is shifted towards higher temperatures.

Figure 10 .
Figure 10.(Color online) The magnetic susceptibility is plotted as a function of the temperature for N = 4 and N = 8 cases.We choose B = 1.The peak in the susceptibility shifts towards lower temperatures as we increase the system size.A similar peak in the magnetic susceptibility at finite temperatures was observed for B = ℘ = 0 and N = 24 in[50] and was interpreted as the result of a competition between antiferromagnetic and ferromagnetic correlations in the system.

Figure 11 .
Figure 11.(Color online) The electric susceptibility as a function of the temperature for the following values of the parameters B = 1 and for fixed values of the electric field: blue solid line ℘ = 1, black triangular line ℘ = 10, red cross line ℘ = 20.Comparing this result to Fig.4 we see that with the increase of the threshold temperature of pair correlations τ 2 the maximum of the electric susceptibility drifts towards higher temperatures.T c ≈ 24, T c ≈ 45.

Figure 12 .
Figure 12. (Color online) The electric susceptibility as a function of the temperature for N = 4 and N = 8 cases.we set B = 1.Note, the electric susceptibility increases as the system size increases.

Figure 13 .
Figure 13.(Color online) The efficiency of the quantum Otto cycle as a function of the modulation of the electric field amplitude, for the following values of the parameters B = 1, T H = 30, T L = 10.The inset corresponds to the quantum adiabatic case.As evident from the figure in the case of a thermodynamic adiabatic quantum Otto cycle, depending on the values of the electric field, the maximal efficiency reaches the value η = 0.75 which is higher than the maximal efficiency of the Carnot cycle η C = 1 − T L /T H = 0.66.The efficiency of the quantum adiabatic Otto cycle (inset) is slightly below the efficiency of the thermodynamic adiabatic Otto cycle.

Figure 14 .
Figure 14.(Color online) The efficiency of the quantum Otto cycle as a function of the pair entanglement, for the following parameters: ℘ 1 = 3.5, B = 1.We find that a strong entanglement in the system is related to an enhanced cycle efficiency.

Figure 15 .
Figure 15.(Color online) The efficiency for ℘ 1 = 3.5 and ℘ = 4, ℘ = 5, ℘ = 10 is plotted as a function of system size N .We find that the efficiency jumps abruptly for N = 3.For N > 4 no significant variation in the efficiency is observed when increasing N .This is we observed for all cases.

Figure 16 .
Figure 16.(Color online) Contour plot of entropy as a function of the modulation of the electric field amplitude and the temperature. .

Figure 17 .
Figure 17.(Color online) The semiclassical efficiency λ as a function of the electric field ℘ and the temperature difference between the hot and the cold baths ∆T = T H − T L for T L = 100, ℘ 1 = 0.5 plotted using(20).The semiclassical efficiency is more sensitive to the values of the electric field than to the temperature difference ∆T .