Molecular dynamics simulation of twinning in devitrite, Na₂Ca₃Si₆O₁₆

Abstract

Reports of Type II twins are quite rare for most crystal structures. When they do occur, they are usually one of a number of possible twinning modes observed in a particular material. However, for the triclinic phase devitrite, Na₂Ca₃Si₆O₁₆, which nucleates from commercial soda−lime−silica float glass subjected to suitable heat treatments, the only reported twinning mode to date is a Type II twinning mode. In this study, this Type II twinning mode is first examined by molecular dynamics simulation to determine the lowest energy configuration of perfect twin boundaries for the twin mode. This is then compared with the lowest energy configurations of perfect twin boundaries found for six possible Type I twinning modes for devitrite for which the formal deformation twinning shear is less than 0.6. The most favourable twin plane configuration for the Type II twinning crystallography is shown to produce reasonably low twin boundary energies and sensible predictions for the optimum locations of the twin plane, K ₁, and the [1 0 0] rotation axis, η ₁, about which the 180° Type II twinning operation takes place. By comparison, all the Type I twinning modes were found to have very energetically unstable atomic configurations, and for each of these twinning modes, the lowest energy configurations found all led to high effective K ₁ twin boundary energies relative to perfect crystal. These results therefore provide a rationale for the experimental observation of the particular Type II twinning mode seen in devitrite.

Keywords:

• boundary
• devitrite
• molecular dynamics
• simulation
• twinning

Acknowledgements

We would like to thank the National University of Defense Technology at Changsha in China for the award of an overseas Academic Visitor scholarship to Dr Li during the course of this research.