next up previous contents
Next: Choice of Parameters Up: Optimising the function Previous: Optimising the function

Which Part of tex2html_wrap_inline6049 to Optimise?

As described in chapter gif, the many-body wavefunction can be written in the following form;

  equation2691

The trial wavefunction in Eq.(gif) is commonly referred to as the Hartree-Fock-Jastrow-Chi wavefunction. All three parts of the wavefunction, i.e. the HF Slater determinant of one-electron orbitals, the Jastrow factor, and the Chi function, would benefit from some form of optimisation of parameters, and so deciding which part to optimise is a matter of deciding which part should yield the greatest improvement in the quality of tex2html_wrap_inline6049 .

The Chi function has previously been constructed according to a scheme introduced by Fahy [26], using Eq.(gif). Of the three parts of the Hartree-Fock-Jastrow-Chi trial wavefunction, this function has the weakest theoretical justification for its form. It is constructed on an ad hoc basis by making the assumption that the LDA charge density is reasonably close to the real charge density. The function is then constructed according to

  equation2707

where tex2html_wrap_inline7065 is the density produced by a VMC calculation using a trial wavefunction, tex2html_wrap_inline6049 , with tex2html_wrap_inline7069 . The overall trial wavefunction should then reproduce a charge density close to the LDA form. Equation (gif) certainly does produce a ţhat is able to reduce the energy of the system (compared to a tex2html_wrap_inline6049 with tex2html_wrap_inline7069 , see figure gif) and produce a charge density that is a reasonable reflection of the LDA charge density (see figure gif). However, there is nothing to suggest that a function constructed using Eq.(gif) is the best possible function either in terms of total energy or charge density. Indeed, there is little hard evidence [66] of exactly how accurate the LDA charge density is. It therefore seemed sensible to start optimising tex2html_wrap_inline6049 by optimising .

  
Figure: Charge density along the Ge-Ge bond for different functions.


next up previous contents
Next: Choice of Parameters Up: Optimising the function Previous: Optimising the function

Andrew Williamson
Tue Nov 19 17:11:34 GMT 1996