Interfacial behavior from the atomic blueprint: Machine learning-guided design of spatially functionalized α-SiO₂...

α-Quartz () surfaces functionalized with hydroxyl () and methyl () groups provide a versatile platform for controlling interfacial properties critical to applications such as catalysis, protective coatings, and energy conversion. The arrangement of these functional groups strongly influences interfacial interactions at solid-liquid interfaces, highlighting their relevance to colloid and interface science. However, conventional models often treat surface functionalization as spatially homogeneous, overlooking the atomic-scale organization of surface groups. We hypothesize that this spatial distribution, beyond overall composition, plays a decisive role in governing surface stability and interfacial behavior. To test this hypothesis, we employ a multi-scale simulation workflow combining density functional theory (DFT), ab initio molecular dynamics (AIMD), and machine-learned force fields (MLFFs). This approach allows us to explore a range of spatial patterns of / functionalization on the α-quartz (0001) surface. We evaluate the impact of spatial arrangements on mixing energy, hydrogen bonding networks, and vibrational properties with high accuracy and robustness. Our results reveal that spatial patterning strongly influences surface stability and interfacial structure. A thermodynamically favored unpaired configuration emerges near 67% substitution, where isolated groups form secondary hydrogen bonds through reorientation toward subsurface oxygen atoms. This rearrangement induces a characteristic blue shift in stretching frequencies, indicating weaker H-bonding. These effects are absent in clustered arrangements. By establishing a clear link between functional group patterning and interfacial behavior, our work uncovers the underlying mechanisms to guide and accelerate the rational design of silica-based materials and coatings, directly relevant to colloid and interface science.