next up previous contents index
Next: Non-selfconsistent Harris-Foulkes functional Up: Theoretical Background Previous: Residual minimization scheme, direct   Contents   Index

N.B. This document is no longer maintained, please visit our wiki.

Wrap-around errors -- convolutions

In this section we will discuss wrap around errors. Wrap around errors arise if the FFT meshes are not sufficiently large. It can be shown that no errors exist if the FFT meshes contain all $ G$ vectors up to $ 2 G_{\rm cut}$.

It can be shown that the charge density contains components up to $ 2 G_{\rm cut}$, where $ 2 G_{\rm cut}$ is the 'longest plane' wave in the basis set:

The wavefunction is defined as

$\displaystyle \vspace*{1mm}
\vert \phi_{n{\bf k}}\rangle = \sum_G C_{{\bf G}n{\bf k}}\vert {\bf k}+{\bf G}\rangle,

in real space it is given by

$\displaystyle \vspace*{1mm}
\langle{\bf r}\vert \phi_{n{\bf k}}\rangle = \sum_...
...{\Omega^{1/2}} \sum_{\bf G}e^{i({\bf k}+{\bf G}){\bf r}} C_{{\bf G}n{\bf k}}.

Using Fast Fourier transformations one can define

$\displaystyle \vspace*{1mm} C_{{\bf r}n{\bf k}}= \sum_{{\bf G}} C_{{\bf G}n{\bf...
...= \frac{1}{N_{\rm FFT}} \sum_{{\bf r}} C_{{\bf r}n{\bf k}}e^{-i{\bf G}{\bf r}}.$ (7.3)

Therefore the wavefunction can be written in real space as

$\displaystyle \vspace*{1mm} \langle{\bf r}\vert \phi_{n{\bf k}}\rangle = \phi_{n{\bf k}}(r) = \frac{1}{\Omega^{1/2}} C_{{\bf r}n{\bf k}}e^{i{\bf k}{\bf r}}.$ (7.4)

The charge density is simply given by

$\displaystyle \vspace*{1mm} \rho^{\rm ps}_{\bf r}\equiv \langle{\bf r}\vert \rh...
...um_{\bf k}w_{\bf k}\sum_n f_{n{\bf k}}\phi_{n{\bf k}}(r) \phi^*_{n{\bf k}}(r) ,$ (7.5)

in the reciprocal mesh it can be written as

$\displaystyle \vspace*{1mm} \rho^{\rm ps}_{\bf G}\equiv \frac{1}{\Omega} \int \...
...frac{1}{N_{\rm FFT}} \sum_{{\bf r}} \rho^{\rm ps}_{\bf r}e^{-i {\bf G}{\bf r}}.$ (7.6)

Inserting $ \rho^{\rm ps}$ from equation (7.5) and $ C_{{\bf r}n{\bf k}}$ from (7.3) it is very easy to show that $ \rho^{\rm ps}_{\bf r}$ contains Fourier-components up to $ 2 G_{\rm cut}$.

Generally it can be shown that a the convolution $ f_r=f^1_r f^2_r$ of two 'functions' $ f^1_r$ with Fourier-components up to $ G_1$ and $ f^2_r$ with Fourier-components up to $ G_2$ contains Fourier-components up to $ G_1+G_2$.

The property of the convolution comes once again into play, when the action of the Hamiltonian onto a wavefunction is calculated. The action of the local-potential is given by

$\displaystyle a_{\bf r}= V_{\bf r}C_{{\bf r}n{\bf k}}

Only the components $ a_{\bf G}$ with $ \vert {\bf G}\vert < G_{\rm cut}$ are taken into account (see section 7.1: $ a_{\bf G}$ is added to the wavefunction during the iterative refinement of the wavefunctions $ C_{{\bf G}n{\bf k}}$, and $ C_{{\bf G}n{\bf k}}$ contains only components up to $ G_{\rm cut}$). From the previous theorem we see that $ a_{\bf r}$ contains components up to $ 3 G_{\rm cut}$ ($ V_{\bf r}$ contains components up to $ 2 G_{\rm cut}$).
Figure 5: The small sphere contains all plane waves included in the basis set $ G<G_{\rm cut}$. The charge density contains components up to $ 2 G_{\rm cut}$ (second sphere), and the acceleration $ a$ components up to $ 3 G_{\rm cut}$, which are reflected in (third sphere) because of the finite size of the FFT-mesh. Nevertheless the components $ a_{\bf G}$ with $ \vert {\bf G}\vert < G_{\rm cut}$ are correct i.e. the small sphere does not intersect with the third large sphere
If the FFT-mesh contains all components up to $ 2 G_{\rm cut}$ the resulting wrap-around error is once again 0. This can be easily seen in Fig. 5.

N.B. Requests for support are to be addressed to: