A smooth, nonsingular, and faithful discretization scheme for polarizable continuum models: The switching/Gaussian approach
- 28 December 2010
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 133 (24), 244111
- https://doi.org/10.1063/1.3511297
Abstract
Polarizable continuum models (PCMs) are a widely used family of implicit solvent models based on reaction-field theory and boundary-element discretization of the solute/continuum interface. An often overlooked aspect of these theories is that discretization of the interface typically does not afford a continuous potential energy surface for the solute. In addition, we show that discretization can lead to numerical singularities and violations of exact variational conditions. To fix these problems, we introduce the switching/Gaussian (SWIG) method, a discretization scheme that overcomes several longstanding problems with PCMs. Our approach generalizes a procedure introduced by York and Karplus [J. Phys. Chem. A 103, 11060 (1999)], extending it beyond the conductor-like screening model. Comparison to other purportedly smooth PCM implementations reveals certain artifacts in these alternative approaches, which are avoided using the SWIG methodology. The versatility of our approach is demonstrated via geometry optimizations, vibrational frequency calculations, and molecular dynamics simulations, for solutes described using quantum mechanics and molecular mechanics.Keywords
This publication has 50 references indexed in Scilit:
- Smooth Solvation Method for d-Orbital Semiempirical Calculations of Biological Reactions. 1. ImplementationThe Journal of Physical Chemistry B, 2005
- An integrated effective fragment—polarizable continuum approach to solvation: Theory and application to glycineThe Journal of Chemical Physics, 2002
- Electrostatics of nanosystems: Application to microtubules and the ribosomeProceedings of the National Academy of Sciences of the United States of America, 2001
- Comment on “Reaction field treatment of charge penetration” [J. Chem. Phys. 112, 5558 (2000)]The Journal of Chemical Physics, 2001
- Solvent Effects on Structure and Reaction Mechanism: A Theoretical Study of [2 + 2] Polar Cycloaddition between Ketene and ImineThe Journal of Physical Chemistry B, 1998
- A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectricsThe Journal of Chemical Physics, 1997
- Analytical derivatives for molecular solutes. II. Hartree–Fock energy first and second derivatives with respect to nuclear coordinatesThe Journal of Chemical Physics, 1994
- A simple method for the calculation of hydration enthalpies of polar molecules with arbitrary shapesThe Journal of Physical Chemistry, 1987
- Electric Moments of Molecules in LiquidsJournal of the American Chemical Society, 1936
- Volumen und Hydratationswärme der IonenThe European Physical Journal A, 1920