Quantum Heat Engines with Singular Interactions

: By harnessing quantum phenomena, quantum devices have the potential to outperform their classical counterparts. Here, we examine using wave function symmetry as a resource to enhance the performance of a quantum Otto engine. Previous work has shown that a bosonic working medium can yield better performance than a fermionic medium. We expand upon this work by incorporating a singular interaction that allows the effective symmetry to be tuned between the bosonic and fermionic limits. In this framework, the particles can be treated as anyons subject to Haldane’s generalized exclusion statistics. Solving the dynamics analytically using the framework of “statistical anyons”, we explore the interplay between interparticle interactions and wave function symmetry on engine performance. powered by energetic coherence.


Introduction
Thermodynamics was originally developed as a physical theory for the purpose of optimizing the performance of large-scale devices, namely steam engines [1]. Despite these practically-focused origins, thermodynamics has proven enormously successful in formulating universal statements, such as the second law. With the growing prominence of quantum technologies, the field of quantum thermodynamics has emerged to understand how the framework of thermodynamics can be extended to quantum systems [2]. One of the principal goals of quantum thermodynamics is to discover how quantum features, such as entanglement, superposition, and coherence, can be best leveraged to optimize the performance of quantum devices.
A notable feature of quantum particles is that they are truly indistinguishable. To account for this, the wave function of a multiparticle state is constructed out of the symmetric (for bosons) or antisymmetric (for fermions) superposition of the single particle states. This wave function symmetry has physical consequences in terms of state occupancy, with any number of bosons allowed to occupy the same quantum state while fermions are restricted to single occupancy-the famous Pauli exclusion principle. The superposition arising from the symmetrization requirements also leads to interference effects that manifest as "exchange forces" in the form of an effective attraction between bosons and effective repulsion between fermions [53]. Wave function symmetry leads to quantum modifications of thermodynamic behavior, allowing for more work to be extracted from indistinguishable particles through mixing [54], or as the working medium of a quantum heat engine [24,[55][56][57]. In [57], we showed that for a working medium of two non-interacting identical particles, a harmonic quantum Otto engine exhibits enhanced performance if the particles are bosons and reduced performance if they are fermions. In this paper we expand upon these results by introducing an interaction proportional to the inverse square of the interparticle distance in addition to the standard harmonic potential. This model is often referred to as the singular [58] or isotonic oscillator [59].
This potential is of particular interest as it provides the basis for the Calogero-Sutherland model [60,61], a system known to host generalized exclusion statistics (GES) anyons [62,63]. In the framework of GES, the Pauli exclusion principle is generalized to allow for a continuum of state occupancy, from the single occupancy allowed to fermions up to the infinite state occupancy allowed to bosons [62]. In the Calogero-Sutherland model this anyonic behavior arises from tuning the interparticle interaction strength, effectively interpolating between the bosonic and fermionic exchange forces [63].
We determine a closed form for the thermal state density matrix for two particles in the singular oscillator potential. By mapping to an equivalent system of "statistical anyons" [24] consisting of a statistical mixture of bosons and fermions, we determine the time-dependent internal energy while varying the oscillator potential. We then analyze a quantum Otto cycle, demonstrating performance that interpolates between the bosonic and fermionic of [57] as the interaction strength parameter is changed.

The Two-Particle Singular Oscillator
We begin by outlining the model, notation, and previous results from the literature that are central to our analysis. We use the following Hamiltonian for two particles in a singular oscillator potential [63,64], where the parameter ν quantifies the strength of the interparticle interaction. Following Calogero's original approach [60], we rewrite this Hamiltonian in the center of mass and relative coordinate frame using the momentum coordinate transformations P = p 1 + p 1 , p = (p 1 − p 1 )/2 and the position coordinate transformations X = ( With this transformation, we can solve the dynamics of the center of mass and relative coordinates separately, where M = 2m and µ = m/2. We see that the center of mass Hamilton is identical to that of an unperturbed harmonic oscillator for a single particle of mass 2m, while the relative motion Hamiltonian is identical to that of a singular oscillator with a single particle of mass m/2. Note that this approach is very commonly applied when solving the dynamics of classical interacting oscillator systems [65]. The eigenfunctions and eigenenergies of Equation (2) are well known, and those of Equation (3) have been determined directly [59,60], using operator methods [66] and as a generalization of the Morse potential [58]. The eigenfunctions are, with the corresponding eigenenergies, where n = 0, 1, 2, ... and L α n (x) is the generalized Laguerre polynomial [67]. For indistinguishable quantum particles, the total wave function must remain symmetric under particle exchange for bosons, and antisymmetric for fermions. In the center of mass and relative coordinate framework, the symmetry condition is satisfied by ψ rel (x) = ψ rel (−x) for bosons and ψ rel (x) = −ψ rel (−x) for fermions [61]. Examining Equation (4), we see that the parity depends on the value of the interaction strength parameter ν. Noting the relationship between the generalized Laguerre and Hermite polynomials [68], along with the Gamma function identity [67], we see that for ν = 0, Equation (4) becomes, Similarly, for ν = 1, Equation (4) becomes, For ν = 0 and ν = 1, we see that the interaction potential vanishes, reducing Equation (1) to the Hamiltonian of two particles in a pure harmonic potential. We note that Equations (8) and (9) correspond to the eigenfunctions of the relative motion Hamiltonian for two bosons and two fermions in a pure harmonic potential, respectively. The restriction of Equation (8) to even and Equation (9) to odd Hermite polynomials ensures that the proper exchange symmetry conditions are satisfied. This demonstrates that interacting particles in the singular oscillator potential can be treated as noninteracting anyons in a harmonic potential obeying generalized exclusion statistics [62], with ν as the parameter that controls the nature of the particle statistics. This behavior was first established by Murthy and Shankar in the context of the thermodynamics of the Calogero-Sutherland model [63].

Singular Oscillator Thermal State
The equilibrium thermodynamic behavior of the two-particle singular oscillator can be determined from the canonical partition function, Noting H = H CM + H rel , the partition function can be split into the product of individual partition functions for the center of mass and relative motion. The center of mass partition function can be found straightforwardly in the energy representation using the harmonic oscillator eigenenergies, We can find the relative partition function using Equation (5), The total canonical partition function is then, Equation (13) can be identified as a product of bosonic and fermionic contributions. To this end, note that the partition functions for two bosons and two fermions in a harmonic potential are [24], and, Therefore we can express the singular oscillator partition function (13) as, This agrees with the partition function determined by Murthy and Shankar for the Calogero-Sutherland model [63] and that of a statistical mixture of bosons and fermions in a harmonic potential, using the framework of statistical anyons [24].
A fuller thermodynamic picture arises from the equilibrium thermal density matrix. This can be determined in position representation using, The density matrix for the center of mass motion is given by the known thermal state of the quantum harmonic oscillator [69], The density matrix for the relative motion can be found in closed form by combining Equations (4) and (17) and applying the Hardy-Hille formula [70], where I α (x) is the modified Bessel function of the first kind [67]. This yields, The total thermal state density matrix in the position representation can then be determined from the simple Cartesian product of ρ CM and ρ rel , Equations (18) and (20), respectively,

Singular Oscillator Dynamics
With the equilibrium behavior of the two-particle singular oscillator established, we next consider the evolution of the system for a time-dependent parameterization of the frequency, ω = ω(t). As in the case of the thermal state, we will treat the dynamics of the relative and center of mass motion separately.
The dynamics of the center of mass coordinate, described by Equation (2), are that of the well-studied time-dependent parametric harmonic oscillator [71]. The time-evolved state can be determined using the path integral formulation by applying the appropriate propagator, where U CM (X, X 0 ; t) is given by [71], Z t and W t are time-dependent solutions to the classical harmonic oscillator equation of motion,Z with initial conditions Z 0 = 0,Ż 0 = 1 and W 0 = 1,Ẇ 0 = 0. The dynamics of the relative coordinate, described by Equation (3), are that of the time-dependent parametric singular oscillator. The corresponding propagator for this Hamiltonian has also been determined exactly [72][73][74][75], where J ν (x) denotes the Bessel function of the first kind [67] and Z t and W t are the same as in Equation (24). As in the center of mass case, the time-evolved density matrix for the relative motion is given by, While the integrals in Equation (22) can be carried out analytically [76][77][78], the product of Bessel functions in Equation (26) makes determining a closed-form expression for general values of ν difficult. To determine an analytical expression for the time-evolved state, we instead turn to the framework of statistical anyons [24].
In [24], we showed that the behavior of a pair of particles in a non-interacting statistical mixture of boson and fermion pairs is, on average, fully equivalent to that of interacting particles obeying generalized exclusion statistics. In this framework, the generalized exclusion statistics parameter is equivalent to the probability that a given pair of particles in the statistical mixture are fermions, p F [24]. In the context of this work, this means that we can exactly map the behavior of particles in a singular oscillator potential to the average behavior of a statistical mixture of bosons and fermions in a harmonic oscillator potential. Using this mapping, the singular oscillator density matrix is given by, where ρ har (x 1 , x 2 , y 1 , y 2 ; t) and ρ har (x 1 , x 2 , y 1 , y 2 ; t) are the time-evolved states for two noninteracting bosons and fermions in a harmonic potential, respectively. The full expressions were determined in [57] and are provided in Appendix A for completeness. Note that in Equation (22), we have switched back into the individual particle frame. Recalling that the generalized exclusion statistics parameter for the singular oscillator is given by the interaction strength [63], we have the relation p F = ν.
The time-dependent average energy can be determined in the usual manner, Using Equation (27), this becomes, In [57], we found H and, where ω 0 is the frequency at t = 0 and Q * is a dimensionless parameter that measures the degree of adiabaticity of the evolution [71], For a perfectly adiabatic stroke, Q * = 1, and in general, Q * ≥ 1 [42,71]. By combining Equations (29)-(31), we find, Note that H t = (ω t /ω 0 )Q * H 0 , where H 0 is the internal energy of the equilibrium thermal state, which can be determined from Equation (13) using, This expression for H t agrees with the result of [55], which determined the timedependent average energy of the singular oscillator using the scale invariance of the potential.

The Quantum Otto Cycle
The quantum Otto cycle, like its classical counterpart, consists of four strokes: (1) isentropic compression, (2) isochoric heating, (3) isentropic expansion, and (4) isochoric cooling. For our working medium of two particles confined within a singular oscillator, the isentropic strokes consist of increasing the oscillator frequency for compression and decreasing it for expansion. The cycle is illustrated graphically in Figure 1. During an isentropic stroke, the state of the working medium evolves unitarily with constant von Neumann entropy, thus fulfilling the isentropic condition. This is notably different from the classical Otto cycle in which the isentropic strokes are thermodynamically adiabatic and reversible. For the quantum cycle, the finite time unitary strokes are not generally adiabatic (i.e., in accordance with the quantum adiabatic theorem) and thus subject to "quantum friction" in the form of nonadiabatic excitations to other energy states [2].
The isochoric strokes consist of coupling the working medium to a hot (cold) bath to increase (decrease) the energy of the system while holding the oscillator frequency constant. We make the standard assumption that the thermalization time is sufficiently short such that the working medium has achieved a state of thermal equilibrium with the bath by the end of each isochoric strokes, which removes the need to explicitly model the system-bath interaction [12,15,24,29,32,42,57,79,80]. Note that, unlike the classical quasistatic Otto cycle, the quantum Otto cycle is fundamentally irreversible [2]. The working medium is only in a state of thermal equilibrium with the cold and hot baths at points A and C, respectively. The nonadiabatic nature of the isentropic strokes will drive the system away from these equilibrium states during the rest of the cycle. The thermalization of the resulting out-ofequilibrium states during the isochoric strokes is then the source of the irreversibility of the cycle [2].      z v r z K w Q Q n 7 C i w d F v P o 7 3 v w b J 8 k e N L G g o a j q p r s r S g U 3 1 v e / v c L a + s b m V n G 7 t L O 7 t 3 9 Q P j x q G p V p h g 2 m h N L t i B o U P M G G 5 V Z g O 9 V I Z S S w F Y 1 u Z 3 7 r C b X h K n m w 4 x R D S Q c J j z m j 1 k n t r p I 4 o L 2 g V 6 7 4 V X 8 O s k q C n F Q g R 7 1 X / u r 2 F c s k J p Y J a k w n 8 F M b T q i 2 n A m c l r q Z w Z S y E R 1 g x 9 G E S j T h Z H 7 v l J w 5 p U 9 i p V 0 l l s z V 3 x M T K o 0 Z y 8 h 1 S m q H Z t m b i f 9 5 n c z G 1 + G E J 2 l m M W G L R X E m i F V k 9 j z p c 4 3 M i r E j l G n u b i V s S D V l 1 k V U c i E E y y + v k u Z F N f C r w f 1 l p X a T x 1 G E E z i F c w j g C m p w B 3 V o A A M B z / A K b 9 6 j 9 + K 9 e x + L 1 o K X z x z D H 3 i f P 7 m H j 7 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " F 9 / v h L s Q P x m T f 6 O 6 t s Y L U D F 6 1 y 0 = " > A A A B 7 3 i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K o M e g F 4 8 R z A O S J c x O e p M h M z v r z K w Q Q n 7 C i w d F v P o 7 3 v w b J 8 k e N L G g o a j q p r s r S g U 3 1 v e / v c L a + s b m V n G 7 t L O 7 t 3 9 Q P j x q G p V p h g 2 m h N L t i B o U P M G G 5 V Z g O 9 V I Z S S w F Y 1 u Z 3 7 r C b X h K n m w 4 x R D S Q c J j z m j 1 k n t r p I 4 o L 2 g V 6 7 4 V X 8 O s k q C n F Q g R 7 1 X / u r 2 F c s k J p Y J a k w n 8 F M b T q i 2 n A m c l r q Z w Z S y E R 1 g x 9 G E S j T h Z H 7 v l J w 5 p U 9 i p V 0 l l s z V 3 x M T K o 0 Z y 8 h 1 S m q H Z t m b i f 9 5 n c z G 1 + G E J 2 l m M W G L R X E m i F V k 9 j z p c 4 3 M i r E j l G n u b i V s S D V l 1 k V U c i E E y y + v k u Z F N f C r w f 1 l p X a T x 1 G E E z i F c w j g C m p w B 3 V o A A M B z / A K b 9 6 j 9 + K 9 e x + L 1 o K X z x z D H 3 i f P 7 m H j 7 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " F 9 / v h L s Q P x m T f 6 O 6 t s Y L U D F 6 1 y 0 = " > A A A B 7 3 i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K o M e g F 4 8 R z A O S J c x O e p M h M z v r z K w Q Q n 7 C i w d F v P o 7 3 v w b J 8 k e N L G g o a j q p r s r S g U 3 1 v e / v c L a + s b m V n G 7 t L O 7 t 3 9 Q P j x q G p V p h g 2 m h N L t i B o U P M G G 5 V Z g O 9 V I Z S S w F Y 1 u Z 3 7 r C b X h K n m w 4 x R D S Q c J j z m j 1 k n t r p I 4 o L 2 g V 6 7 4 V X 8 O s k q C n F Q g R 7 1 X / u r 2 F c s k J p Y J a k w n 8 F M b T q i 2 n A m c l r q Z w Z S y E R 1 g x 9 G E S j T h Z H 7 v l J w 5 p U 9 i p V 0 l l s z V 3 x M T K o 0 Z y 8 h 1 S m q H Z t m b i f 9 5 n c z G 1 + G E J 2 l m M W G L R X E m i F V k 9 j z p c 4 3 M i r E j l G n u b i V s S D V l 1 k V U c i E E y y + v k u Z F N f C r w f 1 l p X a T x 1 G E E z i F c w j g C m p w B 3 V o A A M B z / A K b 9 6 j 9 + K 9 e x + L 1 o K X z x z D H 3 i f P 7 m H j 7 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " F 9 / v h L s Q P x m T f 6 O 6 t s Y L U D F 6 1 y 0 = " > A A A B 7 3 i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K o M e g F 4 8 R z A O S J c x O e p M h M z v r z K w Q Q n 7 C i w d F v P o 7 3 v w b J 8 k e N L G g o a j q p r s r S g U 3 1 v e / v c L a + s b m V n G 7 t L O 7 t 3 9 Q P j x q G p V p h g 2 m h N L t i B o U P M G G 5 V Z g O 9 V I Z S S w F Y 1 u Z 3 7 r C b X h K n m w 4 x R D S Q c J j z m j 1 k n t r p I 4 o L 2 g V 6 7 4 V X 8 O s k q C n F Q g R 7 1 X / u r 2 F c s k J p Y J a k w n 8 F M b T q i 2 n A m c l r q Z w Z S y E R 1 g x 9 G E S j T h Z H 7 v l J w 5 p U 9 i p V 0 l l s z V 3 x M T K o 0 Z y 8 h 1 S m q H Z t m b i f 9 5 n c z G 1 + G E J 2 l m M W G L R X E m i F V k 9 j z p c 4 3 M i r E j l G n u b i V s S D V l 1 k V U c i E E y y + v k u Z F N f C r w f 1 l p X a T x 1 G E E z i F c w j g C m p w B 3 V o A A M B z / A K b 9 6 j 9 + K 9 e x + L 1 o K X z x z D H 3 i f P 7 m H j 7 8 = < / l a t e x i t > o K K q 6 6 e 6 K U s G N 9 f 1 v b 2 1 9 Y 3 N r u 7 B T 3 N 3 b P z g s H R 0 3 j c o 0 w w Z T Q u l 2 R A 0 K n m D D c i u w n W q k M h L Y i k a 3 M 7 / 1 h N p w l T z Y c Y q h p I O E x 5 x R 6 6 R 2 V 0 k c 0 F 6 1 V y r 7 F X 8 O s k q C n J Q h R 7 1 X + u r 2 F c s k J p Y J a k w n 8 F M b T q i 2 n A m c F r u Z w Z S y E R 1 g x 9 G E S j T h Z H 7 v l J w 7 p U 9 i p V 0 l l s z V 3 x M T K o 0 Z y 8 h 1 S m q H Z t m b i f 9 5 n c z G 1 + G E J 2 l m M W G L R X E m i F V k 9 j z p c 4 3 M i r E j l G n u b i V s S D V l 1 k V U d C E E y y + v k m a 1 E v i V 4 P 6 y X L v J 4 y j A K Z z B B Q R w B T W 4 g z o 0 g I G A Z 3 i F N + / R e / H e v Y 9 F 6 5 q X z 5 z A H 3 i f P 7 s L j 8 A = < / l a t e x i t > ! < l a t e x i t s h a 1 _ b a s e 6 4 = " q U g 8 i X I 0 1 z U D f 6 d o e a e Q G g C g J 1 c = " > A A A B 7 X i c b V D L S g N B E O y N r x h f U Y 9 e F o P g K e y K o M e g F 4 8 R z A O S J c x O e p M x 8 1 h m Z o U Q 8 g 9 e P C j i 1 f / x 5 t 8 4 S f a g i Q U N R V U 3 3 V 1 x y p m x Q f D t F d b W N z a 3 i t u l n d 2 9 / Y P y 4 V H T q E x T b F D F l W 7 H x C B n E h u W W Y 7 t V C M R M c d W P L q d + a 0 n 1 I Y p + W D H K U a C D C R L G C X W S c 2 u E j g g v X I l q A Z z + K s k z E k F c t R 7 5 a 9 u X 9 F M o L S U E 2 M 6 Y Z D a a E K 0 Z Z T j t N T N D K a E j s g A O 4 5 K I t B E k / m 1 U / / M K X 0 / U d q V t P 5 c / T 0 x I c K Y s Y h d p y B 2 a J a 9 m f i f 1 8 l s c h 1 N m E w z i 5 I u F i U Z 9 6 3 y Z 6 / 7 f a a R W j 5 2 h F D N 3 K 0 + H R J N q H U B l V w I 4 f L L q 6 R 5 U Q 2 D a n h / W a n d 5 H E U 4 Q R O 4 R x C u I I a 3 E E d G k D h E Z 7 h F d 4 8 5 b 1 4 7 9 7 H o r X g 5 T P H 8 A f e 5 w + Q M Y 8 b < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " q U g 8 i X I 0 1 z U D f 6 d o e a e Q G g C g J 1 c = " > A A A B 7 X i c b V D L S g N B E O y N r x h f U Y 9 e F o P g K e y K o M e g F 4 8 R z A O S J c x O e p M x 8 1 h m Z o U Q 8 g 9 e P C j i 1 f / x 5 t 8 4 S f a g i Q U N R V U 3 3 V 1 x y p m x Q f D t F d b W N z a 3 i t u l n d 2 9 / Y P y 4 V H T q E x T b F D F l W 7 H x C B n E h u W W Y 7 t V C M R M c d W P L q d + a 0 n 1 I Y p + W D H K U a C D C R L G C X W S c 2 u E j g g v X I l q A Z z + K s k z E k F c t R 7 5 a 9 u X 9 F M o L S U E 2 M 6 Y Z D a a E K 0 Z Z T j t N T N D K a E j s g A O 4 5 K I t B E k / m 1 U / / M K X 0 / U d q V t P 5 c / T 0 x I c K Y s Y h d p y B 2 a J a 9 m f i f 1 8 l s c h 1 N m E w z i 5 I u F i U Z 9 6 3 y Z 6 / 7 f a a R W j 5 2 h F D N 3 K 0 + H R J N q H U B l V w I 4 f L L q 6 R 5 U Q 2 D a n h / W a n d 5 H E U 4 Q R O 4 R x C u I I a 3 E E d G k D h E Z 7 h F d 4 8 5 b 1 4 7 9 7 H o r X g 5 T P H 8 A f e 5 w + Q M Y 8 b < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " q U g 8 i X I 0 1 z U D f 6 d o e a e Q G g C g J 1 c = " > A A A B 7 X i c b V D L S g N B E O y N r x h f U Y 9 e F o P g K e y K o M e g F 4 8 R z A O S J c x O e p M x 8 1 h m Z o U Q 8 g 9 e P C j i 1 f / x 5 t 8 4 S f a g i Q U N R V U 3 3 V 1 x y p m x Q f D t F d b W N z a 3 i t u l n d 2 9 / Y P y 4 V H T q E x T b F D F l W 7 H x C B n E h u W W Y 7 t V C M R M c d W P L q d + a 0 n 1 I Y p + W D H K U a C D C R L G C X W S c 2 u E j g g v X I l q A Z z + K s k z E k F c t R 7 5 a 9 u X 9 F M o L S U E 2 M 6 Y Z D a a E K 0 Z Z T j t N T N D K a E j s g A O 4 5 K I t B E k / m 1 U / / M K X 0 / U d q V t P 5 c / T 0 x I c K Y s Y h d p y B 2 a J a 9 m f i f 1 8 l s c h 1 N m E w z i 5 I u F i U Z 9 6 3 y Z 6 / 7 f a a R W j 5 2 h F D N 3 K 0 + H R J N q H U B l V w I 4 f L L q 6 R 5 U Q 2 D a n h / W a n d 5 H E U 4 Q R O 4 R x C u I I a 3 E E d G k D h E Z 7 h F d 4 8 5 b 1 4 7 9 7 H o r X g 5 T P H 8 A f e 5 w + Q M Y 8 b < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " q U g 8 i X I 0 1 z U D f 6 d o e a e Q G g C g J 1 c = " > A A A B 7 X i c b V D L S g N B E O y N r x h f U Y 9 e F o P g K e y K o M e g F 4 8 R z A O S J c x O e p M x 8 1 h m Z o U Q 8 g 9 e P C j i 1 f / x 5 t 8 4 S f a g i Q U N R V U 3 3 V 1 x y p m x Q f D t F d b W N z a 3 i t u l n d 2 9 / Y P y 4 V H T q E x T b F D F l W 7 H x C B n E h u W W Y 7 t V C M R M c d W P L q d + a 0 n 1 I Y p + W D H K U a C D C R L G C X W S c 2 u E j g g v X I l q A Z z + K s k z E k F c t R 7 5 a 9 u X 9 F M o L S U E 2 M 6 Y Z D a a E K 0 Z Z T j t N T N D K a E j s g A O 4 5 K I t B E k / m 1 U / / M K X 0 / U d q V t P 5 c / T 0 x I c K Y s Y h d p y B 2 a J a 9 m f i f 1 8 l s c h 1 N m E w z i 5 I u F i U Z 9 6 3 y Z 6 / 7 f a a R W j 5 2 h F D N 3 K 0 + H R J N q H U B l V w I 4 f L L q 6 R 5 U Q 2 D a n h / W a n d 5 H E U 4 Q R O 4 R x C u I I a 3 E E d G k D h E Z 7 h F d 4 8  We consider a linear driving protocol for duration τ, such that the time dependence of the frequency is given by,

x M Z i 4 w Z R 4 H t j B g O z a I 3 F f / z O i m G 1 9 1 M q C R F U H z 2 U Z h K i j G d R k L 7 Q g N H O b a E c S 3 s r p Q P m W Y c b X B F G 4 K 3 e P I y a V 5 U P L f i 3 V 2 W q z f z O A r k m J y S c + K R K 1 I l N V I n D c L J I 3 k m r + T N e X J e n H f n Y 9 a 6 4 s x n j s g f O J 8 / T t C Y c A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " J + q p F 2 C m y 3 W I i D M P K 9 A 1 B v l t 0 Q g = " > A A A C B X i c b V A 9 S w N B E N 3 z M 8 a v q K U W i 0 G w C n c i a B m 0 S R n B f E A u h L 3 N X L J k b + / Y n R P C k c b G v 2 J j o Y i t / 8 H O f + P m k k I T H w y 8 f W + G n X l B I o V B 1 / 1 2 V l b X 1 j c 2 C 1 v F 7 Z 3 d v f 3 S w W H T x K n m 0 O C x j H U 7 Y A a k U N B A g R L a i Q Y W B R J a w e h 2 6 r c e Q B s R q 3 s c J 9 C N 2 E C J U H C G V u q V T n w J I f q S q Y E E W q O + F o M h + j p / 9 0 p l t + L m o M v E m 5 M y m a P e K 3 3 5 / Z i n E S j k k h n T 8 d w E u x n T K L i E S d F P D S S M j 9 g A O p Y q F o H p Z v k V E 3 p m l T 4 N Y 2 1 L I c 3 V 3 x M Z i 4 w Z R 4 H t j B g O z a I 3 F f / z O i m G 1 9 1 M q C R F U H z 2 U Z h K i j G d R k L 7 Q g N H O b a E c S 3 s r p Q P m W Y c b X B F G 4 K 3 e P I y a V 5 U P L f i 3 V 2 W q z f z O A r k m J y S c + K R K 1 I l N V I n D c L J I 3 k m r + T N e X J e n H f n Y 9 a 6 4 s x n j s g f O J 8 / T t C Y c A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " J + q p F 2 C m y 3 W I i D M P K 9 A 1 B v l t 0 Q g = " > A A A C B X i c b V A 9 S w N B E N 3 z M 8 a v q K U W i 0 G w C n c i a B m 0 S R n B f E A u h L 3 N X L J k b + / Y n R P C k c b G v 2 J j o Y i t / 8 H O f + P m k k I T H w y 8 f W + G n X l B I o V B 1 / 1 2 V l b X 1 j c 2 C 1 v F 7 Z 3 d v f 3 S w W H T x K n m 0 O C x j H U 7 Y A a k U N B A g R L a i Q Y W B R J a w e h 2 6 r c e Q B s R q 3 s c J 9 C N 2 E C J U H C G V u q V T n w J I f q S q Y E E W q O + F o M h + j p / 9 0 p l t + L m o M v E m 5 M y m a P e K 3 3 5 / Z i n E S j k k h n T 8 d w E u x n T K L i E S d F P D S S M j 9 g A O p Y q F o H p Z v k V E 3 p m l T 4 N Y 2 1 L I c 3 V 3 x M Z i 4 w Z R 4 H t j B g O z a I 3 F f / z O i m G 1 9 1 M q C R F U H z 2 U Z h K i j G d R k L 7 Q g N H O b a E c S 3 s r p Q P m W Y c b X B F G 4 K 3 e P I y a V 5 U P L f i 3 V 2 W q z f z O A r k m J y S c + K R K 1 I l N V I n D c L J I 3 k m r + T N e X J e n H f n Y 9 a 6 4 s x n j s g f O J 8 / T t C Y c A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " J + q p F 2 C m y 3 W I i D M P K 9 A 1 B v l t 0 Q g = " > A A A C B X i c b V A 9 S w N B E N 3 z M 8 a v q K U W i 0 G w C n c i a B m 0 S R n B f E A u h L 3 N X L J k b + / Y n R P C k c b G v 2 J j o Y i t / 8 H O f + P m k k I T H w y 8 f W + G n X l B I o V B 1 / 1 2 V l b X 1 j c 2 C 1 v F 7 Z 3 d v f 3 S w W H T x K n m 0 O C x j H U 7 Y A a k U N B A g R L a i Q Y W B R J a w e h 2 6 r c e Q B s R q 3 s c J 9 C N 2 E C J U H C G V u q V T n w J I f q S q Y E E W q O + F o M h + j p / 9 0 p l t + L m o M v E m 5 M y m a P e K 3 3 5 / Z i n E S j k k h n T 8 d w E u x n T K L i E S d F P D S S M j 9 g A O p Y q F o H p Z v k V E 3 p m l T 4 N Y 2 1 L I c 3 V 3 x M Z i 4 w Z R 4 H t j B g O z a I 3 F f / z O i m G 1 9 1 M q C R F U H z 2 U Z h K i j G d R k L 7 Q g N H O b a E c S 3 s r p Q P m W Y c b X B F G 4 K 3 e P I y a V 5 U P L f i 3 V 2 W q z f z O A r k m J y S c + K R K 1 I l N V I n D c L J I 3 k m r + T N e X J e n H f n Y 9 a 6 4 s x n j s g f O J 8 / T t C Y c A = = < / l a t e x i t >
for the compression stroke and, for the expansion stroke. Note that for this driving protocol, Equation (24) can be solved analytically and yields solutions in terms of the Airy functions [76].
To determine the average heat and work exchanged during each stroke, we use the following method. We determine the internal energy at the beginning of the compression stroke (point A in the cycle) from a thermal equilibrium state with the cold bath. Then using Equation (33), we determine the internal energy at the end of the compression stroke (point B). The change in internal energy between A and B gives the average work done on the system (W 1 ). We can find the the internal energy at point C using the thermal equilibrium state with the hot bath. The change in internal energy between C and B gives the average heat exchanged with the hot bath (Q 2 ). Again, applying Equation (33), we find the internal energy at the end of the expansion stroke (point D). The change in internal energy between C and D gives the average work done by the system (W 2 ). Finally, we can then use the difference in internal energy between points D and A to find the average heat exchanged with the cold bath (Q 4 ).
With the average work and heat for each stroke in hand, we can then calculate the efficiency from the ratio of the total work to the heat input and the power from the ratio of total work to the cycle time, We consider isentropic strokes of equal duration, τ 1 = τ 3 = τ. As we do not explicitly model the system-bath interaction, we represent the duration of the isochoric strokes as a multiplicative factor of the isentropic stroke duration. Thus, we can represent the total cycle time as τ cyc = 2γτ. The full expressions for the efficiency and power are cumbersome and detailed in Appendix B.
In Figure 2, we plot the engine efficiency as a function of the ratio of bath temperatures. Confirming the results of [57], we see that the efficiency is greatest in the non-interacting bosonic limit of ν = 0 and least in the non-interacting fermionic limit of ν = 1. Between these limits, we observe that increasing the interaction strength from zero to one interpolates the efficiency smoothly between the bosonic and fermionic bounds. In Figure 3, we plot the engine efficiency as a function of the stroke time, τ. Again, we see that ν = 0 provides the greatest efficiency and ν = 1 the worst, with intermediate values of ν falling between these limits. We see that increasing the stroke time results in increasing efficiency, approaching the bound of η = 1 − ω 0 /ω τ achieved in the limit of perfectly adiabatic strokes (Q * = 1). We note that in the limit of long stroke times the efficiency converges to this limit for all values of ν. This indicates that the influence of the interaction on the engine performance is a fundamentally nonequilibrium effect. The oscillatory behavior of the efficiency arises from the form of Q * for the linear protocol.
Next, we examine the power output of the engine as a function of the ratio of bath temperatures, as shown in Figure 4. As with the efficiency, we see that the bosonic limit (ν = 0) gives the greatest power output, and the fermionic limit (ν = 1) gives the least. Intermediate values of ν fall between these bounds. In Figure 5, we plot the power as a function of the stroke time. As in the case of the efficiency, we see that the shift in the power output arising from the interaction vanishes as the stroke time increases.   Examining both efficiency and power plots, we observe that the zeroes for efficiency and power all occur at finite values of the bath temperature ratio or stroke time. These zeroes mark the transition to the parameter regimes where the cycle no longer satisfies the conditions, where we use the convention that work or heat flowing into the system is positive. These are known as the positive work conditions and must be met for the cycle to function as an engine. In the parameter regimes where the positive work conditions are not satisfied, efficiency and power are no longer meaningful metrics of performance. In these regimes, the cycle functions instead as a heater, accelerator, or refrigerator. In a heater, work is put into the system to induce heat flow into both baths, ( Q 2 < 0, Q 4 < 0, W tot > 0). In an accelerator, work is put into the system to enhance the flow of heat from the hot to the cold bath ( Q 2 > 0, Q 4 < 0, W tot > 0). Lastly, in a refrigerator, work is put into the system to induce heat flow from the cold to the hot bath ( Q 2 < 0, Q 4 > 0, W tot > 0). We note that the values of the bath temperature ratio and stroke time for which the system transitions out of the engine regime vary for different values of ν. To further explore this relationship, in Figure 6, we show the parameter space under which the cycle functions as each type of thermal machine for different interaction strengths. We see that for short stroke times, ν = 0 displays the smallest heater and accelerator regimes, and the largest engine regime, while ν = 0 displays the largest heater and accelerator regimes, and the smallest engine regime. Intermediate values of ν fall between these limits. In general, we see that as cycle time increases, the heater and accelerator regimes vanish, as do the differences in the parameter space for the different values of ν.

Discussion
In this work, we have analyzed the performance of a quantum Otto engine with a working medium of two interacting particles in a singular oscillator potential, including the efficiency, power output, and parameter regimes where the cycle functions as different types of thermal machines. We have determined the full time-dependent dynamics using the framework of statistical anyons and shown that by changing the strength of the interparticle interaction, we can interpolate between the performance of bosonic and fermionic working mediums. We see that the impact of the singular interaction on the engine performance vanishes for long stroke times, as the engine approaches fully adiabatic performance, indicating that the performance differences are a fundamentally nonequilibrium phenomena. Consistent with the results of [57], we have found that the best performance in regards to both efficiency and power arises in the limit of ν = 0, in which the interparticle interaction vanishes and the particles behave as ideal bosons in a pure harmonic potential. While interacting particles yield reduced performance in comparison to ideal bosons, they yield enhanced performance in comparison to ideal fermions. As performance depends directly on the interaction, such an engine may have other applications as well, for instance using differences in performance to measure the generalized exclusion statistics parameter of an anyonic system. By incorporating interparticle interactions, we extend our exploration of the impact of wave function symmetry on heat engine performance to a more physically realistic model. The inverse square interaction potential in particular arises in the case of electron-dipole interactions [81,82] and, as previously noted, as the basis of the Calogero-Sutherland model. In general, solving the dynamics of interacting systems is significantly more challenging than in the idealized case. Here, we demonstrate the usefulness of the statistical anyon framework in overcoming this challenge, as it provides a simple and novel method for mapping the dynamics of interacting particles to that of a mixture of ideal bosons and fermions.
Introducing a time-dependent interaction strength would open up the possibility of optimizing performance by varying the behavior of the particles between that of bosons and fermions during the engine cycle. Engine power output could also be optimized by minimizing the duration of the compression and expansion strokes while suppressing nonadiabatic excitations [83], which can be accomplished using shortcuts to adiabaticity [84,85]. However, for an accurate performance assessment in such cases, it is important that the cost of implementing the shortcut is also accounted for [29][30][31]35,[86][87][88]. We leave these as topics to be explored in future work.
The Calogero-Sutherland model described by Equation (1) has received extensive study in the context of generalized exclusion statistics, and due to its direct connection to spin-chain models, such as the Haldane-Shastry chain [89][90][91][92], which provide a promising route for experimental realization in trapped atom systems [93][94][95][96]. Such systems may provide a means of experimentally implementing a singular quantum heat engine. Furthermore, generalized exclusion statistics anyons can be used to replicate the thermodynamic behavior of fractional exchange statistics anyons that manifest in two-dimensional systems [24]. These anyons are highly sought after due to their close relationship to non-ableian anyons, which may be used to implement fault-tolerant topological quantum computers [97]. However, direct observation and manipulation of fractional exchange statistics anyons is extremely difficult [98]. Using the well-established framework of heat engines to understand the thermodynamic behavior of systems that display intermediate statistics, such as the singular oscillator examined here, may lead to better methods of detection and control for fractional exchange statistics anyons.

Conflicts of Interest:
The authors declare no conflict of interest.