Premium
Atomic radii in molecules for use in a polarizable force field
Author(s) -
Swart Marcel,
Van Duijnen Piet TH.
Publication year - 2011
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/qua.22855
Subject(s) - chemistry , polarizability , intermolecular force , dipole , ab initio , atomic physics , molecule , coupled cluster , density functional theory , field (mathematics) , molecular physics , computational chemistry , physics , mathematics , organic chemistry , pure mathematics
We report here the results for an ab initio approach to obtain the parameters needed for molecular simulations using a polarizable force field. These parameters consist of the atomic charges, polarizabilities, and radii. The former two are readily obtained using methods reported previously (van Duijnen and Swart, J Phys Chem A 1998, 102 , 2399; Swart et al. J Comput Chem 2001, 22 , 79), whereas here we report a new approach for obtaining atomic second‐order radii (SOR), which is based on second‐order atomic moments in scaled Voronoi cells. These parameters are obtained from quantum‐chemistry calculations on the monomers, and used without further adaptation directly for intermolecular interactions. The approach works very well as shown here for four dimers, where high‐level coupled cluster with singles and doubles, and perturbative triples (CCSD(T)) and density functional theory (DFT) Swart‐Solà‐Bickelhaupt functional including Grimme's dispersion correction (SSB‐D) reference data are available for comparison. The energy surfaces for the three methods are very similar, which is also the case for the interaction between a water molecule with either a chloride anion or a sodium cation. These latter systems had previously been used to criticize Thole's damped point‐dipole method, but here we show that with the correct use of the method, it is perfectly able to describe the intermolecular interactions. This is most obvious for the induced dipole moment as function of the chloride–oxygen distance, where the direct (discrete) reaction field results are virtually indistinguishable from those obtained at CCSD(T)/aug‐ cc ‐pVTZ. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010