Infrared Spectroscopy and Hydrogen-Bond Dynamics of Liquid Water from Centroid Molecular Dynamics with an Ab Initio-Based Force Field

Atomic forces for geometry-dependent point multipole and gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry-dependent multipole models. In this study, atomic force expressions for geometry-dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives partial derivative D(m'm)(l)/partial derivative Omega. The force equations can be applied to electrostatic models based on atomic point multipoles or gaussian multipole charge density. Hydrogen-bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry-dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry-dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry-dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry-dependent multipole models.

Infrared and Raman line shapes for ice Ih. I. Dilute HOD in H(2)O and D(2)O

Vibrational spectroscopy of ice Ih provides information about structure, dynamics, and vibrational coupling in this important substance. Vibrational spectra are simplified for HOD in either H(2)O or D(2)O, as in these instances the OD or OH stretch, respectively, functions as a local chromophore. As a first step in providing a theoretical treatment of the vibrational spectroscopy for the fully coupled system (H(2)O or D(2)O), herein we calculate the infrared and Raman spectra for the isotopically substituted systems. The calculation involves a classical molecular dynamics simulation using a new water model, an initial proton-disordered ice configuration, and ab initio based transition frequency, dipole, and polarizability maps. Our theoretical results are in reasonable agreement with experiment, and from our results we provide molecular and physical interpretations for the spectral features.