A reactive molecular dynamics simulation of the silica-water interface.

We report our study of a silica-water interface using reactive molecular dynamics. This first-of-its-kind simulation achieves length and time scales required to investigate the detailed chemistry of the system. Our molecular dynamics approach is based on the ReaxFF force field of van Duin et al. [J. Phys. Chem. A 107, 3803 (2003)]. The specific ReaxFF implementation (SERIALREAX) and force fields are first validated on structural properties of pure silica and water systems. Chemical reactions between reactive water and dangling bonds on a freshly cut silica surface are analyzed by studying changing chemical composition at the interface. In our simulations, reactions involving silanol groups reach chemical equilibrium in approximately 250 ps. It is observed that water molecules penetrate a silica film through a proton-transfer process we call "hydrogen hopping," which is similar to the Grotthuss mechanism. In this process, hydrogen atoms pass through the film by associating and dissociating with oxygen atoms within bulk silica, as opposed to diffusion of intact water molecules. The effective diffusion constant for this process, taken to be that of hydrogen atoms within silica, is calculated to be 1.68 x 10(-6) cm(2)/s. Polarization of water molecules in proximity of the silica surface is also observed. The subsequent alignment of dipoles leads to an electric potential difference of approximately 10.5 V between the silica slab and water.

[1]  Ananth Grama,et al.  Reactive Molecular Dynamics: Numerical Methods and Algorithmic Techniques , 2012, SIAM J. Sci. Comput..

[2]  F. Mauri,et al.  First-Principles Nuclear Magnetic Resonance Structural Analysis of Vitreous Silica , 2009 .

[3]  S. Sen,et al.  Intermediate range order in vitreous silica from a partial structure factor analysis , 2008 .

[4]  Mark S. Gordon,et al.  An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules. , 2008, The journal of physical chemistry. A.

[5]  M. Menziani,et al.  FFSiOH: a New Force Field for Silica Polymorphs and Their Hydroxylated Surfaces Based on Periodic B3LYP Calculations , 2008 .

[6]  S. Garofalini,et al.  Dissociative Chemisorption of Water onto Silica Surfaces and Formation of Hydronium Ions , 2008 .

[7]  A. V. van Duin,et al.  ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. , 2008, The journal of physical chemistry. A.

[8]  S. J. Singer,et al.  Model for the water-amorphous silica interface: the undissociated surface. , 2007, The journal of physical chemistry. B.

[9]  Susan B. Sinnott,et al.  Charge optimized many-body potential for the Si/SiO2 system , 2007 .

[10]  K. Schulten,et al.  Water-silica force field for simulating nanodevices. , 2006, The journal of physical chemistry. B.

[11]  A. Jansen,et al.  Mechanism of oligomerization reactions of silica. , 2006, The journal of physical chemistry. B.

[12]  Paolo Cignoni,et al.  Ambient Occlusion and Edge Cueing for Enhancing Real Time Molecular Visualization , 2006, IEEE Transactions on Visualization and Computer Graphics.

[13]  Jincheng Du,et al.  Molecular Dynamics Simulation of the Structure and Hydroxylation of Silica Glass Surfaces , 2005 .

[14]  Akira Takada,et al.  Molecular dynamics simulations of vitreous silica structures , 2004 .

[15]  M. Tomozawa,et al.  Interaction of water diffusion and hydrogen diffusion in silica glass at high temperatures , 2004 .

[16]  Kemin Wang,et al.  Bionanotechnology based on silica nanoparticles , 2004, Medicinal research reviews.

[17]  Mao-Hua Du,et al.  Hydrolysis of a two-membered silica ring on the amorphous silica surface. , 2004, The Journal of chemical physics.

[18]  Mao-Hua Du,et al.  Water–silica surface interactions: A combined quantum-classical molecular dynamic study of energetics and reaction pathways , 2003 .

[19]  A. V. Duin,et al.  ReaxFFSiO Reactive Force Field for Silicon and Silicon Oxide Systems , 2003 .

[20]  N. Kuzuu,et al.  Molecular dynamics study on surface structure of a-SiO2 by charge equilibration method , 2003 .

[21]  S T Pantelides,et al.  Reactions and diffusion of water and oxygen molecules in amorphous SiO2. , 2002, Physical review letters.

[22]  K. Binder,et al.  Classical and ab-initio molecular dynamic simulation of an amorphous silica surface , 2001, cond-mat/0108481.

[23]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[24]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

[25]  W. Tan,et al.  Biochemically functionalized silica nanoparticles. , 2001, The Analyst.

[26]  L. T. Zhuravlev.  The surface chemistry of amorphous silica. Zhuravlev model , 2000 .

[27]  Dean R. Haeffner,et al.  Electron distribution in water , 2000 .

[28]  E. Papirer.  Adsorption on Silica Surfaces , 2000 .

[29]  N. T. Huff,et al.  Factors affecting molecular dynamics simulated vitreous silica structures , 1999 .

[30]  D. Tildesley,et al.  The free surface of water: molecular orientation, surface potential and nonlinear susceptibility , 1997 .

[31]  Michiel Sprik,et al.  Ab initio molecular dynamics simulation of liquid water: Comparison of three gradient‐corrected density functionals , 1996 .

[32]  Binder,et al.  Cooling-rate effects in amorphous silica: A computer-simulation study. , 1996, Physical review. B, Condensed matter.

[33]  Minoru Tomozawa,et al.  An infrared spectroscopic study of water-related species in silica glasses , 1996 .

[34]  T. Lubensky,et al.  Principles of condensed matter physics , 1995 .

[35]  A. Soper.  Orientational correlation function for molecular liquids: The case of liquid water , 1994 .

[36]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[37]  K. Eisenthal,et al.  Polarization of water molecules at a charged interface: second harmonic studies of the silica/water interface , 1992 .

[38]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[39]  Rino,et al.  Interaction potential for SiO2: A molecular-dynamics study of structural correlations. , 1990, Physical review. B, Condensed matter.

[40]  A. Wright,et al.  Neutron scattering from vitreous silica IV. Time-of-flight diffraction☆ , 1990 .

[41]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[42]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[43]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[44]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[45]  W. J. Orville-Thomas,et al.  Studies of intermolecular interactions by matrix isolation vibrational spectroscopy: Self-association of water , 1980 .

[46]  Kazimierz Krynicki,et al.  Pressure and temperature dependence of self-diffusion in water , 1978 .

[47]  D. G. Pinatti,et al.  A refinement of the structure of vitreous silica , 1975 .

[48]  V. Gold,et al.  Proton-Transfer Reactions , 1975 .

[49]  Laurence S. Rothman,et al.  Dipole moment of water from Stark measurements of H2O, HDO, and D2O , 1973 .

[50]  R. Mozzi,et al.  The structure of vitreous silica , 1969 .

[51]  K. Kuchitsu,et al.  Estimation of Anharmonic Potential Constants. I. Linear XY2 Molecules , 1965 .

[52]  R. Mcdonald.  Surface Functionality of Amorphous Silica by Infrared Spectroscopy , 1958 .

[53]  W. S. Benedict,et al.  Rotation‐Vibration Spectra of Deuterated Water Vapor , 1956 .