Fermi resonance in OH-stretch vibrational spectroscopy of liquid water and the water hexamer

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.

H-Glass Supported Hybrid Gold Nano-Islands for Visible-Light-Driven Hydrogen Evolution

Flat panel reactors, coated with photocatalytic materials, offer a sustainable approach for the commercial production of hydrogen (H2) with zero carbon footprint. Despite this, achieving high solar-to-hydrogen (STH) conversion efficiency with these reactors is still a significant challenge due to the low utilization efficiency of solar light and rapid charge recombination. Herein, hybrid gold nano-islands (HGNIs) are developed on transparent glass support to improve the STH efficiency. Plasmonic HGNIs are grown on an in-house developed active glass sheet composed of sodium aluminum phosphosilicate oxide glass (H-glass) using the thermal dewetting method at 550 °C under an ambient atmosphere. HGNIs with various oxidation states (Au0, Au+, and Au-) and multiple interfaces are obtained due to the diffusion of the elements from the glass structure, which also facilitates the lifetime of the hot electron to be ≈2.94 ps. H-glass-supported HGNIs demonstrate significant STH conversion efficiency of 0.6%, without any sacrificial agents, via water dissociation. This study unveils the specific role of H-glass-supported HGNIs in facilitating light-driven chemical conversions, offering new avenues for the development of high-performance photocatalysts in various chemical conversion reactions for large-scale commercial applications.

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.

Path Integral Simulations of Condensed-Phase Vibrational Spectroscopy

Recent theoretical and algorithmic developments have improved the accuracy with which path integral dynamics methods can include nuclear quantum effects in simulations of condensed-phase vibrational spectra. Such methods are now understood to be approximations to the delocalized classical Matsubara dynamics of smooth Feynman paths, which dominate the dynamics of systems such as liquid water at room temperature. Focusing mainly on simulations of liquid water and hexagonal ice, we explain how the recently developed quasicentroid molecular dynamics (QCMD), fast-QCMD, and temperature-elevated path integral coarse-graining simulations (Te PIGS) methods generate classical dynamics on potentials of mean force obtained by averaging over quantum thermal fluctuations. These new methods give very close agreement with one another, and the Te PIGS method has recently yielded excellent agreement with experimentally measured vibrational spectra for liquid water, ice, and the liquid-air interface. We also discuss the limitations of such methods.

Monitoring water harvesting in metal-organic frameworks, one water molecule at a time

Metal-organic frameworks (MOFs) have gained prominence as potential materials for atmospheric water harvesting, a vital solution for arid regions and areas experiencing severe water shortages. However, the molecular factors influencing the performance of MOFs in capturing water from the air remain elusive. Among all MOFs, Ni2X2BTDD (X = F, Cl, Br) stands out as a promising water harvester due to its ability to adsorb substantial amounts of water at low relative humidity (RH). Here, we use advanced molecular dynamics simulations carried out with the state-of-the-art MB-pol data-driven many-body potential to monitor water adsorption in the three Ni2X2BTDD variants as a function of RH. Our simulations reveal that the type of halide atom in the three Ni2X2BTDD frameworks significantly influences the corresponding molecular mechanisms of water adsorption: while water molecules form strong hydrogen bonds with the fluoride atoms in Ni2F2BTDD, they tend to form hydrogen bonds with the nitrogen atoms of the triazolate linkers in Ni2Cl2BTDD and Ni2Br2BTDD. Importantly, the large size of the bromide atoms reduces the void volume in the Ni2Br2BTDD pores, which enable water molecules to initiate an extended hydrogen-bond network at lower RH. These findings not only underscore the prospect for precisely tuning structural and chemical modifications of the frameworks to optimize their interaction with water, but also highlight the predictive power of simulations with the MB-pol data-driven many-body potential. By providing a realistic description of water under different thermodynamic conditions and environments, these simulations yield unique, molecular-level insights that can guide the design and optimization of energy-efficient water harvesting materials.

The strength of the OH-bend/OH-stretch Fermi resonance in small water clusters

A novel Raman jet-spectrometer is used to study the Fermi resonance between the OH bending overtone and OH stretching fundamental in small cyclic water clusters (H2O)n with n = 3, 4, 5. The new setup features a recirculating vacuum system which reduces the gas consumption by 2 to 3 orders of magnitude and enables long-term measurements of very weak Raman signals. Raman spectra measured from highly diluted expansions with unprecedented signal-to-noise ratio are presented and cluster-specific intensity ratios and effective coupling constants are derived using Markov-Chain Monte-Carlo methods, yielding a high probability for an almost "perfect" resonance for the tetramer and pentamer, i.e. a close frequency match of bend overtone and stretch fundamental with intensity ratios close to 1, but a larger coupling constant for the trimer, with best estimates close to W5 ≲ 50 cm-1 < W4 ≲ 60 cm-1 < W3 ≈ 65 cm-1.

Harmonic and anharmonic studies on THz spectra of two vanillin polymorphs

Polymorphism commonly exists in organic molecular crystals. The fingerprint features in low-frequency vibrational range are important information reflecting different intermolecular interactions of polymorphs. Interpreting these features is very helpful to understand vibrational property of polymorphs and reveal the thermodynamic stability. In this work, the low-frequency vibrations of form I and II of vanillin are investigated using terahertz time-domain spectroscopy. Static DFT calculation and ab initio molecular dynamics (AIMD) are employed to interpret their low-frequency vibrations of both forms in harmonic and anharmonic ways, respectively. Their low-frequency vibration characteristics in harmonic calculations are discussed, and anharmonic mode couplings between OH bond stretch and the stretching and bending motion of hydrogen bonds are uncovered. Moreover, the thermodynamic energies including electronic potential energy and vibrational/kinetic energy arising from nuclear motions are calculated. The result reveals that the stability order of the two forms is mainly dependent on their electric potential energy difference.

Theoretical and experimental OD-stretch vibrational spectroscopy of heavy water

In view of the current situation in which the OD-stretch vibrational spectra have been scarcely computed with non-polarizable rigid D2O models, we investigate the IR and Raman spectra of D2O by using a newly-reported model TIP4P/2005-HW. From the comparison between the calculations and experimental data, we find the excellent performance of TIP4P/2005-HW for vibrational spectroscopy of D2O in the same manner as TIP4P/2005 for H2O, although one may still conveniently employ an alternative method that regards OH as putative OD to calculate the OD-stretch spectra with similar quality from TIP4P/2005 trajectories. We also demonstrate that the appropriate setting for the spectral simulation of D2O under the time-averaging approximation reflects the slower dynamics (i.e., slower motion of translation and rotation due to the heavier mass and stronger hydrogen bond) of D2O than H2O. Moreover, we show from the theoretical calculations that the established interpretation of the OH-stretch spectra of H2O is finely applicable to the OD-stretch of D2O.

Correlations between the Structures and Spectra of Protonated Water Clusters

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. We find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation-dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation-dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

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.

Microhydration of Tertiary Amines: Robust Resonances in Red-Shifted Water

Tertiary amines are strong hydrogen bond acceptors. When a water molecule donates one of the OH groups, its in-phase stretching vibration wavenumber is decreased to such an extent that it comes close to the water bending overtone. This gives rise to anharmonic phenomena such as classical Fermi resonance, resonance with multiple-quantum dark states, or combination transitions with low-frequency intermolecular modes. These effects, which contribute to the spectral breadth of room-temperature hydrogen-bonded amine complexes, are disentangled by Fourier transform infrared spectroscopy in pulsed supersonic slit jet expansions. Monohydrates of the amines quinuclidine, N-methylpyrrolidine, N-methylpiperidine, and dimethylcyclohexylamine exhibit systematic mode coupling signatures. These suggest relatively fast energy flow out of the excited OH stretching fundamental into intra- and intermolecular degrees of freedom of the hydrogen-bonded water molecule. Trimeric complexes are spectroscopically separated from the amine monohydrates.

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.

Exploring the Origins of the Intensity of the OH Stretch-HOH Bend Combination Band in Water

While the intensity of the OH stretching fundamental transition is strongly correlated to hydrogen-bond strength, the intensity of the corresponding transition to the state with one quantum of excitation in both the OH stretching and HOH bending vibrations in the same water molecule shows a much weaker sensitivity to the hydrogen-bonding environment. The origins of this difference are explored through analyses of the contributions of terms in the expansion of the dipole moment to the calculated intensity. It is found that the leading contribution to the stretch-bend intensity involves the second derivative of the dipole moment with respect to the OH bond length and HOH angle. While this is not surprising, the insensitivity of this derivative to the hydrogen-bonding environment is unexpected. Possible contributions of mode mixing are also explored. While mode mixing leads to splittings of the energies of nearly degenerate excited states, it does not result in significant changes in the sum of the intensities of these transitions. Analysis of changes in the partial charges on the hydrogen atoms upon displacement of the HOH angles shows that these charges generally increase with increasing HOH angle. This effect is partially canceled by a decrease in the charge of the hydrogen atom when a hydrogen bond is broken. The extent of this cancellation increases with the hydrogen bond strength, which is reflected in the observed insensitivity of the intensity of the stretch-bend transition to hydrogen-bond strength.

Combinational Vibration Modes in H2O/HDO/D2O Mixtures Detected Thanks to the Superior Sensitivity of Femtosecond Stimulated Raman Scattering

Overtones and combinational modes frequently play essential roles in ultrafast vibrational energy relaxation in liquid water. However, these modes are very weak and often overlap with fundamental modes, particularly in isotopologues mixtures. We measured VV and HV Raman spectra of H₂O and D₂O mixtures with femtosecond stimulated Raman scattering (FSRS) and compared the results with calculated spectra. Specifically, we observed the mode at around 1850 cm^-1 and assigned it to H-O-D bend + rocking libration. Second, we found that the H-O-D bend overtone band and the OD stretch + rocking libration combination band contribute to the band located between 2850 and 3050 cm^-1. Furthermore, we assigned the broad band located between 4000 and 4200 cm^-1 to be composed of combinational modes of high-frequency OH stretching modes with predominantly twisting and rocking librations. These results should help in a proper interpretation of Raman spectra of aqueous systems as well as in the identification of vibrational relaxation pathways in isotopically diluted water.

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

Is Unified Understanding of Vibrational Coupling of Water Possible? Hyper-Raman Measurement and Machine Learning Spectra

The impact of the vibrational coupling of the OH stretch mode on the spectra differs significantly between IR and Raman spectra of water. Unified understanding of the vibrational couplings is not yet achieved. By using a different class of vibrational spectroscopy, hyper-Raman (HR) spectroscopy, together with machine-learning-assisted HR spectra calculation, we examine the impact of the vibrational couplings of water through the comparison of isotopically diluted H₂O and pure H₂O. We found that the isotopic dilution reduces the HR bandwidths, but the impact of the vibrational coupling is smaller than in the IR and parallel-polarized Raman. Machine learning HR spectra indicate that the intermolecular coupling plays a major role in broadening the bandwidth, while the intramolecular coupling is negligibly small, which is consistent with the IR and Raman spectra. Our result clearly demonstrates a limited impact of the intramolecular vibration, independent of the selection rules of vibrational spectroscopies.

Design of an Electrostatic Frequency Map for the NH Stretch of the Protein Backbone and Application to Chiral Sum Frequency Generation Spectroscopy

We develop an electrostatic map for the vibrational NH stretch (amide A) of the protein backbone with a focus on vibrational chiral sum frequency generation spectroscopy (chiral SFG). Chiral SFG has been used to characterize protein secondary structure at interfaces using the NH stretch and to investigate chiral water superstructures around proteins using the OH stretch. Interpretation of spectra has been complicated because the NH stretch and OH stretch overlap spectrally. Although an electrostatic map for water OH developed by Skinner and co-workers was used previously to calculate the chiral SFG response of water structures around proteins, a map for protein NH that is directly responsive to biological complexity has yet to be developed. Here, we develop such a map, linking the local electric field to vibrational frequencies and transition dipoles. We apply the map to two protein systems and achieve much better agreement with experiment than was possible in our previous studies. We show that couplings between NH and OH vibrations are crucial to the line shape, which informs the interpretation of chiral SFG spectra, and that the chiral NH stretch response is sensitive to small differences in structure. This work increases the utility of the NH stretch in biomolecular spectroscopy.

Hydrogen-Bond Restructuring of Water-in-Salt Electrolyte Confined in Ti3C2Tx MXene Monitored by Operando Infrared Spectroscopy

Highly concentrated water-in-salt aqueous electrolytes exhibit a wider potential window compared to conventional, dilute aqueous electrolytes. Coupled with MXenes, a family of two-dimensional transition metal carbides and nitrides with impressive charge storage capabilities, water-in-salt electrolytes present a potential candidate to replace flammable and toxic organic solvents in electrochemical energy storage devices. A new charge storage mechanism was recently discovered during electrochemical cycling of Ti₃C₂T_x MXene electrodes in lithium-based water-in-salt electrolytes, attributed to intercalation and deintercalation of solvated Li⁺ ions at anodic potentials. Nevertheless, direct evidence of the state of Li⁺ solvation during cycling is still missing. Here, we investigate the hydrogen bonding of water intercalated between MXene layers during electrochemical cycling in a water-in-salt electrolyte with operando infrared spectroscopy. The hydrogen-bonding state of the confined water was found to change significantly as a function of potential and the concentration of Li⁺ ions in the interlayer space. This study provides fundamentally new insights into the electrolyte structural changes while intercalating Li⁺ in the MXene interlayer space.

