Surprises from quenches in long-range interacting systems: Temperature inversion and cooling

# Surprises from quenches in long-range interacting systems: Temperature inversion and cooling

## Abstract

What happens when one of the parameters governing the dynamics of a long-range interacting system of particles in thermal equilibrium is abruptly changed (quenched) to a different value? While a short-range system, under the same conditions, will relax in time to a new thermal equilibrium with a uniform temperature across the system, a long-range system shows a fast relaxation to a nonequilibrium quasistationary state (QSS). The lifetime of such an off-equilibrium state diverges with the system size, and the temperature is non-uniform across the system. Quite surprisingly, the density profile in the QSS obtained after the quench is anticorrelated with the temperature profile in space, thus exhibiting the phenomenon of temperature inversion: denser regions are colder than sparser ones. We illustrate with extensive molecular dynamics simulations the ubiquity of this scenario in a prototypical long-range interacting system subject to a variety of quenching protocols, and in a model that mimics an experimental setup of atoms interacting with light in an optical cavity. We further demonstrate how a procedure of iterative quenching combined with filtering out the high-energy particles in the system may be employed to cool the system. Temperature inversion is observed in nature in some astrophysical settings; our results imply that such a phenomenon should be observable, and could even be exploitable to advantage, also in controlled laboratory experiments.

###### pacs:
05.70.Ln, 52.65.Ff, 37.30.+i

Max-Planck-Institut für Physik Komplexer Systeme, Nöthnitzer Straße 38, D-01187 Dresden, Germany
Dipartimento di Fisica e Astronomia and CSDC, Università di Firenze, and INFN, sezione di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, I-50125 Firenze, Italy Keywords: Long-range interactions, Vlasov equation, Quasi-stationary states, Temperature inversion

## 1 Introduction

Long-range-interacting systems abound in nature [1, 2, 3]. Interactions are called long-range when the pair potential energy decays asymptotically with the interparticle distance as , with in spatial dimensions, as it happens for instance in Coulomb and gravitational interactions. One striking feature of many-particle long-range interacting systems is that they are typically found in states that are out of thermodynamic equilibrium. This is at variance with systems with short-range interactions: while nonequilibrium states in short-range systems require an enduring external forcing to counteract the tendency of the system to reach due to “collisional” effects a thermal equilibrium state, an isolated long-range interacting system evolving while starting far from equilibrium will remain stuck in an out-of-equilibrium state for very long times. Such times grow with the number of degrees of freedom, so that they may become larger than any experimentally accessible timescale. For instance, the time needed for a galaxy to reach an “equilibrium” state3 is estimated to be several orders of magnitude larger than the age of the universe [4, 5]. This property is universal, in the sense that it is shared by all systems with long-range interactions, regardless of the details of the interaction potential 4. It is understood as a consequence of the fact that the kinetic equation governing the evolution of the single-particle distribution function , where and are the canonically conjugated positions and momenta, respectively, is the Vlasov equation (often referred to as the collisionless Boltzmann equation, mainly in the astrophysical literature), up to a timescale when collisional effects can no longer be neglected, with diverging with [1, 3]. Such an equation has infinitely many stationary solutions: among these, the stable ones are identified with the observed non-equilibrium states, often referred to as quasi-stationary states (QSSs), and the thermal equilibrium state is just one out of infinitely many.

The overall qualitative picture of the dynamical evolution of an isolated long-range system starting far from thermal equilibrium is the following: After a transient (referred to as violent relaxation after Lynden-Bell [6]) whose lifetime does not depend on , the system sets into a QSS and stays there for a time of the order of . For times , when the Vlasov description no longer applies, the QSS is no longer stationary, and the system eventually evolves to a thermal equilibrium state. The QSS the system relaxes to after the violent relaxation depends on the initial conditions, and no general theory is available yet to predict it. A completely general statistical approach was proposed by Lynden-Bell [6], but it rests on the hypothesis of complete mixing of the Vlasov dynamics that only rarely holds, hence, gives reasonably accurate predictions only in very particular cases [3, 7]. Other methods have been proposed since then, which give good predictions for simple models only for special classes of initial conditions, namely, for single-level initial distributions (the so-called “water-bag” distributions) [3, 7] or when the relaxation towards the QSS is nearly adiabatic [7, 8, 9, 10].

