arxiv: v1 [cond-mat.quant-gas] 5 Apr 2017

Transcription

1 Phase transton of ultracold atoms mmersed n a BEC vortex lattce arxv: v1 [cond-mat.uant-gas] 5 Apr 017 R. H. Chavgur, 1 T. Comparn, V.. Bagnato, 1 and M. A. Caracanhas 1 1 Insttuto de Físca de ão Carlos, Unversdade de ão Paulo, C.P. 369, ão Carlos, P, , Brazl Laboratore de Physue tatstue, École Normale upéreure/pl Research Unversty, UPMC, Unversté Pars Dderot, CNR, 4 rue Lhomond, Pars, France (Dated: February 3, 018) We nvestgate the uantum phases of ultracold atoms trapped n a vortex lattce usng a mxture of two bosonc speces (A and B), n the presence of an artfcal gauge feld. Heavy atoms of speces B are confned n the array of vortces generated n speces A, and they are descrbed through a Bose- Hubbard model. In contrast to the optcal-lattce setups, the vortex lattce has an ntrnsc dynamcs, gven by ts Tkachenko modes. Includng these uantum fluctuatons n the effectve model for B atoms yelds an extended Bose-Hubbard model, wth an addtonal phonon -medated long-range attracton. The ground-state phase dagram of ths model s computed through a varatonal ansatz and the uantum Monte Carlo technue. When compared wth the ordnary Bose-Hubbard case, the long-range nteratomc attracton causes a shft and reszng of the Mott-nsulator regons. Fnally, we dscuss the expermental feasblty of the proposed scheme, whch reles on the proper choce of the atomc speces and on a large control of physcal parameters, lke the scatterng lengths and the vortcty. PAC numbers: Kk, De, k I. INTRODUCTION Two-component condensate systems had an enormous mpact n the feld of ultracold atoms, especally snce the expermental realzaton of a Bose-Ensten condensate (BEC) of fermons. Ths conssts of effectve bosonc molecules, the tghtly bound Cooper pars, formed after the fermons be sympathetcally cooled by another bosonc or fermonc atomc speces. The expermental observaton of the BEC-BC crossover wth ths molecular BEC [1] proved once more the relevance of ultracold atoms as a powerful tool to test condensed-matter models, whch can be studed n a hghly controllable envronment [, 3]. The two-speces BEC of bosonc atoms also has rch physcs to be explored. An mportant experment performed n Cornell s group, wth dfferent 87 Rb hyperfne states, addressed the statc propertes of bnary mxtures, ther relatve phase coherence and ther dynamcs [4]. The same group was able to nucleate vortces n ths system [5], and, more recently, to produce the superposton of vortex-lattce BECs [6]. The vortex confguraton n BEC was predcted by Feynman [7], who suggested that a superflud can rotate when perced by an array of uantzed sngulartes or vortces. In 1969, Tkachenko proposed that a vortex lattce n a superflud would support transverse elastcmodes [8]. He showed that a trangular lattce has the lowest energy of all smple lattces wth one vortex per unt cell, and that t s stable for all normal modes. These predctons were expermentally realzed wth dlute BEC gases n 001 [9], and followed by the observaton of vortex-lne oscllatons, that s, the low-energy Tkachenko modes [10]. Once establshed, the vortex-lattce confguraton proved to be stable, ncludng ts normal modes, for a range of rotaton freuences close to, but below that of the external trap confnng potental [11, 1]. In a prevous work [13], we consdered a neutral mpurty mmersed n a vortex lattce, nteractng wth the Tkachenko modes. We addressed the shallow-lattce regme, assumng a uadratc energy dsperson for the mpurty speces, and takng the contnuous lmt for ts momentum. In ths work, an analogous system s studed n the tght-bndng regme, where heavy atoms (mpurtes) are strongly trapped n the stes of the vortex lattce formed by condensed bosons of another speces. Ths vortex-lattce setup should be compared to the technue of trappng ultracold atoms n optcal lattces generated by laser beams. The study of bosonc atoms n optcal lattces led to the breakthrough observaton of the uantum phase transton between the superflud and Mott-nsulatng phases [14], well descrbed by the Bose- Hubbard (BH) model for lattce bosons [15]. Important results obtaned wth optcal lattces rely on the fact that these are rgd (.e., do not support phonons [16]) and free of defects. To better smulate the condensed-matter models for sold-state crystals, however, a recent trend has been to ntroduce artfcal dynamcs n these lght crystals [17 19]. In our proposed scheme, the dynamcs emerges naturally through the normal-modes exctatons of the vortex lattce. The effectve Hamltonan for the lattce-confned mpurtes corresponds to a polaronc Bose-Hubbard (BH) model, wth parameters that are modfed by the lattce dynamcs, notably by the ntroducton of a long-range phonon - medated nteracton. The generalzaton of BH models n the ultracold atoms context s connected to the use of atoms wth strong dpole moment, for whch the strength of nteratomc nteracton decays slowly wth the dstance [0]. In ths case, the system s descrbed through an extended Bose- Hubbard (EBH) model, wth the addtonal long-range nteracton generatng a rch varety of phases. For the

2 one-dmensonal case, the phase dagram ncludes the pecular Haldane-nsulator phase [1]. In two dmensons, the predcted phases range from densty-wave to supersoldty [, 3]. Ref. [4], n partcular, reports the expermental realzaton of a EBH model n three dmensons. The vortex-lattce setup proposed here also realzes a EBH models, due to the effectve long-range attracton generated by the lattce dynamcs. mlarly to optcal lattces, the hgh and ndependent control of the system parameters would allow one to explore the uantum phases of mpurtes n a vortex lattce. The relevant expermental methods nclude the use of Feshbach resonances to tune atomc scatterng lengths [5] and of the selectve absorpton-mage technue to characterze the uantum-state confguraton [3, 6, 7]. Ths paper s structured as follows: ec. II presents the dervaton of an effectve Bose-Hubbard model for mpurty atoms trapped by the vortex lattce of the other atomc speces. The phase dagram of ths model s determned through the uantum Monte Carlo technue, and compared to the exstng results (cf. ec. III). In ec. IV, we gve a beyond mean-feld treatment for speces A, whch leads to the EBH model for mpurty atoms, ncludng long-range nteractons. The phase dagram of ths dynamcal model s analyzed n ec. V. In ec. VI, we descrbe the physcal parameters relevant for an expermental realzaton of our proposal, and the man conclusons are reported n VII. II. PHYICAL MODEL We consder a two-component mxture of speces A and B, n a uas-d geometry. The vortex lattce n speces A s excted by the artfcal-gauge-feld technue [8], where the nternal atomc structure s carefully engneered by the optcal potentals to produce Berry phases and nucleate vortces [9]. The man dea s to replace the rotatng trap mechansm by usng an artfcal vector potental A, whch selectvely couples to A atoms and does not affect speces B. The Hamltonan H of the system s the sum of the followng three terms: [ H A = d A(r)) r ˆψ A (r)( ˆψA (r) + m A + ˆψ A (r)v ext(r) ˆψ A (r) + g ( A ˆψ A (r) ˆψ ) ] A (r), [ H B = d r ˆψ B (r)( ) ˆψB (r) + m B + ˆψ B (r)v ext(r) ˆψ B (r) + g ( B ˆψ B (r) ˆψ ) ] B (r), H AB = g AB d r ˆψ A (r) ˆψ B (r) ˆψ A (r) ˆψ B (r), (1) where speces {A, B} s descrbed by the creaton (destructon) feld-operator ˆψ (r) ( ˆψ (r)) at the two-dmensonal poston r = (x, y). The strength of ntra and nter-speces repulsve contact nteractons are gven by g = π a /m l z and g AB = π a AB /m AB l AB z, respectvely, where m AB = m A m B /(m A + m B ) s the reduced mass, and a and a AB are the ntra and nter-speces s-wave scatterng lengths. The length scales l z = /(m ω z ) and l AB z = /(mab ω z ) are functons of the harmonc confnng freuency ω z n the transverse drecton. The vector potental A determnes the vortcty Ω for the BEC A along the z axs through Ω = A /m A. At a crtcal vortcty, the resdual confnng potental for A tends to zero [30]. For speces B, the effect of the slowly-varyng external trap potental V ext (r) may be ncluded through the Local Densty Approxmaton. Wthn a mean-feld treatment of H A, the vortexlattce wave functon for speces A s bult as the lnear combnaton of degenerate lowest-landau-levels solutons of the rotatonal Gross-Ptaevsk euaton. Ths yelds an Abrkosov lattce of vortces, encoded n the wave functon ψ A (r) = n A ϕ A (r), where n A = N A / s the average atomc densty (wth beng the surface area). The r-dependent factor reads [30, 31] ϕ A (r) = (υ) 1/4 ϑ 1 (ζ πυ, ρ) exp ζ ζ, () where ζ represents the complex varable (x + ıy)/lz A, ρ = exp(ıπτ), τ = u + ıυ, u = 1/ and υ = 3/. The trangular lattce of vortces s formed by the zeros of the Jacob theta functon ϑ 1. For the stablty of ths vortex lattce we assume a vanshng temperature, T = 0. Then, n spte of not havng a phase-coherent system, as part of the atoms A are outsde the condensate state, we stll have well establshed vortex-lattce densty profle [3]. Moreover, we consder the mean-feld uantum-hall regme [31], wth the number of vortces N V well below the number of atoms n A, that s, a hgh rato ν = N A /N V 1. In ths regme, nether uantum nor thermal fluctuatons affect the vortexlattce stablty. We also assume N B N V N A, that allow us to dsregard the effects of the dlute speces on the stablty of the vortex lattce. In the above scenaro, we can apply the Bogolubov transformaton n the feld operator ˆψ A to derve the exctatons around the vortex-lattce fundamental state, that s, to nclude the uantum fluctuatons beyond the meanfeld approxmaton for A: ˆψ A = ψ A + δ ˆψ A. We consder the grand-canoncal formalsm, wth chemcal potentals µ for the speces A and B, and rewrte the total system Hamltonan n E. (1) as an expanson n powers of δ ˆψ A. To the second-order, ths expanson reads K = H B + H (0) AB µ ˆN B B }{{} K B (1) + HAB }{{} H nt () + H A µ ˆN A A }{{} KA BOG (3) where ˆN A and ˆN B are the number operators for speces A and B, and where the superscrpts n H () A, H(0) AB and

3 3 H (1) AB ndcate the order of the expanson. Consderng the valdty of the mean-feld vortex-lattce soluton for speces A, the coeffcent of the frst order term H (1) A s zero. We truncate the nteracton term H AB to frst order, and we do not show here the mean-feld contrbuton to the energy of speces A. In ths secton we only consder the 0th order term n E. (3), K B, n whch speces A appears as an effectve mean-feld potental for speces B. The nter-speces nteracton term reads H (0) AB = d r V A (r) ˆψ B (r) ˆψ B (r), where V A (r) = n A g AB ϕ A (r) wth ϕ A (r) gven by E. (), consttutes the statc lattce potental seen by B atoms. In the dlute regme for speces B (.e., for N B N A ), the repulsve nterspeces nteracton causes the localzaton of B atoms at the vortex-core postons, whch s energetcally favorable as the densty of A atoms vansh there. Ths effect s at the the bass of the dervaton of an effectve Bose-Hubbard model for B atoms, as detaled below. We focus on the tght-bndng regme, where the lattce depth V 0 = n A g AB s much larger than the recol energy E r = /(m B ξ ). E r s the natural energy-scale for B atoms trapped n a vortex core, wth radus approxmately eual to the healng length ξ(= / m A n A g A ) of the speces A BEC. The parameter Γ LLL = n A g A / Ω s assocated to the lowest-landau-level constrant (Γ LLL < 1, cf. Ref. [11]), and t connects wth the vortex-lattce densty n V through n V 1/πd, where d = /(m A Ω) s the nter-vortex separaton [7, 1]. We assume that the energy-level spacng between the Bloch bands s large compared to the relevant energes of processes nvolvng B atoms, whch are then restrcted to the lowest-energy band. The sngle-band assumpton s especally vald n the tght-bndng regme (V 0 /E r 1) consdered n ths work. Ths allows us to expand the feld operator for speces B n terms of Bloch wave functons: ˆψ B (r) = k Φ k (r) ˆb k, (4) where ˆb k destroys a partcle n a uas-momentum state k. The Wanner functon ϕ B (r j ), defned through Φ k (r) = (1/ N V ) j ϕ B(r j ) e ık Rj, s localzed at the vortex stes R j (r j = r R j ) and s normalzed to one. Expandng the feld operator as a sum of Wanner functons n each lattce ste, and consderng only nearestneghbor hoppng and on-ste nteracton, K B can be rewrtten as a BH Hamltonan [15]: K B = J,j ˆb ˆb j + U ˆn (ˆn 1) µ B ˆn, (5) where ˆb = (1/ N V ) k eık Rˆb k, ˆn = ˆb ˆb, and, j denotes nearest-neghbor pars. The hoppng coeffcent and on-ste repulson strength are gven respectvely by J = d r ϕ B(r ) [ ] + g AB n A ϕ A (r) ϕ B (r j ), m B (6) and U = g B d r ϕ B (r) 4. (7) U can be determned by assumng a Gaussan functon for ϕ B (r) (= B 0 e r /l 0 ) [33], wth d r ϕ B (r) = 1 and wdth l 0 = ξ/ m B V 0 (harmonc approach for the vortex-core densty profle [34]). Ths gves U = g B /πl 0. On the other hand, the Gaussan ansatz leads to a poor approxmaton for J, whch s better determned through the soluton of a 1D Matheu euaton [35]. In analogy wth the case of optcal lattces, the resultng expresson for J reads J = 4 ( ) 3 ( ) V0 4 V0 E r exp. (8) π E r Usng the vortex-lattce potental depth and recol energy n the expressons for J and U, we obtan U J = 1 ( ) 1/4 a B aa m aab m AB B a l B e A m AB. (9) z a AB m B In contrast wth optcal-lattce setups, where U/J s controlled by the laser-beam propertes and atomc scatterng length, n the vortex lattce ths rato s connected to the atomc propertes through the nter- and ntra-speces scatterng lengths (a A, a B, and a AB ), whch can be ndependently controlled wth unform magnetc felds va the Feshbach-resonance technue [5]. The mappng n E. (9) opens the possblty of usng the vortex-lattce setup descrbed n ths paper to explore the BH phase dagram on a trangular lattce, as wll be shown n ec. III. III. E r BOE-HUBBARD PHAE DIAGRAM At zero temperature, the Bose-Hubbard model n E. (5) features a phase transton between a Mott nsulator (MI), wth an nteger number of atoms per ste, and a phase-coherent superflud (F) phase [36]. The phase dagram has been characterzed n detal, for several geometres and dmensonaltes (see Ref. [37] for a revew). For the two-dmensonal trangular lattce, one can compute the phase boundares at several levels of approxmaton. The smplest case conssts n a mean-feld approach [15]. By decouplng the knetc part of E. (5) through ˆb ˆb j ˆb ˆb j + ˆb ˆb j ˆb ˆb j, the Hamltonan becomes a sum of sngle-ste terms: K B = K() B. Assumng a real and homogeneous local order parameter Ψ = ˆb, each term reads K () B = zjψ(ˆb +ˆb )+zjψ + U ˆn (ˆn 1) µ B ˆn, (10) where z s the number of neghbors per ste (z = 6, for a trangular lattce). The frst term n E. (10) s then treated through second-order perturbaton theory

4 4 [36, 38]. We expand the resultng energy spectrum n powers of the order parameter Ψ, and apply Landau crteron to dentfy the phase transton [38]. The boundary between the superflud phase and the Mott nsulator wth fllng factor g s gven by U(g 1) zj U zuj(g + 1) + z µ B = ± J, (11) where the + ( ) sgn refers to the upper (lower) boundary of the Mott-nsulator lobe (cf. Fg. 1). µb/j F Mean feld Techmann et al., 010 QMC (ths work) MI U/J Fgure 1. Phase boundary between the Mott-nsulator lobe wth fllng factor g = 1 (MI) and the superflud phase (F), on a two-dmensonal trangular lattce. Dfferent boundares are obtaned through mean-feld theory (blue sold lne, see E. (11)), process-chan approach (dashed red lne, from Ref. [39]), and the uantum Monte Carlo technue (green crcles, cf. Appendx B). The mean-feld phase dagram s known to underestmate the area of the Mott-nsulator lobes. More accurate boundares are obtaned through the dagrammatc process-chan approach, where a perturbaton seres n J s computed up to hgh order. On the suare lattce, ths technue yelds a phase dagram n extremely good agreement wth unbased uantum Monte Carlo (QMC) results [40]. Results of the process-chan approach are also avalable on the trangular lattce [39], for whch the tp of the Mott-nsulator lobe wth fllng g = 1 s at (U/J, µ B /J) crt (6.6, 10.) cf. Fg. 1. We also study the trangular-lattce BH model through the worm-algorthm QMC technue [41, 4]. We extract the crtcal ponts through the fnte-sze-scalng analyss of the superflud densty (cf. Appendx B), and the numercal results are n good agreement wth the process-chan phase boundary (see Fg. 1). In ec. V, the same QMC algorthm s used to characterze the extended Bose-Hubbard model. IV. VORTEX-LATTICE DYNAMIC The model derved and characterzed n ectons II and III does not nclude a pecular aspect of the vortex-lattce physcs, namely the ntrnsc dynamcs determned by ts normal modes. The ncluson of these modes modfes the hoppng ampltude and the nteractons for B atoms. In addton to the trappng mechansm whch keeps B atoms n the vortex lattce of A, we also consder the scatterng of the B atoms by the Tkachenko modes of speces A. In the followng, we derve the effectve BH Hamltonan that follows from the ncluson of the lattce vbratons,.e., the uantum fluctuatons beyond the mean-feld vortex-lattce soluton n E. (). Accordng to E. (3), the frst-order term n powers of δ ˆψ A yelds [13] H nt = g AB d r [δ ˆψ A ˆψ B ˆψ B ψ A + ψa ˆψ B ˆψ B δ ˆψ ] A. (1) The second-order term,.e., the Hamltonan KA BOG for speces A, s dagonalzed through the Bogolubov-mode expanson δ ˆψ A (r) = 1 [ u (r)â v (r)â ]. (13) Ths expanson leads to KA BOG = ɛ â â, plus a constant term that only shfts the mean-feld chemcal potental. Ths expresson s obtaned for specfc values of u and v, and t ncludes the operator â (â ) whch creates (annhlates) a Tkachenko-mode exctaton wth wave number and energy dsperson ɛ [31]. In analogy wth the dervaton of E. (5), we expand ˆψ B n E. (1) n terms of localzed Wanner functons, whch yelds K B + H nt = J,j where Λ j = Λ j = + g AB na ˆb ˆb j + U,j ˆn (ˆn 1) µ B ˆn + [ Λ j j â + Λ â ] ˆb ˆb j, (14) d r [ ϕ A(r)u (r) ϕ A (r)v (r) ] ϕ B(r )ϕ B (r j ), d r [ ϕ A (r)u (r) ϕ A(r)v (r) ] ϕ B(r )ϕ B (r j ). (15) Λ j s suppressed for large separatons between stes and j, and t can be approxmately set to zero for R R j larger than d (nearest-neghbor approach) cf. Appendx A. Here, however, we only consder the case of = j, that s, we also neglect the nearest-neghbor pars (, j) n the sum n the last term of E. (14) (thus neglectng nduced-tunnelng terms). The next step s to apply a untary transformaton n the total Hamltonan gven by E. (3), to cancel the

5 5 nteracton term between mpurtes and lattce modes. Ths wll be ncorporated n the coeffcents of the mpurty BH Hamltonan, provdng an effectve polaronc model [43, 44]. We consder H = e U He U = H + [U, H] + 1 [U, [U, H]] +..., (16)! wth U = 1,j 1 ɛ [ α,j â α,j â ] ˆnj. (17) The transformed Hamltonan depends on how the mpurty and lattce-modes operators are modfed. For the real-space mpurty operator, we have e Uˆb e U = ˆb ˆX, (18) wth ˆX = eŷ and Ŷ = (1/ ) ( α, â α, ) â, whle the momentum-space lattce-mode operator transforms as e U â e U = â 1 e ı R α, ˆn. (19) By replacng the felds n E. (3) wth the transformed ones from Es. (18) and (19), and choosng α, = (g AB na /ɛ )Λ to exactly cancel the mpurty-lattcemodes nteracton, we obtan K eff B = J ˆb ˆb j + Ũ ˆn (ˆn 1) µ B ˆn,j V,j ˆn ˆn j. (0),j Above, we neglected retardaton effects, assumng that the lattce exctatons nstantaneously follow the moton of a heavy B mpurty. Also, to obtan E. (0), we traced out the lattce degrees of freedom from the transformed K eff B BOG Hamltonan = K A + K B + H nt ph, wth ph = N and N beng a number state of the lattce modes [45]. The effectve BH parameters for the trapped speces B n the presence of the lattce modes read J = J ˆX ˆX j ph = [ = J exp g AB n A Ũ = U n Ag AB µ B = µ B + n Ag AB ] Λ (1 e ı (Rj R) ), ɛ (1) Λ ɛ, () Λ ɛ. (3) Besdes changes of the Bose-Hubbard-model parameters (J, U, and µ B ), the lattce dynamcs also nduces an attractve long-range potental between the mpurtes, whch s medated by the lattce modes. Ths corresponds to the last term n E. (0), where V,j = n Ag AB ( Λ ) Λ jj ɛ. (4) Based on the expressons for the Bogolubov coeffcents u and v at small, we can estmate the strength V,j of the sotropc long-range attracton between atoms on neghborng stes and j, at a lattce dstance d. By usng the uadratc Tkachenko-mode dsperson ɛ [13], we fnd (see Appendx A for a detaled dervaton) V,j e d /l 0 d. (5) The strong suppresson of V,j at large dstance justfes consderng a truncated model, n whch V,j = V when and j are nearest neghbors, and V,j = 0 otherwse. The parameter V > 0 uantfes the strength of the nearestneghbor attracton, and the departure from the ordnary BH model. V. EXTENDED BOE-HUBBARD PHAE DIAGRAM In ths secton, we characterze the phase dagram of the resultng EBH model descrbed n ec. IV by E. (0). We frst treat the case of vanshng hoppng parameter J (atomc lmt), and we assume zv < Ũ to guarantee the system stablty. For larger values of the nearestneghbor attracton V, t s energetcally favorable to add an nfnte number of partcles n the system. The smplest ground state ansatz s a homogeneous state wth g partcles per ste: Ψ g = 1 ) g (ˆb g! 0, (6) where runs over all lattce stes. By mnmzng Ψ g KB eff Ψ g wth respect to g, we obtan the groundstate fllng factor. The boundary between the regons wth g and g + 1 atoms per ste s gven by µ (g,g+1) B = gũ g + 1 zv, (7) and the resultng phase dagram s represented n Fg.. The Mott-nsulator regons of the ordnary BH phase dagram (V = 0) reman present for all values of V, but ther sze and poston are modfed. The sze of each Mott regon decreases as V ncreases, wth the extenson of the g-th regon along the µ B axs beng eual to Ũ zv. Moreover, the phase boundares have a negatve slope as a functon of V. Ths shft follows from the fact that the nearest-neghbor attracton acts as an addtonal chemcal potental, favorng the addton of more atoms n the lattce.

6 6 µb/ũ 1 0 g = 3 g > 8 g = g = 1 g = zv/ũ Fgure. Phase dagram of the EBH model (cf. E. (0)) n the atomc lmt ( J = 0), formed by Mott-nsulator regons wth nteger fllng factor g. Phase boundares are obtaned through E. (7) (regons wth fllng g > 8 are merged together, for clarty). The system s unstable for zv/ũ > 1. The ansatz Ψ g cannot reproduce nhomogeneous states, lke the charge-densty-wave dentfed for repulsve nearest-neghbor nteractons. More general ansätze may be used, that nclude nhomogeneous state. For nstance one can consder the state Ψ g1,g,g 3, wth nteger fllng factors g 1, g, and g 3 on the three sub-lattces represented n Fg. 3. The numercal mnmzaton of the varatonal energy Ψ g1,g,g 3 KB eff Ψ g 1,g,g 3 leads to g 1 = g = g 3, mplyng that the ground state falls n the class of homogeneous states cf. Es. (6) and (7). Note that, for repulsve nearest-neghbor nteractons (not treated n ths work), the trpartte ansatz Ψ g1,g,g 3 would also produce ground states wth fractonal fllng factors (e.g., wth g 1 = 1 and g = g 3 = 0). Fgure 3. Partton of a trangular lattce nto three sublattces, represented wth dfferent symbols. For J = 0, the EBH Hamltonan s dagonal n the bass of states wth fxed fllng factors on all stes, so that fndng the ground state for a L L fnte lattce corresponds to the optmzaton of a functon of L varables (that s, the local occupaton numbers). We address ths multdmensonal optmzaton problem through the smulated-annealng algorthm [46] for extended systems (L 1). The numercal results confrm that the ground state s homogeneous, wth fllng factor determned by E. (7). The two varatonal ansätze and the drectoptmzaton method employed n the atomc lmt cannot be drectly generalzed to the J > 0 case, where the Hamltonan also ncludes off-dagonal terms. To determne the phase dagram n ths regon (see Fg. 4), we compute the average fllng factor (.e., the densty ρ) and superflud fracton ρ s /ρ through QMC smulatons of large systems at low temperature, effectvely probng the ground state. At small J/ Ũ, we observe a drect transton between MI states wth dfferent fllng factors (cf. Fg. 4(a)), wth the phase boundares gven approxmately by E. (7). For larger values of the hoppng coeffcent, the transton towards the superflud phase s sgnaled by the superflud fracton acurng a fnte value. The poston of ths transton, for V > 0, s well captured by a shfted verson of the V = 0 phase boundary. If we denote the crtcal hoppng by J/Ũ = f(g, µ B/Ũ, V/Ũ), then Fg. 4(b) suggests that the smple relaton ( f g, µ B Ũ, Ṽ ) ( = f g, µ ) B gzv, 0 (8) U Ũ holds close to the tps of MI lobes, that s, for largeenough J/Ũ. The rght-hand sde n E. (8) can be obtaned from the process-chan results [39]. VI. EXPERIMENTAL CONIDERATION We now dscuss the expermental feasblty of our proposal, startng from the crtera on the choce of the mxture. The frst reured condton s that m B m A, snce the lattce potental felt by the mpurtes (speces B) must be deep enough that they can be trapped by the vortex cores. Ths condton s also fulflled wth g AB suffcently large. For speces A, we consder bosonc atoms for whch the formaton of hghly ordered vortex lattces has been observed expermentally [1], namely the alkals Na, L, and Rb. We lmt our dscusson to the commonly explored mxture of 3 Na and 87 Rb, to realze the vortex lattce and heavy mpurtes, respectvely. The BEC lfetme s lmted by three-body losses. Reurng that t exceeds a few tens of seconds, we choose n A (10 0 /m 3 ) l z, wth an effectve BEC thckness l z (ω z π 5 khz). Assumng that the bare (non-resonant) scatterng length nteracton s a A = 60a 0, where a 0 s the Bohr radus, the chemcal potental reads µ A = 1.4 khz, ( 70 nk). To form a vortex lattce that s suffcently large, homogeneous and stable, we consder smlar values for the rotaton and the nteracton energes n A g A / Ω 1. At ths crtcal array vortcty Ω, the resdual radal trap vanshes n E. (1) (the harmonc oscllator length of V ext (r) concdes wth the magnetc length l = /(m A Ω)) and one has tghtly packed trangular vortex-lattce geometry, wth lattce parameter d = l 1.6 µm and stes of sze ξ 0.4 µ.

7 7 µb/ũ Fllng factor (a) J/Ũ µb/ũ uperflud fracton (b) J/Ũ Fgure 4. Average fllng factor ρ (a) and superflud fracton ρ s/ρ (b) for the Hamltonan n E. (0), wth nearest-neghbor attracton V/Ũ = Color code: QMC data for a lattce of stes wth perodc boundary condtons, at temperature T = Ũ/0 (smulaton are performed at ponts marked by gray dots). Dashed whte lne: Phase boundares of the Mott nsulators wth g = 1,, 3, obtaned by shftng the V = 0 data from Ref. [39] cf. E. (8). For speces B, the recol energy E r of an mpurty localzed n the vortex core has to be much smaller than the potental barrer of magntude V 0. Takng nto account all consderatons elaborated above, we fnd a reasonable value V 0 /E r 1 for a AB = 40a 0. It s mportant to stress that wth our chosen parameters, the vortex-lattce lfetme can be of the order of several seconds, whle the characterstc tme assocated wth the tunnellng J n E. (8) s approxmately 0. s. In addton, we assume that the occupaton n the stes are low enough that the states of the mpurtes are accurately descrbed usng the lowest-band Wanner functons. Remndng that the densty of stes n V s constrant by n A n the hgh fllngfactor regme, ν = N A /N V 1, the mpurtes have a neglgble nfluence on the lattce bosons, due to the dluteness of speces B (n B n V ). We can then estmate the characterstc energes for the Bose-Hubbard coeffcents of speces B n E. (5). For the hoppng parameter we have J µ A. In the two-speces vortex-lattce setup, the rato U/J can be tuned by changng the scatterng length a B, as shown by E. (9). To access the regon of the phase dagram that corresponds to the MI-F transton, assumng the untary occupaton of the lattce stes, we consder a crt B 77a 0, whch gves U/J 6.6. The extended Bose-Hubbard (EBH) n E. (0), ncludes an attractve long-range potental wth magntude V µ A (cf. E. (A1) n Appendx A). Importantly, we have V/Ũ 0.05, whch obeys the stablty condton for the ground-state soluton of the EBH, zv/ũ < 1. Usng a 7 L- 133 Cs mxture would reure a magnetc feld on the order of 850 G to render the scatterng length of 7 L postve and suffcently large. An advantage of ths mxture would be the hgh value of the rato V 0 /E r whch could be obtaned wthout sgnfcantly ncreasng a AB. Workng wth a more massve mpurty than 87 Rb would allow one to realze the MI-F transton usng smaller values for the scatterng lengths. The dsadvantage, however, would be the hgher three-body loss rate estmated for 7 L [47], reducng lfetme for the vortex lattce. VII. CONCLUION In the present work, we propose a setup to realze a BH model wth a mxture of ultracold atomc gases n the presence of an effectve rotaton (namely, an artfcal gauge-feld). mlarly to the typcal optcal-lattce setups, the tunablty of the physcal parameters of the system allows to explore a wde range of regmes for the effectve BH model. The vortex-lattce Tkachenko modes, n partcular, modfes the BH parameters and ntroduces an addtonal long-range attracton. The EBH model wth long-range repulson (stemmng from strong dpolar nteratomc nteractons) has been studed n detal, and t led to the predcton of addtonal phases whch are not present n the ordnary BH case, lke the densty-wave and supersold phases n two dmensons [48, 49]. For the attractve EBH model n E. (0), we compute the phase dagram through the uantum Monte Carlo technue. The Mott-nsulator regons of the BH case reman n the extended model, but ther sze and poston are modfed by the long-range attracton. A recent publcaton [50] treats a smlar topc, that s, the phase dagram of mpurtes n a vortex lattce. However, Ref. [50] concerns the vortex-lattce deformatons caused by the strong nteracton wth the mpurty, and how the hgher occupaton of the stes affects the BH parameters. In ths work, n contrast, we studed

8 8 the weak-couplng lmt, focussng on the nteracton of the mpurty wth the vortex-lattce degrees of freedom. After establshng the stablty condtons for a vortex lattce n the presence of multple mpurtes, we explored the effects of the lattce dynamcs on the confned speces. Treatng the dlute system by means of an effectve polaronc Hamltonan, we showed how t allows to go beyond the present studes wth atoms trapped n statc optcal lattces. Besdes beng part of an unusual BH class, our proposed model s also an nterestng expermental proposal n the context of Bose-Ensten condensate mxtures, snce t reures the applcaton of the most recent and successful technues n the ultracold atoms feld: The advances n coolng mechansm to produce bnary condensates [51], the Feshbach-resonance technue to control the nteracton parameters, and artfcal gauge-felds to selectvely nucleate vortex lattces n one of the atomc speces [9]. The possble uantum phases can be characterzed through the spatal nose correlatons n the absorpton mage of the free expanded atomc cloud [5]. mlar technues can be used to observe sgnatures of more exotc states, such as spn luds [53]. Fermons could also be consdered n ths same framework. Usng two hyperfne states, nterestng effects are expected to arse n a Mott-nsulator phase of the pseudospn fermons, wth the nterplay between the trangularlattce geometry and spn orderng frustraton [16, 54]. In the optcal lattce, however, the study of uantum spn models and of strongly-correlated magnetc phases, as the spn-lud phase, has been lmted by the hgh temperature of the fermons n the lattce, whch s stll hgher than that reured to observe exchange-drven spn orderng. In our case, n contrast, the BEC vortex-lattce background can act lke a reservor, wth the coolng of the trapped fermons comng from the creaton of exctatons n ths reservor. Ths favorable scenaro was proposed n the context of dsspatve Hubbard models [55], where the strong control over many parameters, as nteractons between atoms, allows one to manpulate the system-reservor couplng. ACKNOWLEDGMENT We thank Marco D Lberto for nsghtful comments, and Martn Holthaus for provdng the data from Ref. [39]. Ths work s supported by CNP and FAPEP. [1] W. Zwerger, n The BC-BEC Crossover and the Untary Ferm Gas (Lecture Notes n Physcs) 836, prnger (011). [] I. Bloch, J. Dalbard, and. Nascmbène, Nat. Phys. 8, 67 (01). [3] O. Mandel, M. Grener, A. Wdera, T. Rom, T. W. Hänsch, and I. Bloch, Nature 45, 937 (003). [4] C. J. Myatt, E. A. Burt, R. W. Ghrst, E. A. Cornell, and C. E. Weman, Phys. Rev. Lett. 78, 586 (1997). [5] M. R. Matthews, B. P. Anderson, P. C. Haljan, D.. Hall, C. E. Weman, and E. A. Cornell, Phys. Rev. Lett. 83, 498 (1999). [6] V. chwekhard, I. Coddngton, P. Engels,. Tung, and E. A. Cornell, Phys. Rev. Lett. 93, (004). [7] R. P. Feynman, Prog. Low Temp. Phys. 1, 17 (1955). [8] V. K. Tkachenko, J. Exp. Theor. Phys. 56, 1763 (1969). [9] J. R. Abo-haeer, C. Raman, J. M. Vogels, and W. Ketterle, cence 9, 476 (001). [10] I. Coddngton, P. Engels, V. chwekhard, and E. A. Cornell, Phys. Rev. Lett. 91, (003). [11] V. chwekhard, I. Coddngton, P. Engels, V. P. Mogendorff, and E. A. Cornell, Phys. Rev. Lett. 9, (004). [1] A. Fetter, Rev. Mod. Phys. 81, 647 (009). [13] M. A. Caracanhas, V.. Bagnato, and R. G. Perera, Phys. Rev. Lett. 111, (013). [14] M. Grener, M. O. Mandel, T. Esslnger, T. Hänsch, and I. Bloch, Nature 415, 39 (00). [15] D. Jaksch, C. Bruder, J. I. Crac, C. W. Gardner, and P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). [16] M. Lewensten, A. anpera, and V. Ahufnger, Ultracold Atoms n Optcal Lattces: mulatng uantum manybody systems, Oxford Unversty Press (01). [17] A. Gressner, A. J. Daley,. R. Clark, D. Jaksch, and P. Zoller, New J. Phys. 9, 44 (007). [18] A. Klen, M. Bruderer,. R. Clark, and D. Jaksch, New J. Phys. 9, 411 (007). [19] M. Bruderer, A. Klen,. R. Clark, and D. Jaksch, Phys. Rev. A 76, (R) (007). [0] K. Góral, L. antos, and M. Lewensten, Phys. Rev. Lett. 88, (00). [1] E. G. Dalla Torre, E. Berg, and E. Altman, Phys. Rev. Lett 97, (006). [] T. Lahaye, C. Menott, L. antos, M. Lewensten, and T. Pfau, Rep. Prog. Phys. 7, (009). [3] C. Trefzger, C. Menott, B. Capogrosso-ansone, and M. Lewensten, J. Phys. B: At. Mol. Opt. Phys (011). [4]. Baer, M. J. Mark, D. Petter, K. Akawa, L. Chomaz, Z. Ca, M. Baranov, P. Zoller, and F. Ferlano, cence 35, 01 (016). [5] C. Chn, R. Grmm, P. Julenne, and E. Tesnga, Rev. Mod. Phys 8, 15 (010). [6] R. Landg, L. Hruby, N. Dogra, M. Landn, R. Mottl, T. Donner, and T. Esslnger, Nature 53, 476 (016). [7] J. T. tewart, J. P. Gaebler, and D.. Jn, Nature 454, 744 (008). [8] J. Dalbard, F. Gerber, G. Juzelūnas, and P. Öhberg, Rev. Mod. Phys. 83, 153 (011). [9] I. B. pelman, Phys. Rev. A 79, (009). [30]. I. Matveenko, D. Kovrzhn,. Ouvry, and G. V. hlyapnkov, Phys. Rev. A 80, (009). [31]. I. Matveenko, and G. V. hlyapnkov, Phys. Rev. A 83, (011).

