To determine the optimal choice of for diamond-structure germanium, a series of LDA calculations were performed using -point meshes with different offsets from the origin[33]. The offset from the origin determines the value of of the Slater determinant constructed from the orbitals. The results of these calculations are illustrated in figure . The results show the convergence of the total energy with supercell size for four different values of . For the -point ( =0) the energy converges slowly as the size of the supercell is increased. A far more rapid convergence is observed when L-point sampling is used. This means the -point grid is offset from the origin by , where and the are the primitive reciprocal lattice translation vectors of the supercell. The best convergence is obtained when using an offset equal to the Baldereschi mean value point[53] of the supercell BZ.The Baldereschi mean point was not chosen for the QMC calculations, as wavefunctions at that -point are necessarily complex and the use of complex arithmetic would slow down the code considerably. It is possible to construct a real wavefunction using Baldereschi-point states by taking a linear combination of and its complex conjugate, , as both and have the same energy eigenvalue. Unfortunately, this combination is no longer a Bloch function, and this results in a local energy function that does not have the periodicity of the supercell.
Although it would be too computationally expensive to repeat all the LDA calculations from figure within the QMC formalism, VMC and DMC calculations for supercells with n = 2 and 3 [33] show that QMC calculations follow the LDA trend very closely. The calculations described in chapters and use Slater determinants where the one-electron orbitals have -points chosen on a mesh offset from the origin at the L-point. This corresponds to Monkhorst-Pack (MP)[54] sampling for n even supercells and better than MP sampling for n odd supercells.
Figure: Convergence of the total
energy with simulation cell size for different -point
sampling schemes. The black line shows the results of
sampling. The green line shows sampling, where the are
the primitive reciprocal lattice translation vectors of the
supercell. The red line shows sampling, and
the blue line shows sampling, where
is the Baldereschi mean value point.