Molecular dynamics study of orientational order and rotational melting in clusters of TeF6

Molecular dynamics simulations of the behavior of molecules in crystalline clusters of TeF₆ were carried out on systems of 100, 150, 250, and 350 molecules. Several diagnostic functions were applied to investigate whether rotational melting occurred before translational melting. These functions included the coefficient of rotational diffusionD _θ(T), the “orientational Lindemann index” δ_θ(T), the “orientational angular distribution function”Q(θ,T), and the “orientational pair-correlation function”g _θ(r, T). All indicators implied that rotational melting occurred before translational melting, that it began with the outermost molecules, and that its onset for smaller clusters was at lower temperatures than for larger clusters. Results also showed that the rotational transition coincided with the transition from a lower symmetry phase (monoclinic) to cubic, a phenomenon that had been noted by others to occur with some regularity for systems of globular molecules.