The energies obtained were eV (LDA), eV (VMC with LDA orbitals), eV (VMC with natural orbitals), eV (DMC with LDA orbitals), and eV (DMC with natural orbitals). The VMC wave function appears to show a very slight improvement with natural orbitals compared with LDA orbitals. However, to within statistical accuracy, the DMC energies obtained with LDA and natural orbitals are the same. This indicates that the nodal surfaces given by the LDA and natural orbitals are of the same quality.
This was investigated by performing DMC calculations in silicon with an fcc simulation cell containing 16 atoms. The smaller simulation cell enabled a large number of independent configurations to be obtained rapidly. Wave functions expanded in a basis of atom-centered Gaussians were obtained from the HF and DFT code CRYSTAL95. [130] Special care was taken to ensure that the LDA and HF calculations were done in equivalent ways to try and eliminate any bias in the comparison. A basis set of four uncontracted functions and one polarization function per pseudo-atom was optimized separately for each calculation. The quality of the basis set is high - to obtain the same energy within a plane wave calculation would require a basis set cutoff of 12.5 Ry. The same non-local LDA pseudopotential was used as in the earlier calculations. The Ewald interaction was used in the many-body Hamiltonian to avoid any charge-density dependence on the interaction. Approximately walker moves obtained DMC total energies and statistical accuracies of -107.488(3) eV per atom and -107.464(3) eV per atom for the LDA and HF guiding wave functions respectively.
The walker energies were approximately normally distributed. Using a conventional students t-test, [19] the 95% confidence interval on the difference in energies obtained was eV per atom, showing that for this system it is very likely that the DMC energy from a determinant of LDA orbitals is lower than that from a determinant of HF orbitals. Therefore, for this system, a determinant of LDA orbitals has a marginally better nodal structure than a determinant of HF orbitals.