An environmental pseudopotential approach to molecular interactions: Implementation in MOLPRO

We present the implementation into the MOLPRO package of a model for the interaction of a central system with its surrounding environment. The properties of a target system enclosed by a noncovalently bound environment or solvent are modeled as those of a system embedded into the effective pseudopotential arising from the exact electrostatic Coulomb potential and the approximated exchange–repulsion potential. For the latter we use the charge‐density overlap model, which relates the exchange–repulsion interaction energy between two species with the overlap of their ground‐state electron charge densities. The solutions of the modified Hartree–Fock equations for the target system are obtained self‐consistently. This way the exchange–induction effects arising from the converged electron–charge density of the embedded system are implicitly included. Inclusion of the correlation effects is provided by the use of post‐Hartree–Fock and density–functional techniques available in the MOLPRO package. The computational and conceptual advantages provided by this approach are shown in the calculation of the dipole polarizabilities of halide and chalcogenide anions in different environments. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 2075–2082, 2003

[1]  J. G. Snijders,et al.  Improved density functional theory results for frequency‐dependent polarizabilities, by the use of an exchange‐correlation potential with correct asymptotic behavior , 1996 .

[2]  R. Wheatley,et al.  A systematic model potential for Li+-H2O , 1995 .

[3]  P. Fowler,et al.  In-crystal ionic polarizabilities derived by combining experimental and ab initio results , 1985, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[4]  S. Price,et al.  An overlap model for estimating the anisotropy of repulsion , 1990 .

[5]  P. Fowler,et al.  Overlap model and ab initio cluster calculations of polarisabilities of ions in solids , 1999 .

[6]  N. Govind,et al.  Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems , 1999 .

[7]  Sarah L. Price,et al.  A systematic intermolecular potential method applied to chlorine , 1990 .

[8]  R. Wheatley,et al.  Intermolecular potential for the interaction of helium with ammonia , 2001 .

[9]  Robert Moszynski,et al.  Perturbation Theory Approach to Intermolecular Potential Energy Surfaces of van der Waals Complexes , 1994 .

[10]  P. Kollman.  Advances and Continuing Challenges in Achieving Realistic and Predictive Simulations of the Properties of Organic and Biological Molecules , 1996 .

[11]  J. Toennies,et al.  A reassessment of electron-gas density-functional calculations of short range repulsive potentials for rare gas atoms , 1994 .

[12]  Martin Karplus,et al.  Calculation of ground and excited state potential surfaces of conjugated molecules. I. Formulation and parametrization , 1972 .

[13]  S. Price,et al.  Use of molecular overlap to predict intermolecular repulsion in N ··· H—O hydrogen bonds , 1998 .

[14]  Keiji Morokuma,et al.  A TEST OF THE NEW INTEGRATED MO + MM (IMOMM) METHOD FOR THE CONFORMATIONAL ENERGY OF ETHANE AND N-BUTANE , 1996 .

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

[16]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[17]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[18]  Peter J. Knowles,et al.  Improved radial grids for quadrature in molecular density‐functional calculations , 1996 .

[19]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[20]  Jackson,et al.  Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. , 1992, Physical review. B, Condensed matter.

[21]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

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

[23]  Carole Van Caillie,et al.  Static and dynamic polarisabilities, Cauchy coefficients and their anisotropies: an evaluation of DFT functionals , 2000 .

[24]  C. W. Murray,et al.  Quadrature schemes for integrals of density functional theory , 1993 .

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

[26]  W. Kohn,et al.  An energy functional for surfaces , 2001 .

[27]  R. Wheatley,et al.  Scaled perturbation theory for weak intermolecular forces: the helium dimer , 2002 .

[28]  R. Gordon,et al.  Theory for the Forces between Closed‐Shell Atoms and Molecules , 1972 .

[29]  T. Wesołowski.  Density Functional Theory with approximate kinetic energy functionals applied to hydrogen bonds , 1997 .

[30]  M. Masamura.  Structures, Energetics, and Spectra of Cl-(H2O)nClusters,n= 1-6: Ab Initio Study , 2002 .

[31]  Elena Akhmatskaya,et al.  Monte Carlo simulation of F−(H2O)4 using an ab initio potential , 1999 .

[32]  Roger D. Amos,et al.  Static and dynamic polarisabilities, Cauchy coefficients and their anisotropies: a comparison of standard methods , 1998 .

[33]  S. Canuto.  Extreme electron correlation effects on the electric properties of atomic anions , 1994 .

[34]  K. Noda,et al.  Repulsive potentials for Cl−–R and Br−–R (R=He, Ne, and Ar) derived from beam experiments , 1976 .

[35]  T. Truong,et al.  Embedded density functional approach for calculations of adsorption on ionic crystals , 1996 .

[36]  Orlando Tapia,et al.  Self-consistent reaction field theory of solvent effects , 1975 .

[37]  Hans Peter Lüthi,et al.  Interaction energies of van der Waals and hydrogen bonded systems calculated using density functional theory: Assessing the PW91 model , 2001 .

[38]  Martin Head-Gordon,et al.  A method for two-electron Gaussian integral and integral derivative evaluation using recurrence relations , 1988 .

[39]  V. Barone,et al.  Accurate static polarizabilities by density functional theory: assessment of the PBE0 model , 1999 .

[40]  Nathaniel O. J. Malcolm,et al.  Cooperative effects in the structuring of fluoride water clusters: Ab initio hybrid quantum mechanical/molecular mechanical model incorporating polarizable fluctuating charge solvent , 1998 .

[41]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[42]  R. Wheatley,et al.  Application of the overlap model to calculating correlated exchange energies , 2000 .

[43]  Roland Lindh,et al.  The reduced multiplication scheme of the Rys quadrature and new recurrence relations for auxiliary function based two‐electron integral evaluation , 1991 .

[44]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[45]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[46]  P. Fowler,et al.  Overlap-model and ab initio cluster calculations of ion properties in distorted environments , 2001 .

[47]  Keiji Morokuma,et al.  Why do molecules interact? The origin of electron donor-acceptor complexes, hydrogen bonding and proton affinity , 1977 .

[48]  Application of density functional theory to calculation of in-crystal anionic polarizability , 1999 .

[49]  M. Vincent,et al.  Binding Energy of F(H2O)- and the Simulation of Fluoride Water Clusters Using a Hybrid QM/MM (Fluctuating Charge) Potential , 1999 .

[50]  Eduardo V. Ludeña,et al.  Energy Density Functional Theory of Many-Electron Systems , 1990 .

[51]  Richard J Wheatley,et al.  Intermolecular potential and second virial coefficient of the water-hydrogen complex. , 2002, The Journal of chemical physics.

[52]  A. Warshel,et al.  Frozen density functional approach for ab initio calculations of solvated molecules , 1993 .

[53]  Cortona,et al.  Self-consistently determined properties of solids without band-structure calculations. , 1991, Physical review. B, Condensed matter.

[54]  B. A. Hess,et al.  TDMP2 calculation of dynamic multipole polarizabilities and dispersion coefficients for the halogen anions F−, Cl−, Br− and I− , 1998 .

[55]  R. Wheatley.  An overlap model for exchange-induction: application to alkali halides , 1998 .

[56]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[57]  Allan H. Harvey,et al.  Intermolecular potentials and second virial coefficients of the water–neon and water–argon complexes , 2002 .

[58]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[59]  Y. Kim,et al.  Dependence of the closed-shell repulsive interaction on the overlap of the electron densities , 1981 .

[60]  R. Wheatley.  The solvation of sodium ions in water clusters: intermolecular potentials for Na+-H2O and H2O-H2O , 1996 .

[61]  Jongseob Kim,et al.  Structures, energetics, and spectra of fluoride–water clusters F−(H2O)n, n=1–6: Ab initio study , 1999 .

[62]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .