arxiv: v1 [cond-mat.quant-gas] 14 Oct 2014

Transcription

1 September, Random-Phase-Approxmaton Exctaton Spectra for Bose-Hubbard Models Jamshd Morad Kurdestany Jawaharlal Nehru Centre For Advanced Scentfc Research, Jakkur, Bangalore, Inda and Department of Physcs and Astronomy, Unversty of Mssour, Columba, MO 5-7, USA arv:.v [cond-mat.uant-gas] Oct Ramesh V. Pa Department of Physcs, Goa Unversty, alegao Plateau, Goa, Inda Rahul Pandt Centre for Condensed Matter heory, Department of Physcs, Indan Insttute of Scence, Bangalore 5, Inda. (Dated: September, ) We obtan the exctaton spectra of the followng three generalzed Bose-Hubbard (BH) models: () a two-speces generalzaton of the spnless BH model, () a sngle-speces, spn- BH model, and () the extended Bose-Hubbard model (EBH) for spnless nteractng bosons of one speces. In all the phases of these models we provde a unfed treatment of random-phase-approxmaton (RPA) exctaton spectra. hese spectra have gaps n all the MI phases and gaps n the DW phases n the EBH model; they are gapless n all the SF phases n these models and n the SS phases n the EBH model. We obtan the dependence of gaps and (b) the sound velocty u s on the parameters of these models and examne and u s as these systems go through phase transtons. At the SF-MI transtons n the spn- BH model, u s goes to zero contnuously (dscontnuously) for MI phases wth an odd (even) number of bosons per ste; ths s consstent wth the natures of these transton n mean-feld theory. In the SF phases of these models, our exctaton spectra agree ualtatvely, at weak couplngs, wth those that can be obtaned from Gross-Ptaevsk-type models. We compare the results of our work wth earler studes of related models and dscuss mplcatons for experments. PACS numbers: 7.5.De, 7.5.Fg,.75.Kk, 5..Jp I. INRODUCION Cold-atom systems, such as spn-polarzed 7 Rb n traps, have provded us excellent laboratores for studes of uantum phase transtons [, ]. One example of such a transton s the superflud (SF) to bosonc Mott-nsulator (MI) transton [, ]. hs SF-MI transton, frst predcted by mean-feld studes [5, ] and later nvestgated n Monte-Carlo smulatons [7] of the Bose-Hubbard (BH) model, has been obtaned n experments n systems of nteractng bosons n optcal lattces [,, ]. Expermental realzatons have also been found for generalzed BH models. For nstance, recent experments [ ] on a degenerate mxture of two types of bosons, namely, 7 Rb and K, n a three-dmensonal optcal lattce, have yelded a laboratory realzaton of a BH model wth two speces of nteractng bosons, whch has been studed theoretcally [ 9] and by Monte Carlo smulatons []. Systems of cold alkal atoms, wth nuclear spn jamshd@physcs.sc.ernet.n Also at Centre for Condensed Matter heory, Department of Physcs, Indan Insttute of Scence, Bangalore 5, Inda rvpa@ungoa.ac.n rahul@physcs.sc.ernet.n; Also at Jawaharlal Nehru Centre For Advanced Scentfc Research, Jakkur, Bangalore, Inda I = / and a hyperfne spn F =, such as, Na, 9 K, and 7 Rb, npurelyoptcal(and notmagnetc)traps, can lead to realzatons of the spn- Bose-Hubbard model, wth spnor condensates [ ], whch can exhbt polar and ferromagnetc superflud phases n addton to MI phases. A dpolar condensate of 5 Cr atoms [5] has been obtaned; to model ths we must nclude long-range nteractons [] n the BH model, n addton to the repulsve nteracton between bosons on the same lattce ste. he frst step n ths drecton can be taken by studyng the extended-bose-hubbard (EBH) model, whch has SF, MI, densty-wave (DW), and supersold (SS) phases [, 7 9]. he study of Ref.[] has ntroduced a smple, transparent, mean-feld theory, whch yelds the phases of the BH model; ths study has also developed a randomphase-approxmaton (RPA) calculaton, whch bulds upon ther mean-feld theory, to obtan the exctaton spectra n the phases of the BH model. We generalze such RPA calculatons so that they can be used for the other bosonc models, namely, the extended-bose- Hubbard (EBH) model, the spn- Bose-Hubbard model, and the Bose-Hubbard model wth two types of bosons. We then use ths RPA to obtan the exctaton spectra n all the phases of these generalzed BH models, the meanfeld theores for whch we have dscussed n Refs.[, 9]. Ourmangoalsto provdeaunfed treatment ofrpa exctaton spectra n the generalzed BH models we have

2 mentoned above. here have been a few studes of exctaton spectra n some of the phases of these models; not all of them are formulated n the way we descrbe here. We dscuss the relaton of our work wth other studes n the last Secton of ths paper. InaddtontoprovdngaunfedtreatmentofRPAexctaton spectra for the three BH models above, our study yelds several nterestng results, whch we summarze below, before we proceed to the detals of our work. Our RPA yelds exctaton spectra, whch have gaps n all the MI phases, n all these BH models, and gaps n the DW phases n the EBH model. hese spectra are gapless n all the SF phases n these models and n the SS phases n the EBHmodel. We obtan the dependence of gaps and (b) the sound velocty u s on the parameters of these models. In partcular, we examne and u s as these systems go through phase transtons. We fnd, e.g., that, at the SF-MI transtons n the spn- BH model, u s goes to zero contnuously (dscontnuously) for MI phases wth an odd (even) number of bosons per ste; ths s consstent wth the natures of these transton (contnuous or dscontnuous) n the mean-feld theory for the spn- BH model []. In the SF phases of these models, our exctaton spectra agree ualtatvely, at weak couplngs, wth those that can be obtaned from Gross-Ptaevsk-type models. For example, our exctaton spectra are ualtatvely smlar to those obtaned for ferromagnetc and polar superfluds n [], whch uses a spn- generalzaton of the Gross-Ptaevsk euaton. he remanng part of ths paper s organzed as follows. InSec.IIwedefnethemodelswestudy, gvetheelements of the mean-feld theory that we use to obtan the phases of these models, and then show how to develop, and then close at the level of the RPA, the euatons of moton for the Green functons for all these models. In Sec. III we present our plots of RPA exctaton spectra for representatve parameter values n these BH models. he concludng Secton contans a dscusson of our results and a comparson of these wth earler studes of such exctaton spectra. II. MODELS, MEAN-FIELD HEORY, AND HE RANDOM-PHASE APPROIMAION In the frst subsecton IIA below we defne the BH model wth two speces of bosons, the spn- BH model, and the EBH model. he next subsecton IIB outlnes our mean-feld theores for these models, prncpally to establsh notatons that are reured for the development of the RPA, whch we present n the last subsecton IIC. Detals of our mean-feld theores have been dscussed n [,, 9]. A. Models A Bose-Hubbard (BH) model, wth two types of bosons, s defned by the followng Hamltonan: H z = t a (a z a j +h.c.) t b (b z b j +h.c.) + U a z + U ab z <,j> <,j> ˆn a (ˆn a )+ U b z ˆn aˆn b µ a z ˆn b (ˆn b ) ˆn a µ b z ˆn b ; () the frst and second terms represent, respectvely, the hoppng of bosons of types a and b between the nearestneghbor pars of stes <,j >, wth hoppng ampltudes t a and t b ; here a, a, and ˆn a a a and b, b, and ˆn b b b are, respectvely, boson creaton, annhlaton, andnumberoperatorsatthestesofad-dmensonalhypercubc lattce (we present exctaton spectra for d = ) for the two bosonc speces. For smplcty, we restrct ourselves to the case t a = t b = t, and, to set the scale of energy, we use zt =, where z = d s the nearest-neghbor coordnaton number. he thrd and fourth terms account for the onste nteractons of bosons of a gven type, wth energes U a and U b, respectvely, whereas the ffth term, wth energy U ab, arses because of the onste nteractons between bosons of types a and b. We have two chemcal-potental terms, µ a and µ b, whch control, respectvely, the total number of bosons of speces a and b. We also study the followng spn- BH Hamltonan[] on a d dmensonal hypercubc lattce wth stes : H zt = <,j>,σ + U zt (a,σ a j,σ +h.c.)+ U zt ˆn (ˆn ) ( F ˆn ) µ ˆn, () zt where spn- bosons hop between the nearest-neghbor pars of stes <,j > wth ampltude t, the spn ndex σ can be,,, a,σ and a,σ are, respectvely, steand spn-dependent boson creaton and annhlaton operators, and the number operator ˆn σ a,σ a,σ; the total number operator at ste s ˆn σ ˆn,σ, and F = σ,σ a,σf σ,σ a,σ, wth F σ,σ standard spn- matrces. hs model () ncludes, n addton to the onste repulson U, an energy U, for nonzero spn confguratons on a ste, whch arses from the dfference between the scatterng lengths for S = and S = channels []. he chemcal potental µ controls the total number of bosons. he Hamltonan for the EBH model s H zt = (a z a j +h.c.)+ U ˆn (ˆn ) zt + V zt <,j> <,j> ˆn ˆn j µ ˆn, () zt

