Cooling strongly self-organized particles using adiabatic demagnetization

Abstract

This study examines the dynamics of polarizable particles, coupled to a lossy cavity mode, that are transversally driven by a laser. The analysis is performed in a regime where the cavity linewidth exceeds the recoil frequency by several orders of magnitude. Using a two-stage cooling protocol, we show that the particles’ kinetic energy can be reduced to the recoil energy. This cooling protocol relies in its first stage on a high laser power such that the particles cool into a strongly self-organized pattern. This can be seen as a strongly magnetized state. In a second stage, we adiabatically ramp down the laser intensity such that the particles’ kinetic energy is transferred to their potential energy and the particles are “demagnetized”. In this second stage, we optimize the ramping speed, which needs to be fast enough to avoid unwanted heating and slow enough such that the dynamics remains approximately adiabatic.

1 Introduction

The realization of quantum technologies (Ladd et al., 2010; Acín et al., 2018; Barzanjeh et al., 2022) based on polarizable particles such as atoms, ions, molecules, and nanoparticles relies on the precise control of their motional degrees of freedom. One important step for achieving full control of these particles is to reduce their residual motion. A key technique to achieve this is laser cooling (Wineland and Itano, 1979; Chu, 1998; Wieman et al., 1999; Cohen-Tannoudji, 1998; Phillips, 1998; Stenholm, 1986; Metcalf and Van der Straten, 1999), which can be used to achieve temperatures that leave particles close to their zero-point motion. The basic principle behind laser cooling is the enhanced absorption rate of laser photons that lower the particle’s momentum. Subsequently, incoherent scattering of a photon from the particle into free space results in lower kinetic energy of the particle. Despite the great success of laser cooling, one major problem is that it typically relies on closed transitions and the atomic species at hand. This hinders the universal application of conventional laser cooling techniques to more complex systems such as molecules or nanoparticles.

A good candidate for overcoming this problem is cavity cooling, where a particles’ motion is cooled by coherent scattering of laser photons (Horak et al., 1997; Vuletić and Chu, 2000; Domokos et al., 2001; Domokos and Ritsch, 2002; Black et al., 2003; Maunz et al., 2004; Morigi et al., 2007; Schleier-Smith et al., 2011; Wolke et al., 2012; Hosseini et al., 2017). Here, the particles’ kinetic energy is carried away by the scattered cavity photons while the internal state of the particles’ remain almost unaltered. The simplest form of cavity cooling requires driving with a laser frequency that is red-detuned with respect to the cavity resonance. That way the cavity promotes the emission of blue-shifted photons which leaves the atoms in average at lower energy. In such a setup the minimum temperature is typically bounded by the linewidth of the cavity (Domokos et al., 2001). Cavity cooling of single atoms (Maunz et al., 2004) and collective cooling (Black et al., 2003; Hosseini et al., 2017) have been realized in experimental labs. Since cavity cooling does not rely on incoherent scattering from a specific internal state, it is has been proposed for cooling molecules (Morigi et al., 2007) and experimentally realized for cooling nanoparticles (Asenbaum et al., 2013; Delić et al., 2019). Although sub-recoil cooling has been achieved experimentally (Wolke et al., 2012) in an optical cavity with very narrow linewidth, the limit set by the cavity linewidth usually lies well above the recoil limit.

This paper investigates a situation where the cavity linewidth is orders of magnitude larger than the recoil frequency, which is, for instance, the case for the experiment described in Hosseini et al. (2017) but also in several other experiments. We demonstrate that it is theoretically possible to achieve temperatures that are of the order of a single recoil by using a combination of cavity cooling and adiabatic control of optomechanical forces. The key ingredient is that the scattered photons, besides cooling, also mediate collective interactions which allow the particles to self-organize (Domokos and Ritsch, 2002; Asbóth et al., 2005). Self-organization occurs if the driving-laser power exceeds a threshold determined by the cavity parameters and the temperature of the particles. Here, the particles spontaneously form a pattern with a spacing that is determined by the wavelength of the light and allows for the constructive interference of scattered photons. Atomic self-organization has been observed with ultra-cold bosons (Baumann et al., 2010), thermal atoms (Arnold et al., 2012), and ultra-cold fermions (Wu et al., 2023; Helson et al., 2023). The formation of a self-organized pattern can be described as a ferromagnetic phase of a long-range interacting system where the collectively scattered light field can be understood as a parameter that measures the magnetization of the atomic ensemble (Schütz et al., 2015).

This study aims to present a protocol which can lower the kinetic energy of polarizable particles close to the recoil limit, even if the cavity linewidth is orders of magnitude wider. It thus proposes a two-stage cooling protocol that uses both cavity cooling and self-organization to bring the particles to a final kinetic energy that is of the order of the recoil energy. The first stage uses collective cavity cooling of many particles with high laser power. The final temperature of the particles is here mostly determined by the cavity linewidth while the particles form a strongly self-organized (magnetized) pattern (Figure 1a). For these parameters, while the atoms possess a rather high kinetic energy, they are tightly confined in space in a pattern which supports constructive interference of scattered laser photons. In the second stage, laser power is slowly decreased such that the magnetization of the particles is adiabatically decreased (Figure 1c). Like the magnetocaloric effect and the principle of adiabatic demagnetization (Tishin and Spichkin, 2003), this results in a decrease of the magnetization of the particles and simultaneously lowers their kinetic energy (Figure 1d). In contrast, however, we do not ramp an external magnetic field but the laser driving amplitude which effectively reduces the particle–particle interactions. This principle is also related to so-called release–retrap or adiabatic trap relaxation protocols which are common in optical lattices and are used to achieve low temperatures and high phase-space densities (DePue et al., 1999; Hu et al., 2017). In such protocols, the particles are cooled in tightly confined trapping potentials while a subsequent adiabatic trap relaxation lowers kinetic energy even further. In contrast, however, the demagnetization presented here is of a collective nature and comes from strong cavity-mediated atom–atom interactions instead of deep laser trapping. This is important, since it allows us to perform the demagnetization fast enough such that cavity shot noise does not significantly heat the system while decreasing the driving-laser power. At the end of this ramp, particles reach a final temperature that can be orders of magnitude lower than that of conventional cavity cooling while the particles reach a spatially homogeneous state (Figure 1b).

FIGURE 1

This paper is structured as follows. Section 2 introduces the semiclassical equations that are used to simulate the system. Furthermore, we show analytical predictions for final kinetic energy following an adiabatic ramp. In Section 3, we analyze the effects of dissipation and show the actual proposed cooling protocol. Conclusions are drawn in Section 4, and Supplementary Appendix A provides details of the calculations in Section 2.

2 Physical setup

We consider a setup of transversally driven polarizable particles with mass inside a single-mode cavity. The particles are driven far off-resonant with detuning between laser frequency and transition frequency such that spontaneous emission and the population of the excited state can be neglected. The laser light is thus coherently scattered with rate by the particles into the cavity. Here, is the Rabi frequency of the driving-laser field and is the vacuum Rabi frequency of the cavity. We assume that the coupling between the particles and the cavity is proportional to the mode function where denotes the wave number of the cavity mode. The laser frequency is red-detuned to the frequency of the single resonator mode with detuning . Furthermore, the cavity mode loses photons at rate . In what follows, we discard effects of the dynamical stark shift . This is possible if and are much larger than .

2.1 Semiclassical description

We now present a semiclassical description of the particles’ center of mass motion and the cavity field. The coupled equations for the motion of the particles with position and momentum , and the real and imaginary part of the cavity field and evolve according to the following stochastic differential equations (Domokos et al., 2001)

and . The noise terms have vanishing first moments, , while the second moments fulfill , , and . Furthermore, the order parameter or magnetization is defined by

Equations 1a, b, c, and d describe the driven-dissipative dynamics of the particles that couple to a dissipative cavity mode.

To better understand the forces that are mediated by the cavity, it is useful to eliminate the cavity degrees of freedom from the dynamics. Here, we work in the limit where and , with the recoil frequency. This implies that the cavity degrees evolve much faster and can be adiabatically eliminated (Schütz et al., 2013). Here, is the Doppler-width and is the single-particle momentum width. Working in this regime allows us to simplify Equations 1a, b, c, and d by calculating the adiabatic stationary state of the cavity field. This is done by formally integrating the differential equations for and . The adiabatic solution is given by

