## Strongly correlated growth of Rydberg aggregates in a vapor cell

A. Urvoy, 1, ∗ F. Ripka, 1 I. Lesanovsky, 2 D. Booth, 3 J.P. Shaffer, 3 T. Pfau, 1 and R. L¨ ow 1, † 1 5. Physikalisches Institut and Center for Integrated Quantum Science and Technology, Universit¨t a Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany 2 School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, UK 3 Homer L. Dodge Department of Physics and Astronomy, The University of Oklahoma, 440 West Brooks Street, Norman, Oklahoma 73019, USA

(Dated: March 10, 2015)

The observation of strongly interacting many-body phenomena in atomic gases typically requires ultracold samples. Here we show that the strong interaction potentials between Rydberg atoms enable the observation of many-body effects in an atomic vapor, even at room temperature. We excite Rydberg atoms in cesium vapor and observe in real-time an out-of-equilibrium excitation dynamics that is consistent with an aggregation mechanism. The experimental observations show qualitative and quantitative agreement with a microscopic theoretical model. Numerical simulations reveal that the strongly correlated growth of the emerging aggregates is reminiscent of soft-matter type systems.

PACS numbers: 32.80.Ee, 34.20.Cf, 03.65.Yz, 61.43.Hv

Due to their exaggerated properties, Rydberg atoms find applications in various research fields ranging from cavity QED [1], quantum information [2], quantum optics [3-7], microwave sensing [8] to molecular physics [9]. One particular area of research employs the strong interactions between Rydberg atoms to create strongly interacting many-body quantum systems for quantum simulation [10-12], quantum phase transitions [13] and the realization of correlated or spatially ordered states [1418]. In any system, gaseous, liquid, glass or solid, spatial correlations can only arise if interactions are present. In Rydberg gases these correlations have recently been revealed by the direct imaging of the resonant excitation blockade effect [17, 19]. In our experiment an initially nearly ideal gas of thermal atoms at room temperature is excited into a strongly interacting Rydberg state. The Rydberg excitations show correlated many-body dynamics that shares similarities with that of soft-matter systems.

We consider the situation of a dense atomic gas, where dense means that the mean inter-particle distance is smaller than the typical interaction distance. The interplay between these strong interactions and an offresonant excitation laser results in a non-equilibrium dynamics with pronounced spatial correlations (related is the antiblockade in ultracold atoms [20, 21]). At a certain distance r fac from a Rydberg atom the interaction shifts the Rydberg states in resonance with the detuned excitation light fields, as depicted in Fig. 1(a). The subsequent resonant excitations are restricted to specific distances and therefore result in the formation of so-called aggregates that feature non-trivial spatial correlations [22-25]. This so-called facilitated excitation occurs predominantly at the boundary of aggregates [see Fig. 1(b)], similarly to seeded nucleation. The number of Rydberg excitations grows increasingly faster as the size of the aggregates gets larger. After a certain time a (quasi) steady state is reached and the excitation saturates because the sample is maximally filled with Rydberg atoms. The spatial correlations expected for this steady state are depicted in Fig. 1(b). The density-density correlation function sketched here shows an enhanced probability for finding two atoms separated by the facilitation radius r fac . There is no apparent long-range order since the isotropic facilitation mechanism does not favor the creation of crystalline arrangements. The existence of such Rydberg aggregates has recently been observed in ultracold systems [24, 25]. In these studies the Rydberg number distributions were extracted and matched to the results of empirical rate equation models, proving the existence of small aggregates with up to 8 excitations.

In this Letter we investigate a complementary regime and study the dynamics of the aggregation process in an atomic vapor at room temperature. Central parameters such as the dephasing rate and the atomic density are up to 2 orders of magnitude larger than in previous work [24, 25]. To quantify the non-equilibrium dynamics we perform a systematic study of the characteristic time scale that underlies the aggregation process. We find that the functional dependence of this time scale on various system parameters can be approximated by power laws. The corresponding exponents agree well with those obtained from a theoretical aggregation model. Further evidence for the aggregation dynamics is found through a quantitative evaluation of the number of Rydberg excitations and the resulting mean inter-particle distance which is consistent with the existence of a facilitation radius.

We perform our experiments in a gas of cesium atoms at room temperature or above. The cesium atoms are confined in a 220 µ m thick glass cell [see Fig. 2(b)]. The atomic density N g ranges from 10 µ m -3 to 500 µ m -3 ,

FIG. 1. Principle of the aggregation. (a) Interaction induced level shift. An already excited Rydberg atom (red dot) at the origin produces an energy shift for the neighboring atoms. At the facilitation radius r fac , the atoms are exactly shifted in resonance with the excitation laser, red-detuned by ∆. The grey shaded area symbolizes the excitation bandwidth, a Gaussian whose width is given by the dephasing rate γ . (b) Left: Typical time evolution of the Rydberg population and of the transmission change of the exciting laser [experimental signal, which is proportional to the derivative of the Rydberg density, see text and Fig. 2 for details] during the aggregation process. Right: Snapshots of the aggregation at times A, B and C, illustrated for a frozen gas in 2D. Red dots correspond to Rydberg atoms. Blue dots correspond to resonant atoms. The dark grey line shows the 'resonant shell' in which the excitation of Rydberg atoms is facilitated. Inset: Sketch of the density-density correlation function expected for the steady state [corresponding to (C)].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-9e68635c-f291-4c0e-9cd4-9b74f562ba4c.jpg)

which is large compared to previous experiments on interacting Rydberg atoms in thermal vapor [26-28] and nonBEC cold atomic gases [24, 25, 29] where N &lt; g ∼ 10 µ m -3 . We excite cesium atoms from the ground state 6S 1 / 2 to a Rydberg state n S 1 / 2 off-resonantly via the intermediate state 7P 3 / 2 [see Fig. 2(a)]. The two-photon transition is red-detuned by ∆ with respect to the unperturbed Rydberg state. The lower transition is driven by a cw laser at 455 nm, which also serves as the probing field. For the upper transition a laser at ∼ 1070 nm is pulsed on the nanosecond scale using a fast Pockels cell. The effective Rabi frequency of the two-photon transition is then given by Ω = Ω 455 Ω 1070 / (2 δ 7P ), where Ω 455 and Ω 1070 are the respective single transition Rabi frequencies and δ 7P = 1 5 GHz is the detuning to the intermedi-. ate state [30]. We measure the transmission change of the lower 455 nm laser after suddenly switching on the upper

FIG. 2. (a) Schematic energy level diagram for the excitation to Rydberg states of cesium atoms. (b) Sketch of the experimental setup. (c) Transmission signals at different atom number densities for the 32S state. The two-photon detuning is ∆ ≈ 2 π × -2200 MHz. The two-photon Rabi frequency is Ω = 2 π × 100 MHz. The grey-shaded area in the background represents the temporal pulse shape of the infrared laser. The time delay t sat used for the quantitative analysis is also depicted. The black curve (upper panel) is the expected transmission signal (not to scale) for non-interacting atoms (following the approach from Ref. [31]).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-3ee759c6-b67b-4001-a7dd-23b1856b0970.jpg)

1070 nm laser on. Note that this signal is directly proportional to the time derivative of the Rydberg population [30]. Typical transmission traces are shown in Fig. 2(c) for various atom number densities. Immediately after the 1070 nm laser is switched on (at t = 0) atoms are transferred to the Rydberg state, and for each excited Rydberg atom one photon is removed from the probe resulting in a decrease of the transmission. One striking feature is the large red detuning of up to 10 GHz at which these signals are observed, in contrast to previous work on coherent excitation in a more dilute gas [27, 31]. Line broadenings and shifts to the blue have been observed with S-states in cold atoms [24], and Rydberg aggregates have been observed only at blue detunings [24, 25], as expected from the repulsive Rydberg-Rydberg interaction potential for most Rydberg S-states [32]. However at large atomic densities, the mean inter-particle distance is so small that the perturbative character of the van-der-Waals potential is not valid anymore. Level crossings between

pair states now have to be considered. For Rydberg Sstates in cesium a strong resonant dipole-quadrupole interaction with neighboring pair-states occurs at 0 55 . µ m for n = 32 [30, 33]. Initially dipole-forbidden attractive pair states acquire some n Sn S admixture at short inter-atomic distances due to the dipole-quadrupole interaction. Therefore the Rydberg-Rydberg interaction for S-states in cesium also features an effective attractive component, allowing for us to observe an excitation signal at large red detunings.

A first indication for the collective nature of the excitation dynamics is provided by the data in Fig. 2. Here we consider the characteristic time scale of the excitation t sat , i.e. the time at which the transmission reaches its minimum [as shown in Fig. 2(c)]. In the non-interacting case this time scale is given by 1 / ∆ ∼ 0 1 ns (black line). . In the experiment however the excitation signal is slow ( t sat glyph[greatermuch] 1 / ∆) and density dependent, which is only possible if many-body effects are included. Pulse propagation effects can be excluded due to the small optical depth at the large detunings applied.

To analyze the dynamics, we systematically varied all accessible parameters independently: ∆, Ω, N g and n ∗ , where n ∗ = n -δ is the effective principal quantum number of the Rydberg state, δ being the quantum defect. Only the dephasing rate γ , determined by the Doppler shifts, the velocity of the atoms and the steep slope of the pair-state potentials, cannot easily be tuned over a significant range (for the derivation of the actual dephasing rate, see Supplemental Material [30]). The basic experimental settings are the following: ∆ ≈ 2 π ×-2200 MHz, Ω = 2 π × 100 MHz, N g = 36 µ m -3 and n ∗ = 27 95 . (32S state), and are all but one kept constant. First, by fitting the values of the time scale of the excitation t sat against ∆ we obtain a power law as t sat ∝ | ∆ | a with a = 1 99(14). . In the same way, we find that t sat ∝ Ω b with b = -2 10(5). . The power law for N g is t sat ∝ N c g with c = -1 09(6). . Finally we performed density scans at three different principal quantum numbers (32S, 34S and 36S states) and were able to extract a scaling behavior as N -1 g ( n ∗ ) d with d = -4 6(5). . These results are summarized in Table I.

By integrating the transmission signal over time we obtain an estimate of the actual number of Rydberg atoms present in the sample. At the time when the Rydberg excitation starts to saturate the Rydberg number in the excitation volume reaches approximately 50000 atoms. The corresponding mean inter-atomic distance is 〈 r 〉 = 0 65 . µ m, leading to approximative values for the effective van-der-Waals or dipole-dipole strength of interaction as C α = ∆ ×〈 〉 r α [22] with C 6 = 2 π ×-109 MHz · µ m and 6 C 3 = 2 π ×-405 MHz · µ m (at ∆ 3 ≈ 2 π ×-1500 MHz). We use these values as an input for simulations described in the following.

We compared these findings to the results of the Rydberg aggregation model adapted from Ref. [22]. This

|     | ∆ a        | Ω b         | N g c       | n ∗ d       |
|-----|------------|-------------|-------------|-------------|
| Exp | 1 . 99(14) | - 2 . 10(5) | - 1 . 09(6) | - 4 . 6(5)  |
| vdW | 2 . 15(1)  | - 2         | - 0 . 86(1) | - 4 . 72(5) |
| dd  | 2 . 38(0)  | - 2         | - 0 . 74(1) | - 2 . 95(2) |

a Range: -4.5 . . . -1.3 GHz

b Range: 75 . . . 175 MHz

c Range: 28 . . . 101 µ m - 3

d Range: 32S - 34S - 36S

TABLE I. Power laws for the aggregation delay t sat . In the experiment ('Exp'), the parameters are ∆ ≈ 2 π ×-2200 MHz, Ω = 2 π × 100 MHz, N g = 36 µ m -3 and n = 32, or varied over the specified range. The power laws labeled as 'vdW' and 'dd' were extracted from simulations in an ensemble of randomly and uniformly distributed atoms with van-der-Waals (ensemble of 10 3 atoms) dipoledipole interaction (ensemble of 15 3 atoms) respectively, using the parameters from the experiment and the extracted interaction strengths (see text). The dephasing rate was also extracted from the experiment [30]. All uncertainties represent the confidence region of the fits. No uncertainty is given for the exponent of Ω in the theory, fixed at -2 as by the assumption of strong dephasing.

model describes the dynamics of Rydberg aggregation within a system of two-level atoms in the presence of strong dephasing. This central assumption of strong dephasing applies to the experimental situation discussed here since the motion-induced dephasing is larger than the Rabi frequency ( γ &gt; ∼ 2 π × 0 5 GHz [30]). . The reduction to a two-level system is well justified in our case, although the commonly used adiabatic elimination of the intermediate state typically breaks down for quantitative considerations of the coherent evolution [27, 31]. Here however, because of the large dephasing rates due to the atomic motion we do not observe coherent evolution, and because of the large detunings we expect only less than 1 5% population of the intermediate state. . We performed simulations for an ensemble of randomly distributed atoms in 3D, assuming a pure van-der-Waals ( V vdW = C /r 6 6 ) or pure dipole-dipole ( V dd = C /r 3 3 ) attractive interaction potential. In fact the actual experimental potential curves have a more complex spatial dependence, due to strong state mixing [30]. The latter makes the actual microscopic description of the excitation dynamics rather intricate since a two-level description of an atom might not be sufficient for covering all possible pair excitation channels. Therefore the van-derWaals and dipole-dipole potential that are used in our two-level model should be regarded as the limiting approximations of the actual experimental situation. We evaluate the time scale of the excitation t sat defined by the transmission minimum, as in the experiment. All results are listed in Table I. We find an excellent agreement between the experimental and theoretical exponents. In particular the dependence of t sat on approximately Ω -2 indicates that a strong dephasing effect is present in our

FIG. 3. Density plots showing the transmission change as a function of time (vertical axis, from top to bottom) and detuning. The parameters are as in Table I. (a) Simulated data with pure van-der-Waals interaction. The Rabi frequency was rescaled with a factor 0 43. . (b) Simulated data with pure dipole-dipole interaction (for 15 3 atoms). The Rabi frequency was rescaled with a factor 0 26. . (c) Experimental data. The blue (resp. green) line shows the position of the absorption maximum t sat in the simulated data with van-der-Waals (resp. dipole-dipole) interaction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-1b6e061f-c286-41b1-aa84-1be89c570b94.jpg)

system [22]. Moreover the power laws vary only slightly, independently of whether pure van-der-Waals or pure dipole-dipole interaction are used. Our experimental results are compatible with both potentials, suggesting a rather weak dependence on the actual shape of the potential and that the basic mechanism of aggregation only relies on the existence of a facilitation radius.

In Fig. 3(a) [resp. (b)] we show the transmission change extracted from the model with van-der-Waals [resp. dipole-dipole] interaction and in Fig. 3(c) from the experiment. The quantitative agreement for the time scale t sat between theory and experiment is excellent. Due to the complex pair state potentials and strong statemixing that are crucial here, an adaptation of the Rabi frequency has to be performed in order to use the twolevel model (i.e. a model with a single Rydberg state). In our simulations the Rabi frequency is reduced by a factor of 0.38 (for van-der-Waals interaction) and 0.23 (for dipole-dipole interaction) compared to the experiment. These factors are compatible with the n Sn S admixture that one can extract from the pair-state potentials [30]. Note that with dipole-dipole interaction there is a dependence on the system size due to the long-range character of the 1 /r 3 potential, hence the different rescaling factors applied in Fig. 3(a) and 3(b).

The comparison of the experimental results to those of the model show that facilitated excitation is occurring, which implies that an aggregation dynamics takes place. The frozen gas approximation which is used to model the aggregation in [22] holds strictly for time intervals of about 0.3 ns, the transit time of the atoms through the 'resonant shell' (see Figure 1). During this time aggregates consisting of up to 4 atoms are formed [30]. Their spatial correlations will then quickly decay to those of a non-interacting ideal gas. Note, however, that even after that time the aggregation of excitations continues at the boundary of the uncorrelated aggregates.

Despite this effect the agreement between the experimental values of t sat and those obtained from the simple theoretical model is in general good [see Fig. 3]. Yet the time evolution differs as the experimental signal appears much more narrow in time than in the simulations. We attribute this discrepancy to the inherent differences between the microscopics of the experimental system and the theoretical model (e.g. the complex potentials and the atomic motion). Interestingly, when considering the temporal width of the signal, the simulations with dipole-dipole interaction seem to agree best with the experiment, which is realistic considering that one component of the actual interaction potentials is of pure dipole-dipole nature. Moreover, the model does not account for the fact that the excitation of isolated, i.e. not facilitated, atoms is in fact coherent.

Another feature of the experimental data that is not reproduced by the theory is the small overshoot in transmission just before the system reaches the steady state [for instance at 20 ns for N g = 88 µ m -3 in Fig. 2(a)]. Remnants of coherent Rabi oscillations can be ruled out given the large dephasing rates and large detunings with respect to the single-atom resonance. Moreover the time scale of the overshoot, e.g. the time between the two maxima, does not show a specific dependency on the detuning, as is the case for Rabi flopping. Explaining this feature will require further study and a more sophisticated theory. Note that ion effects can be neglected. There are no ions in the sample prior to the excitation, otherwise they would act as aggregation seeds [35] at t = 0 and change the signal dramatically. Moreover plasma formation would occur only once Rydberg atoms have been excited [34] and therefore cannot explain our Rydberg excitation signal. In spite of the low collisional ionization rate ( ∼ 40 MHz [34]) it might happen that in the course of the experiment single ions are created. However single ions will not disturb the aggregation dynamics, because of the similarity between the ion-Rydberg and the Rydberg-Rydberg potentials [35].

In conclusion, we have shown a facilitated excitation process of Rydberg atoms in thermal vapor. We have characterized the dynamics by measuring power laws for all relevant experimental parameters. These power laws, as well as absolute numbers for the dynamics and Rydberg density are in excellent agreement with those from a model for Rydberg aggregation. These results were obtained in the framework of thermal vapors, but the general ideas also apply to experiments in ultracold atoms. In the ultracold regime it would be interesting to study how quantum effects influence the aggregation time scale when the strength of dephasing noise is systematically lowered.

The authors would like to thank D. Peter, S. Weber, S. Hofferberth and J.P. Garrahan for fruitful discussions,

as well as H. K¨bler for proof reading. u This project was supported by the Carl-Zeiss-Stiftung and the ERC under contract number 267100. I.L. acknowledges funding from the European Research Council under the European Union's Seventh Framework Programme (FP/20072013) / ERC Grant Agreement n. 335266 (ESCQUMA) and the EU-FET Grant No. 512862 (HAIRS). D.B. and J.P.S. acknowledge financial support from grants NSF(PHY-1205392) and DARPA(60181-PH-DRP) for this work.

- ∗ Electronic address: a.urvoy@physik.uni-stuttgart.de
- † Electronic address: r.loew@physik.uni-stuttgart.de
- [1] S. Haroche, Rev. Mod. Phys. 85 , 1083 (2013).
- [2] M. Saffman, T. G. Walker, and K. Mølmer, Rev. Mod. Phys. 82 , 2313 (2010).
- [3] Y. O. Dudin and A. Kuzmich, Science 336 , 887 (2012).
- [4] T. Peyronel, O. Firstenberg, Q.-Y. Liang, S. Hofferberth, A. V. Gorshkov, T. Pohl, M. D. Lukin, and V. Vuleti´, c Nature 488 , 57 (2012).
- [5] D. Maxwell, D.J. Szwer, D. Paredes-Barato, H. Busche, J.D. Pritchard, A. Gauguet, K.J. Weatherill, M.P.A. Jones, and C.S. Adams, Phys. Rev. Lett. 110 , 103001 (2013).
- [6] H. Gorniaczyk, C. Tresp, J. Schmidt, H. Fedder, and S. Hofferberth, Phys. Rev. Lett. 113 , 053601 (2014).
- [7] D. Tiarks, S. Baur, K. Schneider, S. D¨ urr, and G. Rempe, Phys. Rev. Lett. 113 , 053602 (2014).
- [8] J. A. Sedlacek, A. Schwettmann, H. K¨ ubler, R. L¨ ow, T. Pfau, and J. P. Shaffer, Nature Physics 8 , 819 (2012).
- [9] V. Bendkowsky, B. Butscher, J. Nipper, J. P. Shaffer, R. L¨w, o and T. Pfau, Nature 458 , 1005 (2009).
- [10] H. Weimer, M. M¨ller, I. Lesanovsky, P. Zoller, u and H. P. B¨ uchler, Nature Physics 6 , 382 (2010).
- [11] G. G¨ unter, H. Schempp, M. Robert-de Saint-Vincent, V. Gavryusev, S. Helmrich, C. S. Hofmann, S. Whitlock, and M. Weidem¨ uller, Science 342 , 954 (2013).
- [12] S. Ravets, H. Labuhn, D. Barredo, L. B´guin, T. Lahaye, e and A. Browaeys, Nature Physics 10 , 914 (2014).
- [13] R. L¨w, H. Weimer, U. Krohn, R. Heidemann, V. Bendo kowsky, B. Butscher, H.P. B¨chler, u and T. Pfau, Phys. Rev. A 80 , 033422 (2009).
- [14] T. Pohl, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 104 , 043002 (2010).
- [15] F. Cinti, P. Jain, M. Boninsegni, A. Micheli, P. Zoller, and G. Pupillo, Phys. Rev. Lett. 105 , 135301 (2010).
- [16] R. M. W. van Bijnen, S. Smit, K. A. H. van Leeuwen, E. J. D. Vredenbregt, and S. J. J. M. F. Kokkelmans, J. Phys. B 44 , 184008 (2011).
- [17] P. Schauß, M. Cheneau, M. Endres, T. Fukuhara, S. Hild, A. Omran, T. Pohl, C. Gross, S. Kuhr, and I. Bloch, Nature 491 , 87 (2012).
- [18] P. Schauß, J. Zeiher, T. Fukuhara, S. Hild, M. Cheneau, T. Macr` ı, T. Pohl, I. Bloch, and C. Gross, arXiv:1404.0980.
- [19] A. Schwarzkopf, D. A. Anderson, N. Thaicharoen, and G. Raithel, Phys. Rev. A 88 , 061406 (2013).
- [20] C. Ates, T. Pohl, T. Pattard, and J.M. Rost, Phys. Rev. Lett. 98 , 023002 (2007).
- [21] T. Amthor, C. Giese, C. S. Hofmann, and M. Weidem¨ uller, Phys. Rev. Lett. 104 , 013001 (2010).
- [22] I. Lesanovsky and J. P. Garrahan, Phys. Rev. A 90 , 011603(R) (2014); M. Marcuzzi, J. Schick, B. Olmos, and I. Lesanovsky, J. Phys. A 47 , 482001 (2014).
- [23] M. G¨ arttner, K. P. Heeg, T. Gasenzer, and J. Evers, Phys. Rev. A 88 , 043410 (2013).
- [24] H. Schempp, G. G¨ unter, M. Robert-de-Saint-Vincent, C. S. Hofmann, D. Breyel, A. Komnik, D. W. Sch¨nleber, o M. G¨ arttner, J. Evers, S. Whitlock, and M. Weidem¨ller, u Phys. Rev. Lett. 112 , 013002 (2014).
- [25] N. Malossi, M.M. Valado, S. Scotto, P. Huillery, P. Pillet, D. Ciampini, E. Arimondo, and O. Morsch, Phys. Rev. Lett. 113 , 023006 (2014).
- [26] J. M. Raimond, G. Vitrant, and S. Haroche, J. Phys. B 14 , L655 (1981).
- [27] T. Baluktsian, B. Huber, R. L¨w, o and T. Pfau, Phys. Rev. Lett. 110 , 123001 (2013).
- [28] C. Carr, R. Ritter, C. G. Wade, C. S. Adams, and K. J. Weatherill, Phys. Rev. Lett. 111 , 113901 (2013).
- [29] R. Heidemann, U. Raitzsch, V. Bendkowsky, B. Butscher, R. L¨ ow, L. Santos, and T. Pfau, Phys. Rev. Lett. 99 , 163601 (2007).
- [30] See Supplemental Material at [URL] for details on the experimental procedure, the adiabatic elimination, the theoretical model, the pair-state potentials, the dephasing mechanism and the size of the aggregates..
- [31] B. Huber, T. Baluktsian, M. Schlagm¨ uller, A. K¨ olle, H. K¨ ubler, R. L¨w, o and T. Pfau, Phys. Rev. Lett. 107 , 243001 (2011).
- [32] K. Singer, J. Stanojevic, M. Weidem¨ller, u and R. Cˆ ot´, e J. Phys. B 38 , S295 (2005).
- [33] A. Schwettmann, J. Crawford, K.R. Overstreet, and J.P. Shaffer, Phys. Rev. A 74 , 020701 (2006).
- [34] G. Vitrant, J. M. Raimond, M. Gross, and S. Haroche, J. Phys. B 15 , L49 (1982).
- [35] The Stark-shift of an ion on the 32S state is C /r 4 4 with C 4 = 2 π ×-213 MHz · µ m , which is an attractive poten4 tial with comparable strength to the Rydberg-Rydberg potentials. Therefore facilitation would similarly occur around an ion, if created.
- [36] A. Urvoy, C. Carr, R. Ritter, C. S. Adams, K. J. Weatherill, and R. L¨w, J. Phys. B o 46 , 245001 (2013).
- [37] P. Siddons, C. S. Adams, C. Ge, and I. G. Hughes, J. Phys. B 41 , 155004 (2008).

## SUPPLEMENTAL MATERIAL

## Experimental details

The 455 nm laser is frequency-stabilized with a blue detuning of δ 7P = 1 5 . GHz to the 6S 1 / 2 , F = 4 → 7P 3 / 2 , F = 5 transition (see Fig. 2). The 1070 nm laser is ′ scanned over the two-photon resonance to the Rydberg state | n S 1 / 2 〉 . Both lasers have an estimated linewidth below 5 MHz. The frequency of the 1070 nm laser is calibrated for each measurement using a Fabry-P´rot ine terferometer and additionally by an EIT-signal [36] to fix the origin of the frequency axis. The 455 nm laser typically has a power of 3 mW. The 1070 nm laser has a power

of 15 W and passes through a Pockels cell, allowing to switch the power of this laser for the experiment as fast as 1 5 ns with a repetition rate of 10 kHz. . The two-photon Rabi frequency is Ω = 2 π × 0 05 . . . . 0 5 GHz, the detuning . to the Rydberg state is ∆ = 2 π ×-1 . . . -10 GHz.

The glass cell is home-made and consists of two 1 mmthick quartz optical flats of 5 × 5 cm 2 separated by a 220 µ m spacer and sealed at the edge. A glass tube was connected to the cell, filled with cesium under vacuum and sealed off. This glass tube serves as a reservoir. The temperature of the cell is kept constant at 200 ◦ C to prevent the metallic cesium from condensing, whereas the temperature of the reservoir is varied between 70 ◦ C and 150 C to tune the atom number density. ◦ The atom number density in the cell is determined for each measurement by performing absorption spectroscopy on the D2-line of cesium [37]. Both beams are linearly polarized, overlapped in a counter-propagating configuration and focused inside the cell to a beam waist of approx. 15 µ m. After passing through the cell the blue beam is focused on a pinhole in order to select the central part of radius 6.25 µ m inside the cell. We verified that the resulting imaged volume is indeed almost cylindrical inside the cell. The change in transmission is detected using a fast amplified silicon photodiode (Femto HSAX-S-1G4-Si). During the measurement, the transmission change of the 455 nm laser is monitored and averaged 300 times by a fast oscilloscope. We ensured that we operate in the linear response regime of the photodiode and measured the conversion efficiency of the detection to be approx. 850 V/W. This allows us to convert the transmission change into the real number of Rydberg atoms that are excited, assuming that no decays and losses of Rydberg atoms occur. This assumption is valid for short times as the time spent by an atom inside the excitation volume is on the order of 20 ns.

## Approximation to 2-level atoms

We reduce our 3-level system to a 2-level system by making use of the adiabatic elimination of the intermediate state. The basic principles of this approximation and the relevant parameters are shown in Fig. S1. The density matrix for the 3-level system is ρ = ( ρ ij ) i,j = g e r { , , } , and we define the density matrix for the effective 2-level system as ˆ = (ˆ ρ ρ ij ) i,j = g r { , } . The populations of states | g 〉 and | r 〉 coincide between both descriptions, i.e. ρ gg = ˆ ρ gg and ρ rr = ˆ ρ rr . The approximation relies on assuming that the intermediate state remains unpopulated ( ∂ ρ t ee = 0 and ρ ee | t =0 = 0). The resulting identity from the optical Bloch equations without decays is

$$0 = \Omega _ { g e } \, \text{Im} ( \rho _ { g e } ) - \Omega _ { e r } \, \text{Im} ( \rho _ { e r } ) \quad \text{(S1)} \quad \text{def}$$

FIG. S1. Left panel: 3-level ladder system with the useful definitions. The transition between the levels | i 〉 and | j 〉 is driven with the Rabi frequency Ω ij . The detuning to the intermediate state is δ . Right panel: Approximation to a 2-level system after adiabatic elimination of the intermediate state | e . 〉 The effective 2-level Rabi frequency is Ω eff = Ω ge Ω er / (2 | δ | ). For simplicity reasons the detuning to state | r 〉 and the arising light shifts are not shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-e2f613d9-fe6a-479a-a2c3-ce2c4778c46a.jpg)

Using ∂ ρ t rr = Ω er Im( ρ er ) one obtains

$$\text{Im} ( \rho _ { \text{ge} } ) = \frac { 1 } { \Omega _ { \text{ge} } } \partial _ { t } \rho _ { \text{rr} } \text{ \quad \ \ } ( S 2 )$$

which means that the absorption on the lower transition is proportional to the excitation rate to the Rydberg state | r . 〉 For the 2-level system we also obtain from the optical Bloch equations ∂ ρ t ˆ rr = Ω eff Im(ˆ ρ gr ) and therefore

$$\text{Im} ( \rho _ { \text{gc} } ) = \frac { \Omega _ { \text{eff} } } { \Omega _ { \text{gc} } } \text{Im} ( \hat { \rho } _ { \text{gr} } ) = \frac { \Omega _ { \text{err} } } { 2 | \delta | } \text{Im} ( \hat { \rho } _ { \text{gr} } ) \quad ( \text{S3} )$$

This last relation shows the link between the 2-level model and what is experimentally measured.

## S-state potentials in cesium

At short inter-atomic distances, i.e. below 1 µ m, the pair-state interaction potentials for Rydberg S-states become very complex, as shown in Fig. S2 for the 32S state. We will focus on this 32S state, but the situation is similar for other principal quantum numbers. Neighboring n ′ Pn ′′ D pair-states interact with the 32S-32S state with weak but resonant dipole-quadrupole interaction, leading to avoided crossings and state mixing. This means that the n ′ Pn ′′ D pair-states, which are dipoleforbidden for laser excitation from the 7P 3 / 2 state in the non-interacting case, carry some admixture ε 32S-32S of the 32S-32S state. Therefore it is possible to excite pair-states at negative detunings, where the detuning is defined with respect to the unperturbed 32S-32S state,

ε

FIG. S2. Top panel: Density plot of the 32S-32S admixture ε 32S-32S versus the inter-atomic distance. The molecular quantum number here is M = 0, and interactions up to the quadrupole-quadrupole order are included [33]. The green (resp. red) line depicts the extrapolated van-der-Waals (resp. dipole-dipole) pair-state potential. Bottom panel: 32S32S admixture ε 32S-32S versus the inter-atomic distance at a -2 GHz potential energy (depicted as the blue dash-dot line in the upper panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-81b475c8-33a4-4dd5-92de-8e68b8d08195.jpg)

which would not be the case with purely repulsive vander-Waals interaction.

These pair state potentials are actually only valid for an excitation close to an atom in the 32S state. Once such a pair is created, it is mainly a n ′ Pn ′′ D pair, with an amplitude of 1 -ε 2 32S-32S &gt; 0 99. . For any subsequent excitation around this pair, one therefore has to consider the pair-state interaction potentials for the 32Sn ′ P and 32Sn ′′ Dpair states. For an S-P pair state the interaction potential is purely of the form C /r 3 3 , since the dipoledipole interaction of an S-P pair with its permutation PS is resonant, and the symmetric (resp. antisymmetric) linear combination exhibits attractive (resp. repulsive) interaction. The corresponding interaction strength is C S-P 3 ≈ -290 MHz · µ m , which is in excellent agreement 3 with the value of C 3 = 2 π ×-258 MHz · µ m extracted 3 from the experimental results. For an S-D pair state the interaction is of van-der-Waals character and repulsive ( C S-D 6 ≈ 11 MHz · µ m ). 6 The n ′′ D component therefore plays no role in our case of aggregation with red detunings.

The admixture ε 32S-32S is the rescaling factor for the Rabi frequency when exciting a n ′ Pn ′′ D pair. Because of the nature of the aggregates, which enclose previously excited Rydberg atoms, at most every second Rydberg excitation is created via this process. The other exci- tations occur either as a direct off-resonant excitation from the ground state, or close to a n ′ Pn ′′ D pair. The Rabi frequency does not need rescaling in both cases. An estimate for the overall rescaling factor is therefore √ ε 32S-32S ∼ 0 3, which is close to the rescaling factors of . 0.38 (for van-der-Waals interaction) and 0.23 (for dipoledipole interaction) needed to match the experimental data to the results from the simulation.

In the whole treatment, we neglected the distance dependency of ε 32S-32S . In order to justify this approximation, let us consider just two pair states | 1 〉 and | 2 〉 with energies E 1 ( r ) = C /r 6 6 and E 2 ( r ) = E 0 -C /r 6 6 , and a small dipole quadrupole interaction W r ( ) = U/r 4 between the two pair states. The resulting admixture of pair state | 1 〉 at small r is ε 1 ≈ W r / E ( ) | 2 ( r ) -E 1 ( r ) | ∝ r 2 if we consider E 0 to be small, and the pair state energy is E ′ 1 ( r ) = ∆ ≈ -C /r 6 6 . Therefore ε 1 ∝ ∆ 1 / 3 under these rather crude approximations. Using the same argument as before we obtain Ω ∝ ∆ 1 / 6 which is a sufficiently weak dependence to be neglected.

## Dephasing rate

The first contribution to dephasing in the experiment is the Doppler effect, characterized by the two-photon Doppler width γ D = | k 455 -k 1070 | √ 8 ln 2 k B T m . Here m is the atomic mass of cesium, T is the temperature and k 455 and k 1070 the wavenumbers of the two lasers, from which we consider the difference because the lasers are counterpropagating. At 200 C we have ◦ γ D = 2 π × 512 MHz.

The other important source of dephasing arises from the short transit time during which the atoms are in the resonant shell (or facilitation region) of width ∆ . r In first approximation, only the velocity component perpendicular to the shell v ⊥ determines the transit time, as depicted in Fig.S3. Therefore we define this motional dephasing rate as

$$\gamma _ { m } = \frac { \langle \mathbf v _ { \perp } \rangle } { \Delta r }$$

where 〈 v ⊥ 〉 = 1 √ 3 √ 8 k B T πm is the one dimensional mean atomic velocity. For a temperature of 200 C, this 1D ◦ mean velocity is 〈 v ⊥ 〉 = 158 m s . -1 . Since both dephasing mechanisms are gaussian, the total dephasing rate is defined as γ = √ ( γ D ) 2 +( γ m 2 ) .

Moreover ∆ r is related to the total dephasing rate and on the interaction potential (see Fig.1). First assuming Rydberg interaction of the van-der-Waals (vdW) type, we rename the relevant quantities as ∆ r vdW , γ vdW m (motional dephasing rate) and γ vdW (total dephasing rate). The width of the resonant shell is also given by [22]

$$\Delta r ^ { \text{vdW} } \approx \frac { 1 } { 3 } \left ( \frac { C _ { 6 } } { \Delta } \right ) ^ { \frac { 1 } { 6 } } \left ( \frac { \gamma ^ { \text{vdW} } } { 2 | \Delta | } \right ) \quad \quad ( \text{S5} )$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0039v2-figure-d040baa2-ebf9-4dc6-9735-ddda8c466e25.jpg)

FIG. S3. Top: Sketched snapshots of the excitation process of a 3-atom aggregate. Red (resp. blue) dots represent Rydberg (resp. ground state) atoms. The grey lines show the resonant shell of width ∆ r for each atom, within which the excitation is facilitated. The velocity component v ⊥ perpendicular to the resonant shell is shown for an atom inside the shell. Bottom: Probabilities P n ( t f ) of creating an n -atom aggregate in an interval t f and typical number of n -atom aggregates η n for N fac = 5000 facilitating atoms, similarly to the experimental situation.

| n           |       1 |       2 |       3 |       4 |       5 |
|-------------|---------|---------|---------|---------|---------|
| P n ( t f ) |    0.91 |   0.084 |  0.0061 | 0.00039 | 2.3e-05 |
| η n         | 4740    | 436     | 32      | 2       | 0.12    |

By combining equations (S4) and (S5) we obtain the following self consistent equation for the motional dephasing rate γ vdW m :

$$\frac { \langle \nu _ { \perp } \rangle } { \gamma _ { m } ^ { \text{vdW} } } = \frac { 1 } { 3 } \left ( \frac { C _ { 6 } } { \Delta } \right ) ^ { \frac { 1 } { 6 } } \left ( \frac { \sqrt { \left ( \gamma _ { D } \right ) ^ { 2 } + \left ( \gamma _ { m } ^ { \text{vdW} } \right ) ^ { 2 } } } { 2 | \Delta | } \right ) \quad ( \text{S} ) \quad \text{ans} \ \frac { \text{and} \ c } { \text{prob} }$$

Solving equation (S6) yields γ vdW m = 2 π × 493 MHz, so that the total dephasing rate is γ vdW = 2 π × 711 MHz.

For the case of pure dipole-dipole (dd superscript) interaction, equation (S5) changes to

$$\Delta r ^ { \text{dd} } \approx \frac { 2 } { 3 } \left ( \frac { C _ { 3 } } { \Delta } \right ) ^ { \frac { 3 } { 3 } } \left ( \frac { \gamma ^ { \text{dd} } } { 2 | \Delta | } \right ) \quad \quad ( \text{S7} ) \quad \text{In 1}$$

and the dephasing rate can be estimated to γ dd = 2 π × 591 MHz with γ dd m = 2 π × 296 MHz. γ vdW and γ dd are the values that were used in the simulations.

## Aggregate size

Since we are not in the frozen gas regime, aggregates are only spatially correlated like the ones shown in the inset of Fig. 1(b) for a short period of time. At longer times the atomic motion will destroy the spatial correlations between the atoms. This will, however, have no impact on the fact that Rydberg excitations are facilitated by already existing one. In the following we will estimate the size of the 'frozen aggregates' [Fig. 1(b)]. The time scale over which the ensemble can be considered frozen is given the transit time of the atoms through the resonant shell of width ∆ r , given by t f = ( γ m ) -1 ≈ 0 32 ns. . During this time interval the ensemble can be considered as a frozen gas for our purposes. We estimate the size of the aggregates that can be formed during this time in a very simplified model. The parameters are the same as in Fig. 2(c) with the extracted values of γ and C 6 : ∆ ≈ 2 π × -2200 MHz, Ω = 2 π × 0 4 . × 100 MHz, N g = 88 µ m -3 , C 6 = 2 π × -109 MHz · µ m 6 and γ vdW = 2 π × 711 MHz. 0 4 . is the rescaling factor for Ω discussed in the main text.

Let us first consider one individual Rydberg atom which acts as the seed for an aggregate. In the resonant shell (sphere of radius r fac = 0 61 . µ m and width ∆ = 51 nm, see Eq. (S4)) around this first atom reside r ν res = N g × 4 πr 2 fac ∆ r = 20 7 atoms in average. . Moreover the size of the resonant shell grows with the number of Rydberg atoms. For this we approximate an n -atom aggregate to a sphere (whose volume is occupied by n Rydberg atoms separated by the facilitation radius). Then the entire resonant shell contains n 2 / 3 × ν res atoms. For each ground state atom in the shell, the time constant for the excitation is τ 0 = γ/ Ω 2 ≈ 71 ns, and therefore the excitation rate of the ( n +1)-th atom is Γ n +1 = n 2 / 3 × ν res /τ 0 . If we simplify the problem and consider that the atoms are excited sequentially, the probabilities P n ( ) that an aggregate formed by t n Rydberg atoms at the time t per pre-existing Rydberg atom obeys the following rate equation:

$$\partial _ { n } P _ { n } & = \Gamma _ { n - 1 } P _ { n - 1 } - \Gamma _ { n } P _ { n } \\ & = \frac { \nu _ { \text{res} } } { \tau _ { 0 } } \left ( ( n - 1 ) ^ { 2 / 3 } P _ { n - 1 } - n ^ { 2 / 3 } P _ { n } \right ) \quad \text{(S8)}$$

In the table in Figure S3 we show the solutions of these equations for aggregates consisting of up to 5 atoms at the time t f .

The number of n -atom aggregates that are excited in the whole excitation volume during t f is given by η n = N fac × P n ( t f ), where N fac is the number of Rydberg atoms which can act as a nucleation grains. As the excitation occur at the boundary of the already excited Rydberg ensemble N fac is only a fraction of the ∼ 50000 Rydberg atoms. The values of η n are shown in Figure S3, with N fac = 5000 (corresponding to the outer shell of a sphere with 50000 Rydberg atoms). For this Rydberg atom number the largest spatially correlated aggregates that we create in our experiment can thus be estimated to consist of approximately 4 atoms.
A. Urvoy F. Ripka I. Lesanovsky D. Booth J. P. Shaffer T. Pfau R. Löw
Strongly correlated growth of Rydberg aggregates in a vapor cell
The observation of strongly interacting many-body phenomena in atomic gases
typically requires ultracold samples. Here we show that the strong interaction
potentials between Rydberg atoms enable the observation of many-body effects in
an atomic vapor, even at room temperature. We excite Rydberg atoms in cesium
vapor and observe in real-time an out-of-equilibrium excitation dynamics that
is consistent with an aggregation mechanism. The experimental observations show
qualitative and quantitative agreement with a microscopic theoretical model.
Numerical simulations reveal that the strongly correlated growth of the
emerging aggregates is reminiscent of soft-matter type systems.