Perturbation approach for computing infrared spectra of the local mode of probe molecules

Rui Jie Xue, Adam Grofe, He Yin, Zexing Qu, Jiali Gao, Hui Li

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

Original languageEnglish (US)
Pages (from-to)191-201
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume13
Issue number1
DOIs
StatePublished - Jan 1 2017

Fingerprint

Vibrational spectra
Stretching
Molecular dynamics
infrared spectra
Infrared radiation
perturbation
Molecules
probes
Hydrogen bonds
molecular dynamics
molecules
Potential energy surfaces
Biological systems
Chromophores
electrostatics
Wave functions
Acetone
Coulomb interactions
Computational efficiency
Electrostatics

Cite this

Perturbation approach for computing infrared spectra of the local mode of probe molecules. / Xue, Rui Jie; Grofe, Adam; Yin, He; Qu, Zexing; Gao, Jiali; Li, Hui.

In: Journal of Chemical Theory and Computation, Vol. 13, No. 1, 01.01.2017, p. 191-201.

Research output: Contribution to journalArticle

Xue, Rui Jie ; Grofe, Adam ; Yin, He ; Qu, Zexing ; Gao, Jiali ; Li, Hui. / Perturbation approach for computing infrared spectra of the local mode of probe molecules. In: Journal of Chemical Theory and Computation. 2017 ; Vol. 13, No. 1. pp. 191-201.
@article{5ace1ca96c494203a595d8502a432c73,
title = "Perturbation approach for computing infrared spectra of the local mode of probe molecules",
abstract = "Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.",
author = "Xue, {Rui Jie} and Adam Grofe and He Yin and Zexing Qu and Jiali Gao and Hui Li",
year = "2017",
month = "1",
day = "1",
doi = "10.1021/acs.jctc.6b00733",
language = "English (US)",
volume = "13",
pages = "191--201",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "1",

}

TY - JOUR

T1 - Perturbation approach for computing infrared spectra of the local mode of probe molecules

AU - Xue, Rui Jie

AU - Grofe, Adam

AU - Yin, He

AU - Qu, Zexing

AU - Gao, Jiali

AU - Li, Hui

PY - 2017/1/1

Y1 - 2017/1/1

N2 - Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

AB - Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

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

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

U2 - 10.1021/acs.jctc.6b00733

DO - 10.1021/acs.jctc.6b00733

M3 - Article

VL - 13

SP - 191

EP - 201

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 1

ER -