Orientational Distribution of Free O-H Groups of Interfacial Water is Exponential

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.

Experimental and simulation-based characterization of surfactant adsorption layers at fluid interfaces

Adsorption of surfactants to fluid interfaces occurs in numerous technological and daily-life contexts. The coverage at the interface and other properties of the formed adsorption layers determine the performance of a surfactant with regard to the desired application. Given the importance of these applications, there is a great demand for the comprehensive characterization and understanding of surfactant adsorption layers. In this review, we provide an overview of suitable experimental and simulation-based techniques and review the literature in which they were used for the investigation of surfactant adsorption layers. We come to the conclusion that, while these techniques have been successfully applied to investigate Langmuir monolayers of water-insoluble surfactants, their application to the study of Gibbs adsorption layers of water-soluble surfactants has not been fully exploited. Finally, we emphasize the great potential of these methods in providing a deeper understanding of the behavior of soluble surfactants at interfaces, which is crucial for optimizing their performance in various applications.

Sum rule comparison of narrowband and broadband sum frequency generation spectra and comparison with theory

Earlier sum frequency generation (SFG) experiments involve one infrared and one visible laser, and a measurement of the intensity of the response, yielding data on the surface sensitive properties of the sample. Recently, both the real and imaginary components of the susceptibility were measured in two different sets of experiments. In one set, a broadband infrared laser was used, permitting observations at very short times, while in another set the infrared laser was narrowband, permitting higher spectral resolution. The differences in the spectrum obtained by the two will be most evident in studying narrow absorption bands, e.g., the band due to dangling OH bonds at a water interface. The direct comparisons in the integrated amplitude (sum rule) of the imaginary part of the dangling OH bond region differ by a factor of 3. Due to variations in experimental setup and data processing, corrections were made for the quartz reference, Fresnel factors, and the incident visible laser wavelength. After the corrections, the agreement differs now by the factors of 1.1 within broadband and narrowband groups and the two groups now differ by a factor of 1.5. The 1.5 factor may arise from the extra heating of the more powerful broadband laser system on the water surface. The convolution from the narrowband SFG spectrum to the broadband SFG spectrum is also investigated and it does not affect the sum rule. Theory and narrowband experiments are compared using the sum rule and agree to a factor of 1.3 with no adjustable parameters.

Spontaneous Appearance of Triiodide Covering the Topmost Layer of the Iodide Solution Interface Without Photo-Oxidation

Ions containing iodine atoms at the vapor-aqueous solution interfaces critically affect aerosol growth and atmospheric chemistry due to their complex chemical nature and multivalency. While the surface propensity of iodide ions has been intensely discussed in the context of the Hofmeister series, the stability of various ions containing iodine atoms at the vapor-water interface has been debated. Here, we combine surface-specific sum-frequency generation (SFG) vibrational spectroscopy with ab initio molecular dynamics simulations to examine the extent to which iodide ions cover the aqueous surface. The SFG probe of the free O-D stretch mode of heavy water indicates that the free O-D group density decreases drastically at the interface when the bulk NaI concentration exceeds ∼2 M. The decrease in the free O-D group density is attributed to the spontaneous appearance of triiodide that covers the topmost interface rather than to the surface adsorption of iodide. This finding demonstrates that iodide is not surface-active, yet the highly surface-active triiodide is generated spontaneously at the water-air interface, even under dark and oxygen-free conditions. Our study provides an important first step toward clarifying iodine chemistry and pathways for aerosol formation.

Appraisal of TIP4P-type models at water surface

In view of the current situation in which non-polarizable rigid water models have been scarcely examined against surface-specific properties, we appraise TIP4P-type models at the liquid water surface on the basis of heterodyne-detected sum frequency generation (HD-SFG) spectroscopy. We find in the HD-SFG spectrum of the water surface that the peak frequency of the hydrogen-bonded OH band, the half width at half maximum of the hydrogen-bonded OH band, and the full width at half maximum of the free OH band are best reproduced by TIP4P, TIP4P/Ew, and TIP4P/Ice, respectively, whereas it is already well known that TIP4P/2005 best reproduces the surface tension. These TIP4P-type models perform better at the water surface in terms of the present appraisal items than some polarizable models in the literature.

Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models

