Two-dimensional supersolidity in a dipolar quantum gas Matthew A. Norcia,1, ∗ Claudia Politi,1, 2, ∗ Lauritz Klaus,1, 2 Elena Poli,2 Maximilian Sohmen,1, 2 Manfred J. Mark,1, 2 Russell Bisset,2 Luis Santos,3 and Francesca Ferlaino1, 2, y 1Institut f¨urQuantenoptik und Quanteninformation, Osterreichische¨ Akademie der Wissenschaften, Innsbruck, Austria 2Institut f¨urExperimentalphysik, Universit¨atInnsbruck, Austria 3Institut f¨urTheoretische Physik, Leibniz Universit¨atHannover, Germany (Dated: February 11, 2021) Supersolidity | a quantum-mechanical phenomenon characterized by the presence of both su- perfluidity and crystalline order | was initially envisioned in the context of bulk solid helium, as a possible answer to the question of whether a solid could have superfluid properties [1{5]. While supersolidity has not been observed in solid helium (despite much effort)[6], ultracold atomic gases have provided a fundamentally new approach, recently enabling the observation and study of super- solids with dipolar atoms [7{16]. However, unlike the proposed phenomena in helium, these gaseous systems have so far only shown supersolidity along a single direction. By crossing a structural phase transition similar to those occurring in ionic chains [17{20], quantum wires [21, 22], and theoreti- cally in chains of individual dipolar particles [23, 24], we demonstrate the extension of supersolid properties into two dimensions, providing an important step closer to the bulk situation envisioned in helium. This opens the possibility of studying rich excitation properties [25{28], including vortex formation [29{31], as well as ground-state phases with varied geometrical structure [7, 32] in a highly flexible and controllable system. Ultracold atoms have recently offered a fundamentally stantially higher atom number N and relatively strong re- new direction for the creation of supersolids | rather pulsive contact interactions between atoms. This leads to than looking for superfluid properties in a solid system the formation of large numbers of loosely bound droplets, like 4He, ultracold atoms allow one to induce a crys- enabling us to establish phase coherence in two dimen- talline structure in a gaseous superfluid, a system which sions. In our system, the repulsive dipolar interactions provides far greater opportunity for control and obser- between droplets facilitate a structural transition from vation. This new perspective has enabled supersolid a linear to a two-dimensional array, analogous to the properties to be observed in systems with spin-orbit cou- Coulomb-interaction-mediated structural phase transi- pling [33] or long-range cavity-mediated interactions [34], tions observed with ions [17{20]. Unlike ions however, though in these cases the crystalline structure is exter- our droplets are compressible and result from the spon- nally imposed, yielding an incompressible state. In con- taneous formation of a density wave, allowing for dynam- trast, dipolar quantum gases of highly magnetic atoms ical variation in both droplet number and size. Further, can spontaneously form crystalline structure due to in- the exchange of particles between droplets enables the trinsic interactions [11{13], allowing for a supersolid with spontaneous synchronization of the internal phase of each both crystalline and superfluid excitations [14{16]. In droplet across the system, and the associated superfluid these demonstrations, supersolid properties have only excitations [14{16]. arXiv:2102.05555v1 [cond-mat.quant-gas] 10 Feb 2021 been observed along a single dimension, as a linear chain Dipolar quantum gases exhibit a rich set of ground- of phase-coherent \droplets", i.e. regions of high density and excited-state phenomena due to the competition connected by low-density bridges of condensed atoms, between many energetic contributions. These include confined within an elongated optical trap. mean-field interactions of both contact and dipolar na- The extension of supersolidity into two dimensions is a ture, quantum fluctuations, and external confinement, key step towards creating an ultracold gas supersolid that parameterized by potentially anisotropic trapping fre- is closer to the states envisioned in solid helium. Com- quencies fx;y;z. Such systems can be described with pared to previous studies of incoherent two-dimensional great accuracy by using an extended Gross{Pitaevskii dipolar droplet crystals [8, 35], we work with both a sub- equation (eGPE) [36{39]. Even a fine variation of the 2 a. 10 droplets, immersed in a low-density coherent \halo" that establishes phase-coherence across the system. As the trap becomes more round, the initially linear chain of droplets acquires greater transverse structure, eventually forming a zig-zag state consisting of two offset linear ar- 0 rays. x (µm) Although the eGPE has remarkable predictive power, full simulations in three dimensions are numerically intensive, making a global survey of the array properties -10 as a function of our experimental parameters difficult. -5 0 5 -5 0 5 -5 0 5 To overcome this limitation, we employ a variational y (µm) y (µm) y (µm) ansatz that captures the key behavior of the system, and allows us to disentangle the competing energetic contri- b. 104 butions. In this approach, we describe an array of ND 7 2D9 PND droplets by the wavefunction (r) = j=1 j(r), where 2D8 the j{th droplet is assumed to be of the form: j(r) / rρ,j rz;j p 1 jρ−ρj j 1 jz−zj j 2D Nj exp − exp − , in- 6 1D6 7 2 σρ,j 2 σz;j N terpolating between a Gaussian and a flat-top profile characteristic of quantum droplets [40]. For a given total 5 number of atoms N and droplet number ND, energy 2D6 minimization provides the atom number Nj in each droplet, as well as their widths σρ(z);j, exponents rρ(z);j, 1D5 Atom number 4 and positions ρ = (x ; y ). Repeating this energy 2D5 j j j minimization as a function of ND gives the optimal number of droplets. This model provides a good quali- 3 1D tative description of the overall phase diagram (Fig.1b), 4 1D 3 revealing that the interplay between intra-droplet 0.34 0.36 0.38 0.4 physics and inter-droplet interaction results in a rich Trap aspect ratio t landscape of structural transitions as a function of the atom number and the trap aspect ratio αt = fx=fy. FIG. 1. Calculated phases of dipolar droplet array. a. Several trends are immediately visible from the phase In-trap ground-state density profiles calculated using eGPE diagram. Larger N and higher αt generally produce for atom numbers N 2 [3:3; 4:4; 5:8] × 104 in the droplets states with larger numbers of droplets. Further, as with ions, a large number of droplets favors a 2D configuration, and trap aspect ratios αt = fx=fy 2 [0:33; 0:35; 0:39] (left while tighter transverse confinement (small αt) favors 1D to right). The scattering length a = 88 a0, where a0 is the [17{20]. A transition from 1D to 2D is thus expected Bohr radius. Green dots depict the droplet positions obtained when moving towards larger N or to higher αt. In stark contrast to the case of ions, the number of droplets typi- from the variational model, assuming the same N and droplet cally increases across the 1D to 2D transition, implying a number ND as the eGPE. Stars connect to experimentally first-order nature, while only narrow regions in the phase observed density profiles in Fig.2b. b. Phase diagram, ob- diagram may allow for a 1D-to-2D transition at constant droplet number. tained from our variational model, as a function of N and αt The variational results are in excellent agreement with for fx = 33 Hz, fz = 167 Hz. Linear (two-dimensional) phases our eGPE numerics, in terms of predicting the qualitative with ND droplets are labelled as 1DND (2DND ). structure of droplet array patterns, as shown in Fig.1a. Slight discrepancies exist between the two theories re- garding the predicted droplet positions and the location strength of these energetic contributions can lead to dra- of the 1D-to-2D transition. This is likely because of the matic qualitative changes in the state of the system, for presence of the halo in the eGPE simulation (and pre- example enabling a transition from a uniform conden- sumably in the experiment), visible in Fig.1a, which is sate to a supersolid, or in our present case, from a linear not accounted for in the variational model. This halo ap- supersolid to a two-dimensional one. pears to accumulate at the ends of the trap, pushing the Fig.1a shows ground-state density profiles calculated droplets toward the trap center and likely increasing the across this transition using the eGPE at zero temper- effective trap aspect ratio experienced by the droplets. ature. These profiles feature arrays of high-density To explore the 1D to 2D transition experimentally, we 3 a. c. d. z 0.44 B y x 0.42 0.4 0.38 atomic aspect ratio 0.36 0.34 b. 0.32 2 0.3 4 6 8 10 0.28 Number of droplets lower transverse confinement FIG. 2. Linear to zig-zag transition in an anisotropic trap. a. We confine and condense dipolar 164Dy atoms within an anisotropic optical dipole trap (ODT) formed by the intersection of two laser beams. By tuning the aspect ratio of the trap in the x-y plane (αt), perpendicular to an applied magnetic field B, we induce a transition between linear and zig-zag configurations of droplets. b. Single-trial images of the in-trap density profile of atoms at different αt, showing structural transition from linear to zig-zag states, as well as an increase in droplet number for higher αt. Stars indicate values αt and N corresponding to the eGPE calculations of Fig.1a.