3 where t s the ampltude for a boson to hop from ste to ts nearest-neghbor ste j, z s the nearest-neghbor coordnaton number, <, j > are nearest-neghbor pars of stes, h.c. denotes the Hermtan conjugate, a, a, and ˆn a a are, respectvely, boson creaton, annhlaton, and number operators at the ste, the repulsve potental between bosons on the same ste s U, the chemcal potental µ controls the total number of bosons, and V s the repulsve nteracton between bosons on nearest-neghbor stes. o make a detaled comparson of our results wth experments, the parameters of the BH models must be related to expermental ones []. For the smple BH model U ths s done as follows: zt = π a s V z a exp( E r ), where E r s the recol energy, V the strength of the lattce potental, a s (= 5.5 nm for 7 Rb) the s-wave scatterng coeffcent, a = λ/ the optcal-lattce constant, and λ = 5nm the wavelengthofthe laserused to createthe optcal-lattce; typcally V E r. If we use ths expermental parametrzaton, we scale all the energes by E r. (In ths paper, we set zt =,.e., we measure all energes n unts of zt.) A two-speces BH model, has been realzed n an optcal lattce by usng ellptcally polarzed lght. By changng the polarzaton angle t s possble to shft the lattces wth respect to each other, and, thereby, control the nteractons between the two speces of bosons and also ther hoppng ampltudes (see Refs. [ ] for detals). he spn-dependent term n the spn- BH model follows from the dfference between the scatterng lengths a and a, for S = and S = channels [], respectvely. hese lengths yeld U = π (a + a )/M and U = π (a a )/M, where M s the mass of the atom []. If we consder Na, a = 5.7a B and a = 9.a B, where a B s the Bohr radus, so U > ; n contrast, for 7 Rb, a = (7 ± )a B and a = (±)a B, so U can be negatve. For the EBH case, the relaton of our parameters to those n dpolar bose systems [5, ] s not straghtforward because of long-range nteractons, but we can use the followng estmates: t = w (r r )[ m +V l (r)]w(r r j )d r, () where and j are nearest-neghbor stes, w are Wanner functons, and V l (r) = α=x,y,z V α cos (k α α) s the optcal-lattce potental wth wavevector k. Furthermore, U = U = w(r r ) V nt (r r ) w(r r ) d r d r (5) and V = U <j> = w(r r ) V nt (r r ) w(r r j ) d r d r, () wth V nt = D cos θ r r + π a s M δ(r r ), (7) where D s the dpole moment, a s s the s-wave scatterng constant, and M s the mass of the atom. he s-wave scatterng constant of Chromum s a( 5 Cr) = (7 ± 9)a and a( 5 Cr) = ( ± 5)a, where a =.5 nm []. B. Mean-feld theory We use a homogeneous mean-feld theory for these models because we do not nclude a uadratc confnng potental. (Wth such a potental, we must use an nhomogeneous verson of ths mean-feld theory [, 9].) Conventonal mean-feld theores ntroduce a decouplng approxmaton that reduces a model wth nteractng bosons or fermons to an effectve, nonnteractng problem, whch can be solved easly because the effectve Hamltonan s uadratc n boson or fermon operators. By contrast, the mean-feld theores of [,, 9] decouple the hoppng terms, n the BH models defned above; these hoppng terms are uadratc n boson operators, so, after ths decouplng, we obtan effectve, one-ste Hamltonans, whch can be dagonalzed numercally. For the two-speces BH model E.(), our mean-feld theory [9] obtans an effectve, one-ste problem by decouplng the two hoppng terms as follows: a a j a a j +a a j a a j ; b b j b b j +b b j b b j ; () the superflud order parameters for the ste for bosons of types a and b are ψ a a and ψ b b, respectvely. he approxmaton E.() can now be used to wrte the Hamltonan() as a sum over sngle-ste, meanfeld Hamltonans H MF that are gven below: H MF zt = U a zt ˆn a(ˆn a ) µ a zt ˆn a (φ a a +φ aa )+ψaφ a U b + zt ˆn b(ˆn b ) µ ztˆn b b (φ b b +φ b b )+ψb φ b + U ab z ˆn aˆn b. (9) Here, φ a z δ ψ a+δ and φ b z δ ψ b+δ, where δ labels the nearest neghbors of the ste. For the homogeneous case, the onste chemcal potentals are µ a and µ b, for all, so ρ a = ρ a, ρ b = ρ b, ψ a = ψ a, and ψ b = ψ b are ndependent of. he mean-feld theory for the spn- BH model follows along smlar lnes []: a,σ a j,σ a,σ a j,σ +a,σ a j,σ a,σ a j,σ ()

