Harmonic and Anharmonic Properties of Diamond Structure Crystals with Application to the Calculation of the Thermal Expansion of Silicon

Silicon has interesting harmonic and anharmonic properties such as the low lying transverse acoustic modes at the X and L points of the Brillouin zone, negative Gruneisen parameters, negative thermal expansion and anomalous acoustic attenuation. In an attempt to understand these properties, a lattice dynamical model employing long range, nonlocal, dipole-dipole interactions was developed. Analytic expression for the Gruneisen parameters of several modes are presented. These expressions explain how the negative Gruneisen parameters arise. This model is applied to the calculation of the thermal expansion of silicon from 5K to 1700K. The thermoelastic contribution to the acoustic attenuation of silicon is computed from 1 to 300 K. Strong attenuation anomalies associated with negative thermal expansion are found in the vicinity of 17K and 125K.