Stabilizing the false vacuum: Mott skyrmions
Márton Kanász-Nagy, Balázs Dóra, Eugene A. Demler, Gergely Zaránd
aa r X i v : . [ c ond - m a t . qu a n t - g a s ] M a r Stabilizing the false vacuum: Mott skyrmions
M. Kan´asz-Nagy , B. D´ora , , E. A. Demler and G. Zar´and BME-MTA Exotic Quantum Phases Research Group,Budapest University of Technology and Economics, Budapest 1521, Hungary Department of Physics, Budapest University of Technology and Economics, Budapest 1521, Hungary and Department of Physics, Harvard University, Cambridge, MA 02138, U.S.A
Topological excitations keep fascinating physicists since many decades. While individual vorticesand solitons emerge and have been observed in many areas of physics, their most intriguing higherdimensional topological relatives, skyrmions (smooth, topologically stable textures) and magneticmonopoles – emerging almost necessarily in any grand unified theory and responsible for chargequantization – remained mostly elusive. Here we propose that loading a three-component nematicsuperfluid such as Na into a deep optical lattice and thereby creating an insulating core, one cancreate topologically stable skyrmion textures and investigate their properties in detail. We showfurthermore that the spectrum of the excitations of the superfluid and their quantum numberschange dramatically in the presence of the skyrmion, and they reflect the presence of a trappedmonopole, as imposed by the skyrmion’s topology.
PACS numbers: 03.65.Vf, 03.75.Lm, 37.10.Jk, 67.85.De, 67.85.Fg
Topological excitations and defects emerge and play akey role in almost any area of physics where spontaneoussymmetry breaking occurs. Domain walls, besides beingthe necessary ingredients of magnetic recording media,also emerge in cosmology, string theory [1], and consti-tute basic quasiparticles in 1+1 dimensional field theo-ries [2]. Line defects such as vortices belong today toour basic understanding of superfluidity and supercon-ductivity [3–5], and dislocations have a central impacton the elastic properties of materials and on melting andfracture [6].Topological point defects and excitations such asmonopoles [7–9] and skyrmions may be even more fasci-nating than vortices and domain walls. Skyrmions, orig-inally proposed to describe hadronic particles [10, 11],emerge as smooth, localized, and topologically protectedtextures of some vector field (order parameter). Theirgauge theoretical counterparts, the ’t Hooft-Polyakovmonopoles appear to be necessary ingredients of almostany grand unified theories, and are needed to explaincharge quantization. While planar or line-like topolog-ical defects are abundant in condensed matter, observ-ing topological point defects in the laboratory proved tobe very hard, since only peculiar order parameters sup-port their existence. Moreover, as shown by Derrick [12],skyrmions are generally doomed to shrink to a point orto expand to infinity and vanish in a homogeneous sys-tem. It is in fact only very recently that the existence ofmonopoles has been convincingly demonstrated in spinice materials such as Dy Ti O [13, 14], and the sponta-neous formation of a skyrmion lattice textures has beenreported in certain magnetic materials [15].The advent of ultracold atoms opened new perspec-tives to creating and manipulating individual skyrmionsand monopoles. It has been noticed very early that cer-tain spin F = 1 multicomponent superfluids such as Na or Rb support magnetic phases [1, 2, 17, 19, 20], which– combined with an inhomogeneous trap geometry [21] –were argued to give rise to stable skyrmion and monopoleconfigurations [17, 22]. In a nematic superfluid suchas Na, in particular, the superfluid order parameterΨ( r ) = (Ψ x , Ψ y , Ψ z ) takes a remarkably simple form,Ψ( r ) = ˆ u ( r ) p ̺ ( r ) e iϕ ( r ) , (1)with ̺ denoting the superfluid density, ϕ the superfluidphase, and ˆ u a real unit vector. From equation (1) we canidentify the topological structure of the nematic order pa-rameter space as (cid:0) S × U(1) (cid:1) / Z , with the unit sphere S corresponding to the orientation of the vector ˆ u andthe U(1) symmetry associated with the phase degree offreedom, ϕ . The curious Z factorization is a consequenceof the fact that phase changes ϕ → ϕ + π are equivalentto flipping the orientation of the vector, ˆ u → − ˆ u . Due tothis structure, topologically nontrivial and thus topolog-ically stable hedgehog-like field configurations ˆ u ( r ) exist(see Fig. 1) and can give rise to skyrmion and monopolestructures in two and three dimensions, respectively.Unfortunately however, topological excitationstrapped in usual cold atomic setups turned out to bemuch more unstable than initially thought; monopolesand skyrmions generically slip out of the trap [22] orthey simply gradually unwind and disappear withina short time [23]. Therefore – so far – only unstableskyrmion configurations have been imprinted andobserved experimentally [23, 24].Here we propose to stabilize the skyrmion states ge-ometrically by generating a non-superfluid core at thecenter of a trapped nematic superfluid. We suggest toachieve this by placing the nematic superfluid into a deepoptical lattice, and thus driving the atoms at the centerinto a Mott insulating state (see Fig. 1). In this way a closed two-dimensional superfluid shell is created, which FIG. 1.
Schematic structure of the Mott skyrmion. ( a ) A Mott insulator (MI) core is surrounded by a superfluid(SF) shell with a skyrmion spin structure. Arrows denote theorientation of the nematic order parameter, ˆ u ( r ) in equation(1). The MI core stabilizes the skyrmion by keeping it fromdrifting out from the trap. ( b ) Schematic finite temperaturephase diagram of strongly interacting bosons in an opticallattice. The red arrow shows the chemical potential as onemoves from the center of the trapped skyrmion towards itsedge. The MI has no magnetic structure at the temperaturesconsidered. – unlike open shell configurations – supports topologi-cally stable skyrmions, anchored by the Mott insulatingcore. We compute the free energy of this strongly inter-acting ’Mott skyrmion’ system numerically, and demon-strate that the skyrmion texture is indeed stable.In the skyrmion configuration the superfluid order pa-rameter acquires a non-trivial, topologically protectedtexture, generated virtually by a monopole at the cen-ter of the trap. As predicted by Jackiw and Rebbiand Hasenfratz and ’t Hooft [25, 26], the presence ofa monopole can influence the quantum numbers of theexcitations around the monopole – and turns bosonic ex-citations to fermions, a phenomenon termed ’spin fromisospin’ mechanism. As we show by detailed calculations,somewhat similarly to the ’spin from isospin’ mechanism,the presence of the monopole changes drastically the ex-citation spectrum of the superfluid (see also [27]), andremoves a pseudospin quantum number, present in the’skyrmionless’ ground state. Given the exceptional sta-bility of the ’Mott skyrmion’, the predicted change in theexcitation spectrum should be experimentally accessible.We describe a balanced mixture of interacting spin F =1 bosons using the lattice Hamiltonian H kin + P r H loc , r with the kinetic and local parts defined as H kin = − J X h r , r ′ i b † r α b r ′ α , (2a) H loc , r = − µ ( r ) n r + U n r : + U ~F r : . (2b)Here the operators b † r α create a boson of spin component α ( α = x, y, z ) at the lattice site r , and n r = P α b † r α b r α and ~F r = P α,β b † r α ~F αβ b r β denote their density and mag-netic moments, respectively. The F j stand for the usualangular momentum matrices in the α = x, y, z basis, F jβγ = − i ǫ jβγ , and : . . . : refers to normal ordering. Thehopping J sets the kinetic energy of the bosons, whilethe effect of trapping potential V ( r ) = mω r / µ ( r ) ≡ µ − V ( r ). The (normal ordered) in-teraction term U describes the strong repulsion betweenlattice-confined bosons, while the second, much weakerinteraction term U accounts for the magnetic interac-tion between them. It is this second term, U >
0, whichfor Na forces the superfluid order parameter Ψ α ∝ h b α i to stay within the nematic phase, f r ≡ Ψ † r ~F Ψ r / | Ψ r | ≡ ̺ r ≡ | Ψ r | becomes finite.In the following we focus our attention onto the regime zJ/U ≈ . z = 6 the number of nearest neighbors.Here, increasing the chemical potential at the center ofthe trap beyond some critical value (or equivalently, mak-ing the trap tighter), the density at the center increasesand finally reaches the first Mott lobe (see Fig. 1b).For Na, in particular, we estimate U ≈
250 nK and U ≈ zJ ≈
50 nK,as shown in Supplementary Note 1. We also assume thatthe temperature is already low enough (
T < zJ ) to form-ing a superfluid around the Mott core with typical radius R , albeit it is still higher than the magnetic ordering tem-perature of the Mott insulating core, T > T C ∼ J /U ,so that the interplay of magnetic ordering in the Mottcore and superfluidity can be ignored. RESULTS
Stable skyrmion configuration
First, to verify the stability of the skyrmion, we in-troduced local order parameter fields by performing aHubbard-Stratonovich transformation, b r → Ψ r ≈ h b ( r ) i (see Methods), and traced out the original boson fieldsnumerically to obtain the free energy functional F ( { Ψ r } ) ≈ − Ja X r , r ′ ,α Ψ r α ∆ rr ′ Ψ r ′ α (3)+ X r F loc (cid:0) ̺ r , f r , µ ( r ) , T (cid:1) . Here a denotes the lattice constant, and ∆ rr ′ stands forthe discrete Laplace operator. The first term describesthe stiffness of the superfluid order parameter, while thesecond term incorporates the effect of the interaction aswell as that of the confining potential, and has been de-termined numerically for each lattice point (see Meth-ods). Its structure follows from the obvious O (3) rota-tional symmetry and U (1) gauge symmetry of the under-lying Hamiltonian, equation (2).To find the minimum of F ( { Ψ r } ) we used the imagi- FIG. 2.
Inner structure of the skyrmion in the ( x , z ) plane. ( a ) In-trap SF densities of the | +1 i ( |− i ) bosonsform a vortex (antivortex) around the equator, whereas thatof the | i condensate in ( b ) creates a dark soliton at the poles.( c , d ) show in-trap particle densities. SF order of one of thespin components leads to a local increase of the component’sparticle density at the expense of those of the other two com-ponents, leading to a specific density structure characterizingthe skyrmion (bottom). This structure gets significantly morepronounced at lower temperatures. [Physical parameters ofthe plot: T /U = 0 . , U /U = 0 . , zJ/U = 0 .
18, chem-ical potential in the middle µ mid /U = 0 .
36, and at the edge µ edge /U = − . nary time equations of motions, − ∂ τ Ψ r α = δFδ Ψ r α , − ∂ τ Ψ r α = δFδ Ψ r α . (4)The dynamics generated by equation (4) drives the fieldconfiguration { Ψ r } towards the minima of the free energyfunctional. In particular, for appropriate parameters,starting from a configuration with a skyrmion textureimprinted, the field Ψ r is found to relax to a configura-tion Ψ r ≈ e iϕ p ̺ ( r ) ˆ r , with ˆ r denoting the radial unit vector. We verified byadding a random component to the initial field configura-tion that this final state is indeed a robust local minimumof the free energy, as anticipated.Fig. 2 displays the superfluid and total densities acrossthe trap in the usual hyperfine spin basis, F z = ± F z = 0, where the amplitudes of the various superfluidcomponents read as Ψ ± ( r ) = e iϕ p ̺ ( r )(ˆ x ± i ˆ y ) / √
2, andΨ ( r ) = e iϕ p ̺ ( r )ˆ z . The superfluid density is clearly suppressed at the center of the trap, where the stabiliz-ing Mott insulating core is formed, and it lives on a twodimensional shell around this core. The components Ψ ± form vortices around the equator, while the Ψ compo-nent behaves as a dark soliton, localized at the north andsouth poles. The total density of the components of thesuperfluid is also distorted and reflects the structure ofthe order parameters; the density of h b i ∼ h b z i is elon-gated along the z -axis, while that of the other two spincomponents, h b ± i ∼ h b x i ± i h b y i is squeezed along it (seeFig. 2 bottom). Creation
A possible way to create the skyrmion is to imprintdiabatically a vortex, an antivortex and a dark solitoninto the three spin components [28], and then stabilizethe vortex state by turning on a deep optical optical lat-tice. Starting with a superfluid with all atoms in the |− i state, as a first step, a fraction of the atoms couldbe transferred into a vortex state in component | i us-ing a so-called Λ transition. This is possible by applyinga diabatic pulse of a pair of counter-propagating σ − , σ + Raman beams of first order Laguerre-Gaussian ( LG − )and Gaussian density profiles, respectively [24]. This vor-tex could then be transferred to state | i by a simple RF π -shift. As a next step, the creation of an antivortex incomponent | i could be achieved by changing the chiral-ity of the Laguerre-Gaussian beam from LG − to LG +1 .Finally, another laser, perpendicular to the quantizationaxis, would be used to imprint the dark soliton into theremaining atoms in state |− i .An alternative and maybe even simpler way to createthe skyrmion could be to imprint three dark solitons inthe x , y , and z directions, respectively, and then mixingthem using an RF π/ Excitation spectrum
To study the excitation spectrum of the superfluidshell, we constructed a two-dimensional effective fieldtheory for the superfluid order parameter, ψ ( r ) by as-suming a thin superfluid shell of radius R . Neglectingthe radial motion of the condensate, we can describe thesuperfluid by the Lagrange density L = i ψ ∂ t ψ + ψ (cid:0) ∆ m + ˜ µ (cid:1) ψ − g | ψ | − g ψ ~F ψ ) , (5)generating the following equations of motion for the orderparameter field, i ∂ t ψ = (cid:0) − ∆ m − ˜ µ + g | ψ | (cid:1) ψ + g ( ψ ~F ψ ) · ~F ψ. (6)Here ∆ denotes the two dimensional Laplace operatoron the sphere, ˜ µ the effective chemical potential of thesuperfluid shell, and and g and g stand for the effec-tive couplings. All of these parameters depend on the FIG. 3.
Excitation spectrum.
The left (right) panel visual-izes the low energy part of the excitation spectrum above thetrivial (skyrmion) ground states in units of ω = 1 / ( mRξ ),with the magnetic healing length ξ = 1 / √ mg ρ . Due tothe rotational symmetry of the trivial state around ˆz , eachexcited state has a 2-fold spin degeneracy in addition to the(2 l + 1)-fold orbital degeneracy. In contrast, spin degeneracysplits, and only orbital (rotational) degeneracies survive inthe skyrmion sector. The trivial and skyrmion states exhibitdifferent number of zero modes (Goldstone modes) as well.Apart from the phase degree of freedom, only two zero modesexist in the trivial case, since rotations around ˆz in configu-ration space leave this state invariant. In the skyrmion state,however, the number of zero modes increases by one, sincerotations around all three spin axes provide a zero mode ontop of the phase mode. precise width of the superfluid shell as well as on thelattice parameters. We estimated them from our latticecomputations, as explained in more detail in Supplemen-tary Note 3.The excitation spectrum of the condensate is obtainedby linearizing equation (6) around the ground state fieldconfiguration, and then solving the resulting coupled dif-ferential equations. Equivalently, we can treat the field ψ as a quantum field, and obtain the corresponding Bogoli-ubov spectrum of the condensate (see Methods and Sup-plementary Note 4). For a uniform, ’skyrmionless’ con-figuration, ψ ∝ √ ̺ ˆz , we find that the density (phase)and spin excitations decouple and the spectrum can beobtained analytically, similarly to the case of a spatiallyhomogeneous systems [29]. In the limit of large trap radiicompared to the superfluid and magnetic healing lengths,defined through ξ − ≡ m g ̺ and ξ − ≡ m g ̺ , we ob-tain the spectrum ω ph , l ≈ mRξ p l ( l + 1) , ω sp ,l ≈ mRξ p l ( l + 1) , with l = 0 , , .. the angular momentum quantum num-ber. Since g ≪ g and thus ξ ≫ ξ , spin excitationsdominate the low energy excitation spectrum of the con-densate. For a spherical trap every excited state in the spin sector has a (2 l + 1) × l + 1)-fold degeneracy is due to spherical symmetryand is accidental in the sense that it is lifted once thetrap is distorted, and spherical symmetry is broken. Theother, two-fold degeneracy is, however, a consequence ofthe residual O (2) symmetry of the vector part of the or-der parameter, and can be interpreted as a spin degen-eracy of the Bogoliubov quasiparticles in the spin sector.Notice that this spin degeneracy is absent in the phasesector, where excitations have only a (2 l + 1)-fold angu-lar momentum degeneracy for a spherical trap. In ad-dition to the finite energy excitations, three zero-energyexcitations (Goldstone modes) are found with quantumnumbers l = 0, corresponding to global phase and spinrotations, respectively (see Fig. 3).Excitations in the skyrmion sector are more compli-cated, since spin and density fluctuations couple to eachother due to the spatial winding of the skyrmion tex-ture. Furthermore, the spherical symmetry is found tobe spontaneously broken, and the condensate slightly ex-tends into a randomly selected direction (see Supplemen-tary Note 4). Here we find four skyrmionic zero-energyGoldstone modes and three more excitations of almostzero energy, associated with the spontaneous symme-try breaking. The rest of the excitations have energies ω ∼ / ( m R ξ ), and they come in groups of (2 l + 1) al-most degenerate excitations, split into l states of degen-eracy 2 and a non-degenerate state, as induced by thespontaneous cylindrical distortion of the superfluid. Theskyrmion’s excitation spectrum is sketched in Fig. 3. No-tice that the structure of the excitation spectrum is com-pletely different from that of the ’skyrmionless’ sector,the spin degeneracy of the spin excitations completelydisappears.The above mentioned coupling between density andspin fluctuations leads to a clear fingerprint of theskyrmion state in lattice modulation experiments. Mod-ulation of the atom tunneling, J , excites oscillations inthe amplitude of the superfluid order parameter, whichare coupled to the low energy spin excitations of theskyrmion. Thus, even though these modulations do notcouple directly to spin degrees of freedom, the topologicalwinding of the skyrmion texture gives rise to an indirectcoupling to spin excitations, leading in effect to ’spin-orbit coupling’. Specifically we consider modulation ofthe tunneling along one axis by varying the intensity ofthe optical lattice. This corresponds to perturbation inthe l = 0 and l = 2 angular momentum channels, and ex-cites the lower l = 2 branch of excitations of the skyrmiontexture in Fig. 3 (see Supplementary Note 5). In contrastto the skyrmion case, spin and density fluctuations arecompletely decoupled in the trivial sector, and modula-tion of the tunneling cannot excite any of the low energyspin modes of this state. Similarly, modulation of thetrapping potential along one direction leads to excita-tions of the low energy spectrum of the skyrmion state,but no excitation of the trivial configuration. We thusfind, that the trivial and skyrmion configurations areclearly distinguishable through analyzing the low energyspectrum of modulation experiments.For typical parameters we find that the energies of thesuperfluid excitation are in the order of 10 Hz. It wouldbe essentially impossible to study these excitations in pre-vious unstable skyrmion configurations, however, the ex-treme stability of our ’Mott skyrmion’, should now allowto access them and to study their dynamics. Detection
The skyrmion texture can be detected in many ways.Although the change in the density of the componentsis moderate, their in trap density difference is ratherpronounced and is clearly detectable through absorptionimaging (see Fig. 4).The skyrmion texture can also be easily detectedthrough time-of-flight (ToF) measurements, imaging themomentum distribution of the trapped atoms [30]. TheToF image of the atoms with spin component α at time t is approximately proportional to [31] n ToF α ∝ C α (cid:16) k = m r t (cid:17) , with C α ( k ) denoting the Fourier transform of the cor-relation function h b † r α b r ′ α i , and is approximately givenby C α ( k ) ≈ (cid:12)(cid:12)P r h Ψ r α i e i kr (cid:12)(cid:12) + const . As we show inFig. 5, the ToF image consists of Bragg peaks —- fin-gerprints of the underlying optical lattice. Each Braggpeak, however, displays a fine structure at a momentum FIG. 4.
Difference of in-trap absorption images of thecomponents |± i and | i , taken along the y axis. Due tothe non-trivial SF configurations of the skyrmion, |± i ( | i )bosons have higher densities along the equator (poles), seealso Fig. 2. [Densities are shown as percentages of the largestvalue of the absorption image of component | i . Physicalparameters: identical to the ones in Fig. 2.] FIG. 5. Time of flight characteristics of the skyrmion. ( a ) Schematic picture of time of flight (ToF) peaks locatedat the reciprocal lattice vectors in momentum space, with anadditional fine structure ( ∼ /R ) reflecting the spatial SF cor-relations of the skyrmion of spatial extent R . ( b ) Structure ofthe doughnut (double peak) shaped ToF peaks of component |± i ( | i ), on the left (right), taken along the z ( y ) axes, top(bottom). [Color code and axes: arbitrary but identical units.Physical parameters: identical to those in Fig. 2.] scale k ∼ /R , characteristic of the skyrmion texture ofthe superfluid. The ToF image of the m = 0 component,e.g., displays a circular structure when imaged from the z direction, while a clear double peak structure shouldbe detected under imaging from the x or y directions. Acknowledgement
Illuminating discussions with ´A. Nagy are gratefullyacknowledged. This research has been supported bythe Hungarian Scientific Research Funds Nos. K101244,K105149, CNK80991 and by the Bolyai Program of theHungarian Academy of Sciences. E. A. D. acknowl-edges support through the DOE (FG02-97ER25308), theHarvard-MIT CUA, the ARO-MURI on Atomtronics,and the ARO MURI Quism program.
METHODS
Free energy.
The partition function is given by Z = R D [ b † , b ] e − S [ b † ,b ] , where the action reads as S [ b † , b ] = Z dτ X r α b † r α ∂ τ b r α + H kin + X r H loc , r , with H kin and H loc , r defined in equation (2). The hop-ping term is decoupled by a Hubbard-Stratonovich trans-formation [32] after introducing the superfluid order pa-rameter Ψ as R D [Ψ † , Ψ] e − P rr ′ Ψ † r Jz I − rr ′ Ψ r ′ e P r Jz (Ψ † r b r + b † r Ψ r ) = e P rr ′ b † r J rr ′ b r ′ , (7)with the hopping term expressed as J rr ′ = Jz I rr ′ , with z = 6 the number of nearest neighbors. Its inversecan be rewritten for small kinetic energies using I − rr ′ ≈ δ rr ′ − a z ∆ rr ′ . Within the saddle point approximation,the functional integral in b can be carried out exactly,leading to the free energy F [Ψ † , Ψ] = − T log Z [Ψ † , Ψ] inequation (3) with F loc (cid:0) ρ r , f r , µ ( r ) , T (cid:1) = Jz̺ r − T log Tr b ( e − ( H loc , r − Jz ( b † r Ψ r +h . c . ) ) /T ) . (8) Numerical minimization of the free energy.
Weuse a modified version of the imaginary time minimiza-tion algorithm implemented in [33]. The imaginary timedynamics of the fields in equation (4) lead to a continuousdecrease in the free energy ∂F∂τ = X r δFδ Ψ r ∂ Ψ r ∂τ + δFδ Ψ † r ∂ Ψ † r ∂τ = − X r (cid:12)(cid:12)(cid:12)(cid:12) δFδ Ψ r (cid:12)(cid:12)(cid:12)(cid:12) < . We separate the kinetic part from the remaining on-sitecontributions F = F kin + P r F loc , r . F kin is calculated us-ing a fast Fourier transform, whereas we use numerical in-terpolation in parameter space to calculate F loc , r at eachsite. For additional details, see Supplementary Note 2. Excitation spectrum.
We analyze excitationsaround the trivial ( ψ t = √ ρ t ˆz ) and skyrmion ( ψ s = √ ρ s ˆr ) configurations, assuming uniform ground statedensities on the sphere. Here ˆz and ˆr denote unit vectorsin the z and radial directions, respectively. We deter-mine the two dimensional superfluid densities by usingthe saddle point approximation for the two dimensionaleffective Lagrangian δ L /δψ = 0, yielding ρ t = µ/g and ρ s = ( µ − /mR ) /g , with R the radius of the sphere.In the latter case, the chemical potential becomes renor-malized by − /mR due to the curvature of the groundstate, leading to the depletion of the superfluid.In the trivial state ψ t the fluctuations parallel ( δψ k )and perpendicular ( δψ ⊥ ) to the ground state decouple, and their equations of motion take on the simple form i∂ t δψ t k = − ∆ m δψ t k + g ρ t ( δψ t k + δψ t k ) , (9a) i∂ t δψ t ⊥ = − ∆ m δψ t ⊥ + g ρ t ( δψ t ⊥ − δψ t ⊥ ) . (9b)The interaction term g sets the velocity of the fluctu-ations of the superfluid phase and density, whereas thespin interaction g determines the velocity of the per-pendicular fluctuations. In the skyrmion state, ψ s , how-ever, the kinetic part in equation (9) acquires geometricalterms due to the curvature of the ground state configu-ration. The Laplace operator is replaced by∆ → D + 2 R , (10)with 1 /mR a curvature term shifting the kinetic energyof fluctuations, and the covariant derivative D = ∇ + i A , defined using the non-Abelian vector-potential A ,mixing the components δψ ⊥ ↔ δψ k . These geometricalterms shift the spectrum and lead to the splitting of theexcitation energies. Further details of this calculation arepresented in Supplementary Note 4. [1] Pospelov, M. et al. Detecting domain walls of axion-like models using terrestrial experiments.
Phys. Rev. Lett.
Solitons and Instantons: An Introductionto Solitons and Instantons in Quantum Field Theory (El-sevier Science B. V., 1982).[3] Tinkham, M.
Introduction to Superconductivity (McGraw-Hill Book Co., New York, 1975).[4] Bloch, I., Dalibard, J. & Zwerger, W. Many-body physicswith ultracold gases.
Rev. Mod. Phys.
The Universe in a Helium Droplet. (Claren-don Press, 2003).[6] Kleman, M. & Friedel, J. Disclinations, dislocations, andcontinuous defects: A reappraisal.
Rev. Mod. Phys.
Proc. Roy. Soc. Lond.
A133,
Nucl. Phys. B
JETP Lett
Nucl. Phys.
Proc. Roy.Soc. Lond.
A260,
J. Math. Phys. et al. Dirac strings and magneticmonopoles in the spin ice Dy Ti O . Science
Nature [15] M¨uhlbauer, S. et al.
Skyrmion lattice in a chiral magnet.
Science
Phys.Rev. Lett.
742 (1998).[17] Al Khawaja, U. & Stoof, H. T. C. Skyrmions in a ferro-magnetic BoseEinstein condensate.
Nature
Phys.Rev. Lett.
Phys. Rev. Lett.
Phys. Rev.Lett.
Phys. Rev. Lett.
Phys.Rev. Lett.
Phys. Rev. Lett.
Phys. Rev. Lett.
Phys. Rev. Lett.
Phys. Rev. Lett.
Phys. Rev. Lett.
Phys. Rep.
Nature
Nat. Phys. Phys. Rev. B
Vortex Dynamics in Bose-Einstein Condensates (Ph.D. Thesis, 2000).
SUPPLEMENTARY INFORMATIONSUPPLEMENTARY NOTE 1
In this section we provide estimates of the parame-ters of the lattice Hamiltonian in equation (2). Theinteraction of spin-1 bosons is determined by their s-wave scattering lengths a and a , with a F denotingthe scattering length in the total hyperfine spin F chan-nel. In case of antiferromagnetic interactions, such asin Na considered in this paper, a > a > U = π a +2 a λ ( V /E R ) / E R , with the V depth of the op-tical lattice, and the recoil energy of the optical latticeof wavelength λ given by E R = h mλ [2]. The hopping,measured in units of U , is given by JU = 3 π / λa + 2 a e − √ V /E R , (S1)and can be easily controlled in an experiment by modify-ing the lattice depth to reach the superfluid-Mott insula-tor transition [2, 3]. The magnitude of the on-site spin in-teractions is related to U by the scattering lengths [2, 4], U U = a − a a + 2 a . (S2)As a result, the interaction term U is suppressed, and isapproximately U ≈ . U in case of Na [3, 4].The parameters used throughout this paper are set us-ing the scattering length of Na, a = 2 .
75 nm, and awavelength λ = 594 nm for the optical lattice as used inthe experiment in Ref. 3. In the zJ/U ≈ . U ≈
250 nK, U ≈ zJ ≈
50 nK, with z = 6 the number ofnearest neighbors. Considering 10 atoms in the trap,the radius of the skyrmion is approximately R = 30 a ≈ µ m, with the lattice constant a = λ/ ≈ . µ m. SUPPLEMENTARY NOTE 2
Here we analyze the structure of the local part of thefree energy, F loc , r , in equations (3, 8), and discuss itsnumerical evaluation. Due to the O (3) symmetry of theHamiltonian in equation (2), F loc , r can be written as afunction of the two rotation-invariant quantities of the F = 1 spin sector: the superfluid density ( ̺ r ), and themagnetic moment (f r = | f r | ). In our numerical simula-tions, in order to evaluate equation (8), we truncate theHilbert space at 5 particles per site, and carry out thetrace numerically. Fig. S6 shows plots of F loc , r in caseof nematic ( U >
0) and ferromagnetic ( U <
0) inter-actions. Nematic condensates favor zero magnetization,f r ≡
0, and thus, the structure of the nematic ground
FIG. S6.
Local part of the free energy. ( a ) and ( b )show F loc ( ̺, f) in case of nematic ( U >
0) and ferromagnetic( U <
0) interactions, respectively. The dots indicate theminima of the free energy, favoring a non-magnetized (fullymagnetized) superfluid in the nematic (ferromagnetic) case.[Physical parameters of the plot:
T /U = 0 . U /U =0 . zJ/U = 0 . µ/U = 0 . state configuration space is ( S × U (1)) / Z , as shownin the main text. Ferromagnetic condensates, on theother hand, are fully magnetized f r ≡
1, and thus, theirground state configuration space is SO (3). This topolog-ical structure, however, holds no topologically protected’Mott skyrmion’ configurations [1]. SUPPLEMENTARY NOTE 3
In order to get an order of magnitude estimate of theexcitation energies of the skyrmion, we need to estimatethe parameters of the two-dimensional effective model inequation (5). Therefore, in this part of the Supplemen- tary Information, we relate these parameters to those ofthe lattice Hamiltonian in equation (2). Note, that theaccuracy of this estimation does not influence the ratio ofthe excitation energies of the skyrmion and of the trivialsector, shown in Fig. 3.We approximate the superfluid shell of the skyrmion bya two-dimensional slab of thickness N ≈
10 lattice sitesin the z direction, and we describe it with the action S = Z dt X r α ib r α ∂ t b r α − H kin + X r H loc , r ! , (S3)where the Hamiltonian has been defined in equation (2).By assuming a constant, time-independent profile in the z -direction for low-energy excitations, we can approxi-mate the action, using the two-dimensional effective La-grangian in equation (5), S ≈ N Z dt Z d r L [ ψ, ψ ] , (S4)where we have introduced the continuum fields ψ α through the substitution b r α /a → ψ α ( r ), using the lat-tice constant a . Thus, the parameters of the Lagrangianare given by m = 1 / (2 Ja ), ˜ µ = µ + 6 Jz , g = U a and g = U a . In the weak coupling limit zJ ≫ U , U this Lagrangian density can be used to describe the ex-citation spectrum in the saddle point approximation, asshown in the main text. In the strongly interacting limitquantum corrections renormalize the parameters of theLagrangian density. We assume, however, that these ef-fects do not change the order of magnitude of excitationenergies significantly, and therefore we use the bare pa-rameters above to estimate the latter.The low energy excitations are of the order E =1 / ( mRξ ) = p g ρ/mR both in the skyrmion and inthe trivial configuration. Assuming a superfluid density ρ = 0 . /a , we estimate E ≈ √ U JR/a ≈ Na system with the lattice parameters given inSupplementary Note 1. Given the increased stability ofthe ’Mott skyrmion’ considered here, these frequenciesshould be in the measurable range.
SUPPLEMENTARY NOTE 4
In what follows, we determine and compare the Bogoli-ubov excitation spectra of the trivial and the skyrmionconfigurations. The excitations were analyzed using theeffective Lagrange density in equation (2), considering athin superfluid shell of radius R . Assuming sphericallysymmetric ground state both in the trivial ( ψ t = √ ρ t ˆz )and in the skyrmion configuration ( ψ s = √ ρ s ˆr ), we de-termined the two-dimensional superfluid densities usingthe saddle point equation δ L /δψ = 0. This equationyields ρ t = µ/g and ρ s = (cid:0) µ − / ( mR ) (cid:1) /g , respec-tively. In the skyrmion case, the chemical potential getsrenormalized due to the curvature of the ground state.This curvature effect leads to the depletion of the super-fluid density and affects the excitation spectrum as well.In the trivial configuration, phase and spin excita-tions associated with the fluctuations parallel ( δψ t k ) andperpendicular ( δψ t ⊥ ) to the ground state ψ t , decou-ple in leading order. The fluctuation part of the La-grangian, expanded up to quadratic order, reads δ L = iδψ ∂ t δψ − H , with the Hamiltonian density defined as H t = δψ t (cid:18) − ∆ m δψ t (cid:19) + g ρ t | δψ t k | + δψ t k + δψ t k + g ρ t | δψ t ⊥ | − δψ t ⊥ + δψ t ⊥ ! . (S6)The Bogoliubov excitation energies can be obtained bytreating the above Hamiltonian quantum mechanically,or, equivalently, by determining the eigenvalues of theequations of motions of the fields i∂ t δψ t k = − ∆ m δψ t k + g ρ t (cid:16) δψ t k + δψ t k (cid:17) ,i∂ t δψ t ⊥ = − ∆ m δψ t ⊥ + g ρ t (cid:0) δψ t ⊥ − δψ t ⊥ (cid:1) . These equations can be easily solved by expanding thefluctuations in terms of spherical harmonics, yielding theeigenfrequencies ω ph ,l = s(cid:18) l ( l + 1)2 mR + g ρ t (cid:19) − ( g ρ t ) ,ω sp ,l = s(cid:18) l ( l + 1)2 mR + g ρ t (cid:19) − ( g ρ t ) , with the angular momentum quantum number takingvalues l = 0 , , . . . . For a spherical trap every excita-tion in the spin sector has a (2 l + 1) × l + 1)-fold degen-erate. Although for a non-spherical trap the (2 l + 1)-fold orbital degeneracy is removed, the 2-fold degeneracyof spin modes, a consequence of spontaneous symmetrybreaking, remains. We find three zero-energy excitations(Goldstone modes) with l = 0 quantum numbers, cor-responding to phase fluctuations and rotations of theground state around the x and y axes. (Rotations aroundthe z axis leave the ground state invariant, therefore, theydo not give additional zero modes.) In the limit of largetrap radii compared to the superfluid ( ξ = 1 / √ mρg )and magnetic healing lengths ( ξ = 1 / √ mρg ) the exci- tation energies become ω ph ,l ≈ mRξ p l ( l + 1) ,ω sp ,l ≈ mRξ p l ( l + 1) . Since g ≪ g , and thus ξ ≫ ξ , the low energy spec-trum is dominated by spin excitations.In the skyrmion case, the topological structure of theground state modifies the excitation spectrum signifi-cantly. The kinetic term of the Hamiltonian in equa-tion (S6) acquires a curvature term ∆ → ∆ + 2 /R dueto the non-trivial spatial structure of the ground state.The equations of motion, therefore, do not decouple andcan only be described by the combined equation i∂ t δψ s = − (cid:18) ∆ m + 1 mR (cid:19) δψ s + g ρ s ( δψ s k + δψ s k )+ g ρ s ( δψ s ⊥ − δψ s ⊥ ) . (S7)Notice that the action of the seemingly harmless Lapla-cian is very non-trivial: it mixes parallel and perpendic-ular fluctuations ( δψ s k ↔ δψ s ⊥ ) due to the skyrmion’sgeometric structure, which can also be described by intro-ducing non-Abelian vector potentials, as shown in equa-tion (10). The excitation energies can be most conve-niently found by expanding the fields in the (orthonor-mal) basis of vector spherical harmonic functions [5] Y lm ( r ) = ˆr Y lm ( r ) , Ψ lm ( r ) = r ∇ Y lm ( r ) / p l ( l + 1) , Φ lm ( r ) = ˆr × Ψ lm ( r ) , defined using the spherical harmonics Y lm of angular mo-mentum quantum numbers l and m . Due to their vecto-rial nature, vector spherical functions form a representa-tion of the total angular momentum operators ~J = ~L + ~F with quantum numbers ( j, m J ) = ( l, m ), where the op-erators ~J account for simultaneous spatial ( ~L ) and spin( ~F ) rotations.As can be seen from the formulas above, the vectorfunctions Y lm , defined for all l ≥
0, always point in theradial direction; therefore, they span the parallel fluc-tuations δψ s k . In particular, the function Y ∝ ψ s corresponds to the skyrmion configuration itself, andthus, the fluctuation of the corresponding expansion co-efficient describes the global phase fluctuations of theskyrmion. Perpendicular fluctuations, on the other hand,are spanned by the fields Ψ lm and Φ lm , which are definedfor l = 1 , , . . . angular momenta.Since the Laplacian leaves the Φ -sector invariant, − ∆ Φ lm = l ( l + 1) R Φ lm , (S8)excitations in this sector decouple from the ( Y , Ψ )-fluctuations, and the corresponding (2 l + 1)-fold degen-0erate excitation energies can be derived analytically, ω Φ ,l = s(cid:18) l ( l + 1) − mR + g ρ s (cid:19) − ( g ρ s ) . (S9)In case of large trap radii, R ≫ ξ , ξ , these become ω Φ ,l ≈ mRξ p l ( l + 1) − . (S10)Specifically, for l = 1 angular momenta, we find threezero energy modes corresponding to the rotations of theskyrmion around the x , y and z axes in parameter space.Therefore, together with the global phase fluctuations inthe Y subspace, there are four Goldstone modes in theskyrmion sector. We thus find an increased number ofGoldstone modes as compared to the trivial sector, dueto the topological winding of the skyrmion.The excitation energies of the ( Y , Ψ )-sector are morecomplicated, since the Laplacian is non-diagonal in thesefields, − ∆ (cid:18) Y lm Ψ lm (cid:19) = 1 R (cid:18) l ( l + 1) + 2 − p l ( l + 1) − p l ( l + 1) l ( l + 1) (cid:19) (cid:18) Y lm Ψ lm (cid:19) , thereby mixing parallel and perpendicular fluctuations.The excitation energies are given by the eigenvalues ofthe Bogoliubov-Hamiltonian matrix H YΨ l = 12 mR (cid:18) Ω l Λ m − Λ m − Ω l (cid:19) , (S11)defined using the matrices Ω l = (cid:18) l ( l + 1) + √ R/ξ − p l ( l + 1) − p l ( l + 1) l ( l + 1) − √ R/ξ (cid:19) and Λ m = ( − m √ R (cid:18) /ξ − /ξ (cid:19) . In a spherically symmetric trap there are two branchesof excitation energies for all l = 1 , , . . . angular mo-menta, both being (2 l + 1) degenerate. In the g ≪ g limit the lower of these branches approaches the energiesof the corresponding ω Φ ,l ∼ / ( mRξ ) spin excitations,whereas the other branch, describing mainly phase exci-tations, stays at large energies ∼ / ( mRξ ).An investigation of the l = 1 excitations reveals a weakinstability of the spherically symmetric ground state ψ s towards a slight uniaxial deformation, as we verifiedthrough detailed numerical simulations. This sponta-neous symmetry breaking does not influence the num-ber of Goldstone modes protected by symmetry, how-ever, as indicated in Fig. 3, it slightly splits the non-zeroenergy excitations due to the breaking of the O (3) ro-tational symmetry of the ground state to O (2). In par-ticular, the lower branch of the l = 1 excitations in the ( Y , Ψ )-sector split in a 3 → (2 + 1)-manner and theirenergies become extremely close to zero. No such insta-bility has been observed in our three-dimensional latticesimulations, though the numerical accuracy of the lattermay have been insufficient to detect this small symmetrybreaking. SUPPLEMENTARY NOTE 5
In this section, we analyze the low energy absorptionspectrum of the skyrmion in a lattice modulation exper-iment, in which atom tunneling along one axis is mod-ulated by periodically varying the depth of the opticallattice. Specifically, modulations along the z axis corre-spond to a variation of the z -hopping parameter in equa-tion (2a). In terms of the two-dimensional effective modelof excitations in equation (5), this corresponds to a ∂ z perturbation operator, as can be seen from the discussionbelow equation (S4). This term has spin F = 0, and it isa linear combination of the tensor operators T , = ( ∂ x + ∂ y + ∂ z ) / √ ,T , = ( ∂ x + ∂ y − ∂ z ) / √ , with angular momentum quantum numbers ( l, m ) =(0 ,
0) and (2 , ψ ∝ Y ∝ ˆ r [6].Since T , and T , are derivative operators, they com-mute with the angular momentum L = − ∆ , and theywill not mix the subspace ( Y lm , Ψ lm ) vector sphericalharmonics with Φ lm functions, the latter forming aneigenspace of L (see equation (S8)). Further selectionrules follow from the rotational symmetries of the pertur-bation operators under spatial and spin rotations, dueto the Wigner-Eckart theorem [7]. Working in the ba-sis of total angular momentum quantum numbers, it canbe easily shown that the only non-vanishing matrix ele-ments describing excitations of the spherically symmetricskyrmion ground state are h Y , | T , | Y , i = − r , h Ψ , | T , | Y , i = 1 √ , h Y , | T , | Y , i = − √ . Therefore, modulations of the atom tunneling along the z axis can only create ( l, m ) = (2 ,
0) excitations, andonly in the ( Y , Ψ ) sector. These correspond to a highenergy density excitation, and a small energy spin ex-citation, the latter being shown in the lower branch ofthe l = 2 levels in Fig. 3. Such low energy levels can-not be excited in the trivial configuration, whose lattice1modulation spectrum contains only high energy densityfluctuations, to linear order. Thus, the presence of sucha low energy excitation peak in the modulation spectrumis an unambiguous fingerprint of the skyrmion texture. [1] Ho T.-L. Spinor Bose Condensates in Optical Traps. Phys.Rev. Lett. , 742-745 (1998).[2] Demler, E. & Zhou, F. Spinor Bosonic Atoms in OpticalLattices: Symmetry Breaking and Fractionalization. Phys.Rev. Lett. , 163001 (2002).[3] Xu, K. et al. , Sodium Bose-Einstein condensates in an optical lattice. Phys. Rev. A , 043604 (2005).[4] Stenger, J., Inouye, S., Stamper-Kurn, D. M., Miesner,H.-J., Chikkatur, A. P. & Ketterle, W. Spin domains inground-state Bose-Einstein condensates. Nature , 345-348 (1998).[5] Hill, E. L. The theory of Vector Spherical Harmonics.
Am.J. Phys. , 211-214 (1954).[6] Spontaneous symmetry breaking amounts to only ∼ Naexperiments. Thus, its effect on the excitation spectrumwill be negligible, and we will assume the skyrmion to bespherically symmetric in the following discussion.[7] Sakurai, J. J.