27.11.2012 Views

Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...

Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...

Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2. The <strong>Inhomogeneous</strong> <strong>Bogoliubov</strong> Hamiltonian<br />

+k<br />

0<br />

−k<br />

+θ0<br />

−θ0<br />

−k 0 +k<br />

k<br />

−k 0 +k<br />

(a) |δΨk ′|2<br />

(b) Im(γk ′)<br />

Figure 2.7.: Fourier analysis <strong>of</strong> the stationary scattering state obtained by numerical integration<br />

<strong>of</strong> the GP equation. (a) Density plot |δΨk ′| 2 . The components <strong>of</strong> the scattered<br />

wave are distributed on the elastic circle |k ′ | = k with characteristic nodes at ±θ0.<br />

(b) Scattering state in terms <strong>of</strong> the <strong>Bogoliubov</strong> amplitude Im(γk ′). This figure is more<br />

distinct than (a), because the interference <strong>of</strong> <strong>Bogoliubov</strong> waves with ±k ′ is eliminated.<br />

The sign change across the scattering node θ0 is clearly visible. For the numerical parameters<br />

V (0) = 0.25µ, kξ = 1, kr0 = 0.5, scattering is suppressed around θ0 ≈ π/3.<br />

smaller than the angular k-space resolution 1/kL. Then, the last contribution<br />

with negative A(kξ, θ) is the forward scattering element, which can<br />

be absorbed by shifting the origin <strong>of</strong> the single-particle energy ɛk, and we<br />

recover the Hamiltonian for the potential scattering <strong>of</strong> free matter waves<br />

[22, 23].<br />

Finally, we insert the elastic scattering amplitude (2.58) into equation (2.56)<br />

and find for the two-dimensional scattering problem 3<br />

Im[γ (s) (θ)/γ0] = − 1<br />

4µξ<br />

2.4.4. Numerical verification<br />

+k<br />

0<br />

−k<br />

+θ0<br />

−θ0<br />

kξ<br />

1 + k 2 ξ 2 A(kξ, θ) V 2k sin(θ/2) . (2.61)<br />

The predicted elastic scattering amplitude (2.61) is an ideal opportunity for<br />

a numerical check <strong>of</strong> the formalism so far. The numerical simulation follows<br />

the experimental proposal from above, recall figure 2.5. In the numerics, the<br />

time-dependent Gross-Pitaevskii equation (2.16) is integrated. It relies neither<br />

on the linearization <strong>of</strong> small excitation amplitudes, nor on perturbation<br />

3 Note that in two dimensions the system size L does not enter explicitly in the formula. In general,<br />

44<br />

there is a factor L 2−d<br />

d �<br />

2 − on the right hand side. The Fourier transforms fk = L 2 d −ik·r d re f(r) <strong>of</strong> the<br />

quantities f = γ (0) , γ (s) d − , V scale differently with L: Vk ∝ L 2 , γ (0) ∝ L d<br />

d 1−<br />

2 , and γ (s) ∝ L<br />

k<br />

2 .

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

Saved successfully!

Ooh no, something went wrong!