Reaction pathways for quartz dissolution determined by statistical and graphical analysis of macroscopic experimental data

In light of recent work on the reactivity of specific sites on large (hydr)oxo-molecules and the evolution of surface topography during dissolution, we examined the ability to extract molecular-scale reaction pathways from macroscopic dissolution and surface charge measurements of powdered minerals using an approach that involved regression of multiple datasets and statistical graphical analysis of model fits. The test case (far-from-equilibrium quartz dissolution from 25 to 300 °C, pH 1-12, in solutions with [Na⁺] ? 0.5 M) avoids the objections to this goal raised in these recent studies. The strategy was used to assess several mechanistic rate laws, and was more powerful in distinguishing between models than the statistical approaches employed previously. The best-fit model included three mechanisms—two involving hydrolysis of Si centers by H₂O next to neutral (>Si-OH⁰) and deprotonated (>Si-O⁻) silanol groups, and one involving hydrolysis of Si centers by OH⁻. The model rate law is