Both Inter- and Intramolecular Coupling of O-H Groups Determine the Vibrational Response of the Water/Air Interface

How Thick is the Air-Water Interface?─A Direct Experimental Measurement of the Decay Length of the Interfacial Structural Anisotropy

The air-water interface is a highly prevalent phase boundary impacting many natural and artificial processes. The significance of this interface arises from the unique properties of water molecules within the interfacial region, with a crucial parameter being the thickness of its structural anisotropy, or "healing depth". This quantity has been extensively assessed by various simulations which have converged to a prediction of a remarkably short length of ∼6 Å. Despite the absence of any direct experimental measurement of this quantity, this predicted value has surprisingly become widely accepted as fact. Using an advancement in nonlinear vibrational spectroscopy, we provide the first measurement of this thickness and, indeed, find it to be ∼6-8 Å, finally confirming the prior predictions. Lastly, by combining the experimental results with depth-dependent second-order spectra calculated from ab initio parametrized molecular dynamics simulations, which are also in excellent agreement with this experimental result, we shed light on this surprisingly short correlation length of molecular orientations at the interface.

De Novo Design and Characterization of Amphiphilic Peptides with Basic Side Chains for Tailored Interfacial Chemistries

Lysine-leucine (LK) peptides have been used as model systems and platforms for 2D material design for decades. LK peptides are amphiphilic sequences designed to bind and fold at hydrophobic surfaces through hydrophobic leucine side chains and hydrophilic lysine side chains extending into the aqueous subphase. The hydrophobic periodicity of the sequence dictates the secondary structure at the interface. This robust design makes them ideal candidates for controlling interfacial chemistry. This study presents the de novo design and characterization of two novel peptides: LRα14 and LHα14, which substitute lysine with arginine and histidine, respectively, in the helical LKα14 sequence. This modification is intended to expand the LK peptide platform to a new basic interfacial chemistry. We explore the stability of the new LRα14 and LHα14 designs with respect to changes in pH and salt concentration in bulk solution and at the interface using circular dichroism (UV-CD) and vibrational sum-frequency generation spectroscopy, respectively. Notably, the structural stability of the peptides remains unaffected across a wide range of pH and ionic strength values. At the same time, the variation of side-chain chemistry leads to a wide spectrum of interfacial water structures. By extension of the LK platform to include arginine and histidine, this study broadens the toolbox for designing tailored interfacial chemistries with applications in material and biomedical sciences.

Fully Quantum Simulation of Polaritonic Vibrational Spectra of Large Cavity-Molecule System

The formation of molecular vibrational polaritons, arising from the interplay between molecular vibrations and infrared cavity modes, is a quantum phenomenon necessitating accurate quantum dynamical simulations. Here, we introduce the cavity vibrational self-consistent field/virtual state configuration interaction method, enabling quantum simulation of the vibrational spectra of many-molecule systems within the optical cavity. Focusing on the representative (H2O)21 system, we showcase this parameter-free quantum approach's ability to capture both linear and nonlinear vibrational spectral features. Our findings highlight the growing prominence of molecular couplings among OH stretches and bending excited bands with increased light-matter interaction, revealing distinctive nonlinear spectral features induced by vibrational strong coupling.

Surface stratification determines the interfacial water structure of simple electrolyte solutions

The distribution of ions at the air/water interface plays a decisive role in many natural processes. Several studies have reported that larger ions tend to be surface-active, implying ions are located on top of the water surface, thereby inducing electric fields that determine the interfacial water structure. Here we challenge this view by combining surface-specific heterodyne-detected vibrational sum-frequency generation with neural network-assisted ab initio molecular dynamics simulations. Our results show that ions in typical electrolyte solutions are, in fact, located in a subsurface region, leading to a stratification of such interfaces into two distinctive water layers. The outermost surface is ion-depleted, and the subsurface layer is ion-enriched. This surface stratification is a key element in explaining the ion-induced water reorganization at the outermost air/water interface.

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.

Role of Water during the Early Stages of Iron Oxyhydroxide Formation by a Bacterial Iron Nucleator

Understanding the nucleation of iron oxides and the underlying hydrolysis of aqueous iron species is still challenging, and molecular-level insights into the orchestrated response of water, especially at the hydrolysis interface, are lacking. We follow iron(III) hydrolysis in the presence of a synthetic bacterial iron nucleator, which is a magnetosome membrane specific peptide, by using a constant pH titration technique. Three distinct hydrolysis regimes were identified. Interface-selective sum frequency generation (SFG) spectroscopy was used to probe the interfacial reaction and water in direct contact with the peptide. SFG data reveal that iron(III) species react quickly with interfacial peptides while continuously enhancing water alignment into the later stages of hydrolysis. The gradually aligning water molecules are associated with initially promoted (regimes I and II) and later suppressed (regime III) hydrolysis after the saturation of water alignment has occurred until regime II. These interfacial insights are crucial for understanding the early stage of iron oxide biomineralization.

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.

