Hydrogen-Bonding Interactions at the Vapor/Water Interface Investigated by Vibrational Sum-Frequency Spectroscopy of HOD/H2O/D2O Mixtures and Molecular Dynamics Simulations

Revisiting the Thickness of the Air-Water Interface from Two Extremes of Interface Hydrogen Bond Dynamics

The air-water interface plays a crucial role in many aspects of science because of its unique properties, such as a two-dimensional hydrogen bond (HB) network and completely different HB dynamics compared to bulk water. However, accurately determining the boundary of interfacial and bulk water, that is, the thickness of the air-water interface, still challenges experimentalists. Various simulation-based methods have been developed to estimate the thickness, converging on a range of approximately 3-10 (Å). In this study, we introduce a novel approach, grounded in density functional theory-based molecular dynamics and deep potential molecular dynamics simulations, to measure the air-water interface thickness, offering a different perspective based on prior research. To capture realistic HB dynamics in the air-water interface, two extreme scenarios of the interface HB dynamics are obtained: one underestimates the interface HB dynamics, while the other overestimates it. Surprisingly, our results suggest that the interface HB dynamics in both scenarios converges as the thickness of the air-water interface increases to 4 (Å). This convergence point, indicative of the realistic interface thickness, is also validated by our calculation of anisotropic decay of OH stretch and the free OH dynamics at the air-water interface.

A generalized van't Hoff relation for the temperature dependence of complex-valued nonlinear spectra

The temperature dependence of spectra can reveal important insights into the structural and dynamical behavior of the system being probed. In the case of linear spectra, this has been exploited to investigate the thermodynamic driving forces governing the spectral response. Indeed, the temperature derivative of a spectrum can be used to obtain effective energetic and entropic profiles as a function of the measured frequency. The former can further be used to predict the temperature-dependent spectrum via a van't Hoff relation. However, these approaches are not directly applicable to nonlinear, complex-valued spectra, such as vibrational sum-frequency generation (SFG) or two-dimensional infrared (2D-IR) photon echo spectra. Here, we show how the energetic and entropic driving forces governing such nonlinear spectra can be determined and used within a generalized van't Hoff relation to predict their temperature dependence. The central idea is to allow the underlying energetic profiles to themselves be complex-valued. We illustrate this approach for 2D-IR spectra of water and SFG spectra of the air-water interface and demonstrate the accuracy of the generalized van't Hoff relationship and its implications for the origin of temperature-dependent spectral changes.

Nuclear quantum effects on the vibrational dynamics of the water-air interface

We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.

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.

An electrochemical perspective on the interfacial width between two immiscible liquid phases

Molecular dynamics simulations and vibrational sum-frequency spectroscopy are historically the main techniques applied to the description of the molecular structure and dynamics of the immiscible liquid/liquid interface. A molecular sharpness is estimated for oil/water interfaces, with an interfacial width that extends from hundreds of Å to 1 nm. However, electrochemical studies have elucidated a deeper liquid/liquid interface on the order of several micrometers. The breaking down of single-entity electrochemistry to simpler systems and the combination of high-resolution microscopies is confirming a larger extension of the interface. What can be the role of the electrochemist in clarifying this fundamental question? We try to give a suggestion at the end of a brief historical overview of the liquid/liquid interface studies.

Direct observation of bicarbonate and water reduction on gold: understanding the potential dependent proton source during hydrogen evolution