In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.

Structural evolution of water-in-propylene carbonate mixtures revealed by polarized Raman spectroscopy and molecular dynamics

The liquid structure of systems wherein water is limited in concentration or through geometry is of great interest in various fields such as biology, materials science, and electrochemistry. Here, we present a combined polarized Raman and molecular dynamics investigation of the structural changes that occur as water is added incrementally to propylene carbonate (PC), a polar, aprotic solvent that is important in lithium-ion batteries. Polarized Raman spectra of PC solutions were collected for water mole fractions 0.003 ≤ χwater ≤ 0.296, which encompasses the solubility range of water in PC. The novel approach taken herein provides additional hydrogen bond and solvation characterization of this system that has not been achievable in previous studies. Analysis of the polarized carbonyl Raman band in conjunction with simulations demonstrated that the bulk structure of the solvent remained unperturbed upon the addition of water. Experimental spectra in the O-H stretching region were decomposed through Gaussian fitting into sub-bands and comparison to studies of dilute HOD in D2O. With the aid of simulations, we identified these different bands as water arrangements having different degrees of hydrogen bonding. The observed water structure within PC indicates that water tends to self-aggregate, forming a hydrogen bond network that is distinctly different from the bulk and dependent on concentration. For example, at moderate concentrations, the most likely aggregate structures are chains of water molecules, each with two hydrogen bonds.

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

Experimental and Thermodynamic Viewpoints on Claims of a Spontaneous H2O2 Formation at the Air-Water Interface

Recent claims of the spontaneous H2O2 formation at the air-water interface of water microdroplets have sparked debates on its feasibility. New results from different research groups have provided more insight into these claims, but conclusive proofs are still far from realized. In this Perspective, thermodynamic viewpoints, potential experiments, and theoretical approaches are presented as references for future studies. We suggest that future work should seek for H2 byproduct as indirect evidence to confirm the feasibility of this phenomenon. Examining potential energy surfaces for H2O2 formation reaction when moving from the bulk to the interface under the influence of the local electric fields is also critical to establish this phenomenon.

Deciphering the Properties of Nanoconfined Aqueous Solutions by Vibrational Sum Frequency Generation Spectroscopy

When confined between walls at nanometer distances, water exhibits surprisingly different properties with reference to bare interfacial water. Based on computer simulations, we demonstrate how vibrational sum frequency generation (VSFG) spectroscopy can be used-even with very mild symmetry breaking-to discriminate multilayer water in wide slit pores from both bilayer and monolayer water confined within molecularly narrow pores. Applying the technique, the VSFG lineshapes of monolayer, bilayer, and multilayer water are found to differ in characteristic ways, which is explained by their distinct density stratifications giving rise to different H-bonding patterns in the respective solvation layers.

On the Propensity of Excess Hydroxide Ions at the Alcohol Monolayer-Water Interface

Atomistic molecular dynamics simulations have been employed to study the self-ion (H⁺ and OH^-) distribution at the interface between long-chain C₁₆-OH alcohol (cetyl alcohol) monolayer and water. It is well known that the free air-water interface is acidic due to accumulation of the hydronium (H₃O⁺) ions at the interface. In the present study, we have observed that contrary to the air-water interface, at the long-chain alcohol monolayer-water interface, it is the hydroxide (OH^-) ion, not the hydronium ion (H₃O⁺) that gets accumulated. By calculating the potential of mean forces, it is confirmed that there is extra stabilization for the OH^- ions at the interface relative to the bulk, but no such stabilization is observed for the H₃O⁺ ions. By analyzing the interaction of the self-ions with other constituents in the medium, it is clearly shown that the favorable interaction of the OH^- ions with the alcoholic -OH groups stabilizes this ion at the interface. By calculating coordination numbers of the self-ions it is observed that around 50% water neighbors are substituted by alcoholic -OH in case of the hydroxide ion at the interface, whereas in the case of hydronium ions, only 15% water neighbors are substituted by the alcoholic -OH. The most interesting observation about the local structure and H-bonding pattern is that the hydroxide ion acts solely as the H-bond acceptor, but the hydronium ion acts only as the H-bond donor.

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.