4 and H MF zt = U ztˆn (ˆn )+ U zt ( F ˆn ) ztˆn µ (φ,σ a,σ +φ,σ a,σ)+ ψ,σ φ,σ. () σ σ Here, we use the followng superflud order parameters: ψ,σ a,σ a,σ ; () and φ,σ z δ ψ (+δ),σ, where δ labels the z nearest neghbors of the ste ; recall, furthermore, that σ can assume the values,,, and ˆn σ a,σ a,σ, ˆn σ ˆn,σ, and F = σ,σ a,σ F σ,σ a,σ, wth F σ,σ standard spn- matrces. For the homogeneous case, the chemcal potental µ does not depend on, so ρ,σ = ρ σ and ψ,σ = ψ σ are also ndependent of. In the EBH model wth V >, we decouple the frst and thrd terms of E.() to obtan an effectve one-ste problem [], whch neglects uadratc devatons from eulbrum values (denoted by angular brackets). he two approxmatons we use are as follows: a a j a a j +a a j a a j ; ˆn ˆn j ˆn ˆn j + ˆn ˆn j ˆn ˆn j ; () the superflud order parameter and the local densty for the ste are, respectvely, ψ a and ρ ˆn. he approxmaton () can now be used to wrte the Hamltonan () as a sum over sngle-ste, mean-feld Hamltonans H MF as follows: H MF H MF, H MF U zt ztˆn (ˆn ) ztˆn µ (φ a +φ a ) + (ψ φ +ψ φ )+ V t (ˆn ρ ρ ρ ), () where the superscrpt MF stands for mean-feld, and φ z δ ψ +δ, ρ z δ ρ +δ, and δ labels the z nearest neghbors of the ste. We can have densty-wave (DW) and supersold (SS) phases, so our order parameters should allow for such phases. In the hypercubc lattces we consder, there are two sublattces A and B. Eachste onthe A(B)sublattcehasz nearestneghbors, each one of whch belongs to the B (A) sublattce; thus, ψ = ψ A and ρ = ρ A, f A, andψ = ψ B and ρ = ρ B, f B, whereas φ = ψ B and ρ = ρ B, f A, and φ = ψ A and ρ = ρ A, f B. he chemcal potentals, whch are conjugateto ρ A and ρ B, respectvely, areµ A, f A, and µ B, f B; smlarly, we can defne creaton, annhlaton, and number operators for each sublattce. he mean-feld Hamltonan () can now be wrtten as HAB MF HMF A +HMF B, (5) where H MF A zt H MF B zt (a A ψ B +a A ψ B)+ (ψ Aψ B +ψ Aψ B ) + V t (ˆn Aρ B ρ A ρ B )+ U ztˆn A(ˆn A ) µ A zt ˆn A; () (a B ψ A +a B ψ A)+ (ψ Bψ A +ψ Bψ A ) + V t (ˆn Bρ A ρ B ρ A )+ U ztˆn B(ˆn B ) µ B zt ˆn B. (7) C. Random-Phase Approxmaton (RPA) Exctaton Spectra We now present a systematc method for developng the RPA for generalzed BH models of the types dscussed above. Our RPA s based on that of Ref. [], for the smple BH model, whch s n turn a generalzaton of the work of Ref. [5] for the spn model. Such RPA calculatons use the mean-feld order parameters defned above and the egenstates of our mean-feld Hamltonans of Es.(9), (), and (5). We begn by defnng the projecton operators [5] L αα = α α, where s the ste ndex and α are the egenstates of the mean-feld Hamltonan. Any sngle-ste operator Ô can be expressed as Ô αα α Ô α L αα. In the next three subsectons, we obtan, explctly, the euatons of moton for the Green functons for the three generalzed BH models we consder and show how to close them n the RPA (see the Appendx for detals).. RPA Exctaton Spectra for the two-speces BH model If we use Ô αα α Ô α L αα, E.() becomes (see the Appendx) H = Vαα L αα αα,ββ L αα Lj ββ, αα where and <,j>,αα,ββ j V αα = (µ a + U a ) α ˆn a α U a α ˆn a α () +(µ b + U b ) α ˆn b α U b α ˆn b α U ab α ˆn aˆn b α (9) j αα,ββ = α â α jβ â jβ + α â α jβ â jβ + α ˆb α jβ ˆb jβ + α ˆb α jβ ˆb jβ. ()

5 5 In the RPA, averages of products of operators are replaced by products of ther averages (at fnte temperature we have thermal averages, but we restrct ourselves to zero temperature here), so the euatons of moton for the Green functons, for ths BH model wth two speces of bosons, namely, G j αα,ββ (t,t ) = θ(t t ) [L αα (t);lj ββ (t )] () become, after Fourer transforms over space and tme, the followng lnear euatons, whch can be nverted easly: (ω ω α +ω α )G αα,ββ (,ω)+ () P αα αα,νµ()g,ββ (,ω) = π P αα δ αβ δ βα, where P αα = L αα L α α (we have suppressed the ste ndces, as we are usng the homogeneous mean-feld theory) and ω α = V αα + β αα,ββ( = ) L ββ, and ω and are, respectvely, the freuency and the wave vector. he propagator G j αα,ββ (t,t ) s the ampltude for the j th ste to flp between states β and β at t gven that the th ste has flpped between states α and α at t; and L αα represents the probablty that the ste s n the state α (we have suppressed ste ndces). he poles of these Green functons gve the dfferent branches of the exctaton spectra n the RPA for dfferent phases n the BH model wth two speces of bosons.. RPA Exctaton Spectra for the Spn- BH Model o obtan the RPA exctaton spectrum for the spn- BH model we use Ô αα α Ô α L αα, so E.() becomes (see the Appendx) H = Vαα L αα αα where and j αα,ββ L αα Lj ββ, <,j>,αα,ββ V αα = (µ+ U ) α ˆn α U α ˆn α () U α ( F ˆn) α () j αα,ββ = α â σ α jβ â σ jβ + α â σ α jβ â σ jβ, (5) where ˆn σ ˆn σ, ˆn σ a σ a σ, F = σ,σ a σ F σ,σ a σ, wth F σ,σ standard spn- matrces, and α are the mean-feld egenstates. Agan, by replacng averages of products of operators by products of ther averages, we obtan the Green functons, for ths spn- BH model, namely, G j αα,ββ (t,t ) = θ(t t ) [L αα (t);lj ββ (t )] () and the RPA euaton for ts Fourer transforms: (ω ω α +ω α )G αα,ββ (,ω) +P αα αα,νµ()g,ββ (,ω) = π P αα δ αβ δ βα, (7) where P αα = L αα L α α, and ω α = V αα + β αα,ββ( = ) L ββ. he poles of these Greens functons gve the dfferent branches of the exctaton spectra n the RPA for dfferent phases n the spn- BH model.. RPA Exctaton Spectrum for the Extended Bose-Hubbard Model o obtan the RPA exctaton spectrum for the EBH model, we wrte E.() as (see the Appendx) H = Vαα L αα αα where and j αα,ββ L αα Lj ββ, <,j>,αα,ββ () V αα = (µ+ U ) α ˆn α U α ˆn α (9) j αα,ββ = α â α jβ â jβ + α â α jβ â jβ V α ˆn α jβ ˆn jβ. () hegreenfunctons, forths EBHmodel, wththe two sublattces A and B, are G j αα,ββ (t,t ) = θ(t t ) [L αα (t);lj ββ (t )]. () In the RPA, the euatons ofmoton for ts Fourertransforms are (ω ω α +ω α )G AA αα,ββ (,ω) +P αα αα,νµ()g BA,ββ (,ω) = π P αα δ αβ δ βα () (ω ω α +ω α )G BA αα,ββ (,ω) +P αα αα,νµ()g AA,ββ (,ω) =. () he poles of these Green functons gve the dfferent branches of the RPA exctaton spectra for the dfferent phases of the EBH model. For each value of the freuency ω and the wave vector, we obtan the egenvalues of a (N s ) (N s ) matrx that multples the Green functon, where N s s the dmenson of the number-state bass that we use [,, 7]. hus, we obtan the RPA exctaton spectra for these generalzed BH models.

