Defect-induced supersolidity with soft-core bosons F. Cinti,1, 2, ∗ T. Macr`ı,1 W. Lechner,3 G. Pupillo,4 and T. Pohl1 1Max Planck Institute for the Physics of Complex Systems, N¨othnitzer Str. 38, 01187 Dresden, Germany 2National Institute for Theoretical Physics (NITheP), Stellenbosch 7600, South Africa 3IQOQI and Institute for Theoretical Physics, University of Innsbruck, Austria 4IPCMS (UMR 7504) and ISIS (UMR 7006), Universit´ede Strasbourg and CNRS, Strasbourg, France More than 40 years ago, Andreev, Lifshitz, and Chester suggested the possible existence of a peculiar solid phase of matter, the microscopic constituents of which can flow superfluidly without resistance due to the formation of zero-point defects in the ground state of self-assembled crystals. Yet, a physical system where this mechanism is unambiguously established remains to be found, both experimentally and theoretically. Here we investigate the zero-temperature phase diagram of two- dimensional bosons with finite-range soft-core interactions. For low particle densities, the system is show to feature a solid phase in which zero-point vacancies emerge spontaneously and give rice to superfluid flow of particles through the crystal. This provides the first example of defects-induced, continuous-space supersolidity consistent with the Andreev-Lifshitz-Chester scenario. Spontaneous symmetry breaking is a focal principle of condensed matter physics { yet, simultaneous break- ing of fundamentally different symmetries represents a rare phenomenon. A prime example is the so-called su- persolid phase [1], which displays both crystalline and superfluid properties, that is, the simultaneous breaking of continuous translational and global gauge symmetry. The first mentioning of such a state goes back to Gross [2], who predicted the possibility of a density-modulated superfluid phase of weakly interacting bosons described by a classical field. Later, Andreev, Lifshitz [3], and Chester (ALC) [4] conjectured a microscopic mechanism for strongly interacting systems based on two key as- sumptions: first, that the ground state of a bosonic crys- FIG. 1: Crystallization due to singular and soft-core tal contains defects such as vacancies and interstitials, interactions. (a) For singular potentials with pure power- and, second, that these defects can delocalize, thereby, law repulsion, indicated by the blue areas, particles assemble giving rise to superfluidity. However, the physical real- into a commensurate and insulating solid. (b) Upon remov- izability of this scenario has since remained under active ing the singularity, particles can cluster and, under proper debate [5]. conditions, form an incommensurate crystal with more parti- cles than lattice sites. Defect delocalization due to inter-site In 2004, torsional oscillator experiments [6, 7] provided tunneling can then promote a finite superfluid response of the 4 first suggestive evidence for superfluidity in solid He self-assembled crystalline ground state. through a rapid drop of the resonant oscillation period below a critical temperature, viewed indicative for su- perfluid decoupling of a fraction of the He crystal. This finding has sparked a host of new experimental activity stiffening of bulk solid He, and later found no signature [8{18], that, however, challenged the original interpre- of superfluidity on avoiding this effect [25]. As a result, tation and pointed out several artefacts causing a non- there now seems to be consistent experimental and the- supersolid origin of the observations. Theoretical work oretical evidence for the absence of the long-sought su- has established that crystal incommensurability is a nec- persolid phase in He. The mere existence of continuous- essary condition for superfluidity [19] and that zero-point arXiv:1302.4576v2 [cond-mat.quant-gas] 6 Feb 2014 space supersolidity induced by zero-point defects thus re- defects in ground state solid He are prevented by a large mains an open question. activation energy [20, 21]. In addition, Boninsegni et al. In this Article, we show that bosonic particles in- [21] Rota and Boronat [22], Ma et al. [23], and Lech- teracting via soft-core potentials (see Fig.1) provide a ner and Dellago [24] have shown that point-like defects prototype system for addressing this question. Using experience an effective attraction that results in defect- exact numerical techniques, we determine the underlying clustering and phase separation, ruling out the possibility zero-temperature phase diagram, which reveals the emer- of defect-induced supersolidity [21] as in the ALC sce- gence of defect-induced supersolidity in the vicinity of nario. Several experiments [11, 14{16] have shown that commensurate solid phases, as conjectured by ALC [3, 4]. the original observations were caused by shear modulus 2 RESULTS Supersolidity with Soft-core bosons. We consider a two-dimensional ensemble of N bosons with density ρ, interacting via a pair potential of the type V0 V = γ γ : (1) r + Rc γ This interaction approaches a constant value V0=Rc as the inter-particle distance, r, decreases below the soft- core distance Rc, and drops to zero for r > Rc. The limiting case γ ! 1 yields the soft-disc model [26], while γ = 3 and γ = 6 correspond to soft-core dipole-dipole [27] and van der Waals [28] interactions that can be realized with ultracold atoms [28, 29] or polar molecules [30, 31]. Here, we focus on the latter case (γ = 6), for which the Hamiltonian reads FIG. 2: Zero-temperature phase diagram of two- N N X r2 X U dimensional soft-core bosons. The phase diagram dis- H^ = − i + ; (2) 2 1 + r6 plays the emergence of superfluid (SF) and different solid (NS) i=1 i<j ij and supersolid (SS) phases for varying interaction strength U and density ρ. The density on the left y-axis has been 2 2 where the units of length and energy are Rc and ~ =mRc , scaled by the soft-core radius Rc. The rightp axis gives the 2 respectively, and m denotes the particle mass. In these density in units of the inverse area, A = 3(1:6Rc) =2, of units, the zero-temperature physics is controlled by the the unit cell of the high-density solid phase, corresponding to 2 4 the lattice site occupation N=Ns for a given number of parti- dimensionless interaction strength U = mV0=(~ Rc ) and the dimensionless density R2ρ. cles and lattice sites, N and Ns, respectively. For Aρ & 1:5, c the grey region labeled as NS corresponds to a cluster crystal Particles with soft-core interactions have been studied with N=Ns > 1, as indicated by the grey scale. Supersolid previously in the field of soft condensed matter physics phases with different occupation numbers are found between 2 [32{34], in the classical high-temperature regime. One of two hyperbolas, defined by Rc ρU = const: (dotted lines). At the main findings has been that pair potentials with a high densities (Aρ & 3:5) they can be understood in terms negative Fourier component [32] favor the formation of of density-modulated superfluids. In contrast, superfluidity particle clusters, which in turn can crystallize to form a within the low-density supersolid lobes emerges from delocal- so-called cluster-crystal. In the quantum domain, the- ized zero-point defects according to the ALC scenario. The horizontal error bars represent statistical uncertainties and oretical work has so far focused on the regime of weak uncertainties due to the finite stepping of U. interactions and high particle densities [27, 28, 35{38], which was shown to be well described by mean-field cal- culations [39, 40]. In this limit, one finds strongly mod- p ulated superfluid states [2, 26{28] with broken transla- distancesr ¯ = 1= πρ > Rc, the physics is dominated tional symmetry in the form of a density-wave. by the long-range tail of the interaction potential, 6 In the following we investigate strong coupling V ∼ 1=r . For a fixed interaction strength U & 35, domain where correlations and quantum fluctuations we thus find a first order liquid-solid quantum phase are expected to become important. We employ path transition with increasing Aρ, consistent with previous integral Monte Carlo simulations to determine the work on bosons with power-law interactions [30, 41{43]. ground state properties of the Hamiltonian eq.(2) (see In particular, the location of the liquid-solid phase Methods section). The obtained phase diagram, shown transition for very low densities coincides with that in Fig.2, reveals a rich spectrum of phases with varying for pure van der Waals interactions. With increasing interaction strength and density. density, however, the average inter-particle spacing,r ¯, approaches the soft-core radius Rc and drops to values 2 Small particle densities. At small densities Rc ρ . 0:5 for which equation (1) strongly deviates from pure we find two phases: a superfluid and an insulating tri- ∼ 1=r6 interactions and levels off below the turning 1=6 angular crystal composed of singly occupied sites, that point Ro = (5=7) Rc . As a result of the decreasing is, where the number of lattice sites, Ns, equals the repulsive inter-particle forces, the crystal melts again for particle number N. The observed lobe structure of increasing densities. As indicated in Fig.2, we indeed this crystalline region is readily understood by noticing find a re-entrant superfluid at particle densities for that at very low densities, that is, large inter-particle whichr ¯ < Ro. 3 FIG. 3: Characterization of the different phases across the phase diagram. Panels (a)-(f) show the radial pair 2 correlation function g2(r) and the particle density distribution n(r) for ρR0U = 32 and different values of Aρ, indicated in (g). Panel (g) displays the superfluid fraction fs as a function of Aρ, covering all phases indicated in Fig.2. Global superfluidity is not well defined in the phase-separation region (green area) and therefore omitted.