13.07.2015 Views

PhD thesis in English

PhD thesis in English

PhD thesis in English

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

separate treatment.5. BEC excitation by modulation of scatter<strong>in</strong>g lengthIn Ref. [116] it was predicted that a harmonic modulation of the scatter<strong>in</strong>g lengthleads to the creation of Faraday patterns <strong>in</strong> BEC, i.e. density waves. Up to now,Faraday patterns have been experimentally <strong>in</strong>duced by harmonic modulation of thetransverse conf<strong>in</strong>ement strength [111], which is studied analytically and numerically<strong>in</strong> Ref. [129]. Here we focus only on the nonl<strong>in</strong>ear properties of low-ly<strong>in</strong>g collectivemodes and do not consider possible excitations of Faraday patterns.In order to obta<strong>in</strong> analytical <strong>in</strong>sight <strong>in</strong>to the condensate dynamics <strong>in</strong>duced bythe harmonic modulation of the s-wave scatter<strong>in</strong>g length described by Eq. (5.7), weuse the Gaussian variational approximation approximation. We consider an axiallysymmetric BEC, (λ x = λ y = 1, u x = u y ≡ u ρ ), excited by modulation of the<strong>in</strong>teraction strength, which preserves the axial symmetry of the condensate dur<strong>in</strong>gits time evolution. For this reason, we use a simplified axially-symmetric form ofEqs. (5.3):ü ρ (t) + u ρ (t) − 1u ρ (t) 3 −P(t)u ρ (t) 3 u z (t)= 0 , (5.11)ü z (t) + λ 2 zu z (t) − 1u z (t) − P(t)= 0 , (5.12)3 u ρ (t) 2 u z (t) 2which we will refer to as a Gaussian approximation.To estimate the accuracy of the Gaussian approximation for describ<strong>in</strong>g the dynamics<strong>in</strong>duced by the harmonic modulation of the <strong>in</strong>teraction strength, we compareits solution with an exact numerical solution of the GP equation. In Fig. 5.5, weplot the result<strong>in</strong>g time-dependent axial and radial condensate widths ρ rms (t) andz rms (t), calculated as root mean square valuesρ rms (t) =z rms (t) =√√2π2π∫ ∞−∞∫ ∞−∞dzdz∫ ∞0∫ ∞0ρ dρ |ψ(ρ, z, t)| 2 ρ 2 , (5.13)ρ dρ |ψ(ρ, z, t)| 2 z 2 , (5.14)of the solution of the GP equation, as well as numerical solutions of Eqs. (5.11)and (5.12). We assume that <strong>in</strong>itially the condensate is <strong>in</strong> the ground state. In thevariational description, this translates <strong>in</strong>to <strong>in</strong>itial conditions u ρ (0) = u ρ0 , ˙u ρ (0) = 0,u z (0) = u z0 , ˙u z (0) = 0, where u ρ0 and u z0 are the time-<strong>in</strong>dependent solutions115

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!