Strong, long-range interactions between atoms in high-lying Rydberg states can suppress multiple Rydberg excitations within a micron-sized trapping volume and yield sizable Rydberg level shifts at larger distances. Ensembles of atoms in optical microtraps then form Rydberg superatoms with collectively enhanced transition rates to the singly excited state. These superatoms can represent mesoscopic, strongly interacting spins. We study a regular array of such effective spins driven by a laser field tuned to compensate the interaction-induced level shifts between neighboring superatoms. During the initial transient, a few excited superatoms seed a cascade of resonantly facilitated excitation of large clusters of superatoms. Due to spontaneous decay, the system then relaxes to the steady state having nearly universal Rydberg excitation density . This state is characterized by highly non-trivial equilibrium dynamics of quasi-particles—excitation holes in the lattice of Rydberg excited superatoms. We derive an effective many-body model that accounts for hole mobility as well as continuous creation and annihilation of holes upon collisions with each other. We find that holes exhibit a nearly incompressible liquid phase with highly sub-Poissonian number statistics and finite-range density–density correlations.
We study many-body correlations in the stationary state of a continuously driven, strongly interacting dissipative system. Specifically, we examine resonant optical excitations of Rydberg states of atoms interacting via long-range dipole-dipole and van der Waals potentials employing numerical and analytical techniques. Collection of atoms within a blockade distance form a ``superatom'' that can accommodate at most one Rydberg excitation. The superatom excitation probability saturates to
21 for coherently driven atoms but is significantly higher in the presence of dephasing, approaching unity as the number of atoms increases. Using the exact numerical solution of the density-matrix equations for a small system, we demonstrate that strong dephasing of the optically driven dipoles renders the many-body problem amenable to semiclassical Monte Carlo simulations. We employ the Monte Carlo algorithm for a large number of atoms and find that in the steady state of a uniformly driven, extended one-dimensional system, the saturation of superatoms leads to quasicrystallization of Rydberg excitations whose correlations exhibit damped spatial oscillations. We show that the behavior of the system under the van der Waals interaction potential can be approximated by a rate-equations model based on a ``hard-rod'' interatomic potential, and by solving it we obtain the period and correlation length for the density wave of Rydberg excitations.