Signatures of coherent vibrational energy transfer in IR and Raman line shapes for liquid water

The isosbestic point in the Raman spectra of the hydration shell

Isosbestic point is often observed in a series of spectra, but their interpretation is still controversial, such as whether the continuum model can produce an isosbestic point. In order to answer this question, the Raman spectra of hydration shell with continuous distribution structure in different ionic aqueous solutions were separated by Raman ratio spectra, and an isosbestic point was successfully observed. Our experimental results show that the continuum model can indeed produce the isosbestic point. In order to deepen the understanding of the isosbestic point, we calculate the first moment of the Raman spectra and conduct molecular dynamics (MD) simulations. Both experimental and theoretical findings indicate that elevated temperatures lead to increased disorder among water molecules within the hydration shell.

Modeling Vibrational Sum Frequency Generation Spectra of Interfacial Water on a Gold Surface: The Role of the Fermi Resonance

Studying the hydrogen bonding structure of H2O at the metal-water interface is a highly complex yet fascinating endeavor. The intricate interactions and diverse orientations of water molecules on metal surfaces with varying potentials pose a significant challenge in elucidating the coupling between O-H stretching and H-O-H bending modes. In this study, we employed DFT-MD simulation to explore how the orientation of interfacial water molecules changes with the applied potential on the Au(111) surface. Based on the surface-specific velocity-velocity correlation function (ssVVCF) formula, we calculated vibrational sum frequency generation (VSFG) spectra for the O-H stretches. We found that three assigned peaks (∼3300, ∼3450, and 3650 cm-1) shifted toward lower frequencies as the potential moved toward more negative values. Our results align remarkably well with experimental Raman spectroscopy data. Notably, our VSFG analysis revealed a significant change in the VSFG spectra of the hydrogen-bonded O-H groups (∼3300 cm-1), switching from a negative to a positive sign with decreasing potential. This alteration suggests a substantial change in the orientation of these low-frequency O-H groups owing to their increased interactions with the Au surface. In contrast, the orientations of both the high-frequency O-H groups (∼3450 cm-1) and the dangling O-H groups (∼3650 cm-1) remained unaffected by the applied potentials. Furthermore, our analysis of the decomposed vibrational density of states (VDOS) for the H-O-H bending mode uncovered the coupling between the H-O-H bending and O-H stretching vibrations, known as the Fermi resonance. Our work suggests that the H-O-H bending vibration becomes restricted when water molecules transition from the ″one-H-down″ to the ″two-H-down″ conformation, leading to a redshift in the O-H stretching vibration through the Fermi resonance. By constructing the VSFG and decomposed VDOS spectra, we gained valuable insights into the structural changes that Raman spectra alone cannot fully interpret. Specifically, our analysis revealed the critical role of the Fermi resonance effect in shaping the spectroscopic signature of interfacial water molecules on the Au(111) surface.

Probing the structure of water in individual living cells

Water regulates or even governs a wide range of biological processes. Despite its fundamental importance, surprisingly little is known about the structure of intracellular water. Herein we employ a Raman micro-spectroscopy technique to uncover the composition, abundance and vibrational spectra of intracellular water in individual living cells. In three different cell types, we show a small but consistent population (~3%) of non-bulk-like water. It exhibits a weakened hydrogen-bonded network and a more disordered tetrahedral structure. We attribute this population to biointerfacial water located in the vicinity of biomolecules. Moreover, our whole-cell modeling suggests that all soluble (globular) proteins inside cells are surrounded by, on average, one full molecular layer (about 2.6 Angstrom) of biointerfacial water. Furthermore, relative invariance of biointerfacial water is observed among different single cells. Overall, our study not only opens up experimental possibilities of interrogating water structure in vivo but also provides insights into water in life.

On selection rules in two-dimensional terahertz-infrared-visible spectroscopy

