CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
CASINO manual - Theory of Condensed Matter
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
and the Dirac delta function is<br />
δ(r − r ′ ) = 1 ∑<br />
e −iG·(r−r′) . (284)<br />
Ω<br />
The 3N-dimensional Fourier transform <strong>of</strong> the wave function is<br />
∫<br />
1<br />
˜Ψ(k 1 , . . . , k N ) =<br />
Ω N Ψ(r 1 , . . . , r N ) e ik1·r1 . . . e ik N ·r N<br />
dr 1 . . . dr N (285)<br />
∑<br />
Ψ(r 1 , . . . , r N ) = ˜Ψ(k1 , . . . , k N ) e −ik1·r1 . . . e −ik N ·r N<br />
. (286)<br />
k 1,...,k N<br />
G<br />
33.1.2 Translational averaging<br />
Consider a function f(r, r ′ ) whose Fourier transform is<br />
˜f(k, k ′ ) = 1 Ω 2 ∫<br />
f(r, r ′ ) e ik·r e ik′·r ′ dr dr ′ . (287)<br />
The translational average <strong>of</strong> f(r, r ′ ) is<br />
f T (r) = 1 Ω<br />
∫<br />
f(r ′′ , r ′′′ ) δ(r ′′ − r ′′′ − r) dr ′′ dr ′′′ . (288)<br />
The Fourier transform <strong>of</strong> the translational average is<br />
˜f T (k) = 1 ∫<br />
f T (r) e ik·r dr (289)<br />
Ω<br />
= 1 Ω 2 ˜f(k, −k). (290)<br />
Therefore the Fourier transform <strong>of</strong> f(r, r ′ ) and the Fourier transform <strong>of</strong> its translational average are<br />
related by<br />
˜f(k, −k) = Ω 2 ˜f T (k). (291)<br />
The operations <strong>of</strong> translational averaging and Fourier transforming commute.<br />
The function f(r, r ′ ) has the periodicity <strong>of</strong> the simulation cell,<br />
f(r, r ′ + R) = f(r, r ′ ), (292)<br />
where R is a translation vector <strong>of</strong> the simulation cell lattice. Equation (292) implies that f T (r) is<br />
also periodic,<br />
f T (r + R) = f T (r). (293)<br />
The Fourier transforms <strong>of</strong> f T (r) and f(r, r ′ ) are therefore nonzero only on the reciprocal lattice vectors<br />
<strong>of</strong> the simulation cell lattice.<br />
33.1.3 Rotational averaging<br />
The rotational average <strong>of</strong> a function is<br />
f(r) =<br />
∫<br />
1<br />
4πr 2<br />
f(r ′ ) δ(|r ′ | − r) dr ′ . (294)<br />
The rotationally averaged function f(r) is no longer periodic. The Fourier transform <strong>of</strong> the rotational<br />
average is<br />
˜f(k) = 1 ∫<br />
f(r) e ik·r dr (295)<br />
Ω<br />
= 1 Ω<br />
∫ ∞<br />
0<br />
f(r) 4πr 2 sin(kr)<br />
kr<br />
dr. (296)<br />
183