Thermodynamics, Diffusion, and Structure of Mg2SiO4 forsterite grain boundaries from atomistic modeling

Most materials are not single crystals but consist of crystalline grains of various sizes, misorientated with respect to each other and joint by grain boundaries. The latter influence many of the material properties. For instance, grain boundaries are short circuits for diffusion and thus they strongly influence transport properties of materials such as electrical conductivity, or mineral growth rates, creep, or phase transform. Olivine is a major component of the Earth's upper mantle and therefore it is of considerable importance to study its physical and thermodynamic polycrystalline properties. In the present study, we have used molecular dynamics simulations to model thermodynamics, self-diffusion and structure of a series of [100] symmetric tilt grain boundaries in forsterite. The interactions between the atoms are modeled by an advanced ionic interaction potential (Jahn and Madden, 2007). The parameters of the potential are fitted to ab initio results. The model was optimized for the Ca-Mg-Al-Si-O system and shows good transferability in a wide range of pressures, temperatures, and compositions. Thermodynamics and structure were simulated at ambient conditions, and self-diffusion coefficients were determined at ambient pressure and temperatures of 1250, 1500, 1750, and 2000 K. We find that the energy and excess free volume of the grain boundaries in forsterite depend significantly on the misorientation angle of the grain boundary. One of our modeled structures is compared with an high-resolution transmission electron micrograph (HRTEM) (Heinemann et al., 2005). We relate our findings to previous studies of grain boundaries in ionic materials and in metals. For small misorientation angles (up to 22.1°), grain boundary structures consist of an array of c-edge dislocations as suggested by Heinemann et al. (2005) and their energies can be readily fit with the Read-Shockley dislocation model for grain boundaries. For high misorientation angles (32.1° and 60.8°), the cores of dislocations overlap and form repeated structural units. Similar to energies and excess free volumes, the self-diffusion coefficients of Mg and O depend significantly on the misorientation angle of the grain boundaries and they are well fitted with Arrhenius law. We compare our results to MgO grain boundary diffusion in forsterite derived from reaction rim growth experiments (Gardés and Heinrich, 2010).