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.
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