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.
LIST OF FIGURES<br />
5.3 The LDA energy per electron as a function of the number of grid<br />
points. Cells of size w ∼105 and w ∼210 are compared, both for<br />
finite systems containing 360 electrons and for systems infinite in the<br />
xy-direction. Note that the two vertical scales are offset. The slab<br />
width is 17.64248; the density parameter is r s = 2.07. The length of<br />
the cell in the z-direction varies slightly as a function of the number<br />
of grid points; this is to ensure optimal sampling of the important<br />
slab region. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77<br />
5.4 The dependence of the LDA energy per electron on the slab width;<br />
for this system, r s = 2.07, the number of electrons is 360 and the cell<br />
size in the z-direction is around 105. . . . . . . . . . . . . . . . . . . 78<br />
5.5 The surface energy as a function of slab width, calculated using different<br />
values of the bulk energy. The middle curve uses the Perdew-<br />
Wang parameterisation [69] of Ceperley and Alder’s early QMC results<br />
[12], which was the functional employed in the slab calculation.<br />
The upper curve is from the Perdew-Zunger [70] paremeterisation of<br />
the same results; the lower curve is from the Perdew-Zunger parameterisation<br />
of Ortiz and Ballone’s later QMC calculations [63]. . . . . 79<br />
5.6 The components of the LDA surface energy as a function of the slab<br />
width. All the surface energies are given in mHa bohr −2 . . . . . . . . 80<br />
5.7 The distribution of configuration energies, before and after variance<br />
minimisation. 10000 configurations have been included; the Jastrow<br />
factor was of the form introduced in section 3.4.1 and described in<br />
detail in section 6.3. Optimisation of both u and χ terms was allowed;<br />
The minimisation was performed without reweighting. . . . . . . . . . 81<br />
6.1 The total VMC energy per electron as a function of system size,<br />
comparing the Ewald and MPC interactions. The exchange energy<br />
given by equation (6.25) is also plotted (note the different vertical<br />
scale). Energies, as usual, are in Ha. . . . . . . . . . . . . . . . . . . 90<br />
9