Premium
Using electronic polarization from the internal continuum (EPIC) for intermolecular interactions
Author(s) -
Truchon JeanFrançois,
Nicholl's Anthony,
Grant J. Andrew,
Iftimie Radu I.,
Roux Benoît,
Bayly Christopher I.
Publication year - 2010
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.21369
Subject(s) - polarizability , intermolecular force , ab initio , chemistry , intramolecular force , ab initio quantum chemistry methods , dielectric , polarization (electrochemistry) , water dimer , atom (system on chip) , atomic physics , molecular physics , computational chemistry , molecule , physics , quantum mechanics , hydrogen bond , organic chemistry , stereochemistry , computer science , embedded system
Recently, the vacuum‐phase molecular polarizability tensor of various molecules has been accurately modeled (Truchon et al., J Chem Theory Comput 2008, 4, 1480) with an intramolecular continuum dielectric model. This preliminary study showed that electronic polarization can be accurately modeled when combined with appropriate dielectric constants and atomic radii. In this article, using the parameters developed to reproduce ab initio quantum mechanical (QM) molecular polarizability tensors, we extend the application of the “electronic polarization from internal continuu” (EPIC) approach to intermolecular interactions. We first derive a dielectric‐adapted least‐square‐fit procedure similar to RESP, called DRESP, to generate atomic partial charges based on a fit to a QM abinitio electrostatic potential (ESP). We also outline a procedure to adapt any existing charge model to EPIC. The ability of this to reproduce local polarization, as opposed to uniform polarization, is also examined leading to an induced ESP relative root mean square deviation of 1%, relative to ab initio , when averaged over 37 molecules including aromatics and alkanes. The advantage of using a continuum model as opposed to an atom‐centered polarizable potential is illustrated with a symmetrically perturbed atom and benzene. We apply EPIC to a cation‐π binding system formed by an atomic cation and benzene and show that the EPIC approach can accurately account for the induction energy. Finally, this article shows that the ab initio electrostatic component in the difficult case of the H‐bonded 4‐pyridone dimer, a highly polar and polarized interaction, is well reproduced without adjusting the vacuum‐phase parameters. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010