Quantum Torque on a Non-Reciprocal Body out of Thermal Equilibrium and Induced by a Magnetic Field of Arbitrary Strength

A stationary body that is out of thermal equilibrium with its environment, and for which the electric susceptibility is non-reciprocal, experiences a quantum torque. This arises from the spatially non-symmetric electrical response of the body to its interaction with the non-equilibrium thermal fluctuations of the electromagnetic field: the non-equilibrium nature of the thermal field fluctuations results in a net energy flow through the body, and the spatially non-symmetric nature of the electrical response of the body to its interaction with these field fluctuations causes that energy flow to be transformed into a rotational motion. We establish an exact, closed-form, analytical expression for this torque in the case that the environment is the vacuum and the material of the body is described by a damped oscillator model, where the non-reciprocal nature of the electric susceptibility is induced by an external magnetic field, as for magneto-optical media. We also generalise this expression to the context in which the body is slowly rotating. By exploring the high-temperature expansion of the torque, we are able to identify the separate contributions from the continuous spectral distribution of the non-reciprocal electric susceptibility, and from the resonance modes. In particular, we find that the torque persists in the limiting case of zero damping parameter, due to the contribution of the resonance modes. We also consider the low-temperature expansion of the torque. This work extends our previous consideration of this model to an external magnetic field of arbitrary strength, thereby including non-linear magnetic field effects.

In Ref. [13], we considered the quantum torque on a body made of nonreciprocal material that is out of thermal equilibrium with its environment, an expression for which was first obtained in Ref. [10], and subsequently in Ref. [12]; see also Ref. [11].Broadly, such a torque arises as follows: the nonequilibrium nature of the thermal fluctuations of the electromagnetic field results in a net energy flow through the body, and the spatially non-symmetric electrical response of the body to its interaction with these field fluctuations causes that energy flow to be transformed into a rotational motion.Generally, for a single body, the existence of a quantum torque appears to require both that the system is out of thermal equilibrium and that it exhibits a broken spatial symmetry [4].
In our previous analysis, the environment was the vacuum and the material of the body was described by a damped oscillator model, where the nonreciprocal nature of the electric susceptibility was induced by an external magnetic field, as in the case of magneto-optical media [9][10][11].There, we considered only contributions to first order in the magnetic field strength, which provides a good numerical approximation where the corresponding cyclotron frequency is much less than the damping parameter for the oscillator model.Here, we extend our previous analysis to consideration of an external magnetic field of arbitrary strength, thereby including non-linear magnetic field effects.
The present paper therefore has quite a narrow focus.It is primarily concerned with investigation of the structure of the quantum torque for such a model for arbitrary external magnetic field strength, and, in particular, its spectral decomposition.Accordingly, it has a slightly more mathematical flavour than our previous analysis in Ref. [13].
The remainder of this paper is structured as follows.In Section 2, we briefly review the damped oscillator model for the electromagnetic response of the constituents of the material of the body, and the corresponding spatially nonsymmetric electric susceptibility that is induced by an external magnetic field.In Section 3, we consider the non-equilibrium quantum torque on a stationary such body, and establish an exact, closed-form, analytical expression for this torque that is valid for arbitrary external magnetic field strength.The corresponding torque for a slowly rotating such body is also considered.In Section 4, we explore the high-temperature expansion of the quantum torque for a stationary non-reciprocal body, which enables the identification of the separate contributions from the continuous spectral distribution of the non-reciprocal electric susceptibility of the body, and from the resonance modes.In particular, we find that, in the limit of vanishing damping parameter for the oscillator model, the contribution of the resonance modes results in a persistent quantum torque.The corresponding low-temperature expansion is considered in Section 5 and brief conclusions are presented in Section 6.
We use Heaviside-Lorentz electromagnetic units, and set ℏ = c = 1, except where numerical values are presented.

