“…Due to the discontinuity of the bulk α-quartz structure at the quartz-water interface, the under-coordinated surface oxygen (O) and silicon (Si) atoms are capped with hydrogen (H) or hydroxyl (OH) groups as in previous studies. 30,56,57 Simulations are conducted in the NVT ensemble with a 1.0 fs time step, and the long-range coulombic (electrostatic) interactions are resolved by the particle-particle and particle-mesh (PPPM) method 58 and temperature is controlled by a canonical sampling thermostat that uses global velocity rescaling with Hamiltonian dynamics. 59 Regarding the force field, ClayFF, parameterized by Cygan et al 60 for describing the interaction of clays with water, is adopted to describe the quartz-water system in our simulations and the flexible single point charge (SPC) 61,62 water model is used in this forcefield.…”