4. Mean-field description of an <strong>in</strong>teract<strong>in</strong>g BECthe limit T → 0, however it is very often used for all temperatures. The topic isexplored <strong>in</strong> detail <strong>in</strong> Ref. [84] where the consequences of the approximation are thoroughlydiscussed. Additionally, we neglect the possible depletion of the condensateat the zero temperature, i.e. the depletion that arises due to <strong>in</strong>teractions, which isa reasonable approximation <strong>in</strong> the case of a weakly <strong>in</strong>teract<strong>in</strong>g gas.F<strong>in</strong>ally, after implement<strong>in</strong>g all the described steps, we arrive at the f<strong>in</strong>ite-temperatureHF description of a bosonic gas, which is given by the follow<strong>in</strong>g system ofequations:[ ∂]∂τ − 22M △ + V (⃗r) + g|ψ(⃗r, τ)|2 + 2 gh(⃗r, τ;⃗r, τ) ψ(⃗r, τ) = µψ(⃗r, τ) , (4.9)h(⃗r, τ;⃗r, τ) = ∑ ψ ⃗k (⃗r)ψ ∗ 1⃗ k(⃗r)e β(E ⃗ −µ) k − 1 , (4.10)⃗ k][− 22M △ + V (⃗r) + 2 g|ψ(⃗r, τ)|2 + 2 gh(⃗r, τ;⃗r, τ) ψ ⃗k (⃗r) = E ⃗k ψ ⃗k (⃗r) , (4.11)where ψ ⃗k (⃗r) are effective s<strong>in</strong>gle-particle wave-functions, and E ⃗k are the correspond<strong>in</strong>geigenvalues. More details on the derivation can be found <strong>in</strong> Ref. [78]. Althoughformally the HF equations depend on the imag<strong>in</strong>ary time τ, physically is only relevantthe equilibrium case, when the macroscopic wave-function of the condensateψ(⃗r, τ) does not depend on τ anymore (∂ψ(⃗r, τ)/∂τ = 0), but only on the position⃗r. The above set of equations has to be solved self-consistently, tak<strong>in</strong>g <strong>in</strong>to accountthat the total number of particles is fixed to N, and leads to the solution that consistsof the effective s<strong>in</strong>gle-particle eigenfunctions ψ ⃗k and eigenvalues E ⃗k , the Hartreefunction h, and the condensate wave-function ψ. From Eq. (4.10) we immediatelysee the physical <strong>in</strong>terpretation of the Hartree function h, which represents the densityof the thermal cloud, n th (⃗r). Effectively, with<strong>in</strong> the mean-field description, thegas of bosons is split <strong>in</strong>to the condensate and thermal component. We note the closeanalogy with the non<strong>in</strong>teract<strong>in</strong>g gas description presented <strong>in</strong> Chapter 1, with theimportant exception that the two components now mutually <strong>in</strong>teract. By vary<strong>in</strong>gthe total number of particles N and the temperature T, a complete N −T phase diagramcan be explored. Before consider<strong>in</strong>g a BEC phase transition <strong>in</strong> the mean-fieldapproximation, we first present the zero-temperature limit of the HF approximation.87
4. Mean-field description of an <strong>in</strong>teract<strong>in</strong>g BEC4.1 Gross-Pitaevskii equationIn the zero-temperature limit, we can neglect the thermal cloud, and setlim n th(⃗r) = 0 , (4.12)β→∞which stems directly from Eq. (4.10). In this case, Eq. (4.11) becomes irelevant,while from Eq. (4.9) we f<strong>in</strong>d the time-<strong>in</strong>dependent Gross-Pitaevskii (GP) equation[13, 14, 34] for the order parameter:[− 22M △ + V (⃗r) + g|ψ(⃗r)|2 ]ψ(⃗r) = µψ(⃗r) . (4.13)Effectively, <strong>in</strong> the mean-field approximation at zero temperature, we assume thatall atoms occupy the same state ψ(⃗r), which we denote as the condensate wavefunction.Note that already by neglect<strong>in</strong>g anomalous averages we have disregardeda depletion of the condensate at zero temperature that arises due to <strong>in</strong>teractions.However, it turns out that this is a reasonable approximation <strong>in</strong> the wide range ofexperimental parameters for a weakly <strong>in</strong>teract<strong>in</strong>g gas. On the left-hand side of theGP equation (4.13) we have a k<strong>in</strong>etic energy term, an external trap potential V (⃗r),and a nonl<strong>in</strong>ear term orig<strong>in</strong>at<strong>in</strong>g from the mean-field <strong>in</strong>teraction between the atoms.The GP equation belongs to the class of nonl<strong>in</strong>ear Schröd<strong>in</strong>ger equations, which areextensively studied also <strong>in</strong> the field of nonl<strong>in</strong>ear optics [85, 86].Now, let us analyze solutions of the GP equation. To beg<strong>in</strong> with, we note thatthe non<strong>in</strong>teract<strong>in</strong>g limit is straightforwardly reproduced: for g → 0, the condensatewavefunction is the ground state of the external potential V (⃗r), and the value of thechemical potential is equal to the ground-state energy. In the limit of strong repulsive<strong>in</strong>teractions (for a large number of atoms, for example), the term correspond<strong>in</strong>gto the k<strong>in</strong>etic energy can be safely neglected. In that case we f<strong>in</strong>d an algebraicstationary solution|ψ(⃗r)| 2 = 1 (µ − V (⃗r)) θ(µ − V (⃗r)), (4.14)gwhich is the well-known Thomas-Fermi (TF) solution [19]. The value of the chemicalpotential µ is determ<strong>in</strong>ed as usual, by fix<strong>in</strong>g the total number of atoms <strong>in</strong> the systemto N. In particular, for the harmonic trap, the solution for the non<strong>in</strong>teract<strong>in</strong>g case isa Gaussian, while the TF solution yields a parabolic profile. It is easy to understand88
- 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 and 106: Chapter 4Mean-field description of
- Page 107: 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ć