The electrochemical conversion of CO2 represents a promising way to simultaneously reduce CO2 emissions and store chemical energy. However, the competition between CO2 reduction (CO2R) and the H2 evolution reaction (HER) hinders the efficient conversion of CO2 in aqueous solution. In water, CO2 is in dynamic equilibrium with H2CO3, HCO3 -, and CO3 2-. While CO2 and its associated carbonate species represent carbon sources for CO2R, recent studies by Koper and co-workers indicate that H2CO3 and HCO3 - also act as proton sources during HER (J. Am. Chem. Soc. 2020, 142, 4154-4161, ACS Catal. 2021, 11, 4936-4945, J. Catal. 2022, 405, 346-354), which can favorably compete with water at certain potentials. However, accurately distinguishing between competing reaction mechanisms as a function of potential requires direct observation of the non-equilibrium product distribution present at the electrode/electrolyte interface. In this study, we employ vibrational sum frequency generation (VSFG) spectroscopy to directly probe the interfacial species produced during competing HER/CO2R on Au electrodes. The vibrational spectra at the Ar-purged Na2SO4 solution/Au interface, where only HER occurs, show a strong peak around 3650 cm-1, which appears at the HER onset potential and is assigned to OH-. Notably, this species is absent for the CO2-purged Na2SO4 solution/gold interface; instead, a peak around 3400 cm-1 appears at catalytic potential, which is assigned to CO3 2- in the electrochemical double layer. These spectral reporters allow us to differentiate between HER mechanisms based on water reduction (OH- product) and HCO3 - reduction (CO3 2- product). Monitoring the relative intensities of these features as a function of potential in NaHCO3 electrolyte reveals that the proton donor switches from HCO3 - at low overpotential to H2O at higher overpotential. This work represents the first direct detection of OH- on a metal electrode produced during HER and provides important insights into the surface reactions that mediate selectivity between HER and CO2R in aqueous solution.

Solvents and Stabilization in Ionic Liquid Films

We report the interfacial structures and chemical environments of ionic liquid films as a function of dilution with molecular solvents and over a range of film thicknesses (a few micrometers). Data from spectroscopic ellipsometry and infrared spectroscopy measurements show differences between films comprised of neat ionic liquids, as well as films comprised of ionic liquids diluted with two molecular solvents (water and acetonitrile). While the water-diluted IL films follow thickness trends predicted by the Landau-Levich model, neat IL and IL/MeCN films deviate significantly from predicted behaviors. Specifically, these film thicknesses are far greater than the predicted values, suggesting enhanced intermolecular interactions or other non-Newtonian behaviors not captured by the theory. We correlate film thicknesses with trends in the infrared intensity profiles across film thicknesses and IL-solvent dilution conditions and interpret the changes from expected behaviors as varying amounts of the film volume existing in isotropic (bulk) vs anisotropic (interfacial) states. The hydrogen bonding network of water-diluted ionic liquids is implicated in the agreement of this system with the Landau-Levich model's thickness predictions.

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.

Ultrafast vibrational dynamics of the free OD at the air/water interface: Negligible isotopic dilution effect but large isotope substitution effect

Vibrational relaxation dynamics of the OH stretch of water at the air/water interface has been a subject of intensive research, facilitated by recent developments in ultrafast interface-selective nonlinear spectroscopy. However, a reliable determination of the vibrational relaxation dynamics in the OD stretch region at the air/D₂O interface has not been yet achieved. Here, we report a study of the vibrational relaxation of the free OD carried out by time-resolved heterodyne-detected vibrational sum frequency generation spectroscopy. The results obtained with the aid of singular value decomposition analysis indicate that the vibrational relaxation (T₁) time of the free OD at the air/D₂O interface and air/isotopically diluted water (HOD-H₂O) interfaces show no detectable isotopic dilution effect within the experimental error, as in the case of the free OH in the OH stretch region. Thus, it is concluded that the relaxation of the excited free OH/OD predominantly proceeds with their reorientation, negating a major contribution of the intramolecular energy transfer. It is also shown that the T₁ time of the free OD is substantially longer than that of the free OH, further supporting the reorientation relaxation mechanism. The large difference in the T₁ time between the free OD and the free OH (factor of ∼2) may indicate the nuclear quantum effect on the diffusive reorientation of the free OD/OH because this difference is significantly larger than the value expected for a classical rotational motion.

Polarization-Dependent Sum-Frequency Generation Spectroscopy for Ångstrom-Scale Depth Profiling of Molecules at Interfaces

The three-dimensional spatial distribution of molecules at soft matter interfaces is crucial for processes ranging from membrane biophysics to atmospheric chemistry. While several techniques can access surface composition, obtaining information on the depth distribution is challenging. We develop a noninvasive, polarization-resolved, surface-specific sum-frequency generation spectroscopy providing quantitative depth information. We demonstrate the technique on formic acid molecules at the air-water interface. With increasing molar fraction from 2.5% to 10%, the formic acid molecules shift, on average, ∼0.9  Å into the bulk. The consistency with the simulation data manifests that the technique allows for probing the Ångstrom-scale depth profile.

The Solvation-Induced Onsager Reaction Field Rather than the Double-Layer Field Controls CO2 Reduction on Gold

The selectivity and activity of the carbon dioxide reduction (CO₂R) reaction are sensitive functions of the electrolyte cation. By measuring the vibrational Stark shift of in situ-generated CO on Au in the presence of alkali cations, we quantify the total electric field present at catalytic active sites and deconvolute this field into contributions from (1) the electrochemical Stern layer and (2) the Onsager (or solvation-induced) reaction field. Contrary to recent theoretical reports, the CO₂R kinetics does not depend on the Stern field but instead is closely correlated with the strength of the Onsager reaction field. These results show that in the presence of adsorbed (bent) CO₂, the Onsager field greatly exceeds the Stern field and is primarily responsible for CO₂ activation. Additional measurements of the cation-dependent water spectra using vibrational sum frequency generation spectroscopy show that interfacial solvation strongly influences the CO₂R activity. These combined results confirm that the cation-dependent interfacial water structure and its associated electric field must be explicitly considered for accurate understanding of CO₂R reaction kinetics.

Protonation State of Dopamine Neurotransmitter at the Aqueous Interface: Vibrational Sum Frequency Generation Spectroscopy Study

Dopamine is an important amine-based chemical neurotransmitter whose protonated state plays a crucial role in the recognition process. Understanding the structure and protonated state of dopamine at the aqueous interface is desired as the diffusion as well as binding of dopamine with the receptors take place frequently in the aqueous interface region. Vibrational sum frequency generation (VSFG) study of the OH stretch of water at the air/water interface in the presence of dopamine is performed and compared with its analog, phenylethylamine, and catechol. The VSFG data suggest that, unlike the bulk case, the population of the deprotonated amine group of dopamine is higher at the aqueous interface. This study suggests that the structure of dopamine at the aqueous interface is different from the bulk which may be useful in understanding the recognition process of dopamine in the interfacial region.

Bulk phase behaviour vs interface adsorption: Effects of anions and isotopes on β-lactoglobulin (BLG) interactions

HYPOTHESIS

Protein adsorption is highly relevant in numerous applications ranging from food processing to medical implants. In this context, it is important to gain a deeper understanding of protein-protein and protein-surface interactions. Thus, the focus of this investigation is on the interplay of bulk properties and surface properties on protein adsorption. It was hypothesised that the type of solvent and ions in solution should significantly influence the protein's bulk and interface behaviour, which has been observed in literature and previous work for other net negatively charged, globular proteins such as bovine serum albumin (BSA).

EXPERIMENTS

The phase behaviour of β-lactoglobulin (BLG) with lanthanum chloride (LaCl₃) and iodide (LaI₃) in normal water H₂O(l) and heavy water (D₂O(l)) was established via optical microscopy and ultraviolet-visible spectroscopy. The formation of an adsorption layer and its properties such as thickness, density, structure, and hydration was investigated via neutron reflectivity, quartz-crystal microbalance with dissipation, and infra-red measurements.

FINDINGS

β-lactoglobulin does not show significant anion-induced or isotope-induced effects - neither in bulk nor at the solid-liquid interface, which deviates strongly from the behaviour of bovine serum albumin. We also provide a comprehensive discussion and comparison of protein-specific bulk and interface behaviour between bovine serum albumin and β-lactoglobulin dependent on anion, cation, solvent, and substrate properties. These findings pave the way for understanding the transition from adsorption to crystallisation.

The molecular structure of the surface of water-ethanol mixtures

Mixtures of water and alcohol exhibit an excess surface concentration of alcohol as a result of the amphiphilic nature of the alcohol molecule, which has important consequences for the physico-chemical properties of water-alcohol mixtures. Here we use a combination of intensity vibrational sum-frequency generation (VSFG) spectroscopy, heterodyne-detected VSFG (HD-VSFG), and core-level photoelectron spectroscopy (PES) to investigate the molecular properties of water-ethanol mixtures at the air-liquid interface. We find that increasing the ethanol concentration up to a molar fraction (MF) of 0.1 leads to a steep increase of the surface density of the ethanol molecules, and an increased ordering of the ethanol molecules at the surface. When the ethanol concentration is further increased, the surface density of ethanol remains more or less constant, while the orientation of the ethanol molecules becomes increasingly disordered. The used techniques of PES and VSFG provide complementary information on the density and orientation of ethanol molecules at the surface of water, thus providing new information on the molecular-scale properties of the surface of water-alcohol mixtures over a wide range of compositions. This information is invaluable in understanding the chemical and physical properties of water-alcohol mixtures.

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.

Associated molecular liquids at the graphene monolayer interface

We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.

Protein Denaturation, Zero Entropy Temperature, and the Structure of Water around Hydrophobic and Amphiphilic Solutes

The hydrophobic effect plays a key role in many chemical and biological processes, including protein folding. Nonetheless, a comprehensive picture of the effect of temperature on hydrophobic hydration and protein denaturation remains elusive. Here, we study the effect of temperature on the hydration of model hydrophobic and amphiphilic solutes, through molecular dynamics, aiming at getting insight on the singular behavior of water, concerning the zero-entropy temperature, T_S, and entropy convergence, T_S^*, also observed for some proteins, upon denaturation. We show that, similar to hydrocarbons, polar amphiphilic solutes exhibit a T_S, although strongly dependent on solute-water interactions, opposite to hydrocarbons. Further, the temperature dependence of the hydration entropy, normalized by the solvent accessible surface area, is shown to be nearly solute size independent for hydrophobic but not for amphiphilic solutes, for similar reasons. These results are further discussed in the light of information theory (IT) and the structure of water around hydrophobic groups. The latter shows that the tetrahedral enhancement of some water molecules around hydrophobic groups, associated with the reduction of water defects, leads to the strengthening of the weakest hydrogen bonds, relative to bulk water. In addition, a larger tetrahedrality is found in low density water populations, demonstrating that pure water has encoded structural information, similar to that associated with hydrophobic hydration. The reversal of the hydration entropy dependence on the solute size, above T_S^*, is also analyzed and shown to be associated with a greater loss of water molecules exhibiting enhanced orientational order, in the coordination sphere of large solutes. Finally, the source of the differences between Kauzmann's "hydrocarbon model" on protein denaturation and hydrophobic hydration is discussed, with relatively large amphiphilic hydrocarbons seemingly displaying a more similar behavior to some globular proteins than aliphatic hydrocarbons.

"On-The-Fly" Calculation of the Vibrational Sum-Frequency Generation Spectrum at the Air-Water Interface

In the present work, we provide an electronic structure based method for the "on-the-fly" determination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polarizability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

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.

Strong Hydration at the Poly(ethylene glycol) Brush/Albumin Solution Interface

Albumin molecules are extensively used as biocompatible coatings, and poly(ethylene glycol) (PEG) materials are widely used for antifouling. PEG materials have excellent antifouling property because of their strong surface hydration. Our previous research indicates that hydration at the PEG/bovine serum albumin solution interface is stronger than that at the PEG/water interface. This research shows that this observation is general for different types of albumin molecules. Different albumins including bovine, porcine, rat, rabbit, and sheep serum albumins were studied in this research. It was found that the hydration at the PEG methacrylate (pOEGMA)/albumin solution interface is always stronger than that at the pOEGMA/water interface. Here, we define "strong interfacial hydration" as "ordered strongly hydrogen-bonded interfacial water". We believe that such a strong hydration is because of the strong hydration on the albumin surface, leading to its biocompatible property. All of the albumin molecules demonstrated stronger hydration on the pOEGMA surface compared to other protein molecules such as lysozyme and fibrinogen. The strong hydration on albumin molecules is related to the high surface coverage of glutamic acid and lysine with similar amounts.

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.

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.

Comparative Adsorption of Acetone on Water and Ice Surfaces

Small organic molecules on ice and water surfaces are ubiquitous in nature and play a crucial role in many environmentally relevant processes. Herein, we combine surface‐specific vibrational spectroscopy and a controllable flow cell apparatus to investigate the molecular adsorption of acetone onto the basal plane of single‐crystalline hexagonal ice with a large surface area. By comparing the adsorption of acetone on the ice/air and the water/air interface, we observed two different types of acetone adsorption, as apparent from the different responses of both the free O−H and the hydrogen‐bonded network vibrations for ice and liquid water. Adsorption on ice occurs preferentially through interactions with the free OH group, while the interaction of acetone with the surface of liquid water appears less specific.

Ab initio Modeling of the Vibrational Sum-Frequency Generation Spectrum of Interfacial Water

Understanding the structural and dynamical features of interfacial water is of greatest interest in physics, chemistry, biology, and materials science. Vibrational sum-frequency generation (SFG) spectroscopy, which is sensitive to the molecular orientation and dynamics on the surfaces or at the interfaces, allows one to study a wide variety of interfacial systems. The structural and dynamical features of interfacial water at the air/water interface have been extensively investigated by SFG spectroscopy. However, the interpretations of the spectroscopic features have been under intense debate. Here, we report a simulated SFG spectrum of the air/water interface based on ab initio molecular dynamics simulations, which covers the OH stretching, bending, and libration modes of interfacial water. Quantitative agreement between our present simulations and the most recent experimental studies ensures that ab initio simulations predict unbiased structural features and electrical properties of interfacial systems. By utilizing the kinetic energy spectral density (KESD) analysis to decompose the simulated spectra, the spectroscopic features can then be assigned to specific hydrogen-bonding configurations of interfacial water molecules.

Calcium Ions Affect Water Molecular Structures Surrounding an Oligonucleotide Duplex as Revealed by Sum Frequency Generation Vibrational Spectroscopy

The solvation of DNA in water facilitates the formation of a hydration layer surrounding it, thus stabilizing the DNA duplex in the biological aqueous environment. In this study, via using the lipid bilayer as a soft substrate to accommodate the duplex oligonucleotide, the structure of the water layer surrounding the oligonucleotide was detected under the perturbation of the calcium ions (Ca²⁺) with chiral and achiral sum frequency generation (SFG) vibrational spectroscopy. With increasing Ca²⁺ concentration, both the chiral and achiral water vibrational signals had similar concentration-dependent changes, i.e., an initial decreasing phase followed by an increasing phase. However, when the Ca²⁺ concentrations were adjusted to within the range comparable to those in the human serum, the chiral water vibrational signals remained nearly unchanged, whereas the achiral water vibrational signals still changed as a function of the Ca²⁺ concentration. Therefore, the current experimental result supports the possible protection function of the chiral hydration layer against the Ca²⁺ ions, which generally exist in the cell sap and play important roles in many biological functions.

Neat Water-Vapor Interface: Proton Continuum and the Nonresonant Background

Whether the surface of neat water is "acidic" or "basic" remains an active and controversial field of research. Most of the experimental evidence supporting the preferential adsorption of H₃O⁺ ions stems from nonlinear optical spectroscopy methods typically carried out at extreme pH conditions (pH < 1). Here, we use vibrational sum frequency spectroscopy (VSFS) to target the "proton continuum", an unexplored frequency range characteristic of hydrated protons and hydroxide ions. The VSFS spectra of neat water show a broad and nonzero signal intensity between 1700 and 3000 cm^-1 in the three different polarization combinations examined. By comparing the SF response of water with that from dilute HCl and NaOH aqueous solutions, we conclude the intensity does not originate from either adsorbed H₃O⁺ or OH^- ions. Contributions from the nonresonant background are then critically considered by comparing the experimental results with many-body molecular dynamics (MB-MD) simulated spectra.

Combined Sum Frequency Generation and Thin Liquid Film Study of the Specific Effect of Monovalent Cations on the Interfacial Water Structure

Some salts have been recently shown to decrease the sum frequency generation (SFG) intensity of the hydrogen-bonded water molecules, but a quantitative explanation is still awaited. Here, we report a similar trend for the chloride salts of monovalent cations, that is, LiCl, NaCl, and CsCl, at low concentrations. Specifically, we revealed not only the specific adsorption of cations at the water surface but also the concentration-dependent effect of ions on the SFG response of the interfacial water molecules. Our thin-film pressure balance (TFPB) measurements (stabilized by 10 mM of methyl isobutyl carbinol) enabled the determination of the surface potential that governs the surface electric field affecting interfacial water dipoles. The use of the special alcohol also enabled us to identify a remarkable specific screening effect of cations on the surface potential. We explained the concentration dependency by considering the direct ion-water interactions and water reorientation under the influence of surface electric field as the two main contributors to the overall SFG signal of the hydrogen-bonded water molecules. Although the former was dominant only at the low-concentration range, the effect of the latter intensified with increasing salt concentration, leading to the recovery of the band intensity at medium concentrations. We discussed the likelihood of a correlation between the effect of ions on reorientation dynamics of water molecules and the broad-band intensity drop in the SFG spectra of salt solutions. We proposed a mechanism for the cation-specific effect through the formation of an ionic capacitance at the solution surface. It explains how cations could impart the ion specificity while they are traditionally believed to be repelled from the interfacial region. The electrical potential of this capacitance varies with the charge separation and ion density at the interface. The charge separation being controlled by the polarizability difference between anions and cations was identified using the SFG response of the interfacial water molecules as an indirect probe. The ion density being affected by the absolute polarizability of ions was tracked through the measurement of the surface potentials and Debye-Hückel lengths using the TFPB technique.

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.

Structure from Dynamics: Vibrational Dynamics of Interfacial Water as a Probe of Aqueous Heterogeneity

The structural heterogeneity of water at various interfaces can be revealed by time-resolved sum-frequency generation spectroscopy. The vibrational dynamics of the O–H stretch vibration of interfacial water can reflect structural variations. Specifically, the vibrational lifetime is typically found to increase with increasing frequency of the O–H stretch vibration, which can report on the hydrogen-bonding heterogeneity of water. We compare and contrast vibrational dynamics of water in contact with various surfaces, including vapor, biomolecules, and solid interfaces. The results reveal that variations in the vibrational lifetime with vibrational frequency are very typical, and can frequently be accounted for by the bulk-like heterogeneous response of interfacial water. Specific interfaces exist, however, for which the behavior is less straightforward. These insights into the heterogeneity of interfacial water thus obtained contribute to a better understanding of complex phenomena taking place at aqueous interfaces, such as photocatalytic reactions and protein folding.

Oriented chiral water wires in artificial transmembrane channels

Aquaporins (AQPs) feature highly selective water transport through cell membranes, where the dipolar orientation of structured water wires spanning the AQP pore is of considerable importance for the selective translocation of water over ions. We recently discovered that water permeability through artificial water channels formed by stacked imidazole I-quartet superstructures increases when the channel water molecules are highly organized. Correlating water structure with molecular transport is essential for understanding the underlying mechanisms of (fast) water translocation and channel selectivity. Chirality adds another factor enabling unique dipolar oriented water structures. We show that water molecules exhibit a dipolar oriented wire structure within chiral I-quartet water channels both in the solid state and embedded in supported lipid bilayer membranes (SLBs). X-ray single-crystal structures show that crystallographic water wires exhibit dipolar orientation, which is unique for chiral I-quartets. The integration of I-quartets into SLBs was monitored with a quartz crystal microbalance with dissipation, quantizing the amount of channel water molecules. Nonlinear sum-frequency generation vibrational spectroscopy demonstrates the first experimental observation of dipolar oriented water structures within artificial water channels inserted in bilayer membranes. Confirmation of the ordered confined water is obtained via molecular simulations, which provide quantitative measures of hydrogen bond strength, connectivity, and the stability of their dipolar alignment in a membrane environment. Together, uncovering the interplay between the dipolar aligned water structure and water transport through the self-assembled I-quartets is critical to understanding the behavior of natural membrane channels and will accelerate the systematic discovery for developing artificial water channels for water desalting.

DNA's Chiral Spine of Hydration

The iconic helical structure of DNA is stabilized by the solvation environment, where a change in the hydration state can lead to dramatic changes to the DNA structure. X-ray diffraction experiments at cryogenic temperatures have shown crystallographic water molecules in the minor groove of DNA, which has led to the notion of a spine of hydration of DNA. Here, chiral nonlinear vibrational spectroscopy of two DNA sequences shows that not only do such structural water molecules exist in solution at ambient conditions but that they form a chiral superstructure: a chiral spine of hydration. This is the first observation of a chiral water superstructure templated by a biomolecule. While the biological relevance of a chiral spine of hydration is unknown, the method provides a direct way to interrogate the properties of the hydration environment of DNA and water in biological settings without the use of labels.

A general method for determining molecular interfaces and layers

A general and direct computational scheme to locate the surface separating arbitrarily shaped domains made up of molecules (or any other particles) has been developed and is described and illustrated for several, both artificial and physical examples. The proposed scheme consists of two modules: (i) triangulation and (ii) assignment of simplices to domains. Three different triangulation methods are employed, viz., the Delaunay triangulation, regular triangulation, and quasi-triangulation. In the triangulated system, the assignment step is carried out in two different ways, one based on the characteristic metric of a particular triangulation procedure and the other on the concept of a touching sphere. Some of the combinations of the triangulation and assignment steps lead to methods already used by others to find interfacial or surface molecules, namely the alpha-shape-based method of Usabiaga nad Duque [Phys. Rev. E 79 (2009) 046709] and GITIM of Sega et al. [J. Chem. Phys. 138 (2013) 044110]. The resulting surface is defined not only as a discrete set of particles, but it is build up of facets of the triangulation forming a broken line in two dimensions or a polyhedral surface in three dimensions. Individual molecular layers are identified in a very straightforward manner, starting with the interfacial layer itself and proceeding into the interior of the phase. The proposed scheme is illustrated first by identifying border molecules of pre-sampled domains of several shapes in a plane and then applied to five physically meaningful examples: thin films, near critical water, liquid water slab in an electric field, liquid water at a solid wall, and water at condition of electric-field-induced jetting. Performance of the considered methods is critically assessed. Treatment of domains forming percolating clusters through periodic boundary conditions is also described along with the determination of their periodicity and dimensionality.

On the Assignment of the Vibrational Spectrum of the Water Bend at the Air/Water Interface

We previously reported the spectrum of the water bend vibrational mode (ν₂) at the air/water interface measured using sum-frequency generation (SFG). Here, we present experimental evidence to aid the assignment of the ν₂ spectral features to H-bonded classes of interfacial water, which is in general agreement with two recent independently published theoretical studies. The dispersive line shape shows an apparent frequency shift between SSP and PPP polarization combinations (SFG-visible-infrared). This is naturally explained as an interference effect between the negative (1630 cm^-1) and positive (1662 cm^-1) peaks corresponding to "free-OH" and "H-bonded" species, respectively, which have different orientations and thus different amplitudes in SSP and PPP spectra. A surfactant monolayer of sodium dodecyl sulfate (SDS) was used to suppress the free OH species at the surface, and the corresponding SFG spectral changes indicate that these water molecules with one of the hydrogens pointing up into the air phase contribute to the negative peak at 1630 cm^-1.

Heavy snow: IR spectroscopy of isotope mixed crystalline water ice

Temperature and isotopic dependence of simulated and experimental spectra shed light on the vibrational modes of crystalline water ice.

Fast and slow dynamics and the local structure of liquid and supercooled water next to a hydrophobic amino acid

We study, through molecular dynamics simulations, the structure and orientational dynamics of water next to a blocked hydrophobic amino acid, valine (Val), above and below the freezing point of water.

A QM:MM model for the interaction of DNA nucleotides with carbon nanotubes

Hybrid materials formed by DNA and carbon nanotubes (CNTs) have shown very interesting properties, but their simulation in solution using quantum mechanical approaches is still a challenge in the computational chemistry community. In this paper, we developed a QM:MM model to study the interactions between charged DNA nucleotides and carbon nanotubes in solution. All four types of DNA nucleotides were taken to interact with two CNTs of similar diameter but different chiralities: (4,4) and (7,0). The nucleotides and CNTs were treated at the QM level, while added water and neutralizing ions were modeled at the MM level. ONIOM simulations were performed at the (M06-2X/6-31G(d):Amber) level for the hybrids, as well as for individually solvated CNTs and nucleotides, which allowed us to evaluate the energy of binding. Our binding energy (BE) values range from 146.60 to 503.43 kJ mol(-1), indicating strong physisorption of nucleotides on CNTs. The relatively large BE, compared with past studies on nucleobase-CNT binding in a vacuum, could be due to the larger size of nucleotides compared with nucleobases, the charges on the nucleotides, and the inclusion of solution which causes the release of water molecules upon hybridization.

