Next: Charge density
Up: Total energy using the
Previous: Kinetic energy
For clarity of notation, we first rewrite the nonlocal pseudopotential of equation (14) as

(28) 
which defines the projectors
.
Robust real space methods to calculate the nonlocal pseudopotential energy exist in the context of traditional planewave DFT through the work of KingSmith et al. [19]. We manage to avoid entirely the complications involved in their method through use of the FFT box.
We need to calculate matrix elements such as
which is just a matter of computing quantities like
. There will only be a contribution to a particular matrix element
if there is at least one projector that overlaps with both
and
(see Figure 3). We begin with the radial part of the projectors,
, on a reciprocal space, linear, radial grid, up to arbitrarily high wavevector. Given that the overlap condition with the NGWFs is satisfied, we then continue to calculate each overlap in turn. For example, to compute
, we envisage a real space FFT box that contains the NGWF
and the projector
. The reciprocal representation of the projector is interpolated onto the corresponding reciprocal space FFT box using

(30) 
where
is the position vector of atom , and are real spherical harmonics.
Figure 3:
An example of a typical contribution to the nonlocal potential matrix element
. The overlap of the projector, , with
is calculated using the FFT box with solid outline, and that with
is done using the FFT box with dashed outline.

Performing a FFT on the projector, we obtain it in the real space FFT box, project it into the simulation cell, and take the overlap with the NGWF by summation over the grid points that lie within the localisation region of
. In terms of the projection operator, , this process is represented by

(31) 
Since the FFTs are performed only on a restricted region of the simulation cell, namely the FFT box, and because the NGWFs and atom cores are strictly localised, the cost of calculation of the nonlocal pseudopotential matrix in this way scales as .
Next: Charge density
Up: Total energy using the
Previous: Kinetic energy
Peter D. Haynes
20021029