Next Article in Journal
Optical Lines of Ru21+ to Ru24+ Ions
Next Article in Special Issue
Rydberg Interaction-Induced Distortion of the Autler–Townes Spectra in Cold Lithium Atoms
Previous Article in Journal
Heavy Ion Beam Probing Diagnostics on the TUMAN-3M Tokamak for Study Plasma Potential and Electric Fields in New Operational Regimes
Previous Article in Special Issue
Nonlinear Dynamics in Isotropic and Anisotropic Magneto-Optical Traps
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonequilibrium Steady State in a Large Magneto-Optical Trap

1
Université Côte d’Azur, CNRS, Institut de Physique de Nice, 06560 Valbonne, France
2
Department of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA
3
MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
4
Center for Complex Quantum Systems, Department of Physics and Astronomy, Aarhus University, Aarhus C, 8000 Aarhus, Denmark
*
Author to whom correspondence should be addressed.
Atoms 2022, 10(4), 153; https://doi.org/10.3390/atoms10040153
Submission received: 6 November 2022 / Revised: 29 November 2022 / Accepted: 8 December 2022 / Published: 15 December 2022
(This article belongs to the Special Issue Cold and Rydberg Atoms for Quantum Technologies)

Abstract

:
Considering light-mediated long-range interactions between cold atoms in a magneto-optical trap (MOT), we present numerical evidence of a nonequilibrium steady state (NESS) for sufficiently large number of atoms (> 10 8 ). This state manifests itself as the appearance of an anisotropic distribution of velocity when a MOT approaches the threshold beyond which self-oscillating instabilities occur. Our three-dimensional (3D) spatiotemporal model with nonlocal spatial dependencies stemming from the interatomic interactions has recently been compared successfully to predict different instability thresholds and regimes in experiments with rubidium atoms. The behavior of the NESS is studied as a function of the main MOT parameters, including its spatiotemporal characteristics.

1. Introduction

A magneto-optical trap (MOT) is known worldwide as the staple technology for producing cold and confined ensembles of neutral atoms. Conceptually simple, this device is essential in the preparation stages of highly involved tasks, e.g., Rydberg atom quantum computing [1,2] and Bose-Einstein condensation [3,4]. Additionally, the limit of a large atom number N presents a wide array of interesting physics, such as random lasing [5,6], Anderson localization [7,8], subradiance [9,10,11], superradiance [12,13,14] as well as spontaneous oscillations that bear similarities to plasma and stellar phenomena [15,16]. The last example falls into the category of MOT instabilities, which are, by their dynamic nature, out of thermal equilibrium. In the stable MOT regime, however, the thermodynamical considerations can be complex, requiring handling with particular care.
Stable MOTs governed by single-atom physics normally possess isotropic Gaussian velocity distributions and, thus, are deemed in thermal equilibrium. Exceptions occur in cases with Sisyphus cooling, where the distributions resemble double Gaussians [17]. For conservative traps (such as magnetic traps or optical dipole traps, used in, e.g., Bose-Einstein condensation), the evaporative cooling processes depend on short-range collisional interactions (in particular, van der Waals interactions); after few collisions, the atoms thermalize and acquire isotropic Gaussian velocity distributions, while anisotropically trapped [18]. In the case studied in this paper, light-mediated long-range interactions are considered in a stable MOT containing many atoms. Compared to the case with short-range interactions, where detailed analytical results of the velocity distribution can, in principle, be obtained [18], the collective forces stemming from long-range interactions (see below) pose a different challenge, making it unclear whether analytical modelling is viable. We invoke here numerical simulations and predict a velocity distribution typical of a nonequilbrium steady state (NESS), which we observe to act as a precursor to the unstable (self-oscillating) regime. In general, such a state can appear in a non-periodically driven system with broken reversibility and is characterized by the existence of nonzero net currents and a nonequilibrium probability distribution [19]. States of this kind are abundant in nature, existing in processes as distant as biological ones [20,21].
The many-atom physics in a MOT become important typically for N > 10 4 atoms [22]. The collective forces that appear include, e.g., the shadow force [23] and the rescattering force [22]. The former one is compressive and is due to an imbalance of the beam intensities in the cloud, caused by their attenuation due to the light’s scattering. The latter force is, on the opposite hand, repulsive and appears as the atoms rescatter the scattered photons. These antagonistic forces are critical for understanding different features induced by multiple-scattering, e.g., the stable-cloud size increase with N [22,24] as well as the occurrence of spontaneous oscillations when the beam detuning is brought close to the atomic resonance [15,25].
This article’s main objective is to numerically study the three-dimensional (3D) particle velocity distributions of the stable regime, including how they are impacted as the unstable regime is approached. At large N, velocity anisotropy is observed, which is a signature of a NESS in a MOT. The NESS behavior is studied with respect to all of the main MOT parameters—the atom number N, the beam detuning δ , the magnetic field gradient B along the strong axis, and the intensity I of a single beam. Moreover, given the 3D nature of the numerical simulations, we are permitted to investigate the full spatiotemporal characteristics of the NESS.
This article continues as follows. In Section 2, the theoretical model employed in the 3D simulations is briefly described as well as the details surrounding its numerical implementation. Then, Section 3 presents the main results, including the identification of the NESS, its behavior versus the MOT parameters and its spatiotemporal characteristics. Finally, in Section 4, we conclude and discuss the future perspectives.

2. Theoretical Approach

Our employed theoretical model and its numerical implementation have been detailed in Ref. [26], and here we provide a brief description. We consider the so-called balanced MOT configuration, where the laser beams are independent and of the same, constant intensities before entering the cloud. The model is based on the hyperfine transition F = 0 F = 1 , with each of the three Zeeman transitions m = 0 m = 1 , 0 , + 1 treated as an independent 2-level system and driven by, respectively, σ , π , σ + polarized light (see Figure 1b in Ref. [26]). The magnetic field (which splits the Zeeman levels) is linearized under the assumption of the atom position being much smaller than the radius of the MOT coils and the separation between them. The main physical effects of the model are (i) the (Doppler) trapping force; (ii) the diffusion stemming from its fluctuations; (iii) the beam intensity attenuation due to the light’s scattering by the atoms; and (iv) the rescattering caused by the scattered photon exchange between the atoms. Working with a F = 0 F = 1 system (as opposed to a 2-level system) permits a proper description of the features related to the magnetic field and light polarization.
An important feature of our model is it being spatially nonlocal. Particularly, this is due to the atomic cross sections (in the expressions of the main effects) depending on the beam intensities, which are subject to attenuation. Note that the attenuation affects all of the main effects, including itself. In the case of the trapping force, the attenuation’s inclusion results in an additional compression, i.e., the shadow force (see Figure 2a in Ref. [26]), which is antagonistic to the rescattering (see Figure 2b in Ref. [26]). The attenuation naturally introduces the beam cross-saturation effect, which is computed using a devised numerical scheme.
Moreover, the modelling is done in 3D, without any assumption of spherical symmetry. As shown and discussed more in Section 3, this allows us to observe stable clouds in different non-centrosymmetric configurations.
The complete system dynamics in our model are described by the following collisionless Vlasov-type kinetic equation for the atomic phase-space density f ( r , v , t ) (refer to Equation (26) in Ref. [26]):
t f + v r f + 1 M v [ F t r ( r , v ) + F r s ( r , v ) ] f 1 M 2 2 v 2 D ( r , v ) f = 0
where M is the atomic mass, F t r is the trapping force, F r s is the rescattering force, and D is the momentum diffusion coefficient; the effects depend implicitly on time as well as on the MOT parameters mentioned earlier—N, δ , B , I. The dependence on N arises from their dependence on I that contains information on the cloud density (for F r s , each particle also contributes to greater repulsion), whereas the dependencies on δ and B arise as these are crucial for a proper description of cooling and confining processes. Finding direct numerical solutions to this equation can be deemed impractical considering both local and nonlocal spatial dependencies. Thus, a more sensible yet an equivalent way of solving the system’s dynamics is used by implementing a numerical produce based on a superparticle [27,28] treatment. As an important remark, note that the Vlasov equation is known to contain NESS solutions [29,30], meaning a MOT can potentially be found in such a state. In Section 3, we argue that a MOT indeed contains a NESS.
Our numerical implementation proceeds as follows. We start by generating a Gaussian cloud composed of 7 × 10 3 superparticles representing a larger amount of real particles (here, 10 4 to 1.5 × 10 10 ). The cloud then evolves under the action of the effects of the model. The analysis is done by randomly selecting cloud images (here, 200, unless specified otherwise) after the initial transient time (see Figure 3 in Ref. [25]).
Note that we here concentrate on the predictions of the stable MOT regime. The unstable (self-oscillating) one has been studied with our numerical model in previous works [16,25] showing qualitative agreements with the experimental observations. Nevertheless, the rationale behind picking particular MOT simulation-parameters is based on our observation, to be elaborated in Section 3, that the MOT contains a unique stable regime (NESS phase) before entering the unstable regime.

3. Main Results

In this section, the main results of our simulations are presented. The 3D velocity distributions of stable MOT clouds are studied, including how they are impacted as the unstable regime is approached. The study is done with respect to all of the main MOT parameters—N, δ , B , I—in this particular order. The spatiotemporal profiles of the 3D velocity distributions are also discussed. The analysis of our observations leads us to argue for the existence of a NESS in a MOT.
We begin by observing Figure 1, showing distributions for the three velocity dimensions at different values of N (the dotted lines are Gaussian fits to the histograms). The shape of the histograms is Gaussian due to the stochastic force (proportional to the diffusion coefficient) being of this kind. When above 10 8 atoms, the distributions broaden and, more remarkably, the distribution of the z axis (strong magnetic field axis) narrows compared to the remaining ones. We relate these observations to the growing antagonism between the shadow force and the rescattering force (many-atom effects) as N is increased. Note that different parts of the large clouds ( N > 10 8 ) possess varying degrees of velocity anisotropy and even have its sign switched when the atoms near the cloud edge are considered. This is illustrated in Figure 2, for N = 10 10 : The cloud core ( 70 % of all particles; r < 2 r R M S ) contributes positively to the anisotropy and is responsible for the enhanced peak in Figure 1, whereas the parts further away either do not contribute ( 2 r R M S < r < 2.25 r R M S ) or contribute negatively ( r > 2.25 r R M S ) and thus diminish the overall degree of anisotropy. The change observed in the distribution widths in all the velocity dimensions (Figure 2) is indicative of spatially inhomogeneous flows, as in, e.g., 3D vortexes. As discussed more when presenting the spatiotemporal velocity profiles, the appearance of a nonequilibrium/anisotropic velocity distribution is a signature of a NESS in our system and vortex-like flows can indeed form.
To proceed quantifying the velocity anisotropy (in whole clouds, as in Figure 1), we use the normalized difference of the root-mean-square (RMS) widths of the fitted Gaussians, σ x σ z σ x + σ z , where σ x σ y , σ z are the corresponding Gaussian RMS widths. Figure 3 shows the velocity anisotropy plotted versus N for two different δ values. Above 10 8 atoms, a transition occurs from isotropic to anisotropic velocity, and the anisotropy grows as N is increased but at a slower rate when | δ | is greater (at | δ | / Γ = 3.5 versus 3.2 ). The transition is correlated with the unstable regime occurring above 10 8 atoms, and the slower growth rate—that for a greater | δ | the instability threshold is farther [25]. Below we present a natural way of extracting the NESS threshold at a fixed N.
We continue by studying the NESS velocity anisotropy as the unstable regime is approached. In Figure 4, we display the behavior of the RMS widths σ α = x , y , z versus δ . As can be observed, the widths first decrease and then increase as | δ | gets lower, with the velocity remaining isotropic up to around the point, where they reach their minima. The NESS threshold occurs roughly at this point. We naturally mark it at the detuning beyond which the cloud shape diverges from ellipsoidal, to be discussed when presenting the NESS spatiotemporal characteristics. The decrease of the widths occurs due to the strengthening of the trap’s friction coefficient, whereas the increase is explained by the influence of the many-atom physics. Indeed, in the single-atom limit, we find a continuous decrease of the widths up to a detuning of Γ , beyond which an increase is observed, in agreement with the Doppler theory [31]. In Figure 4, the anisotropy is observed to grow larger as the instability threshold is approached, which is correlated with the cloud shape changes discussed below. Note that we have added the result for the widths for an unstable cloud close to the instability threshold. We observe in this case the velocity distributions to be Gaussian and accompanying greater anisotropy than in the stable regime. We report that the distributions in the unstable regime are generally not Gaussian (close to the threshold or beyond), as different instability regimes exist [16]. This is, however, beyond the scope of this paper and not discussed further.
In Figure 5, we plot the NESS threshold versus B (dashed line), superimposed with the instability threshold (solid line; taken from Figure 6 in Ref. [25]). The NESS threshold is seen to follow roughly the same qualitative behavior as the instability threshold, which decreases linearly with B , with the slope of 0.14 Γ per G/cm. The negative 1.5 Γ offset from the instability threshold coincides with the cloud changing its shape from ellipsoidal (far-detuned case) to a different non-centrosymmetric configuration. Approaching the instability threshold closer leads to greater deformations, which we discuss when presenting the spatiotemporal profiles of the 3D velocity distributions.
We consider next examining the impact of I on the NESS velocity anisotropy. In Figure 6, the results versus I are displayed for the NESS near the instability threshold ( 0.1 Γ away), i.e., where the anisotropy is close to being maximal (see the example Figure 4). The behavior is observed to be highly non-monotonous and separated into three parts: At low I ( 1 mW/cm 2 ) the anisotropy decreases, then increases at intermediate I (between 1 and 10 mW/cm 2 ) and, finally, decreases again at large I (≥10 mW/cm 2 ). This behavior bears some resemblance to the corresponding instability threshold behavior (see the inset), which we note has been qualitatively verified by experimental results (see Figure 3.34 in Ref. [32]). This resemblance indicates that the anisotropy becomes enhanced the closer the instability threshold is to the resonance. The behaviors are, however, not exactly identical: At small I the anisotropy is small ( 3 %) compared to the intermediate region anisotropy (growing from 1 to 16 %), whereas the thresholds mirror equidistantly in these regions (logarithmic scale), and at large I the drop in the anisotropy (toward 11 %) corresponds to the threshold getting plateaued. Due to the NESS phase’s presence being necessary for the MOT’s transition into the unstable regime, this downward trend in the anisotropy size may indicate that the instability mechanism eventually gets broken for larger intensities.
Finally, we discuss the spatiotemporal properties of the 3D velocity distributions of the stable clouds and visually distinguish between the regular and NESS clouds. In Figure 7, the velocity arrow plots are displayed for the clouds at different detunings from the instability threshold detuning δ t h r , i n (rows), viewed from the side (left column) and the top (right column). Far from the instability threshold ( δ t h r , i n δ = 2.4 Γ ), the cloud takes on a flattened ellipsoid shape, as expected for a MOT having the magnetic field stronger along one of the axes (the z axis); the velocities are clearly disordered as the particles are subject to the effect of diffusion. Whereas this cloud is outside the NESS region (see Figure 5), the other two are inside and possess noticeably different non-centrosymmetric configurations. The first one ( δ t h r , i n δ = 0.9 Γ ) appears triangular (side view) due to being pinched along the beam axes (top view; diagonal beam directions); the velocities start to acquire order as the particles move in a vortex-line manner around the lobes separated by the regions where the pinches occur. Closer to the instability threshold ( δ t h r , i n δ = 0.2 Γ ), the pinches become more pronounced and so does the order of the velocities. Figure 8 displays how the superparticle number in each octant is affected for the NESS region: Increasingly unequal amounts are found in adjacent octants when approaching the unstable regime, with the NESS threshold naturally set at the point before this deviation (at δ t h r , i n δ = 1.6 Γ for the considered gradient). We observe the shape changing to be universal as the instability threshold is approached and accompanying growing velocity anisotropy (see the example Figure 4). This change is indeed due to the growing antagonism between the shadow effect and the rescattering, as the unstable regime is not reached when either of the two effects are removed from the simulations [26]. On an equivalent note, the stable cloud takes a flattened ellipsoid shape when attenuation (causing the shadow effect and thus fuelling the antagonism) is removed (see Figure 3.45 in Ref. [32]).
The observation of the inhomogeneous particle-movement is important and supports the argument that a NESS is present. Indeed, as mentioned in Section 1, such a state can appear in non-periodically driven systems (e.g., our balanced MOT) and is characterized by the existence of nonzero net currents and a nonequilibrium probability distribution (which corresponds to velocity anisotropy [see Figure 1]). Note that although the particles move inhomogeneously, the envelopes of the NESS region clouds nevertheless exhibit as little fluctuation as the regular stable clouds (see Figure 3.27 in Ref. [32]), proving that the definition of cloud stability is respected.
As the last remark, there is a possibility of the NESS behavior being quasi-stationary, which otherwise has been observed in different many-body systems with long-range interactions [33,34]. In such a case, the NESS clouds would eventually evolve into other stable configurations. The timescales for that would, however, need to exceed 5 s, given this is the longest duration of our simulations. For comparison, such timescales are much longer than the damping period of MOT clouds, taking some tens of milliseconds.