6 III. RESULS We present, n the next three subsectons, representatve results from our calculatons of the RPA exctaton spectra for the phases of the BH model wth two speces of bosons, the spn- BH model, and the EBH model. In all our plots, we assume that our BH models are defned on a smple suare lattce, and the wave vector moves along certan hgh-symmetry drectons n the recprocal lattce of the smple suare lattce (see, e.g., Ref. []), n partcular, from Γ ( x =, y = ) to ( x = π, y = ), from to ( x = π, y = π), and then back from to Γ. (b) he results of our RPA calculatons for exctaton spectra n the two-speces BH model are shown by the representatve plots n Fgs. -. Frst we consder the exctatonspectranthemiphases,wthtotaldenstyofbosons ρ = ρ a + ρ b =, and, respectvely, n Fgs. (b)- (d). We also plot ρ versus µ a for U a =, U b =.9U a, U ab =.U a and µ b = µ a n Fg. and mark, wth vertcal lnes, the µ a values for whch the spectra was obtaned n Fgs. (b)-(d). he exctaton spectra n the MI phase wth ρ =, gven n Fg. (b), show four exctatons, one of whch s a hole exctaton and the rest are partcle exctatons (see E. (B) from Appendx for detals). However, for ρ = and ρ =, the spectra, plotted, respectvely, n Fgs. (c) and (d), show four branches; these correspond to two-hole and two-partcle exctatons, whch are the solutons of Es. B and B. For ρ =, two exctatons are almost degenerate for the parameters we have consdered. In all these cases we also obtan dspersonless modes, whch are ndependent of modes; we descrbed these n the Appendx so we do not show them here. (c) a () 5 (d) = 5 () 5 () (b) () 5 () = FIG.. (Color onlne) Plots of the exctaton energy ω() versus the wave number (from Γ to, then to, and fnally to Γ, for a smple suare lattce []) for the BH model wth two speces of bosons. hese plots show representatve exctaton spectra n the MI phases. otal densty ρ versus µ a for U a =, U b =.9U a, µ a = µ b and U ab =.U a. he vertcal lnes n the Mott plateaux represent the µ a values for whch the exctaton spectra were obtaned for (b) ρ =, (c) ρ =, and (d) ρ =, wth all other parameters the same n all cases. For the MI phase wth ρ =, the exctaton spectra consst of one-hole and three-partcle exctatons. In contrast, there are two- hole and partcle exctatons for ρ = and. For the parameters chosen here, the lowest-energy exctaton s degenerate for ρ =. A. RPA Exctaton Spectra for the Bose-Hubbard Model wth two speces of bosons = FIG.. (Color onlne) Plots of exctaton spectra when both types of bosons are n the SF phase and (b) when only one type s n the SF and the other type of boson s n MI phase. he exctaton spectra for the case where both a-type and b-type bosons are n the SF phase are gven n the Fg. whch show two gapless exctatons; lnear n wave number near the pont Γ and exctatons whch have fnte gap at the pont Γ. Smlarly, the exctaton spectra for the case where only b-type bosons are n the SF phase and a-type bosons are n the MI phase are gven n Fg. (b); we see only one gapless exctaton correspondng to the SF phase of b-type bosons; and all the other branches of the exctaton spectra have fnte gaps at the pont Γ. We also show some of the dspersonless spectra. In Fg., we present plots versus µ a = µ b of ρ a (red curve) and ρ b (black curve), for bosons of types a andb, foru a =, U b =.9U a, µ a = µ b andu ab =.U a, wth plateaux at ρ a = and/or ρ b =. In a gapless, SF phase, the speed of sound u s follows [] from the slope of the exctaton spectrum near the pont Γ. In Fg. (b), we show how the speed of sound u s behaves, for compo-

7 7 a, b u S ln ln(u-u C ) () (c) (b) (d) () a 9 a a FIG.. (Color onlne) Plots of ρ a (red lne) and ρ b (black lne), for bosons of types a and b, versus µ a for U a =, U b =.9U a, µ a = µ b, andu ab =.U a; notetheplateauxatnteger values of ρ a and ρ b. (b) he speed of sound u s versus ρ a, for U a =, U b =.9U a, ρ a = ρ b, and U ab =.U; u s s zero n the MI phase and fnte n the SF phase. (c) he gap versus µ a, for U a =, U b =.9U a, µ a = µ b and U ab =.U a; the nset shows the power-law behavor of the gap n the vcnty of the transton at U c =.. nents a (red dashed lne) and b (black dashed lne wth trangles), as functons of µ a, for the same parameters as n Fg. ; as we expect, u s s zero n the MI phases and postve n the SF phases. From the exctaton spectra we can obtan the gap at the pont Γ. In Fg. (c), we show plots of versus µ a for the same parameters as n Fgs. and (b); as we expect, s zero n the SF phases (for bosons oftypes a and b) and fnte n the MI phases. In the nset of Fg. (c) we shows that approaches zero as (U U c ) /, where U c s the crtcal value of U at whch the MI-SF transtons occur at these parameter values. We now consder the frst few low-energy exctaton spectra n the SF phase; we hold ρ =.75, but vary the U ab n Fgs. -(c). In all these fgures U a = 9, U b =.9U a and µ a = µ b. In Fg., U ab /U a =, we see that the gapless mode s almost degenerate; and the next branch n the exctaton spectrum s dspersonless; we do not show hgh-energy spectra here. As we ncrease U ab /U a, Fg. (b) for U ab /U a =. and Fg. (c) for U ab /U a =., the degeneracy of the gapless exctaton spectrum s lfted and the spectra move away from each other. he non-dspersve mode now shows dsperson and moves to hgher energes. In Fg. (d), we llustrate the exctaton spectra for U ab /U a < ; we FIG.. (Color onlne) Plots of the exctaton energy ω() versus the wave number for the BH model wth two speces of bosons, for ρ =.75, U a = 9, U b =.9U a, and U ab =., (b) U ab =.U a, (c) U ab =.U a and (d) U ab =.U a. see agan that the gapless modes move away from each other; and we observe that the exctaton freuences of other modes are dfferent compared to ther values for U ab /U a >. Our results are ualtatvely smlar to those of the Gutzwller-based, exctaton-spectrum study of Ref. [7], whch presents k x = k y scans through the recprocal lattce of a two-dmensonal suare lattce. B. RPA Exctaton Spectra for the Spn- Bose-Hubbard Model We now gve representatve plots of our RPA exctaton spectra for the spn- BH model, n Fgs. 5 -(c) for U /U =, U /U =., and U /U =., whch yeld, respectvely, a pure SF wth no spn nteractons, a polar SF, and a ferromagnetc SF for the parameter values gven n the fgure capton []. he exctaton spectra are ualtatvely dfferent n these three cases. In the SF (Fg. 5 ) phase there are two gapless, degenerate modes, shown n black, whose freuences approach zero uadratcally at small wave numbers (near Γ); there s another mode, shown n red, whose energy goes to zero lnearly at small. All the other exctatons spectra have fnte gap at Γ; some of these spectra are degenerate and some are ndependent; e.g., the spectra shown n green, blue, and maroon are ndependent and have degeneraces,, and, respectvely. he branch shown n navy blue s weakly dependent and s doubly