9 9 [3] J. nova, C. B. Hanna, and A. H. MacDonald, Phys. Rev. Lett. 89, (00). [33] I. Bloch, J. Dalbard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (008). [34] A. Fetter, n Lectures n Theoretcal Physcs vol. XIB, p. 351, Gordon and Breach, New York (1969). [35] W. Zwerger, J. Opt. B: Quantum emclasscal Opt. 5, 9 (003). [36] M. P. A. Fsher, P. B. Wechman, G. Grnsten, and D.. Fsher, Phys. Rev. B 40, 546 (1989). [37] K. V. Kruttsky, Phys. Rep. 607, 1 (016). [38] D. van Oosten, P. van der traten, and H. T. C. toof, Phys. Rev. A 63, (001). [39] N. Techmann, D. Hnrchs, and M. Holthaus, EPL 91, (010). [40] N. Techmann, D. Hnrchs, M. Holthaus, and A. Eckardt, Phys. Rev. B 79, 4515 (009). [41] A. F. Albuuerue et al. (ALP collaboraton), J. of Magn. and Magn. Materals 310, 1187 (007). [4] B. Bauer et al., J. tat. Mech.: Theory Exp., P05001 (011). [43] D. Benjamn, and E. Demler, Phys. Rev. A 89, (014). [44] K. Agarwal, I. Martn, M. Lukn, and E. Demler, Phys. Rev. B 87, (013). [45] ee secton 4.31 and 4.3 of G. D. Mahan, Many-Partcle Physcs Plenum, New York, 3rd (000). [46]. Krkpatrck, C. D. Gelatt Jr. and M. P. Vecch, cence, 0, 671 (1983). [47] N. Gross, Z. hotan,. Kokkelmans, and L. Khaykovch, Phys. Rev. Lett. 103, 1630 (009). [48] C. Trefzger, C. Menott, and M. Lewensten, Phys. Rev. A 78, (008). [49] J.-Y. Gan, Y.-C. Wen, and Y. Yu, Phys. Rev. B 75, (007). [50] T. H. Johnson, Y. Yuan, W. Bao,. R. Clark, C. Foot, and D. Jaksch, Phys. Rev. Lett. 116, 4040 (016). [51] C. C. Tannoudj, and D. G. Odeln, n Advances n Atomc Physcs: An Overvew, World centfc (011). [5] E. Altman, E. Demler, and M. D. Lukn, Phys. Rev. A 70, (004). [53] R. Bach, and K. Rz ażewsk, Phys. Rev. Lett. 9, (004). [54] D. Gref, T. Uehlnger, G. Jotzu, L. Tarruell, and T. Esslnger, cence 340, 1307 (013). [55] A. Gressner, A. J. Daley,. R. Clark, D. Jaksch, and P. Zoller, New J. Phys (007). [56] N. V. Prokof ev, B. V. vstunov, and I.. Tuptsyn, Phys. Lett. A 38, 53 (1998). [57] J. Šmakov, and E. & ørensen, Phys. Rev. Lett. 95, (005). Appendx A: Medated potental Based on the Bogolubov transformaton of the vortex-lattce [31], we can determne the profle of the long-range effectve potental gven by E. (4) n the man text. In partcular, ts asymptotc behavor can be estmated analytcally wth the low-energy (long-wavelength) Tkachenko modes contrbuton, as wll be shown. As studed before [13], for l 1, wth the magnetc length l related to the nter-vortex separaton through d = l, the gapless Tkachenko modes have dsperson relaton ɛ /M, wth M = 1 κ Ω η n A g A m A and the lattce constant κ = and η = In the low-energy lmt, we have u (r) ϕ A (r) c 1 e ı r and v (r) ϕ A (r) c e ı r. The small value of the momentum allows us to expand (c 1 c ) 1 η 1/4 ( l). Accordng to the Extended Bose-Hubbard Hamltonan (E. (0)), we can determne the profle of the medated potental from E. (15), assumng agan the gaussan functon for ϕ B (r) (= B 0 e r /l 0 ) V,j = V (d) n Ag AB [ 1 η 1/4 ( l) = g AB g A 1 κ B 0 4 e d = g AB g A 4π κ 1 d e d ɛ /l 0 π/l 0 ] [ [ d d r B 0 e r d /l 0 e r /l 0 rdre r /l 0 ] π ] dθe rd cos θ/l 0 0 /l 0. (A1) Ths result justfes the restrcton of the potental range to nearest-neghborng stes. Appendx B: Detals of the uantum Monte Carlo calculatons The uantum Monte Carlo smulatons n ths work make use of the worm algorthm [56], as mplemented n the ALP lbrares [41, 4]. For the bosonc models n Es. (5) and (0), ths represents an unbased method to compute observables as the densty or energy, at fnte temperature T and for a fnte number of stes L. The phase boundary showed n Fg. 1 s obtaned through the fnte-sze scalng technue [36], by fndng the crtcal hoppng parameter J crt for fxed U and µ B. Away from the tp of the Mott-nsulator lobe, a generc phase transton

10 s expected, wth dynamcal and correlaton-length crtcal exponents eual to and 1/, respectvely [36, 37]. Under ths assumpton, and for J close to J crt, the superflud stffness ρ s satsfes ( ) L Jcrt J ρ s = F L 1, J crt L, (B1) T where F s a unversal functon. Thus J crt can be obtaned by plottng L ρ s as a functon of J for several lnear szes L, at fxed T L, and by extractng the common ntersecton pont of these lnes (see for nstance Ref. [57]). When rescaled as n E. (B1), superflud-stffness lnes correspondng to dfferent values of L collapse onto a sngle curve, confrmng the correct estmate of the crtcal value J crt see the example n Fg L z ρs L = 0 L = 16 L = L 1/ν [J/U (J/U) crt ] Fgure 5. Data collapse of the superflud stffness ρ s: Data for dfferent lnear szes L (see legend) are rescaled as n E. (B1), wth z = and ν = 1/ beng the dynamcal and correlaton-length crtcal exponents. Wth the product T L z kept eual to 0.5, dfferent lnes collapse onto a sngle unversal curve. Data are shown for chemcal potental µ B/U = 0.645, wth (J/U) crt =