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.

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

condensate gets already fragmented, due to rare high potential peaks from<br />

the exponential tail in (3.10), see also figure 4.8. This is reflected in the<br />

positive deviation <strong>of</strong> the mean values and also in the increased width <strong>of</strong> the<br />

distribution. To sum things up: Beyond-Born effects are clearly visible, but<br />

the Born result remains useful in a rather wide range. The correction stays<br />

negative, even for nearly fragmented condensates.<br />

4.3.3. Non-condensed fraction<br />

Even in the ground state, there are particles that are not in the Gross-<br />

Pitaevskii condensate function Φ, but in excited states (subsection 2.5.6).<br />

In order to verify the preconditions for treating the one-dimensional disordered<br />

problem on grounds <strong>of</strong> the Gross-Pitaevskii equation, we compute<br />

the fraction <strong>of</strong> non-condensed particles. We start with the homogeneous<br />

(quasi)condensate and evaluate the sum (2.87). In the thermodynamic limit,<br />

this sum diverges in one dimension. Nonetheless, we can compute the sum<br />

for systems <strong>of</strong> finite size. For the system size L = 200ξ, used in the previous<br />

subsections, we find nnc|L=200ξ = 0.3488/ξ. Note that the non-condensed<br />

fraction diverges logarithmically with L: for L = 105ξ, for example, we find<br />

nnc|L=10 5 ξ<br />

= 1.046/ξ.<br />

We need to relate the non-condensed density to the total density n1D<br />

<strong>of</strong> the <strong>Bose</strong> gas. This non-condensed fraction scales like (n1Dξ) −1 , i.e. the<br />

average particle spacing in relation to the healing length [55, chapter 17].<br />

The numerical procedure for computing the non-condensed density in<br />

presence <strong>of</strong> disorder is as follows. In the same manner as in the Gross-<br />

Pitaevskii integration (subsection 4.3.1), the system is discretized into l<br />

points, the disorder potential is generated, and the Gross-Pitaevskii ground<br />

state is computed. Then, the 2l × 2l real-space matrix (2.68) is set up and<br />

diagonalized numerically. As the matrix is not symmetric, the eigenvectors<br />

(uν(r), vν(r)) are not pairwisely orthogonal. Instead, they satisfy the biorthogonality<br />

relation (2.74). The non-condensed density (2.85) is then<br />

evaluated from the relevant modes with positive frequency (section 2.5). In<br />

figure 4.11, results for the parameters used in the previous subsections are<br />

shown. The initial increase is quadratic in the reduced disorder strength<br />

V0/µ. Strong deviations occur at V0 � 0.1µ. At the disorder strength<br />

V0/µ = 0.1, the fraction <strong>of</strong> non-condensed atoms is increased by about<br />

25%, i.e. the non-condensed fraction is still <strong>of</strong> the same order as in the<br />

homogeneous case.<br />

In conclusion, the 1D validity condition for Gross-Pitaevskii theory in<br />

the homogeneous case is nnc/n1D ∝ (ξn1D) −1 ≪ 1 [55, chapter 17] (with a<br />

100

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

Saved successfully!

Ooh no, something went wrong!