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

Jonggu Jeon, Joon Hyung Lim, Seongheun Kim, Heejae Kim, Minhaeng Cho

Research output: Contribution to journalArticle

8 Citations (Scopus)

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.

Original languageEnglish
Pages (from-to)5356-5367
Number of pages12
JournalJournal of Physical Chemistry A
Volume119
Issue number21
DOIs
Publication statusPublished - 2015 May 28

Fingerprint

System theory
molecular relaxation
Kinetic energy
kinetic energy
Water
Fourier transformation
spectral resolution
Relaxation time
Energy transfer
water
Molecular dynamics
relaxation time
energy transfer
molecular dynamics
Energy Transfer
Molecular Dynamics Simulation
Spectral density
Spectral resolution
Computer simulation
Computational methods

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

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. / Jeon, Jonggu; Lim, Joon Hyung; Kim, Seongheun; Kim, Heejae; Cho, Minhaeng.

In: Journal of Physical Chemistry A, Vol. 119, No. 21, 28.05.2015, p. 5356-5367.

Research output: Contribution to journalArticle

@article{b031c24accb7465e813a1858015150a8,
title = "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",
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.",
author = "Jonggu Jeon and Lim, {Joon Hyung} and Seongheun Kim and Heejae Kim and Minhaeng Cho",
year = "2015",
month = "5",
day = "28",
doi = "10.1021/jp510157y",
language = "English",
volume = "119",
pages = "5356--5367",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "21",

}

TY - JOUR

T1 - Simultaneous spectral and temporal analyses of kinetic energies in nonequilibrium systems

T2 - Theory and application to vibrational relaxation of O-D stretch mode of HOD in water

AU - Jeon, Jonggu

AU - Lim, Joon Hyung

AU - Kim, Seongheun

AU - Kim, Heejae

AU - Cho, Minhaeng

PY - 2015/5/28

Y1 - 2015/5/28

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=84930681551&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84930681551&partnerID=8YFLogxK

U2 - 10.1021/jp510157y

DO - 10.1021/jp510157y

M3 - Article

AN - SCOPUS:84930681551

VL - 119

SP - 5356

EP - 5367

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 21

ER -