Abstract:Rate maps quantify the variability of surface rates during the dissolution of crystalline matter. Recently, highly spatially resolved rate maps revealed the existence of rhythmic pulses of the material flux from the crystal surface. The mechanism underpinning this behavior is not yet understood but the two potential influencing factors are surface-controlled or transport-controlled conditions that may govern the resulting pulsating reaction kinetics in the system. In this study, we apply two numerical methods to identify the dominating mechanism of pulsating dissolution. First, the influence of solute transport is simulated using a reactive transport model, which yields dissolution rates and concentration distribution due to the flow field, i.e., extrinsic reactivity. Second, the influence of intrinsic surface reactivity inherent to the material is simulated using kinetic Monte Carlo (KMC) simulations, where the distribution of reactive sites over time is observed on an atomic scale. Local dissolution rates can be reproduced with the reactive transport model, but no kinetically effective periodic changes of the concentration gradients can be observed and no new dissolution pulses are generated. Control via periodic changes in extrinsic reactivity is thus ruled out as a governing mechanism. In contrast, pulsating dissolution is clearly observed in the KMC simulation, leading to the conclusion that the self-assembly of varying reactive surface building blocks causes the pulsating dissolution. The simulation results suggest that etch pits generate regions with a high number of steps of single crystal layers and, consequently, with sites of increased reactivity at periodic intervals. These steps move over the crystal surface outward from the pit center and are followed by a region with a low number of surface steps. Only when the steps reach a certain spacing to the center new steps are generated at the transition between the hollow core and the etch pit. This self-assembly is observed as a fundamental behavior of crystalline dissolution in KMC models, without any additional parametrization. It thus represents a fundamental mechanism in crystal dissolution and should be considered for an accurate atomic understanding of crystal dissolution.

On the Mechanism of the Dissolution of Quartz and Silica in Aqueous Solutions

The Dissociated Amorphous Silica Surface: Model Development and Evaluation.

Structure of sodium sulphate aqueous solution/quartz interface: A molecular dynamics simulation

The Water-Amorphous Silica Interface: Analysis Of The Stern Layer And Surface Conduction

Dissolution of Crystallites: Surface Energetic Control and Size Effects.

Static and Dynamic Properties of the Electrical Double Layer Near Amorphous Silica: Relevance for Device Design

A New Understanding of Demineralization: the Dynamics of Brushite Dissolution

Mechanism of Dissolution of Sparingly Soluble Electrolytes.

Crumbling Crystals: On the Dissolution Mechanism of NaCl in Water

Stepwise dissolution of silica surface in alkaline solution revealed by molecular modeling

Investigating the quartz (1010)/water interface using classical and ab initio molecular dynamics

Analysis of Spherical Crystal Dissolution in the Solution

Ab initio mechanism revealing for tricalcium silicate dissolution

Unravelling the dissolution dynamics of silicate minerals by deep learning molecular dynamics simulation: A case of dicalcium silicate

A kinetic Monte Carlo study of the C 3 S dissolution mechanism

Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach

Pulsating dissolution of crystalline matter: A surface-controlled process

Mechanism for the effects of surface chemical composition and crystal face on the wettability of α-quartz surface

Crystal face dependent wettability of α-quartz: Elucidation by time-of-flight secondary ion mass spectrometry techniques combined with molecular dynamics.