4. Mean-field description of an <strong>in</strong>teract<strong>in</strong>g BECthe vic<strong>in</strong>ity of the BEC phase transition. We focus on harmonically trapped bosonsfor two reasons: the trapp<strong>in</strong>g is always present <strong>in</strong> the experiments, and also thecalculation of the condensation temperature of the homogenous system turned out tobe a notorious problem, debated <strong>in</strong> many ways for several decades [79]. We scrut<strong>in</strong>yand compare several widely used implementations of the HF approximation with theemphasis on their advantages and drawbacks. The <strong>in</strong>terest <strong>in</strong> the topic <strong>in</strong>creases asnew experiments have reported the observation of beyond mean-field effects [80, 81],which should be <strong>in</strong>corporated <strong>in</strong>to the exist<strong>in</strong>g models.The partition function of the system <strong>in</strong> the grand canonical ensemble is given byZ(β) = Tr e −β(Ĥ−µ ˆN) , (4.2)and can be rewritten as a bosonic functional <strong>in</strong>tegral <strong>in</strong> the imag<strong>in</strong>ary time [32]:∮Z(β) =∮DΨDΨ ∗ e −A E[Ψ(⃗r,τ),Ψ ∗ (⃗r,τ)]/ , (4.3)where Ψ(⃗r, τ) and Ψ ∗ (⃗r, τ) are periodic functions, with the period β:Ψ(⃗r, τ) = Ψ(⃗r, τ + β) , Ψ ∗ (⃗r, τ) = Ψ ∗ (⃗r, τ + β) . (4.4)For the Hamiltonian (4.1) the Euclidean action A E is given byA E [Ψ(⃗r, τ), Ψ ∗ (⃗r, τ)] =∫ β0+ g 2∫dτ∫ β0(d⃗r Ψ ∗ (⃗r, τ) ∂)∂τ − 22M △ + V (⃗r) − µ Ψ(⃗r, τ)∫d⃗r Ψ ∗ (⃗r, τ)Ψ(⃗r, τ)Ψ ∗ (⃗r, τ)Ψ(⃗r, τ) . (4.5)dτA presence of the <strong>in</strong>teract<strong>in</strong>g Ψ 4 term <strong>in</strong> the action makes the calculation of the partitionfunction analytically <strong>in</strong>tractable, and to proceed further we apply the standardmean-field approach. In order to study Bose-E<strong>in</strong>ste<strong>in</strong> condensation, accord<strong>in</strong>g tothe Bogoliubov prescription (1.27), we first decompose the field Ψ <strong>in</strong>to the orderparameter ψ(⃗r, τ), which corresponds to the macroscopic condensate wave-function,and fluctuations δψ(⃗r, τ):Ψ(⃗r, τ) = ψ(⃗r, τ) + δψ(⃗r, τ) . (4.6)In the functional formalism this represents a change of variables, and the action now85
4. Mean-field description of an <strong>in</strong>teract<strong>in</strong>g BECconta<strong>in</strong>s terms up to the 4 th order <strong>in</strong> δψ that we should <strong>in</strong>tegrate over. Withoutfurther approximations, however, this problem is equivalent to the orig<strong>in</strong>al functional<strong>in</strong>tegral. Numerous approximation techniques were developed to treat terms of the3 rd and 4 th order <strong>in</strong> different approximative ways [19, 82]. In the Hartree-Fock-Bogoliubov approach, we approximate the 4 th order term asδψ ∗ δψ δψ ∗ δψ ≈4〈δψ ∗ δψ〉δψ ∗ δψ + 〈δψ ∗ δψ ∗ 〉δψδψ + 〈δψδψ〉δψ ∗ δψ ∗−2〈δψ ∗ δψ〉〈δψ ∗ δψ〉 − 〈δψδψ〉〈δψ ∗ δψ ∗ 〉. (4.7)In accordance with the previous decomposition, we <strong>in</strong>troduce auxiliary functionsh(⃗r, τ;⃗r, τ) = 〈δψ ∗ (⃗r, τ)δψ(⃗r, τ)〉 ,f(⃗r, τ;⃗r ′ , τ ′ ) = 〈δψ ∗ (⃗r, τ)δψ(⃗r ′ , τ ′ )〉 ,b(⃗r, τ;⃗r ′ , τ ′ ) = 〈δψ(⃗r, τ)δψ(⃗r ′ , τ ′ )〉 ,which are denoted as Hartree, Fock and Bogoliubov term, respectively. At themoment, the <strong>in</strong>troduced average values are purely formal, but can be later def<strong>in</strong>edso as to make the complete procedure self-consistent. In the case of the contact<strong>in</strong>teraction (1.21), the Hartree and Fock terms yield equal contributions, hence afactor of 4 <strong>in</strong> front of the correspond<strong>in</strong>g term <strong>in</strong> Eq. (4.7).After apply<strong>in</strong>g the mean-field approximation (4.7), the action A E becomes quadratic<strong>in</strong> δψ, and now the functional <strong>in</strong>tegrations of the Gaussian <strong>in</strong>tegrals can be explicitlyperformed. The f<strong>in</strong>al result for the partition function can be written <strong>in</strong> theformZ(β) = e −β Γ eff[ψ,ψ ∗ ,h,f,b] , (4.8)where Γ eff is the effective action, def<strong>in</strong>ed as a functional of five arguments: ψ(⃗r, τ),ψ ∗ (⃗r, τ), h(⃗r, τ;⃗r, τ), f(⃗r, τ;⃗r ′ , τ ′ ) and b(⃗r, τ;⃗r ′ , τ ′ ). They are determ<strong>in</strong>ed by extremiz<strong>in</strong>gthe effective action Γ eff [ψ, ψ ∗ , h, f, b] with respect to each of them:δΓ effδψ = 0, δΓ effδψ ∗= 0, δΓ effδh = 0, δΓ effδf = 0, δΓ effδb= 0.In the quest for the simplest mean-field description of an <strong>in</strong>homogeneous BEC,we will make another simplification by neglect<strong>in</strong>g Bogoliubov terms, i.e. anomalouscorrelations b(⃗r, τ;⃗r ′ , τ ′ ), as discussed <strong>in</strong> Ref. [83]. This assumption is justified <strong>in</strong>86
- Page 1 and 2:
UNIVERSITY OF BELGRADEFACULTY OF PH
- Page 3 and 4:
Thesis advisor, Committee member:Dr
- Page 6 and 7:
lence for Computer Modeling of Comp
- Page 8 and 9:
dobijanje kondenzata odabrani su at
- Page 10 and 11:
Uticaj slabih interakcija na fenome
- Page 12 and 13:
Abstract of the doctoral dissertati
- Page 14 and 15:
highly accurate information on ener
- Page 16 and 17:
Keywords: cold quantum gases, Bose-
- Page 18 and 19:
CONTENTS3.4.2 Time-of-flight graphs
- Page 20 and 21:
NomenclatureRoman Symbolsagk BLMNn(
- Page 22 and 23:
Chapter 1Introduction1.1 ForewordTh
- Page 24 and 25:
ently explored to illustrate the ve
- Page 26 and 27:
Summations in the last expression c
- Page 28 and 29:
Figure 1.1: The hallmark of the Bos
- Page 30 and 31:
we discuss in some detail the exper
- Page 32 and 33:
In the first papers [3, 4], the TOF
- Page 34 and 35:
where a BG is the off-resonant scat
- Page 36 and 37:
system given by( ) ǫ Bog ⃗k =
- Page 38 and 39:
Having the efficient numerical meth
- Page 40 and 41:
Chapter 2Properties of quantum syst
- Page 42 and 43:
2. Diagonalization of Transition Am
- Page 44 and 45:
2. Diagonalization of Transition Am
- Page 46 and 47:
2. Diagonalization of Transition Am
- Page 48 and 49:
2. Diagonalization of Transition Am
- Page 50 and 51:
2. Diagonalization of Transition Am
- Page 52 and 53:
2. Diagonalization of Transition Am
- Page 54 and 55:
2. Diagonalization of Transition Am
- Page 56 and 57: 2. Diagonalization of Transition Am
- Page 58 and 59: 2. Diagonalization of Transition Am
- Page 60 and 61: 2. Diagonalization of Transition Am
- Page 62 and 63: 2. Diagonalization of Transition Am
- Page 64 and 65: 2. Diagonalization of Transition Am
- Page 66 and 67: ||2. Diagonalization of Transition
- Page 68 and 69: 2. Diagonalization of Transition Am
- Page 70 and 71: 2. Diagonalization of Transition Am
- Page 72 and 73: 2. Diagonalization of Transition Am
- Page 74: 2. Diagonalization of Transition Am
- Page 77 and 78: 2. Diagonalization of Transition Am
- Page 79 and 80: magnetic field ⃗ B = 2M ⃗ Ω.3.
- Page 81 and 82: 3. Rotating ideal BECof the rotatio
- Page 83 and 84: 3. Rotating ideal BEC3.1 Numerical
- Page 85 and 86: 3. Rotating ideal BEC350300SC appro
- Page 87 and 88: 3. Rotating ideal BEC3.2 Finite num
- Page 89 and 90: 3. Rotating ideal BECN - N 03.0·10
- Page 91 and 92: 3. Rotating ideal BEC3.3 Global pro
- Page 93 and 94: 3. Rotating ideal BECever, it does
- Page 95 and 96: 3. Rotating ideal BECT c [nK]120110
- Page 97 and 98: 3. Rotating ideal BEC3.533210 3 κ
- Page 99 and 100: 3. Rotating ideal BECn(x, y)10·10
- Page 101 and 102: 3. Rotating ideal BECpancy numbers
- Page 103 and 104: 3. Rotating ideal BEC6·10 4 0 0.02
- Page 105: Chapter 4Mean-field description of
- Page 109 and 110: 4. Mean-field description of an int
- Page 111 and 112: 4. Mean-field description of an int
- Page 113 and 114: 4. Mean-field description of an int
- Page 115 and 116: 4. Mean-field description of an int
- Page 117 and 118: 4. Mean-field description of an int
- Page 119 and 120: 4. Mean-field description of an int
- Page 121 and 122: 4. Mean-field description of an int
- Page 123 and 124: 4. Mean-field description of an int
- Page 125 and 126: Chapter 5Nonlinear BEC dynamics by
- Page 127 and 128: 5. BEC excitation by modulation of
- Page 129 and 130: 5. BEC excitation by modulation of
- Page 131 and 132: 5. BEC excitation by modulation of
- Page 133 and 134: 5. BEC excitation by modulation of
- Page 135 and 136: 5. BEC excitation by modulation of
- Page 137 and 138: 5. BEC excitation by modulation of
- Page 139 and 140: 5. BEC excitation by modulation of
- Page 141 and 142: 5. BEC excitation by modulation of
- Page 143 and 144: 5. BEC excitation by modulation of
- Page 145 and 146: 5. BEC excitation by modulation of
- Page 147 and 148: 5. BEC excitation by modulation of
- Page 149 and 150: where n = 1, 2, 3, . . . is an inte
- Page 151 and 152: 5. BEC excitation by modulation of
- Page 153 and 154: 5. BEC excitation by modulation of
- Page 155 and 156: Chapter 6SummarySince the first exp
- Page 157 and 158:
short-range interactions.The excita
- Page 159 and 160:
of the ground state. Imaginary-time
- Page 161 and 162:
(A.13). With another boundary condi
- Page 163 and 164:
Appendix B Time-dependent variation
- Page 165 and 166:
List of papers by Ivana VidanovićT
- Page 167 and 168:
References[1] S. N. Bose, Plancks g
- Page 169 and 170:
REFERENCES[21] W. Ketterle, D. S. D
- Page 171 and 172:
REFERENCES[45] A. Bogojević, A. Ba
- Page 173 and 174:
REFERENCES[70] M. R. Matthews, B. P
- Page 175 and 176:
REFERENCES[94] M.-O. Mewes, M. R. A
- Page 177 and 178:
REFERENCES[116] K. Staliunas, S. Lo
- Page 179 and 180:
CURRICULUM VITAE - Ivana Vidanović