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

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

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

4. Disorder—Results and Limiting Cases<br />

contributions in (4.27). The positive part comes from the W (2) contribution,<br />

whereas the negative parts enter via (4.26) and stems from the VG0V part.<br />

The first step in evaluating equation (4.27) is to compute the angular integral<br />

over cos 2 (β), where β is taken as polar angle in d-dimensional spherical<br />

coordinates. The angular average <strong>of</strong> cos 2 (β) decreases with the dimension<br />

and is found as 1/d, cf. discussion before (4.12). Thus, in higher dimensions<br />

the negative contribution to the integral (4.27) loses weight with respect to<br />

the positive contribution. The remaining radial integral reads<br />

�<br />

� 2 2 q ξ<br />

�<br />

. (4.28)<br />

ΛN = 2 Sd<br />

(2π) d<br />

dq q d−1 d Cd(qσ)<br />

σ<br />

(2 + q2ξ2 ) 2<br />

2 + q2 1<br />

−<br />

ξ2 d<br />

This result is plotted in figure 4.7(a) for speckle disorder in d = 1, 2, 3<br />

dimensions.<br />

Limits<br />

The limit ( σ<br />

ξ → ∞)k=0 touches the limit (kσ → 0)ξ=0 <strong>of</strong> the previous section<br />

in the lower right corner ξ ≪ σ ≪ k −1 <strong>of</strong> the parameter space represented<br />

in figure 4.1. Thus, the result (4.19a) is recovered.<br />

In the other limit <strong>of</strong> δ-correlated disorder, new results are found: For<br />

ξ ≫ σ, the smoothing factor (2 + q 2 ξ 2 ) −2 in (4.28) is sharply peaked at<br />

q = 0 with a width <strong>of</strong> ξ −1 , so that the correlator Cd(qσ) can be evaluated<br />

as Cd(0). Substituting y = qξ and integrating the first part in the bracket<br />

<strong>of</strong> (4.28) by parts, we find<br />

ΛN(k = 0, σ = 0) = 2 Sd<br />

(2π) d<br />

σdCd(0) ξd = σd Cd(0)<br />

ξ d<br />

� d<br />

4<br />

� � ∞<br />

1<br />

−<br />

d 0<br />

, d = 1<br />

dy<br />

⎧<br />

⎪⎨ −<br />

×<br />

⎪⎩<br />

3<br />

16 √ 2<br />

0, d = 2<br />

, d = 3.<br />

+ 5<br />

48 √ 2π<br />

y d−1<br />

(2 + y 2 ) 2<br />

(4.29)<br />

Remarkably, the correction is negative in one dimension, vanishes in two<br />

dimensions and is positive in three dimensions. In the scaling with V0/µ<br />

fixed, used so far, the correction vanishes for σ → 0, figure 4.7(a). The limit<br />

<strong>of</strong> uncorrelated disorder becomes well-defined, when Pd(0) = V 2<br />

0 σ d Cd(0) is<br />

kept fixed while σ is decreased, figure 4.7(b). The 1D result will be studied in<br />

detail by means <strong>of</strong> a numerical integration <strong>of</strong> the Gross-Pitaevskii equation<br />

in section 4.3. As expected, the 3D result coincides exactly with the findings<br />

<strong>of</strong> Giorgini, Pitaevskii and Stringari [78], which have been confirmed by<br />

Lopatin and Vinokur [75] and by Falco et al. [80], but have been contradicted<br />

by Yukalov et al. [82].<br />

94

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

Saved successfully!

Ooh no, something went wrong!