Dynamic Charge Equilibration‐Morse stretch force field: Application to energetics of pure silica zeolites

We present the Dynamic Charge Equilibration (DQEq) method for a self‐consistent treatment of charge transfer in force field modeling, where atomic charges are designed to reproduce electrostatic potentials calculated quantum mechanically. Force fields coupled with DQEq allow charges to readjust as geometry changes in classical simulations, using appropriate algorithms for periodic boundary conditions. The full electrostatic energy functional is used to derive the corresponding forces and the second derivatives (hessian) for vibrational calculations. Using DQEq electrostatics, we develop a simple nonbond force field for simulation of silica molecular sieves, where nonelectrostatic interactions are described by two‐body Morse stretch terms. Energy minimization calculations with the new force field yield accurate unit cell geometries for siliceous zeolites. Relative enthalpies with respect to quartz and third‐law entropies calculated from harmonic vibrational analysis agree very well with available calorimetric data: calculated SiO2 enthalpies relative to α‐quartz are within 2 kJ/mol and entropies at 298 K are within 3 J/mol K of the respective experimental values. Contributions from the zero point energy and vibrational degrees of freedom were found to be only about 1 kJ/mol for the free energy of mutual transformations between microporous silica polymorphs. The approach presented here can be applied to interfaces and other oxides as well and it is suitable for development of force fields for accurate modeling of geometry and energetics of microporous and mesoporous materials, while providing a realistic description of electrostatic fields near surfaces and inside pores of adsorbents and catalysts. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1507–1514, 2002

[1]  Mark E. Davis,et al.  Entropy of Pure-Silica Molecular Sieves , 2001 .

[2]  W. Goddard,et al.  The MS-Q Force Field for Clay Minerals: Application to Oil Production , 2001 .

[3]  J. Gale,et al.  Atomistic simulation of the crystal structures and bulk moduli of TiO , 2001 .

[4]  Mark E. Davis,et al.  Thermochemistry of Pure-Silica Zeolites , 2000 .

[5]  W. Goddard,et al.  Morse Stretch Potential Charge Equilibrium Force Field for Ceramics: Application to the Quartz-Stishovite Phase Transition and to Silica Glass , 1999 .

[6]  E. Vogt,et al.  Semi-empirical atomic charges for use in computational chemistry of molecular sieves , 1998 .

[7]  Jianfen Fan,et al.  Extending and simplifying the electronegativity equalization method , 1998 .

[8]  Christopher W. Jones,et al.  Synthesis, Characterization, and Structure Solution of CIT-5, a New, High-Silica, Extra-Large-Pore Molecular Sieve , 1998 .

[9]  K. Knight,et al.  Calibration of excess thermodynamic properties and elastic constant variations associated with the alpha beta phase transition in quartz , 1998 .

[10]  P. Lightfoot,et al.  Synthesis and Structure of ITQ‐3, the First Pure Silica Polymorph with a Two‐Dimensional System of Straight Eight‐Ring Channels , 1997 .

[11]  A. Rappé,et al.  Toward an Understanding of Zeolite Y as a Cracking Catalyst with the Use of Periodic Charge Equilibration , 1996 .

[12]  T. Ohsuna,et al.  Direct Observation of “Pure MEL Type” Zeolite , 1996 .

[13]  Jansen,et al.  Structure-stability relationships for all-silica structures. , 1995, Physical review. B, Condensed matter.

[14]  A. Cheetham,et al.  A Synchrotron X-ray Diffraction, Neutron Diffraction, 29Si MAS-NMR, and Computational Study of the Siliceous Form of Zeolite Ferrierite , 1994 .

[15]  B. Berne,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994, chem-ph/9406002.

[16]  R. Downs,et al.  The pressure behavior of alpha cristobalite , 1994 .

[17]  Mark E. Davis,et al.  Thermochemical study of the stability of frameworks in high silica zeolites , 1993 .

[18]  H. Gies,et al.  The structure of zeolite ZSM-23 (MTT) refined from synchrotron X-ray powder data , 1993 .

[19]  A. Cheetham,et al.  Powder Neutron Diffraction and 29Si MAS NMR Studies of Siliceous Zeolite-Y , 1993 .

[20]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[21]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

[22]  J. Parise,et al.  Elasticity of α-Cristobalite: A Silicon Dioxide with a Negative Poisson's Ratio , 1992, Science.

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

[24]  Jan K. Labanowski,et al.  Density Functional Methods in Chemistry , 1991 .

[25]  H. Gies,et al.  Crystal structure of silica-ZSM-12 by the combined use of high-resolution solid-state MAS NMR spectroscopy and synchrotron x-ray powder diffraction , 1990 .

[26]  J. C. Jansen,et al.  The monoclinic framework structure of zeolite H-ZSM-5. Comparison with the orthorhombic framework of as-synthesized ZSM-5 , 1990 .

[27]  S. L. Lawton,et al.  The Framework Topology of ZSM-18, a Novel Zeolite Containing Rings of Three (Si,Al)-O Species , 1990, Science.

[28]  N. Karasawa,et al.  Acceleration of convergence for lattice sums , 1989 .

[29]  Aoki,et al.  First-principles interatomic potential of silica applied to molecular dynamics. , 1988, Physical review letters.

[30]  C. Catlow,et al.  Computer Simulation Studies of Zeolite Structure , 1988 .

[31]  J. C. Jansen,et al.  On the location and disorder of the tetrapropylammonium (TPA) ion in zeolite ZSM-5 with improved framework accuracy , 1987 .

[32]  G. Artioli,et al.  Crystal structure of coesite, a high-pressure form of SiO/sub 2/, at 15 and 298 K from single-crystal neutron and x-ray diffraction data: test of bonding models , 1987 .

[33]  J. Faber,et al.  Crystal structure of low cristobalite at 10, 293, and 473 K: Variation of framework geometry with temperature , 1985 .

[34]  A. Navrotsky,et al.  The quartz-coesite-stishovite transformations: new calorimetric measurements and calculation of phase diagrams , 1984 .

[35]  P. Richet,et al.  Thermodynamic properties of quartz, cristobalite and amorphous SiO2: drop calorimetry measurements between 1000 and 1800 K and a review from 0 to 2000 K , 1982 .

[36]  J. D. Jorgensen,et al.  Crystal structure and thermal expansion of α‐quartz SiO2 at low temperatures , 1982 .

[37]  R. N. Thurston,et al.  Elastic Moduli of Quartz versus Hydrostatic Pressure at 25 and-195.8C , 1965 .

[38]  C. C. J. Roothaan,et al.  A Study of Two‐Center Integrals Useful in Calculations on Molecular Structure. I , 1951 .

[39]  Nathan Rosen,et al.  Calculation of Interaction between Atoms with s-Electrons , 1931 .

[40]  J. Mintmire,et al.  Efficient parallel algorithms for molecular dynamics simulations using variable charge transfer electrostatic potentials , 2000 .

[41]  A. Corma,et al.  The role of the electrostatic potential, electric field and electric field gradient on the acidity of AFI and CHA zeotypes , 2000 .

[42]  C. Baerlocher,et al.  Synthesis and Structural Characterization of MWW Type Zeolite ITQ-1, the Pure Silica Analog of MCM-22 and SSZ-25 , 1998 .

[43]  Konstantin S. Smirnov,et al.  Consistent implementation of the electronegativity equalization method in molecular mechanics and molecular dynamics , 1996 .

[44]  I. Ohno.  Temperature Variation of Elastic Properties of α- Quartz up to the α-β Transition , 1995 .

[45]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

[46]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[47]  R. T. Sanderson.  Chemical Bonds and Bond Energy , 1976 .