Simulations of the Quartz(1011)/Water Interface: A Comparison of Classical Force Fields, Ab Initio Molecular Dynamics, and X-ray Reflectivity Experiments

Citations Over TimeTop 10% of 2011 papers

Abstract

Classical molecular dynamics (CMD) simulations of the (101̅1) surface of quartz interacting with bulk liquid water are performed using three different classical force fields, Lopes et al., ClayFF, and CHARMM water contact angle (CWCA), and compared to ab initio molecular dynamics (AIMD) and X-ray reflectivity (XR) results. The axial densities of the water and surface atoms normal to the surface are calculated and compared to previous XR experiments. Favorable agreement is shown for all the force fields with respect to the position of the water atoms. Analyses such as the radial distribution functions between water and hydroxyl atoms and the average cosine of the angle between the water dipole vector and the normal of the surface are also calculated for each force field. Significant differences are found between the different force fields from such analyses, indicating differing descriptions of the structured water in the near vicinity of the surface. AIMD simulations are also performed to obtain the water and hydroxyl structure for comparison among the predictions of the three classical force fields to better understand which force field is most accurate. It is shown that ClayFF exhibits the best agreement with the AIMD simulations for water hydroxyl radial distribution functions, suggesting that ClayFF treats the hydrogen bonding more accurately.

Related Papers

• → Polarizable Mean-Field Model of Water for Biological Simulations with AMBER and CHARMM Force Fields(2012)118 cited
• → Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D(2021)44 cited
• → Molecular Dynamics Simulations of Proteins: Can the Explicit Water Model Be Varied?(2007)71 cited
• → The quest for the best nonpolarizable water model from the adaptive force matching method(2010)59 cited
• → Derivation of MM2 force-field parameters for polychlorophosphazenes (NPCl2)n from ab initio data(1994)5 cited