Although predicting the full distribution function in a QSS resulting from a generic initial condition is a formidable task, especially for spatially inhomogeneous QSSs whose stability properties are particularly difficult to analyze, one may ask whether the QSSs that are reached in certain situations do share some common features. One of the hallmarks of thermal equilibrium is a uniform temperature throughout the system, hence, also in inhomogeneous equilibrium states, the temperature will be constant. Conversely, in a nonequilibrium state, the temperature may well be nonuniform. Hence, one may ask whether some properties of the temperature profile in a QSS share a certain degree of universality. A physically relevant case is that of a system prepared in a (spatially inhomogeneous) thermal equilibrium state and then brought out of equilibrium by means of a (non-small) sudden perturbation, acting for a very short time. As shown in recent papers [11, 12], in this case, the typical outcome is somewhat surprising and counterintuitive, and is referred to as temperature inversion. The latter effect implies that the density and temperature profiles are anticorrelated: sparser regions of the system are hotter than denser ones. Temperature inversions are observed in nature in astrophysical settings, the most famous example being the solar corona, where the temperature grows from thousands to millions of Kelvin while going from the photosphere to the sparser external regions of the corona [13]. Other examples of temperature inversions have been observed, for instance, in molecular clouds [14]. As argued in [11, 12, 15], temperature inversion should not be a peculiarity of systems in extreme conditions, as astrophysical systems typically are, but should rather be a generic property of long-range interacting systems relaxing to a QSS after having been brought out of (spatially inhomogeneous) thermal equilibrium by a perturbation. In [11], temperature inversion was observed as a result of perturbations of a thermal state of a prototypical model with long-range interactions, the Hamiltonian mean-field (HMF) model [16], and of a two-dimensional self-gravitating system. Temperature inversion in the latter system is thoroughly investigated in [17] in connection with filamentary structures in galactic molecular clouds. Partial temperature inversions were also observed in QSSs of two-dimensional self-gravitating systems whose initial conditions were of the water-bag type [18]. In one-dimensional self-gravitating systems, temperature inversions in QSSs were observed, which gradually disappear during the slow relaxation of the system from the QSS to thermal equilibrium [19]. Temperature inversion was recently demonstrated to occur in the nonequilibrium stationary state of a class of mean-field systems involving rotators subject to quenched disordered external drive and dissipation [20]. A physical picture of the origin of temperature inversion in a generic long-range interacting system, based on the interplay between a mechanism originally proposed to explain the temperature profile in the solar corona (referred to as velocity filtration [21, 22, 23, 24]) and the interaction of the system particles with the time-dependent mean field, was also suggested [11].

If the phenomenon of temperature inversion is somewhat universal in long-range interacting systems abruptly brought out of thermal equilibrium, one may wonder whether it could be observed also in a controlled laboratory experiment. In the present paper, we aim at taking a step forward towards a positive answer to this question. First of all, using molecular dynamics (MD) simulations, we demonstrate that in the HMF model, with both attractive and repulsive interactions, temperature inversion is typically observed when the system is brought out of thermal equilibrium by means of an abrupt change (from now on referred to as a quench) of a parameter controlling the dynamics of the system, like a coupling constant or an external field. To this end, we employ a variety of different quench protocols. Our studies are relevant to experiments because quench protocols can be performed with experimental setups, especially in systems of (cold) atoms; moreover, an HMF model with repulsive interactions and an external confining field can be seen as a toy model of trapped ions, while an HMF model with attractive interactions does share many features with atoms interacting with light in a cavity. We then elaborate on this point by demonstrating that at least one of the considered quench protocols can be applied to another model, which mimics more closely an experimental setup of atoms interacting with a standing electromagnetic wave in an optical cavity [25, 26, 27], again yielding temperature inversion. It is worth noting that the quenching protocol applicable to the latter model does correspond to changing a parameter that can be actually tuned in an experiment, that is, the intensity of an external laser pump. We also give evidence that the physical picture proposed in [11] does apply to temperature inversions produced by quenches. Finally, we show that an iterative quenching protocol, combined with filtering out the high-energy particles localized in the sparser external regions, can be used to cool the system.

The paper is organized as follows: in Sec. 2, we introduce the HMF model, and study the temperature inversion produced by quenching the external field in the case of repulsive interactions (§2.1.1) as well as that produced in the case with attractive interactions by quenching the external field (§2.2.1) and the coupling constant (§2.2.2). In Sec. 3, we discuss a model of atoms interacting with light in an optical cavity, with emphasis on the dissipationless limit (§3.1). We then report on temperature inversion as observed in the latter system after quenching the coupling constant (§3.1.1), and discuss its possible observation in a laboratory experiment (§3.1.2). Section 4 is devoted to demonstrating that an iterative application of the previously considered quenching protocols can cool the systems. Finally, conclusions are drawn and open issues are discussed in Sec. 5.

## 2 Temperature inversion in the Hamiltonian mean-field (HMF) model

The prototypical model of long-range interacting systems that we shall consider in this paper is the so-called Hamiltonian mean field (HMF) model [16]. The model comprises a system of globally-coupled point particles of unit mass moving on a circle. The Hamiltonian of the system, in presence of an external field of strength , is given by

 H=N∑i=1p2i2+J2NN∑i,j=1[1−cos(ϑi−ϑj)]−hN∑i=1cosϑi. (1)

Here, is the angular coordinate of the -th particle on the circle, while is the conjugated momentum, see Fig. 1(a). The coupling constant can be either positive or negative, defining respectively the ferromagnetic (F) and the antiferromagnetic (AF) version of the model. The reason for the use of terms borrowed from the physics of magnetic systems is that the HMF model can also be seen as a system of planar () spins with mean-field couplings (that is, defined on a complete graph with links of the same strength). For , therefore, the Hamiltonian (1) is -invariant. Ferromagnetic (respectively, antiferromagnetic) interactions in the magnetic picture correspond to attractive (respectively, repulsive) interactions in the particle interpretation. The time evolution of the system is governed by the Hamilton equations derived from the Hamiltonian (1):

 ⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩dϑidt=pi,dpidt=J(mycosϑi−mxsinϑi)−hsinϑi, (2)

