The nonlocality of the pseudopotential only extends over the region occupied by the core of the atom. As the core region is relatively small, it should be possible to deal efficiently with the nonlocality of the pseudopotential by working in real space. Only a small number of operations should be required to project the angular momentum components of each wavefunction in the core of each atom. Furthermore, the number of operations needed to project the angular momentum components of a single wavefunction around a single atom in real space is independent of the size of the system, leading to better scaling with system size than in the reciprocal space projection scheme. The computational cost of the real-space implementation increases with the square of the number of atoms, compared with the cubic relationship that exists for reciprocal space projection methods (Payne et al., 1992).
CASTEP uses the scheme developed originally for norm-conserving potentials (King-Smith et al., 1991) and subsequently generalized for the case of ultrasoft pseudopotentials. The reciprocal-to-real space transformation involves the following steps:
Eq. CASTEP 22
where Vl(r) is the original nonlocal projector and jl(qr) is the spherical Bessel function.
Eq. CASTEP 23
The object of this step is to achieve the best possible convergence of the modified potential in real space.
The modified real-space projectors are no longer identically zero outside Rc, although the minimization in Eq. CASTEP 23 ensures the best possible convergence compatible with the removal of the wrap around error. A quantitative measure of the error introduced by the transformation is given by:
Eq. CASTEP 24
This error function should be small (less than 0.001) for all wavevectors, q, and can be inspected in the CASTEP output file.
One should bear in mind that the use of real-space scheme represents an approximation and that more accurate results are obtained with the reciprocal-space version. The accuracy of the real-space transformation is controlled by the value of the cutoff radius, R0. The default set of these values in CASTEP is designed to be for accuracy rather than performance. It is advisable to carry out independent tests on small systems before embarking on large-scale calculations with drastically reduced R0 values.
The real space transformation of the pseudopotential depends on the cutoff energy for plane wave basis set, and thus it is not possible to precalculate real-space projectors in the same way as the reciprocal space potentials (which are provided in CASTEP database).