The evaluation of non-local pseudopotentials within the QMC solid codes is dealt with in detail in chapter . Optimising wavefunctions with respect to these non-local pseudopotentials is discussed in chapter . In chapters and , DMC calculations involving a non-local pseudopotential to describe silicon are reported. The use of such non-local pseudopotentials within DMC is a problem that is yet to be completely resolved. Here we discuss the approximations used later in the thesis to try and deal with this problem.
Consider Eq.(), in which a pseudopotential is used to describe the valence electrons in a system;
We can divide the full valence Hamiltonian, , into a local and non-local part
where includes the kinetic energy, the local part of the pseudopotential and the Coulomb interaction, and W includes the non-local part of the pseudopotential. The diffusion-drift equation for f, Eq.(), can be split into two terms representing the local and non-local parts of the Hamiltonian
The first three terms on the left hand side of Eq.() can be interpreted as a local diffusion, drifting and branching process. However, the fourth term represents the operation of the non-local pseudopotential on the unknown wavefunction, , producing non-local branching. In the DMC calculations described in chapters and , we have used the ``locality approximation'' which was introduced by Christiansen[38, 39, 40] and recently applied by Mitas[41, 42] to the problem of the unknown wavefunction in Eq.(). In this approximation, the non-local pseudopotential acts not on the unknown wavefunction, , but on the guiding wavefunction, ,
One is then free to impose the fixed-node approximation in the same way as before. In fact, when using the above local model potential, it is important that the fixed node approximation is applied as this steers the random walk away from the nodes of , where there will generally be divergences in the local model potential. Without the fixed node approximation, these divergences would cause large fluctuations in the population of diffusing particles.
It should be noted that within the ``locality approximation'', it is no longer generally true that the DMC estimate of the energy is an upper bound to the true groundstate energy. It must always be less than the VMC energy and it has been shown, [42] that the DMC energy converges quadratically to the exact groundstate energy as the guiding function approaches the true groundstate wavefunction.