where are the components of the magnetization vector , whose magnitude measures the amount of clustering of particles on the circle. A uniform (respectively, non-uniform) distribution of particles on the circle implies the value (respectively, ).

Let us denote by the time-dependent single-particle distribution function, which is the probability density to have a particle with angular coordinate and momentum at time . In thermal equilibrium5 at temperature , the distribution has the form [3]

 feq(ϑ,p)∝exp{−1T[p22−J(meqxcosϑ+meqysinϑ)−hcosϑ]}, (3)

where the equilibrium magnetization components are obtained self-consistently as . In Eq. (3), we have set the Boltzmann constant to unity, as we shall do throughout the paper. In the special case , that is, for distribution of particles symmetric around that we shall focus on in the following, the self-consistent relation for the magnetization may be written as [1, 3]

 meq=I1[β(Jmeq+h)]I0[β(Jmeq+h)] , (4)

where , and denotes the modified Bessel function of order . For , the system in thermal equilibrium is characterized by an inhomogeneous distribution of particles on the circle, with clustering around . In the AF case, a nonzero external field is necessary to have an inhomogeneous thermal state, while in the F case, a clustered equilibrium state is attained even with , via a spontaneous breaking of the symmetry by the attractive interaction between the particles, provided the total energy density is smaller than a critical value , see Fig. 1(b); the corresponding critical temperature is [1, 3]. The F-HMF system with in fact shows a continuous phase transition as a function of temperature, with decreasing continuously from unity at to zero at , while remaining zero at higher temperatures [1, 3].

Let us note that the Hamiltonian (1) can be obtained by taking any long-range interacting system, restricting to one spatial dimension, expanding the potential energy in a Fourier series, and retaining just the first Fourier mode [7, 11]. The AF-HMF model can be seen to represent a one-component Coulomb system, while the F-HMF is a simplified description of a self-gravitating system. As anticipated in the Introduction, and as we shall see further in Sec. 3, the F-HMF model turns out to be very closely related to a different system, i.e., a system of atoms in an optical cavity.

### 2.1 The antiferromagnetic case

Let us now discuss how quenching one of the parameters of the Hamiltonian (1) when the system is in an inhomogeneous thermal equilibrium state produces a non-equilibrium state exhibiting temperature inversion. To start with, we consider the AF case with .

#### Quenching the external field

We consider particles and prepare the system in an inhomogeneous equilibrium state () at an initial value of the field and temperature , by sampling independently for every particle the coordinate and the momentum from the distribution (3). We take the coupling to be . In such a state, while the local density , defined as

 n(ϑ)≡∫dpf(ϑ,p,t), (5)

is non-uniform, corresponding to a magnetization , the local temperature, defined as

 T(ϑ)≡∫dpp2f(ϑ,p,t)n(ϑ), (6)

is nevertheless uniform as a function of . We study the time evolution of the system starting from this equilibrium state by performing a MD simulation that involves integrating6 the equations of motion (2). At , we instantaneously quench the field to . Soon after the quench, the magnetization of the system starts oscillating; the oscillations damp out in time, and eventually the system relaxes to a QSS with at the new value of the field. The latter fact is expected since the value of is smaller than in the initial state. The time evolution of the magnetization for a time window around the time of the quench is shown in Fig. 2.

For times larger than those shown in Fig. 2, the magnetization stays essentially constant with extremely small fluctuations, signaling that the system has relaxed to a QSS. The density and temperature profiles as a function of in the QSS are plotted in Fig. 3. It is apparent that and are anticorrelated. The system exhibits a clear temperature inversion: denser regions are colder than sparser ones. We note that the average temperature

 ⟨T⟩≡∫dϑT(ϑ)n(ϑ) (7)

has the value in the QSS, lower than the initial value . We shall come back to this point in §4.

Apart from the reduction of the average temperature, these results are qualitatively very similar to those obtained in [11] for a ferromagnetic HMF model. In that work, however, the simulation protocol was different, because the system was initiated in equilibrium with and then an external field was applied for a very short time. In Ref. [11], a physical picture was put forward to explain the emergence of temperature inversion, and we argue that it applies also in the present case (and, as we shall discuss below, to all the cases discussed in the present paper). Without entering into details, let us recall the main points of this explanation. After the quench (or the perturbation), a collective oscillation (i.e., a wave) sets in, as witnessed by the time evolution of the magnetization (Fig. 2). The oscillation damps out in time, but the system is conservative with no dissipative mechanism being present. Hence, the energy lost by the wave must be acquired by some of the particles. This may happen via Landau damping7, an ubiquitous phenomenon in systems with long-range interactions: particles whose velocity is not too far from the phase velocity of the wave interact strongly with the wave itself. The net result is that the nearly resonant particles gain kinetic energy from the wave, so that the momentum distribution develops small peaks or shoulders and becomes different from a thermal distribution. In particular, close to the resonances, there will be suprathermal regions, i.e., regions where is larger than the initial thermal distribution. Now, the mechanism of velocity filtration comes into play as follows. The density is maximum at and decreases for larger and smaller ’s, reaching its minimum around . The net effect of the interaction between the particles and the external field is to create an effective force field pushing the particles towards , so that the particles have to climb the potential energy well to reach the sparser regions of the system; particles with larger kinetic energies will do that more easily than “colder” particles, and any suprathermal region present in the velocity distribution will be magnified in regions where the particle density is lower. This is precisely what happens in our case, as shown in Fig. 4, where the velocity distribution function measured at a position where the particle density is smaller, , see panel (a), is plotted together with the same function measured at the position of maximum density, , see panel (b). Magnified suprathermal regions in are apparent, and are responsible for the fact that the variance of the velocity distribution is here larger than in the denser regions of the system.

We performed other numerical experiments with the same quench protocol by varying the initial temperature and the difference between the initial and final values of the field , and obtained qualitatively similar results. No fine-tuning of the parameters is necessary to observe temperature inversion. We repeated some of the numerical experiments with particles, observing no differences with respect to the case.

In principle, one can quench also the coupling constant , while keeping the field fixed. However, as stated before, as far as laboratory systems are concerned, the AF-HMF system can be seen as a toy model of a system of equal electric charges in a confining field, i.e., ions in a trap. In trapped ions, quenching corresponds to quenching the strength of the confining trap, which in principle is an experimentally feasible protocol. Conversely, quenching is equivalent to quenching the charge of the ions, which does not seem very physical, so that we do not study this kind of quench protocol for the AF case.

### 2.2 The ferromagnetic case

Let us now turn to the ferromagnetic (F-HMF) case, that is, . Without loss of generality, we take the coupling to be .

#### Quenching the external field

We repeat for the F-HMF system the numerical experiment we performed with the AF system. We prepare the system of particles in inhomogeneous thermal equilibrium at temperature and with a field , now corresponding to , and evolve the system up to , when we quench the external field to . As in the AF case, the magnetization starts oscillating (Fig. 5), the oscillation damps out in time, and the system settles in a QSS with temperature inversion (Fig. 6).

Figures 5 and 6 are very similar to the corresponding ones for the AF system, that is, Figs. 2 and 3, respectively. The same is true for the momentum distributions, shown in Fig. 7, which is strikingly similar to Fig. 4.

Hence, also in this case, the physical picture based on wave-particle interaction and velocity filtration applies. Consistent with our physical picture, the sign of the interactions is irrelevant for field quenches, as long as the field is sufficiently strong to produce the dominant effects of shaping the density profile and providing the effective potential well that filters the velocities of the particles, thereby producing temperature inversion. Again, we repeated our numerical experiments for different values of the fields, for different initial temperatures and for systems with , observing no qualitative differences. We note that similar to the AF case, the average temperature in the QSS for the F case is smaller than the initial temperature, and one has .

#### Quenching the coupling constant

As we have previously discussed, the F-HMF model admits magnetized (collapsed) thermal equilibrium states also for , provided . Moreover, as we shall see in Sec. 3, the F-HMF model also admits an interpretation in terms of an atomic system, where the coupling constant is tunable in an experiment. We thus turn to analyze what happens if we prepare an F-HMF system in thermal equilibrium with nonzero magnetization and , and then quench the coupling constant. We prepare the system of particles in inhomogeneous thermal equilibrium at temperature and coupling , by sampling independently for every particle the coordinate and the momentum from the distribution (3). As before, we evolve the system until , when we instantaneously quench the coupling to . As shown in Figs. 8, 9 and 10, the behavior of the system is similar to the case when the external field was quenched; although the effect is quantitatively less dramatic, still there is a clear temperature inversion in the QSS, where the average temperature has the value , lower than the initial value of . Note that in contrast to the case when we quenched the strength of the field, here the damping of the oscillations in with time is not quite complete (see Fig. 8); indeed, oscillations are sustained also at long times, though the amplitude of oscillation is small ( at , compared to the amplitude subsequent to the quench).

The physical picture based on wave-particle interactions and velocity filtration holds also in this case. We performed other numerical experiments with larger ’s and different values of and, as before, we can conclude that no fine tuning is needed to produce temperature inversion by means of a quench protocol, acting on either or on .

Let us now discuss how the F-HMF model is related to a system of atoms interacting with a standing electromagnetic wave in an optical cavity, and how a quench of the coupling constant could be performed in a controlled way in a laboratory experiment.

## 3 Temperature inversion in a system of atoms in an optical cavity

Atoms interacting with a single-mode standing electromagnetic wave due to light trapped in a high-finesse optical cavity are subject to an interparticle interaction that is long-ranged owing to multiple coherent scattering of photons by the atoms into the wave mode [25, 26, 27]. The system serves as a unique platform to study long-range interactions, and in particular mean-field interactions, under tunable experimental conditions. The typical setup is sketched in Fig. 11, which also shows optical pumping by a transverse laser of intensity to counter the inevitable cavity losses quantified by the cavity linewidth (the lifetime of a photon in the cavity being ). We follow Ref. [25, 26, 27], and refer the reader to these works and to references therein for all the technical details that we skip here. Let us then consider a system of identical atoms of mass . As far as the interactions with the electromagnetic field is concerned, each atom may be regarded as a two-level system, where the transition frequency between the two levels is . If the atoms are confined in one dimension along the cavity axis (taken to be the -axis), and is the wavenumber of the standing wave, the sum of the electric field amplitudes coherently scattered by the atoms at time depends on their instantaneous positions , and is proportional to the quantity

 Θ≡1NN∑j=1cos(kxj), (8)

