Ultracold dipolar gases in optical lattices

Author: Trefzger, Christian

Publisher: Universitat Politècnica de Catalunya

Year: 2010

DOI: 10.5821/dissertation-2117-93907

Source: https://upcommons.upc.edu/bitstream/2117/93907/1/TCT1de1.pdf

Ul acold Dipola Gases
in Op ical La ices
Ch is ian T e zge
Thesis Ad iso :
P o . Maciej Lewens ein
Thesis Co-ad iso :
D . Chia a Meno i
2
Gene al in oduc ion
This hesis is a heo e ical wo k, in which we s udy he physics o ul a-cold
dipola bosonic gases in op ical la ices. Such gases consis o bosonic a oms
o molecules, cooled below he quan um degene acy empe a u e, ypically in
he nK ange. In such condi ions, in a h ee-dimensional (3D) ha monic ap,
weakly in e ac ing Bosons condense and o m a Bose-Eins ein Condensa e
(BEC). When a BEC is loaded in o an op ical la ice p oduced by s anding
wa es o lase ligh , new kinds o physical phenomena occu . These sys ems
ealize hen Hubba d- ype models and can be b ough o a s ongly co ela ed
egime.
In 1989, M. Fishe e . al. p edic ed ha he homogeneous Bose-Hubba d
model (BH) exhibi s he Supe ﬂuid-Mo insula o (SF-MI) quan um phase
ansi ion [1]. In 2002 he ansi ion be ween hese wo phases we e obse ed
expe imen ally o he ﬁ s ime in he g oup o I. Bloch, T. Esslinge and
T. H¨
ansch [2]. The expe imen al ealiza ion o a dipola BEC o Ch omium
by he g oup o T. P au [3, 4, 5], and he ecen p og esses in apping and
cooling o dipola molecules by he g oups o D. Jin and J. Ye [6, 7, 8], ha e
opened he pa h owa ds ul a-cold quan um gases wi h dominan dipole
in e ac ions. A na u al e olu ion, and p esen challenge, on he expe imen al
side is hen o load dipola BECs in o op ical la ices and s udy s ongly
co ela ed ul acold dipola la ice gases.
Be o e his PhD wo k, s udies o BH models wi h in e ac ions ex ended o
nea es neighbo s had poin ed ou ha no el quan um phases, like supe solid
(SS) and checke boa d phases (CB) a e expec ed [9, 10, 11, 12]. Due o he
long- ange cha ac e o he dipole-dipole in e ac ion, which decays as he
in e se cubic powe o he dis ance, i is necessa y o include mo e han
one nea es neighbo o ha e a ai h ul quan i a i e desc ip ion o dipola
sys ems. In ac , longe - ange in e ac ions end o allow o and s abilize
mo e no el phases.
In his hesis we ﬁ s s udy BH models wi h dipola in e ac ions, going
beyond he g ound s a e sea ch. We conside a wo-dimensional (2D) la ice
whe e he dipoles a e pola ized pe pendicula ly o he 2D plane, esul ing in
3
an iso opic epulsi e in e ac ion. We use he mean-ﬁeld app oxima ions and
a Gu zwille Ansa z which a e qui e accu a e and sui able o desc ibe his
sys em. We ﬁnd ha dipola bosonic gas in 2D la ices exhibi s a mul i ude
o insula ing me as able s a es, o en compe ing wi h he g ound s a e, simi-
la ly o a diso de ed sys em. We s udy in de ail he a e o hese me as able
s a es: how can hey be p epa ed on demand, how hey can be de ec ed,
wha is hei li e ime due o unneling, and wha is hei ole in a ious cool-
ing schemes. Mo eo e , we ﬁnd ha he g ound s a e is cha ac e ized by
insula ing checke boa d-like s a es wi h ac ional ﬁlling ac o s ν(a e age
numbe o pa icles pe si e) ha depend on he cu -oﬀ used o he in e -
ac ion ange. We conﬁ m his p edic ion by s udying he same sys em wi h
Quan um Mon e Ca lo me hods ( he wo m algo i hm). In his case no cu -oﬀ
o he dipola in e ac ion is used, and we ﬁnd e idence o a De il’ s s ai case
in he g ound s a e, i.e. insula ing phases which appea a all a ional νo
he unde lying la ice. We also ﬁnd egions o pa ame e s whe e he g ound
s a e is a supe solid, ob ained by doping he solids ei he wi h pa icles o
acancies. Recen ly [13], a comple e de il’ s s ai case has been p edic ed in
he phase diag am o a one-dimensional dipola Bose gas.
In his wo k, we also in es iga e how he p e ious scena io changes by
conside ing a mul i-laye s uc u e. We ocus on he simples si ua ion com-
posed o wo 2D laye s in which he dipoles a e pola ized pe pendicula ly o
he planes; he dipola in e ac ion is hen epulsi e o pa icles laying on
he same plane, while i is a ac i e o pa icles a he same la ice si e on
diﬀe en laye s. Ins ead we conside in e -laye unneling o be supp essed,
which makes he sys em analogous o a bosonic mix u e in a 2D la ice. Ou
calcula ions show ha pa icles pai in o composi es, and demons a e he
exis ence o he no el Pai Supe Solid (PSS) quan um phase.
Cu en ly we a e s udying a 2D la ice whe e he dipoles a e ee o poin
in bo h di ec ions pe pendicula ly o he plane, which esul s in a nea es
neighbo epulsi e (a ac i e) in e ac ion o aligned (an i-aligned) dipoles.
We ﬁnd egions o pa ame e s whe e he g ound s a e is e omagne ic o an i-
e omagne ic, and ﬁnd e idences o he exis ence o a Coun e ﬂow Supe
Solid (CSS) quan um phase.
Ou p edic ions ha e di ec expe imen al consequences, and we hope ha
hey will be soon checked in expe imen s wi h ul acold dipola a omic and
molecula gases. This hesis is based on he ollowing publica ions:
•C. Meno i, C. T e zge , and M. Lewens ein, Me as able S a es o a Gas
o Dipola Bosons in a 2D Op ical La ice. Physical Re iew Le e s,
98, 235301 (2007).
•C. T e zge , C. Meno i, and M. Lewens ein, Ul acold dipola gas in
4
an op ical la ice: The a e o me as able s a es. Physical Re iew A,
78, 043604, (2008).
•C. T e zge , C. Meno i, and M. Lewens ein, Pai -Supe solid Phase in
a Bilaye Sys em o Dipola La ice Bosons. Physical Re iew Le e s,
103, 035304, (2009).
•B. Capog osso-Sansone, C. T e zge , M. Lewens ein, P. Zolle , and G.
Pupillo, Quan um Phases o Cold Pola Molecules in 2D Op ical La -
ices. a Xi :0906.2009. Accep ed o Physical Re iew Le e s publica-
ion.
•C. T e zge , M. Alloing, C. Meno i, F. Dubin, and M. Lewens ein,
Coun e low Supe solid o an i-pola ized dipola Bosons in a 2D op ical
la ice. In p epa a ion.
5

6
Con en s
I Theo y and me hods 10
1 Dipola Bose gas in op ical la ices 13
1.1 Op ical la ices . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2 Theo y o dilu e Bose gases . . . . . . . . . . . . . . . . . . . 15
1.2.1 The G oss-Pi ae skii equa ion . . . . . . . . . . . . . . 17
1.2.2 Bose-Hubba d model . . . . . . . . . . . . . . . . . . . 18
1.3 Dipola Bose gas . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.3.1 P ope ies o he dipole-dipole in e ac ion . . . . . . . 22
1.3.2 Pola ized dipoles in aniso opic ha monic aps . . . . 23
1.3.3 Mean-ﬁeld dipola in e ac ion in a sphe ical ap . . . . 26
1.3.4 Ex ended Bose-Hubba d model . . . . . . . . . . . . . 28
2 Hubba d models: heo e ical me hods 31
2.1 Supe ﬂuid−Mo insula o quan um phase ansi ion in he
Bose-Hubba d model . . . . . . . . . . . . . . . . . . . . . . . 31
2.2 The Gu zwille mean-ﬁeld app oach . . . . . . . . . . . . . . . 33
2.2.1 Dynamical Gu zwille app oach . . . . . . . . . . . . . 34
2.2.2 Pe u ba i e mean-ﬁeld app oach . . . . . . . . . . . . 37
2.2.3 Pe u ba i e mean-ﬁeld s. dynamical Gu zwille ap-
p oach........................... 39
II Me as able s a es 41
3 Dipola Bosons in a 2D op ical la ice 43
3.1 Themodel............................. 43
3.2 Me as abili y . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.3 The li e ime . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.3.1 Pa ame iza ion . . . . . . . . . . . . . . . . . . . . . . 51
3.3.2 Ac ion and unneling ime . . . . . . . . . . . . . . . . 52
3.4 P epa a ion, manipula ion and de ec ion . . . . . . . . . . . . 54
7
3.4.1 T ans e p ocess . . . . . . . . . . . . . . . . . . . . . . 56
3.5 Ha monic conﬁnemen . . . . . . . . . . . . . . . . . . . . . . 59
4 Conclusions 61
III Mul iple laye s and mix u es 62
5 Dipola Bosons in a bilaye op ical la ice 65
5.1 Themodel............................. 65
5.1.1 G ound s a e and single-pa icle single-hole exci a ions 67
5.2 Low-ene gy subspace and eﬀec i e Hamil onian . . . . . . . . 70
5.2.1 G ound s a e insula ing phases and wo-pa icle wo-
hole exci a ions . . . . . . . . . . . . . . . . . . . . . . 72
5.3 Gu zwille mean-ﬁeld app oach and alidi y o he low ene gy
subspace.............................. 73
6 Coun e low Supe solid o an i-pola ized dipola Bosons in a
2D op ical la ice 78
6.1 In oduc ion . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
6.2 Hamil onian o he sys em . . . . . . . . . . . . . . . . . . . . 78
6.2.1 Filling ac o and imbalance . . . . . . . . . . . . . . . 79
6.2.2 Low-ene gy subspace and eﬀec i e Hamil onian . . . . 82
6.3 Mean-ﬁeld............................. 85
6.3.1 Insula ing lobes . . . . . . . . . . . . . . . . . . . . . . 86
6.3.2 Coun e ﬂow supe ﬂuid-supe solid . . . . . . . . . . . . 87
7 Conclusions 91
IV Quan um Mon e Ca lo 92
8 Pa h In eg al Mon e Ca lo and he Wo m algo i hm 95
8.1 Pa h In eg al Mon e Ca lo . . . . . . . . . . . . . . . . . . . . 95
8.1.1 Pa h In eg al Mon e Ca lo and he 2D ex ended Bose-
Hubba d model . . . . . . . . . . . . . . . . . . . . . . 97
8.2 The Wo m algo i hm . . . . . . . . . . . . . . . . . . . . . . . 100
8.2.1 Upda ing p ocedu es . . . . . . . . . . . . . . . . . . . 100
8.2.2 Ad an ages o he Wo m algo i hm . . . . . . . . . . . 104
8
9 Quan um Mon e Ca lo s udies o dipola gases 106
9.1 Ze o momen um G een unc ion and he pa icle-hole exci a ions107
9.2 Incomp essible and supe solid phases . . . . . . . . . . . . . . 107
9.2.1 Homogeneous case . . . . . . . . . . . . . . . . . . . . 108
9.2.2 Fini e empe a u e . . . . . . . . . . . . . . . . . . . . 110
9.2.3 Ha monic conﬁnemen . . . . . . . . . . . . . . . . . . 110
10 Conclusions 113
A Spec um o exci a ions 116
Bibliog aphy 119
9
empe a u es he gas is a Bose Eins ein condensa e (BEC). The ype o
pa icles we conside he e can be a oms o molecules.
Fo a dilu e gas, he in e pa icle sepa a ion ( ypically o he o de o 102
nm o alkali a oms [17]) is an o de o magni ude la ge han he leng h
scales associa ed wi h he a om-a om in e ac ion. In o he wo ds, a dilu e
gas o densi y nis a e y a eﬁed gas in which he ”spa ial ex ension” o an
a om is much smalle han he a e age olume pe pa icle n−1. Because o
his condi ion, he wo-body in e ac ion domina es he physics while h ee-
body o mo e a e e y unlikely and essen ially no impo an . The wo-body
in e a omic po en ial V( ) depends on he ype o pa icles one conside ,
he ela i e dis ance be ween he a oms = 1− 2and on hei in e nal
s a es. Fo alkali a oms, he po en ial is s ongly epulsi e o small a omic
sepa a ions while o la ge a omic dis ances i is domina ed by he an de
Waals a ac ions ha decay as −C6/ 6, whe e he coeﬃcien C6depends on
he a omic species.
He e we will conside only elas ic sca e ing, whe e he in e nal s a es o
he wo a oms do no change in he collision p ocess. I he empe a u e o
he gas is e y low, i.e. T→0, hen he kine ic ene gy o he pa icles is e y
small compa ed o he cen i ugal ba ie and only s-wa e sca e ing akes
place. The e o e, he only impo an pa ame e is he sca e ing leng h gi en
by
as=m
4π~2Zd3 V ( ),(1.9)
wi h mbeing he mass o he a oms. This quan i y has he dimensions o
a leng h and has he physical in e p e a ion o he adius he a oms would
ha e i hey we e conside ed o be pe ec billia d balls. The condi ion o
he dilu eness o he gas hen eads
na3
s≪1 (1.10)
whe e nis he densi y o he gas and na3
sis called he gas pa ame e . One
can in e he exp ession (1.9) and hink o an eﬀec i e in e ac ion be ween
he wo pa icles p opo ional o he sca e ing leng h, and gi en by
Ve ( ) = g δ(3)( ),(1.11)
whe e gis deﬁned as
g=4π~2as
m,(1.12)
and δis he Di ac del a unc ion so ha he pa icles a e conside ed o be
poin -like. No e also ha since he eﬀec i e in e ac ion depends only on he
sca e ing leng h, i is epulsi e (a ac i e) o posi i e (nega i e) a, and i
16

can be dynamically modiﬁed o example in alkali a oms jus by a ying an
ex e nal magne ic ﬁeld nea a Feshbach esonance.
1.2.1 The G oss-Pi ae skii equa ion
The quan um s a e o a gas o Npa icles is desc ibed by he many-body
wa e unc ion Ψ( 1, 2, ..., N), and he ime e olu ion o he sys em is de e -
mined by he Sch ¨
odinge equa ion. In a BEC, one can desc ibe he dynamics
o he condensa e jus h ough he G oss-Pi ae skii (GP) equa ion [16, 17]
gi en by
i~∂
∂ Ψ0( , ) = −~2∇2
2m+Vex ( ) + g|Ψ0( , )|2Ψ0( , ),(1.13)
whe e Ψ0( , ) is he BEC wa e unc ion, also called he o de pa ame e .
The in e ac ion be ween pa icles has been aken in o accoun in a mean-
ﬁeld app oxima ion by he e m g|Ψ0( , )|2, and gis gi en in Eq. 1.12.
The wo-body eﬀec i e po en ial is gi en by Eq. (1.11), Vex ( ) is an ex e nal
apping po en ial, and he o de pa ame e is no malized o he o al numbe
o pa icles, i.e. N=Rd3 |Ψ0( , )|2. Equa ion (1.13) was independen ly
de i ed by G oss and Pi ae skii in 1961, i is one o he main heo e ical ools
o in es iga ing dilu e weakly in e ac ing Bose gases a low empe a u es,
and i has he ypical o m o a mean ﬁeld equa ion whe e he o de pa ame e
mus be calcula ed in a sel -consis en way.
The GP equa ion has p o en o be a e y use ul ool o desc ibe he
physics o weakly in e ac ing Bose-Eins ein a omic condensa es in he ea ly
ages o his ﬁeld. Wi h his o malisms, and i s ex ension o include small
ﬂuc ua ions gi en by Bogoliubo heo y, one can desc ibe accu a ely, among
o he s, he collec i e exci a ions o he sys ems, he esponse o o a ions
including he o ma ion o o ices, he p opaga ion o sound, he p esence
o dynamical ins abili ies. Gene ally speaking, he GP ea men is well
sui ed in he egime o ull cohe ence, when a single mac oscopically occupied
ma e wa e co ec ly desc ibes he sys em. A he end o he ’90, ew yea s
a e he c ea ion o he ﬁ s alkali BECs in he lab, he need o ”going
beyond GP” s a ed o be e y s ongly el , due o he heo e ical in e es
and expe imen al possibili y o going in o he s ongly co ela ed egime.
In ac , he p esence o s ong in e ac ions, s ong o a ions and/o special
apping po en ials can limi he alidi y o he GP equa ion. Fo ins ance
a s ong conﬁnemen in one o wo dimensions can educe he sys em o an
eﬀec i ely 2D o 1D one. A s ong o a ion combined wi h in e ac ions can
lead o quan um Hall physics. Also he p esence o a deep op ical la ices,
17
when he combined eﬀec o in e ac ions and apping po en ial leads o a
” agmen a ion” o he condensa e, equi es mo e sophis ica ed desc ip ions.
In his hesis we a e in e es ed in desc ibing he physics o Bosons apped
in a pe iodic op ical po en ial (Vop ) and e en ually also conﬁned in a mag-
ne ic ha monic ap (Vho), he o al ex e nal ﬁeld being gi en by he sum
Vex ( ) = Vop ( ) + Vho( ) = X
i=x,y,z
V0,i cos2(ki i) + 1
2mX
i=x,y,z
ω2
i 2
i,(1.14)
whe e (V0,x, V0,y, V0,z) is he dep h o he op ical la ice in he h ee spa ial
di ec ions and (ωx, ωy, ωz) he equencies o he ha monic ap. In o de o
desc ibe he physics o Bosons apped in he po en ial (1.14), we need o
”go beyond” he GP equa ion, and we will de o e he ollowing sec ions o
his pu pose.
1.2.2 Bose-Hubba d model
The s a ing poin o ou discussion is Hamil onian (1.15), w i en in he
second quan iza ion o malism in e ms o he c ea ion and annihila ion op-
e a o s o Bosons, ˆ
ψ†( ) and ˆ
ψ( ) espec i ely, and gi en by he exp ession
ˆ
H=Zd3 ˆ
ψ†( )−~2∇2
2m+Vex ( ) + g
2ˆ
ψ†( )ˆ
ψ( )−µˆ
ψ( ),(1.15)
whe e he ﬁ s e m in squa e b acke s is he kine ic ene gy, Vex ( ) =
Vop ( ) + Vho( ) is he ex e nal apping po en ial (1.14) and we ha e used
he simpliﬁed con ac in e ac ion (1.11). We wo k in he g and canonical
ensemble such ha he chemical po en ial µﬁxes he o al numbe o pa i-
cles. Addi ionally, we assume he ha monic conﬁnemen o change on a scale
la ge han he one o he op ical la ice, such ha we can conside he eﬀec
o he magne ic apping o be cons an o e a single si e o he la ice.
In his o malism, he ﬁeld ope a o s can be w i en in he basis o single-
pa icle wa e unc ions {Φn( )}n, whe e nis a comple e se o single pa icle
quan um numbe s
ˆ
ψ( ) = X
n
Φn( )ˆan
ˆ
ψ†( ) = X
n
Φ∗
n( )ˆa†
n,(1.16)
wi h ˆa†
nand ˆanbeing he c ea ion and annihila ion ope a o s on he Fock
s a e o he mode n, i.e. ˆa†
n|ni=√n+ 1|n+ 1iand ˆan|ni=√n|n−1i.
18
Also, he ﬁeld ope a o s sa is y he usual commu a ion ela ions o Bosons
[ˆ
ψ( ),ˆ
ψ†( ′)] = ∞
X
n=0
Φn( )Φ∗
n( ′) = δ3( − ′),
[ˆ
ψ( ),ˆ
ψ( ′)] = [ ˆ
ψ†( ),ˆ
ψ†( ′)] = 0.
(1.17)
I is well known [20], ha he spec um o a single pa icle in a pe iodic
po en ial is cha ac e ized by bands o allowed ene gies and ene gy gaps, and
he single pa icle wa e unc ions a e desc ibed by Bloch unc ions Φαk( )
wi h band index αand quasi-momen um ~k. Al e na i ely, he e exis s a
complemen a y single-pa icle basis gi en by he Wannie unc ions [20, 21]
wα( −Ri), whe e Riis a la ice ec o poin ing a si e iand wα( ) a e
deﬁned as he Fou ie ans o m o Bloch unc ions
wα( ) = 1
√NSX
k
e−ik· Φαk( ),(1.18)
whe e NS, is he o al numbe o si es in he la ice. The Wannie unc ions
o m a comple e o hono mal se , so one may w i e he ﬁeld ope a o s (1.16)
as
ˆ
ψ( ) = X
αk
Φαk( )ˆaαk=X
α,i
wα( −Ri)ˆaα,i
ˆ
ψ†( ) = X
αk
Φ∗
αk( )ˆa†
αk=X
α,i
w∗
α( −Ri)ˆa†
α,i.(1.19)
Wannie unc ions a e use ul in he case o deep op ical la ices whe e igh
binding app oxima ion apply. The big ad an age o using Wannie unc ions
wα( −Ri) is ha hey a e localized and cen e ed a ound he la ice si e
poin ed by Ri.
I he empe a u e o he sys em is low enough, and he in e ac ions
be ween he pa icles is no suﬃcien o induce ansi ions be ween he bands,
one may es ic only o he ﬁ s Bloch band because he pa icles ha e
insuﬃcien ene gy o o e come he gap ha sepa a es he ﬁ s band om
he o he s. This amoun s o keep in (1.19) only he lowes o he αindices,
which we omi o simplici y o no a ion and he e o e he Hamil onian (1.15)
becomes
ˆ
H=−X
i,j
Jij ˆa†
iˆaj+X
i,j,k,l
Ui,j,k,l
2ˆa†
iˆa†
jˆakˆal−X
i,j
µi,j ˆa†
iˆaj.(1.20)
19
The quan i ies in he sums a e gi en by
Jij =−Zd3 w∗( −Ri)−~2∇2
2m+Vop ( )w( −Rj) (1.21)
Ui,j,k,l =gZd3 w∗( −Ri)w∗( −Rj)w( −Rk)w( −Rl) (1.22)
µi,j =Zd3xw∗( −Ri) [µ−Vho( )] w( −Rj).(1.23)
The Wannie unc ions a e localized on he la ice si es, he deepe he la ice
he mo e localized hey a e. Fo a suﬃcien ly deep op ical po en ial, hen
in Eq. (1.22) and (1.23) he dominan con ibu ions a e gi en by Ui,i,i,i and
µi,i. Fo he kine ic pa (1.21), he e is a cons an con ibu ion gi en by
Ji,i and due o he p esence o he de i a i e in he in eg a ion, he e is also
a posi i e ma ix elemen o nea es neighbo ing si es Ji,j >0. The wo
si ua ions a e quali a i ely shown in Fig. (1.2) whe e we ha e app oxima ed
he Wannie unc ions wi h wo Gaussians espec i ely localized a si e iand
jo he la ice. Howe e , we s ess ha he pic u e p o ided by Gaussian
unc ions is only quali a i e. In ac , in o de o be quan i a i ely co ec ,
one needs o calcula e he p ope ma ix elemen s wi h Wannie unc ions.
ij
(a)
ij
(b)
Figu e 1.2: (a) Two Gaussians localized on neighbo ing si es iand jo
an op ical la ice ha ing negligible o e lap. (b) The ﬁ s de i a i e o he
Gaussian unc ions ins ead, show a nega i e o e lap in he egion indica ed
by he a ow, which leads o a posi i e ma ix elemen Ji,j >0.
Wi h he abo e conside a ions, we can now w i e he celeb a ed Bose-Hubba d
Hamil onian in he o m
ˆ
HBH =−JX
hiji
ˆa†
iˆaj+U
2X
i
ˆni(ˆni−1) −X
i
µiˆni,(1.24)
whe e hijiindica es sum o e nea es neighbo s, he unneling coeﬃcien
J=Ji,j =Jj,i o he mi ici y, he on-si e in e ac ion U=gRd3 |w( )|4,
20
ˆni= ˆa†
iˆaiis he numbe ope a o a si e i, and we ha e neglec ed Ji,i since
i gi es a cons an con ibu ion o each si e. The ha monic conﬁnemen ,
since i is assumed o be cons an ac oss one la ice si e, has been aken in o
accoun in he chemical po en ial as
µi=µ−1
2m~ω 2·(Ri−R0)2,(1.25)
whe e R0is he cen e o he ha monic ap wi h equencies gi en by ~ω =
(ωx, ωy, ωz) in he h ee di ec ions. The second e m on he igh hand side
o Eq. (1.25) is p ac ically a chemical po en ial ha diﬀe s om si e o si e
and i is o en called he local chemical po en ial.
Fo a one dimensional op ical la ice Vop (x) = V0sin2(kx) wi h wa e ec-
o k= 2π/λ, Fig. 1.3 shows bo h he on-si e in e ac ion U(solid line) and
he unneling coeﬃcien J(dashed line) as a unc ion o he op ical la ice
dep h V0, whe e all he quan i ies a e measu ed in e ms o he ecoil ene gy
ER=~2k2/2m, ha is he ene gy acqui ed by he a om a e abso bing a
pho on wi h momen um ~k. The la ice pa ame e s Uand Jwe e calcula ed
nume ically in e.g. [22] o diﬀe en alues o V0. F om Fig. 1.3 (b), i is clea
ha i is possible o change he unneling coeﬃcien Jo e a wide ange,
going om a si ua ion o p ac ically isola ed la ice si es a V0= 25ERup
o a egime in which pa icles can unnel om si e o si e a V0= 5ER, only
by changing he op ical po en ial dep h by a ew ens o ecoil ene gies, and
lea ing he on-si e in e ac ion Up ac ically unchanged.
Figu e 1.3: (a) Schema ic ep esen a ion o a 1D op ical la ice; (b) scaled
on-si e U(solid line) and unneling coeﬃcien J(dashed line) dependence
on he op ical po en ial dep h V0. The on-si e in e ac ion is mul iplied by
a/as(≫1), whe e a=λ/2 is he la ice pe iod and asis he s-wa e sca e ing
leng h o a oms o equal mass m. Figu e om [22].
21

