Using electronic polarization from the internal continuum (EPIC) for intermolecular interactions.

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-pi 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.

Accurate Molecular Polarizabilities Based on Continuum Electrostatics.

A novel approach for representing the intramolecular polarizability as a continuum dielectric is introduced to account for molecular electronic polarization. It is shown, using a finite-difference solution to the Poisson equation, that the Electronic Polarization from Internal Continuum (EPIC) model yields accurate gas-phase molecular polarizability tensors for a test set of 98 challenging molecules composed of heteroaromatics, alkanes and diatomics. The electronic polarization originates from a high intramolecular dielectric that produces polarizabilities consistent with B3LYP/aug-cc-pVTZ and experimental values when surrounded by vacuum dielectric. In contrast to other approaches to model electronic polarization, this simple model avoids the polarizability catastrophe and accurately calculates molecular anisotropy with the use of very few fitted parameters and without resorting to auxiliary sites or anisotropic atomic centers. On average, the unsigned error in the average polarizability and anisotropy compared to B3LYP are 2% and 5%, respectively. The correlation between the polarizability components from B3LYP and this approach lead to a R2 of 0.990 and a slope of 0.999. Even the F2 anisotropy, shown to be a difficult case for existing polarizability models, can be reproduced within 2% error. In addition to providing new parameters for a rapid method directly applicable to the calculation of polarizabilities, this work extends the widely used Poisson equation to areas where accurate molecular polarizabilities matter.