An accurate classical simulation of a two-dimensional vibrational spectrum

OD stretch spectrum of a hydrated HOD molecule

Jonggu Jeon, Minhaeng Cho

Research output: Contribution to journalArticle

10 Citations (Scopus)

Abstract

An accurate computational method for the classical simulation of the two-dimensional vibrational spectra is presented. The method refines our previous computational method for the third order vibrational response function in the classical limit, and it enables capturing the diagonal elongation and its waiting time (T) dependence widely observed in experimental two-dimensional infrared (2D IR) spectra of intramolecular modes. The improvement is achieved by a series of new developments including (i) a block algorithm for the stability matrix computation, (ii) new equations of motion for the position-perturbed molecular dynamics (MD) trajectory, and (iii) enhanced sampling efficiency by exploiting the time-reversal invariance of MD trajectories. The method is applied to the simulation of 2D IR spectra of the OD stretch mode in a hydrated HOD molecule, employing a hybrid quantum mechanical/molecular mechanical force field. The simulated spectra exhibit diagonal elongation of the 2D IR signal at small T, reflecting the correlation of individual transitions among the inhomogeneously broadened ensemble. The slopes of the nodal lines of the elongated signals are found to decay with a time scale of 1.6 ps as T increases, in reasonable agreement with the frequency correlation decay time of 1.2 ps. The amplitudes of the positive and negative peaks also decay as T increases, due to vibrational population relaxation and molecular rotation. The peak positions tend to blue shift with increasing T, reflecting the different relaxation rates of the strongly and weakly solvated HOD species. These results indicate that the present method can reliably predict the waiting-time-dependent changes of 2D IR spectra of a single vibrational chromophore in solution.

Original languageEnglish
Pages (from-to)8148-8161
Number of pages14
JournalJournal of Physical Chemistry B
Volume118
Issue number28
DOIs
Publication statusPublished - 2014 Jul 17

Fingerprint

Vibrational spectra
Computational methods
vibrational spectra
Molecular dynamics
Elongation
Trajectories
Molecules
Stiffness matrix
Chromophores
Invariance
Equations of motion
molecules
simulation
elongation
Sampling
Infrared radiation
decay
trajectories
molecular dynamics
molecular rotation

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Materials Chemistry
  • Surfaces, Coatings and Films

Cite this

An accurate classical simulation of a two-dimensional vibrational spectrum : OD stretch spectrum of a hydrated HOD molecule. / Jeon, Jonggu; Cho, Minhaeng.

In: Journal of Physical Chemistry B, Vol. 118, No. 28, 17.07.2014, p. 8148-8161.

Research output: Contribution to journalArticle

@article{407c5979c1a24df6b44c493d9a218a7d,
title = "An accurate classical simulation of a two-dimensional vibrational spectrum: OD stretch spectrum of a hydrated HOD molecule",
abstract = "An accurate computational method for the classical simulation of the two-dimensional vibrational spectra is presented. The method refines our previous computational method for the third order vibrational response function in the classical limit, and it enables capturing the diagonal elongation and its waiting time (T) dependence widely observed in experimental two-dimensional infrared (2D IR) spectra of intramolecular modes. The improvement is achieved by a series of new developments including (i) a block algorithm for the stability matrix computation, (ii) new equations of motion for the position-perturbed molecular dynamics (MD) trajectory, and (iii) enhanced sampling efficiency by exploiting the time-reversal invariance of MD trajectories. The method is applied to the simulation of 2D IR spectra of the OD stretch mode in a hydrated HOD molecule, employing a hybrid quantum mechanical/molecular mechanical force field. The simulated spectra exhibit diagonal elongation of the 2D IR signal at small T, reflecting the correlation of individual transitions among the inhomogeneously broadened ensemble. The slopes of the nodal lines of the elongated signals are found to decay with a time scale of 1.6 ps as T increases, in reasonable agreement with the frequency correlation decay time of 1.2 ps. The amplitudes of the positive and negative peaks also decay as T increases, due to vibrational population relaxation and molecular rotation. The peak positions tend to blue shift with increasing T, reflecting the different relaxation rates of the strongly and weakly solvated HOD species. These results indicate that the present method can reliably predict the waiting-time-dependent changes of 2D IR spectra of a single vibrational chromophore in solution.",
author = "Jonggu Jeon and Minhaeng Cho",
year = "2014",
month = "7",
day = "17",
doi = "10.1021/jp501182d",
language = "English",
volume = "118",
pages = "8148--8161",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "28",

}

TY - JOUR

T1 - An accurate classical simulation of a two-dimensional vibrational spectrum

T2 - OD stretch spectrum of a hydrated HOD molecule

AU - Jeon, Jonggu

AU - Cho, Minhaeng

PY - 2014/7/17

Y1 - 2014/7/17

N2 - An accurate computational method for the classical simulation of the two-dimensional vibrational spectra is presented. The method refines our previous computational method for the third order vibrational response function in the classical limit, and it enables capturing the diagonal elongation and its waiting time (T) dependence widely observed in experimental two-dimensional infrared (2D IR) spectra of intramolecular modes. The improvement is achieved by a series of new developments including (i) a block algorithm for the stability matrix computation, (ii) new equations of motion for the position-perturbed molecular dynamics (MD) trajectory, and (iii) enhanced sampling efficiency by exploiting the time-reversal invariance of MD trajectories. The method is applied to the simulation of 2D IR spectra of the OD stretch mode in a hydrated HOD molecule, employing a hybrid quantum mechanical/molecular mechanical force field. The simulated spectra exhibit diagonal elongation of the 2D IR signal at small T, reflecting the correlation of individual transitions among the inhomogeneously broadened ensemble. The slopes of the nodal lines of the elongated signals are found to decay with a time scale of 1.6 ps as T increases, in reasonable agreement with the frequency correlation decay time of 1.2 ps. The amplitudes of the positive and negative peaks also decay as T increases, due to vibrational population relaxation and molecular rotation. The peak positions tend to blue shift with increasing T, reflecting the different relaxation rates of the strongly and weakly solvated HOD species. These results indicate that the present method can reliably predict the waiting-time-dependent changes of 2D IR spectra of a single vibrational chromophore in solution.

AB - An accurate computational method for the classical simulation of the two-dimensional vibrational spectra is presented. The method refines our previous computational method for the third order vibrational response function in the classical limit, and it enables capturing the diagonal elongation and its waiting time (T) dependence widely observed in experimental two-dimensional infrared (2D IR) spectra of intramolecular modes. The improvement is achieved by a series of new developments including (i) a block algorithm for the stability matrix computation, (ii) new equations of motion for the position-perturbed molecular dynamics (MD) trajectory, and (iii) enhanced sampling efficiency by exploiting the time-reversal invariance of MD trajectories. The method is applied to the simulation of 2D IR spectra of the OD stretch mode in a hydrated HOD molecule, employing a hybrid quantum mechanical/molecular mechanical force field. The simulated spectra exhibit diagonal elongation of the 2D IR signal at small T, reflecting the correlation of individual transitions among the inhomogeneously broadened ensemble. The slopes of the nodal lines of the elongated signals are found to decay with a time scale of 1.6 ps as T increases, in reasonable agreement with the frequency correlation decay time of 1.2 ps. The amplitudes of the positive and negative peaks also decay as T increases, due to vibrational population relaxation and molecular rotation. The peak positions tend to blue shift with increasing T, reflecting the different relaxation rates of the strongly and weakly solvated HOD species. These results indicate that the present method can reliably predict the waiting-time-dependent changes of 2D IR spectra of a single vibrational chromophore in solution.

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

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

U2 - 10.1021/jp501182d

DO - 10.1021/jp501182d

M3 - Article

VL - 118

SP - 8148

EP - 8161

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 28

ER -