5. BEC excitation by modulation of scatter<strong>in</strong>g lengthwhere P = √ 2/πNa av /l denotes the average <strong>in</strong>teraction strength, Q = √ 2/πNδ a /lis a modulation amplitude, and Ω represents the modulation or driv<strong>in</strong>g frequency.The dimensionless experimental parameters from Ref. [15] have the valuesP = 15, Q = 10, λ z = 0.021 , (5.9)correspond<strong>in</strong>g to the highly elongated trap with the modulated but always positive(repell<strong>in</strong>g) <strong>in</strong>teraction. In the experiment, the oscillations of the condensate sizewere observed by <strong>in</strong>-situ phase-contrast imag<strong>in</strong>g, presented <strong>in</strong> Fig. 5.4.Figure 5.4: Oscillations of the BEC cloud, presented by a series of density profilestaken at equidistant time step of 15 ms. Results are taken from Ref. [15].The <strong>in</strong>terpretation of the experimental data was based on the analytical resultsfor the frequencies of the low-ly<strong>in</strong>g collective modes obta<strong>in</strong>ed from the l<strong>in</strong>earizedform of the Gaussian approximation. For experimental data, Eq. (5.5) yields thefollow<strong>in</strong>g values for the frequencies of the quadrupole and the breath<strong>in</strong>g mode:ω Q0 = 0.035375 , ω B0 = 2.00002 . (5.10)The external trap was stationary, thus prevent<strong>in</strong>g excitations of the dipole (Kohn)mode, correspond<strong>in</strong>g to the center-of-mass motion. For the specific set of experimentalparameters basically only the quadrupole oscillation mode was excited <strong>in</strong> thisway and resonances located at the quadrupole frequency and its second harmonicwere observed.There are several advantages of such an experimental scheme: for <strong>in</strong>stance, <strong>in</strong>future experiments with multi-species BEC, a s<strong>in</strong>gle component could be <strong>in</strong>dividuallyexcited <strong>in</strong> this way, while the excitation of other components would occur only<strong>in</strong>directly.113
5. BEC excitation by modulation of scatter<strong>in</strong>g lengthHowever, due to the nonl<strong>in</strong>ear form of the underly<strong>in</strong>g GP equation, we expectnonl<strong>in</strong>earity-<strong>in</strong>duced shifts <strong>in</strong> the frequencies of low-ly<strong>in</strong>g modes compared to thevalues obta<strong>in</strong>ed from Eq. (5.5), calculated us<strong>in</strong>g the l<strong>in</strong>ear stability analysis. Inparticular, <strong>in</strong> the case of a close match<strong>in</strong>g of the driv<strong>in</strong>g frequency Ω and one ofthe BEC eigenmodes, we expect resonances, i.e. large amplitude oscillations of thecondensate size. When this happens, the role of the nonl<strong>in</strong>ear terms becomes crucialand nonl<strong>in</strong>ear phenomena become dom<strong>in</strong>ant. Furthermore, we emphasize thatoscillations with very small amplitudes, which occur <strong>in</strong> the l<strong>in</strong>ear regime, are experimentallyhard to observe. On the other hand, very large amplitude oscillationslead to a fragmentation of the condensate [15, 123]. Thus, the case <strong>in</strong> between is ofthe ma<strong>in</strong> experimental <strong>in</strong>terest and represents our ma<strong>in</strong> objective, as we discuss <strong>in</strong>the next section.5.3 Harmonic modulation of the s-wave scatter<strong>in</strong>g length:theoretical frameworkTo study nonl<strong>in</strong>ear BEC dynamics, we use an approach that is complementary tothe recent theoretical considerations [113, 114, 115, 116] of a BEC with harmonicallymodulated <strong>in</strong>teraction. In Ref. [114] the real-time dynamics of a sphericallysymmetric BEC was numerically studied and analytically expla<strong>in</strong>ed us<strong>in</strong>g the resonantBogoliubov-Mitropolsky method [124]. On the other hand, <strong>in</strong> our approach <strong>in</strong>order to discern <strong>in</strong>duced dynamical features, we look directly at the excitation spectrumobta<strong>in</strong>ed from a Fourier transform of the time-dependent condensate width.From this type of numerical analysis we f<strong>in</strong>d characteristic nonl<strong>in</strong>ear properties:higher harmonic generation, nonl<strong>in</strong>ear mode coupl<strong>in</strong>g, and significant shifts <strong>in</strong> thefrequencies of collective modes with respect to their l<strong>in</strong>ear response counterparts.In addition, we work out an analytic perturbative theory with respect to the modulationamplitude, capable of captur<strong>in</strong>g many of the mentioned nonl<strong>in</strong>ear effectsobta<strong>in</strong>ed numerically.Nonl<strong>in</strong>earity-<strong>in</strong>duced frequency shifts were considered previously <strong>in</strong> Ref. [117] forthe case of bosonic collective modes excited by modulation of the trapp<strong>in</strong>g potential,and also <strong>in</strong> Ref. [125] for the case of a superfluid Fermi gas <strong>in</strong> the BCS-BEC crossover.Our analytical approach is based on the Po<strong>in</strong>caré-L<strong>in</strong>dstedt method [126, 127, 128,124], <strong>in</strong> the same spirit as presented <strong>in</strong> Refs. [117, 125, 126]. However, the harmonicmodulation of the <strong>in</strong>teraction strength <strong>in</strong>troduces additional features that require a114
- 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 and 130: 5. BEC excitation by modulation of
- Page 131 and 132: 5. BEC excitation by modulation of
- Page 133: 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ć