Modeling hydroxylated nanosilica: Testing the performance of ReaxFF and FFSiOH force fields

We analyze the performance of the FFSiOH force field and two parameterisations of the ReaxFF force field for modeling hydroxylated nanoscale silica (SiO₂). Such nanosystems are fundamental in numerous aspects of geochemistry and astrochemistry and also play a key role during the hydrothermal synthesis of technologically important nanoporous silicas (e.g., catalysts, absorbents, and coatings). We consider four aspects: structure, relative energies, vibrational spectra, and hydroxylation energies, and compare the results with those from density functional calculations employing a newly defined dataset (HND: Hydroxylated Nanosilica Dataset). The HND consists of three sets of (SiO₂)₁₆(H₂O)_N nanoparticles (NPs), each with a different degree of hydroxylation and each containing between 23 and 26 distinct isomers and conformers. We also make all HND reference data openly available. We further consider hydroxylated silica NPs of composition (SiO₂)_M(H₂O)_N with M = 4, 8, 16, and 24 and infinite surface slabs of amorphous silica, both with variable hydroxylation. For energetics, both ReaxFF and FFSiOH perform well for NPs with an intermediate degree of hydroxylation. For increased hydroxylation, the performance of FFSiOH begins to significantly decline. Conversely, for the lower degree of hydroxylation both parameterisations of ReaxFF do not perform well. For vibrational frequencies, FFSiOH performs particularly well and significantly better than ReaxFF. This feature also opens the door to inexpensively calculating Gibbs free energies of the hydroxylated nanosilica systems in order to efficiently correct density functional theory calculated electronic energies. We also show how some small changes to FFSiOH could improve its performance for higher degrees of hydroxylation.