Using this result in Equations 1a, b, c, and d results in

with

We emphasize that is positive since we assumed , which allows for self-organization and cavity cooling (Asbóth et al., 2005). The dynamics given by Equation 5a, b can be rewritten using an effective Hamiltonian

with . The term is a long-range interaction potential which tries to maximize the value of . The latter is the order parameter or magnetization and is used to distinguish between the self-organized and the spatially homogeneous phase. In this context, the values of can be seen as a continuous magnetization for each atom which takes values between and . In the spatially homogeneous or paramagnetic phase, takes random values between and such that . In the self-organized or ferromagnetic phase, the particles form a pattern with a periodicity that is determined by the wavelength such that , meaning that the individual spins fulfil either all or all . In an experiment, the magnetization can be detected from the cavity output. This can be seen by finding the stationary state of Equations 1c, d that can be used to calculate the intra-cavity photon number

where we used Equations 3, 4 and the average runs over different initializations and trajectories.

After adiabatic elimination of the cavity degrees of freedom, we derived a dynamical description from a classical Hamiltonian. This implies that Equations 5a, b conserve the mean energy for a time-independent interaction strength . The description by means of Hamiltonian dynamics is, however, only true on a timescale where dissipative effects can be discarded (Schütz et al., 2015; Schütz et al., 2013; Schütz and Morigi, 2014; Jäger et al., 2016; Schütz et al., 2016).

In the following, we are interested in changing very slowly such that the particles evolve mainly adiabatically but do so sufficiently fast such that dissipative effects are negligible.

2.2 Adiabatic ramp of the interaction strength

We assume that the distribution function of the particles is given by a thermal state which can be seen as the stationary state of the system reached after sufficiently long times. This state is given by

with single-particle kinetic energy

and partition function . Note that and are explicitly time-dependent. The expectation value is here defined by with integrals and and for an arbitrary function of the atomic positions and momenta.

We assume a time dependent and, in particular, that the temporal change of is sufficiently slow such that the particles remain in a thermal state. With this assumption, a dynamical equation for the kinetic energy can be derived in the following.

Using Equations 5a, b and Equation 7, we obtain the dynamical evolution of the single-particle kinetic energy

Furthermore, we may write

with

Then, using , we can rewrite Equation 8 as

and the integration of this equation leads towhere with . The latter can be solved using integration by parts to obtainusing the notation . Defining and performing the limit with we obtain the resultwhere is the th modified Bessel function and describes the stable solution of the equationSee Supplementary Appendix A for a detailed derivation. The value for calculated from Equation 10 is the mean magnetization of the particles for the given value of . Thus, Equation 9 connects the magnetization before and after the ramp with the kinetic energy before and after the ramp. We now discuss how this result can be used to lower the kinetic energy of the particles.

In Figure 2a we plotted as a function of . It can be seen that is 0 for (paramagnetic phase) and increases for (ferromagnetic phase) while it tends to 1 in . Note that only the positive solution has been shown, but there is also the solution . This transition from spatially homogeneous (paramagnetic) to self-organized (ferromagnetic) has been described as a “phase transition.”

FIGURE 2

In the spatially homogeneous phase, for , the quotient of the kinetic energies in Equation 9 is always 1. This implies that any adiabatic change within the spatially homogeneous region will, to good approximation, not affect the kinetic energy.

However, when we assume that the coupling strength is initialized such that the particles are in the self-organized phase and ramped to a value where the particles are distributed spatially homogeneously—that is, and —we obtain

This result of the right-hand side of Equation 11 in Figure 2b is shown as black solid line. It is a monotonous decreasing function with . Therefore, we conclude that a potentially very low kinetic energy can be reached by starting the ramp from a high coupling strength . In this regime, for , we obtain the asymptotic resultwhere is the Euler number. This shows that the ratio of the kinetic energies is proportional to and is plotted as a dashed gray line in Figure 2b.

We now discuss how this principle might be applicable to the driven-dissipative dynamics of particles in a cavity.

3 Cooling protocol

