Next: Appendix B: From the Up: Nonorthogonal generalized Wannier function Previous: acknowledgments

# Appendix A: Delta functions

In this paper, whenever we refer to delta functions' we will assume a periodic and bandwidth limited version of the Dirac delta functions. These delta functions are three-dimensional versions of the impulse functions' that are common in signal processing applications of FFTs [53]. In electronic structure, similar functions have been used as mesh delta functions' in the exact finite difference method' of Hoshi et al. [54] and in recent studies of their possible application when we consider the limit of an infinite simulation cell [55].

In our derivations we will assume that we have a simulation cell of any symmetry, which in general is a parallelepiped defined by its primitive lattice vectors , and . In this simulation cell we define a regular grid with an odd number of points , , and in every direction (the adaptation of our results to the case of even numbers of points is straightforward). Therefore point of this regular grid is defined as

 (31)

with , etc.

Bandwidth limited delta functions centred at points of the regular grid are defined as

 (32)

where , and are the primitive reciprocal lattice vectors of the simulation cell. Plane-waves whose wavevector is a linear combination of these reciprocal lattice vectors have periodicity compatible with the simulation cell and therefore so do our delta functions, or any other function expanded in terms of these plane-waves. These periodic bandwidth limited delta functions are our basis set. A plot of a two-dimensional version of one of these delta functions is shown in Figure 7. It is obvious from (32) that the delta functions are real-valued everywhere in space. They are not normalised to unity but they are normalised to the grid point volume ( is the volume of the simulation cell)
 (33)

Their value at grid points is equal to one when the grid point coincides with the centre of the function and zero for all other grid points
 (34)

The delta functions act as Dirac delta functions with the added effect of filtering out any plane-wave components that are not part of them. For example, if is a function periodic with the periodicity of the simulation cell but not bandwidth limited, it can be expressed in terms of its discrete Fourier transform (plane-wave) expansion

 (35)

where is the volume of the simulation cell.

It is straightforward to show that the projection of onto is

We define here to be the bandwidth limited version of the function , limited to the same frequency components as .

As the NGWFs are linear combinations of the delta functions according to (7), the result of equation (36) is very important since it leads to the following relation

 (36)

which means that the integral in the lefthand side of the above equation is exactly equal to a discrete summation of values on the grid, provided we use the bandwidth limited version of .

As a corollary we observe that the delta functions are an orthogonal set since

 (37)

We also need to define the fine grid delta functions (here the are just grid point indices for the fine grid, they are not related to any Cartesian coordinates). These functions are the analogues of the delta functions we have just described that would be obtained if we doubled the minimum and maximum values that their wavevectors can take. Consequently, they have the same periodicity but they correspond to a grid with twice the number of points in every direction, i.e. , and points. They are defined by

 (38)

As expected, the fine grid delta functions also satisfy an equation similar to (36)
 (39)

where is again a bandwidth limited version of but this time it is limited to contain any of the plane-waves that constitute rather than . It is easy to verify that any function that can be written as a sum of products of pairs of delta functions can also be written as a fine grid delta function expansion. We define and use the fine grid delta functions because of this `product rule' property.

Next: Appendix B: From the Up: Nonorthogonal generalized Wannier function Previous: acknowledgments
Peter D. Haynes 2002-10-31