Polarizable intermolecular potentials for water and benzene interacting with halide and metal ions

A complete derivation of polarizable intermolecular potentials based on high-level, gas-phase quantum-mechanical calculations is proposed. The importance of appreciable accuracy together with inherent simplicity represents a significant endeavor when enhancement of existing force fields for biological systems is sought. Toward this end, symmetry-adapted perturbation theory can provide an expansion of the total interaction energy into physically meaningful, for example, electrostatic, induction and van der Waals terms. Each contribution can be readily compared with its counterpart in classical force fields. Since the complexity of the different intermolecular terms cannot be fully embraced using a minimalist description, it is necessary to resort to polyvalent expressions capable of encapsulating overlooked contributions from the quantum-mechanical expansion. This choice results in consistent force field components that reflect the underlying physical principles of the phenomena. This simplified potential energy function is detailed, and definitive guidelines are drawn. As a proof of concept, the methodology is illustrated through a series of test cases that include the interaction of water and benzene with halide and metal ions. In each case considered, the total energy is reproduced accurately over a range of biologically relevant distances.