In previous work [33, 50] VMC and DMC calculations were performed on diamond-structure germanium using the techniques described in chapter , but variance minimisation was not used to optimise the function in those calculations. The VMC calculations described in Refs.[33, 50] were repeated using the optimised functions described above. A fcc simulation cell of diamond structure germanium containing 16 atoms was studied. The single particle orbitals were obtained from a LDA calculation using a plane-wave basis set containing all plane waves up to a kinetic energy cutoff of 40 Ry. The calculations were performed at a non-zero wavevector of the simulation cell Brillouin zone using the techniques described in chapter and Refs.[33, 50].
The first set of calculations used a function containing 6 stars of vectors, whose coefficients had been optimised by the variance minimisation procedure, described in this chapter. The second set of optimisations of were performed again with 6 stars of vectors and with a further 6 parameters in the g(r) function (5 for to and B). 10,000 configurations were used in both the optimisations. The configurations for the second optimisation were generated from a VMC calculation using the best from the first optimisation. Values of L=3.0, 3.5, 4.0, 5.0 were used for the range of the g(r) function. The results yielding the lowest variance of the local energy are illustrated in figure .
Figure: Reduction in energy from optimising the function and adding g(r) to the function.
This figure shows that the introduction of a function according to the Fahy prescription, lowers the total energy by 0.15 eV per atom. Optimising the values of the coefficients in this function using the variance minimisation scheme lowers the energy by a further 0.17 eV per atom. The introduction of g(r) reduces the VMC energy by at most a further 0.05 eV per atom. The variance of the energy is only reduced by a few percent during the whole optimisation process.