4. Conclusions

In this paper, we numerically predicted the existence of a NESS in a MOT containing large amount of atoms ( > 10 8 ). This stable-MOT state was found to emerge before the MOT’s transition into the unstable (self-oscillating) regime. The NESS clouds appear in different non-centrosymmetric configurations, where the particles move in a vortex-line manner that becomes more pronounced and accompanying greater velocity anisotropy as the unstable regime is approached. This enhancement of the features results from the growth of antagonism between two collective forces—the shadow force (stemming from attenuation) and the rescattering force. Close to the instability threshold ( 0.1 Γ away), the velocity anisotropy exceeds 10 % for a wide range of beam intensities, making it possible to be detected in realistic experimental conditions. The anisotropy measurements can be further enhanced by probing the core of the clouds (refer to Figure 2). Recoil-induced resonance (RIR) spectroscopy could be employed in the velocity measurements (due to the large cloud size) [35,36], whereas the shape measurements could be achieved by cloud slicing and subsequent image reconstruction from the slices. A more systematic theoretical study of the vortex-like flows of particles might also be an interesting direction for future investigations. In this regard, conditions should arise for which the flows become turbulent, resulting in unstable cloud motion.

Author Contributions

Conceptualization, M.G. and R.K.; methodology, M.G. and R.K.; software, M.G.; data curation, M.G. and Y.-C.Z.; investigation, M.G.; validation, T.P., G.L. and R.K.; writing—original draft, M.G.; writing—review & editing, M.G., Y.-C.Z., G.L. and R.K.; visualization, M.G.; supervision, G.L. and R.K. All authors have read and agreed to the published version of the manuscript.

Funding

M.G., G.L. and R.K. acknowledge the financial support from the European training network ColOpt, which is funded by the European Union (EU) Horizon 2020 programme under the Marie Sklodowska-Curie action (Grant No. 721465). G.L. and R.K. further acknowledge the financial support from the European project ANDLICA, ERC advanced grant agreement (Grant No. 832219). Y.-C.Z. acknowledges the financial support from the National Key Research and Development Program Project of China (Grant No. 2021YFA1401700), the National Nature Science Foundation of China (Grant No. 12104359), and Xi’an Jiaotong University through the “Young Top Talents Support Plan” and the Basic Research Funding (Grant No. xtr042021012). T.P. acknowledges the financial support from the Danish National Research Foundation through the Niels Bohr Professorship and the Center of Excellence “CCQ” (Grant No. DNRF156).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saffman, M.; Walker, T.G.; Mølmer, K. Quantum information with Rydberg atoms. Rev. Mod. Phys. 2010, 82, 2313. [Google Scholar] [CrossRef] [Green Version]
  2. Isenhower, L.; Urban, E.; Zhang, X.L.; Gill, A.T.; Henage, T.; Johnson, T.A.; Walker, T.G.; Saffman, M. Demonstration of a neutral atom controlled-NOT quantum gate. Phys. Rev. Lett. 2010, 104, 010503. [Google Scholar] [CrossRef]
  3. Cornell, E.A.; Wieman, C.E. Nobel Lecture: Bose-Einstein condensation in a dilute gas, the first 70 years and some recent experiments. Rev. Mod. Phys. 2002, 74, 875. [Google Scholar] [CrossRef] [Green Version]
  4. Ketterle, W. Nobel lecture: When atoms behave as waves: Bose-Einstein condensation and the atom laser. Rev. Mod. Phys. 2002, 74, 1131. [Google Scholar] [CrossRef] [Green Version]
  5. Baudouin, Q.; Mercadier, N.; Guarrera, V.; Guerin, W.; Kaiser, R. A cold-atom random laser. Nat. Phys. 2013, 9, 357. [Google Scholar] [CrossRef] [Green Version]
  6. Gerasimov, L.V.; Kuprianov, D.V.; Havey, M.D. Random lasing in an inhomogeneous and disordered system of cold atoms. Opt. Spectrosc. 2015, 119, 377. [Google Scholar] [CrossRef] [Green Version]
  7. Cottier, F.; Cipris, A.; Bachelard, R.; Kaiser, R. Microscopic and macroscopic signatures of 3D Anderson localization of light. Phys. Rev. Lett. 2019, 123, 083401. [Google Scholar] [CrossRef] [Green Version]
  8. Cobus, L.A.; Maret, G.; Aubry, A. Crossover from renormalized to conventional diffusion near the three-dimensional Anderson localization transition for light. Phys. Rev. B 2022, 106, 014208. [Google Scholar] [CrossRef]
  9. Guerin, W.; Araújo, M.O.; Kaiser, R. Subradiance in a large cloud of cold atoms. Phys. Rev. Lett. 2016, 116, 083601. [Google Scholar] [CrossRef] [Green Version]
  10. Ferioli, G.; Glicenstein, A.; Henriet, L.; Ferrier-Barbut, I.; Browaeys, A. Storage and release of subradiant excitations in a dense atomic cloud. Phys. Rev. X 2021, 11, 021031. [Google Scholar] [CrossRef]
  11. Gold, D.C.; Huft, P.; Young, C.; Safari, A.; Walker, T.G.; Saffman, M.; Yavuz, D.D. Spatial coherence of light in collective spontaneous emission. PRX Quantum 2022, 3, 010338. [Google Scholar] [CrossRef]
  12. Araújo, M.O.; Krešić, I.; Kaiser, R.; Guerin, W. Superradiance in a large and dilute cloud of cold atoms in the linear-optics regime. Phys. Rev. Lett. 2016, 117, 073002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ferioli, G.; Glicenstein, A.; Robicheaux, F.; Sutherland, R.T.; Browaeys, A.; Ferrier-Barbut, I. Laser-driven superradiant ensembles of two-level atoms near Dicke regime. Phys. Rev. Lett. 2021, 127, 243602. [Google Scholar] [CrossRef]
  14. Pennetta, R.; Blaha, M.; Johnson, A.; Lechner, D.; Schneeweiss, P.; Volz, J.; Rauschenbeutel, A. Collective radiative dynamics of an ensemble of cold atoms coupled to an optical waveguide. Phys. Rev. Lett. 2022, 128, 073601. [Google Scholar] [CrossRef] [PubMed]
  15. Labeyrie, G.; Michaud, F.; Kaiser, R. Self-sustained oscillations in a large magneto-optical trap. Phys. Rev. Lett. 2006, 96, 023003. [Google Scholar] [CrossRef] [Green Version]
  16. Gaudesius, M.; Zhang, Y.-C.; Pohl, T.; Kaiser, R.; Labeyrie, G. Phase diagram of spatiotemporal instabilities in a large magneto-optical trap. Phys. Rev. A 2021, 103, L041101. [Google Scholar] [CrossRef]
  17. Jersblad, J.; Ellmann, H.; Støchkel, K.; Kastberg, A.; Sanchez-Palencia, L.; Kaiser, R. Non-Gaussian velocity distributions in optical lattices. Phys. Rev. A 2004, 69, 013410. [Google Scholar] [CrossRef] [Green Version]
  18. Bohn, J.L.; Jin, D.S. Differential scattering and rethermalization in ultracold dipolar gases. Phys. Rev. A 2014, 89, 022702. [Google Scholar] [CrossRef] [Green Version]
  19. Busiello, D.M.; Jarzynski, C.; Raz, O. Similarities and differences between non-equilibrium steady states and time-periodic driving in diffusive systems. New J. Phys. 2018, 20, 093015. [Google Scholar] [CrossRef]
  20. Qian, H. Open-system nonequilibrium steady state: Statistical thermodynamics, fluctuations, and chemical oscillations. J. Phys. Chem. B 2006, 110, 31. [Google Scholar] [CrossRef]
  21. Zhang, D.; Ouyang, Q. Nonequilibrium thermodynamics in biochemical systems and its application. Entropy 2021, 23, 271. [Google Scholar] [CrossRef] [PubMed]
  22. Sesko, D.W.; Walker, T.G.; Wieman, C.E. Behavior of neutral atoms in a spontaneous force trap. J. Opt. Soc. Am. B 1991, 8, 946. [Google Scholar] [CrossRef]
  23. Dalibard, J. Laser cooling of an optically thick gas: The simplest radiation pressure trap? Opt. Commun. 1988, 68, 203. [Google Scholar] [CrossRef]
  24. Camara, A.; Kaiser, R.; Labeyrie, G. Scaling behavior of a very large magneto-optical trap. Phys. Rev. A 2014, 90, 063404. [Google Scholar] [CrossRef] [Green Version]
  25. Gaudesius, M.; Kaiser, R.; Labeyrie, G.; Zhang, Y.-C.; Pohl, T. Instability threshold in a large balanced magneto-optical trap. Phys. Rev. A 2020, 101, 053626. [Google Scholar] [CrossRef]
  26. Gaudesius, M.; Zhang, Y.-C.; Pohl, T.; Kaiser, R.; Labeyrie, G. Three-dimensional simulations of spatiotemporal instabilities in a magneto-optical trap. Phys. Rev. A 2022, 105, 013112. [Google Scholar] [CrossRef]
  27. Birsdall, C.K.; Langdon, A.B. Plasma Physics via Computer Simulation; CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  28. Arber, T.D.; Bennett, K.; Brady, C.S.; Lawrence-Douglas, A.; Ramsay, M.G.; Sircombe, N.J.; Gillies, P.; Evans, R.G.; Schmitz, H.; Bell, A.R.; et al. Contemporary Particle-In-Cell approach to laser-plasma modelling. Plasma Phys. Control. Fusion 2015, 57, 113001. [Google Scholar] [CrossRef]
  29. Yamaguchi, Y.Y.; Barré, J.; Bouchet, F.; Dauxois, T.; Ruffo, S. Stability criteria of the Vlasov equation and quasi-stationary states of the HMF model. Phys. A Stat. Mech. App. 2004, 337, 36. [Google Scholar] [CrossRef] [Green Version]
  30. Farias, C.A.F.; Pakter, R.; Levin, Y. Entropy production and Vlasov equation for self-gravitating systems. J. Phys. A Math. Theor. 2018, 51, 494002. [Google Scholar] [CrossRef]
  31. Chang, R.; Hoendervanger, A.L.; Bouton, Q.; Fang, Y.; Klafka, T.; Audo, K.; Aspect, A.; Westbrook, C.I.; Clément, D. Three-dimensional laser cooling at the Doppler limit. Phys. Rev. A 2014, 90, 063407. [Google Scholar] [CrossRef]
  32. Gaudesius, M. Self-Oscillating Clouds in Magneto-Optical Traps. Ph.D. Thesis, Université Côte d’Azur, Nice, France, 2021. Available online: https://tel.archives-ouvertes.fr/tel-03273416 (accessed on 12 December 2022).
  33. Campa, A.; Dauxois, T.; Ruffo, S. Statistical mechanics and dynamics of solvable models with long-range interactions. Phys. Rep. 2009, 480, 57. [Google Scholar] [CrossRef] [Green Version]
  34. Gabrielli, A.; Joyce, M.; Marcos, B. Quasistationary states and the range of pair interactions. Phys. Rev. Lett. 2010, 105, 210602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Meacher, D.R.; Boiron, D.; Metcalf, H.; Salomon, C.; Grynberg, G. Method for velocimetry of cold atoms. Phys. Rev. A 1994, 50, R1992(R). [Google Scholar] [CrossRef] [PubMed]
  36. Vengalattore, M.; Prentiss, M. Recoil-induced resonances in the high-gain regime. Phys. Rev. A 2005, 72, 021401(R). [Google Scholar] [CrossRef]
