Lattice Constants

The lattice constant of a solid, , corresponds to the size of the conventional unit cell length at the equilibrium volume, and is obtained computationally by minimising the total energy as a function of cell volume. Experimental lattice constants are usually obtained from low temperature X-ray diffraction measurements, and extrapolated to zero Kelvin.

The calculation of theoretical lattice constants is straightforward
for cubic systems - single-point energy calculations are performed
at several different volumes using the same kinetic energy cutoff and
-point sampling, and the results are fitted to the Murnaghan equation
of state [110],

Here is the equilibrium bulk modulus which is defined in Sec. 3.3.2, and is the derivative of the bulk modulus. Usually to calculated energy-volume points that span the equilibrium volume by approximately are sufficient for an accurate determination of . An example of a Murnaghan interpolation is shown in Fig. 3.2 for Si obtained using the HCTH functional.

Table 3.1 presents the calculated lattice constants obtained with the LDA, PW91 and HCTH functionals. As expected, the LDA uniformly underestimates with a mean error of 0.06 Å, while PW91 uniformly overestimates them, with a mean error of 0.03 Å. HCTH performs worse than even the LDA with a systematic overestimation of 0.08 Å. For HCTH there is a clear correlation between the accuracy of the lattice constant and the number of occurrences of the constituent atoms in the G fitting data used to determine the functional. For C, Si and Ge, the lattice constant errors are 0.01 Å (), 0.07 Å (), and 0.14 Å () respectively - the number of systems in the fitting data containing carbon, silicon, and germanium atoms are 19, 7, and 0 respectively. The errors for the aluminium semiconductors AlN, AlP, and AlAs are Å (), Å (), and Å () with 8, 4, and 0 occurrences of nitrogen, phosphorus, and arsenic atoms in the fitting data. The errors for the gallium semiconductors GaN, GaP, and GaAs are Å (), Å (), and Å () respectively. The fitting data contains one system with aluminium but none with gallium, which is consistent with the lattice constants for GaN and GaAs being less accurate than AlN and AlAs.

The significant overestimation of lattice constants made by HCTH is
consistent with the findings of Kurth * et al.* [108].
Their unit cell volumes for Si, Ge, and GaAs correspond to lattice
constants of , , and Å respectively, which are
close to the values shown in Table 3.1. It should be noted
that the degradation in performance for heavier systems is also evident
in molecular calculations [98]. A subset of the molecules in the
HCTH fitting data have well-known experimental data. The subset includes
28 molecules containing just first-row atoms and 12 containing second-row
atoms. HCTH overestimates the bond lengths of these two sets of systems
by an average of 0.008 Å and 0.025 Å respectively [98].