Molecular dynamics simulation of nonlinear spectroscopies of intermolecular motions in liquid water

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water

A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.

2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions

The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.

Low-frequency anharmonic couplings in bromoform revealed from 2D Raman-THz spectroscopy: from the liquid to the crystalline phase

Two-dimensional (2D) Raman-THz spectroscopy in the frequency up to 7 THz has been applied to study the crystalline beta-phase of bromoform (CHBr3). As for liquid CHBr3, cross peaks are observed, which however sharpen up in the crystalline sample and split into assignable sub-contributions. In the Raman dimension, the frequency positions of these cross peaks coincide with the intramolecular bending modes of the CHBr3 molecules, and in the THz dimension with the IR active lattice modes of the crystal. This work expands the applicability of this new 2D spectroscopic technique to solid samples at cryogenic temperatures. Furthermore, it provides new experimental evidence that the cross peaks indeed originate from the coupling between intra- and intermolecular vibrational modes.

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

Surface Chemistry Can Unlock Drivers of Surface Stability of SARS-CoV-2 in a Variety of Environmental Conditions

The surface stability and resulting transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), specifically in indoor environments, have been identified as a potential pandemic challenge requiring investigation. This novel virus can be found on various surfaces in contaminated sites such as clinical places; however, the behavior and molecular interactions of the virus with respect to the surfaces are poorly understood. Regarding this, the virus adsorption onto solid surfaces can play a critical role in transmission and survival in various environments. In this article, we first give an overview of existing knowledge concerning viral spread, molecular structure of SARS-CoV-2, and the virus surface stability is presented. Then, we highlight potential drivers of the SARS-CoV-2 surface adsorption and stability in various environmental conditions. This theoretical analysis shows that different surface and environmental conditions including temperature, humidity, and pH are crucial considerations in building fundamental understanding of the virus transmission and thereby improving safety practices.

A Feynman diagram description of the 2D-Raman-THz response of amorphous ice

The 2D-Raman-THz response in all possible time-orderings (Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman) of amorphous water ice is calculated in two ways: from atomistic molecular dynamics simulations and with the help of a Feynman diagram model, the latter of which power-expands the potential energy surface and the dipole and polarizability surfaces up to leading order. Comparing both results allows one to dissect the 2D-Raman-THz response into contributions from mechanical anharmonicity, as well as electrical dipole and polarizability anharmonicities. Mechanical anharmonicity dominates the 2D-Raman-THz response of the hydrogen-bond stretching and hydrogen-bond bending bands of water, and dipole anharmonicity dominates that of the librational band, while the contribution of polarizability anharmonicity is comparably weak. A distinct echo of the hydrogen-bond stretching band is observed for the THz-Raman-THz pulse sequence, again dominated by mechanical anharmonicity. A peculiar mechanism is discussed, which is based on the coupling between the many normal modes within the hydrogen-bond stretching band and which will inevitably generate such an echo for an amorphous structure.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

Delocalization and stretch-bend mixing of the HOH bend in liquid water

Liquid water's rich sub-picosecond vibrational dynamics arise from the interplay of different high- and low-frequency modes evolving in a strong yet fluctuating hydrogen bond network. Recent studies of the OH stretching excitations of H₂O indicate that they are delocalized over several molecules, raising questions about whether the bending vibrations are similarly delocalized. In this paper, we take advantage of an improved 50 fs time-resolution and broadband infrared (IR) spectroscopy to interrogate the 2D IR lineshape and spectral dynamics of the HOH bending vibration of liquid H₂O. Indications of strong bend-stretch coupling are observed in early time 2D IR spectra through a broad excited state absorption that extends from 1500 cm^-1 to beyond 1900 cm^-1, which corresponds to transitions from the bend to the bend overtone and OH stretching band between 3150 and 3550 cm^-1. Pump-probe measurements reveal a fast 180 fs vibrational relaxation time, which results in a hot-ground state spectrum that is the same as observed for water IR excitation at any other frequency. The fastest dynamical time scale is 80 fs for the polarization anisotropy decay, providing evidence for the delocalized or excitonic character of the bend. Normal mode analysis conducted on water clusters extracted from molecular dynamics simulations corroborate significant stretch-bend mixing and indicate delocalization of δ_HOH on 2-7 water molecules.