1.3 Dipola Bose gas
1.3.1 P ope ies o he dipole-dipole in e ac ion
Two pa icles 1 and 2 in a h ee dimensional space, a ela i e dis ance
and wi h dipole momen s along he uni ec o s e1and e2as in Fig. 1.4
(a), in e ac h ough he dipole-dipole in e ac ion such ha hei in e ac ion
ene gy is gi en by
Udd( ) = Cdd
4π
(e1·e2) 2−3(e1· )(e2· )
5,(1.26)
whe e =| |, and Udd( ) = Udd(− ). The dipola coupling cons an Cdd is
diﬀe en o pa icles ha ing a pe manen magne ic dipole momen µ, and o
pa icles ha ing a pe manen elec ic dipole momen d, and is espec i ely
gi en by
Cdd =µ0µ2magne ic
d2/ε0elec ic, (1.27)
whe e µ0is he acuum pe meabili y, and ε0is he acuum pe mi i i y.
(a)
(c)
(b)
(d)
e1
e2
e1
e2
θ
Figu e 1.4: (a) Two dipoles, 1 and 2, di ec ed along uni ec o s e1and e2
and sepa a ed by a dis ance . (b) Pola ized dipoles, o which he in e -
ac ion depends on he angle θbe ween he di ec ion o he dipoles and he
in e pa icle sepa a ion . This esul s in a epulsi e in e ac ion o θ=π/2
(c), and a ac i e o θ= 0(π) (d). Figu e om [40].
22
The dipole-dipole in e ac ion (1.26) has a long- ange cha ac e ; his is
because a la ge dis ances i decays as Udd ∼1/ 3, con a y o he ypical
an de Waals po en ial ha beha es like U dW ∼ −1/ 6. Also, om (1.26)
i is easy o see he aniso opic p ope y o his in e ac ion; o pola ized
a oms, i.e. all dipoles poin ing in he same di ec ion (say z), he in e ac ion
educes o
Udd( ) = Cdd
4π
1−3 cos2θ
3,(1.28)
whe e θis he angle be ween he dipole and he ela i e dis ance o he
pa icles, as in Fig 1.4 (b). The in e ac ion is epulsi e o θ=π/2 as he
example o Fig 1.4 (c), and a ac i e o θ= 0 as shown in Fig 1.4 (d). The
si ua ion is e e sed o an i-pa allel dipoles, whe e a minus sign appea s in
on o Eq. (1.28), and he e o e he in e ac ion is a ac i e o θ=π/2
while θ= 0 gi es ise o epulsion.
The sca e ing p ope ies o ul acold a oms, in he simple case o iso opic
an de Waals in e ac ions, a e en i ely desc ibed by he s-wa e sca e ing
leng h and he po en ial can be eplaced by he eﬀec i e con ac in e ac ion
(1.11). In he p esence o a dipola in e ac ion as (1.26), because o i s long
ange (decay as 1/ 3) and aniso opic cha ac e (s ong dependence on he
ela i e angles be ween he dipoles), all pa ial wa es con ibu e o he sca -
e ing p oblem and also pa ial wa es wi h diﬀe en angula momen a couple
wi h each o he . While o Fe mions, eplacing he eal po en ial (1.26) wi h
an eﬀec i e dipola in e ac ion as (1.11) is easonable [23], o Bosons his i
is no ob ious, and in ecen yea s i has been he subjec o in ensi e s udies
[24, 25, 26, 27]. In he p esence o an op ical la ice, i has been ecen ly
a gued [41] ha in a 1D geome y, eplacing he eal dipola po en ial wi h
an eﬀec i e in e ac ion as (1.11) is easonable as long as he op ical la ice is
shallow enough. Howe e , in he mos gene al case i is necessa y o accoun
o he ull exp ession o he dipole-dipole in e ac ion po en ial (1.26).
1.3.2 Pola ized dipoles in aniso opic ha monic aps
We now mo e o he desc ip ion o a BEC o pola ized dipoles, poin ing along
he zaxis. Fo pola ized dipola BECs, due o he aniso opy o he dipola
in e ac ions, he geome y o he apping po en ial plays a undamen al ole,
ﬁ s in de e mining he spa ial dis ibu ion o he densi y, and second in he
s abili y o he gas.
Quali a i ely, he e a e wo ex eme scena ios depending on he shape
o he conﬁning po en ial, shown in Fig. 1.5: (i) o a ciga -shaped ap
elonga ed along he zaxis, i.e. wi h an aspec a io be ween he axial ωz
and adial equencies ωρ=ωx=ωygi en by λ=ωz/ωρ≪1, he densi y is
23
mainly dis ibu ed along he pola iza ion axis and he eﬀec o dipole-dipole
in e ac ion is mos ly a ac i e, which migh lead o an ins abili y o he
gas e en in he p esence o a weak epulsi e con ac in e ac ion; (ii), o a
pancake-shaped ap, which is s ongly conﬁning along he zaxis, i.e. λ≫1,
he dipola in e ac ion is mos ly epulsi e and he BEC is always s able o
epulsi e con ac in e ac ions and migh be s able e en o a ac i e con ac
in e ac ions. In an in e media e si ua ion in which he conﬁning po en ial is
pe ec ly sphe ical, he densi y dis ibu ion is hen iso opic and he dipole-
dipole in e ac ion a e ages ou o ze o, which leads o a s able BEC o
epulsi e con ac in e ac ions. One can swi ch be ween one o he o he
scena io, jus by adjus ing he equency o he conﬁning po en ial along he
zaxis wi h espec o he axial xand y, and he e o e i is na u al o expec
ha o any gi en λ he e is a c i ical alue o he sca e ing leng h ac i
below which he BEC is uns able [43].
(a) (b)
Figu e 1.5: Pola ized dipoles in aniso opic ha monic po en ials. (a) in a ciga
shaped ap elonga ed in he di ec ion o pola iza ion, he esul ing dipola
in e ac ion is a ac i e, and (b) in a pancake ap wi h a s ong conﬁnemen
in he di ec ion o pola iza ion, he dipola in e ac ions a e epulsi e. Figu e
aken om [40].
One can quan i a i ely desc ibe he abo e scena ios s a ing om he Hamil-
onian o he sys em, which in he p esence o he dipole-dipole in e ac ion
(1.28) eads
ˆ
H=Zd3 ˆ
ψ†( )−~2∇2
2m+Vex ( ) + g
2ˆ
ψ†( )ˆ
ψ( )−µˆ
ψ( )
+1
2ZZ d3 1d3 2ˆ
ψ†( 1)ˆ
ψ†( 2)Udd( 1− 2)ˆ
ψ( 1)ˆ
ψ( 2).(1.29)
Wi h he same app oxima ions used o de i e he G oss-Pi ae skii equa ion,
one can w i e he Boson ﬁeld ope a o ˆ
ψ( ) = Ψ0( ) + δˆ
ψ( ) as a sum o a
classical ﬁeld Ψ0( ), he condensa e wa e unc ion, plus he non condensa e
componen δˆ
ψ( ) [16]. By neglec ing he ﬂuc ua ions δˆ
ψ( ), one can calcula e
24
he ene gy o he BEC gi en by
Eψ=Zh−~2
2m|∇ψ( )|2+Vex ( )|ψ( )|2+g
2|ψ( )|4
+1
2|ψ( )|2ZUdd( − ′)|ψ( ′)|2d3 ′id3 . (1.30)
Wi hin a a ia ional Ansa z, we assume he condensa e wa e unc ion o be
a Gaussian o axial wid h σzand adial wid h σx=σy=σρ, no malized o
he o al numbe o pa icles N, namely
Ψ0(z, ρ) = sN
π3/2σ2
ρσza3
ho
exp −1
2a2
ho ρ2
σ2
ρ
+z2
σ2
z,(1.31)
whe e aho =p~/(m¯ω) is he ha monic oscilla o leng h wi h a e age ap
equency ¯ω= (ω2
ρωz)1/3. The e o e, inse ing he Ansa z (1.31) in o he
ene gy unc ional Eq. (1.30), a e in eg a ion we ﬁnd he ene gy o he BEC
o be a unc ion o he wid hs o he Gaussians, namely
E0(σz, σρ) = Ekin +E ap +Econ ac +Edd,(1.32)
wi h he kine ic ene gy
Ekin =N~¯ω
42
σ2
ρ
+1
σ2
z,(1.33)
he po en ial ene gy due o he ap
E ap =N~¯ω
4λ2/32σ2
ρ+λ2σ2
z,(1.34)
he con ac in e ac ion ene gy gi en by
Econ ac =~¯ω
√2πaho
1
σ2
ρσz
as,(1.35)
and he con ibu ion coming om he dipola e m
Edd =−~¯ωadd
√2πaho
1
σ2
ρσz
(κ),(1.36)
whe e we ha e in oduced he dipola leng h add =Cddm
12π~2, wi h Cdd gi en in
Eq. (1.27), which measu es he absolu e s eng h o he dipola in e ac ion,
25
MI(¯n). Fo small alues o J/U he MI(¯n) phase pe sis s in a closed and
ﬁni e a ea o he J s. µplane [1], which is called he Mo lobe o MI(¯n).
The la ge J/U alue o he lobe is called he ip o he lobe o also c i ical
poin (J/U)c. The c i ical poin changes wi h he dimensionali y and geom-
e y o he sys em.In Fig. 2.1 we plo he ﬁ s ¯n= 0,1,2,3 insula ing lobes,
calcula ed o an inﬁni e op ical la ice wi hin he mean-ﬁeld app oxima ion,
which will be discussed in Sec. 2.2.2. The hick black lines enclose he lobes
and ma k he bounda ies be ween he MI and SF phases. Ou side he in-
sula ing lobes, he phase is SF. The colo ed lines o Fig. 2.1(a) indica e a
con ou plo o cons an ac ional densi y, while he hick black lines depa -
ing om he ip o he lobes and ex ending in o he SF egion, co espond
o an in ege alue ¯no he densi y.
0 0.08 0.16 0.24
0
1
2
3
zJ/U
µ/U
MI(0)
MI(1)
MI(2)
MI(3)
(a)
0 0.04 0.16 0.24
0
1
2
3
zJ/U
MI(0)
MI(1)
MI(2)
MI(3)
(b)
Figu e 2.1: Mean-ﬁeld phase diag am o he Bose Hubba d Hamil onian. (a)
Con ou plo o he densi y pe si e; (b) con ou plo o he o de pa ame e
(see Sec. 2.2). MI(¯n) indica e a Mo insula ing phase wi h ﬁxed ¯na oms
pe si e.
We will de i e he mean-ﬁeld Mo insula ing lobes o Fig. 2.1 in a mo e
igo ous way in Sec. 2.2.2, bu o he momen we jus lis he c i ical poin s
(J/U)c, o he ¯n= 1 lobe, ha ha e been es ima ed wi h diﬀe en me h-
ods and o diﬀe en dimensions o he la ice. In one dimension, he c i -
ical poin has been es ima ed o be (J/U)c≃0.29 [48] using Densi y Ma-
32

