Proton-containing defects at forsterite {010} tilt grain boundaries and stepped surfaces

Atomistic simulation techniques are used to investigate the effect of proto n-containing defects on the structure and stabilities of a range of grain b oundaries of forsterite. We study two series of stepped {010} tilt boundari es that are at 90 degrees to each other: one with the {100} plane as step w all and the other with the {001} plane. Each series consist of several grai n boundaries with increasing terrace area (i.e., decreasing boundary angle) . The ratios of boundary and surface energies gamma(b)/gamma(s), and gamma( b)/gamma({010}) with boundary angle show maxima at a boundary angle phi = s imilar to 30 degrees and minima at phi = 0 degrees and similar to 60 degree s. The adhesion energies of the two series show a minimum at low boundary a ngle (phi = 20-30 degrees), indicating that there is an optimum size for th e (010) terrace area, where the relative stabilities of grain boundary and related surface make separation of the boundary into the free surfaces ener getically least expensive. Dissociative adsorption of water molecules in the bulk crystal is preferred at the M2 site, but is an endothermic process with a calculated hydration energy of +119 kJ/mol. Hydration of the grain boundaries on the other hand is energetically favorable with hydration energies tending toward -80 kJ/mo l compared with a surface hydration energy at the planar {010} surface of - 90 kJ/mol. We also investigated the adsorption of protons at cation defects , by modeling the process of replacement of Mg ions in the bulk and along t he grain boundaries by two H+ ions each. Replacement of Mg2+ in the bulk cr ystal occurs preferentially at the MI site and is calculated to be exotherm ic with a replacement energy of -78 kJ/mol; Mg2+ replacement along both ser ies of grain boundaries is also an exothermic process. The replacement ener gies tend to a constant value of approximately -210 kJ/mol with increasing terrace area. (cf., the calculated value for the planar {010} surface of -2 30.4 kJ/mol).