Dynamics of interstitial hydrogen molecules in crystalline silicon

The static and dynamic properties of interstitial H 2 , HD and D 2 molecules in crystalline silicon are obtained from ab initio molecular-dynamics simulations with atomic-like basis sets. The static ( T = 0) calculations agree with those of most other authors: the centre of mass (CM) of H 2 is at the tetrahedral interstitial (T) site, the molecule is a nearly-free rotator, and the activation energy for diffusion is 0.90 eV. However, these results fail to explain a number of experimental observations, such as why H 2 is infrared (IR) active, why the expected ortho / para splitting is not present, why the symmetry is C 1 , why the piezospectroscopic tensors of H 2 and D 2 are identical or why the exposure to an H/D mix results in a single HD line which is not only at the wrong place but also much weaker than expected. In the present work, we extend the static calculations to include the constant-temperature dynamics for H 2 in Si. At T >0 K, the CM of the molecule no longer remains at the T site. Instead, H 2 `bounces' off the walls of its tetrahedral cage and exchanges energy with the host crystal. The average position of the CM is away from the T site along ⟨100⟩. Under uniaxial stress, the CM shifts off that axis and the molecule has C 1 symmetry. The H-H stretch frequency calculated from the Fourier transform of the v-v autocorrelation function is close to the measured one. Since the potential energy experienced by H 2 in Si near the T site is very flat, we argue that H 2 should be a nearly free quantum mechanical rotator. Up to room temperature, only the j = 0 and j = 1 rotational states are occupied, H 2 resembles a sphere rather than a dumbbell, the symmetry is determined by the position of the CM and HD is equivalent to DH in any symmetry. The rapid motion of the CM implies that an ortho-to-para transition will occur if a large magnetic moment is nearby. Several candidates are proposed. Since nuclear quantum effects are not included in our calculations, we cannot address the possibility that the observed vibrational spectrum of H 2 results from a tunnelling excitation as proposed by Stoneham.