Unravelling the interfacial water structure at the photocatalyst strontium titanate by sum frequency generation spectroscopy

The direct conversion of solar energy to hydrogen is considered as a possible method to produce carbon neutral hydrogen fuel. The mechanism of photocatalytic water splitting involves the chemical breakdown of water and re-assembly into hydrogen and oxygen at the interface of a photocatalyst. The selection rules of a suitable material are well established, but the fundamental understanding of the mechanisms, occurring at the interface between the catalyst and the water, remains missing. Using surface specific sum frequency generation spectroscopy, we present here characterisation of the interface between water and the photocatalyst strontium titanate (SrTiO3). We monitor the OH-stretching vibrations present at the interface. Their variations of intensities and frequencies as functions of isotopic dilution, pH and salt concentration provide information about the nature of the hydrogen bonding environment. We observe the presence of water molecules that flip their orientation at pH 5 indicating the point of zero charge of the SrTiO3 layer. These water molecules are oriented with their hydrogen away from the surface when the pH of the solutions is below 5 and pointing towards the surface when the pH is higher than 5. Besides, water molecules donating a H-bond to probably surface TiOH groups are observed at all pH.

How Does pH Affect the Adsorption of Human Serum Protein in the Presence of Hydrophobic and Hydrophilic Nanoparticles at Air-Water and Lipid-Water Interfaces?

This study investigates interaction between hydrophilic (11-mercaptoundecanoic acid (MUA)) and hydrophobic (1-undecanethiol (UDT)) gold nanoparticles (GNPs) with human serum albumin (HSA) protein on air-water and lipid-water interfaces at pH 3 and 7. Vibrational sum frequency generation (VSFG) spectroscopy is used to analyze changes in the intensity of interfacial water molecules and the C-H group of the protein. At the air-water interface, the hydrophobic interaction between the HSA protein and hydrophobic GNPs at pH 3 leads to their accumulation at the interface, resulting in an increased C-H intensity of the protein with a slight decrease in water intensity. Whereas, at pH 7, where the negative charge of the protein results in the reduced surface activity of the HSA compared to pH 3, the interaction between alkyl chain of the hydrophobic GNPs and alkyl group of the protein results in the adsorption of the protein-capped GNPs at the interface. This leads to an increased intensity of the C-H group of protein and water molecules. However, negatively charged hydrophilic GNPs do not induce significant changes in the interfacial water structure or the C-H group of the protein due to the electrostatic force of repulsion with the negatively charged HSA at pH 7. In contrast, at the lipid-water interface, both hydrophobic and hydrophilic GNPs interact with HSA protein, causing disordering of interfacial water molecules at pH 3 and ordering at pH 7. Interestingly, similar behavior of the protein with both types of GNPs results in comparable ordering/disordering at the interface depending on the pH of solution. Furthermore, the VSFG results obtained with the deuterated lipid suggest that changes in ordering and disorder occur due to increased protein adsorption in the presence of GNPs, causing alterations in the membrane structure. These findings give a better understanding of the mechanisms that govern protein-nanoparticle interaction and their consequential effects on the structure, function, and behavior of molecules at the biological membrane interface, which is crucial for developing safe and effective nanoparticle-based therapeutics.

Fully First-Principles Surface Spectroscopy with Machine Learning

Our current understanding of the structure and dynamics of aqueous interfaces at the molecular level has grown substantially due to the continuous development of surface-specific spectroscopies, such as vibrational sum-frequency generation (VSFG). As in other vibrational spectroscopies, we must turn to atomistic simulations to extract all of the information encoded in the VSFG spectra. The high computational cost associated with existing methods means that they have limitations in representing systems with complex electronic structure or in achieving statistical convergence. In this work, we combine high-dimensional neural network interatomic potentials and symmetry-adapted Gaussian process regression to overcome these constraints. We show that it is possible to model VSFG signals with fully ab initio accuracy using machine learning and illustrate the versatility of our approach on the water/air interface. Our strategy allows us to identify the main sources of theoretical inaccuracy and establish a clear pathway toward the modeling of surface-sensitive spectroscopy of complex interfaces.

Inconsistent hydrogen bond-mediated vibrational coupling of amide I