On the Fresnel factor correction of sum-frequency generation spectra of interfacial water

Insights into the microscopic structure of aqueous interfaces are essential for understanding the chemical and physical processes on the water surface, including chemical synthesis, atmospheric chemistry, and events in biomolecular systems. These aqueous interfaces have been probed by heterodyne-detected sum-frequency generation (HD-SFG) spectroscopy. To obtain the molecular response from the measured HD-SFG spectra, one needs to correct the measured ssp spectra for local electromagnetic field effects at the interface due to a spatially varying dielectric function. This so-called Fresnel factor correction can change the inferred response substantially, and different ways of performing this correction lead to different conclusions about the interfacial water response. Here, we compare the simulated and experimental spectra at the air/water interface. We use three previously developed models to compare the experiment with theory: an advanced approach taking into account the detailed inhomogeneous interfacial dielectric profile and the Lorentz and slab models to approximate the interfacial dielectric function. Using the advanced model, we obtain an excellent quantitative agreement between theory and experiment, in both spectral shape and amplitude. Remarkably, we find that for the Fresnel factor correction of the ssp spectra, the Lorentz model for the interfacial dielectric function is equally accurate in the hydrogen (H)-bonded region of the response, while the slab model underestimates this response significantly. The Lorentz model, thus, provides a straightforward method to obtain the molecular response from the measured spectra of aqueous interfaces in the H-bonded region.

Hyper-mobile Water and Raman 2900 cm-1 Peak Band of Water Observed around Backbone Phosphates of Double Stranded DNA by High-Resolution Spectroscopies and MD Structural Feature Analysis of Water

High-resolution measurements of microwave dielectric relaxation and Raman spectroscopies for waters in double-stranded (ds) 10-mer DNA solution revealed the presence of hyper-mobile water (HMW) and a marked OH stretching band appearing in the range from 2500 to 3100 cm^-1, here called the LA band, at the low wavenumber tail of the major OH stretching band of water. Quantitation of the Raman scattering intensity for ds 10-mer DNA in phosphate or tris(hydroxymethyl)aminomethane (TRIS) buffers showed that the LA band was formed by 2000-3000 water molecules per ds 10-mer DNA, indicating collective OH stretching vibrations of water molecules around the backbone phosphate oxygen atoms. The LA band intensity of ds 10-mer DNA in 10 mM TRIS increased and decreased by 30% with the addition of 2 mM MgCl₂ and 2 mM CaCl₂, respectively. The LA band intensity and the effect of adding Mg(II) or Ca(II) ions to the band intensity were maintained in the presence of 0.14 M KCl; however, the changes induced by the divalent cations were reduced by half. Molecular dynamics calculations of water molecules around the backbone phosphate groups of ds 10-mer DNA indicate the presence of high-density water and broad regions of fluctuating water density, suggesting that they correspond to HMW and the LA band, respectively.

Evolution of Chemical Bonding and Crystalline Swelling-Shrinkage of Montmorillonite upon Temperature Changes Probed by in Situ Fourier Transform Infrared Spectroscopy and X-ray Diffraction

Clay minerals are distributed in Earth's crust and troposphere and in Martian crust where temperature varies. Understanding the changes of chemical bonding and crystalline swelling-shrinkage of montmorillonite (Mnt) upon temperature changes is fundamental for studying its surface reactivity and interaction in specific surroundings. However, such an issue remains poorly understood. Here, in situ high- and low-temperature Fourier transform infrared (HT- and LT-FTIR) spectroscopy and X-ray diffraction (HT- and LT-XRD) were performed to study the evolution of chemical bonding and crystalline swelling-shrinkage of sodium-montmorillonite (NaMnt) upon temperature changes. The FTIR results show that the hydroxyl content in NaMnt decreased when the temperature increased from 20 to 700 °C, while it is independent of temperature from 0 to -150 °C. The formation of hydroxyls at the "broken" layer edges of NaMnt is related to adsorbed water molecules on the surfaces, and its content increased when the particle size of the NaMnt decreased. The water molecules in the interlayer space of NaMnt could bond to the tetrahedral sheet of NaMnt through Si₂O-H₂O bonds. HT- and LT-XRD results reveal that all of those water molecules in NaMnt were removed after heating to 100 °C in the heating chamber. The NaMnt was transformed from a state of monolayer interlayer water molecules at 20 °C to a dehydrated state at 100 °C, and then to a dehydroxylated state at 700 °C. Accordingly, the basal spacings of NaMnt were changed from 1.27 to 0.97 nm and then to 0.96 nm, respectively. When NaMnt was cooled from 20 to -268 °C, a crystalline swelling of NaMnt occurred with an increase of 0.03 nm of basal spacing. This work demonstrates that high/low temperature has a remarkable effect on the hydroxyls and the water molecules in NaMnt, which in turn affects its swelling-shrinkage performance. These findings provide some in-depth understanding of the changes of chemical bonding and crystalline swelling-shrinkage of montmorillonite upon temperature changes and the reasons behind these, which might be helpful for the design of engineering Mnt in high-/low-temperature applications.

Detecting the First Hydration Shell Structure around Biomolecules at Interfaces

Understanding the role of water in biological processes remains a central challenge in the life sciences. Water structures in hydration shells of biomolecules are difficult to study in situ due to overwhelming background from aqueous environments. Biological interfaces introduce additional complexity because biomolecular hydration differs at interfaces compared to bulk solution. Here, we perform experimental and computational studies of chiral sum frequency generation (chiral SFG) spectroscopy to probe chirality transfer from a protein to the surrounding water molecules. This work reveals that chiral SFG probes the first hydration shell around the protein almost exclusively. We explain the selectivity to the first hydration shell in terms of the asymmetry induced by the protein structure and specific protein-water hydrogen-bonding interactions. This work establishes chiral SFG as a powerful technique for studying hydration shell structures around biomolecules at interfaces, presenting new possibilities to address grand research challenges in biology, including the molecular origins of life.

Experimental and Theoretical Heterodyne-Detected Sum Frequency Generation Spectroscopy of Isotopically Pure and Diluted Water Surfaces

The χ⁽²⁾ (second-order nonlinear optical susceptibility) spectrum of the water surface has been a matter of debate for a few decades. Here, we report that we experimentally measured the isotopic dilution dependence of the χ⁽²⁾ spectrum and theoretically reproduced it by employing the quantum/classical mixed approach with a new idea to subtract an artifact. The present theoretical framework allows for clarifying the effects of the intramolecular, intermolecular, and Fermi resonance couplings on the OH-stretch vibrational spectra of water at the surface as well as in the bulk.

Character of the OH Bend-Stretch Combination Band in the Vibrational Spectra of the "Magic" Number H3O+(H2O)20 and D3O+(D2O)20 Cluster Ions

The fundamental transitions that contribute to the diffuse OH stretching spectrum of water are known to increase in width and intensity with increasing red shift from the free OH frequency. In contrast, the profile of the higher-energy combination band involving the OH stretching and the intramolecular HOH bending modes displays a qualitatively different spectral shape with a much faster falloff on the lower-energy side. We elucidate the molecular origin of this difference by analyzing the shapes of the combination bands in the IR spectra of cryogenically cooled H₃O⁺(H₂O)₂₀ and D₃O⁺(D₂O)₂₀ clusters. The difference in the shapes of the bands is traced to differences in the dependence of their transition dipole matrix elements on the hydrogen-bonding environment. The fact that individual transitions across the combination band envelope have similar intensities makes it a useful way to determine the participation of various sites in extended H-bonding networks.

Rapid Water Dynamics Structures the OH-Stretching Spectra of Solitary Water in Ionic Liquids and Dipolar Solvents

In a recent study [ J. Phys. Chem. B 126, 4584 (2022)] we used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands common to all high-polarity solvents were assigned to water in several solvation states. In the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [ J. Chem. Phys. 128, 224511 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra also confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared to the inhomogeneous limit, and that this narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.

Pressure-induced Fermi resonance between fundamental modes in phthalic anhydride

In situhigh-pressure Raman spectra of phthalic anhydride (PA) have been measured up to 16 GPa through diamond anvil cell technique. The results show that all the Raman bands are blue-shifted with the increase of pressure, accompanied by appearance of some new bands. A Fermi resonance phenomenon of the two Raman fundamental modes of PA at 773 cm^-1and 801 cm^-1is proposed at pressures above 6.6 GPa, where a possible first-order phase transition occurs. The pressure-induced changes of Fermi resonance parameters, e.g., intensity ratio, coupling coefficient and frequency gap of unperturbed transition, are discussed.

Time-Dependent Friction Effects on Vibrational Infrared Frequencies and Line Shapes of Liquid Water

From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged nonlinear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits not only adiabatic contributions at and below the vibrational time scales but also much slower nonadiabatic contributions, reflecting homogeneous and inhomogeneous line broadening mechanisms, respectively. Intermolecular interactions in liquid water soften both stretch and bend potentials compared to the gas phase, which by itself would lead to a red-shift of the corresponding vibrational bands. In contrast, nonadiabatic friction contributions cause a spectral blue shift. For the stretch mode, the potential effect dominates, and thus, a significant red shift when going from gas to the liquid phase results. For the bend mode, potential and nonadiabatic friction effects are of comparable magnitude, so that a slight blue shift results, in agreement with well-known but puzzling experimental findings. The observed line broadening is shown to be roughly equally caused by adiabatic and nonadiabatic friction contributions for both the stretch and bend modes in liquid water. Thus, the quantitative analysis of the time-dependent friction that acts on vibrational modes in liquids advances the understanding of infrared vibrational frequencies and line shapes.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

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.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H⁺(H₂O)₄ and H⁺(H₂O)₂₁ agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H⁺(H₂O)₅, which can adopt several low-energy conformations.

Structural Characterization of Protonated Water Clusters Confined in HZSM-5 Zeolites

A molecular description of the structure and behavior of water confined in aluminosilicate zeolite pores is a crucial component for understanding zeolite acid chemistry under hydrous conditions. In this study, we use a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and ab initio molecular dynamics (AIMD) to study H₂O confined in the pores of highly hydrated zeolite HZSM-5 (∼13 and ∼6 equivalents of H₂O per Al atom). The 2D IR spectrum reveals correlations between the vibrations of both terminal and H-bonded O-H groups and the continuum absorption of the excess proton. These data are used to characterize the hydrogen-bonding network within the cluster by quantifying single-, double-, and non-hydrogen-bond donor water molecules. These results are found to be in good agreement with the statistics calculated from an AIMD simulation of an H⁺(H₂O)₈ cluster in HZSM-5. Furthermore, IR spectral assignments to local O-H environments are validated with DFT calculations on clusters drawn from AIMD simulations. The simulations reveal that the excess charge is detached from the zeolite and resides near the more highly coordinated water molecules in the cluster. When they are taken together, these results unambiguously assign the complex IR spectrum of highly hydrated HZSM-5, providing quantitative information on the molecular environments and hydrogen-bonding topology of protonated water clusters under extreme confinement.

Ab initio molecular dynamics study on energy relaxation path of hydrogen-bonded OH vibration in bulk water

The vibrational energy relaxation paths of hydrogen-bonded (H-bonded) OH excited in pure water and in isotopically diluted (deuterated) water are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The present study extends the previous NE-AIMD simulation for the energy relaxation of an excited free OH vibration at an air/water interface [T. Ishiyama, J. Chem. Phys. 154, 104708 (2021)] to the energy relaxation of an excited H-bonded OH vibration in bulk water. The present simulation shows that the excited OH vibration in pure water dissipates its energy on a timescale of 0.1 ps, whereas that in deuterated water relaxes on a timescale of 0.7 ps, consistent with the experimental observations. To decompose these relaxation energies into the components due to intramolecular and intermolecular couplings, constraints are introduced on the vibrational modes except for the target path in the NE-AIMD simulation. In the case of pure water, 80% of the total relaxation is attributed to the pathway due to the resonant intermolecular OH⋯OH stretch coupling, and the remaining 17% and 3% are attributed to intramolecular couplings with the bend overtone and with the conjugate OH stretch, respectively. This result strongly supports a significant role for the Förster transfer mechanism of pure water due to the intermolecular dipole-dipole interactions. In the case of deuterated water, on the other hand, 36% of the total relaxation is due to the intermolecular stretch coupling, and all the remaining 64% arises from coupling with the intramolecular bend overtone.

Coupling of three intramolecular vibration modes of liquid H2O molecules in the framework of the fluctuation theory of hydrogen bonding

Vibrational spectra of liquid water contain a wealth of information about its structure and dynamics but there are no generally acknowledged interpretation of their band profiles, unfortunately. We have tried to calculate them taking into account only the coupling of three intramolecular vibrations in the set of H₂O molecules with different initial oscillator frequencies and intensities. The matter is that each water molecule forms hydrogen bonds of different strengths; thus the OH stretching band spans several hundreds of wave numbers (the fluctuation theory of hydrogen bonding). This distribution overlaps with the similar band of the first overtone of the HOH bending frequencies thus triggering a Fermi resonance between three vibrations. There were some problems causing some simplifications in previous theoretical modeling of vibrational transitions in condensed water. To solve them we extract the statistical distribution of OH frequencies of H₂O molecules directly from the experimental spectra of HOD molecules at the same conditions instead of defining it theoretically. Also the bending overtone is allowed to have non-zero intrinsic intensity when calculating the Fermi resonance. The test calculation of the isotropic component of the Raman spectrum which is most critical for interpretation shows that our algorithm can reproduces the characteristic peculiarities of the experiment. The spectrum consists of three non-Gaussian contours. The overtone of bending vibration, being strengthened by Fermi resonance, makes the greatest contribution to the dominating low-frequency spectrum component that was previously attributed to intermolecular coupling of adjacent OH oscillators by some authors. Further we plan to calculate similar band profiles in the IR and also Raman isotropic and anisotropic spectra of Н₂О and D₂О molecules within a wide temperature range for their quantitative comparison with experiment.

Vibrational couplings and energy transfer pathways of water's bending mode

Coupling between vibrational modes is essential for energy transfer and dissipation in condensed matter. For water, different O-H stretch modes are known to be very strongly coupled both within and between water molecules, leading to ultrafast dissipation and delocalization of vibrational energy. In contrast, the information on the vibrational coupling of the H-O-H bending mode of water is lacking, even though the bending mode is an essential intermediate for the energy relaxation pathway from the stretch mode to the heat bath. By combining static and femtosecond infrared, Raman, and hyper-Raman spectroscopies for isotopically diluted water with ab initio molecular dynamics simulations, we find the vibrational coupling of the bending mode differs significantly from the stretch mode: the intramode intermolecular coupling of the bending mode is very weak, in stark contrast to the stretch mode. Our results elucidate the vibrational energy transfer pathways of water. Specifically, the librational motion is essential for the vibrational energy relaxation and orientational dynamics of H-O-H bending mode.

Vibrational energy transfer in water involves intermolecular coupling of O-H stretching modes, but much less is known about the role of the bending modes. Here the authors, combining static and femtosecond infrared, Raman, and hyper-Raman spectroscopy and ab initio molecular dynamics simulations, provide insight into the energy dynamics of the bend vibrations.

Visualization of the in situ distribution of contents and hydrogen bonding states of cellular level water in apple tissues by confocal Raman microscopy

Raman spectroscopy has been employed for studying the hydrogen bonding states of water molecules for decades, however, Raman imaging data contain thousands of spectra, making it challenging to obtain information on water with different hydrogen bonds. In the current study, a novel method combining confocal Raman microscopy (CRM) imaging with the iterative curve fitting algorithms was developed to determine the distribution of water contents at the cellular level and water states with different hydrogen bonds in apple tissues. Raman imaging data ranging from 2700 to 3800 cm-1 were acquired from whole cells in the apple tissue, which were then decomposed into seven sub-peaks using the fixed-position Gaussian iterative curve fitting (FPGICF) algorithm. The content and hydrogen bonding states of cellular water were calculated as the area sum of the OH stretching vibration and the area ratio of DA-OH over DDAA-OH stretching vibration or the number of hydrogen bonds of each water molecule, respectively. Finally, the area of each sub-peak, the area sum of the OH stretching vibration, and the area ratio of DA-OH over DDAA-OH stretching vibration were used to visualize the distribution of each sub-peak, water contents and water states with different hydrogen bonds, respectively. In addition, it was found that the number of hydrogen bonds of each water molecule could also be considered as a criterion to describe the hydrogen bond states of water in apple tissues. The availability of such information should provide new insights for future study of cellular water in other food materials.

Temperature and Nuclear Quantum Effects on the Stretching Modes of the Water Hexamer

The water hexamer has many low-lying isomers, e.g., ring, book, cage, and prism, shifting from two- to three-dimensional structures. We show that this dimensionality change is accompanied by a drop in the quantum nature of the cluster, as manifested in the red shift of the quantal OH stretching modes as compared with their classical counterparts. We obtain this "nuclear quantum effect" (NQE) as the mean deviation between the OH stretch frequencies from velocity autocorrelation Fourier transforms from classical trajectories on a high-level water potential (MB-pol) as compared with scaled harmonic frequencies from high-level quantum chemistry calculations. With a universal scaling factor, the predicted OH frequencies agree with experiment to a mean absolute deviation ≤10 cm-1, which allows unequivocal isomer assignments. By assuming temperature-independent NQEs, we produce the temperature dependence of the cage isomer OH stretch spectrum below 70 K, where it is the dominant structure. All bands widen and blue-shift with increasing temperature, most conspicuously the reddest mode, which thus constitutes a "vibrational thermometer".

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.

Icelike Vibrational Properties of Strong Hydrogen Bonds in Hydrated Lithium Nitrate

The hydrogen bond network accounts for many of the extraordinary physical properties of liquid water and ice. Its vibrational dynamics are quite complex in their entirety but can be accessed in detail by investigating small groups of only a few water molecules. Here, aqueous salt hydrates turned out to be an exceptional model system for water molecules arranged in well-defined geometrical structures that can be accessed by means of femtosecond spectroscopy of the OH stretching vibration. In this study, we find striking resemblance between the vibrational properties of three water molecules connected via strong hydrogen bonds in the trihydrate of LiNO₃ and those of ordinary ice I_h. As in ice, the vibrations of the hydrate water molecules show ultrafast excited state dynamics that are strongly accelerated when proceeding from deuterated to neat H₂O samples. The latter is analyzed by means of an additional relaxation channel that is due to Fermi resonance between the OH stretching vibration and the bend overtone accompanied by delocalization of the vibration over neighboring water molecules in the H₂O species. Moreover, in the hydrate and ice samples severe spectral broadening is examined when comparing fundamental and excited state absorption bands. Here, proton delocalization along the strong hydrogen bonds is given as a possible underlying mechanism.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

Formation and characterization of protein-based films from yellow pea (Pisum sativum) protein isolate and concentrate for edible applications

This study investigated the properties of films or bioplastics fabricated using a wet processing method from yellow pea protein isolate (YPI) and protein concentrate (YPC) for potential application in food packaging. The wet processing method included mixing the protein with water and glycerol followed by casting and drying the films in a humidity- and temperature-controlled chamber. Whey protein isolate (WPI) and a film from a blend of equal amounts of YPI and WPI, labelled as YPI + WPI, were also studied. Fourier transform-infra red analysis revealed that films from YPI, YPC, WPI and YPI + WPI were formed by protein polymerisation with the plasticiser, glycerol, via hydrophobic and hydrophilic interactions. The protein films had contact angles of <90° demonstrating that they had a hydrophilic surface, with YPC < YPI < YPI + WPI < WPI. The pattern of ultraviolent light transmission of the films was WPI > YPC > YPI + WPI > YPI, whereas the mechanical and thermal resilience of films formulated from YPI, YPC and the protein blend were comparable to the properties of WPI-based films. The findings demonstrate that yellow pea proteins can be used as biomaterials to develop protein and protein-blend films or bioplastics for food packaging and edible applications.

•

Bioplastics were fabricated from yellow pea protein isolate and concentrate, with glycerol.

•

Contact angles of pea protein films indicate more hydrophobic surface than whey protein films.

•

Pea protein films had more surface structure homogeneity and limited light transmission.

•

Pea + whey protein blend did not produce synergistic effects in film property.

•

Film physico-mechanical properties are promising for food packaging application.

Unveiling Heterogeneity of Interfacial Water through the Water Bending Mode

The water bending mode provides a powerful probe of the microscopic structure of bulk aqueous systems because its frequency and spectral line shape are responsive to the intermolecular interactions. Furthermore, interpreting the bending mode response is straightforward, as the intramolecular vibrational coupling is absent. Nevertheless, bending mode has not been used for probing the interfacial water structure, as it has been yet argued that the signal is dominated by bulk effects. Here, through the sum-frequency generation measurement of the water bending mode at the water/air and water/charged lipid interfaces, we demonstrate that the bending mode signal is dominated not by the bulk but by the interface. Subsequently, we disentangle the hydrogen-bonding of water at the water/air interface using the bending mode frequency distribution and find distinct interfacial hydrogen-bonded structures, which can be directly related to the interfacial organization of water. The bending mode thus provides an excellent probe of aqueous interfacial structure.

Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk

The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H₂O)₁₀₀H⁺, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.

Hiding in the Crowd: Spectral Signatures of Overcoordinated Hydrogen-Bond Environments

Molecules with an excess number of hydrogen-bonding partners play a crucial role in fundamental chemical processes, ranging from anomalous diffusion in supercooled water to transport of aqueous proton defects and ordering of water around hydrophobic solutes. Here we show that overcoordinated hydrogen-bond environments can be identified in both the ambient and supercooled regimes of liquid water by combining experimental Raman multivariate curve resolution measurements and machine learning accelerated quantum simulations. In particular, we find that OH groups appearing in spectral regions usually associated with non-hydrogen-bonded species actually correspond to hydrogen bonds formed in overcoordinated environments. We further show that only these species exhibit a turnover in population as a function of temperature, which is robust and persists under both constant pressure and density conditions. This work thus provides a new tool to identify, interpret, and elucidate the spectral signatures of crowded hydrogen-bond networks.

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.