Abstract
The proper description of the electric environment of condensed phases is a critical challenge for force field methods. To test and validate the ability of the CHARMM additive force field to describe the electric environment in aqueous solution combined QM/MM calculations have been used to calculate the vibrational Stark effect (VSE). We utilized a first principles methodology using correlated electronic structure techniques to compute the Stark shift between the gas phase and solvent environments and between two different solvent environments of three VSE probes containing acetonitrile or fluorine functionalities which have been well-characterized experimentally. Reasonable agreement with the experimentally determined Stark shifts is obtained when the MM atoms are described by the CHARMM additive force field, though it is essential to employ an anharmonic correction in the frequency calculation. In addition, the electric field created by the solvent is computed along the CN bond and a theoretical Stark tuning rate is determined for acetonitrile and shown to be in satisfactory agreement with experiment.