TCB Publications - Abstract

Francois Dehez, Janos G. Angyan, Ignacio Soteras Gutierrez, F. Javier Luque, Klaus Schulten, and Christophe Chipot. Modeling induction phenomena in intermolecular interactions with an ab initio force field. Journal of Chemical Theory and Computation, 3:1914-1926, 2007.

DEHE2007 One possible road towards the development of a polarizable potential energy function relies on the use of distributed polarizabilities derived from the induction energy mapped around the molecule. Whereas such polarizable models are expected to reproduce the signature induction energy with an appreciable accuracy, it is far from clear whether they will perform equally well in the context of intermolecular interactions. To address this issue, while pursuing the ultimate goal of a "plug-and-play"-like approach, polarizability models determined quantum mechanically and consisting of atomic isotropic dipole plus charge-flow polarizabilities were combined with the classical, non-polarizable Charmm force field. Performance of the models was probed in the challenging test cases of cation-$\pi$ binding and the association of a divalent calcium ion with water, where induction effects are envisioned to be considerable. Since brute force comparison of the binding energies estimated from the polarizable and the classical Charmm potential energy functions is not justified, the individual electrostatic and induction contributions of the force field were confronted to the corresponding terms of a symmetry-adapted perturbation theory (SAPT) expansion conducted at the MP2/6-311++G(d,p) level of approximation. While the quantum-mechanical and the molecular-mechanical electrostatic and damped induction contributions agree reasonably well, overall reproduction of the binding energies is plagued by an underestimated repulsion that underlines the necessity of de novo parametrization of the classical 6-12 form of the van der Waals potential. Based on the SAPT expansion, new Lennard-Jones parameters were optimized, which, combined with the remainder of the polarizable force field, yield an improved reproduction of the target binding energies.


Request Full Text

Request Paper

Full Name
Email Address
Institution
Type the number eight in the box