Two-dimensional terahertz-infrared-visible (2D TIRV) spectroscopy directly measures the coupling between quantum high-frequency vibrations and classical low-frequency modes of molecular motion. In addition to coupling strength, the signal intensity in 2D TIRV spectroscopy can also depend on the selection rules of the excited transitions. Here, we explore the selection rules in 2D TIRV spectroscopy by studying the coupling between the high-frequency CH3 stretching and low-frequency vibrations of liquid dimethyl sulfoxide (DMSO). Different excitation pathways are addressed using variations in laser pulse timing and different polarizations of exciting pulses and detected signals. The DMSO signals generated via different excitation pathways can be readily distinguished in the spectrum. The intensities of different excitation pathways vary unequally with changes in polarization. We explain how this difference stems from the intensities of polarized and depolarized Raman and hyper-Raman spectra of high-frequency modes. These results apply to various systems and will help design and interpret new 2D TIRV spectroscopy experiments.

Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model

The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer. We find that our model is capable of reproducing the bimodal nature of the main peak, with the shoulder at 3250 cm-1 resulting almost entirely from Fermi resonance. Furthermore, we capture the temperature and polarization dependence of the shoulder, which has proven to be difficult to obtain with previous methods, and analyze the origin of this dependence. We expect our model to be generally useful for understanding and predicting how Raman spectra change under different conditions and with different probe reporters beyond water.

Simulating the isotropic Raman spectra of O-H stretching mode in liquid H2O based on a machine learning potential: the influence of vibrational couplings

In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.

Vibrational Coupling and Hydrogen-Bond Structure of Water in the Extended Hydration Shell of Metal Ions

Intra- and intermolecular vibrational coupling (VC) and hydrogen bonding (H-bonding) of water are sparsely understood in the hydration shell (HS) of a metal ion, though the corresponding knowledge for an anion is quite extensive. This is primarily due to the overwhelming effect of anions on water, which masks the subtle perturbing influence of most of the cations. Using Raman difference spectroscopy with simultaneous curve fitting (Raman-DS-SCF) in combination with isotopic dilution and polarized Raman spectroscopy, we have elucidated the VC and H-bonding of water in the HS of bi- and trivalent metal ions─Mg2+, Ca2+, La3+, Gd3+, Dy3+. Polarized Raman measurement of the HS water with VC "turned on" and "turned off" (using isotopically diluted water, HOD) reveals that water retains the intra- and intermolecular vibrational coupling in the HS of high-charge-density metal ions, which is in stark contrast to that of an anion. Hydration shell spectroscopy in HOD unambiguously shows that the average H-bonding of water becomes stronger in the HS than that of bulk water. The first HS water strongly donates two H-bonds to the second HS water (ν̅max ≈ 3200 cm-1) but weakly accepts a H-bond from the second HS water (ν̅max ≈ 3590 cm-1), which makes the HS water heterogeneous in terms of its H-bond structure. The weakly interacting OH (ν̅max 3585 cm-1 in HOD) red-shifts by ∼ 15 cm-1 while the VC is "turned on" (ν̅max ≈ 3600 cm-1 in H2O), revealing the intramolecular coupling of water in the HS of metal ions.

Dispersion related coupling effects in IR spectra on the example of water and Amide I bands

We discuss coupling effects in infrared spectra which are caused by dispersion and local field effects. The first effect is instigated by changes of the refractive index due to absorption which have an impact on the strength of adjacent absorptions. The second effect is a consequence of the light-induced polarization of one molecule affecting neighboring ones. These coupling effects do not only effect band positions, but also influence relative intensities. They are particularly strong in case of overlapping bands and complicate their deconvolution by band fitting. We investigated the corresponding challenges for the HO-stretching vibrations in water and the Amide I band in proteins. Our findings show that the effects are significant and of high interest for protein and water structure determination. Especially, for the water stretching vibrations we conclude that it is of utmost importance to consider such coupling effects in quantum mechanical calculations of water spectra. Otherwise, progress in understanding band positions and profiles is likely to be hampered. Also, in case of the Amide I band we found a distinct impact of such coupling effects. Accordingly, we strongly recommend consideration of dispersion and local field effects to ensure the possibility of an accurate, quantitative determination of α-helix and β-sheet structures.

Electron Impact with the Liquid-Vapor Interface

Reaction dynamics in the liquid-vapor interface is one of the crucial physical sciences but is still starving for in-depth exploration. It is challenging to selectively detect the interfacial species or the yields of chemical reaction therein, meanwhile shielding or reducing the interference from the vapor and liquid bulk. Mass spectrometry is a straightforward method but is also frustrated in such a selective detection. Using a liquid microjet in combination with a pulsed electron beam, a linear time-of-flight mass spectrometer, and a quadrupole mass filter, we recently innovated time-delayed mass spectrometry for investigations of the liquid-vapor interface. In this Account, we illustrate how this unique method succeeds in disentangling different sources, i.e., the vapor and liquid-vapor interface, of the ionic yields of the electron impacts with a liquid beam of alcohol in vacuum. These achievements are basically attributed to the application of an onion-peeling strategy in the ion detection. Concretely, the microsecond time scale of molecular volatilization can be resolved well by tuning the delay time between the nanosecond pulses of incident electron bunch and ion attractor. First, the specific orientation of the interfacial molecule, i.e., a well-known fact about the hydrophobic hydrocarbon groups pointing outside the liquid surface of alcohol, is validated again. More importantly, the dynamic features of time-delayed mass spectra, in particular, for the ionic yields from the liquid-vapor interface, are rationalized explicitly. Moreover, we demonstrate evidence of in situ molecular dimers in the liquid-vapor interface of 1-propanol. As the first example of electron-induced reaction in the liquid-vapor interface, dimethyl ether can be synthesized in the liquid methanol interface due to local interfacial acidification by high-energy electron impacts. On the contrary, the low energy electron can lead to local basicity through dissociative electron attachment (DEA). Besides the primary low-energy electrons, the low-energy secondary and inelastically scattered electrons in the higher-energy impacts of the primary electrons can also participate in the DEA process. In contrast to the gas- or solid-phase DEAs, that in the liquid-vapor interface shows distinct differences in both the types and efficiencies of anionic products. With these and efforts in the future, we develop a molecular-level understanding of how the chemical reactions happen in the liquid-vapor interface.

Local Ice-like Structure at the Liquid Water Surface

Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.

Isotope effects on the structural transformation and relaxation of deeply supercooled water

We have examined the structure of supercooled liquid D2O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D2O agree remarkably well with those for H2O.

Hydration and Hydrogen Bond Order of Octadecanoic Acid and Octadecanol Films on Water at 21 and 1 °C

The temperature-dependent hydration structure of long-chain fatty acids and alcohols at air-water interfaces has great significance in the fundamental interactions underlying ice nucleation in the atmosphere. We present an integrated theoretical and experimental study of the temperature-dependent vibrational structure and electric field character of the immediate hydration shells of fatty alcohol and acid headgroups. We use a combination of surface-sensitive infrared reflection-absorption spectroscopy (IRRAS), surface potentiometry, and ab initio molecular dynamics simulations to elucidate detailed molecular structures of the octadecanoic acid and octadecanol (stearic acid and stearyl alcohol) headgroup hydration shells at room temperature and near freezing. In experiments, the alcohol at high surface concentration exhibits the largest surface potential; yet we observe a strengthening of the hydrogen-bonding for the solvating water molecules near freezing for both the alcohol and the fatty acid IRRAS experiments. Results reveal that the hydration shells for both compounds screen their polar headgroup dipole moments reducing the surface potential at low surface coverages; at higher surface coverage, the polar headgroups become dehydrated, which reduces the screening, correlating to higher observed surface potential values. Lowering the temperature promotes tighter chain packing and an increase in surface potential. IRRAS reveals that the intra- and intermolecular vibrational coupling mechanisms are highly sensitive to changes in temperature. We find that intramolecular coupling dominates the vibrational relaxation pathways for interfacial water determined by comparing the H₂O and the HOD spectra. Using ab initio molecular dynamics (AIMD) calculations on cluster systems of propanol + 6H₂O and propionic acid + 10H₂O, a spectral decomposition scheme was used to correlate the OH stretching motion with the IRRAS spectral features, revealing the effects of intra- and intermolecular coupling on the spectra. Spectra calculated with AIMD reproduce the red shift and increase in intensity observed in experimental spectra corresponding to the OH stretching region of the first solvation shell. These findings suggest that intra- and intermolecular vibrational couplings strongly impact the OH stretching region at fatty acid and fatty alcohol water interfaces. Overall, results are consistent with ice templating behavior for both the fatty acid and the alcohol, yet the surface potential signature is strongest for the fatty alcohol. These findings develop a better understanding of the complex surface potential and spectral signatures involved in ice templating.

