Supersolid droplet crystal in a dipole-blockaded gas F. Cinti,1 P. Jain,1 M. Boninsegni,1 A. Micheli,2 P. Zoller,2 and G. Pupillo2 1Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2J1 2IQOQI and Institute for Theoretical Physics, University of Innsbruck, 6020 Innsbruck, Austria. (Dated: October 31, 2018) A novel supersolid phase is predicted for an ensemble of Rydberg atoms in the dipole-blockade regime, interacting via a repulsive dipolar potential \softened" at short distances. Using exact numerical techniques, we study the low temperature phase diagram of this system, and observe an intriguing phase consisting of a crystal of mesoscopic superfluid droplets. At low temperature, phase coherence throughout the whole system, and the ensuing bulk superfluidity, are established through tunnelling of identical particles between neighbouring droplets. PACS numbers: 67.80.K-, 32.80.Rm, 67.85.Hj, 67.85.Jk, 67.85.-d, 02.70.Ss The search for novel phases of matter drives much of densates (droplets) organize into a single, global conden- the current research in condensed matter physics. Of sate. Specifically, we consider the following two-body particular interest are phases simultaneously displaying potential: different types of order. A chief example, of great cur- 8 3 rent interest, is the so-called supersolid, namely a phase < D=a if r ≤ a featuring crystalline order, and also capable of sustaining v(r) = ; (1) dissipation-less flow. Attempts to observe experimentally : D=r3 if r > a a supersolid phase of matter, primarily in a crystal of solid helium, have spanned four decades since early the- D being the characteristic strength of the interaction. oretical predictions [1]. The most credible claim of such This kind of interaction potential can be realized with an observation to date [2, 3], has been subjected to in- cold dipole-blockaded Rydberg atoms [6{10]. The pa- depth scrutiny over the past few years, and it seems fair rameters D and a above can be controlled with external to state that agreement is lacking at the present time, as fields [5, 11] (we come back to this point below). to whether experimental findings indeed signal a super- Our system of interest comprises N identical bosons of solid phenomenon [4]. mass m, confined to two dimensions [21]. The many-body A new, fascinating avenue to the observation of super- Hamiltonian is the following (in dimensionless form): solid and other phases of matter not yet observed (or N even thought of), is now opened by advances in cold 1 X X H = − r2 + v(r ) (2) atom physics, providing not only remarkably clean and 2 i ij i=1 i>j controlled experimental systems, but also allowing one to \fashion" artificial inter-particle potentials, not aris- where rij = jri − rjj is the distance between particles i ing in any known condensed matter system. This allows and j, and v is given by Eq. (1). All lengths are expressed one to address a key theoretical question, namely which 2 in terms of the characteristic length r◦ = mD=~ , and we two-body interaction potential(s), if any, can lead to the introduce a dimensionless cutoff Rc = a=r◦ for the poten- occurrence of a supersolid phase in free space (i.e., not tial (1). The system is enclosed in a square cell of area A, on a lattice). with periodic boundary conditions. The particle density is n = N=A, but we shall express our results in terms of In a recent article [5], Henkel et al. have proposed, p 2 the (dimensionless) inter-particle distance rs=1= nr◦. based on a mean-field treatment, that a Bose conden- 3 2 2 The energy scale is ◦ = D=r◦ = ~ =mr◦. sate of particles interacting through an effective poten- The low-temperature phase diagram of such a system tial which flattens off at short distance, might support has been explored by means of first principles numerical a density modulation. In this Letter, we show by first simulations, based on the Continuous-space Worm Al- arXiv:1005.2403v4 [cond-mat.quant-gas] 21 Sep 2010 principle numerical simulations that interaction poten- gorithm [12, 13]. Numerical results shown here pertain tials which combine a long-distance repulsion with a to simulations with a number of particles N varying be- short-distance cutoff, lead in fact to the appearance of a tween 50 and 400, in order to carry out extrapolation novel self-assembled crystalline phase of mesoscopic su- of the results to the thermodynamic limit. Our ground perfluid droplets in a system of bosons. Furthermore, state estimates are obtained as extrapolations of results such a crystal can turn supersolid in the T ! 0 limit, at finite temperature. Details of the simulations are stan- as tunneling of particles across neighbouring droplets dard, as the use of the potential (1) entails no particular takes place, and superfluid phase coherence is established technical difficulty. across the whole system, as individual separate Bose con- In the limit Rc rs, the truncation of the dipolar 2 (a) (b) 1.5 4 (a) 3 1 y g(r) 2 0.5 1 0 0 0.2 0.4 0.6 0.8 0 r (c) (d) 0.2 1.5 (b) s 1 ρ y 0.1 0.5 0 0 10 20 30 40 0 0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 T x x 1.0×10 (c) -2 FIG. 1: Snapshots of a system of bosons interacting via po- 1.0×10 tential (1), at the four different temperatures 200 (a), 20 (b), P -4 1.0 (c) and 0.1 (d), expressed in units of ◦. Points shown are 1.0×10 taken along individual particle world lines. The nominal value of rs in this case is 0.14, whereas the cutoff of the potential -6 1.0×10 (1) is Rc=0.3. 0 20 40 60 80 100 L FIG. 2: (Color online) Results shown are for rs=0.14 and potential at short distances does not play an important Rc=0.3. Temperature is in units of ◦. (a): Pair correla- role, and the low temperature phase diagram of (2) is tion function g(r) at a temperature T = 200 (triangles), 20 that of purely dipolar bosons in two dimensions, investi- (squares), 1.0 (diamonds) and 0.1 (circles). The simulated system comprises N=200 particles. (b): Superfluid density gated previously by several authors [14, 15]. It is known vs. T for systems with N=100 (square), and 200 (diamond) < C that for rs ∼ rs = 0.06 the ground state of the system particles. (c): Frequency of occurrence of permutation cycles > L is a triangular crystal, whereas for rs ∼ rs = 0.08 it is of length L at the same four temperatures reported in panel a uniform superfluid (in the intermediate density range (a). Longer permutation cycles occur at lower temperature. a more complex scenario is predicted [16]). As we show > below, a very different physics sets in when Rc ∼ rs, in the density ranges which correspond to either the crys- talline or superfluid phase in the purely dipolar system. The formation of such droplets is a purely classical Fig. 1 shows typical configurations (i.e., particle world effect, that depends on the flattening off of the repul- lines) produced by Monte Carlo simulations of a system sive inter-particle potential below the cutoff distance. In of bosons interacting via the potential (1), at a nominal fact, a simple estimate of the number Nd of particles per droplet, can be obtained by considering a triangular lat- density corresponding to rs = 0:14, at different temper- atures spanning three orders of magnitude. The value of tice of point-like dipoles, each one of strength / Nd (as it comprises Nd particles), and by minimizing with re- the cutoff Rc in this case is 0.3. At the highest temper- ature, a simple classical gas phase is observed, as shown spect to Nd the potential energy per particle, for a fixed by the pair correlation function g(r), shown in Fig. 2 (a), density. The result is which is just a constant (note that g(r) does not vanish 2 Rc at the origin, owing to the flattening off of the potential Nd = γ (3) r at short distance). As T is decreased, an intriguing ef- s fect takes place, namely particles bunch into mesoscopic where γ ≈ 2:79. Eq. (3) furnishes a fairly accurate esti- droplets, in turn forming a regular (triangular) crystal. mate of Nd for the (wide) range of values of the param- This is shown qualitatively in the snapshots in Fig. 1, but eters rs and Rc explored here. For instance, using the also confirmed quantitatively by the structure of the g(r) parameters of Fig. 1, we find from (3) Nd ≈ 13, which as well (Fig. 2(a)), which displays pronounced, broad agrees quite well with our simulation result. It is worth maxima, as well as well-defined minima, where the func- noting that a similar sort of pattern formation, due to tion approaches zero. We henceforth refer to this phase competing interactions, has been previously established as the droplet-crystal phase. for classical colloidal systems [17, 18]. 3 2 IDC SDC SF (a) 1 L Rc > rs 1 IDC C SF y 0 0.5 C n(k) Rc < rs -1 rs -2 0 FIG. 3: Schematic ground state phase diagram of (2) as a -2 -1 0 1 2 -20 -10 0 10 20 2 function of rs. The superfluid droplet crystal (SDC) is sand- 1 wiched between an insulating droplet crystal (IDC) and a su- (b) < C 1 perfluid (SF).