My PhD thesis - Condensed Matter Theory - Imperial College London
My PhD thesis - Condensed Matter Theory - Imperial College London
My PhD thesis - Condensed Matter Theory - Imperial College London
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 3.<br />
QUANTUM MONTE CARLO METHODS<br />
3.4 Trial wave functions<br />
In VMC, the quality of the trial wave function sets a limit on the accuracy of the<br />
calculation. In fixed-node DMC, this is also true to some extent: the nodes are<br />
determined by the trial wave function, which therefore also determines the fixednode<br />
error. However, the statistical efficiency of a DMC calculation strongly depends<br />
on the quality of the trial wave function, as this determines the importance sampling.<br />
In addition, from the computational point of view, evaluating the trial wave function<br />
normally constitutes the major part of the calculation. The form of the trial wave<br />
function is therefore very important; it must be both accurate and easy to evaluate.<br />
The wave functions normally used in QMC simulations are of the Slater-Jastrow<br />
type, consisting of a Slater determinant multiplied by an exponential Jastrow factor:<br />
Ψ(X) = e J(X) D(X), (3.76)<br />
where D(X) is a determinant of one-electron orbitals, exactly as in equation (2.13).<br />
In fact, the computation is made significantly faster by the use of wave functions<br />
of the form<br />
Ψ(X) = e J(X) D ↑ (R ↑ )D ↓ (R ↓ ). (3.77)<br />
The spin-dependence of the one-electron orbitals in the determinants has been removed,<br />
and the evaluation of the two smaller determinants D ↑ and D ↓ is more<br />
efficient than that of the large determinant D. This function is no longer antisymmetric<br />
on exchange of electrons with opposite spins; however, the expectation value<br />
of any spin-independent operator is unaffected by this alteration.<br />
The single-electron orbitals may be obtained from density-functional theory or<br />
Hartree-Fock calculations. The optimal orbitals in these two mean-field schemes<br />
are usually very similar. The nodal surface of the resulting trial wave function<br />
is completely defined by these orbitals, since the Jastrow factor is never zero; for<br />
DMC calculations, this means the expectation value of the energy is not affected by<br />
including the Jastrow factor. However, the variance of the energy is affected, and<br />
the Jastrow factor is an important part of the trial wave function.<br />
52