Vibrational exciton delocalization precludes the use of infrared intensities as proxies for surfactant accumulation on aqueous surfaces

Surface-sensitive vibrational spectroscopy is a common tool for measuring molecular organization and intermolecular interactions at interfaces. Peak intensity ratios are typically used to extract molecular information from one-dimensional spectra but vibrational coupling between surfactant molecules can manifest as signal depletion in one-dimensional spectra. Through a combination of experiment and theory, we demonstrate the emergence of vibrational exciton delocalization in infrared reflection–absorption spectra of soluble and insoluble surfactants at the air/water interface. Vibrational coupling causes a significant decrease in peak intensities corresponding to C–F vibrational modes of perfluorooctanoic acid molecules. Vibrational excitons also form between arachidic acid surfactants within a compressed monolayer, manifesting as signal reduction of C–H stretching modes. Ionic composition of the aqueous phase impacts surfactant intermolecular distance, thereby modulating vibrational coupling strength between surfactants. Our results serve as a cautionary tale against employing alkyl and fluoroalkyl vibrational peak intensities as proxies for concentration, although such analysis is ubiquitous in interface science.

Coupling between surfactant molecules at the air/water interface bleeds intensity into a diffuse background, such that single-wavelength vibrational intensity is effectively depleted at high surface coverage.

Restructuring of Hydration Shell Water due to Solvent-Shared Ion Pairing (SSIP): A Case Study of Aqueous MgCl2 and LaCl3 Solutions

Hydration of ions plays a crucial role in interionic interactions and associated processes in aqueous media, but selective probing of the hydration shell water is nontrivial. Here, we introduce Raman difference with simultaneous curve fitting (RD-SCF) analysis to extract the OH-stretch spectrum of hydration shell water, not only for the fully hydrated ions (Mg²⁺, La³⁺, and Cl^-) but also for the ion pairs. RD-SCF analyses of diluted MgCl₂ (0.18 M) and LaCl₃ (0.12 M) solutions relative to aqueous NaCl of equivalent Cl^- concentrations provide the OH-stretch spectra of water in the hydration shells of fully hydrated Mg²⁺ and La³⁺ cations relative to that of Na⁺. Integrated intensities of the hydration shell spectra of Mg²⁺ and La³⁺ ions increase linearly with the salt concentration (up to 2.0 M MgCl₂ and 1.3 M LaCl₃), which suggests no contact ion pair (CIP) formation in the MgCl₂ and LaCl₃ solutions. Nevertheless, the band shapes of the cation hydration shell spectra show a growing signature of Cl^--associated water with the rising salt concentration, which is a manifestation of the formation of a solvent-shared ion pair (SSIP). The OH-stretch spectrum of the shared/intervening water in the SSIP, retrieved by second-round RD-SCF analysis (2RD-SCF), shows that the average H-bonding of the shared water is weaker than that of the hydration water of the fully hydrated cation (Mg²⁺ or La³⁺) but stronger than that of the anion (Cl^-). The shared water displays an overall second-order dependence on the concentration of the interacting ions, unveiling 1:1 stoichiometry of the SSIP formed between Mg²⁺ and Cl^- as well as La³⁺ and Cl^-.

Hydration Shell Changes in Surfactant Aggregate Transitions Revealed by Raman-MCR Spectroscopy

Hydration states of many self-assemblies directly relate to their structures and functions. Here, we use Raman multivariate curve resolution (Raman-MCR) assisted by differential scanning calorimetry and nuclear magnetic resonance to explore the hydration properties of aggregates formed by three cationic ammonium surfactants, trimethylene-1,3-bis(dodecyldimethylammonium bromide) (12-3-12(Br)₂), didodecyldimethylammonium bromide (DDAB), and dodecyltrimethylammonium bromide (DTAB). For 12-3-12(Br)₂, the transitions from spherical to rodlike and wormlike micelles lead to about 20% and 60% dehydration and gradually weaken water tetrahedral order and H-bond in hydration shells for both headgroup and hydrophobic chain. As to DDAB, unilamellar vesicles contain two kinds of hydration water species, but multicompartment vesicles exhibit decreased water order and weaker H-bond. DTAB only forms spherical micelles and its hydration structure is similar to that of the 12-3-12(Br)₂ spherical micelles. This work provides a basis to explore the hydration states of complex biological self-assemblies.

Vibrational coupling to hydration shell - Mechanism to performance enhancement of qualitative analysis in NIR spectroscopy of carbohydrates in aqueous environment

Vibrational coupling between carbohydrates and the hydration shell is unveiled as the underlying mechanism that improves wavenumber-selectively the carbohydrate discrimination performance by near-infrared (NIR) spectroscopy. The investigation is based on measurement of six carbohydrates (fructose, glucose, mannose, ribose, xylose and sorbitol) in aqueous solution in different concentration levels (5 mg/L, ~0.03 mmol/dm^-3 and 20 mg/L, ~0.1 mmol/dm^-3). The results of multivariate classification are interpreted by quantum mechanical NIR spectra simulations. The simulation unveils that the phenomenon is vibration-selective and thus wavenumber-selective, and leads to an enhancement of the qualitative information contained in the specific spectral regions. The location of these regions and the related performance correspond fully to the appearance and magnitude of the unveiled cooperative vibration effect.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Raman spectrum and polarizability of liquid water from deep neural networks

Using deep neural networks to model the polarizability and potential energy surfaces, we compute the Raman spectrum of liquid water at several temperatures with ab initio molecular dynamics accuracy.

Electrically induced liquid-liquid phase transition in water at room temperature

In this work we expand on findings previously reported [Wexler et al., Phys. Chem. Chem. Phys., 2016, 18, 16281] on the experimental observation of a phase transition in a hydrogen bonded liquid manifesting in long range dipole-dipole interactions. The studied system, liquid water stressed by an electric field, exhibits collective oscillations brought about through spontaneous breakdown of symmetry. Raman spectroscopy identifies the primary excitation of the emergent phase as transverse optically active phonon-like sidebands that appear on the hydrogen bonded asymmetric stretch mode. The phase transition is observed throughout the entire volume of liquid. The system also exhibits a self-similarity relation between the scattered Raman intensity and the electric field strength which further supports the conclusion that collective behavior persists against thermal disruption. The experimental findings are discussed in terms of a quantum field theory for macroscopic quantum systems.

Hydration-Shell Vibrational Spectroscopy

Hydration-shell vibrational spectroscopy provides an experimental window into solute-induced water structure changes that mediate aqueous folding, binding, and self-assembly. Decomposition of measured Raman and infrared (IR) spectra of aqueous solutions using multivariate curve resolution (MCR) and related methods may be used to obtain solute-correlated spectra revealing solute-induced perturbations of water structure, such as changes in water hydrogen-bond strength, tetrahedral order, and the presence of dangling (non-hydrogen-bonded) OH groups. More generally, vibrational-MCR may be applied to both aqueous and nonaqueous solutions, including multicomponent mixtures, to quantify solvent-mediated interactions between oily, polar, and ionic solutes, in both dilute and crowded fluids. Combining vibrational-MCR with emerging theoretical modeling strategies promises synergetic advances in the predictive understanding of multiscale self-assembly processes of both biological and technological interest.

Probing hydrogen bond strength in liquid water by resonant inelastic X-ray scattering

Local probes of the electronic ground state are essential for understanding hydrogen bonding in aqueous environments. When tuned to the dissociative core-excited state at the O1s pre-edge of water, resonant inelastic X-ray scattering back to the electronic ground state exhibits a long vibrational progression due to ultrafast nuclear dynamics. We show how the coherent evolution of the OH bonds around the core-excited oxygen provides access to high vibrational levels in liquid water. The OH bonds stretch into the long-range part of the potential energy curve, which makes the X-ray probe more sensitive than infra-red spectroscopy to the local environment. We exploit this property to effectively probe hydrogen bond strength via the distribution of intramolecular OH potentials derived from measurements. In contrast, the dynamical splitting in the spectral feature of the lowest valence-excited state arises from the short-range part of the OH potential curve and is rather insensitive to hydrogen bonding.