so that the cavity electric field at time is [25]. Here, is the maximum intracavity-photon number per atom, given by , with being the ratio between the cavity vacuum Rabi frequency and the detuning between the laser and the atomic transition frequency, and being the detuning between the laser and the cavity-mode frequency. It is important to note that is tunable by controlling either the strength of the external laser pump or the detuning . The quantity characterizes the amount of spatial ordering of atoms within the cavity mode, with corresponding to atoms being uniformly distributed and the resulting vanishing of the cavity field, and implying spatial ordering. In particular, corresponds to perfect ordering of atoms to form Bragg gratings and the resulting scattering of photons in phase, thereby maximizing the cavity field, which in turn traps the atoms by means of mechanical forces. Note that the wave number is related to the linear dimension of the cavity through and , where is the wavelength of the standing wave, and .

The dynamics of the system may be conveniently studied by analyzing the time evolution of the -atom phase space distribution at time , with ’s denoting the momenta conjugate to the positions . Treating the cavity field quantum mechanically, and regarding the atoms as classically polarizable particles with semiclassical center-of-mass dynamics, it was shown that the distribution evolves in time according to the Fokker-Planck equation (FPE) [25, 26]

 ∂tfN+{fN,H}=−¯¯¯nΓLfN. (9)

Here, the operator describes the dissipative processes (damping and diffusion) that explicitly depend on the positions and momenta of the atoms8, and , where is the recoil frequency due to collision between an atom and a photon, and is the reduced Planck constant. The Hamiltonian is given by [25, 26]

 H=N∑j=1p2j2m+NℏΔc¯¯¯nΘ2. (10)

This semiclassical limit is valid under the condition of being larger than , and the Hamiltonian describes the conservative dynamical evolution of in the limit of vanishing cavity losses, or for times sufficiently small such that dissipative effects are negligible. The Hamiltonian (10) contains the photon-mediated long-ranged (mean-field) interaction between the atoms encoded in the quantity . Note that the interaction is attractive (respectively, repulsive) when is negative (respectively, positive). On time scales sufficiently long so that the effect of the right-hand-side of Eq. (9) is non-negligible, the mean-field description of the dynamics is no longer valid [25, 27].

### 3.1 The dissipationless limit and the connection with the HMF model

Let us now consider the case of effective attractive interactions between the atoms and the cavity field (), and study the dynamics of the system in the limit in which the effect of the dissipation can be neglected, that is, for sufficiently small times. In this limit, the dynamics of the atoms is conservative and governed by the Hamiltonian given by Eq. (10). The positions of the atoms enter the Hamiltonian only as , so that we may define the phase variables

 ϑj=kxj=2πxjλ , (11)

for . Now, the length of the cavity is times the wavelength , with an integer. Then, setting the origin of the -axis in the center of the cavity, we have , so that on using the periodicity of the cosine function, we can take the phase variables modulo such that . Then, by measuring lengths in units of the reciprocal wavenumber of the cavity standing wave, masses in units of the mass of the atoms , and energies in units of , the Hamiltonian can be rewritten in dimensionless form as

 H=N∑j=1(pϑ)2j2−¯¯¯nNΘ2, (12)

where, in terms of the variables, is now expressed as

 Θ=1NN∑j=1cosϑj. (13)

The ’s in Eq. (12) are the momenta canonically conjugated to the variables.

The similarity between the system with Hamiltonian (12) and an F-HMF model, already noted in [25], is now well apparent, since the expression of the interaction field given in Eq. (13) coincides with the expression of the -component of the magnetization of the HMF model. Indeed, the equations of motion derived from the Hamiltonian (12) read

 ⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩dϑjdt=(pϑ)j,d(pϑ)jdt=−¯¯¯nΘsinϑj, (14)

for . Equations (14) coincide with Eqs. (2) once , , and . Hence, the dynamics of a system of atoms interacting with light in a cavity in the dissipationless limit is equivalent to that of a model that differs from the ferromagnetic HMF model in zero field just for the fact that particles in the former interact only with the -component of the magnetization. Equivalently, the dynamics is equivalent to that of an HMF model with . This means that the thermal equilibrium distribution for atoms in a cavity (still assuming that dissipative effects can be neglected) is still given by Eq. (3), with and .

#### Quenching the coupling constant