Coherent X-rays reveal the influence of cage effects on ultrafast water dynamics

The dynamics of liquid water feature a variety of time scales, ranging from extremely fast ballistic-like thermal motion, to slower molecular diffusion and hydrogen-bond rearrangements. Here, we utilize coherent X-ray pulses to investigate the sub-100 fs equilibrium dynamics of water from ambient conditions down to supercooled temperatures. This novel approach utilizes the inherent capability of X-ray speckle visibility spectroscopy to measure equilibrium intermolecular dynamics with lengthscale selectivity, by measuring oxygen motion in momentum space. The observed decay of the speckle contrast at the first diffraction peak, which reflects tetrahedral coordination, is attributed to motion on a molecular scale within the first 120 fs. Through comparison with molecular dynamics simulations, we conclude that the slowing down upon cooling from 328 K down to 253 K is not due to simple thermal ballistic-like motion, but that cage effects play an important role even on timescales over 25 fs due to hydrogen-bonding.

Nuclear quantum effects on the vibrational dynamics of liquid water

Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm^-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.

Ultrafast Vibrational Echo Spectroscopy of Liquid Water from First-Principles Simulations

Vibrational echo spectroscopy has become a powerful technique to study vibrational spectral diffusion in water and aqueous solutions. The dynamics of vibrational spectral diffusion is intimately related to the hydrogen bond fluctuations in liquid water and other hydrogen bonded liquids. Earlier theoretical calculations of vibrational echo spectroscopy of aqueous systems were based on classical molecular dynamics simulations involving empirical force fields of water. In the current work, we have employed the method of ab initio molecular dynamics simulation to calculate the spectral observables of vibrational echo and two-dimensional infrared (2D-IR) spectroscopy of liquid water at room temperature under Condon and cumulant approximations. The time scales extracted from the temporal decay of the frequency-time correlation function (FTCF), short-time slope of three pulse photon echo (SP3E), dynamic line width (DLW), and the slope of nodal line of 2D-IR spectra are found to be in reasonably close agreement with each other which reinforces the assertion that signatures of FTCF can be captured using three pulse photon echo and 2D-IR spectroscopy.

Ultrafast dynamics of liquid water: frequency fluctuations of the OH stretch and the HOH bend

Frequency fluctuations of the OH stretch and the HOH bend in liquid water are reported from the third-order response function evaluated using the TTM3-F potential for water. The simulated two-dimensional infrared spectra of the OH stretch are similar to previously reported theoretical results. The present study suggests that the frequency fluctuation of the HOH bend is faster than that of the OH stretch. The ultrafast loss of the frequency correlation of the HOH bend is due to the strong couplings with the OH stretch as well as the intermolecular hydrogen bond bend.

Instantaneous normal mode analysis for intermolecular and intramolecular vibrations of water from atomic point of view

By exploiting the instantaneous normal mode (INM) analysis for models of flexible molecules, we investigate intermolecular and intramolecular vibrations of water from the atomic point of view. With two flexible SPC/E models, our investigations include three aspects about their INM spectra, which are separated into the unstable, intermolecular, bending, and stretching bands. First, the O- and H-atom contributions in the four INM bands are calculated and their stable INM spectra are compared with the power spectra of the atomic velocity autocorrelation functions. The unstable and intermolecular bands of the flexible models are also compared with those of the SPC/E model of rigid molecules. Second, we formulate the inverse participation ratio (IPR) of the INMs, respectively, for the O- and H-atom and molecule. With the IPRs, the numbers of the three species participated in the INMs are estimated so that the localization characters of the INMs in each band are studied. Further, by the ratio of the IPR of the H atom to that of the O atom, we explore the number of involved OH bond per molecule participated in the INMs. Third, by classifying simulated molecules into subensembles according to the geometry of their local environments or their H-bond configurations, we examine the local-structure effects on the bending and stretching INM bands. All of our results are verified to be insensible to the definition of H-bond. Our conclusions about the intermolecular and intramolecular vibrations in water are given.

Fluctuations and Relaxation Dynamics of Liquid Water Revealed by Linear and Nonlinear Spectroscopy

