Interaction of Aromatic Molecules with Forsterite: Accuracy of the Periodic DFT-D4 Method

Unlocking the surface chemistry of ionic minerals: a high-throughput pipeline for modeling realistic interfaces

A systematic procedure is introduced for modeling charge-neutral non-polar surfaces of ionic minerals containing polyatomic anions. By integrating distance- and charge-based clustering to identify chemical species within the mineral bulk, our pipeline, PolyCleaver, renders a variety of theoretically viable surface terminations. As a demonstrative example, this approach was applied to forsterite (Mg2SiO4), unveiling a rich interface landscape based on interactions with formaldehyde, a relevant multifaceted molecule, and more particularly in prebiotic chemistry. This high-throughput method, going beyond techniques traditionally applied in the modeling of minerals, offers new insights into the potential catalytic properties of diverse surfaces, enabling a broader exploration of synthetic pathways in complex mineral systems.

Adsorption of Polycyclic Aromatic Hydrocarbons and C60 onto Forsterite: C-H Bond Activation by the Schottky Vacancy

Understanding how to catalytically break the C-H bond of aromatic molecules, such as polycyclic aromatic hydrocarbons (PAHs), is currently a big challenge and a subject of study in catalysis, astrochemistry, and planetary science. In the latter, the study of the breakdown reaction of PAHs on mineral surfaces is important to understand if PAHs are linked to prebiotic molecules in regions of star and planet formation. In this work, we employed a periodic density functional theory along with Grimme's D4 (DFT-D4) approach for studying the adsorption of a sample of PAHs (naphthalene, anthracene, fluoranthene, pyrene, coronene, and benzocoronene) and fullerene on the [010] forsterite surface and its defective surfaces (Fe-doped and Ni-doped surfaces and a MgO-Schottky vacancy) for their implications in catalysis and astrochemistry. On the basis of structural and binding energy analysis, large PAHs and fullerene present stronger adsorption on the pristine, Fe-doped, and Ni-doped forsterite surfaces than small PAHs. On a MgO-Schottky vacancy, parallel adsorption of the PAH leads to the chemisorption process (C-Si and/or C-O bonds), whereas perpendicular orientation of the PAH leads to the catalytic breaking of the aromatic C-H bond via a barrierless reaction. Spin density and charge analysis show that C-H dissociation is promoted by electron donation from the vacancy to the PAH. As a result of the undercoordinated Si and O atoms, the vacancy acts as a Frustrated Lewis Pair (FLP) catalyst. Therefore, a MgO-Schottky vacancy [010] forsterite surface proved to have potential catalytic activity for the activation of C-H bond in aromatic molecules.