The surface roughness, but not the water molecular orientation varies with temperature at the water-air interface

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.

The Role of Sum-Frequency Generation Spectroscopy in Understanding On-Surface Reactions and Dynamics in Atmospheric Model-Systems

Surfaces, both water/air and solid/water, play an important role in mediating a multitude of processes central to atmospheric chemistry, particularly in the aerosol phase. However, the study of both static and dynamic properties of surfaces is highly challenging from an experimental standpoint, leading to a lack of molecular level information about the processes that take place at these systems and how they differ from bulk. One of the few techniques that has been able to capture ultrafast surface phenomena is time-resolved sum-frequency generation (SFG) spectroscopy. Since it is both surface-specific and chemically sensitive, the extension of this spectroscopic technique to the time domain makes it possible to study dynamic processes on the femtosecond time scale. In this Perspective, we will explore recent advances made in the field both in terms of studying energy dissipation as well as chemical reactions and the role the surface geometry plays in these processes.

Theoretical Study of the Two-Dimensional Vibrational Sum Frequency Generation Spectroscopy of the Air-Water Interface at Varying Temperature and Its Connections to the Interfacial Structure and Dynamics

We performed a theoretical study of the temperature variation of two-dimensional vibrational sum frequency generation (2D-VSFG) spectra of the OH stretch modes at air-water interfaces in the mid-IR region. The calculations are performed at four different temperatures from 250 to 325 K by using a combination of techniques involving response function formalism of nonlinear spectroscopy, electronic structure calculations, and molecular dynamics simulations. Also, the calculations are performed for isotopically dilute solutions so that the intra- and intermolecular coupling between the vibrational modes of interest can be ignored. We have established the connections of temperature variation of various frequency- and time-dependent features of the calculated spectra to the changes in the underlying structure and dynamics of the interfaces. The results reveal that interfacial water is dynamically more heterogeneous than bulk water, with three dominant dynamical processes exhibiting their corresponding time-dependent features in the 2D-VSFG spectrum. These are the spectral diffusion of hydrogen-bonded OH groups at the interface, conversion of an initially hydrogen-bonded OH group to a dangling OH which is a stable state for surface water, unlike the bulk water, and the third one, which involves the conversion of an initially free or dangling OH group to its hydrogen-bonded state at the interface. The temporal appearance of the cross peaks corresponding to interconversion of the hydrogen-bonded state to the dangling state or vice versa of an interfacial OH group is found to take place at a slower rate than the dynamics of spectral diffusion of hydrogen-bonded molecules at the interface, which, in turn, is slower than the corresponding spectral diffusion of bulk water molecules. The temperature variation of these dynamic processes can be linked to the decay of appropriate hydrogen-bond and non-hydrogen-bond time correlation functions of interfacial water molecules for the different air-water systems studied in this work.

Microscopic properties of forces from ice solidification interface acting on silica surfaces based on molecular dynamics simulations

The origin of the forces acting on a silica surface from an ice solidification interface was investigated to understand the solidification phenomenon and its impact on nanometer-scale structures using molecular dynamics simulations. The microscopic forces were determined by appropriately averaging the forces acting on the silica wall from the water molecules in time and space; the time evolutions of these microscopic forces during the solidification processes were investigated for three types of silica surfaces. The results indicate that the microscopic forces fluctuate more after the solidification interface makes contact with the wall surface. To visualize the changes in the microscopic forces and hydrogen bonds due to solidification, their differences compared to the liquid state were calculated. When the solidification interface is near the wall, the changes in these microscopic forces and hydrogen bonds due to solidification are correlated. This tendency is more significant for an amorphous wall and a wall with a structure than for a crystalline wall. The changes in the microscopic force depend on the water molecules that behave as acceptors when forming the hydrogen bonds with the wall and on the configuration of the silanol groups on the silica surfaces.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.

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

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

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.

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.

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.

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.

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.

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.

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.

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 Surface of Ice under Equilibrium and Nonequilibrium Conditions

The ice premelt, often called the quasi-liquid layer (QLL), is key for the lubrication of ice, gas uptake by ice, and growth of aerosols. Despite its apparent importance, in-depth understanding of the ice premelt from the microscopic to the macroscopic scale has not been gained. By reviewing data obtained using molecular dynamics (MD) simulations, sum-frequency generation (SFG) spectroscopy, and laser confocal differential interference contrast microscopy (LCM-DIM), we provide a unified view of the experimentally observed variation in quasi-liquid (QL) states. In particular, we disentangle three distinct types of QL states of disordered layers, QL-droplet, and QL-film and discuss their nature.

The topmost ice layer is energetically unstable, as the topmost interfacial H₂O molecules lose a hydrogen bonding partner, generating a disordered layer at the ice–air interface. This disordered layer is homogeneously distributed over the ice surface. The nature of the disordered layer changes over a wide temperature range from −90 °C to the bulk melting point. Combined MD simulations and SFG measurements reveal that the topmost ice surface starts to be disordered around −90 °C through a process that the topmost water molecules with three hydrogen bonds convert to a doubly hydrogen-bonded species. When the temperature is further increased, the second layer starts to become disordered at around −16 °C. This disordering occurs not in a gradual manner, but in a bilayer-by-bilayer manner.

When the temperature reaches −2 °C, more complicated structures, QL-droplet and QL-film, emerge on the top of the ice surface. These QL-droplets and QL-films are inhomogeneously distributed, in contrast to the disordered layer. We show that these QL-droplet and QL-film emerge only under supersaturated/undersaturated vapor pressure conditions, as partial and pseudopartial wetting states, respectively. Experiments with precisely controlled pressure show that, near the water vapor pressure at the vapor-ice equilibrium condition, no QL-droplet and QL-film can be observed, implying that the QL-droplet and QL-film emerge exclusively under nonequilibrium conditions, as opposed to the disordered layers formed under equilibrium conditions.

These findings are connected with many phenomena related to the ice surface. For example, we explain how the disordering of the topmost ice surface governs the slipperiness of the ice surface, allowing for ice skating. Further focus is on the gas uptake mechanism on the ice surface. Finally, we note the unresolved questions and future challenges regarding the ice premelt.

Evaporation of nanoscale water on a uniformly complete wetting surface at different temperatures

The evaporation of nanoscale water films on surfaces affects many processes in nature and industry. Using molecular dynamics (MD) simulations, we show the evaporation of a nanoscale water film on a uniformly complete wetting surface at different temperatures. With the increase in temperature, the growth of the water evaporation rate becomes slow. Analyses show that the hydrogen bond (H-bond) lifetimes and orientational autocorrelation times of the outermost water film decrease slowly with the increase in temperature. Compared to a thicker water film, the H-bond lifetimes and orientational autocorrelation times of a monolayer water film are much slower. This suggests that the lower evaporation rate of the monolayer water film on a uniformly complete wetting surface may be caused by the constriction of the water rotation due to the substrate. This finding may be helpful for controlling nanoscale water evaporation within a certain range of temperatures.

Water at surfaces with tunable surface chemistries

Aqueous interfaces are ubiquitous in natural environments, spanning atmospheric, geological, oceanographic, and biological systems, as well as in technical applications, such as fuel cells and membrane filtration. Where liquid water terminates at a surface, an interfacial region is formed, which exhibits distinct properties from the bulk aqueous phase. The unique properties of water are governed by the hydrogen-bonded network. The chemical and physical properties of the surface dictate the boundary conditions of the bulk hydrogen-bonded network and thus the interfacial properties of the water and any molecules in that region. Understanding the properties of interfacial water requires systematically characterizing the structure and dynamics of interfacial water as a function of the surface chemistry. In this review, we focus on the use of experimental surface-specific spectroscopic methods to understand the properties of interfacial water as a function of surface chemistry. Investigations of the air-water interface, as well as efforts in tuning the properties of the air-water interface by adding solutes or surfactants, are briefly discussed. Buried aqueous interfaces can be accessed with careful selection of spectroscopic technique and sample configuration, further expanding the range of chemical environments that can be probed, including solid inorganic materials, polymers, and water immiscible liquids. Solid substrates can be finely tuned by functionalization with self-assembled monolayers, polymers, or biomolecules. These variables provide a platform for systematically tuning the chemical nature of the interface and examining the resulting water structure. Finally, time-resolved methods to probe the dynamics of interfacial water are briefly summarized before discussing the current status and future directions in studying the structure and dynamics of interfacial water.

Reduced Near-Resonant Vibrational Coupling at the Surfaces of Liquid Water and Ice

We study the resonant interaction of the OH stretch vibrations of water molecules at the surfaces of liquid water and ice using heterodyne-detected sum-frequency generation (HD-SFG) spectroscopy. By studying different isotopic mixtures of H₂O and D₂O, we vary the strength of the interaction, and we monitor the resulting effect on the HD-SFG spectrum of the OH stretch vibrations. We observe that the near-resonant coupling effects are weaker at the surface than in the bulk, for both water and ice, indicating that for both phases of water the OH vibrations are less strongly delocalized at the surface than in the bulk.

Second-order spectral lineshapes from charged interfaces

Second-order nonlinear spectroscopy has proven to be a powerful tool in elucidating key chemical and structural characteristics at a variety of interfaces. However, the presence of interfacial potentials may lead to complications regarding the interpretation of second harmonic and vibrational sum frequency generation responses from charged interfaces due to mixing of absorptive and dispersive contributions. Here, we examine by means of mathematical modeling how this interaction influences second-order spectral lineshapes. We discuss our findings in the context of reported nonlinear optical spectra obtained from charged water/air and solid/liquid interfaces and demonstrate the importance of accounting for the interfacial potential-dependent χ⁽³⁾ term in interpreting lineshapes when seeking molecular information from charged interfaces using second-order spectroscopy.

Charged interfaces are important in chemical systems, but the influence of charge on vibrational sum frequency spectra has only recently been considered. Here the authors show the importance of accounting for the interfacial potential-dependent χ⁽³⁾ term in interpreting spectral lineshapes from charged interfaces.

Excess Hydrogen Bond at the Ice-Vapor Interface around 200 K

Phase-resolved sum-frequency generation measurements combined with molecular dynamics simulations are employed to study the effect of temperature on the molecular arrangement of water on the basal face of ice. The topmost monolayer, interrogated through its nonhydrogen-bonded, free O-H stretch peak, exhibits a maximum in surface H-bond density around 200 K. This maximum results from two competing effects: above 200 K, thermal fluctuations cause the breaking of H bonds; below 200 K, the formation of bulklike crystalline interfacial structures leads to H-bond breaking. Knowledge of the surface structure of ice is critical for understanding reactions occurring on ice surfaces and ice nucleation.

Observation and Identification of a New OH Stretch Vibrational Band at the Surface of Ice

We study the signatures of the OH stretch vibrations at the basal surface of ice using heterodyne-detected sum-frequency generation and molecular dynamics simulations. At 150 K, we observe seven distinct modes in the sum-frequency response, five of which have an analogue in the bulk, and two pure surface-specific modes at higher frequencies (∼3530 and ∼3700 cm^-1). The band at ∼3530 cm^-1 has not been reported previously. Using molecular dynamics simulations, we find that the ∼3530 cm^-1 band contains contributions from OH stretch vibrations of both fully coordinated interfacial water molecules and water molecules with two donor and one acceptor hydrogen bond.

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

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

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

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

Vibrational coupling is relevant not only for dissipation of excess energy after chemical reactions but also for elucidating molecular structure and dynamics. It is particularly important for O-H stretch vibrational spectra of water, for which it is known that in bulk both intra- and intermolecular coupling alter the intensity and line shape of the spectra. In contrast with bulk, the unified picture of the inter/intra-molecular coupling of O-H groups at the water-air interface has been lacking. Here, combining sum-frequency generation experiments and simulation for isotopically diluted water and alcohols, we unveil effects of inter- and intramolecular coupling on the vibrational spectra of interfacial water. Our results show that both inter- and intramolecular coupling contribute to the O-H stretch vibrational response of the neat H₂O surface, with intramolecular coupling generating a double-peak feature, while the intermolecular coupling induces a significant red shift in the O-H stretch response.

Dissecting the Molecular Structure of the Air/Water Interface from Quantum Simulations of the Sum-Frequency Generation Spectrum

The molecular characterization of the air/water interface is a key step in understanding fundamental multiphase phenomena ranging from heterogeneous chemical processes in the atmosphere to the hydration of biomolecules. The apparent simplicity of the air/water interface, however, masks an underlying complexity associated with the dynamic nature of the water hydrogen-bond network that has so far hindered an unambiguous characterization of its microscopic properties. Here, we demonstrate that the application of quantum many-body molecular dynamics, which enables spectroscopically accurate simulations of water from the gas to the condensed phase, leads to a definitive molecular-level picture of the interface region. For the first time, excellent agreement is obtained between the simulated vibrational sum-frequency generation spectrum and the most recent state-of-the-art measurements, without requiring any empirical frequency shift or ad hoc scaling of the spectral intensity. A systematic dissection of the spectral features demonstrates that a rigorous representation of nuclear quantum effects as well as of many-body energy and electrostatic contributions is necessary for a quantitative reproduction of the experimental data. The unprecedented accuracy of the simulations presented here indicates that quantum many-body molecular dynamics can enable predictive studies of aqueous interfaces, which by complementing analogous experimental measurements will provide unique molecular insights into multiphase and heterogeneous processes of relevance in chemistry, biology, materials science, and environmental research.

Molecular Mechanism of Water Evaporation

Evaporation is the process by which water changes from a liquid to a gas or vapor, and is a key step in Earth's water cycle. At the molecular level, evaporation requires breaking at least one very strong intermolecular bond between two water molecules at the interface. Despite the importance of this process the molecular mechanism by which an evaporating water molecule gains sufficient energy to escape from the surface has remained elusive. Here, we show, using molecular dynamics simulations at the water-air interface with polarizable classical force field models, that the high kinetic energy of the evaporated water molecule is enabled by a well-timed making and breaking of hydrogen bonds involving at least three water molecules at the interface, the recoil of which allows one of the molecules to escape. The evaporation of water is thus enabled by concerted, ultrafast hydrogen-bond dynamics of interfacial water, and follows one specific molecular pathway.

Lipid carbonyl groups terminate the hydrogen bond network of membrane-bound water

We present a combined experimental sum-frequency generation (SFG) spectroscopy and ab initio molecular dynamics simulations study to clarify the structure and orientation of water at zwitterionic phosphatidylcholine (PC) lipid and amine N-oxide (AO) surfactant monolayers. Simulated O-H stretch SFG spectra of water show good agreement with the experimental data. The SFG response at the PC interface exhibits positive peaks, whereas both negative and positive bands are present for the similar zwitterionic AO interface. The positive peaks at the water/PC interface are attributed to water interacting with the lipid carbonyl groups, which act as efficient hydrogen bond acceptors. This allows the water hydrogen bond network to reach, with its (up-oriented) O-H groups, into the headgroup of the lipid, a mechanism not available for water underneath the AO surfactant. This highlights the role of the lipid carbonyl group in the interfacial water structure at the membrane interface, namely, stabilizing the water hydrogen bond network.