Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...
Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...
Bogoliubov Excitations of Inhomogeneous Bose-Einstein ...
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 .