In order to apply the results of the previous section, we must first analyze the dissipative effects in the particles’ dynamics. This is done by comparing the results of Equations 5a, b that discard any dissipative effects with the dynamics, including dissipation, in Equations 1a, b, c, d.

3.1 Effects of dissipation

We first initialize the particles in a strongly self-organized thermal state with kinetic energy Then, we ramp exponentially as

for different ramping times . While the choice of an exponential ramp is a technical detail, it allows for a rather fast change of the interaction strength for large values of and slow changes for . This seems a good compromise between being fast and remaining approximately adiabatic. We study the dynamics of the full system, including the cavity degrees of freedom (Equations 1a, b, c, d) and the dynamics where the cavity degrees of freedom are eliminated (Equations 5a, b). Figures 3a, b show the dynamics of the kinetic energy following the ramp for a ramping time of and , respectively. To put these values into actual numbers, we provide an explicit example and use the value from Hosseini et al. (2017) for ¹³³Cs, which results in ramping times and , respectively. The dashed line shows the result using the conservative dynamics (Equation 5a, b), while the solid line represents the full dissipative dynamics (Equation 1a, b, c, d). Both curves for both ramping times show decreased kinetic energy. There is good agreement of both dynamics on short timescales while we observe discrepancies for the longer ramping time. Therefore, we expect that for sufficiently short times, dissipative effects are still negligible while they affect the dynamics on longer timescales. This relies on a timescale separation of dissipative and conservative forces that relies on (i) the number of particles and (ii) the typical timescale separation of motion and cavity relaxation—that is, . This was studied by Jäger et al., (2016) and Schütz et al., (2016) and observed in Wu et al. (2023). In conclusion, this preliminary analysis demonstrates that there must be an optimal ramping time for which the lowest possible temperature can be achieved.

FIGURE 3

While our original assumption was that the ramp is close to adiabatic, we expect this assumption to fail, especially because the system parameters are ramped across a phase transition. An observable to test this is the kurtosis

The kurtosis is for a Gaussian state and deviates from 3 for non-Gaussian states. Figures 3c, d show the kurtosis with the same labeling for the two different ramping times. In Figure 3c, we observe that the kurtosis remains close to 3 for times , while it deviates for longer times as soon as the value of crosses the phase transition line. In Figure 3d we observe the same for the simulations of the conservative dynamics (Equations 5a, b) and times , while the simulation of the full dissipative dynamics (Equations 1a, b, c, d) shows values of on much shorter timescales. This finding supports our claim that the dynamics do not remain adiabatic across the phase transition. In addition, the discrepancies between the conservative and dissipative dynamics predict a dissipation-assisted creation of non-Gaussian states that is closely related to the dissipation-assisted stabilization of non-Gaussian states predicted in Schütz et al. (2016).

We now analyze the dependence of the minimum achievable temperature on the ramping time . Figure 4 compares the values of the final kinetic energies for different ramping times , different particle numbers, and different ratios of . The black line with symbols are calculated using simulations of Equations 1a, b, c, d with (circles), (crosses), and (pluses), where we show in Figure 4a and in Figure 4b. These are realistic values, and the lower value of is close to that realized in Hosseini et al. (2017). For both simulations, the system has been initialized with a pumping strength of and . The results of the simulations predict a local minimum of the kinetic energy in the range , thus showing that there is an optimal ramping time. In general, this optimal ramping time is shorter for smaller particle numbers. In addition, the minimum achievable kinetic energy is larger for smaller particle numbers. This should be due to the timescale separation between the conservative and dissipative forces that becomes larger for increasing particle numbers.

FIGURE 4

We also find that the result of the kinetic energy for (b) appears to be slightly displaced to larger ramping times with respect to the simulations for (a). We expect that this is due to a violation of the adiabaticity criteria, and , that are not fulfilled for smaller values of and and short ramping times in Figure 4b. For completeness, we have also included a simulation of Equations 5a and b that does not include dissipation and noise. The results are visible as gray dashed lines in Figure 4. Those curves are monotonically decreasing, thus showing that noise and dissipation are the origins for the local minima in the kinetic energy in the full simulations. The theoretical minimum of the achievable kinetic energy is shown as a gray dashed–dotted line. It is calculated using Equation 12 and , resulting in . We observe that the simulation without dissipation and noise (gray dashed line) converges to this theoretical minimum in the limit .

