24.10.2014 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!