Implementation of a Morse potential to model hydroxyl behavior in phyllosilicates.

The accurate molecular simulation of many hydrated chemical systems, including clay minerals and other phyllosilicates and their interfaces with aqueous solutions, requires improved classical force field potentials to better describe structure and vibrational behavior. Classical and ab initio molecular dynamics simulations of the bulk structure of pyrophyllite, talc, and Na-montmorillonite clay phases exhibit dissimilar behavior in the hydroxyl stretch region of power spectra derived from atomic trajectories. The classical simulations, using the CLAYFF force field, include either a standard harmonic potential or a new Morse potential parametrized for both dioctahedral and trioctahedral phases for the O-H bond stretch. Comparisons of classical results with experimental values and with ab initio molecular dynamics simulations indicate improvements in the simulation of hydroxyl orientation relative to the clay octahedral sheet and in the O-H bond stretch in the high frequency region of the power spectrum.

[1]  Thomas A. Halgren,et al.  The representation of van der Waals (vdW) interactions in molecular mechanics force fields: potential form, combination rules, and vdW parameters , 1992 .

[2]  L. Rintoul,et al.  Lattice vibrations of montmorillonite: an FT Raman and X-ray diffraction study , 1996 .

[3]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[4]  I. Swainson,et al.  Triclinic muscovite; X-ray diffraction, neutron diffraction and photo-acoustic FTIR spectroscopy , 1998 .

[5]  V. Farmer.  The Infrared spectra of minerals , 1974 .

[6]  William Jones,et al.  On the application of computer simulation techniques to anionic and cationic clays: A materials chemistry perspective , 2006 .

[7]  V. Farmer.  The Layer Silicates , 1974 .

[8]  H. Lischka,et al.  Formation of 2,4–D complexes on montmorillonites – an ab initio molecular dynamics study , 2007 .

[9]  C. I. Sainz-Díaz,et al.  Isomorphous cation substitution in dioctahedral phyllosilicates by means of ab initio quantum mechanical calculations on clusters , 2003 .

[10]  Garrison Sposito,et al.  The surface chemistry of soils , 1984 .

[11]  S. Churakov.  Ab initio study of sorption on pyrophyllite: structure and acidity of the edge sites. , 2006, The journal of physical chemistry. B.

[12]  Randall T. Cygan,et al.  Molecular Models of Hydroxide, Oxyhydroxide, and Clay Phases and the Development of a General Force Field , 2004 .

[13]  S. Churakov.  Structure and dynamics of the water films confined between edges of pyrophyllite : A first principle study , 2007 .

[14]  H. Heinz,et al.  Molecular models and simulations of layered materials , 2009 .

[15]  H. Lischka,et al.  Ab initio calculations of relative stabilities of different structural arrangements in dioctahedral phyllosilicates , 2007 .

[16]  C. Johnston,et al.  Infrared Study of Water Sorption on Na-, Li-, Ca-, and Mg-Exchanged (SWy-1 and SAz-1) Montmorillonite , 2000 .

[17]  C. I. Sainz-Díaz,et al.  Hydrogen bonding and vibrational properties of hydroxy groups in the crystal lattice of dioctahedral clay minerals by means of first principles calculations , 2004 .

[18]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[19]  James P. Larentzos,et al.  An ab Initio and Classical Molecular Dynamics Investigation of the Structural and Vibrational Properties of Talc and Pyrophyllite , 2007 .

[20]  G. S. Premachandra,et al.  Polarized ATR-FTIR Study of Smectite in Aqueous Suspension , 2001 .

[21]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[22]  H. Lischka,et al.  SORPTION OF PHENOXYACETIC ACID HERBICIDES ON THE KAOLINITE MINERAL SURFACE – AN AB INITIO MOLECULAR DYNAMICS SIMULATION , 2004 .

[23]  J. Reimers,et al.  A local mode potential function for the water molecule , 1984 .

[24]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[25]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[26]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[27]  G. Brindley,et al.  Crystal Structures of Clay Minerals and their X-ray Identification , 1982 .

[28]  Free energy of liquid water on the basis of quasichemical theory and ab initio molecular dynamics. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.