Using infrared spectroscopy and density functional theory (DFT) calculations, we scrutinized an amide (dimethylformamide) as a "model" compound to interpret the interactions of amide 1 with different phenol derivatives (para-chlorophenol (PCP) and para-cresol (CP)) as "model guest molecules". We established the involvement of amide I in vibrational coupling with symmetric and asymmetric C[double bond, length as m-dash]C modes of different phenolic derivatives and how their coupling was dependent upon different guest aromatic phenolic compounds. Interestingly, substitution of phenol perturbed the pattern of vibrational coupling with amide I. The symmetric and asymmetric C[double bond, length as m-dash]C modes of PC were coupled significantly with amide 1. For PCP, the symmetric C[double bond, length as m-dash]C mode coupled significantly, but the asymmetric mode coupled negligibly, with amide I. Here, we reveal the nature of vibrational coupling based on the structure of a guest molecule hydrogen-bonded with amide I. Our conclusions could be valuable for depiction of the unusual dynamics of coupled amide-I modes as well as the dependency of vibrational coupling on altered factors.

Vibrational Sum Frequency Generation Spectra of Water-Vapor Interfaces Covered by Alcohols: Effects of Surface Coverage and Coupling between Oscillators

The present study deals with the effects of varying coverage of water surface by alcohols on the vibrational sum frequency generation (VSFG) spectrum of interfacial water. We have considered two different alcohols: Tertiary butyl alcohol (TBA) whose alkyl part is fully branched and stearyl alcohol (STA) which has a long linear alkyl chain with larger hydrophobic surface area than that of TBA. With increase of the alcohol concentration, the hydrogen bonded OH stretch region of the VSFG spectrum is found to change following a regular trend for the STA-water system, whereas non-monotonic variation of the VSFG spectrum is observed for the TBA-water system which can be correlated with the presence of very different interactions of TBA molecules at different concentrations. On increasing the concentration of TBA, the hydrophobic groups get more tilted towards the water phase and significant hydrophobic interactions are introduced at higher concentrations. Whereas, for STA, there is a gradual increase in the hydrophilic interaction. Because of stacking interactions between the long chain alkyl groups, the hydrophobic parts stay outward from the water phase at higher concentrations and a regular change in the VSFG spectrum is observed. We have also presented a computationally efficient scheme to calculate the VSFG spectrum of interfacial systems for coupled oscillators which is expected to be beneficial for the treatment of coupling where the interfacial system size is inherently large.

Surface Properties of Saponin-Chitosan Mixtures

The surface properties of saponin and saponin-chitosan mixtures were analysed as a function of their bulk mixing ratio using vibrational sum-frequency generation (SFG), surface tensiometry and dilational rheology measurements. Our experiments show that saponin-chitosan mixtures present some remarkable properties, such as a strong amphiphilicity of the saponin and high dilational viscoelasticity. We believe this points to the presence of chitosan in the adsorption layer, despite its complete lack of surface activity. We explain this phenomenon by electrostatic interactions between the saponin as an anionic surfactant and chitosan as a polycation, leading to surface-active saponin-chitosan complexes and aggregates. Analysing the SFG intensity of the O-H stretching bands from interfacial water molecules, we found that in the case of pH 3.4 for a mixture consisting of 0.1 g/L saponin and 0.001 g/L chitosan, the adsorption layer was electrically neutral. This conclusion from SFG spectra is corroborated by results from surface tensiometry showing a significant reduction in surface tension and effects on the dilational surface elasticity strictly at saponin/chitosan ratios, where SFG spectra indicate zero net charge at the air-water interface.

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.

Infrared Characterization of Isotopic Analogues of Methanediol in Aqueous Solution

Dissolved methanediol in aqueous solution has been treated as the precursor for the formation of atmospheric formic acid in multiphase environments. In this work, methanediol, CH₂(OH)₂, and its isotopic analogues, CH₂(OD)₂, CD₂(OH)₂, and CD₂(OD)₂, in aqueous solution were prepared by dissolving paraformaldehyde and deuterium-substituted paraformaldehyde powders in H₂O and D₂O under reflux. Their infrared absorption contours of formaldehyde solutions at concentrations of <1 wt % are not dependent on the concentration, mainly referring to the characteristics of the monomeric configuration, and can be categorized into two parts. At wavenumbers >2000 cm^-1, broad bands of moderate strengths were ascribed to the stretching modes of two OH or OD groups, observed at 3000-3700 and 2050-2750 cm^-1, respectively. At wavenumbers of 950-1200 cm^-1, the isotopic analogues of methanediols composed of CH₂ moieties are featured with a singlet strong band at ca. 1030 cm^-1, mainly attributed to the O-C-O stretching modes; the isotopic methanediols containing CD₂ moieties manifested two intense bands at ca. 1100 and 980 cm^-1, majorly enveloping the CD₂ deformation and O-C-O stretching modes. The aforementioned spectral features were assigned on the basis of density functional theory, ωB97XD, with the basis set aug-cc-pVTZ and the solvent effect using the conductor-like polarizable continuum model. In addition, the predicted energetics suggested that the trans-methanediol is more stable than the cis- conformer by ca. 0.62 kcal mol^-1 and majorly contributes to the infrared features. At higher concentrations of CH₂(OH)₂, extra bands at 920 and 1104 cm^-1 appeared and were attributed to the C-O-C stretching modes of the dimeric/polymeric methanediol; that is, HO(CH₂O)_nH, n ≥ 2. These infrared characterizations of the isotopic analogues of methanediols provided suitable detection windows in the relevant atmospheric and aerosol reactions in the laboratory studies.

Water Structure at the Hydrophobic Nanodroplet Surface Revealed by Vibrational Sum Frequency Scattering Using Isotopic Dilution

The water structure at the hydrophobic/water interface is key toward understanding hydrophobicity at the molecular level. Herein, we characterize the hydrogen-bonding network of interfacial water next to sub-micron-sized hydrophobic oil droplets dispersed in water using isotopic dilution vibrational sum frequency scattering (SFS) spectroscopy. The relative intensity of different modes, the frequency shift of the uncoupled O-D spectrum, and a low-frequency shoulder (2395 cm-1) reveal that water forms an overall stronger hydrogen-bonding network next to hydrophobic droplets compared to bulk water and the air/water interface. Half of the spectral width of the oil droplet SFS spectrum is determined by inter- and intramolecular coupling of water molecules. Isotopic dilution also confirms the presence of a broad distribution (ca. 2640-2745 cm-1) of non-water-hydrogen-bonded O-D modes that are red-shifted and broadened compared to similar species at the air/water interface. This band corroborates the presence of charge transfer between water and oil.

Passively Stabilized Phase-Resolved Collinear SFG Spectroscopy Using a Displaced Sagnac Interferometer

Sum-frequency generation (SFG) vibrational spectroscopy is a powerful technique to study interfaces at the molecular level. Phase-resolved SFG (PR-SFG) spectroscopy provides direct information on interfacial molecules' orientation. However, its implementation is technologically demanding: it requires the generation of a local oscillator wave and control of its time delay with sub-fs accuracy. Commonly used noncollinear PR-SFG provides this control naturally but requires very accurate sample height control. Collinear PR-SFG spectroscopy is less demanding regarding sample positioning, but tuning the local oscillator time delay with this beam geometry is challenging. Here, we develop a collinear PR-SFG setup using a displaced Sagnac interferometer. This scheme allows full, independent control of the time delay and intensity of the local oscillator and provides long-time phase stabilization (better than 5° over 12 h) for the measured signal. This approach substantially reduces the complexity of an experimental setup and combines the advantages of collinear and noncollinear PR-SFG techniques.

Effects of Stearyl Alcohol Monolayer on the Structure, Dynamics and Vibrational Sum Frequency Generation Spectroscopy of Interfacial Water

The structure, dynamics and vibrational spectroscopy of water surface covered by a monolayer of stearyl alcohol (STA) are investigated by means of molecular dynamics simulations and vibrational sum frequency generation...

Interfacial Water Structure of Binary Liquid Mixtures Reflects Nonideal Behavior

The evaporation of molecules from water-organic solute binary mixtures is key for both atmospheric and industrial processes such as aerosol formation and distillation. Deviations from ideal evaporation energetics can be assigned to intermolecular interactions in solution, yet evaporation occurs from the interface, and the poorly understood interfacial, rather than the bulk, structure of binary mixtures affects evaporation kinetics. Here we determine the interfacial structure of nonideal binary mixtures of water with methanol, ethanol, and formic acid, by combining surface-specific vibrational spectroscopy with molecular dynamics simulations. We find that the free, dangling OH groups at the interfaces of these differently behaving nonideal mixtures are essentially indistinguishable. In contrast, the ordering of hydrogen-bonded interfacial water molecules differs substantially at these three interfaces. Specifically, the interfacial water molecules become more disordered (ordered) in mixtures with methanol and ethanol (formic acid), showing higher (lower) vapor pressure than that predicted by Raoult's law.

Using Diffusion Monte Carlo Wave Functions to Analyze the Vibrational Spectra of H7O3+ and H9O4. J Phys Chem A 2021;125:7185-7197. [PMID: 34433268 DOI: 10.1021/acs.jpca.1c05025]

An approach for evaluating spectra from ground state probability amplitudes (GSPA) obtained from diffusion Monte Carlo (DMC) simulations is extended to improve the description of excited state energies and allow for coupling among vibrational excited states. This approach is applied to studies of the protonated water trimer and tetramer, and their deuterated analogs. These ions provide models for solvated hydronium, and analysis of these spectra provides insights into spectral signatures of proton transfer in aqueous environments. In this approach, we obtain a separable set of internal coordinates from the DMC ground state probability amplitude. A basis is then developed from products of the DMC ground state wave function and low-order polynomials in these internal coordinates. This approach provides a compact basis in which the Hamiltonian and dipole moment matrix are evaluated and used to obtain the spectrum. The resulting spectra are in good agreement with experiment and in many cases provide comparable agreement to the results obtained using much larger basis sets. In addition, the compact basis allows for interpretation of the spectral features and how they evolve with cluster size and deuteration.