Non-Reciprocal Electric Susceptibility
As in Ref. [13], we employ the following damped oscillator model to describe the non-equilibrium displacement, u(t), at time t, of the particles that constitute the material of a body, in response to the application of a fluctuating electric field, E(t), and a constant external magnetic field, B: where e is the particle charge, m is the particle mass, η is the damping parameter, and ω 0 is the free oscillation frequency.For a metal, ω 0 = 0, and the model then reduces to that of Refs.[9,10]; see also Refs.[6,11,14,15].This dynamical equation is easily solved in the frequency domain to yield the polarisation, P(ω) = neu(ω) ≡ χ(ω)E(ω), which is here expressed in terms of the electric susceptibility, where, without loss of generality, the external magnetic field has been chosen to lie in the z-direction, ω c ≡ e m B is the corresponding cyclotron frequency, and ω 2 p ≡ ne 2 m is the square of the plasma frequency, n being the particle density.For B ̸ = 0, this electric susceptibility is clearly non-symmetric and, moreover, non-reciprocal, that is, Re χ(ω) is non-symmetric.
For a metal, we may set ω 0 = 0 and use the charge and mass of the electron.For gold, the parameter values at room temperature are ω p = 9 eV and η = 0.035 eV [16].For an external magnetic field strength of 1 T, the corresponding cyclotron frequency is ω c ∼ 10 −4 eV, which is much smaller than the damping parameter, so, in practice, keeping only terms linear in ω c provides a good numerical approximation, one that we followed in our previous work [13].Here, however, our interest is primarily in the structure of the torque, and, in particular, its spectral decomposition, rather than its magnitude.We therefore include non-linear magnetic field effects to all orders, and the expressions we derive are exact in ω c .
In the dilute limit, the mean polarisability of the body is given by where the electric susceptibility may, in general, be position-dependent, and the integration extends over the domain of the body.Thus, for a non-reciprocal such body, 1 Re α(ω) is, in general, non-symmetric.
In particular, for a homogeneous such body of volume V that is made of non-reciprocal material with electric susceptibility as in Eq. (2.2) and with ω 0 = 0, which we will henceforth assume, where Alternatively, Eq. (2.4) may be expressed in terms of partial fractions, as that is, (2.9) 3 Quantum Torque on a Non-Reciprocal Body In Ref. [13], we used the Fluctuation-Dissipation Theorem to express the quantum torque on a stationary non-reciprocal body, with mean polarisability α(ω) in the dilute limit, as where β ′ is the inverse temperature of the body and β is the inverse temperature of the environment, which is taken to be the vacuum.See also Refs.[10,12].Thus, from Eq. (2.5), the only non-zero component of the quantum torque for the damped oscillator model with ω 0 = 0 is which becomes that is, where we have employed the integral representation of the digamma function [17][18][19][20], , Re s > 0.
(3.5)Quantum Torque on a Non-Reciprocal Body out of Thermal Equilibrium Keeping only terms to first order in ω c , it is easily verified that Eq. (3.4) agrees with Eq. (4.15) of our earlier work [13].
For a sense of scale, note that for a gold nanosphere of radius 100 nm, with the parameter values given in Section 2, in an external magnetic field of strength 1 T, the magnitude of the coefficient of the log term in Eq. (3.4) is 1.6 × 10 −24 Nm.
The extension to a slowly rotating body was also considered in Ref. [13]; see also Ref. [10].There, the z-component of the quantum torque on a nonreciprocal body undergoing a non-relativistic rotation about the z-axis with angular frequency Ω, and with mean polarisability α(ω) in the dilute limit, was expressed as where ω + ≡ ω + Ω. Keeping only terms to first order in Ω, we may write where τ 0 z is given by Eq. (3.1) and Thus, from Eq. (2.9) for the damped oscillator model with ω 0 = 0, whence, (3.12) Keeping only terms to first order in ω c , it is easily verified that Eq. (3.12) agrees with Eq. (5.10) of our earlier work [13].
The quantum torque on a rotating non-reciprocal body, made of material described by the damped oscillator model with ω 0 = 0, is therefore given, to first order in Ω, by Eq. (3.7), together with Eqs.(3.4) and (3.12).

High-Temperature Expansion
It is instructive to consider the high-temperature expansion of the quantum torque.We do so here only for a stationary body; the same approach could also be applied to a slowly rotating body.
For high temperatures, β, β ′ → 0 + , use of the power series representation [18][19][20], in Eq. (3.4) yields where It is easily seen that ξ k satisfies the reflection formula and the recurrence relation The latter may be used to recursively generate ξ k .For example, the first few ξ k for non-negative k, generated in this way, are as follows: ) ) ) ) ) Alternatively, these expressions may be established directly by using a suitable multiple-angle formula, as in the following, for k ≥ 1: where U n denotes the Chebyshev polynomial of the second kind of order n.
It follows immediately from Eq. (4.5), or from Eqs. (4.13) and (4.4), that, in the limit of vanishing damping parameter, the behaviour of ξ k differs, depending on whether the index, k ∈ Z, is even or odd: As a consequence, only the odd β, β ′ powers in Eq. (4.2) survive in this limit, resulting in lim It may be slightly surprising to note that lim η→0 + τ 0 z ̸ = 0, since it is clear from Eq. (2.4) that if η = 0 then Re α xy (ω) = 0, and it then follows from Eq. (3.1) that τ 0 z = 0.That is, τ 0 z is a discontinuous function of η as η → 0 + .
In fact, this behaviour is due to the existence of resonance modes at ω = ±ω c − iη in the complex-ω plane, corresponding to poles of the integrand in Eq. (3.2).In the limit η → 0 + , these poles approach the real line from below, and the integral in Eq. (3.2) then becomes formally divergent.However, it may be assigned a finite value by invoking the Sokhotski-Plemelj Theorem, which may be written in the form lim where P denotes the corresponding Cauchy Principal Value (CPV) integral.See, for example, Ref. [21], where it is shown that the CPV integral may, equivalently, be expressed as an Analytic Principal Value integral.Indeed, using the partial fraction representation for Re α xy (ω) given in Eq. (2.7), Eq. (3.1) may be alternatively written as where Echoing the decomposition in Eq. (4.16), Eq. (4.17) has the limit lim where the CPV contribution, indicated by P, and the resonance contribution, indicated by R, are respectively given by lim since the integrand in Eq. (4.17) is real in this limit, and It is easily verified, using the exponential generating function of the Bernoulli numbers, that Eq. (4.21) is a more compact equivalent of Eq. (4.15).Thus, the fact that lim η→0 + τ 0 z ̸ = 0 is due entirely to the contribution of the resonance modes.Quantum Torque on a Non-Reciprocal Body out of Thermal Equilibrium The persistence of the quantum torque in the limit of vanishing damping found here is reminiscent of similar behaviour found in our earlier work [22] for a charged particle passing a conducting plate with permittivity governed by the Drude model.There, the transverse magnetic component of the classical electrodynamic friction on the particle was found to persist in the limit of vanishing damping (resistivity) parameter, again due to the presence, in this limit, of a singularity in the integrand for the friction.
Let us now generalise the above to finite η.Using the residue theorem to replace the integral over the real line in Eq. (4.17) by that over the contour C displayed in Figure 1, we may write where the CPV integral over the contour C excludes the semi-circles around the poles at ω = ±ω c − iη, and is evaluated in the limit ε → 0 + , where it is equal to the corresponding integral over the real line, excluding the points ω = ±ω c .

