Effect of the environment on the hydroxyl density of α-quartz (111)

Generation of Amorphous Silica Surfaces with Controlled Roughness

Amorphous silica (a-SiO2) surfaces, when grafted with select metals on the active sites of the functionalized surfaces, can act as useful heterogeneous catalysts. From a molecular modeling perspective, one challenge has been generating a-SiO2 slab models with controllable surface roughness to facilitate the study of the effect of surface morphology on the material properties. Previous computational methods either generate relatively flat surfaces or periodically corrugated surfaces that do not mimic the full range of potential surface roughness of the amorphous silica material. In this work, we present a new method, inspired by the capillary fluctuation theory of interfaces, in which rough silica slabs are generated by cleaving a bulk amorphous sample using a cleaving plane with Fourier components randomly generated from a Gaussian distribution. The width of this Gaussian distribution (and thus the degree of surface roughness) can be tuned by varying the surface roughness parameter α. Using the van Beest, Kramer, and van Santen (BKS) force field, we create a large number of silica slabs using cleaving surfaces of varying roughness (α) and using two different system sizes. These surfaces are then characterized to determine their roughness (mean-squared displacement), density profile, and ring size distribution. This analysis shows a higher concentration of surface defects (under-/overcoordinated atoms and strained rings) as the surface roughness increases. To examine the effect of the roughness on surface reactivity, we re-equilibriate a subset of these slabs using the reactive force field ReaxFF and then expose the slabs to water and observe the formation of surface silanols. We observe that the rougher surfaces exhibit higher silanol concentrations as well as bimodal acidity.

Surface Structure and Stability of Partially Hydroxylated Silica Surfaces

Surface energies of silicates influence crack propagation during brittle fracture and decrease with surface relaxation caused by annealing and hydroxylation. Molecular-level simulations are particularly suited for the investigation of surface processes. In this work, classical MD simulations of silica surfaces are performed with two force fields (ClayFF and ReaxFF) to investigate the effect of force field reactivity on surface structure and energy as a function of surface hydroxylation. An unhydroxylated fracture surface energy of 5.1 J/m² is calculated with the ClayFF force field, and 2.0 J/m² is calculated for the ReaxFF force field. The ClayFF surface energies are consistent with the experimental results from double cantilever beam fracture tests (4.5 J/m²), whereas ReaxFF underestimated these surface energies. Surface relaxation via annealing and hydroxylation was performed by creating a low-energy equilibrium surface. Annealing condensed neighboring siloxane bonds increased the surface connectivity, and decreased the surface energies by 0.2 J/m² for ClayFF and 0.8 J/m² for ReaxFF. Posthydroxylation surface energies decreased further to 4.6 J/m² with the ClayFF force field and to 0.2 J/m² with the ReaxFF force field. Experimental equilibrium surface energies are ∼0.35 J/m², consistent with the ReaxFF force field. Although neither force field was capable of replicating both the fracture and equilibrium surface energies reported from experiment, each was consistent with one of these conditions. Therefore, future computational investigations that rely on accurate surface energy values should consider the surface state of the system and select the appropriate force field.