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