Many efforts have been devoted to elucidating the intra- and intermolecular dynamics of liquid water because of their important roles in many fields of science and engineering. Nonlinear spectroscopy is a powerful tool to investigate the dynamics. Because nonlinear response functions are described by more than one time variable, it is possible to analyze static and dynamic mode couplings. Here we review the intra- and intermolecular dynamics of liquid water revealed by recent linear and nonlinear spectroscopic experiments and computer simulations. In particular, we discuss the population relaxation, anisotropy decay, and spectral diffusion of the intra- and intermolecular motions of water and their temperature dependence, which play important roles in ultrafast dynamics and relaxations in water.

Time-averaging approximation in the interaction picture: absorption line shapes for coupled chromophores with application to liquid water

The time-averaging approximation (TAA), originally developed to calculate vibrational line shapes for coupled chromophores using mixed quantum/classical methods, is reformulated. In the original version of the theory, time averaging was performed for the full one-exciton Hamiltonian, while herein the time averaging is performed on the coupling (off-diagonal) Hamiltonian in the interaction picture. As a result, the influence of the dynamic fluctuations of the transition energies is more accurately described. We compare numerical results of the two versions of the TAA with numerically exact results for the vibrational absorption line shape of the OH stretching modes in neat water. It is shown that the TAA in the interaction picture yields theoretical line shapes that are in better agreement with exact results.

Three-dimensional infrared spectroscopy of isotope-substituted liquid water reveals heterogeneous dynamics

The dynamics of the hydrogen bond network of isotopically substituted liquid water are investigated with a new ultrafast nonlinear vibrational spectroscopy, three-dimensional infrared spectroscopy (3D-IR). The 3D-IR spectroscopy is sensitive to three-point frequency fluctuation correlation functions, and the measurements reveal heterogeneous structural relaxation dynamics. We interpret these results as subensembles of water which do not interconvert on a half picosecond time scale. We connect the experimental results to molecular dynamics (MD) simulations, performing a line shape analysis as well as complex network analysis.

Effects of nonadditive interactions on ion solvation at the water/vapor interface: a molecular dynamics study

The solvation of halide ions at the water/vapor interface is investigated by using molecular dynamics simulations with nonpolarizable molecular mechanical (MM), polarizable MM, and quantum mechanical (QM)/MM methods. The free energy profile of the ion solvation is decomposed into the energy and the entropic contributions along the ion displacement from inside to the surface of water. It is found that the surface affinity of the ion, relative to the bulk value, is determined by a subtle balance between the energetic destabilization and the entropic stabilization with the ion displacement. The amount of energetic destabilization is found to be reduced when nonadditive interactions are included, as in the polarizable MM and QM/MM models. The structure of water around the ion at the interface is also largely modified when the higher order effects are considered. For example, the induced dipole effect enhances the solvation structure around the ion at the interface significantly and thus reduces the amount of entropic stabilization at the interface, relative to in the bulk. It is found that this induced dipole effect causes the slowing in the ion-water hydrogen bond dynamics at the interface. On the other hand, the higher order induced multipole effects in the QM/MM method suppress both the excessive enhancement of the solvation structure and the slowing of the ion-water hydrogen bond dynamics at the interface. The present study demonstrates that not only the induced dipole moment but also the higher order induced multipole moments, which are neglected in standard empirical models, are essential for the correct description of the ion solvation at the water/vapor interface.

Vibrational sum-frequency generation spectroscopy at the water/lipid interface: molecular dynamics simulation study

The sum-frequency generation (SFG) spectrum from the water/[1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine] (DMPC) interface in the OH stretching mode region of water is simulated and shows three spectral peaks which are assigned to different environments. The weak 3590 cm(-1) peak is attributed to a few water molecules coupled to the glycerol backbone of DMCP. The 3470 cm(-1) feature comes from the top water layer adjacent to the hydrophilic headgroup of DMPC. The 3290 cm(-1) peak arises from the near-bulk water nonadjacent to DMPC. The stretching mode corresponding to the 3290 cm(-1) peak is strongly coupled to the neighboring water molecules. In contrast, the 3470 cm(-1) mode is decoupled from the surrounding water molecules, and the orientation of water is governed by DMPC. This decoupling explains the slow relaxation dynamics of water measured in the time-resolved SFG experiment. Despite the similarity of the SFG spectra, the peak origins of water/lipid and water/vapor interfaces are different.