Water Orientation at the Calcite-Water Interface

Mineral-water interfaces play an important role in many natural as well as technological fields. Fundamental properties of these interfaces are governed by the presence of the interfacial water and its specific structure at the surface. Calcite is particularly interesting as a dominant rock-forming mineral in the earth's crust. Here, we combine atomic force microscopy, sum-frequency generation spectroscopy, and molecular dynamics simulations to determine the position and orientation of the water molecules in the hydration layers of the calcite surface with high resolution. While atomic force microscopy provides detailed information about the position of the water molecules at the interface, sum-frequency generation spectroscopy can deduce the orientation of the water molecules. Comparison of the calcite-water interface to the interfaces of magnesite-water, magnesite-ethanol, and calcite-ethanol reveals a comprehensive picture with opposite water orientations in the first and second layer of the interface, which is corroborated by the molecular dynamics simulations.

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.

Conformation-induced vibrational spectral dynamics of hydrogen peroxide and vicinal water molecules

We studied the conformation-induced spectral response of water molecules due to site-specific structural alterations of solvated hydrogen peroxide (H₂O₂) employing DFT-based first principles molecular dynamics (FPMD) simulations. Wavelet transform was used to determine the time-dependent frequencies of the hydroxyls of water molecules and the O-H stretch modes of H₂O₂. Shifts in the vibrational frequency of the hydrogen-bonded hydroxyls inside the solvation shell of H₂O₂ support multiple distinctive hydrogen bonding environments. This paper classifies two distinct hydrogen bond types inside the O-O_W solvation shell of H₂O₂, and the dynamical calculations provide a quantitative estimation of the relative hydrogen bond strength. We ascertain the reason for not observing the escape of water molecules from the hydrogen peroxide hydration shell, unlike the solvation shell of ionic solutions and neutral solutes. Besides, we provide a comprehensive analysis of the spectral shifts in the normalized frequency distribution, the time-dependent decay of frequency-frequency correlation functions, and the hydrogen bond length scale fluctuations. We also quantify the relative contribution of the cisoid and transoid conformers affecting the vibrational spectral signature of the vicinal water molecules. While the transoid conformers promote the hydrogen bonding interactions through the oxygen site (OH_W), the cisoid conformers facilitate hydrogen peroxide-water hydrogen bond formation through the hydrogen site (HO_W). These non-identical hydrogen bond associations stabilize hydrogen peroxide in water.

Role of Water in CaCO3 Biomineralization

Biomineralization occurs in aqueous environments. Despite the ubiquity and relevance of CaCO₃ biomineralization, the role of water in the biomineralization process has remained elusive. Here, we demonstrate that water reorganization accompanies CaCO₃ biomineralization for sea urchin spine generation in a model system. Using surface-specific vibrational spectroscopy, we probe the water at the interface of the spine-associated protein during CaCO₃ mineralization. Our results show that, while the protein structure remains unchanged, the structure of interfacial water is perturbed differently in the presence of both Ca²⁺ and CO₃^2– compared to the addition of only Ca²⁺. This difference is attributed to the condensation of prenucleation mineral species. Our findings are consistent with a nonclassical mineralization pathway for sea urchin spine generation and highlight the importance of protein hydration in biomineralization.

Isolating the Contributions of Specific Network Sites to the Diffuse Vibrational Spectrum of Interfacial Water with Isotopomer-Selective Spectroscopy of Cold Clusters

Decoding the structural information contained in the interfacial vibrational spectrum of water requires understanding how the spectral signatures of individual water molecules respond to their local hydrogen bonding environments. In this study, we isolated the contributions for the five classes of sites that differ according to the number of donor (D) and acceptor (A) hydrogen bonds that characterize each site. These patterns were measured by exploiting the unique properties of the water cluster cage structures formed in the gas phase upon hydration of a series of cations M⁺·(H₂O)_n (M = Li, Na, Cs, NH₄, CH₃NH₃, H₃O, and n = 5, 20-22). This selection of ions was chosen to systematically express the A, AD, AAD, ADD, and AADD hydrogen bonding motifs. The spectral signatures of each site were measured using two-color, IR-IR isotopomer-selective photofragmentation vibrational spectroscopy of the cryogenically cooled, mass selected cluster ions in which a single intact H₂O is introduced without isotopic scrambling, an important advantage afforded by the cluster regime. The resulting patterns provide an unprecedented picture of the intrinsic line shapes and spectral complexities associated with excitation of the individual OH groups, as well as the correlation between the frequencies of the two OH groups on the same water molecule, as a function of network site. The properties of the surrounding water network that govern this frequency map are evaluated by dissecting electronic structure calculations that explore how changes in the nearby network structures, both within and beyond the first hydration shell, affect the local frequency of an OH oscillator. The qualitative trends are recovered with a simple model that correlates the OH frequency with the network-modulated local electron density in the center of the OH bond.

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.

Probing the intermolecular coupled vibrations in a water cluster with inelastic electron tunneling spectroscopy

The hydrogen-bonding networks of water have strong intra- and intermolecular vibrational coupling which influences the energy dissipation and proton transfer in water. Disentangling and quantitative characterization of different coupling effects in water at a single-molecular level still remains a great challenge. Using tip-enhanced inelastic electron tunneling spectroscopy (IETS) based on low-temperature scanning tunneling microscopy, we report the direct quantitative assessment of the intermolecular coupling constants of the OH-stretch vibrational bands of an isolated water tetramer adsorbed on a Au(111)-supported NaCl(001) bilayer film. This is achieved by distinguishing various coupled modes of the H-bonded O-H stretching vibrations through tip-height dependent IET spectra. In contrast, such vibrational coupling is negligible in the half-deuterated water tetramer owing to the large energy mismatch between the OH and OD stretching modes. Not only do these findings advance our understanding on the effects of local environment on the intermolecular vibrational coupling in water, but also open up a new route for vibrational spectroscopic studies of extended H-bonded network at the single-molecular level.

Structure and solvation of confined water and water-ethanol clusters within microporous Brønsted acids and their effects on ethanol dehydration catalysis

Water networks confined within zeolites solvate clustered reactive intermediates and must rearrange to accommodate transition states that differ in size and polarity, with thermodynamic penalties that depend on the shape of the confining environment.

Aqueous-phase reactions within microporous Brønsted acids occur at active centers comprised of water-reactant-clustered hydronium ions, solvated within extended hydrogen-bonded water networks that tend to stabilize reactive intermediates and transition states differently. The effects of these diverse clustered and networked structures were disentangled here by measuring turnover rates of gas-phase ethanol dehydration to diethyl ether (DEE) on H-form zeolites as water pressure was increased to the point of intrapore condensation, causing protons to become solvated in larger clusters that subsequently become solvated by extended hydrogen-bonded water networks, according to in situ IR spectra. Measured first-order rate constants in ethanol quantify the stability of S_N2 transition states that eliminate DEE relative to (C₂H₅OH)(H⁺)(H₂O)_n clusters of increasing molecularity, whose structures were respectively determined using metadynamics and ab initio molecular dynamics simulations. At low water pressures (2–10 kPa H₂O), rate inhibition by water (–1 reaction order) reflects the need to displace one water by ethanol in the cluster en route to the DEE-formation transition state, which resides at the periphery of water–ethanol clusters. At higher water pressures (10–75 kPa H₂O), water–ethanol clusters reach their maximum stable size ((C₂H₅OH)(H⁺)(H₂O)_4–5), and water begins to form extended hydrogen-bonded networks; concomitantly, rate inhibition by water (up to –3 reaction order) becomes stronger than expected from the molecularity of the reaction, reflecting the more extensive disruption of hydrogen bonds at DEE-formation transition states that contain an additional solvated non-polar ethyl group compared to the relevant reactant cluster, as described by non-ideal thermodynamic formalisms of reaction rates. Microporous voids of different hydrophilic binding site density (Beta; varying H⁺ and Si–OH density) and different size and shape (Beta, MFI, TON, CHA, AEI, FAU), influence the relative extents to which intermediates and transition states disrupt their confined water networks, which manifest as different kinetic orders of inhibition at high water pressures. The confinement of water within sub-nanometer spaces influences the structures and dynamics of the complexes and extended networks formed, and in turn their ability to accommodate the evolution in polarity and hydrogen-bonding capacity as reactive intermediates become transition states in Brønsted acid-catalyzed reactions.

Orientation independent vibrational dynamics of lipid-bound interfacial water

Zwitterionic phospholipids are one of the main constituents of biological membranes. The electric field associated with the two opposite headgroup charges aligns water molecules in the headgroup region. Here, we study the role of water alignment on the sub-picosecond vibrational dynamics of lipid-bound water. To this end, we compare the dynamics of oppositely oriented water associated with, respectively, a phosphocholine (PC) headgroup and an inverse-phosphocholine with non-ethylated phosphate groups (CP). We find that the dynamics are independent of the water orientation, implying that the vibrational dynamics report on the local properties of the water molecules.

Decoding the molecular water structure at complex interfaces through surface-specific spectroscopy of the water bending mode

The structure of interfacial water determines atmospheric chemistry, wetting properties of materials, and protein folding. The challenge of investigating the properties of specific interfacial water molecules has frequently been confronted using surface-specific sum-frequency generation (SFG) vibrational spectroscopy using the O-H stretch mode. While perfectly suited for the water-air interface, for complex interfaces, a potential complication arises from the contribution of hydroxyl or amine groups of non-water species present at the surface, such as surface hydroxyls on minerals, or O-H and N-H groups contained in proteins. Here, we present a protocol to extract the hydrogen bond strength selectively of interfacial water, through the water bending mode. The bending mode vibrational frequency distribution provides a new avenue for unveiling the hydrogen bonding structure of interfacial water at complex aqueous interfaces. We demonstrate this method for the water-CaF₂ and water-protein interfaces. For the former, we show that this method can indeed single out water O-H groups from surface hydroxyls, and that with increasing pH, the hydrogen-bonded network of interfacial water strengthens. Furthermore, we unveil enhanced hydrogen bonding of water, compared to bulk water, at the interface with human serum albumin proteins, a prototypical bio-interface.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

Sum frequency generation, calculation of absolute intensities, comparison with experiments, and two-field relaxation-based derivation

The experimental sum frequency generation (SFG) spectrum is the response to an infrared pulse and a visible pulse and is a highly surface-sensitive technique. We treat the surface dangling OH bonds at the air/water interface and focus on the absolute SFG intensities for the resonant terms, a focus that permits insight into the consequences of some approximations. For the polarization combinations, the calculated linewidths for the water interface dangling OH SFG band at 3,700 [Formula: see text] are, as usual, too large, because of the customary neglect of motional narrowing. The integrated spectrum is used to circumvent this problem and justified here using a Kubo-like formalism and theoretical integrated band intensities rather than peak intensities. Only relative SFG intensities are usually reported. The absolute integrated SFG intensities for three polarization combinations for sum frequency, visible, and infrared beams are computed. We use molecular dynamics and the dipole and the polarizability matrix elements obtained from infrared and Raman studies of [Formula: see text]O vapor. The theoretical expressions for two of the absolute susceptibilities contain only a single term and agree with experiment to about a factor of 1.3, with no adjustable parameters. The Fresnel factors are included in that comparison. One of the susceptibilities contains instead four positive and negative terms and agrees less well. The expression for the SFG correlation function is normally derived from a statistical mechanical formulation using a time-evolving density matrix. We show how a derivation based on a two-field relaxation leads to the same final result.

Nature of Excess Hydrated Proton at the Water-Air Interface

Understanding the interfacial molecular structure of acidic aqueous solutions is important in the context of, e.g., atmospheric chemistry, biophysics, and electrochemistry. The hydration of the interfacial proton is necessarily different from that in the bulk, given the lower effective density of water at the interface, but has not yet been elucidated. Here, using surface-specific vibrational spectroscopy, we probe the response of interfacial protons at the water-air interface and reveal the interfacial proton continuum. Combined with spectral calculations based on ab initio molecular dynamics simulations, the proton at the water-air interface is shown to be well-hydrated, despite the limited availability of hydration water, with both Eigen and Zundel structures coexisting at the interface. Notwithstanding the interfacial hydrated proton exhibiting bulk-like structures, a substantial interfacial stabilization by -1.3 ± 0.2 kcal/mol is observed experimentally, in good agreement with our free energy calculations. The surface propensity of the proton can be attributed to the interaction between the hydrated proton and its counterion.

Vibrational mode frequency correction of liquid water in density functional theory molecular dynamics simulations with van der Waals correction

We develop a frequency correction scheme for the stretch and bending modes of liquid water, which substantially improves the prediction of the vibrational spectra.

Capturing intrinsic site-dependent spectral signatures and lifetimes of isolated OH oscillators in extended water networks

The extremely broad infrared spectrum of water in the OH stretching region is a manifestation of how profoundly a water molecule is distorted when embedded in its extended hydrogen-bonding network. Many effects contribute to this breadth in solution at room temperature, which raises the question as to what the spectrum of a single OH oscillator would be in the absence of thermal fluctuations and coupling to nearby OH groups. We report the intrinsic spectral responses of isolated OH oscillators embedded in two cold (~20 K), hydrogen-bonded water cages adopted by the Cs⁺·(HDO)(D₂O)₁₉ and D₃O⁺·(HDO)(D₂O)₁₉ clusters. Most OH oscillators yield single, isolated features that occur with linewidths that increase approximately linearly with their redshifts. Oscillators near 3,400 cm^-1, however, occur with a second feature, which indicates that OH stretch excitation of these molecules drives low-frequency, phonon-type motions of the cage. The excited state lifetimes inferred from the broadening are considered in the context of fluctuations in the local electric fields that are available even at low temperature.

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.

Hydrogen bonding and molecular orientations across thin water films on sapphire

HYPOTHESIS

Water vapor binding to metal oxide surfaces produces thin water films with properties controlled by interactions with surface hydroxo sites. Hydrogen bonding populations vary across films and induce different molecular orientations than at the surface of liquid water. Identifying these differences can open possibilities for tailoring film-mediated catalytic reactions by choice of the supporting metal oxide substrate.

EXPERIMENTS

The (0001) face of a single sapphire (α-Al₂O₃) sample exposed to water vapor and the surface of liquid water were probed by polarization dependent Sum Frequency Generation-Vibration Spectroscopy (SFG-VS). Molecular dynamics (MD) provided insight into the hydrogen bond populations and molecular orientations across films and liquid water.

FINDINGS

SFG-VS revealed a submonolayer film on sapphire exposed to 43% relative humidity (R.H.), and a multilayer film at 78% R.H. Polarization dependent SFG-VS spectra showed that median tilt angles of free OH bonds on the top of films are at ∼43° from the normal of the (0001) face but at 38° on neat liquid water. These values align with MD simulations, which also show that up to 36% of all OH bonds on films are free. This offers new means for understanding how interfacial reactions on sapphire-supported water films could contrast with those involving liquid water.

The Surface Activity of the Hydrated Proton Is Substantially Higher than That of the Hydroxide Ion

The behavior of hydroxide and hydrated protons, the auto-ionization products of water, at surfaces is important for a wide range of applications and disciplines. However, it is unknown at which bulk concentration these ions start to become surface active at the water-air interface. Here, we report changes in the D2 O-air interface in the presence of excess D+ hyd /OD- hyd determined using surface-sensitive vibrational sum-frequency generation (SFG) spectroscopy. The onset of the perturbation of the D2 O surface occurs at a bulk concentration as low as 2.7±0.2 mm D+ hyd . In contrast, a concentration of several hundred mm OD- hyd is required to change the D2 O surface. The hydrated proton is thus orders of magnitude more surface-active than hydroxide at the water-air interface.

Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water-Air Interface

Density functional theory-based molecular dynamics simulations are increasingly being used for simulating aqueous interfaces. Nonetheless, the choice of the appropriate density functional, critically affecting the outcome of the simulation, has remained arbitrary. Here, we assess the performance of various exchange-correlation (XC) functionals, based on the metrics relevant to sum-frequency generation spectroscopy. The structure and dynamics of water at the water-air interface are governed by heterogeneous intermolecular interactions, thereby providing a critical benchmark for XC functionals. We find that the XC functionals constrained by exact functional conditions (revPBE and revPBE0) with the dispersion correction show excellent performance. The poor performance of the empirically optimized density functional (M06-L) indicates the importance of satisfying the exact functional condition. Understanding the performance of different XC functionals can aid in resolving the controversial interpretation of the interfacial water structure and direct the design of novel, improved XC functionals better suited to describing the heterogeneous interactions in condensed phases.

Structure and Dynamics of Water at the Water-Air Interface Using First-Principles Molecular Dynamics Simulations. II. NonLocal vs Empirical van der Waals Corrections

van der Waals (vdW) correction schemes have been recognized to be essential for an accurate description of liquid water in first-principles molecular dynamics simulation. The description of the structure and dynamics of water is governed by the type of the vdW corrections. So far, two vdW correction schemes have been often used: empirical vdW corrections and nonlocal vdW corrections. In this paper, we assess the influence of the empirical vs nonlocal vdW correction schemes on the structure and dynamics of water at the water–air interface. Since the structure of water at the water–air interface is established by a delicate balance of hydrogen bond formation and breaking, the simulation at the water–air interface provides a unique platform to testify as to the heterogeneous interaction of water. We used the metrics [Ohto et al. J. Chem. Theory Comput., 2019, 15, 595−60230468702] which are directly connected with the sum-frequency generation spectroscopic measurement. We find that the overall performance of nonlocal vdW methods is either similar or worse compared to the empirical vdW methods. We also investigated the performance of the optB88-DRSLL functional, which showed slightly less accuracy than the revPBE-D3 method. We conclude that the revPBE-D3 method shows the best performance for describing the interfacial water.

Deconstructing water's diffuse OH stretching vibrational spectrum with cold clusters

The diffuse vibrational envelope displayed by water precludes direct observation of how different hydrogen-bond topologies dictate the spectral response of individual hydroxy group (OH) oscillators. Using cold, isotopically labeled cluster ions, we report the spectral signatures of a single, intact water (H₂O) molecule embedded at various sites in the clathrate-like cage structure adopted by the Cs⁺·(D₂O)₂₀ ion. These patterns reveal the site-dependent correlation between the frequencies of the two OH groups on the same water molecule and establish that the bound OH companion of the free OH group exclusively accounts for bands in the lower-energy region of the spectrum. The observed multiplet structures reveal the homogeneous linewidths of the fundamentals and quantify the anharmonic contributions arising from coupling to both the intramolecular bending and intermolecular soft modes.