Atomic structures and energies of grain boundaries in Mg2SiO4 forsterite from atomistic modeling

Grain boundaries influence many physical and chemical properties of crystalline materials. Here, we perform molecular dynamics simulations to study the structure of a series of [100] symmetric tilt grain boundaries in Mg2SiO4 forsterite. The present results show that grain boundary energies depend significantly on misorientation angle. For small misorientation angles (up to 22A degrees), grain boundary structures consist of an array of partial edge dislocations with Burgers vector associated with stacking faults and their energies can be readily fit with a model which adds the Peach-Koehler equation to the Read-Shockley dislocation model for grain boundaries. The core radius of partial dislocations and the spacing between the partials derived from grain boundary energies show that the transition from low- to high-angle grain boundaries occurs for a misorientation angle between 22A degrees and 32A degrees. For high misorientation angles (32.1A degrees and 60.8A degrees), the cores of dislocations overlap and form repeated structural units. Finally, we use a low energy atomic configuration obtained by molecular dynamics for the misorientation of 12.18A degrees as input to simulate a high-resolution transmission electron microscopy (HRTEM) image. The simulated image is in good agreement with an observed HRTEM image, which indicates the power of the present approach to predict realistic atomic structures of grain boundaries in complex silicates.