Crystal–melt interfaces in Mg2SiO4 at high pressure: structural and energetics insights from first-principles simulations

The interplay between crystal–melt and grain boundary interfaces in partially melted polycrystalline aggregates controls many physical properties of mantle rocks. To understand this process at the fundamental level requires improved knowledge about the interfacial structures and energetics. Here, we report the results of first-principles molecular dynamics simulations of two grain boundaries of (0l1)/[100] type for tilt angles of 30.4° and 49.6° and the corresponding solid–liquid interfaces in Mg₂SiO₄ forsterite at the conditions of the upper mantle. Our analysis of the simulated position time series shows that structural distortions at the solid–liquid interfacial region are stronger than intergranular interfacial distortions. The calculated formation enthalpy of the solid–solid interfaces increases nearly linearly from 1.0 to 1.4 J/m² for the 30.4° tilt and from 0.8 to 1.0 J/m² for the 49.6° tilt with pressure from 0 to 16 GPa at 1500 K, being consistent with the experimental data. The solid–liquid interfacial enthalpy takes comparable values in the range 0.9 to 1.5 J/m² over similar pressure interval. The dihedral angle of the forsterite–melt system estimated using these interfacial enthalpies takes values in the range of 67° to 146°, showing a decreasing trend with pressure. The predicted dihedral angle is found to be generally larger than the measured data for silicate systems, probably caused by compositional differences between the simulation and the measurements.