In an HMF model in the thermodynamic limit , the condition holds for all times if it holds at , that is, if the initial spatial distribution is symmetric around . In a finite system, even starting from a symmetric distribution, will not stay exactly zero for all times, due to fluctuations induced by finite-size effects, but will remain very small. Hence, the numerical experiment reported in Sec. 2.2.2, which was performed while starting from an equilibrium state with , almost directly applies also to the model of atoms in an optical cavity (with ), but for the fact that there was a coupling of the particles to the small fluctuating -component of the magnetization in the former that would be absent in the atom case. It is reasonable to expect that this coupling has only a small effect, since is very large (). To check this, we repeated the numerical experiment using the dynamics given by Eqs. (14), that is, we prepared the system in the same initial condition at and , then evolved the dynamics—now using Eqs. (14)—until when we quenched the coupling to the new value . The results are shown in Figs. 12, 13 and 14, that are essentially indistinguishable from Figs. 8, 9 and 10 of Sec. 2.2.2. There is clear temperature inversion, the momentum distribution exhibits the same features as before, and the average temperature has the value in the QSS, lower than the initial value of .

#### From numerical to laboratory experiments?

What do the results we have shown suggest as regards the possibility of observing temperature inversion in a laboratory experiment with atoms interacting with light in an optical cavity? To answer this question, let us consider the system and the values of the parameters that were already considered in Ref. [25] so that the semiclassical approximation for the atomic dynamics holds. We thus consider a system of atoms, whose mass is kg; the atomic transition is the line, with a wavelength nm. In terms of MHz, the half-width of the line, the cavity linewidth is and the detuning is . The energy unit is thus J. Preparing a system of atoms in a thermal state with a temperature of order unity in these energy units means reaching temperatures of the order of K, which can be achieved with laser cooling techniques. Hence, the initial state we have considered in our simulations can be prepared in a laboratory, and, as already noted, can be tuned in the range we considered (and even in a much wider range, if needed) so that also the quench of is feasible. The crucial point is to understand whether the dissipationless limit corresponding to assuming a Hamiltonian dynamics reasonably describes the dynamics of the atoms over the timescale needed for the relaxation of the system to the QSS with temperature inversion. With our choice of energy, length and mass units, the time unit is fixed to

 τ=√mλ204π2ℏΔc≈10−6 s, (15)

this has to be compared with the timescale that rules the dissipative effects, s. Luckily enough, the simulations reported in Refs. [25, 26] show that the dissipative effects set in on a timescale that is of the order of , while our simulations show that the QSS with temperature inversion is reached after times of the order of , so that there should be ample room for measurements of temperature inversion before it is destroyed by dissipation.

## 4 Iterative cooling via temperature inversion

We have already noted in all the cases considered above that the system (be it an HMF or an atomic system) is colder in the QSS with temperature inversion than it was in the initial thermal equilibrium state. Moreover, the fact that the QSS has temperature inversion means that hotter particles reside in the external parts of the system, where the density is smaller. This suggests that temperature inversion could be exploited to make the system even colder by means of an iterative protocol that goes as follows. We prepare the system in an inhomogeneous thermal state at temperature symmetric about , perform a quench, and let the system relax to the QSS with temperature inversion as before. Now, we filter out hotter particles, that is, all the particles that are located at , with a suitable choice of : the resulting system will be close to isothermal at a temperature . We let the system relax for some time and then perform another quench. The system goes into another QSS with temperature inversion and average temperature . The protocol can be iterated as long as it continues to appreciably cool the system, or as long as there is room for another quench of the parameters.

We now show the results obtained by applying this iterative cooling protocol to the cases studied before.

### 4.1 Cooling the antiferromagnetic HMF model

We start with the AF-HMF model. After performing the quench described in Sec. 2.1.1, whereby starting in equilibrium at , and , the field was quenched to , the system is in a QSS with temperature and density profiles shown in Fig. 3. Now, we filter out the hotter (high-energy) particles, which we choose to be those with ’s lying outside the interval . Subsequently, we instantaneously quench the field from to ; the system eventually settles into a QSS, and the corresponding local density and local temperature profiles are shown in Fig. 15 (a). By further filtering out the high-energy particles and quenching instantaneously the field to , the system goes to another QSS whose density and temperature profiles are reported in Fig. 15 (b). The average temperature as a function of time for the full sequence of events starting with equilibrium at and is shown in Fig. 16. The system cools down from an initial average temperature to a final value . Each small downward step in corresponds to filtering out the high-energy particles, while each large downward step corresponds to settling of the system into a QSS with a lower temperature as a result of quenching of the field to a lower value. The percentage decrease in the number of particles during filtering out the high-energy particles equals about during the first stage of filtration, and about during the second stage. Not surprisingly, the cooling protocol becomes less efficient at subsequent stages.

### 4.2 Cooling the ferromagnetic HMF model and atoms in a cavity

We now apply the same cooling protocol described above to the F-HMF model in an external field, while starting from equilibrium at , and . The results are shown in Figs. 17 and 18. Again, we obtain a cooling from an initial average temperature to a final value . The percentage decrease in the number of particles during filtering out the high-energy particles equals about during the first stage of filtration, and equals about during the second stage.

We now turn to the case of a quench of the coupling constant, still for the F-HMF model but now without external field. The cooling protocol goes as before, but we can fruitfully perform only one stage of filtering instead of two. We start with the equilibrium state at and that is quenched to , yielding the QSS whose density and temperature profiles are shown in Fig. 9. We then filter out the hotter particles, in this case defined as those with ’s lying outside the interval . Subsequently, we instantaneously quench the coupling to . The density and temperature profile in the resulting QSS are shown in Fig. 19, and the time evolution of the average temperature is reported in Fig. 20. Under the full quenching protocol, the system cools down from an initial average temperature to a final value . The percentage decrease in the number of particles during filtering out the high-energy particles equals about .

