Molecular dynamics simulations and instantaneous normal-mode analysis of the vibrational relaxation of the C-H stretching modes of N-methylacetamide-d in liquid deuterated water

Nonequilibrium molecular dynamics (MD) simulations and instantaneous normal mode (INMs) analyses are used to study the vibrational relaxation of the C-H stretching modes (ν(s)(CH?)) of deuterated N-methylacetamide (NMAD) in aqueous (D2O) solution. The INMs are identified unequivocally in terms of the equilibrium normal modes (ENMs), or groups of them, using a restricted version of the recently proposed Min-Cost assignment method. After excitation of the parent ν(s)(CH?) modes with one vibrational quantum, the vibrational energy is shown to dissipate through both intramolecular vibrational redistribution (IVR) and intermolecular vibrational energy transfer (VET). The decay of the vibrational energy of the ν(s)(CH?) modes is well fitted to a triple exponential function, with each characterizing a well-defined stage of the entire relaxation process. The first, and major, relaxation stage corresponds to a coherent ultrashort (τ(rel) = 0.07 ps) energy transfer from the parent ν(s)(CH?) modes to the methyl bending modes δ(CH?), so that the initially excited state rapidly evolves into a mixed stretch-bend state. In the second stage, characterized by a time of 0.92 ps, the vibrational energy flows through IVR to a number of mid-range-energy vibrations of the solute. In the third stage, the vibrational energy accumulated in the excited modes dissipates into the bath through an indirect VET process mediated by lower-energy modes, on a time scale of 10.6 ps. All the specific relaxation channels participating in the whole relaxation process are properly identified. The results from the simulations are finally compared with the recent experimental measurements of the ν(s)(CH?) vibrational energy relaxation in NMAD/D?O(l) reported by Dlott et al. (J. Phys. Chem. A 2009, 113, 75.) using ultrafast infrared-Raman spectroscopy.