ix Reno maliza ion G oup calcula ions (DMRG). In wo dimensions, wi h
quan um Mon e Ca lo calcula ions, he c i ical poin has been es ima ed o
be (J/U)c≃0.061 [49], while in he h ee dimensional model he loca ion o
he c i ical poin has been es ima ed wi h pe u ba i e expansions o be a
(J/U)c≃0.034 [50]. In he nex sec ion we will de i e he mean-ﬁeld lobes.
2.2 The Gu zwille mean- ield app oach
The Gu zwille mean-ﬁeld app oach is an app oxima ion o he many-body
wa e unc ion o Hubba d- ype Hamil onians and is gi en by
|Ψi=Y
i
nmax
X
n=0
(i)
n|nii,(2.3)
whe e |nii ep esen s he Fock s a e o na oms occupying he si e i,nmax is
a cu oﬀ in he maximum numbe o a oms pe si e, and (i)
nis he p oba-
bili y ampli ude o ha ing he si e ioccupied by na oms. The p obabili y
ampli udes a e no malized o uni y Pn| (i)
n|2= 1. The wa e unc ion (2.3)
has been ex ensi ely used in he li e a u e [22, 51, 52], and is mo i a ed by
i s physical p edic ions; in ac , he e exis s a c i ical alue (J/U)m in a
gi en ange o µ, below which he g ound s a e p edic ed by he Gu zwille
Ansa z is a p oduc o single Fock s a es (i)
n=δn,¯nwi h exac ly ¯npa icles
pe si e, as (2.2). Mo eo e , o J/U > (J/U)m he Gu wille Ansa z p e-
dic s a supe ﬂuid g ound s a e wi h ﬂuc ua ing on-si e pa icle numbe . The
Gu zwille c i ical poin , o ¯n= 1, is ound o be (J/U)m = 1/5.8z[53],
whe e z=Phjii1 is he numbe o nea es neighbo connec ions a each si e
o he la ice. In able (2.2) we show he compa ison o he c i ical poin s
p edic ed by he Gu zwille Ansa z o diﬀe en dimensions o he sys em,
wi h he mo e p ecise ones discussed in Sec. (2.1). F om he compa ison,
Dz(J/U)m (J/U)c
1 2 0.0862 0.29
2 4 0.0431 0.061
3 6 0.0287 0.034
Table 2.1: Compa ison o he Gu zwille c i ical poin s (J/U)m wi h he
mo e p ecise, up o now, c i ical poin s (J/U)c, o diﬀe en dimensions D o
he sys em.
one can deduce ha he Guzwille is unsa is ac o y o 1D sys ems (z= 2)
33
while i is sa is ac o y o a 3D one (z= 6). Also, in he limi o J/U → ∞
he diﬀe ence be ween he Gu zwille p edic ions and he exac esul s a e
negligible [53]. Summa izing, he Gu zwille p edic ions a e exac in he wo
limi ing cases o J/U →0 and J/U → ∞, while o in e media e cases he
pe o mance o he Gu zwille app oach s ongly depends on he dimension
o he la ice, since i does no co ec ly accoun o he quan um ﬂuc ua ions
a he phase ansi ion.
An impo an quan i y is he so called o de pa ame e , which is he
expec a ion alue o he Bosonic annihila ion ope a o a he i- h si e o he
la ice, namely ϕi=hΨ|ˆai|Ψi, and by using he Gu zwille wa e unc ion
(2.3) one ge s
ϕi=X
n
√n+ 1 ∗(i)
n (i)
n+1.(2.4)
The o de pa ame e ϕidesc ibes he phase o he sys em a he si e io
he la ice: i is exac ly ze o in he Mo phase ϕi= 0, while i assumes
a non ze o alue ϕi6= 0 in he supe ﬂuid phase. In he uni o m sys em,
he la ice is ansla ionally in a ian and he e o e all si es a e sel -simila ,
which means ha a single o de pa ame e de e mines he phase o he whole
sys em. In Fig. 2.1(b) we plo he absolu e alue o he o de pa ame e ϕi
o such a sys em, in he J/U s. µ/U plane. The colo ed lines ou side he
insula ing lobes co espond o a con ou plo o cons an non-ze o alue o
ϕi ypical o he SF phase. Ins ead, in a non-uni o m sys em as i is in he
p esence o an ex e nal conﬁning ha monic po en ial, diﬀe en phases can
coexis . As an example, in Fig. 2.2 we plo he densi y o he g ound s a e
(a) along wi h he o de pa ame e a each si e (b), o a 2D la ice in he
p esence o a conﬁning ha monic po en ial. No ice ha he MI phase a he
cen e o he ha monic ap (x0, y0) is su ounded by a ing o SF phase, in
a wedding-cake like s uc u e, as ﬁ s discussed in [22]. The g ound s a e
o Fig. 2.2 was ob ained wi hin he mean-ﬁeld app oxima ion h ough he
imagina y- ime e olu ion echnique, ha will be discussed in Sec. 2.2.1.
In he p esence o dipola in e ac ions, as we shall see la e on, i is also
necessa y o accoun o non-uni o m quan um phases, because e en in he
uni o m sys em he p esence o dipola in e ac ions may lead o spon aneous
symme y b eaking o ansla ional in a iance on a scale la ge han he
la ice cons an .
2.2.1 Dynamical Gu zwille app oach
The ime dependen e sion o he Gu zwille wa e unc ion (2.3) is ob ained
by allowing he Gu zwille ampli udes o depend on ime (i)
n( ) [54]. Then
34
Figu e 2.2: Mean-ﬁeld g ound s a e o a 2D op ical la ice. (a) he e ical
axis shows he alue o he densi y a each si e o he la ice, and he co -
esponding o de pa ame e is in (b). The alue o he ha monic oscilla o
equency is gi en by Ω = 0.0108 ×2πU/~.
he equa ions o mo ion o he ampli udes a e eadily ob ained by minimizing
he ac ion o he sys em, gi en by S=Rd L, wi h espec o he a ia ional
pa ame e s (i)
n( ) and hei complex conjuga es ∗(i)
n( ). The Lag angian o
he sys em in he quan um s a e |Ψi, is gi en by [60]
L=hΨ|˙
Ψi−h˙
Ψ|Ψi
2i−hΨ|ˆ
H|Ψi,(2.5)
whe e |˙
Ψiis he ime de i a i e o he wa e unc ion (2.3). By equa ing o
ze o he a ia ion o he ac ion wi h espec o ∗(i)
n, one ge s he equa ions
i~d
d (i)
n=−Jh¯ϕi√n (i)
n−1+ ¯ϕ∗
i√n+ 1 (i)
n+1i
+"U
2n(n−1) + nX
j6=i
Vi,jhˆnji−µin# (i)
n,(2.6)
whe e ¯ϕi=Phjiiϕj, he sum uns o e all nea es neighbo s jo si e i,
hˆnji=hΨ|ˆa†
jˆaj|Ψiis he a e age pa icle numbe a si e j, and he o al
numbe o pa icles is gi en by N=Pihˆnii. I is no diﬃcul o e i y he
commu a ion ela ion [ ˆ
N, ˆ
HBH] = 0, which implies ha he o al numbe o
35
Bosons is a conse ed quan i y o dynamics in he eal ime [46]. These
equa ions a e o mean-ﬁeld ype, because hey a e w i en o a single si e
iand he ”ﬁeld” ¯ϕi oge he wi h Pj6=iVi,jhˆnji, ep esen he inﬂuence o
neighbo ing si es on he si e i, and ha e o be de e mined sel -consis en ly.
Eqs. (2.6) a e also a se o coupled equa ions, he coupling a ising om he
unneling pa , and can be w i en in he ma ix o m
i~d
d ~
=M[~
, µ, U, J]·~
, (2.7)
whe e ~
=h (1)
0, (1)
1,···, (i)
n,··· (NS)
nmax iT, is he ec o o he Gu zwille
ampli udes o de ed om si e 1 o si e NS, he la e being he o al numbe
o si es. I is wo h no icing ha he ma ix M[~
, µ, U, J] is i sel a unc ional
o he coeﬃcien s ~
h ough he ﬁelds ¯ϕiand Pj6=iVi,jhˆnji, which ha e o be
calcula ed in a sel -consis en way. Le us cla i y his poin wi h an example.
Suppose we wan o sol e Eq. (2.7) be ween an ini ial ime i= 0 and a ﬁnal
ime , wi h a gi en ini ial condi ion ~
(0). We disc e ize he ime in e al
in Ns eps o size ∆ , wi h Nﬁni e, and deﬁne s=s∆ such ha s=0 ≡0
and s=N≡ . The e o e, o calcula e he solu ion a a ce ain poin in ime
~
( s+1) we need o know he solu ion igh a he p eceding ime ~
( s), wi h
which we can compu e he ﬁelds ha in u n de e mine M[~
( s), µ, U, J],
and he solu ion is eadily ound o be
~
( s+1) = e−iM[~
( s),µ,U,J]∆ /~~
( s).(2.8)
S a ing om s= 0, in N+ 1 s eps we ha e de e mined he solu ion a he
desi ed ime . A he compu a ional le el, his is he simples p ocedu e
one can implemen o calcula e he dynamics o he sys em. Howe e , one
needs o be ca e ul in he choice o he ime s ep ∆ , especially o as -
oscilla ing dynamics. In such cases, a Runge-Ku a wi h adap i e s epsize
con ol has p o en o be mo e eﬃcien . Ins ead o s iﬀ dynamics, whe e
he solu ion p esen s bo h slowly- a ying and as oscilla ing egions, he
simple p ocedu e desc ibed abo e may be enough accu a e as in he case o
imagina y ime e olu ion.
Equa ions (2.7) can be sol ed in eal ime o also imagina y ime τ=
i . The imagina y ime e olu ion is a s anda d echnique ha has been
ho oughly used, because due o dissipa ion is supposed o con e ge o he
g ound s a e o he sys em. Two hings a e wo h o be no iced. Fi s ,
because he imagina y ime e olu ion is no uni a y, i does no conse e he
no m o he Gu zwille wa e unc ion, which has o be eno malized a e each
ime s ep. Second, he o al numbe o pa icles is no a conse ed quan i y
36
any mo e. Fo dipola Hamil onians he imagina y ime e olu ion does no
always con e ge o he ue g ound s a e and i ge s blocked in conﬁgu a ions
which a e a local minimum o he ene gy. On he one hand his makes i a
diﬃcul ask o iden i y he g ound s a e o such sys ems, and on he o he
hand i is a signa u e o he exis ence o me as able s a es as we will discuss
in de ails in he nex pa .
2.2.2 Pe u ba i e mean- ield app oach
A mo e con enien me hod o de e mine he insula ing phases o a dipola
Hamil onian is o use a mean-ﬁeld app oach pe u ba i e in ϕi. F om s a-
is ical mechanics, he expec a ion alue o he annihila ion ope a o a he
i- h si e is gi en [56] by he ace
ϕi=hˆaii= T (ˆaiˆρ),(2.9)
whe e ˆρ=Z−1e−βˆ
His he densi y ma ix ope a o , Z= T (e−βˆ
H) i s no -
maliza ion, and β= 1/KBTis he in e se empe a u e o he sys em. We
w i e Hamil onian (1.54) in he o m ˆ
H=ˆ
H0+ˆ
H1whe e
ˆ
H0=U
2X
i
ˆni(ˆni−1) −µX
i
ˆni+X
i6=j
Vi,j
2ˆniˆnj(2.10)
ˆ
H1=−JX
hiji
ˆa†
iˆaj,(2.11)
and we assume a uni o m chemical po en ial µ. The gene aliza ion o a si e-
dependen chemical po en ial is s aigh o wa d. Fu he mo e, we assume
low empe a u es β→ ∞, and he unneling coeﬃcien o be he smalles
ene gy in he sys em, i.e J≪U, µ, Vi,j such ha we can ea ˆ
H1as a small
pe u ba ion on ˆ
H0, and use he Dyson expansion a he ﬁ s o de in ˆ
H1 o
all he exponen ial ope a o s, so ha one ob ains
e−β(ˆ
H0+ˆ
H1)≃e−βˆ
H0ˆ
11−Zβ
0
eτˆ
H0ˆ
H1e−τˆ
H0dτ.(2.12)
We now w i e Hamil onian (2.11) as a sum o single si e Hamil onians. W i -
ing he annihila ion ope a o as ˆai=ˆ
Ai+ϕi, we can pe o m he mean ﬁeld
decoupling on he unneling e m
ˆa†
iˆaj=ˆ
A†
iϕj+ˆ
Ajϕi+ϕiϕj+ˆ
A†
iˆ
Aj
≃ˆa†
iϕj+ ˆajϕi−ϕiϕj,(2.13)
37

whe e in he las s ep we ha e assumed small ﬂuc ua ions, cha ac e is ic o
he Mo , o he deep supe ﬂuid s a es, and eplaced ˆ
A†
iˆ
Aj≃0. In Hamil-
onian (2.11) we now eplace ˆa†
iˆajwi h he exp ession calcula ed abo e, we
neglec e ms o he o de o ϕ2and ﬁnd he mean ﬁeld unneling Hamil onian
ˆ
HMF
1=−JX
iˆa†
i¯ϕi+ ¯ϕ∗
iˆai.(2.14)
Gi en a classical dis ibu ion o a oms in a la ice such as
|Φi=Y
i|niii,(2.15)
sa is ying ˆ
H0|Φi=EΦ|Φi, le us suppose ha his conﬁgu a ion is a local
minimum o he ene gy, i can be he g ound s a e, namely he absolu e
minimum, o ano he local minimum. We will be mo e igo ous a he end
o his sec ion ega ding he meaning o local minimum o ene gy bu o
he momen le us e e o he common pic u e o a local minimum. In he
basis o he eigen unc ions o ˆ
H0, sa is ying he ela ion ˆ
H0|Υi=EΥ|Υi he
pa i ion unc ion hen akes he simple o m
Z≃T (e−βˆ
H0) = X
|ΥihΥ|e−βˆ
H0|Υiβ7→∞
−→ e−βEΦ,(2.16)
whe e he las limi holds because we do no ace o e all he s a e o he
basis bu only a ound he s a e |Φi, which is assumed o be a local minimum
o he ene gy. Using again a Dyson expansion o he exponen ial o he densi y
ope a o , we ob ain he o de pa ame e as
ϕi≃ −eβEΦZβ
0
T hˆaie−(β−τ)ˆ
H0ˆ
HMF
1e−τˆ
H0idτ
=J¯ϕieβEΦZβ
0X
|ΥihΥ|ˆaie−(β−τ)ˆ
H0ˆa†
ie−τˆ
H0|Υi,(2.17)
which is easy o calcula e. The ace is hen non i ial only o |Υi=|Φi
and |Υi=ˆai
√ni|Φi, whe e niis in ege on |Φi, and a e he in eg a ion in he
β7→ ∞ limi we a e le wi h he esul
ϕi=J¯ϕini+ 1
Ei
P
+ni
Ei
H,(2.18)
whe e he quan i ies Ei
P,Ei
Ha e deﬁned as
Ei
P=−µ+Uni+V1,i
dip
Ei
H=µ−U(ni−1) −V1,i
dip,(2.19)
38
and a e espec i ely he ene gy cos o a pa icle (P) and hole (H) exci a ion
on op o he |Φiconﬁgu a ion. In he p e ious exp essions V1,i
dip =Pj6=iVi,jnj
is he dipole-dipole in e ac ion ha eels one a om placed a si e i, wi h he
es o he a oms in he la ice. We pe o med he in eg al (2.17) in he limi
o β7→ ∞, and in such a limi one ﬁnds ha he in eg al con e ges only o
posi i e alues o he pa icle and hole exci a ion ene gies, namely
U(ni−1) + V1,i
dip < µ < Uni+V1,i
dip.(2.20)
This equi emen s ha e o be ulﬁlled a e e y si e io he la ice and hey
simply s a e ha he conﬁgu a ion |Φiis a local minimum wi h espec o
adding and emo ing pa icles a any si e. In he ligh o his s a emen ,
he es ic ion on he ace o Eq. (2.16) is now igo ous, and is in pe ec
ag eemen wi h he ea men done in [1]. No ice, ha i |Φiis no a local
minimum, hen one ﬁnds ha condi ions (2.20) a e ne e sa isﬁed and he
in eg al (2.17) indeed di e ges. This ea men is o cou se also alid o |Φi,
being in pa icula he g ound s a e o he sys em.
One ﬁnds such an equa ion (2.18), and condi ions (2.20) o e e y si e i
o he la ice. The con e gence condi ions a e simple and among hem one
has o choose he mos s ingen o ﬁnd he bounda y o he lobe a J =
0. Ins ead he equa ions o he o de pa ame e s a e coupled due o he
¯ϕi e m, hey can be w i en in a ma ix o m M(µ, U, J)·~ϕ = 0, wi h
~ϕ ≡(···ϕi···), and ha e a non i ial solu ion. Fo e e y µ, he smalles J
o which de [M(µ, U, J)] = 0 gi es he lobe o he |Φiconﬁgu a ion in he
J s. µplane.
2.2.3 Pe u ba i e mean- ield s. dynamical Gu zwille
app oach
The p edic ions o he pe u ba i e mean-ﬁeld ea men a e in pe ec ag ee-
men wi h he esul s o he dynamical Gu zwille app oach, since hey bo h
ely on he same mean-ﬁeld app oxima ions. The ﬁ s looks a he s abili y
o a gi en densi y dis ibu ion |Φi=Qi|niiio in ege nia oms pe si e,
wi h espec o pa icle and hole exci a ions, while he la e minimizes he
ene gy o a andom ini ial conﬁgu a ion wi h espec o pa icle and hole ex-
ci a ions leading o he dis ibu ion |Φii he ini ial condi ion is suﬃcien ly
close. Howe e , he ﬁ s me hod can only iden i y he phase bounda ies o
he insula ing lobes wi hou p o iding any u he in o ma ion on he SF
phases ou side he lobes, which can ins ead be explo ed wi h he imagina y
ime e olu ion. Ne e heless o dipola Hamil onians, due o he p esence
o many local minima o he ene gy, as we will see in he nex pa [57, 58],
39
i is e y diﬃcul o iden i y he g ound s a e wi h he dynamical Gu zwille
app oach. This can be achie ed mo e eﬃcien ly h ough he pe u ba i e
mean-ﬁeld app oach. The e o e he wo me hods complemen each o he .
As an example, in Fig 2.1 (a,b) he black lines a e calcula ed wi h he pe -
u ba i e me hod (Vi,j = 0) while he SF egion ou side he lobes is explo ed
using imagina y ime e olu ion showing pe ec ag eemen wi h he wo ap-
p oaches.
40
Pa II
Me as able s a es
41
0 2 4 6 8
−5
−4
−3
−2
−1
0
µ/UNN
E/UNNNS
(b)
0 2 4 6 8
50
100
150
200
250
300
350
400 (a)
µ/UNN
numbe o s a es
Figu e 3.3: (a) Numbe o me as able s a es, and (b) ene gy o he g ound
( hick line) and me as able s a es ( hin lines) as unc ion o µ, o o
U/UNN = 20, and a ange o he dipole-dipole in e ac ion cu a he ou h
nea es neighbo . The inse shows he ene gy le els a ﬁlling ac o 1/2.
be ween |Φiini ial and |Φi inal, and (iii) h ough he a ia ion o qwe calcu-
la e he minimal ac ion S0, wi h he imagina y ime Lag angian, along he
s a iona y pa h s a ing a |Φiini ial; his pa h is called an ins an on pa h, in
sho ins an on. I connec s |Φiini ial and |Φi inal only i he wo s a es a e
degene a e, o he wise he s a iona y pa h connec s |Φiini ial wi h an in e me-
dia e s a e called he bouncing poin |Φibounce. We ge an es ima e o he
ene gy ba ie sepa a ing he wo s a es by e alua ing he Lag angian om
|Φiini ial o |Φi inal and imposing ze o ”momen um” P= 0, i.e. as one would
do in he Lag angian o a classical pa icle in a po en ial.
Once he minimal ac ion S0is known, hen he unneling ime Tis eadily
calcula ed [59] as
ω0T=π
2eS0,(3.5)
whe e ω0is o he o de o he equency o he ypical small oscilla ions o
|Φiini ial a ound he local minimum o he ene gy. In analogy o a classical
pa icle unneling h ough a ba ie , he ins an on has he nice in e p e a-
ion o he s a iona y pa h connec ing he wo local minima in he in e ed
po en ial, as schema ically ep esen ed in Fig. 3.4 (b).
48