Finally, we apply exactly the same cooling protocol to the dynamics of atoms in a cavity, given by Eqs. 14. The results are shown in Figs. 21 and 22, and as expected are very close to that obtained for the F-HMF model: the system under the quenching protocol cools down, from an initial average temperature to a final value . The percentage decrease in the number of particles during filtering out the high-energy particles equals about .

Although the cooling protocol does not yield dramatic results, especially in the last two considered cases, still it seems able to (at least) decrease the temperature of a system by a factor of two. Moreover, it opens up the possibility of a different way of cooling a system with respect to the nowadays commonly used ones.

## 5 Concluding remarks

In this paper, we have demonstrated that quasi-stationary states that an isolated long-range interacting system relaxes to when starting from a spatially inhomogeneous thermal equilibrium and subsequently subject to a quench of one of the parameters of the Hamiltonian generically exhibit the phenomenon of temperature inversion. Using molecular dynamics simulations of a prototypical model, the HMF model, we have shown that different quench protocols lead to temperature inversion, in presence and in absence of an external field, and with both attractive and repulsive interactions: no fine-tuning is necessary in any of the cases. This lends further support to the claim made in [11, 15] that, rather than a peculiarity of some astrophysical systems, temperature inversion should be a “universal” feature of nonequilibrium quasi-stationary states resulting from perturbations of spatially inhomogeneous thermal equilibrium. Moreover, our analysis of the momentum distribution functions supports the interpretation put forward in [11] that temperature inversion arises due to an interplay between wave-particle interaction (the role of the wave being played by the oscillation of the mean field induced by the perturbation or the quench) and velocity filtration. The latter is a mechanism proposed by Scudder [21, 22, 23] to explain the temperature profile of the solar corona, and basically amounts to the statement that a suprathermal velocity distribution becomes broader when climbing a potential well.

The fact that quench protocols are able to produce temperature inversions is of particular relevance in view of observing this phenomenon in controlled laboratory experiments. Exploiting the close connection between the HMF model and a model describing the dynamics of a system of atoms interacting with light in an optical cavity, as noted in Ref. [25], we have argued that a quench of the external laser pump should be able to produce a quasi-stationary state with temperature inversion in a system of, say, Rubidium atoms in a cavity that is initially in thermal equilibrium at temperatures reachable by laser cooling. Such a quasi-stationary state should survive for quite a long time before being destroyed by dissipative effects. The latter claim is made on the assumption that the timescales of dissipation are the same as those observed in [25, 26]. Similar to our study, these works also considered a quench of the coupling constant, but from different initial states and to different final states, and this may affect the dissipation timescales. However, results in [26] seem to indicate that dissipation is inhibited in spatially inhomogeneous states, and this would be an advantage for a quench experiment like those suggested here. An interesting follow-up of the present work, in view of a detailed study of the feasibility of an experiment, would surely be to use the techniques of [25, 26] to study the evolution after quenches like the ones described here under a dynamics that fully takes into account the dissipative effects. Another interesting system to study in view of a possible experiment to detect temperature inversion would be a system of trapped ions, which may be related to the antiferromagnetic HMF model. Work is in progress along this direction [33].

Finally, we have shown that temperature inversion may be exploited to cool the system. The cooling protocol is of iterative type, and takes advantage of the fact that in states with temperature inversion hotter particles are spatially separated from colder ones, so that they can be filtered out. Let us sign off by saying that temperature inversion is one of many fascinating features of the nonequilibrium states of long-range interacting systems that are yet to be unveiled.

## 6 Acknowledgments

We thank P. Di Cintio, H. Landa and A. Trombettoni for fruitful discussions and suggestions. SG acknowledges the support and hospitality of INFN (Italy) and the University of Florence.

### Footnotes

1. Corresponding author; Email: shamikg1@gmail.com
2. Email: lapo.casetti@unifi.it
3. For non-confined self-gravitating systems with a finite mass, the “true” thermodynamic equilibrium state where the velocities obey a Maxwell-Boltzmann distribution is never reached; this estimate refers to the time needed to reach a state where the memory of the initial condition has been completely lost [4, 5].
4. This is true provided there are only long-range interactions; this property does not hold for systems with mixed long- and short-range interactions.
5. Despite long-range interactions implying a non-additive Hamiltonian, in the HMF model in the thermodynamic limit, canonical and microcanonical ensembles are equivalent for any temperature and energy density, so that one can speak of thermal equilibrium without specifying if it is a canonical or a microcanonical equilibrium [1, 3].
6. In all the MD simulations reported in this paper, we used a fourth-order symplectic algorithm [28], with time step . This ensured energy conservation up to a relative fluctuation of .
7. The theory of Landau damping is completely understood for homogeneous states and small perturbations [29, 30]; its extension to inhomogeneous states still has many open issues [31, 32]. However, the basic physical mechanism works in any situation where particles interact with a collective excitation in the system.
8. We omit the explicit expression of because it is not relevant in the context of the present paper; It may be found in Ref. [25, 26].