Im ω
Re ω 1 The contour C used for the evaluation of τ 0 z .The corresponding CPV integral over C excludes the semi-circles around the poles at ω = ±ωc − iη, and is evaluated in the limit ε → 0 + .Thus, we may write where and Again using Eq.(4.22), it is easily verified that Eq. ( 4.26) may be written as which agrees with the odd β, β ′ powers in Eq. (4.2).Also, it is immediate that Eq. (4.26) agrees with Eq. (4.21) in the limit η → 0 + .To evaluate Eq. (4.25), it suffices to again express the integrand in the form seen in Eq. (3.2), but accommodate the CPV by using the high-frequency expansion of Im 1 ω 2 +ξ 2 away from the poles at ω = ±ω c − iη.This results in where we have used the functional equation Quantum Torque on a Non-Reciprocal Body out of Thermal Equilibrium to evaluate the formally divergent integrals by analytic continuation and have interpreted the k = 1 term as the limit, See Ref. [23] for interpretation of the finite part of the divergent integrals that may result from term-by-term integration, and how missing terms in such integration arise from singularities of the integrand.As expected, Eq. (4.28) agrees with the even β, β ′ powers, and logarithms, in Eq. (4.2), and, from Eq. (4.14), with Eq. (4.20) in the limit η → 0 + .So, we may conclude that the continuous spectral distribution of Re α xy (ω) in Eq. (3.1) generates the even β, β ′ powers, and logarithms, in the hightemperature expansion of the quantum torque, Eq. (4.2), while the discrete spectral distribution of Re α xy (ω), arising from the resonance modes, generates the odd β, β ′ powers.
However, it should be noted that, like the power series representation, Eq. (4.1), Eq. (4.2) is an exact expression for the quantum torque, simply one that has been written in a form that is suitable for high temperature expansion.Thus, the evenness and oddness of the two types of contribution above under change of sign of the inverse temperatures is a property that enables closed-form expressions for these two types of contribution to be easily identified.
Evenness or oddness under change of sign of the inverse temperatures is suggestive of similar behaviour under time reversal, since the corresponding thermal Green functions are periodic in imaginary time, with period equal to the relevant inverse temperature.Indeed, the model in Eq. (2.1) is invariant under time reversal if we also reverse the signs of both η and ω c , the latter corresponding to change in the sign of B. It is easily seen from Eq. (4.13) and Eq.(4.4) that, for k ∈ Z, Thus, alternatively, the two types of contribution to the quantum torque may be extracted by writing and (4.39)

Low-Temperature Expansion
The low-temperature expansion of the quantum torque is less interesting, but is still worth considering.Again, for illustration, we do so here only for a stationary body.For low temperatures, β, β ′ → ∞, use of the asymptotic series representation [18,20], (5.2) On the other hand, use of the low-frequency expansion of Im in agreement with the asymptotic expression, Eq. (5.2).So, for low temperatures, τ 0 z may again be decomposed as in Eq. ( 4.24), with the CPV contribution as in Eq. ( 5.3) and the resonance contribution as in Eq. (4.26), but, in this case, the resonance contribution is exponentially damped (that is, the resonance modes are not excited at low temperatures), and so it does not feature in the asymptotic representation, Eq. (5.2).

Conclusion
In this paper, we have considered the quantum torque experienced by a nonreciprocal body that is out of thermal equilibrium with its environment, where the material of the body is described by a damped oscillator model and where the non-reciprocal nature of the electric susceptibility is induced by an external magnetic field of arbitrary strength, thereby extending our previous consideration of this model [13] to include non-linear magnetic field effects.
We have established an exact, closed-form, analytical expression for the quantum torque on a stationary such body, and have generalised this expression to the context in which the body is slowly rotating.
By exploring the high-temperature expansion of the quantum torque on a stationary non-reciprocal body, we have been able to identify the separate contributions to this torque arising from the continuous spectral distribution of the electric susceptibility, and from the resonance modes.These two types of contribution have been found to differ in parity: that from the continuous spectral distribution is even under combined reflection of the inverse temperatures of the body and of its environment, while that from the resonance modes is odd under such reflection.
In particular, we have demonstrated that, while the contribution to the quantum torque due to the continuous spectral distribution vanishes in the limit of vanishing damping parameter, that due to the resonance modes does not, resulting in a persistent quantum torque in this limit.The quantum torque is, therefore, a discontinuous function of the damping parameter in this limit.
We have also considered the low-temperature expansion of the quantum torque on a stationary non-reciprocal body, where the two types of contribution again display the parity characteristics described above, but where that due to the resonance modes is exponentially damped.
There are, of course, limitations to our analysis.In particular, we have considered the mean polarisability of the non-reciprocal body only in the dilute limit, and have assumed that the temperatures of the body and of its environment do not change with time, so that the non-equilibrium nature of the configuration is maintained.These limitations are discussed in our previous work [13], and are not explored further here.
In addition, for simplicity, we have restricted attention only to the case in which the environment is the vacuum.It would be interesting to extend this work to encompass more general environments, such as a dielectric medium, to explore related phenomena, such as angular momentum transfer, and to consider more realistic physical models for nonreciprocal bodies, including magneto-optical media.

. 1 )
Quantum Torque on a Non-Reciprocal Body out of Thermal Equilibrium in Eq. (3.4) yields