Isotope Effects in 2D correlation infrared Spectra of Water: HEOM Analysis of Molecular Dynamics-Based Machine Learning Models

TL;DR

This study employs the HEOM framework to simulate and analyze isotope effects in 2D infrared spectra of water, revealing detailed insights into vibrational dynamics and environmental interactions in H2O and D2O.

Contribution

It introduces a non-perturbative HEOM approach to accurately model 2D IR spectra of water isotopes, capturing complex environmental and anharmonic effects.

Findings

HEOM accurately reproduces experimental 2D spectra.

Differences in isotope effects elucidate vibrational relaxation mechanisms.

Environmental interactions significantly influence spectral features.

Abstract

We model, simulate, and analyze the intramolecular modes of liquid H2O and D2O to elucidate how energy excitation, relaxation, and vibrational dephasing interplay through anharmonic mode-mode coupling. Our analysis employs two-dimensional (2D) correlation spectra, a representative observable in nonlinear infrared vibrational spectroscopy. Accurate reproduction of these 2D spectral profiles requires not only a precise dynamical description of intramolecular vibrations but also an appropriate treatment of thermal environmental effects arising from strong interactions with surrounding molecules, which act as thermal baths. Capturing the essential features of the 2D spectra further demands a non-Markovian, non-perturbative, and nonlinear description of the interactions between intramolecular modes and their baths. To this end, we adopt a hierarchical equations of motion (HEOM) framework to compute the 2D spectra. By comparing the resulting spectra of H2O and D2O, we explore the underlying mechanisms governing their complex energy and phase relaxation dynamics.