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.1 The Jastrow factor<br />
In section 2.3 it was shown that a single-determinant wave function takes account<br />
of exchange but not of correlation; the Jastrow factor allows correlation effects to<br />
be incorporated.<br />
The most important correlations are those involving pairs of electrons. These<br />
are included by having a term of the form<br />
− ∑ i>j<br />
u σi ,σ j<br />
(|r i − r j |) (3.78)<br />
in the Jastrow exponent J(X). Recall that the single-determinant wave function<br />
does nothing to prevent electrons of opposite spin from coming together; this term<br />
keeps these electrons apart, resulting in a significant lowering of energy. Electrons<br />
of like spin are also kept apart more than before, although this affects the energy<br />
less dramatically.<br />
The two-body term of equation (3.78) does not simply keep electrons apart. Both<br />
the long- and short-range behaviour of u are constrained by theoretical arguments.<br />
When two electrons approach each other, the Coulomb energy diverges; for a<br />
wave function to be an eigenstate of Ĥ, this divergence must be cancelled by a<br />
corresponding divergence in the kinetic energy. Such a divergence is produced by<br />
cusps in the wave function: discontinuities in the first derivative with respect to the<br />
distance between the electrons. A full discussion of the cusp conditions is given in<br />
appendix B.<br />
The long-range behaviour of u may be determined by arguments based on the<br />
random phase approximation of Bohm and Pines [9], and is the subject of chapter 7.<br />
A connection is made between the long-range electron-electron correlations and the<br />
long-wavelength density fluctuations known as plasmons; for a homogeneous system,<br />
the resulting u function has the form 1/ω p |r i − r j | in the limit |r i − r j | → ∞, where<br />
ω p = √ 4πn is the plasma frequency.<br />
A function which combines the required short- and long-range behaviour is<br />
u σi σ j<br />
(|r i − r j |) =<br />
1<br />
)<br />
(1 − e −|r i−r j |/F σi σ j<br />
, (3.79)<br />
ω p |r i − r j |<br />
53