3.2 Cooling protocol

We now use this gained insight to minimize the kinetic energy of particles that are initially in a spatially homogeneous configuration. We thus assume that the initial state is a thermal state with temperature that can be reached by cavity cooling for . The actual choice of this state is rather arbitrary but should be sufficiently cold such that the ensemble can be cavity (laser) cooled.

In a first stage, we perform a quench in the driving-laser intensity determined by such that has a value for the system to reach a state well inside the self-organized phase. Over a very long time, the system again reaches a stationary state which is thermal. To be consistent, for large laser intensities, we need to take corrections of the final temperature into account which come from the laser driving power. This final kinetic energy is given in the well-organized regime (Niedenzu et al., 2011; Grießer et al., 2012) bywhere is the effective trapping frequency.

In a second stage, we consider a ramp from back to a value close to 0 such that both magnetization and kinetic energy are adiabatically reduced. If we assume that the system remains adiabatic in a thermal state, the optimum final kinetic energy can be approximated using Equation 12 by

By minimizing Equation 16 with respect to , we find

at an optimum value of

This minimum kinetic energy is of the order of the recoil energy .

We now consider simulations that test this prediction. Following the procedure of the first stage, Figure 5a shows the dynamics of the kinetic energy after a quench from to (Equation 18). Initially, a rapid increase is observed in the magnetization determined by . This can be seen in Figure 5c, which plots the cavity field determined by Equation 6. Since energy is conserved on short timescales, the kinetic energy is also exponentially increasing, reaching a maximum of . On longer timescales, dissipation guides the system toward a thermal state with a temperature given by Equation 15 (horizontal gray dashed line in Figure 5a) with corresponding magnetization (horizontal gray dashed line in Figure 5c). The steady-state magnetization has been calculated using where is the solution of Equation 10. Figure 5e plots the kurtosis (Equation 14) that starts and ends at a value close to , suggesting that both the final and initial state are thermal.

FIGURE 5

Following the second stage, we ramp the coupling strength according to Equation 13 with a ramping time of that we have found to be close to optimal for the choice of parameters in the previous subsection. Figure 5b shows the decrease of the kinetic energy that eventually reaches a value that is of the order of the recoil energy. We observe a final kinetic energy of , whereby Equation 17 predicts a similar value of . During this process, the field intensity and mean magnetization of the system are decreased (Figure 5d), reminiscent of adiabatic demagnetization. We emphasize that the adiabatic process is much faster than cavity cooling in the first stage. This is visible by comparing the time axes in Figures 5a and b. This difference in the timescales comes from the fact that the adiabatic process is collective. The kurtosis, visible in Figure 5f, remains 3 during the ramp until the phase transition line is crossed and the process is no longer adiabatic.

4 Conclusion

In this paper, we have studied the possibility of cooling transversally driven particles inside an optical cavity using a combination of cavity cooling and a protocol reminiscent of adiabatic demagnetization. To analyze the effect of dissipation, we have performed simulations of dissipative and conservative dynamical models for this physical setup. We have shown that the particles can reach kinetic energies comparable to the recoil limit for the parameter choice below the typical limit of cavity cooling. To achieve this final kinetic energy, we have tuned the laser power from a value well above the self-organization threshold to below it. The duration of this ramp is chosen sufficiently long such that it seems to be quasi-adiabatic for the coherent dynamics but sufficiently rapid such that dissipation has only a minor effect on the final kinetic energy.

While the results presented here rely on adiabatically changing the coupling or interaction strength that results in a change of the internal magnetization, we expect that similar physics can be achieved by changing an additional external field that simulates an effective magnetic field. This can, for instance, be accomplished by modulating a laser that directly drives the cavity beside the transversal laser field (Niedenzu et al., 2013).

It is important to emphasize that the cooling stage is crucial to achieving the final kinetic energy as it reduces entropy. However, the precise cooling protocol is rather arbitrary. In particular, the choice of cavity cooling, visible in Figure 5, can be replaced by other schemes, such as a ramp of the interaction strength instead of a quench, or even by other laser cooling mechanisms. Importantly, the sole outcome of this first stage is the preparation of a sufficiently cold and highly self-organized (magnetized) particle ensemble.