### References

1. Campa A, Dauxois T and Ruffo S 2009 Statistical mechanics and dynamics of solvable models with long-range interactions Phys. Rep. 480 57
2. Bouchet F, Gupta S and Mukamel D 2010 Thermodynamics and dynamics of systems with long-range interactions Physica A: Special Issue FPSP XII 389 4389
3. Campa A, Dauxois T, Fanelli D and Ruffo S 2014 Physics of Long-Range Interacting Systems (Oxford: Oxford University Press)
4. Choudouri A R 2010 Astrophysics for Physicists (Cambridge: Cambridge University Press)
5. Binney J and Tremaine S 2008 Galactic Dynamics, 2nd ed. (Princeton: Princeton University Press)
6. Lynden-Bell D 1967 Statistical mechanics of violent relaxation in stellar systems Mon. Not. R. Astron. Soc. 136 101
7. Levin Y, Pakter R, Rizzato F B, Teles T N and Benetti F P C 2014 Nonequilibrium statistical mechanics of systems with long-range interactions Phys. Rep. 535 1
8. Benetti F P C, Ribeiro-Teixeira A C, Pakter R and Levin Y 2014 Nonequilibrium Stationary States of 3D Self-Gravitating Systems Phys. Rev. Lett. 113 100602
9. Leoncini X, Van Der Berg T N and Fanelli D 2009 Out-of-equilibrium solutions in the XY-Hamiltonian Mean-Field model Europhys. Lett. 86 20002
10. De Buyl P, Mukamel D and Ruffo S 2011 Self-consistent inhomogeneous steady states in Hamiltonian mean-field dynamics Phys. Rev. E 84 061151
11. Teles T N, Gupta S, Di Cintio P and Casetti L 2015 Temperature inversion in long-range interacting systems Phys. Rev. E 92 020101(R)
12. Teles T N, Gupta S, Di Cintio P and Casetti L 2016 Reply to ‘Comment on “Temperature inversion in long-range interacting systems”’ Phys. Rev. E 93 066102
13. Golub L and Pasachoff J M 2009 The Solar Corona (Cambridge: Cambridge University Press)
14. Myers P C and Fuller G A 1992 Density structure and star formation in dense cores with thermal and nonthermal motions Astrophys. J. 396 631
15. Casetti L and Gupta S 2014 Velocity filtration and temperature inversion in a system with long-range interactions Eur. Phys. J. B 87 91
16. Antoni M and Ruffo S 1995 Clustering and relaxation in Hamiltonian long-range dynamics Phys. Rev. E 52 2361
17. Di Cintio P, Gupta S and Casetti L 2016 Dissipationless collapse and the origin of temperature inversion in galactic filaments in preparation
18. Teles T N, Levin Y, Pakter R and Rizzato F P 2010 Statistical mechanics of unbound two-dimensional self-gravitating systems J. Stat. Mech.: Theory Exp. P05007
19. Worrakitpoonpon T 2011 Relaxation of isolated self-gravitating systems in one and three dimensions PhD Thesis
20. Campa A, Gupta S and Ruffo S 2015 Nonequilibrium inhomogeneous steady state distribution in disordered, mean-field rotator systems J. Stat. Mech.: Theory Exp. P05011
21. Scudder J 1992 On the causes of temperature change in inhomogeneous low-density astrophysical plasmas Astrophys. J. 398 299
22. Scudder J 1992 Why all stars should possess circumstellar temperature inversions Astrophys. J. 398 319
23. Scudder J 1994 Ion and electron suprathermal tail strengths in the transition region: Support for the velocity filtration model of the corona Astrophys. J. 427 446
24. A primer on velocity filtration is given in the Supplementary material to [11] at http://link.aps.org/supplemental/10.1103/PhysRevE.92.020101
25. Schütz S and Morigi G 2014 Prethermalization of atoms due to photon-mediated long-range interactions Phys. Rev. Lett. 113 203002
26. Schütz S, Jäger S B and Morigi G 2015 Relaxation after quenches in systems with photon-mediated long-range interactions arXiv:1512.05243
27. Jäger S B, Schütz S and Morigi G 2016 Mean-field theory of atomic self-organization in optical cavities arXiv:1603.05148
28. McLachlan R I and Atela P 1992 The accuracy of symplectic integrators Nonlinearity 5 541
29. Choudhouri A R 1998 The Physics of Fluids and Plasmas (Cambridge: Cambridge University Press)
30. Mohout C and Villani C 2011 On Landau damping Acta Mathematica 207 29
31. Barré J, Olivetti A and Yamaguchi Y Y 2010 Dynamics of perturbations around inhomogeneous backgrounds in the HMF model J. Stat. Mech.: Theory Exp. P08002
32. Barré J, Olivetti A and Yamaguchi Y Y 2011 Algebraic damping in the one-dimensional Vlasov equation J. Phys. A: Math. Theor. 44 405502
33. Landa H, Gupta S, Di Cintio P and Casetti L 2016 in preparation
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters