Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well-understood. Harmonic phonons adapted to the specific volume at temperature, quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test the atomistic mechanisms, we performed inelastic
neutron scattering experiments on a single crystal of silicon to measure the changes in lattice dynamics
from 100 to 1500 K. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity, reproduced the measured shifts of individual phonons with temperature, whereas the quasiharmonic approximation typically gave results of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a fortuitous cancellation of
contributions from individual phonons