Figure 1. Histograms (blue; bin width of 25 mm/s) for the three velocity dimensions at different values of N. The dotted lines are Gaussian fits. The appearance of a nonequilibrium/anisotropic velocity distribution ( v z ) above 10 8 atoms is a signature of a nonequilibrium steady state (NESS) in a stable magneto-optical trap (MOT; dashed lines are drawn near the v x distribution max). The remaining MOT parameters used in the simulations are the following: B = 3 G/cm, δ / Γ = 3.2 , I = 5 mW/cm 2 .
Figure 1. Histograms (blue; bin width of 25 mm/s) for the three velocity dimensions at different values of N. The dotted lines are Gaussian fits. The appearance of a nonequilibrium/anisotropic velocity distribution ( v z ) above 10 8 atoms is a signature of a nonequilibrium steady state (NESS) in a stable magneto-optical trap (MOT; dashed lines are drawn near the v x distribution max). The remaining MOT parameters used in the simulations are the following: B = 3 G/cm, δ / Γ = 3.2 , I = 5 mW/cm 2 .
Atoms 10 00153 g001
Figure 2. Histograms (green; bin width of 25 mm/s) for the three velocity dimensions for the bottom row cloud in Figure 1 ( N = 10 10 atoms), at 3 different volumes measured from the center, with r R M S referring to the root-mean-square radius. 1000 cloud images were used (instead of 200). The dotted lines are Gaussian fits, and the dashed lines are drawn near the v x distribution max. In the volume considered in the first row, 70 % of all particles are concentrated, whereas for the second and third rows the concentrations are 15 % each (hence smaller amount of bin counts). The particles in the first volume (the core particles) are responsible for the enhanced v z distribution peak in Figure 1, whereas the remaining ones (the edge particles) diminish it.
Figure 2. Histograms (green; bin width of 25 mm/s) for the three velocity dimensions for the bottom row cloud in Figure 1 ( N = 10 10 atoms), at 3 different volumes measured from the center, with r R M S referring to the root-mean-square radius. 1000 cloud images were used (instead of 200). The dotted lines are Gaussian fits, and the dashed lines are drawn near the v x distribution max. In the volume considered in the first row, 70 % of all particles are concentrated, whereas for the second and third rows the concentrations are 15 % each (hence smaller amount of bin counts). The particles in the first volume (the core particles) are responsible for the enhanced v z distribution peak in Figure 1, whereas the remaining ones (the edge particles) diminish it.
Atoms 10 00153 g002
Figure 3. Velocity anisotropy versus N for two different δ values. See the main text for the corresponding definition and the commentary on the displayed behaviors. Same remaining MOT simulation-parameters as in Figure 1.
Figure 3. Velocity anisotropy versus N for two different δ values. See the main text for the corresponding definition and the commentary on the displayed behaviors. Same remaining MOT simulation-parameters as in Figure 1.
Atoms 10 00153 g003
Figure 4. The root-mean-square (RMS) widths of the fitted Gaussians to velocity histograms, σ α = x , y , z , versus δ . The vertical dashed line represents the threshold beyond which the NESS behavior appears, and the vertical solid line represents the threshold beyond which the MOT becomes unstable (i.e., self-oscillations occur). The instability threshold is taken from Figure 6 in Ref. [25]. Same remaining MOT simulation-parameters as in Figure 1, but N = 1.5 × 10 10 atoms.
Figure 4. The root-mean-square (RMS) widths of the fitted Gaussians to velocity histograms, σ α = x , y , z , versus δ . The vertical dashed line represents the threshold beyond which the NESS behavior appears, and the vertical solid line represents the threshold beyond which the MOT becomes unstable (i.e., self-oscillations occur). The instability threshold is taken from Figure 6 in Ref. [25]. Same remaining MOT simulation-parameters as in Figure 1, but N = 1.5 × 10 10 atoms.
Atoms 10 00153 g004
Figure 5. NESS threshold versus B (dots, connected by a dashed line), superimposed with instability threshold (dots, connected by a solid line). The NESS clouds are contained between the two delineated thresholds, whereas the unstable ones are beyond the latter one. The regular stable clouds occur before the NESS clouds; unlike the latter, they possess isotropic velocity distributions and have ellipsoidal shapes (see text). The instability threshold results are taken from Figure 6 in Ref. [25] (note the same data points being blue circles). Same remaining MOT simulation-parameters as in Figure 4.
Figure 5. NESS threshold versus B (dots, connected by a dashed line), superimposed with instability threshold (dots, connected by a solid line). The NESS clouds are contained between the two delineated thresholds, whereas the unstable ones are beyond the latter one. The regular stable clouds occur before the NESS clouds; unlike the latter, they possess isotropic velocity distributions and have ellipsoidal shapes (see text). The instability threshold results are taken from Figure 6 in Ref. [25] (note the same data points being blue circles). Same remaining MOT simulation-parameters as in Figure 4.
Atoms 10 00153 g005
Figure 6. NESS velocity anisotropy near the instability threshold ( 0.1 Γ away) versus I, displayed in logarithmic scale. The circles are the data points. Same remaining MOT simulation-parameters as in Figure 4. The inset displays the corresponding instability threshold results (dots) showing a resembling behavior; the results for I 5 mW/cm 2 are taken from Figure 3.34 in Ref. [32].
Figure 6. NESS velocity anisotropy near the instability threshold ( 0.1 Γ away) versus I, displayed in logarithmic scale. The circles are the data points. Same remaining MOT simulation-parameters as in Figure 4. The inset displays the corresponding instability threshold results (dots) showing a resembling behavior; the results for I 5 mW/cm 2 are taken from Figure 3.34 in Ref. [32].
Atoms 10 00153 g006
Figure 7. (Color online.) Three-dimensional (3D) velocity arrow plots for the stable clouds at different detunings from the instability threshold detuning δ t h r , i n (rows), viewed from the side (left column) and the top (right column). Single-shot images are displayed. The x and y axes have been rotated by 45 between the side- and top-view images. The diagonals in the top-view images correspond to the directions of two pairs of MOT beams (dashed lines). The clouds in the bottom two rows posses 4 lobes (separated by the beam axes); in each side-view image the lobe 1 is obscured; each top-view image shows that the lobes 1 and 3 are lower in the z plane compared to the same cloud portions in the upper-row top-view image, whereas the lobes 2 and 4 are higher. The remaining MOT parameters used in the simulations are the following: B = 3.6 G/cm, I = 5 mW/cm 2 , N = 1.5 × 10 10 atoms. Note that the shape changing is universal as the instability threshold is approached and accompanies growing velocity anisotropy (see the example Figure 4).
Figure 7. (Color online.) Three-dimensional (3D) velocity arrow plots for the stable clouds at different detunings from the instability threshold detuning δ t h r , i n (rows), viewed from the side (left column) and the top (right column). Single-shot images are displayed. The x and y axes have been rotated by 45 between the side- and top-view images. The diagonals in the top-view images correspond to the directions of two pairs of MOT beams (dashed lines). The clouds in the bottom two rows posses 4 lobes (separated by the beam axes); in each side-view image the lobe 1 is obscured; each top-view image shows that the lobes 1 and 3 are lower in the z plane compared to the same cloud portions in the upper-row top-view image, whereas the lobes 2 and 4 are higher. The remaining MOT parameters used in the simulations are the following: B = 3.6 G/cm, I = 5 mW/cm 2 , N = 1.5 × 10 10 atoms. Note that the shape changing is universal as the instability threshold is approached and accompanies growing velocity anisotropy (see the example Figure 4).
Atoms 10 00153 g007
Figure 8. Evolution of the superparticle number in an octant, for the stable clouds in Figure 7 together with a few showing the transition to the NESS region (from δ t h r , i n δ = 1.6 to 1.4 Γ ). There are eight different-colored solid curves, each corresponding to one of the eight octants in 3D. When all the curves are on the dashed line at 12.5 % , the superparticle amounts in each octant are the same. Beyond the NESS threshold, the amounts in adjacent octants become increasingly unequal. Note that the initial growths in evolutions are transients, and the information contained in such plots is used to determine the NESS threshold.
Figure 8. Evolution of the superparticle number in an octant, for the stable clouds in Figure 7 together with a few showing the transition to the NESS region (from δ t h r , i n δ = 1.6 to 1.4 Γ ). There are eight different-colored solid curves, each corresponding to one of the eight octants in 3D. When all the curves are on the dashed line at 12.5 % , the superparticle amounts in each octant are the same. Beyond the NESS threshold, the amounts in adjacent octants become increasingly unequal. Note that the initial growths in evolutions are transients, and the information contained in such plots is used to determine the NESS threshold.
Atoms 10 00153 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gaudesius, M.; Zhang, Y.-C.; Pohl, T.; Labeyrie, G.; Kaiser, R. Nonequilibrium Steady State in a Large Magneto-Optical Trap. Atoms 2022, 10, 153. https://doi.org/10.3390/atoms10040153

AMA Style

Gaudesius M, Zhang Y-C, Pohl T, Labeyrie G, Kaiser R. Nonequilibrium Steady State in a Large Magneto-Optical Trap. Atoms. 2022; 10(4):153. https://doi.org/10.3390/atoms10040153

Chicago/Turabian Style

Gaudesius, Marius, Yong-Chang Zhang, Thomas Pohl, Guillaume Labeyrie, and Robin Kaiser. 2022. "Nonequilibrium Steady State in a Large Magneto-Optical Trap" Atoms 10, no. 4: 153. https://doi.org/10.3390/atoms10040153

APA Style

Gaudesius, M., Zhang, Y. -C., Pohl, T., Labeyrie, G., & Kaiser, R. (2022). Nonequilibrium Steady State in a Large Magneto-Optical Trap. Atoms, 10(4), 153. https://doi.org/10.3390/atoms10040153

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop