Hybrid quantum/classical simulation of the vibrational relaxation of the bend fundamental in liquid water

The Bending Mode of Water: A Powerful Probe for Hydrogen Bond Structure of Aqueous Systems

Insights into the microscopic structure and dynamics of the water's hydrogen-bonded network are crucial to understand the role of water in biology, atmospheric and geochemical processes, and chemical reactions in aqueous systems. Vibrational spectroscopy of water has provided many such insights, in particular using the O-H stretch mode. In this Perspective, we summarize our recent studies that have revealed that the H-O-H bending mode can be an equally powerful reporter for the microscopic structure of water and provides more direct access to the hydrogen-bonded network than the conventionally studied O-H stretch mode. We discuss the fundamental vibrational properties of the water bending mode, such as the intermolecular vibrational coupling, and its effects on the spectral lineshapes and vibrational dynamics. Several examples of static and ultrafast bending mode spectroscopy illustrate how the water bending mode provides an excellent window on the microscopic structure of both bulk and interfacial water.

A Microscopic Interpretation of Pump-Probe Vibrational Spectroscopy Using Ab Initio Molecular Dynamics

What happens when extra vibrational energy is added to water? Using nonequilibrium molecular dynamics simulations, also including the full electronic structure, and novel descriptors, based on projected vibrational density of states, we are able to follow the flow of excess vibrational energy from the excited stretching and bending modes. We find that the energy relaxation, mostly mediated by a stretching-stretching coupling in the first solvation shell, is highly heterogeneous and strongly depends on the local environment, where a strong hydrogen bond network can transport energy with a time scale of 200 fs, whereas a weaker network can slow down the transport by a factor 2-3.

Molecular origin of the difference in the HOH bend of the IR spectra between liquid water and ice

The intensity of the HOH bend in the infrared (IR) spectrum of ice is significantly smaller than the corresponding one in liquid water. This difference in the IR intensities of the HOH bend in the two systems is investigated using Molecular Dynamics (MD) simulations with the flexible, polarizable, ab initio based TTM3-F model for water, a potential that correctly reproduces the experimentally observed increase of the HOH angle in liquid water and ice from the water monomer value. We have identified two factors that are responsible for the difference in the intensity of the HOH bend in liquid water and ice: (i) the decrease of the intensity of the HOH bend in ice caused by the strong anti-correlation between the permanent dipole moment of a molecule and the induced dipole moment of neighboring hydrogen bond acceptor molecules, and (ii) the weakening of this anti-correlation by the disordered hydrogen bond network in liquid water. The presence of the anti-correlation in ice is further confirmed by ab initio electronic structure calculations of water pentamer clusters extracted from the trajectories of the MD simulations with the TTM3-F potential for ice and liquid water.

Juxtaposing density matrix and classical path-based wave packet dynamics

In many physical, chemical, and biological systems energy and charge transfer processes are of utmost importance. To determine the influence of the environment on these transport processes, equilibrium molecular dynamics simulations become more and more popular. From these simulations, one usually determines the thermal fluctuations of certain energy gaps, which are then either used to perform ensemble-averaged wave packet simulations, also called Ehrenfest dynamics, or to employ a density matrix approach via spectral densities. These two approaches are analyzed through energy gap fluctuations that are generated to correspond to a predetermined spectral density. Subsequently, density matrix and wave packet simulations are compared through population dynamics and absorption spectra for different parameter regimes. Furthermore, a previously proposed approach to enforce the correct long-time behavior in the wave packet simulations is probed and an improvement is proposed.

Hybrid quantum/classical simulations of the vibrational relaxation of the amide I mode of N-methylacetamide in D2O solution

Hybrid quantum/classical molecular dynamics (MD) is applied to simulate the vibrational relaxation (VR) of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution. A novel version of the vibrational molecular dynamics with quantum transitions (MDQT) treatment is developed in which the amide I mode is treated quantum mechanically while the remaining degrees of freedom are treated classically. The instantaneous normal modes of the initially excited NMAD molecule (INM(0)) are used as internal coordinates since they provide a proper initial partition of the system in quantum and classical subsystems. The evolution in time of the energy stored in each individual normal mode is subsequently quantified using the hybrid quantum-classical instantaneous normal modes (INM(t)). The identities of both the INM(0)s and the INM(t)s are tracked using the equilibrium normal modes (ENMs) as templates. The results extracted from the hybrid MDQT simulations show that the quantum treatment of the amide I mode accelerates the whole VR process versus pure classical simulations and gives better agreement with experiments. The relaxation of the amide I mode is found to be essentially an intramolecular vibrational redistribution (IVR) process with little contribution from the solvent, in agreement with previous theoretical and experimental studies. Two well-defined relaxation mechanisms are identified. The faster one accounts for ≈40% of the total vibrational energy that flows through the NMAD molecule and involves the participation of the lowest frequency vibrations as short-life intermediate modes. The second and slower mechanism accounts for the remaining ≈60% of the energy released and is associated to the energy flow through specific mid-range and high-frequency modes.