Hydration Structure of Na+ and K+ Ions in Solution Predicted by Data-Driven Many-Body Potentials

The hydration structure of Na⁺ and K⁺ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na⁺ and K⁺ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

Investigation of Water Evaporation Process at Air/Water Interface using Hofmeister Ions

Evaporation is an interfacial phenomenon in which a water molecule breaks the intermolecular hydrogen (H-) bonds and enters the vapor phase. However, a detailed demonstration of the role of interfacial water structure in the evaporation process is still lacking. Here, we purposefully perturb the H-bonding environment at the air/water interface by introducing kosmotropic (HPO₄^-2, SO₄^-2, and CO₃^-2) and chaotropic ions (NO₃^- and I^-) to determine their influence on the evaporation process. Using time-resolved interferometry on aqueous salt droplets, we found that kosmotropes reduce evaporation, whereas chaotropes accelerate the evaporation process, following the Hofmeister series: HPO₄^-2 < SO₄^-2 < CO₃^-2 < Cl^- < NO₃^- < I^-. To extract deeper molecular-level insights into the observed Hofmeister trend in the evaporation rates, we investigated the air/water interface in the presence of ions using surface-specific sum frequency generation (SFG) vibrational spectroscopy. The SFG vibrational spectra reveal the significant impact of ions on the strength of the H-bonding environment and the orientation of free OH oscillators from ∼36.2 to 48.4° at the air/water interface, where both the effects follow the Hofmeister series. It is established that the slow evaporating water molecules experience a strong H-bonding environment with free OH oscillators tilted away from the surface normal in the presence of kosmotropes. In contrast, the fast evaporating water molecules experience a weak H-bonding environment with free OH oscillators tilted toward the surface normal in the presence of chaotropes at the air/water interface. Our experimental outcomes showcase the complex bonding environment of interfacial water molecules and their decisive role in the evaporation process.

The dielectric function profile across the water interface through surface-specific vibrational spectroscopy and simulations

The dielectric properties of interfacial water on subnanometer length scales govern chemical reactions, carrier transfer, and ion transport at interfaces. Yet, the nature of the interfacial dielectric function has remained under debate as it is challenging to access the interfacial dielectric with subnanometer resolution. Here we use the vibrational response of interfacial water molecules probed using surface-specific sum-frequency generation (SFG) spectra to obtain exquisite depth resolution. Different responses originate from water molecules at different depths and report back on the local interfacial dielectric environment via their spectral amplitudes. From experimental and simulated SFG spectra at the air/water interface, we find that the interfacial dielectric constant changes drastically across an ∼1 Å thin interfacial water region. The strong gradient of the interfacial dielectric constant leads, at charged planar interfaces, to the formation of an electric triple layer that goes beyond the standard double-layer model.

Polarization-Dependent Heterodyne-Detected Sum-Frequency Generation Spectroscopy as a Tool to Explore Surface Molecular Orientation and Ångström-Scale Depth Profiling

Sum-frequency generation (SFG) spectroscopy provides a unique optical probe for interfacial molecules with interface-specificity and molecular specificity. SFG measurements can be further carried out at different polarization combinations, but the target of the polarization-dependent SFG is conventionally limited to investigating the molecular orientation. Here, we explore the possibility of polarization-dependent SFG (PD-SFG) measurements with heterodyne detection (HD-PD-SFG). We stress that HD-PD-SFG enables accurate determination of the peak amplitude, a key factor of the PD-SFG data. Subsequently, we outline that HD-PD-SFG can be used not only for estimating the molecular orientation but also for investigating the interfacial dielectric profile and studying the depth profile of molecules. We further illustrate the variety of combined simulation and PD-SFG studies.

Local Ice-like Structure at the Liquid Water Surface

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

Accurate molecular orientation at interfaces determined by multimode polarization-dependent heterodyne-detected sum-frequency generation spectroscopy via multidimensional orientational distribution function

Many essential processes occur at soft interfaces, from chemical reactions on aqueous aerosols in the atmosphere to biochemical recognition and binding at the surface of cell membranes. The spatial arrangement of molecules specifically at these interfaces is crucial for many of such processes. The accurate determination of the interfacial molecular orientation has been challenging due to the low number of molecules at interfaces and the ambiguity of their orientational distribution. Here, we combine phase- and polarization-resolved sum-frequency generation (SFG) spectroscopy to obtain the molecular orientation at the interface. We extend an exponentially decaying orientational distribution to multiple dimensions, which, in conjunction with multiple SFG datasets obtained from the different vibrational modes, allows us to determine the molecular orientation. We apply this new approach to formic acid molecules at the air–water interface. The inferred orientation of formic acid agrees very well with ab initio molecular dynamics data. The phase-resolved SFG multimode analysis scheme using the multidimensional orientational distribution thus provides a universal approach for obtaining the interfacial molecular orientation.

Ion adsorption stabilizes bulk nanobubbles

The mechanism leading to the extraordinary stability of bulk nanobubbles in aqueous solutions remains an outstanding problem in soft matter, modern surface science, and physical chemistry science. In this work, the stability of bulk nanobubbles in electrolyte solutions under different pH levels and ionic strengths is studied. Nanobubbles are generated via the technique of ultrasonic cavitation, and characterized for size, number concentration and zeta potential under ambient conditions. Experimental results show that nanobubbles can survive in both acidic and basic solutions with pH values far away from the isoelectric point. We attribute the enhanced stability with increasing acidity or alkalinity of the aqueous solutions to the effective accumulation of net charges, regardless of their sign. The kinetic stability of the nanobubbles in various aqueous solutions is evaluated within the classic DLVO framework. Further, by combining a modified Poisson-Boltzmann equation with a modified Langmuir adsorption model, we describe a simple model that captures the influence of ion species and bulk concentration and reproduce the dependence of the nanobubble's surface potential on pH. We also discuss the apparent contradiction between quantitative calculation by ion stabilization model and experimental results. This essentially requires insight into the structure and dynamics of interfacial water on the atomic-scale.

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

Data-driven many-body models enable a quantitative description of chloride hydration from clusters to bulk

We present a new data-driven potential energy function (PEF) describing chloride-water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride-water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride-water and water-water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments.

Ordered Water Layer on the Macroscopically Hydrophobic Fluorinated Polymer Surface and Its Ultrafast Vibrational Dynamics

Hydrophobic-like water monolayers have been predicted at the metal and some polar surfaces by theoretical simulations. However, direct experimental evidence for the presence of this water layer at surfaces, particularly at biomolecule and polymer surfaces, is yet to be validated at room temperature. Here we observe experimentally that an ordered molecular water layer is present at the hydrophobic fluorinated polymer such as polytetrafluoroethylene (PTFE) surface by using sum frequency generation vibrational spectroscopy. The macroscopic hydrophobicity of PTFE surface is actually hydrophilic at the molecular level. The macroscopically hydrophobic character of PTFE is indeed resulting from the hydrophobicity of the ordered two-dimension (2D) water layer, in which cyclic water tetramer structure is found. The water layer at humidity of ≤40% has a vibrational relaxation time of 550 ± 60 fs. The vibrational relaxation time in the frequency range of 3200-3400 cm^-1 shows remarkable difference from the interfacial water at the air/H₂O interface and the lipid/H₂O interface. No discernible frequency dependence of the vibrational relaxation time is observed, indicating the homogeneous dynamics of OH groups in the water layer. These insights into the water layer at the macroscopically hydrophobic surface may contribute to a better understanding of the hydrophobic interaction and interfacial water dynamics.

Assessing the Impact of Solvent Selection on Vibrational Sum-Frequency Scattering Spectroscopy Experiments

The development of vibrational sum-frequency scattering (S-VSF) spectroscopy has opened the door to directly probing nanoparticle surfaces with an interfacial and chemical specificity that was previously reserved for planar interfacial systems. Despite its potential, challenges remain in the application of S-VSF spectroscopy beyond simplified chemical systems. One such challenge includes infrared absorption by an absorptive continuous phase, which will alter the spectral lineshapes within S-VSF spectra. In this study, we investigate how solvent vibrational modes manifest in S-VSF spectra of surfactant stabilized nanoemulsions and demonstrate how corrections for infrared absorption can recover the spectral features of interfacial solvent molecules. We also investigate infrared absorption for systems with the absorptive phase dispersed in a nonabsorptive continuous phase to show that infrared absorption, while reduced, will still impact the S-VSF spectra. These studies are then used to provide practical recommendations for anyone wishing to use S-VSF to study nanoparticle surfaces where absorptive solvents are present.

Reorientation-induced relaxation of free OH at the air/water interface revealed by ultrafast heterodyne-detected nonlinear spectroscopy

The uniqueness of water originates from its three-dimensional hydrogen-bond network, but this hydrogen-bond network is suddenly truncated at the interface and non-hydrogen-bonded OH (free OH) appears. Although this free OH is the most characteristic feature of interfacial water, the molecular-level understanding of its dynamic property is still limited due to the technical difficulty. We study ultrafast vibrational relaxation dynamics of the free OH at the air/water interface using time-resolved heterodyne-detected vibrational sum frequency generation (TR-HD-VSFG) spectroscopy. With the use of singular value decomposition (SVD) analysis, the vibrational relaxation (T₁) times of the free OH at the neat H₂O and isotopically-diluted water interfaces are determined to be 0.87 ± 0.06 ps (neat H₂O), 0.84 ± 0.09 ps (H₂O/HOD/D₂O = 1/2/1), and 0.88 ± 0.16 ps (H₂O/HOD/D₂O = 1/8/16). The absence of the isotope effect on the T₁ time indicates that the main mechanism of the vibrational relaxation of the free OH is reorientation of the topmost water molecules. The determined sub-picosecond T₁ time also suggests that the free OH reorients diffusively without the switching of the hydrogen-bond partner by the topmost water molecule.

Water’s hydrogen-bond network is truncated at hydrophobic interfaces and the dynamics of the resulting free OH groups is not well understood. The authors experimentally show that the main vibrational relaxation mechanism for free OH at the air-water interface is a diffusive molecular reorientation, rather than intramolecular energy transfer.

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.

Charge transfer as a ubiquitous mechanism in determining the negative charge at hydrophobic interfaces

The origin of the apparent negative charge at hydrophobic–water interfaces has fueled debates in the physical chemistry community for decades. The most common interpretation given to explain this observation is that negatively charged hydroxide ions (OH^–) bind strongly to the interfaces. Using first principles calculations of extended air–water and oil–water interfaces, we unravel a mechanism that does not require the presence of OH^–. Small amounts of charge transfer along hydrogen bonds and asymmetries in the hydrogen bond network due to topological defects can lead to the accumulation of negative surface charge at both interfaces. For water near oil, some spillage of electron density into the oil phase is also observed. The computed surface charge densities at both interfaces is approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.015\ {\rm{e}}/{{\rm{nm}}}^{2}$$\end{document}−0.015e∕nm2 in agreement with electrophoretic experiments. We also show, using an energy decomposition analysis, that the electronic origin of this phenomena is rooted in a collective polarization/charge transfer effect.

The accumulation of negative charge at hydrophobic–water interfaces has been a source of debate for a long time. Here the authors use ab initio calculations to show that the charge accumulation at air–water and oil–water interfaces is caused by subtle charge transfer processes.

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.

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.

Hydrogen bonding structure of confined water templated by a metal-organic framework with open metal sites

Water in confinement exhibits properties significantly different from bulk water due to frustration in the hydrogen-bond network induced by interactions with the substrate. Here, we combine infrared spectroscopy and many-body molecular dynamics simulations to probe the structure and dynamics of confined water as a function of relative humidity within a metal-organic framework containing cylindrical pores lined with ordered cobalt open coordination sites. Building upon the agreement between experimental and theoretical spectra, we demonstrate that water at low relative humidity binds initially to open metal sites and subsequently forms disconnected one-dimensional chains of hydrogen-bonded water molecules bridging between cobalt atoms. With increasing relative humidity, these water chains nucleate pore filling, and water molecules occupy the entire pore interior before the relative humidity reaches 30%. Systematic analysis of rotational and translational dynamics indicates heterogeneity in this pore-confined water, with water molecules displaying variable mobility as a function of distance from the interface.

