Simultaneous spectral and temporal analyses of kinetic energies in nonequilibrium systems: theory and application to vibrational relaxation of O-D stretch mode of HOD in water

J Phys Chem A. 2015 May 28;119(21):5356-67. doi: 10.1021/jp510157y. Epub 2014 Dec 22.

Abstract

A time series of kinetic energies (KE) from classical molecular dynamics (MD) simulation contains fundamental information on system dynamics. It can also be analyzed in the frequency domain through Fourier transformation (FT) of velocity correlation functions, providing energy content of different spectral regions. By limiting the FT time span, we have previously shown that spectral resolution of KE evolution is possible in the nonequilibrium situations [Jeon and Cho, J. Chem. Phys. 2011, 135, 214504]. In this paper, we refine the method by employing the concept of instantaneous power spectra, extending it to reflect an instantaneous time-correlation of velocities with those in the future as well as with those in the past, and present a new method to obtain the instantaneous spectral density of KE (iKESD). This approach enables the simultaneous spectral and temporal resolution of KE with unlimited time precision. We discuss the formal and novel properties of the new iKESD approaches and how to optimize computational methods and determine parameters for practical applications. The method is specifically applied to the nonequilibrium MD simulation of vibrational relaxation of the OD stretch mode in a hydrated HOD molecule by employing a hybrid quantum mechanical/molecular mechanical (QM/MM) potential. We directly compare the computational results with the OD band population relaxation time profiles extracted from the IR pump-probe measurements for 5% HOD in water. The calculated iKESD yields the OD bond relaxation time scale ∼30% larger than the experimental value, and this decay is largely frequency-independent if the classical anharmonicity is accounted for. From the integrated iKESD over intra- and intermolecular bands, the major energy transfer pathways were found to involve the HOD bending mode in the subps range, then the internal modes of the solvent until 5 ps after excitation, and eventually the solvent intermolecular modes. Also, strong hydrogen-bonding of HOD is found to significantly hinder the initial intramolecular energy transfer process.