Quantum vibrational state-dependent potentials for classical many-body simulations

In this paper, we present a method for constructing simple state-dependent many-body potentials for quantum vibrations in a classical bath. The approach is based on an adiabatic separation between high-frequency quantum vibrational modes of the solute and the lower frequency classical motion of the solvent, and on a first-order perturbation theory description of the dependence of the quantum energies on bath configuration. In the simplest realization of the method, the delocalized quantum probability density of the vibrational mode is approximated by a sum of two delta functions, with positions and weights chosen to represent the lowest three moments of the exact distribution. Thus, in the many-body description of the system, each atom describing the quantum vibration is represented by a pair of particles. These quantum particles are held in rigid relative position and interact with the bath via potentials the magnitudes of which are modified by the delta-function weights. The resulting approach allows the classical molecular dynamics of molecules in arbitrary quantum vibrational states to be simulated with a little more effort than a purely classical description. The applicability of the method is illustrated in many-body simulations of the dephasing of vibrational superposition states of I(2) in a cryogenic krypton matrix, yielding results in good agreement with experiment.