Regarding the ultimate limits of this cooling protocol, this analysis is performed with semiclassical equations. This means our approach is only valid for kinetic energies that are above the recoil limit. In addition, we have not included the quantum statistics of the particles, which becomes relevant for low temperatures. Including the latter would be an interesting extension of our work since one might expect different distributions for bosons (Baumann et al., 2010) and fermions (Helson et al., 2023; Zhang et al., 2021). In future research, it might also be interesting to use multi-mode cavities that provide more possibilities to tune the interactions and dissipation (Torggler and Ritsch, 2014; Keller et al., 2017; Keller et al., 2018). The study of such systems is not only interesting for advances in laser and cavity cooling but also as a simulator for classical and quantum thermodynamics (Vinjanampathy and Anders, 2016; Niedenzu et al., 2018). In conclusion, the engineering of interactions and dissipation for particles in optical cavities is a versatile tool for quantum technologies and studying new physics.

References

• 1

  AcínA.BlochI.BuhrmanH.CalarcoT.EichlerC.EisertJ.et al (2018). The quantum technologies roadmap: a european community view. New J. Phys.20, 080201. 10.1088/1367-2630/aad1ea

  • CrossRef
  • Google Scholar

• 2

  ArnoldK. J.BadenM. P.BarrettM. D. (2012). Self-organization threshold scaling for thermal atoms coupled to a cavity. Phys. Rev. Lett.109, 153002. 10.1103/physrevlett.109.153002

  • CrossRef
  • Google Scholar

• 3

  AsbóthJ. K.DomokosP.RitschH.VukicsA. (2005). Self-organization of atoms in a cavity field: threshold, bistability, and scaling laws. Phys. Rev. A72, 053417. 10.1103/physreva.72.053417

  • CrossRef
  • Google Scholar

• 4

  AsenbaumP.KuhnS.NimmrichterS.SezerU.ArndtM. (2013). Cavity cooling of free silicon nanoparticles in high vacuum. Nat. Commun.4, 2743. 10.1038/ncomms3743

  • CrossRef
  • Google Scholar

• 5

  BarzanjehS.XuerebA.GröblacherS.PaternostroM.RegalC. A.WeigE. M. (2022). Optomechanics for quantum technologies. Nat. Phys.18, 15–24. 10.1038/s41567-021-01402-0

  • CrossRef
  • Google Scholar

• 6

  BaumannK.GuerlinC.BrenneckeF.EsslingerT. (2010). Dicke quantum phase transition with a superfluid gas in an optical cavity. Nature464, 1301–1306. 10.1038/nature09009

  • CrossRef
  • Google Scholar

• 7

  BlackA. T.ChanH. W.VuletićV. (2003). Observation of collective friction forces due to spatial self-organization of atoms: from Rayleigh to bragg scattering. Phys. Rev. Lett.91, 203001. 10.1103/physrevlett.91.203001

  • CrossRef
  • Google Scholar

• 8

  ChuS. (1998). Nobel lecture: the manipulation of neutral particles. Rev. Mod. Phys.70, 685–706. 10.1103/revmodphys.70.685

  • CrossRef
  • Google Scholar

• 9

  Cohen-TannoudjiC. N. (1998). Nobel lecture: manipulating atoms with photons. Rev. Mod. Phys.70, 707–719. 10.1103/revmodphys.70.707

  • CrossRef
  • Google Scholar

• 10

  DelićU.ReisenbauerM.GrassD.KieselN.VuletićV.AspelmeyerM. (2019). Cavity cooling of a levitated nanosphere by coherent scattering. Phys. Rev. Lett.122, 123602. 10.1103/physrevlett.122.123602

  • CrossRef
  • Google Scholar

• 11

  DePueM. T.McCormickC.WinotoS. L.OliverS.WeissD. S. (1999). Unity occupation of sites in a 3d optical lattice. Phys. Rev. Lett.82, 2262–2265. 10.1103/physrevlett.82.2262

  • CrossRef
  • Google Scholar

