5. BEC excitation by modulation of scatter<strong>in</strong>g lengthfor the phases φ σ and ϕ σ can be solved explicitly <strong>in</strong> terms of the widths u σ andthe center of mass coord<strong>in</strong>ates σ 0 . In this way we obta<strong>in</strong> two sets of ord<strong>in</strong>arysecond-order differential equations that govern the condensate dynamics. A centerof-massmotion is decoupled from the shape oscillations and is determ<strong>in</strong>ed by simpleharmonic oscillator equations, which are <strong>in</strong>dependent of <strong>in</strong>teratomic <strong>in</strong>teractions:¨σ 0 (t) + λ 2 σσ 0 (t) = 0, σ ∈ {x, y, z} . (5.2)On the other hand, the widths of the condensate exhibit non-trivial dynamics givenby a set of coupled nonl<strong>in</strong>ear differential equations:ü σ (t) + λ 2 σu σ (t) − 1u σ (t) − P= 0, σ ∈ {x, y, z} . (5.3)3 u σ (t)u x (t)u y (t)u z (t)In the previous equations and from now on, we use dimensionless notation: wechoose a convenient frequency scale ω (for example, the external trap frequency <strong>in</strong>one of the spatial directions) and express all lengths <strong>in</strong> the units of the characteristicharmonic oscillator length l = √ /Mω, time <strong>in</strong> units of ω −1 and external frequencies<strong>in</strong> units of ω: λ σ = ω σ /ω, σ ∈ {x, y, z}. The dimensionless <strong>in</strong>teraction parameter Pis given by P = g/((2π) 3/2 ωl 3 ) = √ 2/πNa/l.In this approach, the <strong>in</strong>itial partial differential equation (4.15) is approximatedwith the two sets of ord<strong>in</strong>ary differential equations, given by Eqs. (5.2) and (5.3),which allow analytical considerations. The first properties of the condensate thatcan be calculated are the equilibrium widths u x0 , u y0 and u z0 . They are found bysolv<strong>in</strong>g an algebraic system of equations:λ 2 σu σ0 − 1 P− = 0, σ ∈ {x, y, z} . (5.4)u σ0 u σ0 u x0 u y0 u z0The equilibrium widths represent stationary solutions of Eqs. (5.3).Now we turn to the calculation of the frequencies of low-ly<strong>in</strong>g modes. To beg<strong>in</strong>with, from Eqs.(5.2) we read off frequencies that correspond to the center-of-massmotion. These are dipole modes and their frequencies are equal to the external trapfrequencies (for the case of a harmonic trap). Most often, this type of excitationsis created by shift<strong>in</strong>g the trap <strong>in</strong> space. Actually, the well established experimentalprocedure for the precise determ<strong>in</strong>ation of the trap parameters is based on themeasurement of the dipole mode frequencies [21].109
5. BEC excitation by modulation of scatter<strong>in</strong>g lengthNext, by l<strong>in</strong>eariz<strong>in</strong>g the Eqs. (5.3) around the equilibrium widths (5.4), we canobta<strong>in</strong> <strong>in</strong>formation on collective modes related to the BEC shape oscillations. Wewill consider experimentally relevant case of an axially symmetric trap, such thatλ x = λ y = 1, u x0 = u y0 ≡ u ρ0 . Due to the axial symmetry of the consideredsystem, the projection of the angular momentum along the z-axis, L z , is a goodquantum number and we use it for the classification of the modes. By solv<strong>in</strong>gthe correspond<strong>in</strong>g eigenproblem, we f<strong>in</strong>d three different modes depicted <strong>in</strong> Fig. 5.2.Accord<strong>in</strong>g to their symmetry properties, we designate the modes as the |L z | = 2quadrupole mode, the L z = 0 quadrupole mode, and the breath<strong>in</strong>g mode (whichalso corresponds to L z = 0). Eigenvalues of the l<strong>in</strong>earized system of equations yieldthe frequencies of the collective modes. Frequencies of the L z = 0 quadrupole modeω Q0 , and the breath<strong>in</strong>g mode ω B0 are given by:ω B0,Q0 = √ 2[ (1 + λ 2 z − P4u 2 ρ0 u3 z0)√ (± 1 − λ 2 z +P4u 2 ρ0 u3 z0while the frequency of the |L z | = 2 quadrupole mode is given by:ω |Lz|=2Q0=) 2 ( )]1/2 2 P+ 8,4u 3 ρ0 u2 z0(5.5)√P4 − 2 . (5.6)u 4 ρ0u z0As shown <strong>in</strong> Fig. 5.2, the |L z | = 2 quadrupole mode is characterized by out-ofphaseoscillations <strong>in</strong> x and y directions, the |L z | = 0 quadrupole mode exhibitsout-of-phase radial and axial oscillations, while <strong>in</strong>-phase radial and axial oscillationscorrespond to the breath<strong>in</strong>g mode.Figure 5.2: A schematic illustration of the condensate eigenmodes: |L z | = 2quadrupole mode (left), |L z | = 0 quadrupole mode (middle) and |L z | = 0 breath<strong>in</strong>gmode (right).110
- 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 and 108: 4. Mean-field description of an int
- 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: 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ć