The properties of water under confinement are significantly altered with respect to the bulk phase. Here the authors use infrared spectroscopy and many-body molecular dynamics simulations to show the structure and dynamics of confined water as a function of relative humidity within a metal-organic framework.

Time-dependent vibrational sum-frequency generation spectroscopy of the air-water interface

Vibrational sum-frequency generation spectroscopy is a powerful method to study the microscopic structure and dynamics of interfacial systems. Here we demonstrate a simple computational approach to calculate the time-dependent, frequency-resolved vibrational sum-frequency generation spectrum (TD-vSFG) of the air-water interface. Using this approach, we show that at the air-water interface, the transition of water molecules with bonded OH modes to free OH modes occurs at a time scale of $$\sim$$ ~ 3 ps, whereas water molecules with free OH modes rapidly make a transition to a hydrogen-bonded state within $$\sim$$ ~ 2 ps. Furthermore, we also elucidate the origin of the observed differential dynamics based on the time-dependent evolution of water molecules in the different local solvent environments.

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.

Vibrational spectra of dissociatively adsorbed D2O on Al-terminated α-Al2O3(0001) surfaces from ab initio molecular dynamics

Water can adsorb molecularly or dissociatively onto different sites of metal oxide surfaces. These adsorption sites can be disentangled using surface-sensitive vibrational spectroscopy. Here, we model Vibrational Sum Frequency (VSF) spectra for various forms of dissociated, deuterated water on a reconstructed, Al-terminated α-Al₂O₃(0001) surface at submonolayer coverages (the so-called 1-2, 1-4, and 1-4' modes). Using an efficient scheme based on velocity-velocity autocorrelation functions, we go beyond previous normal mode analyses by including anharmonicity, mode coupling, and thermal surface motion in the framework of ab initio molecular dynamics. In this way, we calculate vibrational density of states curves, infrared, and VSF spectra. Comparing computed VSF spectra with measured ones, we find that relative frequencies of resonances are in quite good agreement and linewidths are reasonably well represented, while VSF intensities coincide not well. We argue that intensities are sensitively affected by local interactions and thermal fluctuations, even at such low coverage, while absolute peak positions strongly depend on the choice of the electronic structure method and on the appropriate inclusion of anharmonicity.

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.

The Influence of Sodium Iodide Salt on the Interfacial Properties of Aqueous Methanol Solution by a Combined Molecular Simulation and Sum Frequency Generation Vibrational Spectroscopy Study

Understanding the influence of salt ions on the microscopic properties of liquid interfaces is of both fundamental and practical importance. A large number of previous experimental and theoretical investigations have explored the salt effects on the surfaces of either pure water or neat organic liquid. However, how the salt ions affect the interfacial structures of water/organic liquid mixtures has rarely been studied. Here, the molecular dynamics (MD) simulations and sum frequency generation vibrational spectroscopy (SFG-VS) were carried out to investigate the influence of sodium iodide (NaI) on the air/liquid interfaces of the methanol-water mixtures. The SFG-VS spectral intensities were discovered to increase with the addition of 3 M NaI, while the center frequencies of the C-H stretching vibrations at high methanol concentrations showed a ∼2 cm^-1 blue shift compared with those obtained before adding NaI. The MD results indicated that Na⁺ and I^- can only affect Part I (near the bulk phase) but not Part II (near the gas phase) of the interfacial region. It was also found that the average orientations of interfacial methyl groups were constant and not effectively disturbed by the changes of methanol concentrations or the addition of NaI. It is therefore concluded that the changes of the SFG-VS intensities upon the addition of NaI salts were mainly caused by the increasing number of interfacial methanol molecules. Further analysis showed that the existence of NaI affects the surface tensions more for the interfaces with higher bulk methanol concentrations, which is in agreement with the SFG-VS results. It is noteworthy that the maximum number density of methanol molecules with the net nonzero orientations is reached near the Gibbs dividing surface, the reasons of which are worth further investigating.