• 12

  DomokosP.HorakP.RitschH. (2001). Semiclassical theory of cavity-assisted atom cooling. J. Phys. B34, 187–198. 10.1088/0953-4075/34/2/306

  • CrossRef
  • Google Scholar

• 13

  DomokosP.RitschH. (2002). Collective cooling and self-organization of atoms in a cavity. Phys. Rev. Lett.89, 253003. 10.1103/physrevlett.89.253003

  • CrossRef
  • Google Scholar

• 14

  GrießerT.NiedenzuW.RitschH. (2012). Cooperative self-organization and sympathetic cooling of a multispecies gas in a cavity. New J. Phys.14, 053031. 10.1088/1367-2630/14/5/053031

  • CrossRef
  • Google Scholar

• 15

  HelsonV.ZwettlerT.MivehvarF.ColellaE.RouxK.KonishiH.et al (2023). Density-wave ordering in a unitary fermi gas with photon-mediated interactions. Nature618, 716–720. 10.1038/s41586-023-06018-3

  • CrossRef
  • Google Scholar

• 16

  HorakP.HechenblaiknerG.GheriK. M.StecherH.RitschH. (1997). Cavity-induced atom cooling in the strong coupling regime. Phys. Rev. Lett.79, 4974–4977. 10.1103/physrevlett.79.4974

  • CrossRef
  • Google Scholar

• 17

  HosseiniM.DuanY.BeckK. M.ChenY.-T.VuletićV. (2017). Cavity cooling of many atoms. Phys. Rev. Lett.118, 183601. 10.1103/physrevlett.118.183601

  • CrossRef
  • Google Scholar

• 18

  HuJ.UrvoyA.VendeiroZ.CrépelV.ChenW.VuletićV. (2017). Creation of a bose-condensed gas of ⁸⁷rb by laser cooling. Science358, 1078–1080. 10.1126/science.aan5614

  • CrossRef
  • Google Scholar

• 19

  JägerS. B.SchützS.MorigiG. (2016). Mean-field theory of atomic self-organization in optical cavities. Phys. Rev. A94, 023807. 10.1103/physreva.94.023807

  • CrossRef
  • Google Scholar

• 20

  KellerT.JägerS. B.MorigiG. (2017). Phases of cold atoms interacting via photon-mediated long-range forces. J. Stat. Mech.2017, 064002. 10.1088/1742-5468/aa71d7

  • CrossRef
  • Google Scholar

• 21

  KellerT.TorgglerV.JägerS. B.SchützS.RitschH.MorigiG. (2018). Quenches across the self-organization transition in multimode cavities. New J. Phys.20, 025004. 10.1088/1367-2630/aaa161

  • CrossRef
  • Google Scholar

• 22

  LaddT. D.JelezkoF.LaflammeR.NakamuraY.MonroeC.O’BrienJ. L. (2010). Quantum computers. Nature464, 45–53. 10.1038/nature08812

  • CrossRef
  • Google Scholar

• 23

  MaunzP.PuppeT.SchusterI.SyassenN.PinkseP. W. H.RempeG. (2004). Cavity cooling of a single atom. Nature428, 50–52. 10.1038/nature02387

  • CrossRef
  • Google Scholar

• 24

  MetcalfH. J.Van der StratenP. (1999). Laser cooling and trapping. Springer Science & Business Media.

  • Google Scholar

• 25

  MorigiG.PinkseP. W. H.KowalewskiM.de Vivie-RiedleR. (2007). Cavity cooling of internal molecular motion. Phys. Rev. Lett.99, 073001. 10.1103/physrevlett.99.073001

  • CrossRef
  • Google Scholar

• 26

  NiedenzuW.GrießerT.RitschH. (2011). Kinetic theory of cavity cooling and self-organisation of a cold gas. Europhys. Lett.96, 43001. 10.1209/0295-5075/96/43001

  • CrossRef
  • Google Scholar

• 27

  NiedenzuW.MukherjeeV.GhoshA.KofmanA. G.KurizkiG. (2018). Quantum engine efficiency bound beyond the second law of thermodynamics. Nat. Commun.9, 165. 10.1038/s41467-017-01991-6

  • CrossRef
  • Google Scholar

