Inclusion of the quadrupole moment when describing polarization. The effect of the dipole‐quadrupole polarizability

A method to compute distributed dipole–quadrupole polarizabilities is suggested. The method is based on numerical differentiation of distributed quadrupole moments, using finite field perturbation calculations. It is tested using two different multicenter multipole expansions, and compared with results using polarizabilities obtained via the uncoupled Hartree–Fock approximation. The accuracy of these dipole–quadrupole polarizabilities are tested for different molecules and basis sets, by comparing the induced electrostatic potential of the Hartree–Fock density with the induced electrostatic potential of the polarization models. This is done by perturbing the molecules with an external homogeneous field and with an external dipole. It is found that inclusion of the dipole–quadrupole polarizability significantly improves the accuracy of the response of the molecule to these external perturbations. This suggests that inclusion of higher‐order induced moments can be of importance when improving the description of intermolecular interactions using force fields. © 2008 Wiley Periodicals, Inc. J Comput Chem 2008

[1]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[2]  Herman J. C. Berendsen,et al.  A mobile charge densities in harmonic oscillators (MCDHO) molecular model for numerical simulations: The water-water interaction , 2000 .

[3]  D. Hagberg,et al.  On the solvation of ions in small water droplets. , 2005, The journal of physical chemistry. B.

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

[5]  M. V. Subbotin,et al.  A quantum mechanical polarizable force field for biomolecular interactions , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Samuel Krimm,et al.  A new electrostatic model for molecular mechanics force fields , 2000 .

[7]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[8]  K. Mikkelsen,et al.  An atomic capacitance–polarizability model for the calculation of molecular dipole moments and polarizabilities , 2001 .

[9]  Pär Söderhjelm,et al.  Accuracy of typical approximations in classical models of intermolecular polarization. , 2008, The Journal of chemical physics.

[10]  Roland Lindh,et al.  Local properties of quantum chemical systems: the LoProp approach. , 2004, The Journal of chemical physics.

[11]  A. Buckingham.  Permanent and Induced Molecular Moments and Long‐Range Intermolecular Forces , 2007 .

[12]  Patrick Bultinck,et al.  Fast Calculation of Quantum Chemical Molecular Descriptors from the Electronegativity Equalization Method , 2003, J. Chem. Inf. Comput. Sci..

[13]  L. Renee Olano,et al.  Fluctuating charge normal modes: An algorithm for implementing molecular dynamics simulations with polarizable potentials , 2005, J. Comput. Chem..

[14]  G. Karlström.  Local polarizabilities in molecules, based on ab initio Hartree-Fock calculations , 1982 .

[15]  O Engkvist,et al.  Accurate Intermolecular Potentials Obtained from Molecular Wave Functions: Bridging the Gap between Quantum Chemistry and Molecular Simulations. , 2000, Chemical reviews.

[16]  K. D. Collins,et al.  The Hofmeister effect and the behaviour of water at interfaces , 1985, Quarterly Reviews of Biophysics.

[17]  Zhong-Zhi Yang,et al.  Atom-Bond Electronegativity Equalization Method Fused into Molecular Mechanics. II. A Seven-Site Fluctuating Charge and Flexible Body Water Potential Function for Liquid Water , 2004 .

[18]  B. Roos,et al.  Molcas: a program package for computational chemistry. , 2003 .

[19]  Zhong-Zhi Yang,et al.  Atom-bond electronegativity equalization method fused into molecular mechanics. I. A seven-site fluctuating charge and flexible body water potential function for water clusters. , 2004, The Journal of chemical physics.

[20]  Per-Olof Widmark,et al.  Density matrix averaged atomic natural orbital (ANO) basis sets for correlated molecular wave functions , 1990 .

[21]  Arieh Warshel,et al.  Polarizable Force Fields:  History, Test Cases, and Prospects. , 2007, Journal of chemical theory and computation.

[22]  K. Sundberg,et al.  An atom monopole–dipole interaction model with charge transfer for the treatment of polarizabilities of π‐bonded molecules. , 1978 .

[23]  Henry Margenau,et al.  Theory of intermolecular forces , 1969 .

[24]  J. D. Bernal,et al.  A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions , 1933 .

[25]  Roland Lindh,et al.  Accuracy of distributed multipoles and polarizabilities: Comparison between the LoProp and MpProp models , 2007, J. Comput. Chem..

[26]  Zhong-Zhi Yang,et al.  General atom-bond electronegativity equalization method and its application in prediction of charge distributions in polypeptide , 2000 .

[27]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[28]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[29]  F. Stillinger,et al.  Improved simulation of liquid water by molecular dynamics , 1974 .

[30]  Douglas J. Tobias,et al.  Ions at the Air/Water Interface , 2002 .