Tribochemical mechanism of amorphous silica asperities in aqueous environment: a reactive molecular dynamics study
Evidence and attribution¶
Authority of statements
Prose sections below (Summary, Methods, Findings, etc.) are curated summaries of the publication identified by doi, title, and pdf_path in the front matter above. They are not new primary claims by this wiki.
For definitive numerical values, reaction schemes, and interpretations, use the peer-reviewed article (and optional records under normalized/papers/ when present)—not this page alone.
Summary¶
This Langmuir article reports ReaxFF reactive MD of amorphous silica asperities sliding in liquid water versus phosphoric acid solution, with applied normal loads in the GPa range. In pure water, low normal load tends to keep surfaces from strongly adhering, whereas higher load promotes formation and rupture of interfacial siloxane bridges. In phosphoric acid, acid polymerization yields oligomers under milder contact conditions, while higher load drives more degradative tribochemical pathways that increase wear. The authors emphasize bridging oxygen atoms in the silica network as key to covalent bridge formation, with hydrogen migration weakening bridges and enabling shear-driven rupture.
Methods¶
LAMMPS ReaxFF (papers/Yue_Langmuir_Asperity_2015.pdf, §2.3) models a-SiO₂ asperity on a-SiO₂ substrate in water or phosphoric acid (Fig. 1, §2). 3D PBC cells use frozen substrate bases, rigid asperity drivers, and movable buffer layers flanking the reactive zone. 50 ps equilibration at 300 K precedes sliding. Velocity Verlet, 0.25 fs timestep. Langevin thermostats (100 fs damping) act only on the two buffer slabs (adjacent to the frozen base and below the rigid cap); the interior reacts with deterministic MD—a hybrid thermal setup rather than a single global NVT/NVE label. A uniform normal load on the cap (−z); examples 0.5 and 1.0 GPa (6.28 and 12.56 nN) align with an experimental window cited to Li et al. Lateral sliding at 10 m s⁻¹ (0.1 Å ps⁻¹, +x) records friction. Normal stress is imposed by forces, not a hydrostatic barostat. 300 K via buffer thermostats. No electric field or enhanced sampling.
Force-field training: N/A — established Si/O/H (acid-compatible) ReaxFF parametrizations from §2 citations.
Static QM / DFT: N/A — not the centerpiece.
Findings¶
Tribochemistry tracks friction and wear: bridging oxygens form interfacial covalent bridges while mobile H weakens them and enables shear-driven rupture (abstract, §3). In water, low load stays relatively passive / low adhesion; high load drives siloxane formation–rupture cycles. In phosphoric acid, oligomerization appears at lower load, whereas higher load favors more degradative chemistry and higher wear (abstract). Normal load (0.5 vs 1.0 GPa in the examples) is the main lever between passivation, polymerization, and wear. Discussion notes single-asperity idealization versus multi-contact statistics and third bodies. International work with van Duin-group ReaxFF support—use Langmuir figures for time series.
Limitations¶
- Single-asperity models omit multi-contact statistics and third-body particles from wear debris.
- Force-field limits on pH, ionic strength, and long-range polarization may affect absolute reaction rates.
Relevance to group¶
International collaboration (Tsinghua + Penn State) applying group ReaxFF expertise to nanotribology and silica–water interfaces.
Citations and evidence anchors¶
- Abstract and Sec. 1 in
papers/Yue_Langmuir_Asperity_2015.pdf; DOI:10.1021/la5042663.