• 28

  NiedenzuW.SchützS.HabibianH.MorigiG.RitschH. (2013). Seeding patterns for self-organization of photons and atoms. Phys. Rev. A88, 033830. 10.1103/physreva.88.033830

  • CrossRef
  • Google Scholar

• 29

  PhillipsW. D. (1998). Nobel lecture: laser cooling and trapping of neutral atoms. Rev. Mod. Phys.70, 721–741. 10.1103/revmodphys.70.721

  • CrossRef
  • Google Scholar

• 30

  Schleier-SmithM. H.LerouxI. D.ZhangH.Van CampM. A.VuletićV. (2011). Optomechanical cavity cooling of an atomic ensemble. Phys. Rev. Lett.107, 143005. 10.1103/physrevlett.107.143005

  • CrossRef
  • Google Scholar

• 31

  SchützS.HabibianH.MorigiG. (2013). Cooling of atomic ensembles in optical cavities: semiclassical limit. Phys. Rev. A88, 033427. 10.1103/physreva.88.033427

  • CrossRef
  • Google Scholar

• 32

  SchützS.JägerS. B.MorigiG. (2015). Thermodynamics and dynamics of atomic self-organization in an optical cavity. Phys. Rev. A92, 063808. 10.1103/physreva.92.063808

  • CrossRef
  • Google Scholar

• 33

  SchützS.JägerS. B.MorigiG. (2016). Dissipation-assisted prethermalization in long-range interacting atomic ensembles. Phys. Rev. Lett.117, 083001. 10.1103/physrevlett.117.083001

  • CrossRef
  • Google Scholar

• 34

  SchützS.MorigiG. (2014). Prethermalization of atoms due to photon-mediated long-range interactions. Phys. Rev. Lett.113, 203002. 10.1103/physrevlett.113.203002

  • CrossRef
  • Google Scholar

• 35

  StenholmS. (1986). The semiclassical theory of laser cooling. Rev. Mod. Phys.58, 699–739. 10.1103/revmodphys.58.699

  • CrossRef
  • Google Scholar

• 36

  TishinA.SpichkinY. (2003). The magnetocaloric effect and its applications. 1st ed.Boca Raton: CRC Press.

  • Google Scholar

• 37

  TorgglerV.RitschH. (2014). Adaptive multifrequency light collection by self-ordered mobile scatterers in optical resonators. Optica1, 336. 10.1364/optica.1.000336

  • CrossRef
  • Google Scholar

• 38

  VinjanampathyS.AndersJ. (2016). Quantum thermodynamics. Contemp. Phys.57, 545–579. 10.1080/00107514.2016.1201896

  • CrossRef
  • Google Scholar

• 39

  VuletićV.ChuS. (2000). Laser cooling of atoms, ions, or molecules by coherent scattering. Phys. Rev. Lett.84, 3787–3790. 10.1103/physrevlett.84.3787

  • CrossRef
  • Google Scholar

• 40

  WiemanC. E.PritchardD. E.WinelandD. J. (1999). Atom cooling, trapping, and quantum manipulation. Rev. Mod. Phys.71, S253–S262. 10.1103/revmodphys.71.s253

  • CrossRef
  • Google Scholar

• 41

  WinelandD. J.ItanoW. M. (1979). Laser cooling of atoms. Phys. Rev. A20, 1521–1540. 10.1103/physreva.20.1521

  • CrossRef
  • Google Scholar

• 42

  WolkeM.KlinnerJ.KeßlerH.HemmerichA. (2012). Cavity cooling below the recoil limit. Science337, 75–78. 10.1126/science.1219166

  • CrossRef
  • Google Scholar

• 43

  WuZ.FanJ.ZhangX.QiJ.WuH. (2023). Signatures of prethermalization in a quenched cavity-mediated long-range interacting fermi gas. Phys. Rev. Lett.131, 243401. 10.1103/physrevlett.131.243401

  • CrossRef
  • Google Scholar

• 44

  ZhangX.ChenY.WuZ.WangJ.FanJ.DengS.et al (2021). Observation of a superradiant quantum phase transition in an intracavity degenerate fermi gas. Science373, 1359–1362. 10.1126/science.abd4385

  • CrossRef
  • Google Scholar