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.5. Exact diagonalization <strong>of</strong> the <strong>Bogoliubov</strong> problem<br />

used, instead <strong>of</strong> density δˆn and phase δ ˆϕ. This makes the structure better<br />

comparable with the Schrödinger equation and helps to point out the<br />

differences. The expansion <strong>of</strong> the grand canonical Hamiltonian (2.11) reads<br />

ˆF = 1<br />

�<br />

d<br />

2<br />

d r � δ ˆ Ψ † (r), δ ˆ Ψ(r) � � L G<br />

G∗ � �<br />

δΨ(r) ˆ<br />

L δ ˆ Ψ † �<br />

+ cst., (2.65)<br />

(r)<br />

with the Hermitian operator L = − �2<br />

2m∇2 + V (r) − µ + 2g|Φ(r)| 2 and G =<br />

gΦ(r) 2 . The goal is now to find a transformation to quasiparticles ˆ βν, such<br />

that the <strong>Bogoliubov</strong> Hamiltonian ˆ F takes the diagonal form<br />

ˆF = �<br />

�ων ˆ β † ν ˆ βν. (2.66)<br />

ν<br />

The <strong>Bogoliubov</strong> Hamiltonian was obtained by the expansion around the<br />

saddlepoint. It turns out that the saddlepoint is indeed an energy minimum,<br />

not only in the homogeneous problem (2.3.1), but also in inhomogeneous<br />

systems. Therefore, the spectrum �ων is positive. The transformation mixes<br />

δ ˆ Ψ(r) and δ ˆ Ψ † (r), similar to the homogeneous <strong>Bogoliubov</strong> transformation in<br />

subsection 2.3.1. Furthermore, the inhomogeneous potential and its density<br />

imprint will mix the plane waves, such that the eigenstates will have some<br />

unknown shape. The general transformation to a set <strong>of</strong> eigenmodes reads<br />

δ ˆ Ψ(r) = ��<br />

uν(r) ˆ βν − v ∗ ν(r) ˆ β † �<br />

ν , (2.67)<br />

ν<br />

with a set <strong>of</strong> initially unknown functions uν(r) and vν(r) [106]. Note that<br />

there are different conventions concerning the sign in front <strong>of</strong> v ∗ ν. Here and<br />

e.g. in [54, 105, 106], the minus sign is used, with vk(r) = vke ik·r in the<br />

homogeneous case. In other literature, vν(r) is defined with opposite sign<br />

[55, 79, 85, 86]. The quasiparticle annihilators ˆ βν appearing in equation<br />

(2.67) and their creators obey bosonic commutator relations, as will be<br />

verified below.<br />

The equation <strong>of</strong> motion for δ ˆ Ψ(r) is set up using both forms <strong>of</strong> the Hamiltonian<br />

(2.65) and (2.66). By equating the coefficients, we find the following<br />

equation on the uν(r) and vν(r)<br />

�<br />

L −G<br />

η<br />

−G∗ � �<br />

uν(r)<br />

L vν(r)<br />

� �� �<br />

=:Q<br />

�<br />

= �ων<br />

� uν(r)<br />

vν(r)<br />

�<br />

, (2.68)<br />

where η = diag(1, −1). This equation is known as the <strong>Bogoliubov</strong>-de-Gennes<br />

equation for the <strong>Bogoliubov</strong> eigenmodes [86]. Although both operators Q<br />

49

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

Saved successfully!

Ooh no, something went wrong!