Investigating the Quartz ($$10\overline{10}$$)/Water Interface using Classical and Ab Initio Molecular Dynamics
- Vanderbilt Univ., Nashville, TN (United States). Dept. of Chemical Engineering
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Chemical Sciences Division
- Vanderbilt Univ., Nashville, TN (United States). Dept. of Chemical Engineering; Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Chemical Sciences Division
Two different terminations of the ($$10\overline{10}$$) surface of quartz (α and β) interacting with water are simulated by classical (CMD) (using two different force fields) and ab initio molecular dynamics (AIMD) and compared with previously published X-ray reflectivity (XR) experiments. Radial distribution functions between hydroxyl and water show good agreement between AIMD and CMD using the ClayFF force field for both terminations. The Lopes et al. (Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D. J. Phys. Chem. B 2006, 110, 27822792) force field (LFF), however, underestimates the extent of hydroxyl-water hydrogen bonding. The β termination is found to contain hydroxylhydroxyl hydrogen bonds; the quartz surface hydroxyl $$-$$ hydrogens and oxygens that hydrogen bond with each other exhibit greatly reduced hydrogen bonding to water. Conversely, the hydroxyl hydrogen and oxygens that are not hydrogen bonded to other surface hydroxyls but are connected to those that are show a considerable amount of hydrogen bonding to water. The electron density distribution of an annealed surface of quartz ($$10\overline{10}$$) obtained by XR is in qualitative agreement with electron densities calculated by CMD and AIMD. In all simulation methods, the interfacial water peak appears farther from the surface than observed by XR. Agreement among AIMD, LFF, and XR is observed for the relaxation of the near-surface atoms; however, ClayFF shows a larger discrepancy. Overall, results show that for both terminations of ($$10\overline{10}$$), LFF treats the near-surface structure more accurately whereas ClayFF treats the interfacial water structure more accurately. It is shown that the number of hydroxyl and water hydrogen bonds to the bridging Si$$-$$O$$-$$Si oxygens connecting the surface silica groups to the rest of the crystal is much greater for the α than the β termination. Lastly, it is suggested that this may play a role in the greater resistance to dissolution of the β termination than that of the α termination.
- Research Organization:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22); National Science Foundation (NSF)
- DOE Contract Number:
- AC05-00OR22725
- OSTI ID:
- 1025421
- Journal Information:
- Langmuir, Journal Name: Langmuir Journal Issue: 14 Vol. 27; ISSN 0743-7463
- Publisher:
- American Chemical Society
- Country of Publication:
- United States
- Language:
- English
Similar Records
Simulations of the quart (101-bar1)/water interface: A comparison of classical force fields, ab initi molecular dynamics, and x-ray reflectivity experiments.
Structures of quartz (10{ovr 1}0)- and (10{ovr 1}1)-water interfaces determined by x-ray reflectivity and atomic force microscopy of natural growth surfaces.
Journal Article
·
Fri Dec 31 23:00:00 EST 2010
· Journal of Physical Chemistry C
·
OSTI ID:1021988
Structures of quartz (10{ovr 1}0)- and (10{ovr 1}1)-water interfaces determined by x-ray reflectivity and atomic force microscopy of natural growth surfaces.
Journal Article
·
Sun Sep 01 00:00:00 EDT 2002
· Geochim. Cosmochim. Acta
·
OSTI ID:949484