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.
Keywords: PolyCleaver; Python; interfaces; mineral surfaces.
© Mates-Torres and Rimola 2024.