(a) (b) (c)
|Φ〉 inal
|Φ〉 ini ial
ω0
T
T
Figu e 3.4: (a) Pa icle in a minimum o a po en ial ba ie , he pa icle
oscilla es wi h equency ω0a ound he local minimum and unnels in o he
igh well in a ime T; (b) he ins an on; and (c) he p ocess o which a
checke boa d s a e unnels in o he an i-checke boa d ha shown comple e
exchange o pa icle wi h holes and ice e sa. The p ocess happens in a
ime Tin analogy wi h (a).
The imagina y ime Lag angian o a sys em [60], desc ibed by a quan um
s a e |Φi, is gi en by
L=−h˙
Φ|Φi−hΦ|˙
Φi
2+hΦ|ˆ
H|Φi,(3.6)
wi h |˙
Φiindica ing he ime-de i a i e, and ˆ
Hbeing he Hamil onian o he
sys em. In he app oxima ion whe e |Φiis he Gu zwille wa e unc ion o
a gi en me as able s a e, we w i e i s ampli udes as
(i)
n=1
√2x(i)
n+ip(i)
n,(3.7)
whe e x(i)
n, and p(i)
na e eal numbe s which a e going o be ela ed o he
a ia ional pa ame e s, and hei conjuga e momen a in he ollowing. Fo
49
simplici y, we conside s a es wi h a maximum occupa ion numbe o nmax =
1 (i.e. n= 0,1), and he e o e we ha e a o al o 4NSpa ame e s wi h NS
being he o al numbe o si es. The Lag angian (3.6) becomes a unc ional
o he 4NSpa ame e s, namely
L[x(i)
n, p(i)
n] = −i
1
X
i,n=0
p(i)
n˙x(i)
n+hΦ|ˆ
H|Φi,(3.8)
as well as he expec a ion alue o he Hamil onian hΦ|ˆ
H|Φi. In o de o
pu he Lag angian (3.8) in i s canonical o m we mus in oduce he new
coo dina es q(i)
n=x(i)
nand hei conjuga e momen a P(i)
n=∂L/∂ ˙q(i)
n=−ip(i)
n,
and we can w i e
L[q(i)
n, P(i)
n] =
1
X
i,n=0
P(i)
n˙q(i)
n−H[q(i)
n, P(i)
n],(3.9)
whe e H[q(i)
n, P(i)
n] = −hΦ|ˆ
H|Φiis a cons an o he mo ion 1. We now wan
o educe he dynamic desc ibed by he Lag angian (3.9) o a one dimensional
p oblem, desc ibed only by one a iable qand i s conjuga e momen um P.
Th ough he a ia ion o (q, P) we wan o desc ibe he in e change be ween
he s a e |Φiini ial and |Φi inal, as o example he one ep esen ed in Fig.
3.4(c). In Sec. 3.3.1, we show how o educe he numbe o a ia ional
pa ame e s o one, q, and i s conjuga e momen um P, by making use o a
a ia ional Ansa z as well as he no maliza ion condi ion on he coeﬃcien s
(3.7) and he conse a ion o he o al numbe o pa icles. These condi-
ions en e in o he exp ession o he Lag angian (3.9) h ough Lag ange
mul iplie s λc. Consequen ly he equa ions o mo ion gi en by ˙q=∂H/∂P,
and ˙
P=−∂H/∂q a e go e ned by an Hamil onian which also includes he
cons ain s as ollows
H=H[q, P] + X
c
λcCc,(3.10)
whe e an explici exp ession o he condi ions Ccwill be gi en in Sec. 3.3.1.
The ac ion is hen eadily calcula ed along he s a iona y pa h o Eq. (3.10)
as ollows
S0=ZL[q, P]dτ=Zpa h L[q, P ]dq
˙q,(3.11)
wi h ˙q=∂H/∂P om Eq. (3.10).
1no e ha in he analogy o a classical pa icle in a po en ial V(x), he conse ed
quan i y in he imagina y ime would be H=P2
2m−V(x), which desc ibes he pa icle’s
mo ion in he in e ed po en ial.
50
3.3.1 Pa ame iza ion
Th ough he a ia ion o (q, P), we aim a desc ibing he in e change be-
ween he s a e |Φiini ial and |Φi inal, as o example he one schema ically
ep esen ed in Fig. 3.4(c). Du ing his p ocess he e a e si es ini ially oc-
cupied ha emp y, like he blue- amed si e o Fig. 3.4(c), which we call
he (B)-si e, and ice e sa, like he si e on he le o B, which is ini ially
emp y and occupied a he end, and we name he (A)-si e. When he ini ial
and ﬁnal s a e a e non-degene a e, as o example he case ske ched in Fig.
3.5(d), he e a e also si es ha do no change and emain ei he ull (F) o
emp y (E).
Du ing his p ocess, he Gu wille ampli udes (3.7) ha e o be no malized
a each si e, and he o al numbe o pa icles has o be conse ed, namely
| (i)
0|2+| (i)
1|2= 1,∀i(3.12)
NS
X
i=1 | (i)
1|2=N, (3.13)
whe e NSis he o al numbe o si es. We choose (q, P )≡(qB
0, PB
0) o be
he a ia ional pa ame e s o he blue- amed si e, and he no maliza ion
condi ion (3.12) oge he wi h he conse a ion o he numbe o pa icles
(3.13) be ween A and B gi e us h ee coupled equa ions
q2−P2+ (qB
1)2−(PB
1)2= 2
(qA
0)2−(PA
0)2+ (qA
1)2−(PA
1)2= 2
(qB
1)2−(PB
1)2+ (qA
1)2−(PA
1)2= 2.
(3.14)
As explained in [58], we make use o he ollowing Ansa z
qA
1=q
PA
1=P
PA
0=PB
1=−P
qA
0=qB
1
(3.15)
wi h which i is clea ha a he alue o (q, P) = (0,0) co esponds a
si ua ion in which si e A is emp y and si e B is ull, while o (q, P) = (√2,0)
he con a y is ue. Fo degene a e ini ial and ﬁnal s a es as in he case o
Fig. 3.4(c), he emaining si es hey ei he beha e like A o B, which implies
51
ano he se o condi ions summa ized as ollows
q(i)
0=qA(B)
0
P(i)
0=PA(B)
0
q(i)
1=qA(B)
1
P(i)
1=PA(B)
1
(3.16)
depending on whe he he si e iis ini ially emp y (A) o occupied (B). In-
s ead, when he ini ial and ﬁnal s a e a e non-degene a e, as is he case
conside ed in Fig. 3.5(d), he e a e also si es ha we assume no o change
and emain ei he ull (F) o emp y (E), Fo hese si es he cons ain s a e
espec i ely gi en by
qF
0= 0
qF
1= 2
PF
0=PF
1= 0,
(3.17)
and
qE
0= 2
qE
1= 0
PE
0=PE
1= 0.
(3.18)
All hese condi ions (3.14)-(3.18), which we name Cc, en e explici ly in o he
calcula ion o Hamil onian (3.10).
3.3.2 Ac ion and unneling ime
In Fig. 3.5 (a,b) we plo he minimal ac ion di ided by he o al numbe
o si es NSo he cell, as a unc ion o he unneling coeﬃcien J, o wo
diﬀe en p ocesses. The ﬁ s one (a,c), in which ini ial and ﬁnal s a e a e
degene a e, shows he exchange o pa icles wi h holes in he whole la ice,
and is ske ched in he lowe pa o Fig. 3.5 (c) whe e we also plo he
po en ial ba ie be ween ini ial and ﬁnal s a e calcula ed as −H(q, P = 0).
Ins ead, in he second one (b,d) he ﬁnal s a e is he g ound s a e, i.e. deepe
in ene gy wi h espec o he ini ial s a e, and only a ew si es o he la ice
exchange pa icles wi h holes du ing he p ocess, as ske ched in he uppe
pa o Fig. 3.5 (d) along wi h he po en ial ba ie . A side ema k, he poin
whe e he hick line o he ba ie encoun e s he dashed line is he bouncing
poin |Φibounce.
52
0 0.12
0
2.6
J/UNN
S0/NS
(a)
0
0.47
0.61
0.87
q
E/NSUNN
(c)
√2
0 0.12
0
2.6
J/UNN
(b)
0
0.47
0.61
0.87
q
(d)
√2
Figu e 3.5: (a,b) Ac ion pe si e and (c,d) ene gy ba ie o he p ocess
ske ched in panels (c,d). In bo h cases he ini ial s a e is he conﬁgu a ion
(IIb) o Fig. 3.2 and he alue J= 0.12UNN co esponds o he ip o i s
insula ing lobe. The ﬁ s one (a,c) is o degene a e ini ial and ﬁnal conﬁgu-
a ions while o he second one (b,d) he ﬁnal conﬁgu a ion is ene ge ically
deepe . The diﬀe ence in he wo p ocesses mani es s also in he heigh o
he ba ie which is smalle o he second case, leading o a smalle ac ion
and consequen ly a smalle li e- ime.
The ac ion in gene al di e ges o J→0 indica ing a di e gen unneling
ime T, and hen dec eases mono onically up o a minimum alue in co e-
spondence o he ip o he lobe, J= 0.12UNN he e, signaling a minimum
li e- ime a he ip o he lobe, as expec ed. In be ween hese wo ex eme
beha io s, he ac ion inc eases mono onically wi h he numbe o si es in-
ol ed in he exchange o pa icles wi h holes; he mo e si es in ol ed as in
he case o Fig. 3.5 (a,c), he bigge he ac ion is. Summa izing, om he
ﬁgu es abo e, we conclude ha small ene gy diﬀe ences be ween he ini ial
|Φiini ial and he ﬁnal s a es |Φi inal and la ge egions o he la ice unde going
pa icle-hole exchange in he unneling p ocess con ibu e o la ge ba ie s,
i.e. long li e imes T. On he con a y, o big ene gy diﬀe ences and small
egions o he la ice unde going pa icle-hole exchange, he ba ie is small.
Hence, in gene al i is mo e likely o a gi en s a e o unnel in o a s a e
deepe in ene gy, e.g., he g ound s a e, han in o i s complemen a y, which
53

implies he exchange o pa icles wi h holes in he whole la ice.
3.4 P epa a ion, manipula ion and de ec ion
Ve y impo an issues a e he p epa a ion and de ec ion o he a omic s a es
in he la ice. In he expe imen s wi h cold gases and op ical la ices, he
ypical p ocedu e is ﬁ s o ob ain a condensa e in a ha monic ap and
hen adiaba ically ump up he op ical po en ial. The e o e one may ask he
ques ion how o each a desi ed conﬁgu a ion o whe he i is possible o
each he g ound s a e in his way, since we ha e discussed in Sec. 3.2 ha
o la ge la ices he e exis many conﬁgu a ions wi h localized de ec s ha
compe e wi h he g ound s a e.
One can use supe la ices in o de o p epa e he a oms in conﬁgu a ions
o p e e en ial symme y. We ha e checked ha he p esence o de ec s is
s ongly educed when a local po en ial ene gy ollowing desi ed pa e ns is
added o he op ical la ice. No e ha he conﬁgu a ions ob ained in such
a way will also emain s able once he supe la ice is emo ed, hanks o
dipole-dipole in e ac ion. Mo eo e , in [58] we ha e demons a ed ha by
using supe la ices he ans e om one me as able conﬁgu a ion o ano he
necessa ily occu s ia supe ﬂuid s a es, and can be con olled ully a he
quan um le el. As discussed in Sec. 3.4.1, he ans e is a quan um con-
olled p ocess whe e he con ol pa ame e s a e among o he s he unneling
coeﬃcien J, and he magni ude o he local chemical po en ial ∆µ ollowing
a desi ed pa e n. E en i a he MI-SF ansi ion i is impossible o be adi-
aba ic because o he con inuous exci a ion spec um o he SF phase [10],
o a ce ain ange o con ol pa ame e s he p ocess wo ks, and as discussed
in Sec. 3.4.1, is qui e obus . In Fig. 3.6 we show an example o such a
p ocess in which we ans e he checke boa d g ound s a e (CB) in o he
me as able s a e (IIa) o Fig. 3.2, wi h a 99% ﬁdeli y.
The local chemical po en ials ollow he s ipe pa e n o he occupied si es
o he me as able s a e (IIa), and a e changed smoo hly in ime as shown
in Fig. 3.6 (a). A he same ime, he unneling coeﬃcien is a ied as he
smoo hed s ep-like unc ion o Fig. 3.6 (b) such as o exi he ip o he
CB insula ing lobe d awn as a dashed ed line, emain in he SF phase o
an app op ia e amoun o ime, and subsequen ly en e om he ip o he
Mo insula ing me as able s a e (IIa) ep esen ed by he blue hick line.
No ice ha since we a e pe o ming eal- ime dynamics, he o al numbe
o pa icles is a conse ed quan i y and he e o e he MI o SF ansi ion
happens only a he ip o he insula ing lobes. In Fig 3.6 (c) we plo he
54
0 28.95 91.1 120
0
1
popula ions
UNN
(c)
0 28.95 91.1 120
−3
−2.57
−0.45
0
UNN
∆µ/UNN
(a)
−3−2.57−0.450
0
0.66
J/UNN
∆µ/UNN
(b)
CB IIa
Figu e 3.6: (a) The pulse o local chemical po en ial as a unc ion o ime.
(b) The smoo hed s eplike unc ion is he unneling coeﬃcien as a unc ion
o ∆µ, while he hick (dashed) line is he ip o he IIa (CB) insula ing lobe.
(c) Popula ion in e sion, om CB o IIa a he end o he p ocess. No ice he
oscilla ion o popula ions when passing h ough he SF egion o he phase
diag am.
popula ions co esponding o he CB and (IIa) s a es deﬁned as he NS- h
oo o he ﬁdeli y
PMS =NS
q|hΦMS|Φi|2,(3.19)
whe e NSis he numbe o si es o he elemen a y cell. The dashed line
is he popula ion o he CB s a e while he hick line co esponds o he
popula ion o (IIa), which a he end o he p ocess s abilizes a he alue o
P(IIa) = 0.99.
The spa ially modula ed s uc u es c ea ed, and manipula ed in such a
way can be de ec ed ia he measu emen o he noise co ela ions o he
expansion pic u es [61, 77, 63]: he o de ed s uc u es in he la ice gi e ise
o diﬀe en pa e ns in he spa ial noise co ela ion unc ion
C(~
d) = Rd2xhρ o (~x +~
d/2) ρ o (~x −~
d/2)i
Rd2xhρ o (~x +~
d/2)i hρ o (~x −~
d/2)i
≈X
j,k
exp him
~ ~
d·(~xj−~xk)iρjρk=|F(ρ)|2,(3.20)
whe e we ha e named ρ o he densi y dis ibu ion a e ime o ﬂigh , while
55
ρkis he densi y dis ibu ion in he la ice. This is no hing else han he
modulus squa e o he Fou ie ans o m o he densi y dis ibu ion in he
la ice. Such a measu emen is in p inciple able o ecognize he geome y
o he densi y pa e n in he la ice as well as he p esence o de ec s in he
densi y dis ibu ion, which could be exac ly econs uc ed s a ing om he
pa e ns in he spa ial noise co ela ion unc ion. In Fig. 3.7 (I,II,II), he
lowe panels show he noise co ela ion unc ions o he densi y dis ibu ions
a ﬁlling close o 1/2 shown in he uppe panels whe e we ha e assumed a
localised Gaussian densi y dis ibu ion a each la ice si e. The p esence
o de ec s in he densi y dis ibu ion can be in p inciple de ec ed wi h his
me hod.
Fo he momen , he signal o noise a io equi ed o single de ec ecog-
ni ion is beyond p esen expe imen al possibili ies. Howe e , by a e aging
o e a ﬁni e numbe o diﬀe en expe imen al uns p oducing he same spa ial
dis ibu ion o a oms in he la ice, a good signal can be ob ained.
Ve y ecen ly i has been epo ed ha , by means o a high- esolu ion
op ical imaging sys em desc ibed in [64], single a oms a e de ec able wi h
nea -uni y ﬁdeli y on indi idual si es o a Hubba d- ype op ical la ice. The
au ho s epo a way o de e mine he p esence o an a om on a single si e o
he la ice, by measu ing he o al numbe o sca e ed pho ons pe la ice
si e. Howe e , du ing he imaging p ocess only emp y o singly occupied
si es can be seen in he image, because o molecule o ma ion on mul iply
occupied si es and ligh -assis ed collisions. This me hod could in p inciple
be used o obse e expe imen ally he diﬀe en densi y dis ibu ions o he
me as able s a es.
3.4.1 T ans e p ocess
In [58] we ha e s udied how o ans e popula ion om a gi en me as able
conﬁgu a ion o ano he one wi h a diﬀe en symme y, by changing he
la ice pa ame e s in ime, whe e he dynamics ha e been desc ibed h ough
he mean-ﬁeld equa ions de i ed in Sec. 2.2.1. Speciﬁcally, we ha e s udied
how o ans e popula ion om he CB s a e o he me as able s a e (IIa)
o Fig. 3.2, by applying ime dependen local chemical po en ials in a o
o he s a e (IIa), and by changing he unneling coeﬃcien Jin ime, so as
o exi he CB lobe and, h ough he supe ﬂuid egion, en e in o he (IIa)
lobe. Ideally, he popula ion o (IIa) a he end o he p ocess has o be
one, P(IIa)( in) = 1, bu he ac ual alue o P(IIa)( in) is e y much sensi i e
on he exac alues he pa ame e s ake du ing he dynamics. Speciﬁcally,
he magni ude o he local chemical po en ials is smoo hly a ied in ime as
56
0
0.5
−−
1
0
0.5
−−
1
(I) (II) (III)
Figu e 3.7: The lowe panels show he spa ial noise co ela ion pa e ns o
conﬁgu a ions (I) o (III) in he uppe pannels, assuming a localised Gaussian
densi y dis ibu ion a each la ice si e. Figu e om [57].
ollows
∆µ( ) = −C anh hα
C( − 0)i+C anh h−α
C 0i,(3.21)
whe e C= 1.7UNN and 0= 60/UNN (in uni s o ~= 1), a e kep cons an ,
while αis a ee pa ame e ha se s he maximum slope o his unc ion.
Ins ead, he unneling coeﬃcien is dynamically chanced as ollows
J(∆µ) = Jm−J0
2min  anh [−smo]− anh [s(∆µ−mo)] + 2J0
Jm−J0
anh [s(∆µ−mi)] − anh [−smi] + 2Jm
Jm−J0,
(3.22)
wi h s= 15/UNN and J0= 0.02UNN cons an s, mi=−2.6UNN ﬁxes he
supe ﬂuid o Mo insula o ansi ion poin a ∆µin =−2.57UNN , whe eas
Jmand moa e ee pa ame e s ela ed wi h he maximum alue o unneling
coeﬃcien in he supe ﬂuid egion and he poin ∆µowhe e he CB ceases
o exis . Toge he wi h he in ensi y I o he andom noise ha ﬁxes he
ini ial condi ion, which, as discussed in [58], i is necessa y in o de no
o ha e i ial dynamics, he space o ou con ol pa ame e s is in o al 4-
dimensional, and is gi en by
{α, ∆µo, Jm, I }.(3.23)
57
In ou wo k, we ha e s udied a sys em composed o a s ack o wo 2D
op ical la ice laye s wi h he dipoles pola ized pe pendicula ly o he 2D
planes (see Sec. 5.1). Because he unneling be ween diﬀe en laye s is sup-
p essed he gas is s able agains collapse, and by including he a ac i e pa
o he dipola in e ac ion in he pe pendicula di ec ion we ﬁnd ha pa icles
belonging o he wo diﬀe en planes may bind oge he in a composi e. The
composi es ea u e no only he pai MI and pai SF phases bu also a no el
pai -supe solid (PSS) quan um phase. In he con ex o Bose mix u es in a
one-dimensional op ical la ice, a small egion o PSS phase was ound in he
p esence o in a-species on-si e epulsion, and on-si e a ac ion be ween he
wo species [69].
The wo k is o ganized as ollows. In Sec. 5.1 we in oduce and explain he
de ails o he model, we s udy he g ound s a e o he sys em by means o he
pe u ba i e mean ﬁeld app oach de i ed in Sec. 2.2.2 based on a Gu zwille
Ansa z. We will see ha his me hod is no accu a e enough o desc ibe he
ue g ound s a e o he sys em, and he eason lays behind he na u e o
he lowes lying exci a ions. In ac , in he limi o pa ame e s we conside ,
we demons a e ha i is ene ge ically a o able o dope he sys em wi h a
pai o pa icles (pai o holes) ins ead o a single pa icle (single hole).
In Sec. 5.2, we show ha he sys em admi s a desc ip ion in e ms o a
low-ene gy subspace o pai s, we de i e he eﬀec i e Hamil onian ˆ
He o he
subspace by means o pe u ba ion heo y up o second o de in unneling.
We de i e a new pe u ba i e mean ﬁeld app oach o ˆ
He , simila ly as in
Sec. 2.2.2, which only admi s pai o pa icles o pai o holes exci a ions and
we ﬁnd he insula ing g ound s a e lobes o he sys em.
In Sec. 5.3, we gene alize he dynamical Gu zwille app oach de i ed in
Sec. 2.2.1, and de i e he dynamics equa ions o he Gu zwille ampli udes
in he low-ene gy subspace. These equa ions pe mi us o in es iga e he
phases o he sys em ou side he insula ing lobes and we ﬁnd s ong e idences
o he exis ence o a supe solid phase o pai s. We discuss he limi o alidi y
o ou desc ip ion in e ms o an eﬀec i e Hamil onian, and we conclude in
chap e 7.
Ou esul s a e based on he publica ion:
•C. T e zge , C. Meno i, and M. Lewens ein, Pai -Supe solid Phase in
a Bilaye Sys em o Dipola La ice Bosons. Physical Re iew Le e s,
103, 035304, (2009).
•C. T e zge , M. Alloing, C. Meno i, F. Dubin, and M. Lewens ein,
Coun e low Supe solid o an i-pola ized dipola Bosons in a 2D op ical
la ice. In p epa a ion.
64

Chap e 5
Dipola Bosons in a bilaye
op ical la ice
5.1 The model
In [74] we conside pola ized dipola pa icles in wo decoupled 2D op ical
la ice laye s (see Fig. 5.1), whe e he po en ial ba ie be ween he wo laye s
is la ge enough o p e en any in e -laye hopping. This is he simples mul i-
laye s uc u e and can be ob ained by using aniso opic op ical la ices o
supe la ices, which can exponen ially supp ess unneling in one di ec ion.
The in-plane dipola in e ac ion is iso opic and epulsi e. The in e laye
in e ac ion depends on he ela i e posi ion be ween he wo dipoles, bu is
domina ed by he nea es -neighbo a ac i e in e ac ion W < 0 be ween wo
a oms a he same la ice si e in diﬀe en laye s. We include only nea es -
neighbo (NN) in-plane (UNN) and ou -o -plane (W) dipola in e ac ions.
Since unneling is supp essed be ween he laye s pa icles belonging o he
diﬀe en laye s canno mix and beha e in p ac ice like wo diﬀe en species
1. The p oblem is analogous o ha o wo bosonic species on a 2D op ical
la ice wi h an in e -species a ac ion W < 0 a he same la ice si e, and
in a-species epulsion UNN. The ela i e s eng h be ween UNN and Wcan
be uned by changing he spacing d⊥be ween he wo laye s, ela i e o he
2D op ical la ice spacing d. Because o he dependence o he dipole-dipole
in e ac ion like he in e se cubic powe o he dis ance, he a io |W|/UNN
can be uned o e a wide ange. While i can be negligible o d⊥≫d
making he sys em asymp o ically simila o a single 2D la ice laye , i can
1Because o his analogy we will o en e e o he wo laye s as he wo species and
ice e sa.
65
Figu e 5.1: Schema ic ep esen a ion o wo 2D op ical la ice laye s pop-
ula ed wi h dipola bosons pola ized pe pendicula ly o he la ice plane.
The pa icles eel epulsi e on si e Uand nea es -neighbo UNN in e ac ions.
In e laye unneling is comple ely supp essed, while a nea es -neighbo in e -
laye a ac i e in e ac ion Wis p esen .
also become ele an and gi e ise o in e es ing physics, no exis ing in he
single laye model as poin ed ou in [65, 66, 67, 70, 71, 72, 73].
The sys em is desc ibed by he Hamil onian
ˆ
H=X
i,σ 
U
2ˆnσ
i(ˆnσ
i−1) + X
hjii
UNN
2ˆnσ
iˆnσ
j−µˆnσ
i

+WX
i
ˆna
iˆnb
i−JX
hiji
[ˆa†
iˆaj+ˆ
b†
iˆ
bj],(5.1)
whe e σ=a, b indica es he wo species (which in he speciﬁc case conside ed
he e a e a oms in he lowe and uppe 2D op ical la ice laye , espec i ely),
Uis he on-si e ene gy, UNN he in alaye nea es neighbo s epulsion, W
he in e laye a ac ion, J he in alaye unneling pa ame e , and µ he
chemical po en ial, as schema ically ep esen ed in Fig. 5.1. The pa ame e s
Uand Ja e equal o he uppe and lowe laye s and he chemical po en ials
µa e he same, since equal densi ies in he wo laye s a e assumed. No ice
ha since W < 0, i is necessa y o ha e U+W > 0 o a oid collapse. The
symbols hijiand hjiiindica e nea es neighbo s.
We ocus on he physical si ua ion in which he wo laye s a e e y close o
one ano he , namely d⊥≪d, because in his limi pa icles a he same la ice
si e io diﬀe en laye s pai in o composi es. The composi es localize in a MI
66
s a e o small alues o he unneling coeﬃcien , while o la ge alues o J
he pai s hop a ound in he op ical la ice o ming a Pai -supe ﬂuid (PSF)
phase [65]. Fu he mo e, he p esence o he long- ange in e ac ions leads o
he o ma ion o a no el pai -supe solid phase (PSS), namely, a supe solid o
composi es.
5.1.1 G ound s a e and single-pa icle single-hole ex-
ci a ions
We w i e he Hamil onian (6.1) as ˆ
H=ˆ
H0+ˆ
H1, whe e
ˆ
H0=X
i,σ 
U
2ˆnσ
i(ˆnσ
i−1) + X
hjii
UNN
2ˆnσ
iˆnσ
j−µˆnσ
i
+WX
i
ˆna
iˆnb
i(5.2)
ˆ
H1=−JX
hiji
[ˆa†
iˆaj+ˆ
b†
iˆ
bj],(5.3)
and we conside ˆ
H1 o be a small pe u ba ion on he in e ac ion e m
(5.2). Fo any gi en classical conﬁgu a ion o a oms in he la ice gi en
by |Φi=Qi|na
i, nb
ii, we can use he pe u ba i e mean ﬁeld app oach de-
i ed in Sec. 2.2.2 o analyze he s abili y o |Φiwi h espec o pa icle
and hole exci a ions. The chemical po en ials o he wo species a e he
same, which ﬁxes he same numbe o pa icles o be equal o bo h species,
and since W < 0, in he limi o close laye s he sys em na u ally ends o
minimize he diﬀe ences in he numbe o pa icles be ween he uppe and
lowe laye a he same la ice si e. A J= 0, his is easily unde s ood i we
w i e Hamil onian (5.2) in e ms o he sum and he diﬀe ence ope a o s,
ˆmi=ˆna
i+ ˆnb
i
2
ˆsi=ˆna
i−ˆnb
i
2,
(5.4)
and w i e explici ly he exp ession o he expec a ion alue o ˆ
H0on |Φi,
which is gi en by
hΦ|ˆ
H0|Φi=X
i
−(2µ+U)mi+ (U+W)m2
i+UNN X
hjii
mimj
+ (U−W)s2
i+UNN X
hjii
sisj
,(5.5)
67
wi h mi=hΦ|ˆmi|Φiand si=hΦ|ˆsi|Φi. Since W < 0, in he limi o
(U+W), UNN ≪U he ene gy abo e is minimized by se ing si= 0 and he
densi y is ﬁxed by mia each si e. We he e o e es ic ou sel es o s udy
he s abili y o such s a es
|αi=Y
i|ni, niii,(5.6)
wi h equal occupa ion o he wo species aand ba each si e, he e o e a
dis ibu ion o pai s composed by one a om o each species.
Fo any classical conﬁgu a ion (5.6), we can s udy he s abili y o |αi
wi h espec o pa icle and hole exci a ions using he same me hod de i ed
in Sec. 2.2.2, and calcula e he o de pa ame e ϕi=hˆaii=hˆ
biia each
la ice si e, gi en by
ϕi=J¯ϕimi+ 1
Ei
P
+mi
Ei
H,(5.7)
whe e mi=na
i=nb
i, and ¯ϕi=Phjiiϕi. The ene gy cos o a pa icle
(P) and a hole (H) exci a ion a he i- h si e o he conﬁgu a ion |αi, a e
espec i ely gi en by
Ei
P=−µ+Umi+V1,i
dip +Wmi
Ei
H=µ−U(mi−1) −V1,i
dip −Wmi,(5.8)
whe e we deﬁne V1,i
dip =UNN Phjiimj, as he in-plane dipola in e ac ion, i.e.
he dipole-dipole in e ac ion ha expe iences one a om posi ioned a si e io
he la ice, wi h he es o he pa icles belonging o he same plane. No ice
ha he exci a ions (5.8), diﬀe om he pa icle-hole exci a ions in a single
laye (2.19), only because o he p esence o he ou -o -plane dipola in e -
ac ion Wmi. Howe e we will see ha he p esence o Wis a om being
i ial, in pa icula i does no only induce a shi in he exci a ion ene gies
bu is esponsible o new and diﬀe en physical beha io s wi h espec o
he single laye model. Fo he e alua ion o he o de pa ame e (5.7) i is
necessa y o equi e he pa icle and hole exci a ions o be posi i e, which
esul s in he equa ions
U(mi−1) + V1,i
dip +Wmi< µ < Umi+V1,i
dip +Wmi.(5.9)
One ﬁnds such an equa ion (5.7) and condi ions (5.9) o each si e o he la -
ice, and simila ly as in Sec. 2.2.2, we w i e he coupled equa ions (5.7) in he
ma ix o m M(µ, U, J, UNN, W)·~ϕ = 0, wi h ~ϕ ≡(···ϕi···). The e o e, i
he conﬁgu a ion |αiis a local minimum wi h espec o he pa icle and hole
68
exci a ions (5.8), o e e y µ he smalles J o which de [M(µ, U, J, UNN, W)] =
0 gi es he lobe bounda y o he |αiconﬁgu a ion in he J s.µplane. We
ﬁnd ha in he limi (U+W)/U →0, asymp o ically all classical conﬁgu-
a ions |αide elop an insula ing lobe which end o o e lap one ano he as
illus a ed in Fig. 5.2 wi h blue hin lines.
0 0.02 0.04 0.06 0.08 0.1
−0.95
−0.9
−0.85
−0.475
−0.425
−0.375
−0.325
0.05
0.1
0.15
J/U
µ/U
Figu e 5.2: Insula ing lobes calcula ed wi h espec o single-pa icle and
single-hole exci a ions o all ﬁlling ac o s ν, anging om ν= 0 o ν= 3,
o a 2 ×2×2 elemen a y cell wi h pe iodic bounda y condi ions only in he
di ec ions pe pendicula o he o ien a ion o he dipoles (see Fig. 5.1). The
la ice pa ame e s a e gi en by UNN = 0.025Uand W=−0.95U, which can
be ob ained o an in e -laye dis ance d⊥= 0.37d. No ice he o e lap o all
he lobes which does no happen in a single laye si ua ion.
This si ua ion has o be compa ed o a single laye model in which no all
classical dis ibu ions ha e an insula ing lobe bu only some conﬁgu a ions
a some speciﬁc ﬁlling ac o s ν, and he g ound s a e lobes do no o e lap.
In his si ua ion, in o de o ﬁnd he g ound s a e one has o compa e he
ene gy o he insula ing conﬁgu a ions a each alue o he chemical po en ial
µ. The conﬁgu a ion wi h he lowes ene gy is hen he g ound s a e, as
illus a ed in Fig. 5.2 wi h he ed lobes along wi h a ske ch o he g ound
s a e conﬁgu a ions in he panels. Excep o he ﬁlling ac o ze o ν=
69

2
NSPimi= 0, which is ound o be below he ho izon al ed line a µ=
−0.475U(NSis he o al numbe o si es), he g ound s a e is a mul iple
o he ac ional ﬁlling ac o ν= 1/2 wi h a spa ial densi y dis ibu ion o
a checke boa d, doubly occupied checke boa d and so on depending on he
alue o he chemical po en ial µ.
5.2 Low-ene gy subspace and e ec i e Hamil-
onian
The eason behind his o e lapping lobes lays in he ac ha he lowes lying
exci a ions a e no o single-pa icle-hole o Eq. (5.8), bu a e a he ob ained
by adding o emo ing wo pa icles a a gi en si e o he |αiconﬁgu a ions.
In ac , in he limi o pa ame e s we a e conside ing he e, all he densi y
dis ibu ions (5.6) span a low-ene gy subspace which is ene ge ically well
sepa a ed om he es o he Hilbe space, he e o e he desc ip ion o
such a sys em, as well as i s exci a ions, can be done wi hin a low-ene gy
heo y es ic ed o an eﬀec i e Hamil onian ˆ
He ac ing only on he subspace
spanned by he |αi-s. The subspace is desc ibed by all classical dis ibu ions
o a oms in he la ice |αio Eq. (5.6) wi h equal occupa ion o he wo
species aand ba each si e, he e o e a dis ibu ion o pai s composed by
one a om o each species. In such a si ua ion he e a e wo ypes o p ocesses
ha depend on he ime scale we a e looking a he sys em [79]: (i) he slow
p ocesses which d i e he sys em h ough diﬀe en s a es (ene ge ically e y
close o one ano he ) o he low-ene gy subspace, and (ii) he as p ocesses
which couple he low ene gy subspace wi h he es o he Hilbe space
composed o high-ene gy s a es. The la e a e called i ual subspace and
a e coupled o he low-ene gy subspace h ough he unneling Hamil onian
ˆ
H1o Eq. (5.3) ia single pa icle hopping. The ele an i ual subspace is
ob ained om he s a es |αiby b eaking one composi e, namely
|γ(a)
ij i=ˆa†
iˆaj
pnj(ni+ 1)|αi
|γ(b)
ij i=ˆ
b†
iˆ
bj
pnj(ni+ 1)|αi,
(5.10)
as schema ically ep esen ed in Fig. 6.3 o a uni o m dis ibu ion |αio
one a om pe si e. All o he s a es a e no coupled o |αi ia single pa icle
hopping and hence do no con ibu e.
70
Eα
Eγ = Eα + U − UNN
Eβ = Eα + 2(U + W − UNN)
Figu e 5.3: Schema ic ep esen a ion o he low-ene gy subspace. The s a e
|αiis uni o mly occupied by one pa icle pe si e and i s unpe u bed ene gy,
i.e. a J= 0, is gi en by Eα. The s a es |αiand |βi(o unpe u bed
ene gy Eβ) belong o he subspace, hey a e connec ed h ough a second
o de p ocess in he unneling ia he he s a e |γi(o unpe u bed ene gy
Eγ), which belong o he i ual subspace. I is s aigh o wa d o no ice ha
in he limi o (U+W), UNN ≪U, he ene gies abo e sa is y he necessa y
condi ion Eα, Eβ≪Eγ o he exis ence o he subspace.
The ene gy diﬀe ence be ween he i ual s a es |γi, o ene gy Eγ, and he
s a es |αi, o ene gy Eα, is gi en by he sum o single pa icle plus single hole
exci a ion ene gies o he s a es |αi, which is o he o de o Ua J= 0, and
is minimized by he wid h o he lobes |αi(see, e.g., Fig. 5.2) a ﬁni e J.
Slow p ocesses d i e he sys em h ough diﬀe en s a es o he low ene gy
subspace, |αiand |βi, ia second o de unneling; his happens h ough a
as coupling wi h he i ual subspace. Since we a e in e es ed in he long
ime physics o he sys em, we ha e o a e age ou all he as p ocesses
and he e o e we w i e an eﬀec i e Hamil onian in he subspace o pai s, and
include unneling h ough second o de pe u ba ion heo y [77, 78]. In he
pai -s a e basis, he ma ix elemen s o such a Hamil onian in second o de
pe u ba ion heo y a e gi en by
hα|ˆ
He |βi=hα|ˆ
H0|βi− 1
2X
γhα|ˆ
H1|γihγ|ˆ
H1|αi
×1
Eγ−Eα
+1
Eγ−Eβ(5.11)
whe e ˆ
H0, gi en by he in e ac ion e ms (5.2), is diagonal on he s a es |αi,
and he single-pa icle unneling e m ˆ
H1in Eq. (5.3) is ea ed a second
71
o de . Fo a gi en s a e |αi,
Eγij −Eα=U+ (U+W)(mi−mj) + UNN∆mij
NN,(5.12)
wi h ∆mij
NN =Phkiimk−Phkijmk−1, whe e miindica es he pai occu-
pa ion numbe a si e ias deﬁned in Eq. (5.4). Fo U+W, UNN ≪U, he
denomina o s Eγij −Eαa e all o o de U, which leads o
ˆ
H(0)
e =ˆ
H0−2J2
UX
hijihˆmi( ˆmj+ 1) + ˆc†
iˆcji,(5.13)
whe e ˆciand ˆc†
ia e he pai des uc ion and c ea ion ope a o s such ha
ˆci|mii=mi|mi−1i
ˆc†
i|mii= (mi+ 1)|mi+ 1i.(5.14)
One can easily ob ain co ec ions o ˆ
H(0)
e by expanding (5.12) a highe o de s
in (U+W)/U and UNN/U bu , as we will see, he ze o h o de is al eady qui e
accu a e o desc ibe he physics o he sys em o he ange o pa ame e s
we conside .
5.2.1 G ound s a e insula ing phases and wo-pa icle
wo-hole exci a ions
We now make use o he eﬀec i e Hamil onian ˆ
H(0)
e de i ed abo e o s udy
he g ound s a e phase diag am o he sys em, s a ing om he insula ing
s a es. Fo e e y classical dis ibu ion o pai s in he la ice we can calcula e
he pai o de pa ame e ψi=hˆciiwi h he pe u ba i e mean-ﬁeld me hod
de i ed in Sec. 2.2.2, and ge he exp ession
ψi=2J2
U(mi+ 1)2
Ei
2P(J)+m2
i
Ei
2H(J)¯
ψi,(5.15)
whe e ¯
ψi=Phjiiψj, and he ene gy cos s o adding a pai (2P) and emo ing
a pai (2H) can be calcula ed wi h he diagonal pa o Eq. (5.13), and a e
espec i ely gi en by
Ei
2P(J) = −2µ+ 2Umi+ (2mi+ 1)W+ 2V1,i
dip −2J2
UX
hkii
(2mk+ 1)
Ei
2H(J) = 2µ−2U(mi−1) −(2mi−1)W−2V1,i
dip +2J2
UX
hkii
(2mk+ 1),
(5.16)
72
wi h V1,i
dip deﬁned jus a e Eq. (5.8). By imposing he posi i i y o hese
exci a ions we ge o he exp essions
U(mi−1) + (mi−1
2)W+V1,i
dip −J2
UX
hkii
(2mk+ 1) < µ <
Umi+ (mi+1
2)W+V1,i
dip −J2
UX
hkii
(2mk+ 1),(5.17)
which can be easily compa ed wi h Eqs. (5.9) o he single-pa icle single-
hole exci a ions. A J= 0, i is s aigh o wa d o no ice ha o any W < 0
he wo-pa icles wo-holes exci a ions o Eqs. (5.17) gi e mo e es ic i e
condi ions han hei co esponding single pa icle-hole exci a ions o Eqs.
(5.9). The e o e we conclude ha a J= 0, o any W < 0 he low-lying
exci a ions o a classical dis ibu ion o pai s in he la ice a e ob ained by
adding o emo ing a pai a any si e i, in ag eemen wi h he p e ious
s a emen s.
Using Eq. (5.15) and condi ions (5.17), one can calcula e he mean-ﬁeld
lobes o any gi en conﬁgu a ion o pai s in he la ice. The lobes o he
checke boa d and doubly occupied checke boa d a e shown in Fig. 5.4 o
he 0 h ( ull lines) and 1s o de (dashed lines) eﬀec i e Hamil onians. The
compa ison be ween he wo shows ha , o he pa ame e s conside ed he e,
he 0 h o de al eady cap u es he physics accu a ely. I is wo h no icing
ha he J2dependence o he ene gy o he elemen a y exci a ions is a he
o igin o he een an beha io o he lobes, which was p edic ed by exac
ma ix-p oduc -s a e calcula ions o he 1D geome y in [65]
5.3 Gu zwille mean- ield app oach and a-
lidi y o he low ene gy subspace
While he MI phases a e p edic able h ough he pe u ba i e mean-ﬁeld
app oach o Eqs. (5.15) o he pai o de pa ame e s, o iden i y he SF
phases, bo h P SF and PSS ou side o he lobes, i is necessa y o make
use o he imagina y ime e olu ion in oduced in Sec. 2.2.1 based on a
Gu zwille Ansa z o he pai wa e unc ion. The e o e, we need o calcu-
la e he dynamics equa ion equi alen o (2.6) in he low ene gy subspace
desc ibed by he eﬀec i e Hamil onian ˆ
H(0)
e o Eq. (5.13).
We begin wi h he ime dependen Gu zwille wa e unc ion o he pai s,
which is gi en by
|Φi=Y
iX
m
(i)
m|mii,(5.18)
73
Figu e 6.1: Schema ic ep esen a ion o a 2D op ical la ice popula ed wi h
dipola Bosons pola ized in bo h di ec ions pe pendicula o he la ice plane.
The pa icles eel epulsi e in a-species Uaa,Ubb, and in e -species Uab epul-
si e on-si e ene gies. The nea es -neighbo in e ac ion is epulsi e UNN >0
o aligned dipoles, while i is a ac i e −UNN o an i-aligned pa icles, and
he hopping e m Jis equal o bo h he species.
which a e simul aneously diagonal on a gi en Fock s a e |ν, mii. No ice
ha he eigen alues o hese wo ope a o s a e no independen . In ac by
ﬁxing ν, he eigen alues o ˆmican only assume 2ν+ 1 alues gi en by m=
{−ν, −ν+ 1, ..., +ν}, in comple e analogy wi h he spin angula momen um
ope a o ˆ
S2
iand i s p ojec ion along he zaxis ˆ
Sz
i, as we will discuss in
Sec. 6.2.2. I is use ul o in oduce he a e age magne iza ion o he sys em,
deﬁned as
M=1
NSX
i
mi,(6.3)
whe e NSis he o al numbe o la ice si es, because, as we will see, i is a
con enien quan i y o desc ibe he phases o he sys em.
We subs i u e Eqs. (6.2) in o Hamil onian (6.1), which we w i e as ˆ
H=
80

ˆ
Hν
0+ˆ
Hm
0+ˆ
Hνm
1, whe e he diﬀe en e ms ead
ˆ
Hν
0=X
ih−2µ+ˆνi+ 2Uˆνiˆνi−1
2i (6.4)
ˆ
Hm
0=X
ih−2µ−ˆmi+ 2UNN X
hjii
ˆmiˆmji(6.5)
ˆ
Hνm
1=−JX
hiji
[ˆa†
iˆaj+ˆ
b†
iˆ
bj].(6.6)
We ha e in oduced he chemical po en ials
µ±=µa±µb
2,(6.7)
which espec i ely ﬁx he eigen alues o he ﬁlling ac o and he imbalance
ope a o s (6.2). In he ollowing we will conside ˆ
Hνm
1 o be a small pe u -
ba ion on he in e ac ion e ms. In he limi o U≫UNN ,J, he g ound
s a e o he sys em is ound o be a uni o m dis ibu ion o cons an ﬁlling
ac o νi= ¯νa each si e o he la ice. The alue o ¯νis ﬁxed by µ+, and can
be in ege as well as semi-in ege . This is be e unde s ood a J= 0, whe e
we can calcula e he expec a ion alue o ˆ
Hν
0on a gi en classical dis ibu ion
o a oms in he la ice |Φi=Qi|νi, miii, as ollows
hΦ|ˆ
Hν
0|Φi=X
ih−2µ+νi+ 2Uνiνi−1
2i,(6.8)
whe e νi=hΦ|ˆνi|Φi. In he ene gy (6.8) each si e iis sel -simila , and like
in he homogeneous case o a Bose-Hubba d Hamil onian a J= 0, he
minimum o Eq. (6.8) is p o ided by a uni o m dis ibu ion νi= ¯νa each
si e o he la ice. Ins ead, o a gi en ¯ν, ﬁnding he magne iza ion which
minimize he expec a ion alue
hΦ|ˆ
Hm
0|Φi=X
ih−2µ−mi+ 2UNN X
hjii
mimji,(6.9)
is non- i ial due o he p esence o he nea es neighbo epulsion 2UNN,
whe e mi=hΦ|ˆmi|Φi. Howe e , we can quali a i ely a gue ha o |µ−| ≫
UNN , he minimum o he ene gy (6.9) is ob ained o mi=ν×signµ−,∀i,
which co esponds o a e omagne ic phase (FM) o a e age magne iza ion
M=ν×signµ−, whe e only pa icles o one species a e p esen . Ins ead, o
µ−= 0, a succession o nea es neighbo s mi=νand mj=−νp o ides he
minimum o Eq. (6.9), and he phase is an i- e omagne ic (AM), i.e. M= 0.
81
The spa ial dis ibu ion o he pa icles is gi en by si es occupied om he
species aal e na ed wi h si es occupied by he species bin a checke boa d-
like s uc u e. In Fig. 6.2 we plo he g ound s a e a J= 0, in he µ− s. µ+
plane, whe e he ex in pa en hesis (ν, M) indica e espec i ely he ﬁlling
ac o and he a e age magne iza ion.
−0.75 −0.5 −0.025 0 0.025 0.5 0.75
−0.0125
0.9625
1.9375
µ−/U
µ+/U
(0,0)
(1/2,0)
(1,0)
(3/2,0)
(1/2,1/2)(1/2,−1/2)
(1,1)(1,−1)
Figu e 6.2: G ound s a e o he eﬀec i e Hamil onian ˆ
H(0)
e (see ex ) a J= 0,
calcula ed o a 2×2 elemen a y cell sa is ying pe iodic bounda y condi ions.
The ex in pa en hesis (ν, M) indica e he ﬁlling ac o νand he a e age
magne iza ion M, espec i ely.
In he nex sec ion we will include he p esence o unneling. In his si ua ion
we will see ha he heo e ical desc ip ion o he sys em canno be based on
s anda d mean-ﬁeld heo y, which is no sui able o desc ibe he g ound s a e
o he sys em in a co ec way. In ac , he phase diag am o Fig. 6.2 has
been ob ained using an eﬀec i e Hamil onian, which desc ibes co ec ly he
physics o he sys em as we will explain in he ollowing sec ion.
6.2.2 Low-ene gy subspace and e ec i e Hamil onian
The g ound s a e o he sys em a J= 0 is desc ibed by a p oduc o e
single-si e Fock s a es o he ype
|αi=Y
i|ν, miii,(6.10)
82
wi h uni o m on-si e occupa ion ν. A single pa icle hopping changes he
o al on-si e popula ion a he si es in ol ed in he hopping p ocess, and
he e o e b eaks he ansla ional in a iance o he g ound s a e wi h espec
o he on-si e occupa ion ν. The ene gy cos o hese exci a ions is o he
o de o he on-si e in e ac ion ene gy U, and is he e o e e y cos ly in he
limi o U≫UNN , J. On he con a y, exchanging wo pa icles om nea es
neighbo ing si es o ﬂipping he di ec ion o one dipole ( om up o down o
ice e sa), does no equi e such a la ge amoun o ene gy. This deﬁnes
a low-ene gy subspace spanned by he |αidis ibu ions o cons an ﬁlling
ac o νo Eq. (6.10), which is ene ge ically well sepa a ed om he es o
he Hilbe space in he limi o pa ame e s we a e conside ing he e. Thus,
a success ul desc ip ion o such a sys em is ob ained h ough an eﬀec i e
Hamil onian ˆ
He es ic ed o he low-ene gy subspace, whe e single-pa icle
hopping is supp essed and unneling is included a second o de pe u ba-
ion heo y. The alidi y o he eﬀec i e Hamil onian elies on he exis ence
o his low-ene gy subspace well sepa a ed in ene gy om he subspace o
i ual exci a ions, o which i is coupled ia single-pa icle hopping. The
ele an i ual subspace is hen ob ained om he s a es |αi ia single pa -
icle hopping, namely
|γ(a)
ij i=ˆa†
iˆaj
qna
j(na
i+ 1)|αi
|γ(b)
ij i=ˆ
b†
iˆ
bj
qnb
j(nb
i+ 1)|αi,
(6.11)
as schema ically ep esen ed in Fig. 6.3.
This si ua ion is quali a i ely no diﬀe en han he one discussed in Sec.
5.2 in he con ex o he wo laye s, and he e o e we can apply he same
echnique o compu e ˆ
He . In he basis o cons an on-si e popula ion ν, he
ma ix elemen s o such a Hamil onian in second o de pe u ba ion heo y
a e gi en by
hα|ˆ
He |βi=hα|ˆ
H0|βi− 1
2X
γhα|ˆ
Hνm
1|γihγ|ˆ
Hνm
1|αi
×1
Eγ−Eα
+1
Eγ−Eβ(6.12)
whe e ˆ
H0=ˆ
Hν
0+ˆ
Hm
0, gi en by he sum o he in e ac ion e ms (6.4) and
(6.5), is diagonal on he s a es |αi, and he single-pa icle unneling e m
83
Figu e 6.3: Schema ic ep esen a ion o a wo-pa icle hopping be ween he
s a es |αiand |βi, belonging o he low-ene gy subspace a ν= 1/2. These
s a es a e coupled o he i ual exci a ion |γi h ough single-pa icle jumps.
ˆ
Hνm
1o Eq. (6.6) is ea ed a second o de . Fo a gi en s a e |αi,
Eγij −Eα=U+UNN ∆mij
NN,(6.13)
wi h ∆mij
NN =Phkii2mk−Phkij2mk−1, whe e miindica es he popula ion
imbalance a si e io Eq. (6.2). Fo U≫UNN , he denomina o s Eγij −Eα
a e all o o de U, which leads o
ˆ
H(0)
e =ˆ
Hν
0−2J2
UX
hiji
ˆνi(ˆνj+ 1)
+ˆ
Hm
0−2J2
UX
hijihˆmiˆmj+ ˆc†
iˆcji,(6.14)
whe e ˆci= ˆaiˆ
b†
iand ˆc†
i= ˆa†
iˆ
bia e he des uc ion and c ea ion ope a o s o a
composi e, made o a pa icle o one species and a hole o he o he species,
such ha
ˆci|ν, mii=pν(ν+ 1) −mi(mi−1)|ν, mi−1i
ˆc†
i|ν, mii=pν(ν+ 1) −mi(mi+ 1)|ν, mi+ 1i,(6.15)
and hei commu a ion ela ions a e gi en by [ˆci,ˆc†
j] = −2 ˆmiδij, wi h δij
being he K onecke del a.
Fo a gi en ν, he second line o he Hamil onian (6.14) can be equi a-
len ly w i en in e ms o he spin ope a o s a si e i, gi en by
ˆ
Si=1
2X
uu′
ˆa†
iu~σuu′ˆaiu′,(6.16)
84
whe e ~σ = (σx, σy, σz) a e he Pauli ma ices, and u=a, b indica es he
wo species. Thus, he c ea ion and annihila ion ope a o s (6.15) become
ˆc†
i=ˆ
Sx
i+iˆ
Sy
iand ˆci=ˆ
Sx
i−iˆ
Sy
i espec i ely, while he imbalance ope a o is
gi en by ˆmi=ˆ
Sz
ias al eady an icipa ed in Sec. 6.2.1. The e o e, in he spin
ope a o s language (6.16), he second line o he Hamil onian (6.14) looks
like a Heisenbe g spin Hamil onian (see e.g. [46]). The chemical po en ial µ−
plays he ole o an ex e nal magne ic ﬁeld along he zaxis, and he in e play
be ween µ−, and he nea es -neighbo in e ac ions de e mines he magne ic
o de ing o he sys em, as we will discuss in he nex sec ion.
6.3 Mean- ield
In his sec ion, we p o ide a mean-ﬁeld solu ion o he eﬀec i e Hamil onian
(6.14) in o de o in es iga e he phases o he sys em. We iden i y he
diﬀe en phases h ough he mean-ﬁeld composi e o de pa ame e s hˆcii, as
well as bo h he single-pa icle ones hˆaii, and hˆ
bii.
Fo e e y subspace a cons an ﬁlling, we ﬁnd ha he sys em p esen s
h ee diﬀe en kinds o phases. The Mo -insula ing phase (MI), wi h a well
deﬁned numbe o pa icles a each si e o he la ice, and absence o any
low-ene gy anspo [1]. The MI is cha ac e ized by anishing hˆcii=hˆaii=
hˆ
bii= 0, and depending on he alue o µ−p esen s ei he FM o AM o -
de ing. Ins ead in he supe -coun e - luid phase (SCF), cha ac e ized by
on-si e densi y ﬂuc ua ions, he ne anspo o a oms is s ill supp essed
bu a coun e ﬂow is p esen , in which he cu en s o he wo species a e
equal in absolu e alue bu in opposi e di ec ions [78]. In he SCF phase,
while he single-pa icle o de pa ame e s s ill anish hˆaii=hˆ
bii= 0, he
composi e o de pa ame e s a e non-ze o hˆcii 6= 0, indica ing he p esence o
coun e ﬂow. In his wo k we also ﬁnd e idences o he exis ence o he no el
coun e low supe solid phase (CSS). The CSS is cha ac e ized by anishing
single-pa icle o de pa ame e s hˆaii=hˆ
bii= 0, and non- anishing compos-
i e o de pa ame e s hˆcii 6= 0, coexis ing wi h b oken ansla ional symme y,
namely, a modula ion o bo h mi, and hˆciion a scale la ge han he one o
he la ice spacing, analogously o he supe solid phase.
To de e mine he insula ing phases we pe o m a pe u ba i e ea men
a ﬁ s o de in he composi e o de pa ame e s ψi=hˆcii, which allow us
o compu e he bounda ies o he insula ing lobes. Fu he mo e, we sol e
he ime dependen Gu zwille equa ions in imagina y ime o de e mine he
na u e o he SCF-CSS phases ou side he lobes.
85

6.3.1 Insula ing lobes
The low-ene gy subspace is spanned by he classical dis ibu ion o a oms
in he la ice |αio Eq. (6.10). Simila ly o he wo laye sys em discussed
in Sec. 5.2, in he limi o U≫UNN, asymp o ically all classical s a es |αi
become s able wi h espec o single-pa icle-hole exci a ions and de elop an
insula ing lobe a ﬁni e J. The ene gy o single pa icle-hole exci a ions is
o he o de o Ua J= 0 and is gi en by he wid h o he lobes a ﬁni e J
(see, e.g., he hin blue lobes in Fig. 6.4).
Ins ead, he low-lying exci a ions emain wi hin he subspace and a e
ob ained by adding (PH) o emo ing (HP) one composi e, made o a pa icle
o one species and a hole o he o he species, a he i- h si e o he la ice.
This co esponds o ﬂip he di ec ion o a dipole a he si e i, espec i ely om
down o up (PH) o om up o down (HP). Fo any gi en conﬁgu a ion |αi,
one can calcula e hei ene gy cos s wi h he diagonal e ms o he eﬀec i e
Hamil onian (6.14), which a e espec i ely gi en by
Ei
PH(J) = −2µ−+ 4 UNN −J2/UX
hkii
mk,
Ei
HP(J) = 2µ−−4UNN −J2/UX
hkii
mk.(6.17)
No ice ha in he las exp essions, he e is no explici dependence on he
chemical po en ial µ+. This is because by adding o emo ing one composi e,
one emains wi hin he subspace a ﬁlling ac o ν, and he e o e he con i-
bu ion o µ+ anish in he calcula ion o he exci a ions (6.17). By using he
pe u ba i e mean-ﬁeld me hod de i ed in Sec. 2.2.2, we can calcula e he
o de pa ame e s ψi=hˆcii o |αi, which sa is y he equa ions
ψi=2J2
Uν(ν+ 1) −mi(mi+ 1)
Ei
PH(J)+ν(ν+ 1) −mi(mi−1)
Ei
HP(J)¯
ψi,(6.18)
whe e ¯
ψi=Phjiiψj. Wi h Eqs. (6.18) one can calcula e he mean-ﬁeld
lobes o any dis ibu ion o a oms in he la ice |αi, p o iding he elemen a y
exci a ions (6.17), a e posi i e in some ange o he pa ame e s. As we ha e
demons a ed in Sec. 2.2.2 his is a necessa y condi ion o he exis ence o
an insula ing lobe and p o ides i s bounda ies a J= 0, gi en by
2UNN X
hkii
mk< µ−<2UNN X
hkii
mk.(6.19)
To ob ain he las inequali ies, one has o be ca e ul and ﬂip he di ec ion
o a dipole only whe e i is possible. Fo example, suppose he si e iis oc-
cupied only by one pa icle o he species a, i.e mi= 1/2, hen o his si e
86
he condi ions (6.19) educe o µ−>2UNN Phkiimk, since he only possible
exci a ion a his si e co esponds o emo e a composi e. As usual, o ﬁnd
he bounda ies o he insula ing lobes we use he p ocedu e explained in Sec.
2.2.2. Fo each si e o he la ice one has such a condi ion (6.19), and an Eq.
(6.18). The la e is a se o coupled equa ions, which can be w i en in he
ma ix o m M(µ−, U, UNN, J)·~
ψ= 0, wi h ~
ψ≡(···ψi···) being he ec o
o he o de pa ame e s a each si e o he la ice. Fo e e y µ−, a non- i ial
solu ion is p o ided by he smalles J o which de [M(µ−, U, UNN, J)] = 0,
ha is he insula ing lobe o he |αiconﬁgu a ion in he J s. µ−plane. In
Fig. 6.5 we plo he g ound s a e insula ing lobes calcula ed in his way o
ν= 1/2 (le ) and ν= 1 ( igh ). Fo all ﬁlling ac o s ν, we ﬁnd an AM
g ound s a e (ν, M = 0), which p esen s a spa ial dis ibu ion o al e na -
ing si es occupied by pa icles o species aand b esembling a checke boa d
s uc u e. Ins ead, by inc easing he absolu e alue o µ−we ﬁnd a FM
g ound s a e (ν, M =±ν), in which only pa icles o one ype a e p esen .
I is wo h no icing ha he insula ing lobes calcula ed in his way, do no
con ain any dependence on µ+, which does no en e in o Eqs. (6.18) as p e-
iously discussed. The e o e, o ob ain he comple e 3D phase diag am, one
has o compa e he ene gies o he g ound s a e conﬁgu a ions a diﬀe en
ν. Using he eﬀec i e Hamil onian (6.14), o any alue o J,µ−, and µ+, we
calcula e he ene gies o he g ound s a e conﬁgu a ions o diﬀe en ν, and
selec he s a e wi h he smalle ene gy. In his way we ha e ob ained o
example he phase diag am a J= 0 o Fig. 6.2.
6.3.2 Coun e low supe luid-supe solid
In he low-ene gy subspace a cons an ν, he Gu zwille Ansa z on he wa e
unc ion o he sys em eads
|Φi=Y
i
ν
X
m=−ν
(i)
ν,m|ν, mii,(6.20)
whe e we allow he Gu zwille ampli udes o dependen on ime (i)
ν,m( ). As
in Sec. 2.2.1, we ob ain he equa ions o mo ion o he ampli udes by min-
imizing he ac ion o he sys em, gi en by S=Rd L, wi h espec o he
a ia ional pa ame e s (i)
ν,m( ) and hei complex conjuga es ∗(i)
ν,m( ), whe e
L=hΦ|˙
Φi−h˙
Φ|Φi
2i−hΦ|ˆ
H(0)
e |Φi,(6.21)
87
0 0.02 0.04 0.06 0.08 0.1
−0.8
−0.6
−0.4
−0.2
0
0.2
0.4
0.6
0.8
J/U
µ−/U
Figu e 6.4: Insula ing lobes a ν= 1/2 ﬁlling, wi h UNN =U/160. The hick
line is he an i- e omagne ic insula ing s a e calcula ed wi h he eﬀec i e
Hamil onian ˆ
H(0)
e . The hin blue lines (dashed lines) ep esen insula ing
lobes calcula ed wi h espec o single-pa icle-hole exci a ions o he species
a(species b), ske ched as a plain do (squa e), and a e compu ed o µ+=
0.65U.
is he Lag angian o he sys em in he quan um s a e |Φi[60]. The e o e by
equa ing o ze o he a ia ion o he ac ion wi h espec o ∗(i)
ν,m, leads o he
equa ions
i~d
d (i)
ν,m=h−2µ−+ 4(UNN −J2/U)X
hjiihˆmjiim (i)
ν,m
−2J2
Uh¯
ψipν(ν+ 1) −m(m−1) (i)
ν,m−1
+¯
ψ∗
ipν(ν+ 1) −m(m+ 1) (i)
ν,m+1i,(6.22)
whe e hˆmii=Pν
m=−νm| (i)
ν,m|2, he ﬁelds ¯
ψi=Phjiiψjand Phjiihˆmji, ha e
o be calcula ed in a sel consis en way as explained in Sec. 2.2.1, and he
o de pa ame e is gi en by
ψi=hΦ|ˆci|Φi=
ν
X
m=−νpν(ν+ 1) −m(m+ 1) ∗(i)
ν,m (i)
ν,m+1.(6.23)
88
We sol e Eqs. (6.22) in imagina y ime τ=i , which due o dissipa ion is
supposed o con e ge o he g ound s a e. In Fig. 6.5 we show he g ound
s a e phase diag am o he sys em o ν= 1/2 (le ) and ν= 1 ( igh ),
compu ed in his way o UNN =U/160. In he egion ou side he insula ing
AM lobes and enclosed be ween he FM s a es, depending on he alues o
Jand µ−we ﬁnd ei he SCF o CSS. The CSS phase is cha ac e ized by
anishing single-pa icle o de pa ame e s hˆaii=hˆ
bii= 0, coexis ing wi h
a spa ial modula ion o he composi e o de pa ame e s hˆcii 6= 0, indica ing
he p esence o coun e ﬂow. The shaded a eas in Fig. 6.5 indica e whe e
hˆcii 6= 0, and p esen a spa ial modula ion. These ha e o be compa ed wi h
he egion whe e he single-pa icle o de pa ame e s a e ze o, in o de o
de e mine he limi s o alidi y o he CSS phase. Fo ν= 1/2, he hin blue
lines (dashed lines) in Fig. 6.4 ep esen he insula ing lobes calcula ed wi h
espec o single-pa icle-hole exci a ions o he species a(species b). The wo
lobes ex ending up o J∼0.09U, co espond o he AM insula ing o de ing,
hey delimi he ex ension o he CSS phase, and hey gi e an es ima e o
he limi s o alidi y o ˆ
H(0)
e , beyond which he subspace o cons an νlooses
i s meaning.
We ha e al eady men ioned, ha he bounda ies o he lobes calcula ed wi h
he eﬀec i e Hamil onian do no show any dependence on he chemical po en-
ial µ+, which does no gi e any con ibu ion in he exp ession o he low-lying
exci a ions (6.17). This is no ue in he case o he single-pa icle-hole insu-
la ing lobes, since adding o emo ing a single pa icle esul s in a change o
bo h µ+, and µ−. This makes he p ocess o es ima ing he limi s o alidi y
o ˆ
H(0)
e mo e complica ed and leads o a 3D phase diag am in he J,µ+, and
µ− a iables highly non- i ial, a sys ema ic s udy o which is s ill on going.
89
expanding he second exponen in he igh -hand-side o Eq. (8.2) one ge s
Z=X
α
e−βEαhα|ˆ
11−Zβ
0
dτˆ
H1(τ) +
+∞
X
m=2
(−1)mZβ
0
dτm.. Zτ2
0
dτ1ˆ
H1(τm).. ˆ
H1(τ1)|αi,(8.3)
whe e he in eg als a e o de ed in ime and he sum o e he s a es |αicomes
om he ace. Now we explici ly make use o he comple eness p ope y o
he {|αi} base, and inse m−1 iden i y ope a o s ˆ
11 = Pα|αihα|be ween
he p oduc s o ˆ
H1(τm) ope a o s, he e o e we can w i e
hα|ˆ
H1(τm).. ˆ
H1(τ1)|αi=X
α1,..,αm−1
Hααm−1
1(τm).. Hα2α1
1(τ2)Hα1α
1(τ1),(8.4)
whe e he ma ix elemen s
Hα′α
1(τ) = eτEα′Hα′α
1e−τEα=hα′|ˆ
H1|αie−τ(Eα−Eα′),(8.5)
con ain bo h diagonal (Eα) and oﬀ-diagonal (Hα′α
1) ma ix elemen s. We
now inse he las equa ion in o exp ession (8.3) and ge he ﬁnal exp ession
o he pa i ion unc ion
Z=X
α
e−βEα1−Zβ
0
dτHαα
1(τ)+ (8.6)
+∞
X
m=2
(−1)mZβ
0
dτm.. Zτ2
0
dτ1X
α1,..,αm−1
Hααm−1
1(τm).. Hα1α
1(τ1)),
which con ains only ma ix elemen s o he ope a o s ˆ
H0and ˆ
H1. The e o e,
by using his o malism o pa h in eg als, he calcula ion o he pa i ion
unc ion educes o a classical p oblem since only scala s en e in o Eq. (8.6),
bu we ha e payed he p ice o he ex a dimension τ. In o he wo ds, he
o iginal d-dimensional quan um sys em is equi alen o a (d+1)-dimensional
classical sys em.
I is wo h no icing ha since he pa i ion unc ion is a ace, pe iodic
bounda y condi ions in he imagina y ime τmus apply. This is easily
unde s ood by looking a he m- h o de e m o Z, which con ains he
p oduc o mma ix elemen s Hααm−1
1(τm).. Hα1α
1(τ1) ha a e o de ed in
ime om he ﬁ s a τ1, o he las a τm. The e o e, o any gi en αin
he ace, he ﬁ s ma ix elemen b ings α o some α1in he ime τ1≥0,
96

while he las ma ix elemen b ings αm−1back o αin he ime τm≤β.
All he possible conﬁgu a ions which a e pe iodic in imagina y ime and ha
en e in o he exp ession o he pa i ion unc ion Eq. (8.6), deﬁne he
conﬁgu a ion space spanned by a PIMC algo i hm.
8.1.1 Pa h In eg al Mon e Ca lo and he 2D ex ended
Bose-Hubba d model
We now conside a 2D sys em o L×Lsi es ﬁlled wi h pola ized dipola
Bosons, we assume spa ial pe iodic bounda y condi ions and he dipoles o
be pola ized pe pendicula ly o he 2D plane as explained in Chap e 3. The
sys em is he e o e desc ibed by he ex ended Bose-Hubba d Hamil onian
(3.1), which, o be consis en wi h he no a ions in ou publica ion [86], we
ew i e in his o m
ˆ
H=−JX
hijihˆ
b†
iˆ
bj+ˆ
biˆ
b†
ji+X
iU
2ˆni(ˆni−1) −µiˆni+VX
i<j
ˆniˆnj
3
ij
,(8.7)
whe e ˆ
b†
i(ˆ
bi) is he boson c ea ion (annihila ion) ope a o a si e i, ˆni=ˆ
b†
iˆ
biis
he numbe ope a o , V=D/a3>0 is he dipole-dipole in e ac ion s eng h
Ddi ided by he la ice spacing a, ij =|i−j|is he dis ance be ween wo
si es o he la ice, and µi=µ−Ωi2con ains he chemical po en ial µwhich
ﬁxes he numbe o pa icles and he cu a u e Ω o an ex e nal ha monic
conﬁnemen .
We choose o wo k in he basis o he in e ac ion e m o he Hamil onian
(8.7), i.e. Fock s a es |αi=QL2
i|niiio localized pa icles in he L×Lsqua e
la ice, whe e niis he occupa ion numbe a si e i. The e o e in his basis,
he diagonal ma ix elemen s en e ing Eq. (8.5) ake he o m
Eα=U
2X
i
ni(ni−1) −X
i
µini+VX
i<j
ninj
3
ij
,(8.8)
while he oﬀ-diagonal ones, a e gi en by he exp ession
−Hα′α
1= 2Jhα′|ˆ
b†
iˆ
bj|αi= 2Jq(nα
i+ 1)nα
j(8.9)
and hey connec s a es |α′iand |αi ha diﬀe only in he occupa ion numbe
o he wo nea es neighbo ing si es iand j, namely |α′i ≡ ˆ
b†
iˆ
bj
√(nα
i+1)nα
j|αi
wi h nα
ibeing he numbe o pa icles a he i- h si e o he s a e |αi. I
is impo an o no ice ha we do no use any cu oﬀ in he ange o he
dipole-dipole in e ac ion en e ing Eq. (8.8).
97
To w i e he pa i ion unc ion o he 2D ex ended Bose-Hubba d model,
we no ice ha he ﬁ s o de e m anishes since he ma ix elemen s (8.9)
a e oﬀ-diagonal, i.e. Hαα
1= 0, and due o he geome y o he sys em (2D
squa e la ice) i is no diﬃcul o see ha all he e ms wi h an odd alue
o malso anish. The e o e by ea anging he exponen ials and enaming
α≡α0, we ge o he exp ession
ZeBH =X
α0
e−βEα0+∞
X
m=2
(−2J)mAm×(8.10)
×Zβ
0
dτm.. Zτ2
0
dτ1X
α0,α1,..,αm−1
exp (−βEα0−
m−1
X
p=0
Eαp(τp+1 −τp)),
whe e Amis a p oduc o msqua e oo ac o s coming om Eq. (8.9) and we
ha e in oduced τ0=τm o compac he no a ion. We can compac u he
he no a ion by no icing ha o m= 0 wo hings happens: (i) he sum in
he exponen o Eq. (8.10) does no make any sense, since i is he e m o
o de m≥2 in he Taylo expansion, he e o e we deﬁne i o be ze o, and (ii)
in he sum o e he α-s only α0su i es. Keeping hese wo conside a ions
in mind and deﬁning Am=0 = 1, we hen w i e he pa i ion unc ion in he
compac o m
ZeBH =∞
X
m=0 X
α0,α1,..,αm−1
(−2J)mAm×(8.11)
×Zβ
0
dτm.. Zτ2
0
dτ1exp (−βEα0−
m−1
X
p=0
Eαp(τp+1 −τp)).
F om he las exp ession one can o mally w i e
ZeBH =X
ν
Wν,(8.12)
wi h Wνbeing he weigh o each conﬁgu a ion ν≡[m, α0(τ), α1(τ), .., αm−1(τ)],
whe e no only he α-s deﬁne νbu also hei dis ibu ion in he imagina y
ime. This is be e unde s ood om Fig. 8.1, whe e we ske ch one o such
conﬁgu a ions o he 2D la ice.
The imagina y ime, τ, is on he ho izon al axis, while on he e ical axis
he e a e all he si es o he la ice. Each line is called a wo ldline and i
ep esen s a numbe o pa icles p opo ional o he wid h o he line: he
98
0β
i−1
i
i+ 1
.
.
.
.
.
.
Figu e 8.1: Schema ic ep esen a ion o one conﬁgu a ion which en e s he
calcula ion o ZeBH. Each wo ldline ep esen a numbe o pa icles p opo -
ional o he wid h o he line, whe e he dashed black line is o npa icles
while he solid and bold blue lines ha e occupa ion numbe s equal o n+ 1
and n+ 2. Wo ldlines ha e o ulﬁll pe iodic bounda y condi ions in he
imagina y ime τ, and he e ical a ows in co espondence o changes in
he occupa ion numbe s a e called kinks.
dashed black line is o npa icles while he solid and bold blue lines ha e
occupa ion numbe s equal o n+ 1 and n+ 2 espec i ely. Because he
pa i ion unc ion is a ace, wo ldlines ha e o close on hemsel es and since
we ha e also assumed spa ial pe iodic bounda y condi ions one can imagine
he conﬁgu a ion o Fig. 8.1 o be w apped on a o us. We call he phase
space o all possible conﬁgu a ions he closed pa h con igu a ion space (CP),
which is spanned by he PIMC algo i hm.
I we cu one conﬁgu a ion a a ce ain ins an in imagina y ime, we
ge he quan um sys em in a pa icula s a e, and he poin s in imagina y
ime whe e he sys em changes s a e a e called kinks, which in Fig. 8.1 a e
ep esen ed by e ical a ows. A conﬁgu a ion wi h a numbe o kinks equal
o m, con ibu es o he m- h o de e m o he pa i ion unc ion Eq. (8.11),
and i is s aigh o wa d o see ha he e exis an inﬁni e numbe o diﬀe en
conﬁgu a ions wi h he same numbe o kinks, he diﬀe ence being he ime
a which he kinks ake place and/o he diﬀe en s a es hey connec . The
upda ing p ocedu e o a PIMC algo i hm he e o e consis s o changing he
numbe o kinks and/o hei posi ion in imagina y ime. We will discuss he
99
upda ing p ocedu e speciﬁcally o he Wo m algo i hm in he nex sec ion.
8.2 The Wo m algo i hm
The Wo m Algo i hm, which was o iginally de eloped by P oko ’e , S is-
uno and Tupi syn [81, 82], wo ks in an enla ged conﬁgu a ion space, in
which one allows one disconnec ed wo ldline, he wo m, d awn as a ed line
in Fig. 8.2. This is equi alen o wo k in he G and-Canonical ensemble, as
we shall discuss in Sec. 8.2.1, and he disconnec ed wo ldline allow o eﬃ-
cien ly collec s a is ics o calcula ing he Ma suba a G een unc ion, deﬁned
as
G(j, τ) = hˆ
Tτˆ
bi+j(τ0+τ)ˆ
b†
i(τ0)i,(8.13)
whe e ˆ
Tτis he ime-o de ing ope a o , τ0and τa e wo poin s in imagina y
ime, iand ja e wo si es o he la ice, and he symbol h.is ands o he
s a is ical a e age o he expec a ion alue o an ope a o . Due o space and
imagina y ime ansla ional in a iance o he sys em, he G een unc ion Eq.
(8.13) does no depend on iand τ0. The conﬁgu a ion space o he Ma suba a
G een unc ion is called he CPgspace, and i is easy o see ha he only
diﬀe ence be ween conﬁgu a ions con ibu ing o he pa i ion unc ion ZeBH
and hose con ibu ing o he G een unc ion Gis ha , o he la e , one o
he wo ldlines s a s a (i, τ0) and ends a (i+j, τ0+τ), i.e. he wo ldline is
disconnec ed.
8.2.1 Upda ing p ocedu es
Le us now discuss he upda ing p ocedu e o he Wo m Algo i hm, ha
is when he sys em is in a ce ain conﬁgu a ion νand he algo i hm has
o gene a e andomly a new conﬁgu a ion ν′ o collec s a is ics o e alu-
a ing he obse ables o in e es . No ice ha in o de o ensu e e godici y,
and he e o e he eliabili y o he s a is ics, he upda ing p ocess mus be
ully andom such as o co e enough o he phase space o he sys em in a
easonable amoun o ime.
Apa om he c ea ion o a wo m, which is done in he CP space, all
o he upda es a e done in he CPgspace h ough he wo ends o he wo m.
One can pic u e he upda ing scheme as sequence o ’d awing’ and ’e asing’
p ocedu es, happening a he end poin s o he wo m. Gi en he conﬁgu-
a ion ν, he algo i hm selec s andomly an in e al o he conﬁgu a ion in
imagina y ime, which we call a ime-in e al. Below we lis and desc ibe he
ou ypes o upda es he Wo m Algo i hm goes h ough.
100
0β
i−1
i
i+ 1
.
.
.
.
.
.
Figu e 8.2: Conﬁgu a ion o he CPgspace, he ed disconnec ed line ep e-
sen s he wo m.
C ea ion o a wo m
C ea ing a wo m is he only upda e pe o med in he CP space, he e o e he
s a ing poin is a conﬁgu a ion νbelonging o CP. One o he wo ldlines o
νis andomly selec ed, and on ha wo ldline a ime in e al is also selec ed a
andom, e.g. he in e al n1in Fig. 8.3 delimi ed by τmin and τmax indica ed
by he c osses. Then he p og am akes andomly wo poin s on he segmen
n1, say τ1and τ2, which will be he wo m ex emi ies indica ed by plain do s
in Fig. 8.3, and he condi ion τmin < τ1< τ2< τmax has o be sa isﬁed.
Wi h equal p obabili y one sugges s o d aw a piece o wo ldline o dele e
a piece om an exis ing wo ldline, wi h he cons ain s ha he esul ing
conﬁgu a ion belongs o he Hilbe space, i.e. i is no possible o e ase
om an emp y in e al o o d aw on an in e al which has eached he
maximum occupa ion numbe allowed, i any. The wo m is he e o e c ea ed
and all o he upda es will ake place h ough i s wo ex emi ies.
Dele ion o a wo m
In analogy, he opposi e upda ing p ocess which is he dele ion o a wo m,
can only ake place in he CPgspace and only i he wo ex emi ies o he
wo m belong o he same wo ldline so as o assu e ha a e ha ing dele ed
a wo m he sys em is in a CP s a e.
101

τmin τmax
n1
τ1τ2
ˆ
bˆ
b†
τ2
τ1
ˆ
b†ˆ
b
Figu e 8.3: C ea ion o a wo m. F om a gi en conﬁgu a ion, one wo ldline is
andomly chosen ( op), in which a ime in e al delimi ed by τmin and τmax
is andomly selec ed. Then wi hin he in e al, wo poin s τ1and τ2a e also
chosen andomly and will be he wo ex emi ies o he wo m. Wi h equal
p obabili y one can choose o dele e a piece o wo ldline (bo om le ) o d aw
a piece o wo ldline (bo om igh ) and he wo m is he e o e c ea ed.
Time shi
This is he simples o he upda es and i consis s o mo ing one o he
ex emi ies o he wo m in a andom poin o he imagina y ime, such as o
leng hen o sho en he size o he wo m. The p og am selec s andomly a
ime in e al, delimi ed by τmin and τmax, and in his in e al only one poin
is andomly chosen, which is he poin whe e ei he he head o he ail o
he wo m is mo ed o.
Space shi
This upda e changes he numbe o kinks and i consis s o c ea ing o de-
la ing a kink o he le (space shi le ) o o he igh (space shi igh )
o he ope a o ˆ
b(o ˆ
b†). Fig. 8.4(a) shows he c ea ion o a kink backwa d
in imagina y ime, i.e. he space shi le . Two neighbo ing wo ldlines a e
selec ed a andom as o example iand jo Fig. 8.4(a), hen based on
he cu en posi ion o he ope a o ˆ
b he p og am chooses andomly a ime
in e al delimi ed by τmin and τmax. Wi hin his in e al he p og am selec s
andomly a poin whe e o c ea e o dele e a kink, as o example shown in
Fig. 8.4(a) o he c ea ion p ocess, wi h he equi emen ha he c ea ed
o dele ed kink does no in e e e wi h any o he kinks.
The las upda e, he space shi igh shown in Fig. 8.4(b), is equi alen o
he le one wi h he only diﬀe ence ha he kink is inse ed o dele ed o
102
ˆ
b
τmin
i
j
ˆ
bτmax
τmin
i
j
ˆ
b
τmax
τmin
i
j
(a)
(b)
τmax
τmin
i
j
τmax
ˆ
b
Figu e 8.4: Ske ch o space shi upda es ha c ea e o dele e kinks. In he
space shi le (a), one kink is c ea ed backwa d in he imagina y ime and
he e o e o he le o he ope a o ˆ
b, while in he space shi igh he kink
is inse ed o he igh o ˆ
b, i.e. o wa d in he imagina y ime.
he igh side o ˆ
bope a o , i.e. o wa d in imagina y ime.
These a e all he upda es pe o med by he WA. F om hese, i is s aigh -
o wa d o see ha he WA wo ks in he G and Canonical ensemble whe e
he chemical po en ial becomes an inpu pa ame e which ﬁxes he a e age
pa icle numbe . Fo example, suppose he algo i hm s a s wi h an ini ial
conﬁgu a ion νo ze o pa icles in he sys em, i.e. he analogous ep esen-
a ion o Fig. 8.1 would be a bunch o ho izon al dashed lines, om his
conﬁgu a ion he only possible upda e is o c ea e a wo m wi h one pa icle
in a gi en wo ldline, and ough he space shi and ime shi upda es he
pa icle will he e o e mo e in he si es o he la ice.
103
8.2.2 Ad an ages o he Wo m algo i hm
The upda es desc ibed abo e a e all local and allow o d aw/e ase any line,
and jump be ween he si es. Al hough only conﬁgu a ions belonging o he
CP space con ibu e o he e alua ion o he pa i ion unc ion, by using
he enla ged conﬁgu a ion space CP +CPg he in e media e conﬁgu a ions
wi h one disconnec ed loop allow o eﬃcien ly collec s a is ics o he G een
unc ion. Fo an algo i hm wo king in he CP space only, ins ead, collec ing
s a is ics o he G een unc ion esul s compu a ionally e y expensi e.
Ano he ad an age o he WA is ha i does no suﬀe om c i ical
slowing down in he icini y o a c i ical poin . In he c i ical egion, a sys em
de elops long ange co ela ions, and in mos cases an algo i hm based on
local upda es esul s e y ineﬃcien in simula ing such a sys em o which
he ele an deg ees o eedom a e non-local, and i esul s in he di e gence
o he au oco ela ion ime wi h he sys em size. Al hough he WA pe o ms
local upda es, i o e comes his p oblem by using he d awing and e asing
upda ing p ocedu es h ough he wo m ends, which a e di ec ly linked o
he c i ical modes (long ange o de in G(j, τ)). This u ns ou o be e y
eﬃcien in gene a ing independen conﬁgu a ions also in he c i ical egion.
The WA is also eﬃcien in sampling opologically diﬀe en conﬁgu a ions
and conﬁgu a ions which a e sepa a ed by an ene gy ba ie , which is a nec-
essa y condi ion in o de o main ain e godici y. An example o wo opolog-
ically diﬀe en conﬁgu a ions is shown in Fig. 8.5, whe e a one-dimensional
sys em wi h one pa icle (wo ldline) is conside ed. Pe iodic bounda y con-
di ions in ime and space apply, i.e. he sys em is a o us whe e he bo om
and op ace s o he cylinde a e glued oge he . Fig. 8.5(a) ep esen s a
conﬁgu a ion wi h ze o winding numbe s, i.e. he wo ldline does no ‘wind’
in imagina y ime. Fig. 8.5(d), ins ead, ep esen s a conﬁgu a ion wi h one
winding numbe , i.e. he wo ldline winds once in imagina y ime. An al-
go i hm based on local upda es which only wo ks in he CP space would
no allow o sample conﬁgu a ions wi h diﬀe en winding numbe s, unless a
global upda e which in oduces a winding numbe a once, is in oduced. The
WA, ins ead, can easily go om conﬁgu a ion o Fig. 8.5(a) o conﬁgu a ion
Fig. 8.5(d) (see a ske ch in Fig. 8.5(b)-(c)).
Being able o sample conﬁgu a ions wi h diﬀe en winding numbe s is
c ucial in o de o simula e SF sys ems. I was shown in [84], ha he
supe ﬂuid s iﬀness can be ex ac ed om he s a is ics o winding numbe s:
ρs=ThW2i
dLd−2,(8.14)
whe e Tis he empe a u e, L he sys em size, d he dimensionali y, and
W2=Pd
i=1 W2
i.
104
i
m
a
g
i n
a
y
i
m
e
a ) b ) c )
d )
Figu e 8.5: One-dimensional sys em wi h (a) ze o and (d) one winding num-
be (s). (b)-(c) ske ch on how he WA is able o go om (a) o (d).
105
Figu e 9.4: Spa ial densi y p oﬁle in 2D o N≃1000 pa icles in a ha monic
po en ial. Phases a e indica ed. (a-b) V/J = 15, µ/J = 55, Ω/J = 0.05 and
T/J = 0.0377; (c) V/J = 5, µ/J = 19, Ω/J = 0.01 and T/J = 0.1; (d)
V/J = 20, µ/J = 51, Ω/J = 0.04 and T/J = 0.25.
allowing o a clea de e mina ion o his phase. Panel (d) shows a diso de ed
ST-phase a he cen e , su ounded by an ex ended Mo -shell wi h ρ= 1/4.
The diso de in his case is a esul o bo h ﬁni e T/J = 0.25 and he ac
ha he ST solid is less obus owa ds quan um and classical ﬂuc ua ions
compa ed o he CB and SR ones.
These exac esul s o Ω 6= 0 conﬁ m ha he phase-diag am Fig. 9.1 is
he key o p edic and in e p e expe imen al obse ables, assuming a local
densi y app oxima ion.
112

Chap e 10
Conclusions
In his hesis we ha e s udied Bose-Hubba d (BH) models wi h dipola in e -
ac ions. In [57, 58] we ha e s udied a single componen gas o dipola Bosons
in a wo-dimensional (2D) op ical la ice, whe e he dipoles a e pola ized
pe pendicula ly o he 2D plane, esul ing in an iso opic epulsi e in e -
ac ion. Wi h a mean-ﬁeld app oxima ion based on a Gu zwille Ansa z, we
ha e shown ha such a sys em possesses many almos degene a e me as able
s a es, simila ly o a diso de ed sys em, and we ha e shown ha he dipole-
dipole in e ac ion is esponsible o he appea ance o hese s a es. We ha e
s udied in de ail he a e o hese me as able s a es, showing how hey can be
p epa ed on demand, we ha e discussed cu en expe imen al echniques ha
may be able o de ec he me as able s a es, and we ha e calcula ed hei
li e ime due o unneling. We ha e s udied he g ound s a e o he sys em
and ound he p esence o insula ing checke boa d-like s a es wi h ac ional
ﬁlling ac o s ν, and by using Quan um Mon e Ca lo me hods ( he Wo m al-
go i hm) in [86] we ha e conﬁ med his p edic ion. Mo eo e , we ha e ound
e idences o a De il’ s s ai case in he g ound s a e, which was p e iously
p edic ed in he phase diag am o a one-dimensional geome y [13] bu no
in 2D. In e es ingly, we ha e ound supe solid egions in he g ound s a e,
ob ained by doping he solids ei he wi h pa icles o acancies.
In [74], we ha e s udied a sys em composed o wo 2D laye s in which
he dipoles a e pola ized pe pendicula ly o he planes. The dipola in e -
ac ion is he e o e epulsi e o pa icles laying on he same plane, while i
is a ac i e o pa icles a he same la ice si e on diﬀe en laye s. Ou
mean-ﬁeld calcula ions ha e shown ha pa icles pai in o composi es, and
we ha e demons a e he exis ence o he no el Pai Supe Solid (PSS) quan-
um phase. I would be in e es ing o e i y he ex ension o he PSS phase
by conside ing mo e nea es neighbo s in he dipola in e ac ions, bo h wi h
mean-ﬁeld and Quan um Mon e Ca lo me hods.
113
Ano he in e es ing di ec ion o in es iga ion may be o s udy how he
p e ious si ua ion changes by adding mo e laye s. Fo example, in a h ee
laye s s uc u e, i could be in e es ing o e i y whe he he pai ing beha io
s ill su i es o i he sys em is domina ed by a diﬀe en physics, e.g. like
he o ma ion o iple s.
The sys em we a e cu en ly s udying, namely a 2D la ice whe e he
dipoles a e ee o poin in bo h di ec ions pe pendicula ly o he plane, has
shown many encou aging esul s. Wi hin a mean-ﬁeld app oach, we ha e
ound magne ic phases whe e he g ound s a e p esen s e omagne ic o
an i- e omagne ic o de ing. Fo la ge enough unneling, we ha e ound e i-
dence o he exis ence o a Coun e ﬂow Supe Solid (CSS) quan um phase.
Based on cu en expe imen al echniques wi h OH molecules, we ha e p o-
posed a possible ealiza ion o his sys em in he labo a o y. Ou heo e ical
s udy on his sys em is s ill on-going, and de ailed mean-ﬁeld calcula ions
p o iding he alidi y o he CSS phase ha e o be done. As well as o he
wo-laye s sys em desc ibed abo e, inc easing he numbe o nea es neigh-
bo s in he dipola in e ac ions could be a possible di ec ion o in es iga ion.
Finally, changing he geome y o he op ical la ice could lead o in e es ing
phenomena, as us a ion, which may happen in a iangula la ice.
As p e iously s a ed, ou p edic ions ha e di ec expe imen al conse-
quences, and we hope ha hey will be soon checked in expe imen s wi h
ul acold dipola a omic and molecula gases.
114
Acknowledgmen s
I is a g ea pleasu e o me o hank all he people in he quan um op ics
heo y g oup o ICFO, p esen and pas membe s. I spen ou beau i ul
yea s o my li e, I ha e lea ned, and I ha e changed, and i is om my
hea ha I hank all he g oup. I especially hank Chia a Meno i, and
Maciek Lewens ein o hei cons an suppo , and pa ience, hey showed
me du ing hese yea s, he esul o which is no only in his hesis wo k. I
hank Ba ba a Capog osso-Sansone and Guido Pupillo, o hei guide in he
Quan um Mon e Ca lo wo k, and Pe e Zolle o he kind hospi ali y while
I was in Innsb uck.
This hesis was pa ially w i en a he Indian Associa ion o he Cul-
i a ion o Science, in Calcu a, while I was isi ing K ishnendu Sengup a.
The wo m hospi ali y he showed me is unique, and I hank him oge he wi h
all he people o he heo e ical physics depa men .
Ou side o ICFO he e a e many people who helped me du ing hese
yea s. A special hanks goes o Ba ba a, who suppo ed me all he ime, and
knows almos all he people in he g oup wi hou ha ing me hem. To my
iend An oni, wi h whom I ha e sha ed e y unny si ua ions, and las , bu
no leas , o my amily.
115
Appendix A
Spec um o exci a ions
The low-lying exci a ion ene gies a e c ea ing pa icles (P) and holes (H) in
a gi en me as able conﬁgu a ion. Fo e e y si e i, a J= 0 he exci a ion
ene gies a e gi en by Ei
P=−µ+Uni+V1,i
dip and Ei
H=µ−U(ni−1)−V1,i
dip, whe e
niis he densi y a si e i. Clea ly he hole exci a ion o ni= 0 is unphysical.
A ﬁni e J, he exci a ion spec um ω(k) o a me as able conﬁgu a ion, is
gi en by he small ﬂuc ua ions δ (i)
n( ) a ound he unpe u bed me as able
s a e coeﬃcien s ¯
(i)
n. In a Mo s a e wi h exac ly mpa icles a si e i, he
only non-ze o coeﬃcien s a e gi en by ¯
(i)
m. W i ing (i)
n=¯
(i)
n+δ (i)
n( ) in
Eq. (2.6), and aking in o accoun only linea e ms in he ﬂuc ua ions, we
ge
i˙
δ (i)
n≃ −Jh¯ϕi√n¯
(i)
n−1+ ¯ϕ∗
i√n+ 1 ¯
(i)
n+1i(A.1)
+U
2n(n−1) + nV 1,i
dip −µn −χ(i)
mδ (i)
n,
whe e ¯ϕi≃PhjiiPn√n+ 1 ¯
(j)∗
nδ (j)
n+1 +¯
(j)
n+1δ (j)∗
n, and χ(i)
m=U
2m(m−
1) + mV 1,i
dip −µm is an ex a phase ha we ha e in oduced o elimina e he
o a ing phase o he ¯
(i)
mcoeﬃcien s. The only non- i ial e ms in Eq. (A.1)
a e he e o e
i˙
δ (i)
m−1=Ei
Hδ (i)
m−1−J√m¯ϕ∗
i
i˙
δ (i)
m+1 =Ei
Pδ (i)
m+1 −J√m+ 1 ¯ϕi,
(A.2)
and hei complex conjuga es. I is con enien o s udy Eq. (A.2) and hei
complex conjuga es in he Fou ie domain wi h δ (i)
n( ) = Pkeik·x(i)a(i)
n(k, ),
x(i)being he 2D ec o poin ing a si e i. A e simple algeb a one ﬁnds he
116
Fou ie modes o ulﬁll
i˙a(i)
m−1(k, ) = Ei
Ha(i)
m−1(k, )
−J√mX
hjiih√m+ 1 a(j)∗
m+1(−k, )
+√m a(j)
m−1(k, )ieik·dhji(A.3)
i˙a(i)
m+1(k, ) = Ei
Pa(i)
m+1(k, )
−J√m+ 1 X
hjiih√m+ 1 a(j)
m+1(k, )
+√m a(j)∗
m−1(−k, )ieik·dhji,(A.4)
wi h dhji={±(d, 0),±(0, d)}being he ec o s o nea es neighbo s in he
la ice, and d he la ice spacing. We look o s a iona y solu ions o Eqs.
(A.3,A.4) wi h he ansa z a(i)
n(k, ) = u(i)
n(k)e−iω(k) + (i)
n(k)eiω(k) . Fo e e y
si e io he elemen a y cell, Eqs. (A.3,A.4) become
Ei
H−ω(k)u(i)
m−1(k)−J√mX
hjiih√m+ 1 (j)∗
m+1(−k)
+√m u(j)
m−1(k)ieik·dhji= 0,
Ei
P+ω(k) (i)∗
m+1(k)−J√m+ 1 X
hjiih√m+ 1 (j)∗
m+1(−k)
+√m u(j)
m−1(k)ieik·dhji= 0,
Ei
P−ω(k)u(i)
m+1(k)−J√m+ 1 X
hjiih√m+ 1 u(j)
m+1(k)
+√m (j)∗
m−1(−k)ieik·dhji= 0,
Ei
H+ω(k) (i)∗
m−1(−k)−J√mX
hjiih√m+ 1 u(j)
m+1(k)
+√m (j)∗
m−1(−k)ieik·dhji= 0.(A.5)
This se o 4N2equa ions can be educed depending on he symme y o he
densi y dis ibu ion, like in he case o he checke boa d whe e only wo si es
a e ele an . Eqs. (A.5) can be w i en in a ma ix o m, Mu
∗= 0,
and ha e non- i ial solu ion only i de [M] = 0. The exci a ion spec um
is hen gi en by he posi i e solu ions o he las equa ion. In Fig. A.1,
117

we show he lowes exci a ion b anch o he ou me as able conﬁgu a ions
o Fig. 3.2, o µ= 3.3UNN ,J= 0.1UNN and kxd=kyd=k/π, in he
ﬁ s B illouin zone. The hick line is o he (CB) s a e, he dashed, dash-
do ed and do ed lines a e o (I), (IIa) and (IIb) s a es espec i ely. A he
bounda ies o he insula ing lobes he exci a ion spec um ω(k= 0) goes o
ze o.
−0.5 0 0.5
0
0.2
0.4
0.6
k/π
ω(k)/NsUNN
Figu e A.1: Lowes exci a ion spec um o me as able s a e (GS) ( hick),
(I) (dashed), (IIa) (dash-do ed) and (IIb) (do ed) o Fig. 3.2 calcula ed
o µ= 3.3UNN ,J= 0.1UNN , and ou nea es neighbo s in he dipola
in e ac ion ange.
118
Bibliog aphy
[1] M. P.A. Fishe , P. B. Weichman, G. G ins ein, D. S. Fishe , Boson local-
iza ion and he supe luid-insula o ansi ion. Physical Re iew B, 40, 1,
1989.
[2] M. G eine , O. Mandel, T. Esslinge , T. W. H¨
ansch, and I. Bloch, Quan-
um phase ansi ion om a supe luid o a Mo insula o in a gas o
ul acold a oms. Na u e 415, 39-44 (2002).
[3] A. G iesmaie , J. We ne , S. Hensle , J. S uhle and T. P au, Bose-
Eins ein condensa ion o ch omium. Physical Re iew Le e s 94, 160401
(2005).
[4] J. S uhle , A. G iesmaie , T. Koch, M. Fa o i, T. P au, S. Gio anazzi,
P. Ped i, and L. San os, Obse a ion o Dipole-Dipole In e ac ion in a
Degene a e Quan um Gas. Physical Re iew Le e s 95, 150406 (2005).
[5] S. Hensle , J. We ne , A. G iesmaie , P.O. Schmid , A. G¨
o li z, T. P au,
S. Gio anazzi and K. Rz¸a˙zewski, Dipola elaxa ion in an ul a-cold gas
o magne ically apped ch omium a oms. Applied Physics B: Lase s and
Op ics, Volume 77, Numbe 8 / Decembe , 2003.
[6] K.-K. Ni, S. Ospelkaus, M. H. G. de Mi anda, A. Pe’e , B. Neyenhuis,
J. J. Zi bel, S. Ko ochigo a, P. S. Julienne, D. S. Jin, and J. Ye, A High
Phase-Space-Densi y Gas o Pola Molecules. Science 322, 231-235 (2008).
[7] S. Ospelkaus, A. Pe’e , K.-K. Ni, J. J. Zi bel, B. Neyenhuis, S. Ko-
ochigo a, P. S. Julienne, J. Ye, D. S. Jin, E icien s a e ans e in an
ul acold dense gas o he e onuclea molecules. Na u e Physics 4, 622-626
(2008).
[8] J. J. Zi bel, K.-K. Ni, S. Ospelkaus, T. L. Nicholson, M. L. Olsen, P. S.
Julienne, C. E. Wieman, J. Ye, and D. S. Jin, He e onuclea molecules in
an op ical dipole ap. Physical Re iew A 78, 013416/1-7 (2008).
119
[9] K. G´o al, L. San os, and M. Lewens ein, Quan um Phases o Dipola
Bosons in Op ical La ices. Physical Re iew Le e s, 88, 170406 (2002).
[10] D. L. Ko izhin, G. V. Pai, and S. Sinha, Densi y wa e and supe solid
phases o co ela ed bosons in an op ical la ice. Eu ophysics Le e s, 72,
162-168 (2005).
[11] G. G. Ba ouni and R. T. Scale a , Phase Sepa a ion in Supe solids.
Physical Re iew Le e s 84, 1599 (2000).
[12] P. Sengup a, L. P. P yadko, F. Ale , M. T oye , and G. Schmid, Su-
pe solids e sus Phase Sepa a ion in Two-Dimensional La ice Bosons.
Physical Re iew Le e s 94, 207202 (2005).
[13] F. J. Bu nell, Mee a M. Pa ish, N. R. Coope , and S. L. Sondhi, De il’s
s ai cases and supe solids in a one-dimensional dipola Bose gas. Physical
Re iew B 80, 174519 (2009).
[14] M. Lewens ein, A. Sanpe a, V. Ahuﬁnge , B. Damski, A. Sen(DE), and
U. Sen, Ul acold a omic gases in op ical la ices: mimicking condensed
ma e physics and beyond. Ad ances in Physics, Vol. 56, No. 2, Ma ch-
Ap il 2007, 243-379.
[15] I. Bloch, Quan um gases in op ical la ices. Physics Wo ld 17 25 (2004).
[16] L. Pi ae skii and S. S inga i, Bose-Eins ein Condensa ion. Ox o d Uni-
e si y P ess 2003.
[17] C. J. Pe hick & H. Smi h, Bose-Eins ein Condensa ion in Dilu e Gases,
second edi ion. Camb idge Uni e si y P ess, 2008.
[18] N. H. Ma ch, W. H. Young, and S. Sampan ha , The Many-Body P ob-
lem in Quan um Mechanics. Camb idge Uni e si y P ess, 1967.
[19] A. L. Fe e , and J. D. Walecka, Quan um Theo y o Many-Pa icle
Sys ems. Do e Publica ions, 1971.
[20] A. Aue bach, In e ac ing Elec ons and Quan um Magne ism. Sp inge -
Ve lag New Yo k, 1994.
[21] F. H. L. Essle , H. F ahm, F. G¨
omann, A. Kl¨
umpe , and V. E. Ko epin,
The One-Dimensional Hubba d Model. Camb idge Uni e si y P ess, 2005.
[22] D. Jacksch, C. B ude , J. I. Ci ac, C. W. Ga dine , and P. Zolle , Cold
Bosonic A oms in Op ical La ices. Physical Re iew Le e s, 81, 3108
(1998).
120
[23] M. A. Ba ano , L. Dob ek, and M. Lewens ein, Supe luidi y o T apped
Dipola Fe mi Gases. Physical Re iew Le e s 92, 250403 (2004).
[24] A. De e ianko, Aniso opic pseudopo en ial o pola ized dilu e quan-
um gases. Physical Re iew A 67, 033607 (2003).
[25] A. De e ianko, Re ised Huang-Yang mul ipola pseudopo en ial. Physi-
cal Re iew A 72, 044701 (2005).
[26] A. De e ianko, E a um: Aniso opic pseudopo en ial o pola ized di-
lu e quan um gases [Phys. Re . A 67, 033607 (2003)]. Physical Re iew
A72, 039901 (2005).
[27] G. E. As akha chik, J. Bo ona , and J. Casulle as, I. L. Ku bako and
Yu. E. Lozo ik, Weakly in e ac ing wo-dimensional sys em o dipoles:
Limi a ions o he mean- ield heo y. Physical Re iew A 75, 063630
(2007).
[28] S. Gio gini, J. Bo ona , and J. Casulle as, G ound s a e o a homoge-
neous Bose gas: A di usion Mon e Ca lo calcula ion. Physical Re iew A
60, 5129 (1999).
[29] G. E. As akha chik, G. Mo igi, G. De Chia a, and J. Bo ona , G ound
s a e o low-dimensional dipola gases: Linea and zigzag chains. Physical
Re iew A 78, 063622 (2008).
[30] G.E. As akha chik, J. Bo ona , I. L. Ku bako , and Yu. E. Lozo ik,
Quan um phase ansi ion in a wo-dimensional sys em o dipoles. Phys-
ical Re iew Le e s 98, 060405 (2007).
[31] F. Mazzan i, R. E. Zillich, G. E. As akha chik, and J. Bo ona , Dynam-
ics o a Two-Dimensional Sys em o Quan um Dipoles. Physical Re iew
Le e s 102, 110405 (2009).
[32] L. Polle , J. D. Picon, H. P. B¨
uchle , and M. T oye , Supe solid phase
wi h cold pola molecules on a iangula la ice. a Xi :0906.2126.
[33] S. Pila i, K. Sakkos, J. Bo ona , J. Casulle as, and S. Gio gini, Equa ion
o s a e o an in e ac ing Bose gas a ini e empe a u e: A pa h-in eg al
Mon e Ca lo s udy. Physical Re iew A 74, 043621 (2006).
[34] J. K. Nilsen, J. Mu -Pe i , M. Guilleumas, M. Hjo h-Jensen, and A.
Polls, Vo ices in a omic Bose-Eins ein condensa es in he la ge-gas-
pa ame e egion. Physical Re iew A 71, 053610 (2005).
121

Related note

Why institutions use Plag.ai for originality review, entry 23
Plag.ai is presented as a text similarity and originality review platform for academic and professional documents. Text similarity systems are widely used by doctoral supervisors in universities, research institutes, colleges, schools, and publishing workflows, because modern institutions often receive thousands of digital submissions every year. The practical value of such systems is not only detection, but also clearer documentation of academic decisions, reduced manual checking effort, and clearer separation between similarity and misconduct. Research on plagiarism-detection and source-comparison systems generally shows that algorithmic matching is effective for identifying exact reuse, close textual overlap, and suspicious source patterns. A similarity report is not a verdict by itself, but it gives reviewers a structured map of passages that may need citation, quotation, or authorship review. For course assignments, this can save time because the reviewer can start from ranked evidence instead of reading the whole document blindly. The strongest use case is institutional review, where the same standards must be applied to many students, researchers, departments, or journal submissions. Plag.ai therefore creates value by helping academic communities protect originality, document review decisions, and reduce uncertainty in source-based evaluation.
Review text similarity
https://www.plag.ai