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

Investigation of ion association kinetics in Na2SO4-H2O by excess Raman spectroscopy and 2D correlation Raman spectroscopy

The ion association of salts aqueous solutions have long captivated the attention of researchers within the field of physical chemistry. In this paper, we have performed a comprehensive analysis of ion interactions in sodium sulfate (Na2SO4) aqueous solutions using a combination of high-resolution techniques, including excess (ERS) and two-dimensional correlation (2DCRS) Raman spectroscopy in conjunction with molecular dynamics (MD) calculations. The Raman spectrum shows that two inflection points in the Raman shift of the O-H vibration are observed with the increase in Na2SO4 concentration. Simultaneously, a new peak of the SO42- vibration appears at first inflection point, representing the formation of ion association. Further analysis based on ERS and 2CRS reveals that these two inflection points correspond respectively to the formation of ion pairs (CIPs) and small ion clusters. Importantly, MD simulations confirm the above experimental results. Our study provides evidence for ion association and clustering in aqueous in salt ion aqueous solutions.

Dissecting the hydrogen bond network of water: Charge transfer and nuclear quantum effects

The molecular structure of water is dynamic, with intermolecular (H)-bond interactions being modified by both electronic charge transfer and nuclear quantum effects (NQEs). Electronic charge transfer and NQEs potentially change under acidic / basic conditions, but such details have not been measured. Here, we developed correlated vibrational spectroscopy, a symmetry-based method that distinctively separates interacting from non-interacting molecules in self- and cross-correlation spectra, giving access to previously inaccessible information. We found that OH- donated ~8% more negative charge to the H-bond network of water and H3O+ accepted ~4% less negative charge from the H-bond network of water. D2O had ~9% more H-bonds compared to H2O, and acidic solutions displayed more dominant NQEs than basic ones.

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.

i-PI 3.0: A flexible and efficient framework for advanced atomistic simulations

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.

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.

Learning Electronic Polarizations in Aqueous Systems

The polarization of periodically repeating systems is a discontinuous function of the atomic positions, a fact which seems at first to stymie attempts at their statistical learning. Two approaches to build models for bulk polarizations are compared: one in which a simple point charge model is used to preprocess the raw polarization to give a learning target that is a smooth function of atomic positions and the total polarization is learned as a sum of atom-centered dipoles and one in which instead the average position of Wannier centers around atoms is predicted. For a range of bulk aqueous systems, both of these methods perform perform comparatively well, with the former being slightly better but often requiring an extra effort to find a suitable point charge model. As a challenging test, we also analyze the performance of the models at the air-water interface. In this case, while the Wannier center approach delivers accurate predictions without further modifications, the preprocessing method requires augmentation with information from isolated water molecules to reach similar accuracy. Finally, we present a simple protocol to preprocess the polarizations in a data-driven way using a small number of derivatives calculated at a much lower level of theory, thus overcoming the need to find point charge models without appreciably increasing the computation cost. We believe that the training strategies presented here help the construction of accurate polarization models required for the study of the dielectric properties of realistic complex bulk systems and interfaces with ab initio accuracy.

TMAO perturbs intermolecular vibrational motions of water revealed by low-frequency modes

Trimethylamine N-oxide (TMAO) as a representative natural osmolyte has received much attention because of its unique properties, including enhancement of hydrogen bonding networks in liquid water and stabilization of three-dimensional structures of proteins in living organisms. As a hydrogen bond maker and/or a protein stabilizer, its hydrated structures and orientation dynamics in aqueous solutions have been investigated by various spectroscopic methods. Particularly, distinct from other natural osmolytes, it has been found that TMAO molecules form complexes with water molecules even at low concentrations, showing extraordinarily long lifetimes and much larger effective dipole moments. In this study, we demonstrated that collective motions of water molecules are closely correlated to TMAO molecules, as revealed by the changes of the librational modes observed in hyper-Raman (HR) spectra in the low-frequency region (<1000 cm-1) for the first time. Based on HR spectra of the TMAO solutions at submolar concentrations, we observed that the librational bands originating from water apparently upshift (∼15 cm-1) upon the addition of TMAO molecules. Compared to the OH stretching band of water showing a negligible downshift (<5 cm-1), the librational bands of water are more sensitive to reflect changes in the hydrogen bonding networks in the TMAO solutions, suggesting formation of transient TMAO-water complexes plays an essential role toward surrounding water molecules in perturbing their librational motions. We expect to provide a supplementary approach to understand that water molecules in TMAO aqueous solutions are strongly affected by TMAO molecules, different from other osmolytes.

Correlated libration in liquid water

The libration spectrum of liquid H2O is resolved into an octupolar twisting libration band at 485 cm-1 and dipolar rocking-wagging libration bands at 707 and 743 cm-1 using polarization analysis of the hyper-Raman scattering (HRS) spectrum. Dipole interactions and orientation correlation over distances less than 2 nm account for the 36 cm-1 splitting of the longitudinal and transverse polarized bands of the dipolar rocking-wagging libration mode, while the intensity difference observed for the bands is the result of libration correlation over distances larger than 200 nm. The coupled rock and wag libration in water is similar to libration modes in ice. The libration relaxation time determined from the width of the spectrum is 36-54 fs. Polarization analysis of the HRS spectrum also shows long range correlation for molecular orientation and hindered translation, bending and stretching vibrations in water.

Hyper-Raman spectroscopy of biomolecules

The hyper-Raman scattering (HRS) spectra of biologically significant molecules (D-glucose, L-alanine, L-arabinose, L-tartaric acid) in aqueous solutions are reported. The HRS spectra were measured using a picosecond laser at 532 nm operating at a MHz repetition rate. High signal to noise spectra were collected with a commercial spectrometer and CCD without resonant, nanoparticle, or surface enhancement. The HRS peak frequencies, relative intensities, band assignments, and depolarization ratios are examined. By comparing HRS to Raman scattering (RS) and infrared absorption spectra we verify that the IR-active vibrational modes of the target molecules are observed in HRS spectra but come with substantially different peak intensities. The HRS of the biomolecules as well as water, dimethyl sulfoxide, methanol, and ethanol were deposited into a data repository to support the development of theoretical descriptions of HRS for these molecules. Depositing the spectra in a repository also supports future dual detection RS, HRS microscopes which permit simultaneous high-spatial-resolution vibrational spectroscopy of IR-active and Raman-active vibrational modes.

Bulklike Vibrational Coupling of Surface Water Revealed by Sum-Frequency Generation Spectroscopy

Vibrational coupling between interfacial water molecules is important for energy dissipation after on-water chemistry, yet intensely debated. Here, we quantify the interfacial vibrational coupling strength through the linewidth of surface-specific vibrational spectra of the water's O─H (O─D) stretch region for neat H_{2}O/D_{2}O and their isotopic mixtures. The local-field-effect-corrected experimental SFG spectra reveal that the vibrational coupling between hydrogen-bonded interfacial water O─H groups is comparable to that in bulk water, despite the effective density reduction at the interface.

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.

Characterization of Secondary Structures of Model Polypeptides in Solutions with Hyper-Raman Spectroscopy

Characterization of the secondary structures of two model polypeptides, poly-l-lysine and poly-l-glutamic acid in aqueous solutions has been demonstrated by hyper-Raman (HR) spectroscopy for the first time. Complementary to infrared (IR) and visible Raman spectroscopy, HR spectroscopy gives the amide I, II, and III bands originating from the polypeptide backbones and the CCH₃ symmetric bending mode, enabling us to distinguish different conformations. The α-helix gives the broad and weak amide III band, while the β-sheet and the random coil show similar spectral patterns with different relative intensities between the amide I and II bands. HR spectra from aqueous solutions of the α-helix and the random coil of poly-l-ornithine also possess these spectral features. The HR spectra are analogous to UV resonance Raman (UVRR) spectra, indicating the signal enhancement due to the electronic resonance effect via the π-π* transition. In contrast, the vibrational frequencies of the amide I band in the HR spectra are much higher than those in the IR, visible Raman, and UVRR spectra, suggesting the non-coincidence between HR, IR, and Raman bands. Our finding suggests that HR spectroscopy is promising to provide complementary information on the secondary structures of polypeptides in aqueous solutions as a spectral approach differing from existing vibrational spectroscopic methods.

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.