Estimating pH at the Air/Water Interface with a Confocal Fluorescence Microscope

One way to determine the pH at the air/water interface with a confocal fluorescence microscope has been proposed. The relation between the pH at the air/water interface and that in a bulk solution has been formulated in connection with the adsorption equilibrium and the dissociation equilibrium of the dye adsorbed. Rhodamine B (RhB) is used as a surface-active fluorescent pH probe. The corrected fluorescence spectrum of RhB molecules at the air/water interface with the surface density of 1.0 nmol m(-2) level shows pH-dependent shifts representing an acid-base equilibrium. Two ways to determine the unknown acid-base equilibrium constant of RhB molecules at the air/water interface have been discussed. With surface-tension measurements, the adsorption properties, maximum surface density, and adsorption equilibrium constants were estimated for both cationic and zwitterionic forms of RhB molecules at the air/water interface.

A means to an interface: investigating monoethanolamine behavior at an aqueous surface

The use of amine scrubbers to trap carbon dioxide from flue gas streams is one of the most promising avenues for atmospheric carbon dioxide reduction. However, modifications are necessary to efficiently scale these scrubbers for use in fossil fuel plants. Current advances in tailoring amines for CO2 capture involve improvements of bulk kinetic and thermodynamic parameters, with little consideration to surface chemistry and behavior. Aqueous alkanolamine solutions, such as monoethanolamine (MEA), are currently highly favored sorbents in CO2 post-combustion capture. Although numerous studies have explored MEA-CO2 chemistry at the macroscopic scale, few have investigated the role of the interface in the gas adsorption process. Additionally, as these amines become more industrially ubiquitous, their presence on and the need to understand their behavior at atmospheric and environmental surfaces will increase. This study investigates the surface behavior of monoethanolamine at the vapor/water interface, with particular focus on MEA's surface orientation and footprint. Using vibrational sum frequency spectroscopy, surface tensiometry, and computational techniques, MEA is found to adopt a constrained gauche interfacial conformation with its methylene backbone oriented toward the vapor phase and its functional groups solvated in the bulk solution. Computational and experimental analysis agree well, giving a complete picture with vibrational mode assignments and surface orientation of MEA. These findings can assist in the tailoring of amine structures or to facilitate improvements in engineering design to exploit favorable surface chemistry, as well as to serve as a starting point toward understanding aqueous amine surface behavior relevant to environmental systems.

Bilayer Charge Reversal and Modification of Lipid Organization by Dendrimers as Observed by Sum-Frequency Vibrational Spectroscopy

Polyamidoamine (PAMAM) dendrimers are hyperbranched, nanosized polymers with promising biomedical applications as nanocarriers in targeted drug delivery and gene therapy. For the development of safe dendrimer-based biomedical applications it is necessary to gain an understanding of the detailed mechanism of the interactions of both cationic and anionic dendrimers with cell membranes. To characterize dendrimer-membrane interactions we applied solid-supported lipid bilayers as biomembrane models and utilized infrared-visible sum-frequency vibrational spectroscopy to independently probe the interactions of cationic G5-NH2 and anionic G4.5-COONa dendrimers with the two leaflets of the lipid bilayers. Interaction with both dendrimers led to changes in the interfacial water structure and charge density as evidenced by the changes in the OH band intensities in the sum-frequency spectra of the bilayers. Interaction with the G5-NH2 dendrimer also led to a unique inversion of the sign of the OH-stretch amplitudes, in addition to a decrease in their absolute values. We suggest that the positively charged amino groups on the G5-NH2 dendrimer surface bind to the negatively charged bilayer, while uncompensated positive charges not involved in the binding cause a reversal of the electric field and thus an opposite orientation of the interfacial water molecules. More subtle but nonetheless significant changes were seen in the relative magnitudes of the CH amplitudes. The methyl antisymmetric to symmetric stretch amplitude ratios are altered, implying changes in the tilt angles of the phospholipid alkyl chains. The conformational order of the phospholipid alkyl chains of both leaflets is also influenced by the G5-NH2 dendrimer while G4.5-COONa has no effect on the alkyl chain conformation.