AIPAR: ab initio parametrization of intermolecular potentials for computer simulations

An unambiguous, fully ab initio and automated technique denoted AIPAR (“ab initio parametrization”) implemented in the SJBR program has been proposed to yield intermolecular interaction potentials between polar molecules and water. The AIPAR procedure has been applied to several organic molecules covering a wide range of structure and functional groups, namely methanol, acetone (propanone), methanethiol (methyl mercaptan), imidazole (1,3-diazole), oxazole and furan. The AIPAR-derived sets of parameters compare well with the empirical OPLS ones, mainly when the all-atoms model is employed in the OPLS procedure. Monte Carlo simulations were performed for an aqueous solution of methanol and for an equimolar binary mixture methanol–water using the AIPAR and OPLS parameters. The thermodynamic and geometric results obtained with the parameters obtained with the AIPAR procedure compare favorably with the OPLS simulations, even for the binary mixture, demonstrating the precision, robustness and transferability of the parameters obtained with the AIPAR procedure.Figure Superimposed configurations of water (without the hydrogen atoms) around the methanol molecule obtained with the AGOA procedure.

[1]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[2]  G. Karlstroem,et al.  New intermolecular energy calculation scheme: applications to potential surface and liquid properties of water , 1990 .

[3]  Antonio Rizzo,et al.  Selected topics in ab initio computational chemistry in both very small and very large chemical systems , 1991 .

[4]  William L. Jorgensen,et al.  Ab initio and Monte Carlo calculations for a nucleophilic addition reaction in the gas phase and in aqueous solution , 1986 .

[5]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid alcohols , 1986 .

[6]  William L. Jorgensen,et al.  Temperature and size dependence for Monte Carlo simulations of TIP4P water , 1985 .

[7]  E. Clementi,et al.  Analytical potentials from "ab initio" computations for the interaction between biomolecules. 2. Water with the bases of DNA. , 1977, Journal of the American Chemical Society.

[8]  Harold A. Scheraga,et al.  Analytical intermolecular potential functions from ab initio SCF calculations for hydration of methylamine and methylammonium ion , 1990 .

[9]  William L. Jorgensen,et al.  Intermolecular potential functions and Monte Carlo simulations for liquid sulfur compounds , 1986 .

[10]  Scott F. Smith,et al.  Theoretical examination of the SN2 reaction involving chloride ion and methyl chloride in the gas phase and aqueous solution , 1985 .

[11]  E. Clementi,et al.  Analytical potentials from "ab initio" computations for the interaction between biomolecules. 3. Reliability and transferability of the pair potentials. , 1977, Journal of the American Chemical Society.

[12]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[13]  C. Reichardt.  Solvents and Solvent Effects in Organic Chemistry , 1988 .

[14]  E. Clementi,et al.  Analytical potentials from ab initio computations for the interaction between biomolecules. VII. Polar amino acids and conclusions , 1981 .

[15]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[16]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[17]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[18]  E. Clementi,et al.  Analytical potentials from "ab initio" computations for the interaction between biomolecules. 1. Water with amino acids. , 1977, Journal of the American Chemical Society.

[19]  Ranbir Singh,et al.  J. Mol. Struct. (Theochem) , 1996 .

[20]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[21]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[22]  S. Canuto,et al.  Solvent Effects from a Sequential Monte Carlo - Quantum Mechanical Approach , 1997 .

[23]  F. Sato,et al.  On the Transferability of Force Field ParametersWith an ab Initio Force Field Developed for Sulfonamides , 2003 .

[24]  E. Clementi,et al.  Analytical potentials from ab initio computations for the interaction between biomolecules. V. Formyl‐triglycyl amide and water , 1979 .

[25]  Monte Carlo simulation of the binary liquid mixture water—methanol , 1993 .

[26]  William L. Jorgensen,et al.  Development of an All-Atom Force Field for Heterocycles. Properties of Liquid Pyrrole, Furan, Diazoles, and Oxazoles , 1998 .

[27]  U. Lüning.  Chr. Reichardt: Solvents and Solvent Effects in Organic Chemistry, second, completely revised and enlarged edition, VCH Verlagsgesellschaft, Weinheim, Basel, Cambridge, New York 1988. 534 Seiten, Preis: DM 148,–. , 1989 .

[28]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[29]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[30]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

[31]  W. L. Jorgensen.  OPLS Force Fields , 2002 .

[32]  Claudio N. Cavasotto,et al.  Direct Derivation of van der Waals Force Field Parameters from Quantum Mechanical Interaction Energies , 2003 .