14.09.2014 Views

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

CASINO manual - Theory of Condensed Matter

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!