8 5 5 (b) 5 (c) (b) () () FIG. 5. (Color onlne) Plots of the exctaton energy ω versus the wave number for the spn- BH model: superflud (U =, µ = 5 and U /U = ), (b) polar-superflud (U =, µ = 5 and U /U =.) and (c) ferromagnetcsuperflud (U =, µ = 5 and U /U =.) phases. In and (b) the lowest modes are two gapless, degenerate modes, shown n black, whose energy approaches zero uadratcally at small wave numbers (near Γ). In (c) we see clearly the splttng of the degenerate modes, whch now have a gap and then a uadratc dependence on near Γ. In all three cases, there s another mode, shown n red, whose energy goes to zero lnearly at small. All the hgher exctatons have fnte gaps at Γ and are hghly degenerate. he degeneraces of these modes are partally lfted n the case of polar and ferromagnetc SF phases. FIG.. (Color onlne) Plots of the exctaton freuency ω() versus the wave number for the spn- BH model n the MI phase wth ρ = : U =, µ =.5 and U /U =. and (b) U =, µ =.5 and U /U =.. he MI ground state s degenerate; and these excted states do not couple to the ground state through partcle or hole exctatons; ths gves the the dspersonless mode at ω() =. he mode shown n red s a hole exctaton; all the hgher-energy modes shown are partcle exctatons. degenerate. hese degeneraces are partally lfted when U s fnte. In the polar SF (Fg. 5 (b)), the lowest modes, shown n black, are doubly degenerate and have a uadratc dependence on as n the case of a pure SF. However, the next branch, whch dsplays a lnear dependence on near the Γ pont, as n the case of a pure SF, becomes approxmately dspersonless at hgher value of, and avods level crossng wth the branch shown n green (these acure a dependence). Note that ths spectrum, dspersonless and wth a trple degeneracy n the pure SF case, has ts degeneracy lfted partally: the spectrum shown n navy blue s non-degenerate and the one shown n blue s doubly degenerate; we obtan smlar behavors n hgher-energy exctaton spectra. In the ferromagnetc SF (Fg. 5 (c)), we see clearly the splttng of the degenerate modes, whch now have a gap and then a uadratc dependence on near Γ; the red curve ndcates the gapless mode, whch approaches Γ lnearly. he degeneraces of hgher modes are also lfted n ths case. he ualtatve features of these spectra agree wth those found for the varous SF phases n the spn-, Gross- Ptaevsk study of Ref. []. We present exctaton spectra n the MI and MI phases n Fgs. and 7, respectvely. Here, MI and MI denote, respectvely, Mott phases wth one and two bosons per ste. For MI, we obtan ualtatvely smlar exctaton spectra for U /U > (Fgs. ) and U /U < (Fgs. (b)). For U /U >, the ground state for MI s trply degenerate(appendx); these degenerate states do not couple to each other throughthe creatonof a partcle or hole. hus, the lowest exctaton, ω() =, represented n black n Fg., s doubly degenerate and dspersonless. he non-degenerate hole exctaton s shown n red. he exctatons, above the hole exctaton, are sx-partcle exctatons, two of whch are dspersonless and degenerate. We fnd smlar exctaton spectra for U /U < (Fgs. (b)). For MI, we obtan ualtatvely dfferent exctaton spectra for U /U > (Fgs. 7) and U /U < (Fgs. 7(b)). For U /U >, the MI ground state s non-degenerate(appendx) because of the formaton of a snglet, so there s no dspersonless mode wth ω() =. he frst excted state of our MF hamltonan (E. ) has fve-fold degeneracy; and the ground state does not couple to these states va partcle or hole exctatons. hs gves us the dspersonless mode (black lne). he dspersve mode wth the lowest energy comprses hole exctatons (red curves); ths mode s trply degenerate. he partcle exctatons are also dspersve (green curve,

9 9 () (b) and U /U =.) u s goes to zero contnuously at the phase-sf-mi transton. he natures of these transtons are consstent wth the predctons of the mean-feld theory [] on whch we base our RPA study. We do not nvestgate magnetc orderng and magnetc exctatons n the MI phase of the spn- model because, as has been noted earler, our mean-feld theory does not allow for any magnetc structure n the MI phase n ths model []. C. RPA Exctaton Spectra for the Extended Bose-Hubbard model FIG. 7. (Color onlne) Plots of the exctaton energy ω() versusthewave numberfor thespn-bhmodel wth ρ = : U =, µ =.5 and U /U =. and (b) U =, µ = 7.5 and U /U =.. u S..5 u S U 7 U We now show representatve RPA exctaton spectra n Fgs. 9, (b), (c), (d), and (e) for SF, DW / (see Ref. []), SS, MI, and DW / phases for V =.U and U =, µ =., (b) U =, µ = 9, (c) U =., µ = 5., (d) U =, µ =, and (e) U =, µ = 5. By the symbol DW n/m, we mean a denstywave phase wth n bosons per ste and m stes per unt cell. From these exctaton spectra we obtan u s, whose dependence on U s shown n Fg. 9 (f), as the system moves from the SF to the SS phase and then to the DW / phase, for V =.U, and ρ total = ; the sharp drop n u s at the SS-DW / boundary ndcates a frst-order transton; the nset shows that u s U n the small-u regme n the SF phase. he RPA exctaton spectra have a gap at the Γ pont when the system s n the DW / phase and ths gap goes to zero as the system approaches the SS phase. In Fg. we plot versus U as the system goes from the DW to the SS phase. Fgure (b) show analogous plots for the SF-MI transton. he behavors of ths gap are what we expect on general grounds,.e., t jumps dscontnuosly to zero at the DW / to SS transton, but t goes to zero contnuously, as (U U c ) /, at the MI to SF transton at the crtcal value U c. hs exponent of / s a mean-feld exponent, whch should be modfed by fluctuatons (ths MI-SF transton les n the d-dmensonal Y unversalty class). FIG.. he sound speed u s versus U n the vcnty of the polar-sf-mi transton for µ = and U /U =., wth (ρ = ); u s jumps to zero at ths frst-order transton; by contrast, (see nset for µ = 5 and U /U =.) u s goes to zero contnuously at the SF-MI transton. degeneracy ), except for the one shown n blue (degeneracy ). Hgher-energy exctatons are mult-partcle-hole exctatons and are dspersonless; we do not show these here. In Fg. we plot the sound speed u s versus U n the vcnty of the polar-sf-mi transton for µ = and U /U =., wth (ρ = ); u s jumps to zero at ths frst-order transton; by contrast, (see nset for µ = 5 IV. CONCLUSIONS We have obtaned representatve exctaton spectra n all phases of three generalzed Bose-Hubbard (BH) models, namely, () a two-speces generalzaton of the spnless BH model, () a spn- BH model, and () the extended Bose-Hubbard model (EBH) for spnless nteractng bosons of one speces. Our study uses the random phase approxmaton (RPA), whch t develops, n a unfed way, by startng from mean-feld theores of the type that we have dscussed n Ref. []. hus, t generalzes the RPA studes ntated n Ref. [] for the smple BH model and contnued, e.g., n Refs. [, 7, 9]. Our study yelds a varety of nterestng results that we have descrbed n detal n the prevous Secton. In

10 () (b) (c).. () (d) (e) u S SF u S (f) SF U /... SS DW / 7 9 U 9 U.5 (b) U ln(u-u C ) ln SF to MI 7 U SS to DW/ FIG. 9. (Color onlne) Plots of the exctaton freuency ω() versus the wave number for the EBH model n U =, V =.U, µ =., and ρ total =, (b) U =, V =.U, µ = 9, ρ total =, (c) U =., V =.U, µ = 5., ρ total =, (d) U =, V =.U, µ =, ρ total =, and (e) U =, V =.U, µ = 5, ρ total =. We show the exctaton spectra n the SF, DW/, SS, MI, and DW/ phases, respectvely. In (f) V =.U, ρ total = ; and we show the dependence of u s on U as the system goes from the SF to the SS and then to the DW/ phase; the last transton s clearly frst order; the nset shows the u s U behavor at low U n the SF phase. partcular, our RPA exctaton spectra show clear gaps n MI phases, n all the models above, n the DW phases n the EBH model, and gapless spectra n all SF phases and the SS phase n the EBH model. From these spectra we have obtaned the dependence of gaps and (b) the sound velocty u s on the parameters of these models. We have also nvestgated the behavors of and u s as these systems go through phase transtons. We fnd that, at the polar SF-MI transtons n the spn- BH model, u s goes to zero contnuously(dscontnuously) for MI phases wth an odd (even) number of bosonsper ste; ths s consstent wth the natures of these transton (contnuous or dscontnuous) n the mean-feld theory for the spn- BH model []. In the SF phases of these models, our exctaton spectra agree ualtatvely, at weak couplngs, wth those that can be obtaned from Gross-Ptaevsk-type models. For example, our RPA exctaton spectra are ualtatvely smlar to those obtaned for ferromagnetc and polar superfluds n [], whch uses a spn- generalzaton of the Gross-Ptaevsk euaton. Exctaton spectra, at the level of the RPA, can also FIG.. Plotsofthegap ntherpaexctatonspectrumas afuncton of U near thess-dw/transton and(b) near the SF-MI transton; the nset shows the power-law behavor of the gap n the vcnty of the transton at U c = 5.7. be obtaned by startng wth a Gutzwller-type meanfeld theory. Such studes, carred out, e.g., for two-boson BH models, the spn- BH models, and the EBH model, should be formally euvalent to our RPA study; however, we are not aware of any study that has shown ths n complete detal. We do not know of one, unfed treatment of Gutzwller-type RPA, lke our unfed RPA study for these three models, so a drect comparson of our results wth these Gutzwller-type studes s not completely straghtforward. Reference [] has evaluated the exctaton spectrum of spn- bosonc atoms n a Mott-nsulator phase. Hou, et al. have nvestgated the uantum phase transton of spn- cold bosons n an optcal lattce wth and wthout an external magnetc feld, []; they also obtan Gutzwller-type RPA exctaton spectra. Prevous studes have analyzed exctaton spectra n the two-speces BH model. In the MI phase, the lowest two branches of the exctaton spectrum have gaps; and they correspond to the partcle- and hole-exctaton modes [, ]. In the SF phase, the exctaton spectrum has onegaplessmode and one mode wth agap [5 9]; the gapless mode arses because of oscllatons of the phase of the SF order parameter (ths s the Bogolubov mode); the lowest mode wth the gap arses because of the ampltude mode n the vcnty of the SF-MI transton at nteger fllngs; ths ampltude mode becomes gapless at the crtcal pont. Our study explores the dependence of the exctaton spectra on model parameters n far greater detal than earler studes.

11 Furthermore, the RPA method can be used to obtan the complete ω and dependence of the Green functon as noted, e.g., n Refs. [, 7] for the smple, spnless BH model. From ths, varous propertes can be calculated, e.g., the momentum dstrbuton functon; the latter has been obtaned n the RPA n [9] for the EBH model. We gve a bref descrpton of such Green-functon studes n the Appendx and ther generalzatons for the three models we consder here. Exctaton spectra have been measured expermentally by Bragg-spectroscopy [5 5] and lattce-ampltudemodulaton [57 59] methods; and these measurements havebeen used to characterzesf and MI phases and the transton between them. We hope our work wll lead to expermental measurements of the spectra of elementary exctatons n the dfferent phases n the physcal realzatons of the generalzed BH models we have dscussed above. Momentum-dstrbuton functons can be calculated, at the RPA level, from the Green functons that we have obtaned. he calculaton s especally smple n a Mott phase, as llustrated, e.g., n Ref. [9]. Such momentum-dstrbuton functons can be compared wth those that are measured expermentally as dscussed, e.g., n Refs. [, 9]. he studes of Refs. [7, 9] also menton some of the lmtatons of the RPA n the context of the BH model; e.g., Ref. [7] notes that ths type of of RPA, when appled to the smple BH model, yelds a volaton of the total densty sum rule. Exctaton spectra can also be measured by Quantum Monte Carlo (QMC) smulaton as shown for the smple BH model n Ref. [7] and a contnuaton model n Ref. []. We hope our study wll lead to QMC studes of the exctaton spectral of the generalzed BH models we have dscussed here. the Unversty of Goa and the Indan Insttute of Scence, for hosptalty durng the perod n whch ths paper was beng wrtten. Appendx A: he Random-Phase Approxmaton. We llustrate below how we obtan the RPA euatons for the Green functons for our Bose-Hubbard(BH) models of Es. (), (), and (), respectvely. he RPA euaton of moton, for any of these BH models use the operators L αα = α α, where α s the egenvector for the mean-feld Hamltonan for ste. Any sngle-ste operator can be expressed as Ô αα α Ô α L αα, so the Hamltonan for these BH models can be wrtten as H m = V,m αα L,m αα j,m αα,ββ L,m αα Lj,m ββ, αα j,αα,ββ (A) where the superscrpt m =, s for two speces of bosons, m = s for the spn- BH model, and m = s for the EBH model. V, αα =(µ a + U a ) α ˆn a α U a α ˆn a α +(µ b + U b ) α ˆn b α U b α ˆn b α U ab α ˆn aˆn b α, (A) V. ACKNOWLEDGMENS We thank DS, UGC, CSIR (Inda)(No. ()//EMR-II) for support and K. Sheshadr for useful dscussons. JMK and RVP thank, respectvely, and V, αα = (µ+ U ) α ˆn α U α ˆn α U α ( F ˆn ) α, V, αα = (µ+ U ) α ˆn α U α ˆn α, (A) (A) j, αα,ββ = α â α jβ â j jβ + α â α jβ â j jβ + α ˆb α jβ ˆb j jβ + α ˆb α jβ ˆb j jβ, (A5) j, αα,ββ = σ { α â,σ α jβ â j,σ jβ + α â,σ α jβ â j,σ jβ }, (A) j, αα,ββ = α â α jβ â j jβ + α â α jβ â j jβ V α ˆn α jβ ˆn j jβ. (A7) he propagatorg j,m αα,ββ (t,t ), s the ampltude for the j th ste to flp between the states β and β at t, gven

12 that the th ste has flpped between states α and α at t. We have G j,m αα,ββ (t,t ) = θ(t t ) [L,m αα (t),lj,m ββ (t )] = L,m αα (t) Lj,m ββ (t ). (A) he euaton of moton for ths Green functon, wth t =, s calculated as follows: where d dt Gj,m αα,ββ d dt L,m αα (t) = (t) = dθ(t) dt [ L,m αα (t),lj,m ββ ] () θ(t) d dt [ L,m αα (t),lj,m ββ () ], (A9) [ ] [ H m,l,m αα (t) L,m αα (t),hm] (A) we have set =, and we obtan the second term on the rght-hand sde of E. (A9) d dt [ L,m ] αα (t),lj,m ββ () = L j,m ββ () [ L,m αα (t),hm] [ L,m αα (t),hm] L j,m ββ () (A) By usng E. (A) n the above euaton, substtutng t n E. (A9) and carryng out the RPA, where we replace thermal averages of products of operators by the products of ther thermal averages, to obtan d dt Gj,m αα,ββ (t) = δ(t)δ α βδ αβ Pαα m +(ωm α ωm α )G j,m αα,ββ (t)+pm αα Fourer transformaton over space and tme for m = and m = now yelds (ω ω m α +ω m α )Gm αα,ββ (,ω)+pm αα l l,m αα,νµ Glj,m,ββ (t). m αα,νµ()g m,ββ,ω) = π Pm αα δ αβ δ βα ; and, for m = wth a bpartte hypercubc lattce, wth two sublattces A and B, we get (ω ωα m +ωα m )GAA αα,ββ (,ω)+pm αα m αα,νµ()g BA,ββ (,ω) = π Pm αα δ αβ δ βα, (ω ωα m +ωα m )GBA αα,ββ (,ω)+pm αα m αα,νµ()g AA,ββ (,ω) =, (A5) (A) (A) (A) where αα m,ββ () = ǫ ( j,m αα,ββ + j,m ββ,αα ), ǫ t j=x,y cos( j), Pαα m = Lm αα Lm α α and ωm α = Vαα m + β m αα,ββ ( = ) Lm ββ, and ω and are, respectvely, the freuency and wave vector of the exctaton. he poles of these Green functons gve the dfferent branches of the exctaton spectrum n the RPA. We can relate the Green functons calculated above to the sngle-partcle Green functon g,j (t) = θ(t) [â (t),â j ]. We frst llustrate ths for the smple BH model by followng the treatment of Ref. []. We defne the sngle-partcle Green functons g (,ω) = y α,α y β,β G αα,ββ (,ω), (A) αα, ββ g (,ω) = y α,α y β,β G αα,ββ (,ω), (A7) αα, ββ where y α,α,α â,α and y α,α,α â,α are ndependent of the poston. It s easy to show that euaton of moton for the green functons g (,ω) and g (,ω) can be wrtten as ( ) ( ) ǫ A ǫ A g (,ω) ǫ A ǫ A g (,ω) = ( A A ), (A) where

13 ( ) A A A A = α ( y,α y α, y,αy α, y,α y α, y,α y α, ω +δ +(ω ω α ) ) α ( ) y,α y α, y,α y α, y,α y α, y,αy α,, (A9) ω +δ +(ω α ω ) represents the ground state and the summaton s over all the excted states α. We gve below the analogs of Es.(A)-(A9) for the BH model wth two speces of bosons n Es.(A) and (A) and for the spn- BH model n E.(A). For a- type bosons, we defne the sngle-partcle Green functons g MN (,ω) = y M α,α yn β,β G αα,ββ (,ω), αα, ββ (A) g MN (,ω) = y M α,α y N β,β G αα,ββ (,ω), αα, ββ (A) where the superscrpts M and N can be a or b, for the two-speces BH model, or σ and σ for the spn- BH model. For the two-speces BH model, the euaton of motonforthegreenfunctonsg MN (,ω)andg MN (, ω) can be wrtten as ǫ p A aa ǫ p A ba ǫ p A aa ǫ p A ba ǫ p A ab ǫ p A bb ǫ p A ab ǫ p A bb ǫ p A aa ǫ p A ba ǫ p A aa ǫ p A ba ǫ p A ab ǫ p A bb ǫ p A ab ǫ p A bb ǫ p A bb ǫ p A ab ǫ p A bb ǫ p A ab ǫ p A ba ǫ p A aa ǫ p A ba ǫ p A aa ǫ p A bb ǫ p A ab ǫ p A bb ǫ p A ab ǫ p A ba ǫ p A aa ǫ p A ba ǫ p A aa (,ω) (,ω) (,ω) (,ω) g aa g ba g aa g ba (,ω) (,ω) (,ω) (,ω) g bb g ab g bb g ab = = A aa A ab A aa A ab A bb A ba A bb A ba,. (A) (A) For the spn- BH model we have: ( ǫp A σσ ǫ p σ Aσσ ǫ p A σ σ ǫ p A σ σ ) ( g σσ (,ω) g σ σ (, ω) ) = ( ) A σσ A σ σ. (A) Fnally we can wrte ( ) ( ) ( A MN A MN A MN A MN ) = α y N,α y M α, y N,α ym α, y N,α y M α, y N,α ym α, ω +δ +(ω ω α ) α y M,α yn α, y M,α yn α, y M,α y N α, ym,α y N α, ω +δ +(ω α ω ). (A5) Appendx B: Calculaton of exctaton spectra analytcally n the MI phases for BH models: We calculate the exctaton spectra for the MI phases of BH models analytcally. We begn wth the BH model. In the MI phase, the superflud order parameter s zero, the densty of bosons ρ s eual to an nteger, and the egenstates of the mean-feld Hamltonan are gven by the Fock states,,α,n = n! (a )n vacuum. he presence of P αα L αα L α α n the RPA euaton for the Green functon (we have suppressed the ste ndex ) (ω ω α +ω α )G αα,ββ (,ω)+ (B) P αα αα,νµ()g,ββ (,ω) = π P αα δ αβ δ βα

14 yeld nonzero Green functon G αα,ββ (,ω) only f α or α are ground states. Let us obtan the Green functon for the ground state wth densty of bosons eual to ρ, G αα,α α(,ω) G nn+,n+n (,ω) wth α = n = ρ and α = n +. he hoppng matrx n E. (B) s gven by αα,νµ = ( α a α ν a µ + α a α ν a µ ). (B) Snce α = n = ρ and α = n+, the nonzero term n the above euaton have µ = n, ν = n + or µ = n, ν = n. hus the RPA euaton for the Green functon G n n+,n+ n (,ω) s [ω ω n +ω n+ +(n+)ǫ ]G nn+,n+n (,ω)+ ǫ n(n+)gn n,n+ n (,ω) = π.(b) Smlarly the RPA euaton for the Green functon G n n,n+ n (,ω) s ǫ n(n+)gnn+,n+n (,ω) +[ω ω n +ω n nǫ ]G n n,n+ n (,ω) =. (B) For the MI phase, ω n = µn U n(n ), so by solvng the coupled Es. (B) and (B) we obtan the followng exctaton spectra for the MI phase wth densty eual to ρ: ω ± = [ (µ U(ρ )+ǫ ) ± ǫ ǫ U(ρ+)+U ]; (B5) here the plus (mnus) sgn corresponds to the partcle (hole) exctaton. he MI phases can also dsplay dspersonless exctaton spectra; we llustrate ths by evaluatng the Green functon G αα,α α(,ω) G n n+,n+ n (,ω) wth α = n = ρ and α = n +. It s easly seen from E.(B)that, forthesevaluesofαandα, αα,νµ() = for all values of µ and ν. hus the RPA euaton for G n n+,n+ n (,ω) s gven by (ω ω n +ω n+ )G n n+,n+ n (,ω) = π. (B) hs euaton gves -ndependent exctaton spectra ω = ω n ω n+,.e., mult-partcle (or mult-hole) exctaton spectra can be dspersonless. For the BH model wth two-speces of bosons, the egenstates of the mean-feld Hamltonan for the MI phase, wth even total densty of bosons ρ = n a + n b are gven by α n a n b = na! n b! (a ) na (b ) n b vacuum. (B7) o obtan the partcle(hole) exctaton for the a-type bosons n the MI phase wth densty of bosons for a-type and b-type, respectvely, eual to ρ a and ρ b, we calculate the Green functon G αα,α α(,ω) G a n an a+ (,ω) wth α = n a = ρ a n b = ρ b and α = n a + n b usng the E. (). In ths case, the hoppng matrx element αα,νµ() gven by E. () s non-zero f { µ, ν } eual to { n a n b, n a + n b } or { n a n b, n a n b }. herefore, the RPA euaton for G a n a n a+(,ω) s gven by [ω ω na n b +ω na+ n b +(n a +)ǫ ]G a n an a+ (,ω)+ ǫ na (n a +)G a n a n a (,ω) = π (B) where the Green functon G a n a n a (,ω) G αα,β β(,ω) wth α = n a n b, α = n a n b, β = n a + n b and β = n a n b and ts RPA euaton s gven by ǫ na (n a +)G a n a n a+(,ω) +[ω ω na n b +ω na n b n a ǫ ]G a n a n a (,ω) =. (B9) Here ω nan b = µ a n a Ua n a(n a )+µ b n b U b n b (n b ) U ab n a n b. By solvng coupled Es. (B) and (B9) we obtan the exctaton spectra for the MI phase wth the denstes of a-type and b-type bosons, respectvely, eual to ρ a and ρ b : ω ± = [ (µ a U a (ρ a ) U ab ρ b +ǫ ) ± ǫ ǫ U a (ρ a +)+Ua ] (B) where the plus (mnus) sgn corresponds to the partcle (hole) exctaton. hs expresson has been obtaned earler by usng Gutzwller approxmaton [7]. In the lmt U ab =, ths expresson matches wth that for BH model (E. (B5)). Smlarly, for the partcle (hole) exctatons for b-type partcles we obtan ω ± = [ (µ b U b (ρ b ) U ab ρ a +ǫ ) ± ǫ ǫ U b (ρ b +)+Ub ]. (B) For the MI phase wth total densty of bosons ρ =, the frst few egenstates of the mean-feld Hamltonan are gven by α = C +C, α = C C, α =, α =, α = D + D, and α 5 = D D. Here C, C, D, and D depend on the values of U a, U b and U ab. he ground state α couples to mean-feld egenstates () α va the annhlaton of a-type or b-type bosons () α, α and α 5 va the creaton of a-type or b-type bosons. Startng wth the Green functon G αα,α α the coupled RPA euatons can be wrtten n a matrx form gven by

15 ω ω α,α +α α,α α α α,α α α α,α α α α,α α 5 α α,α α ω ω α,α α α,α α α α,α α α α,α α 5 α α,α α α α,α α ω ω α,α α α,α α α α,α α 5 5 α 5α,α α α 5α,α α α 5α,α α ω ω α5,α α 5α,α α 5 π =, G αα,α α (,ω) G αα,α α (,ω) G αα,α α (,ω) G α5α,α α (,ω) (B) where ω α,α j = ω α ω αj and α α j,α jα = αα j,α jα (). Solutons of ths euaton gve one-hole and three-partcle exctaton spectra. he exctaton spectra of the MI phase of the spn- BH model can also be obtaned n a smlar manner. he egenstates of the mean-feld Hamltonan for the MI phase wth total densty of bosons ρ = n +n +n are gven by α n n n = n!n!n! (a, )n (a, )n (a, )n vacuum. (B) Frst we dscuss the MI phase wth ρ = where we fnd the ground state s trply degenerate and the egenstates are lnear combnatons of Fock states, and. he frst excted state s and s non-degenerate and ts value s ndependent of U. he ground state s coupled to ths state va the annhlaton of a boson, whch gves rse to a hole exctaton n the spectrum. he next excted states have degeneracy of sx f U = ; ths s partally lfted f U s fnte. hese egenstates are lnear combnatons of Fock states,,,, and. he ground state couples to these states va the creaton of a boson; ths gves rse to sx-partcle exctatons n the spectra. hese are the relevant exctatons n the ρ = MI phase. All the hgher energy egenstates have total boson numbermorethantwoandwllnotcoupletothegroundstate through the hoppng matrx and thus gves dspersonless mult-partcle exctaton spectra. For ρ =, the ground state of the spn- mean-feld Hamltonan s non-degenerate n the case of U >. However, t has a degeneracy of fve and sx, respectvely, for U < and U =. hese states are lnear combnatons of the two-boson Fock states,,,, and. he ground state couples to the one-boson Fock states, and va the hoppng matrx; ths yelds three-hole exctatons. Smlarly, the ground state couples to three-boson Fock states (a total of states) va the creaton of bosons. All the hgher exctatons n the mean-feld soluton do not couple to the ground state through the hoppng matrx; and they gve rse to dspersonless mult-partcle exctatons. [] I. Bloch, J. Dalbard, and W. Zwerger, Rev. Mod. Phys., 5 (); M. Lewensten, et al., Adv. n Physcs, 5 (7). [] S. Sachdev, Quantum Phase ranstons (Cambrdge Unversty Press, 999). [] D. Jaksch, et al., Phys. Rev. Lett. (99). [] M. Grener, et al. Nature (London) 5 9 (). [5] M. P. A. Fsher, et al. Phys. Rev. B 5 (99); D. S. Rokhsar and B. G. Kotlar, Phys. Rev. B (99); W. Krauth, M. Caffarel, and J. P. Bouchaud, Phys. Rev. B 5 7 (99). [] K. Sheshadr, et al., Europhys. Lett. 57 (99). [7] W. Krauth, N.rved, andd. Ceperley, Phys. Rev.Lett. 7, 7 (99); N. rved and M. Makvc, Phys. Rev. Lett. 7, 9 (995). [] J. Catan, et al., Phys. Rev. A 77, (R) (). [9] S. rotzky, et al., Scence 9, 95 (). [] A. Wdera, S. rotzky, P. Chenet, S. Föllng, F. Gerber, I. Bloch, V. Grtsev, M. D. Lukn, and E. Demler, Phys. Rev. Lett., (). [] D. M. Weld, P. Medley, H. Myake, D. Hucul, D. E. Prtchard, and W. Ketterle, Phys. Rev. Lett., 5 (9). [] B. Gadway, D. Pertot, R. Remann, and D.Schneble, Phys. Rev. Lett. 5 5 (). [] A. B. Kuklov and B.V. Svstunov Phys. Rev. Lett. 9, (). [] J.-R. Han, Physcs Letters A, (). [5] P. Buonsante, et al., Phys. Rev. Lett., (). [] A. Hu, et al., Phys. Rev. A, 9 (9). [7]. Ozak, I. Danshta and. Nkun arxv:condmat/.7v. [] J.M. Kurdestany, R.V. Pa, and R. Pandt, Ann. Phys. (Berln), () / DOI./andp.7. [9] R.V. Pa, J. M. Kurdestany, K. Sheshadr and R. Pandt, Phys. Rev. B 5, 5 (). []. Rosclde and J. Ignaco Crac, Phys. Rev. Lett. 9, 9 (7). [] See, e.g., H.-J. Mesner, et al., Phys. Rev. Lett., (999). [].L. Ho, Phys. Rev. Lett., 7 (99). [] S. Mukerjee, C. u, and J.E. Moore, Phys. Rev. Lett. 97, (). [] R.V. Pa, K. Sheshadr and R. Pandt, Phys. Rev. B (77)