The Interplay of Structure and Dynamics in the Raman Spectrum of Liquid Water over the Full Frequency and Temperature Range

While many vibrational Raman spectroscopy studies of liquid water have investigated the temperature dependence of the high-frequency O-H stretching region, few have analyzed the changes in the Raman spectrum as a function of temperature over the entire spectral range. Here, we obtain the Raman spectra of water from its melting to boiling point, both experimentally and from simulations using an ab initio-trained machine learning potential. We use these to assign the Raman bands and show that the entire spectrum can be well described as a combination of two temperature-independent spectra. We then assess which spectral regions exhibit strong dependence on the local tetrahedral order in the liquid. Further, this work demonstrates that changes in this structural parameter can be used to elucidate the temperature dependence of the Raman spectrum of liquid water and provides a guide to the Raman features that signal water ordering in more complex aqueous systems.

CO2 Hydration Shell Structure and Transformation

The hydration-shell of CO₂ is characterized using Raman multivariate curve resolution (Raman-MCR) spectroscopy combined with ab initio molecular dynamics (AIMD) vibrational density of states simulations, to validate our assignment of the experimentally observed high-frequency OH band to a weak hydrogen bond between water and CO₂. Our results reveal that while the hydration-shell of CO₂ is highly tetrahedral, it is also occasionally disrupted by the presence of entropically stabilized defects associated with the CO₂-water hydrogen bond. Moreover, we find that the hydration-shell of CO₂ undergoes a temperature-dependent structural transformation to a highly disordered (less tetrahedral) structure, reminiscent of the transformation that takes place at higher temperatures around much larger oily molecules. The biological significance of the CO₂ hydration shell structural transformation is suggested by the fact that it takes place near physiological temperatures.

Structural and electrostatic effects at the surfaces of size- and charge-selected aqueous nanodrops

The effects of ion charge, polarity and size on the surface morphology of size-selected aqueous nanodrops containing a single ion and up to 550 water molecules are investigated with infrared photodissociation (IRPD) spectroscopy and theory.

The effects of ion charge, polarity and size on the surface morphology of size-selected aqueous nanodrops containing a single ion and up to 550 water molecules are investigated with infrared photodissociation (IRPD) spectroscopy and theory. IRPD spectra of M(H₂O)_n where M = La³⁺, Ca²⁺, Na⁺, Li⁺, I^–, SO₄^2– and supporting molecular dynamics simulations indicate that strong interactions between multiply charged ions and water molecules can disrupt optimal hydrogen bonding (H-bonding) at the nanodrop surface. The IRPD spectra also reveal that “free” OH stretching frequencies of surface-bound water molecules are highly sensitive to the ion's identity and the OH bond's local H-bond environment. The measured frequency shifts are qualitatively reproduced by a computationally inexpensive point-charge model that shows the frequency shifts are consistent with a Stark shift from the ion's electric field. For multiply charged cations, pronounced Stark shifting is observed for clusters containing ∼100 or fewer water molecules. This is attributed to ion-induced solvent patterning that extends to the nanodrop surface, and serves as a spectroscopic signature for a cation's ability to influence the H-bond network of water located remotely from the ion. The Stark shifts measured for the larger nanodrops are extrapolated to infinite dilution to obtain the free OH stretching frequency of a surface-bound water molecule at the bulk air–water interface (3696.5–3701.0 cm^–1), well within the relatively wide range of values obtained from SFG measurements. These cluster measurements also indicate that surface curvature effects can influence the free OH stretching frequency, and that even nanodrops without an ion have a surface potential that depends on cluster size.

Dimension of discrete variable representation for mixed quantum/classical computation of three lowest vibrational states of OH stretching in liquid water

Three low-lying vibrational states of molecular systems are responsible for the signals of linear and third-order nonlinear vibrational spectroscopies. Theoretical studies based on mixed quantum/classical calculations provide a powerful way to analyze those experiments. A statistically meaningful result can be obtained from the calculations by solving the vibrational Schrödinger equation over many numbers of molecular configurations. The discrete variable representation (DVR) method is a useful technique to calculate vibrational eigenstates subject to an arbitrary anharmonic potential surface. Considering the large number of molecular configurations over which the DVR calculations are repeated, the calculations are desired to be optimized in balance between the cost and accuracy. We determine a dimension of the DVR method which appears to be optimum for the calculations of the three states of molecular vibrations with anharmonic strengths often found in realistic molecular systems. We apply the numerical technique to calculate the local OH stretching frequencies of liquid water, which are well known to be widely distributed due to the inhomogeneity in molecular configuration, and found that the frequencies of the 0-1 and 1-2 transitions are highly correlated. An empirical relation between the two frequencies is suggested and compared with the experimental data of nonlinear IR spectroscopies.

Non-equilibrium thermodynamics and collective vibrational modes of liquid water in an inhomogeneous electric field

In this experiment liquid water is subject to an inhomogeneous electric field (∇(2)Ea≈ 10(10) V m(2)) using a high voltage (20 kV) point-plane electrode system. Using interferometry it was found that the application of a strong electric field gradient to water generates local changes in the refractive index of the liquid, polarizes the surface and creates a downward moving electro-convective jet. A maximum temperature difference of 1 °C is measured in the immediate vicinity of the point electrode. Raman spectroscopy performed on water reveals an enhancement of the vibrational collective modes (3250 cm(-1)) as well as an increase in the local mode (3490 cm(-1)) energy. This bimodal enhancement indicates that the spectral changes are not due to temperature changes. The intense field gradient thus establishes an excited subpopulation of vibrational oscillators far from thermal equilibrium. Delocalization of the collective vibrational mode spatially expands this excited population beyond the microscale. Hindered rotational freedom due to electric field pinning of molecular dipoles retards the heat flow and generates a chemical potential gradient. These changes are responsible for the observed changes in the refractive index and temperature. It is demonstrated that polar liquids can thus support local non-equilibrium thermodynamic transient states critical to biochemical and environmental processes.

Vibrational Quantum Decoherence in Liquid Water

Traditional descriptions of vibrational energy transfer consider a quantum oscillator interacting with a classical environment. However, a major limitation of this simplified description is the neglect of quantum decoherence induced by the different interactions between two distinct quantum states and their environment, which can strongly affect the predicted energy-transfer rate and vibrational spectra. Here, we use quantum-classical molecular dynamics simulations to determine the vibrational quantum decoherence time for an OH stretch vibration in liquid heavy water. We show that coherence is lost on a sub-100 fs time scale due to the different responses of the first shell neighbors to the ground and excited OH vibrational states. This ultrafast decoherence induces a strong homogeneous contribution to the linear infrared spectrum and suggests that resonant vibrational energy transfer in H2O may be more incoherent than previously thought.

Hydrogen-bonding and vibrational coupling of water in a hydrophobic hydration shell as observed by Raman-MCR and isotopic dilution spectroscopy

Raman multivariate curve resolution (Raman-MCR) spectroscopy reveals the perturbation of vibrational coupling of water in a hydrophobic hydration shell.

Femtosecond coherent anti-Stokes Raman scattering spectroscopy of hydrogen bonded structure in water and aqueous solutions

Femtosecond coherent anti-Stokes Raman scattering (fsCARS) spectroscopy, together with perturbation theory based numerical calculation, is employed to study OH stretching (υOH) of pure water and aqueous lithium chloride solutions. Vibrational OH stretching (υOH) modes of aqueous solutions are Raman-excited by a pair of ultrashort, femtosecond laser pulses, and then probed through inelastic scattering of a third, time-delayed laser field. In order to overcome limited spectral resolution of fsCARS, numerical evaluation of the CARS signal through vibrational wave packet propagation was employed in order to confirm the position of distinctive OH stretching mode that is complicated by intramolecular and intermolecular vibrational coupling. Moreover, in order to come to a microscopic description of the observed CARS spectra for aqueous solutions, we have performed molecular dynamics simulations of aqueous lithium chloride solutions with varying concentrations at ambient conditions. To this end we have analyzed the equilibrium distributions of hydrogen bonds in the first solvation shells of the ions as well as in bulk water and also computed the average number of hydrogen bonds per water molecule. According to our experimental and theoretical results on time evolution of Raman OH stretching band of water, it can be inferred that the dissolved ions mainly influence hydrogen bond strength and structure of water molecules in the first hydration shell, the addition of lithium chloride primarily breaks the tetrahedral hydrogen bonding, promotes formation of the donor hydrogen bonding in water, and slightly increases the amount of free OH bonds.

Determination of Spectroscopic Band Shapes by Second Derivatives, Part II: Infrared Spectra of Liquid Light and Heavy Water

Second derivative and band simulation techniques are used in a synergetic relationship to identify components in the infrared (IR) spectra of liquid light and heavy water. Nine Gaussian components are retrieved in massive OH and OD stretch absorption. In this context, ν1 and ν3 are the principal components along with satellites derived from harmonic and combination bands. The Raman spectrum of light water matches the IR components with intensity variations as expected.

Kubo's Line Shape Function for a Linear-Quadratic Chromophore-Solvent Coupling

An exact, closed-form solution is obtained for the line shape function of an optical transition with the transition frequency depending linearly plus quadratically on a Gaussian coordinate of the thermal bath. The dynamical modulation of the line shape involves two parameters corresponding to the linear and quadratic components of the transition frequency. The increase of the second component results in a non-Gaussian line shape that splits into two Lorenzian lines in the limit of fast modulation.

Local orientational order in liquids revealed by resonant vibrational energy transfer

We demonstrate that local orientational ordering in a liquid can be observed in the decay of the vibrational anisotropy caused by resonant transfer of vibrational excitations between its constituent molecules. We show that the functional form of this decay is determined by the (distribution of) angles between the vibrating bonds of the molecules between which energy transfer occurs, and that the initial drop in the decay reflects the average angle between nearest neighbors. We use this effect to observe the difference in local orientational ordering in the two hydrogen-bonded liquids ethanol and N-methylacetamide.

Non-Gaussian lineshapes and dynamics of time-resolved linear and nonlinear (correlation) spectra

Signatures of nonlinear and non-Gaussian dynamics in time-resolved linear and nonlinear (correlation) 2D spectra are analyzed in a model considering a linear plus quadratic dependence of the spectroscopic transition frequency on a Gaussian nuclear coordinate of the thermal bath (quadratic coupling). This new model is contrasted to the commonly assumed linear dependence of the transition frequency on the medium nuclear coordinates (linear coupling). The linear coupling model predicts equality between the Stokes shift and equilibrium correlation functions of the transition frequency and time-independent spectral width. Both predictions are often violated, and we are asking here the question of whether a nonlinear solvent response and/or non-Gaussian dynamics are required to explain these observations. We find that correlation functions of spectroscopic observables calculated in the quadratic coupling model depend on the chromophore's electronic state and the spectral width gains time dependence, all in violation of the predictions of the linear coupling models. Lineshape functions of 2D spectra are derived assuming Ornstein-Uhlenbeck dynamics of the bath nuclear modes. The model predicts asymmetry of 2D correlation plots and bending of the center line. The latter is often used to extract two-point correlation functions from 2D spectra. The dynamics of the transition frequency are non-Gaussian. However, the effect of non-Gaussian dynamics is limited to the third-order (skewness) time correlation function, without affecting the time correlation functions of higher order. The theory is tested against molecular dynamics simulations of a model polar-polarizable chromophore dissolved in a force field water.

Effects of excluded volume and correlated molecular orientations on Förster resonance energy transfer in liquid water

Förster theory for the survival probability of excited chromophores is generalized to include the effects of excluded volume and orientation correlation in the molecular distribution. An analytical expression for survival probability was derived and written in terms of a few simple elementary functions. Because of the excluded volume, the survival probability exhibits exponential decay at early times and stretched exponential decay at later times. Experimental schemes to determine the size of the molecular excluded volume are suggested. With the present generalization of theory, we analyzed vibrational resonance energy transfer kinetics in neat water. Excluded volume effects prove to be important and slow down the kinetics at early times. The majority of intermolecular resonance energy transfer was found to occur with exponential kinetics, as opposed to the stretched exponential behavior predicted by Förster theory. Quantum yields of intra-molecular vibrational relaxation, intra-, and intermolecular energy transfer were calculated to be 0.413, 0.167, and 0.420, respectively.