Interpretation of the water surface vibrational sum-frequency spectrum

Beyond the "spine of hydration": Chiral SFG spectroscopy detects DNA first hydration shell and base pair structures

Experimental methods capable of selectively probing water at the DNA minor groove, major groove, and phosphate backbone are crucial for understanding how hydration influences DNA structure and function. Chiral-selective sum frequency generation spectroscopy (chiral SFG) is unique among vibrational spectroscopies because it can selectively probe water molecules that form chiral hydration structures around biomolecules. However, interpreting chiral SFG spectra is challenging since both water and the biomolecule can produce chiral SFG signals. Here, we combine experiment and computation to establish a theoretical framework for the rigorous interpretation of chiral SFG spectra of DNA. We demonstrate that chiral SFG detects the N-H stretch of DNA base pairs and the O-H stretch of water, exclusively probing water molecules in the DNA first hydration shell. Our analysis reveals that DNA transfers chirality to water molecules only within the first hydration shell, so they can be probed by chiral SFG spectroscopy. Beyond the first hydration shell, the electric field-induced water structure is symmetric and, therefore, precludes chiral SFG response. Furthermore, we find that chiral SFG can differentiate chiral subpopulations of first hydration shell water molecules at the minor groove, major groove, and phosphate backbone. Our findings challenge the scientific perspective dominant for more than 40 years that the minor groove "spine of hydration" is the only chiral water structure surrounding the DNA double helix. By identifying the molecular origins of the DNA chiral SFG spectrum, we lay a robust experimental and theoretical foundation for applying chiral SFG to explore the chemical and biological physics of DNA hydration.

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

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

Effect of counterions on the structure and dynamics of water near a negatively charged surfactant: a theoretical vibrational sum frequency generation study

Charged aqueous interfaces are of paramount importance in electrochemical, biological and environmental sciences. The properties of aqueous interfaces with ionic surfactants can be influenced by the presence of counterions. Earlier experiments involving vibrational sum frequency generation (VSFG) spectroscopy of aqueous interfaces with negatively charged sodium dodecyl sulfate (Na+DS- or SDS) surfactants revealed that the hydrogen bonding strength of the interfacial water molecules follows a certain order when salts of monovalent and divalent cations are added. It is known that cations do not directly participate in hydrogen bonding with water molecules, rather they only influence the hydrogen bonded network through their electrostatic fields. In the current work, we have simulated the aqueous interfacial systems of sodium dodecyl sulfate in the presence of chloride salts of mono and divalent countercations. The electronic polarization effects on the ions are considered at a mean-field level within the electronic continuum correction model. Our calculations of the VSFG spectra show a blue shift in the presence of added countercations whose origin is traced to different relative contributions of water molecules from the solvation shells of the surfactant headgroups and the remaining water molecules in the presence of countercations. Furthermore, the cations shield the electric fields of the surfactant headgroups, which in turn influences the contributions of water molecules to the total VSFG spectrum. This shielding effect becomes more significant when divalent countercations are present. The dynamics of water molecules is found to be slower at the interface in comparison to the bulk. The interfacial depth dependence of various dynamical quantities shows that the interface is structurally and dynamically more heterogeneous at the microscopic level.

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.

Neural Network-Based Sum-Frequency Generation Spectra of Pure and Acidified Water Interfaces with Air

The affinity of hydronium ions (H3O+) for the air-water interface is a crucial question in environmental chemistry. While sum-frequency generation (SFG) spectroscopy has been instrumental in indicating the preference of H3O+ for the interface, key questions persist regarding the molecular origin of the SFG spectral changes in acidified water. Here we combine nanosecond long neural network (NN) reactive simulations of pure and acidified water slabs with NN predictions of molecular dipoles and polarizabilities to calculate SFG spectra of long reactive trajectories including proton transfer events. Our simulations show that H3O+ ions cause two distinct changes in phase-resolved SFG spectra: first, a low-frequency tail due to the vibrations of H3O+ and its first hydration shell, analogous to the bulk proton continuum, and second, an enhanced hydrogen-bonded band due to the ion-induced static field polarizing molecules in deeper layers. Our calculations confirm that changes in the SFG spectra of acidic solutions are caused by hydronium ions preferentially residing at the interface.

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

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

Theoretical basis for interpreting heterodyne chirality-selective sum frequency generation spectra of water

Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O-H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O-H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O-H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O-H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model β-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra.

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.

Detecting Interplay of Chirality, Water, and Interfaces for Elucidating Biological Functions

Chemists have long been fascinated by chirality, water, and interfaces, making tremendous progress in each research area. However, the chemistry emerging from the interplay of chirality, water, and interfaces has been difficult to study due to technical challenges, creating a barrier to elucidating biological functions at interfaces. Most biopolymers (proteins, DNA, and RNA) fold into macroscopic chiral structures to perform biological functions. Their folding requires water, but water behaves differently at interfaces where the bulk water hydrogen-bonding network terminates. A question arises as to how water molecules rearrange to minimize free energy at interfaces while stabilizing the macroscopic folding of biopolymers to support biological function. This question is central to solving many research challenges, including the molecular origin of biological homochirality, folding and insertion of proteins into cell membranes, and the design of heterogeneous biocatalysts. Researchers can resolve these challenges if they have the theoretical tools to accurately predict molecular behaviors of water and biopolymers at various interfaces. However, developing such tools requires validation by the experimental data. These experimental data are scarce because few physical methods can simultaneously distinguish chiral folding of the biopolymers, separate signals of interfaces from the overwhelming background of bulk solvent, and differentiate water in hydration shells of the polymers from water elsewhere.We recently illustrated these very capacities of chirality-sensitive vibrational sum frequency generation spectroscopy (chiral SFG). While chiral SFG theory dictates that the method is surface-specific under the condition of electronic nonresonance, we show the method can distinguish chiral folding of proteins and DNA and probe water structures in the first hydration shell of proteins at interfaces. Using amide I signals, we observe protein folding into β-sheets without background signals from α-helices and disordered structures at interfaces, thereby demonstrating the effect of 2D crowding on protein folding. Also, chiral SFG signals of C-H stretches are silent from single-stranded DNA, but prominent for canonical antiparallel duplexes as well as noncanonical parallel duplexes at interfaces, allowing for sensing DNA secondary structures and hybridization. In establishing chiral SFG for detecting protein hydration structures, we observe an H218O isotopic shift that reveals water contribution to the chiral SFG spectra. Additionally, the phase of the O-H stretching bands flips when the protein chirality is switched from L to D. These experimental results agree with our simulated chiral SFG spectra of water hydrating the β-sheet protein at the vacuum-water interface. The simulations further reveal that over 90% of the total chiral SFG signal comes from water in the first hydration shell. We conclude that the chiral SFG signals originate from achiral water molecules that assemble around the protein into a chiral supramolecular structure with chirality transferred from the protein. As water O-H stretches can reveal hydrogen-bonding interactions, chiral SFG shows promise in probing the structures and dynamics of water-biopolymer interactions at interfaces. Altogether, our work has created an experimental and computational framework for chiral SFG to elucidate biological functions at interfaces, setting the stage for probing the intricate chemical interplay of chirality, water, and interfaces.

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

Design of an Electrostatic Frequency Map for the NH Stretch of the Protein Backbone and Application to Chiral Sum Frequency Generation Spectroscopy

We develop an electrostatic map for the vibrational NH stretch (amide A) of the protein backbone with a focus on vibrational chiral sum frequency generation spectroscopy (chiral SFG). Chiral SFG has been used to characterize protein secondary structure at interfaces using the NH stretch and to investigate chiral water superstructures around proteins using the OH stretch. Interpretation of spectra has been complicated because the NH stretch and OH stretch overlap spectrally. Although an electrostatic map for water OH developed by Skinner and co-workers was used previously to calculate the chiral SFG response of water structures around proteins, a map for protein NH that is directly responsive to biological complexity has yet to be developed. Here, we develop such a map, linking the local electric field to vibrational frequencies and transition dipoles. We apply the map to two protein systems and achieve much better agreement with experiment than was possible in our previous studies. We show that couplings between NH and OH vibrations are crucial to the line shape, which informs the interpretation of chiral SFG spectra, and that the chiral NH stretch response is sensitive to small differences in structure. This work increases the utility of the NH stretch in biomolecular spectroscopy.

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.

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

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

Detecting the First Hydration Shell Structure around Biomolecules at Interfaces

Understanding the role of water in biological processes remains a central challenge in the life sciences. Water structures in hydration shells of biomolecules are difficult to study in situ due to overwhelming background from aqueous environments. Biological interfaces introduce additional complexity because biomolecular hydration differs at interfaces compared to bulk solution. Here, we perform experimental and computational studies of chiral sum frequency generation (chiral SFG) spectroscopy to probe chirality transfer from a protein to the surrounding water molecules. This work reveals that chiral SFG probes the first hydration shell around the protein almost exclusively. We explain the selectivity to the first hydration shell in terms of the asymmetry induced by the protein structure and specific protein-water hydrogen-bonding interactions. This work establishes chiral SFG as a powerful technique for studying hydration shell structures around biomolecules at interfaces, presenting new possibilities to address grand research challenges in biology, including the molecular origins of life.

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.

Review of High-Frequency Ultrasounds Emulsification Methods and Oil/Water Interfacial Organization in Absence of any Kind of Stabilizer

Emulsions are multiphasic systems composed of at least two immiscible phases. Emulsion formulation can be made by numerous processes such as low-frequency ultrasounds, high-pressure homogenization, microfluidization, as well as membrane emulsification. These processes often need emulsifiers' presence to help formulate emulsions and to stabilize them over time. However, certain emulsifiers, especially chemical stabilizers, are less and less desired in products because of their negative environment and health impacts. Thus, to avoid them, promising processes using high-frequency ultrasounds were developed to formulate and stabilize emulsifier-free emulsions. High-frequency ultrasounds are ultrasounds having frequency greater than 100 kHz. Until now, emulsifier-free emulsions' stability is not fully understood. Some authors suppose that stability is obtained through hydroxide ions' organization at the hydrophobic/water interfaces, which have been mainly demonstrated by macroscopic studies. Whereas other authors, using microscopic studies, or simulation studies, suppose that the hydrophobic/water interfaces would be rather stabilized thanks to hydronium ions. These theories are discussed in this review.

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.

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.

Chiral Sum Frequency Generation Spectroscopy Detects Double-Helix DNA at Interfaces

Many DNA-based technologies involve the immobilization of DNA and therefore require a fundamental understanding of the DNA structure-function relationship at interfaces. We present three immobilization methods compatible with chiral sum frequency generation (SFG) spectroscopy at interfaces. They are the "anchor" method for covalently attaching DNA on a glass surface, the "island" method for dropcasting DNA on solid substrates, and the "buoy" method using a hydrocarbon moiety for localizing DNA at the air-water interface. Although SFG was previously used to probe DNA, the chiral and achiral SFG responses of single-stranded and double-stranded DNA have not been compared systemically. Using the three immobilization methods, we obtain the achiral and chiral C-H stretching spectra. The results introduce four potential applications of chiral SFG. First, chiral SFG gives null response from single-stranded DNA but prominent signals from double-stranded DNA, providing a simple binary readout for label-free detection of DNA hybridization. Second, with heterodyne detection, chiral SFG gives an opposite-signed spectral response useful for distinguishing native (D-) right-handed double helix from non-native (L-) left-handed double helix. Third, chiral SFG captures the aromatic C-H stretching modes of nucleobases that emerge upon hybridization, revealing the power of chiral SFG to probe highly localized molecular structures within DNA. Finally, chiral SFG is sensitive to macroscopic chirality but not local chiral centers and thus can detect not only canonical antiparallel double helix but also other DNA secondary structures, such as a poly-adenine parallel double helix. Our work benchmarks the SFG responses of DNA immobilized by the three distinct methods, building a basis for new chiral SFG applications to solve fundamental and biotechnological problems.

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

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

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

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

Bulk Response of Carboxylic Acid Solutions Observed with Surface Sum-Frequency Generation Spectroscopy

We study the molecular properties of aqueous acetic acid and formic acid solutions with heterodyne-detected vibrational sum-frequency generation spectroscopy (HD-VSFG). For acid concentrations up to ∼5 M, we observe a strong increase of the responses of the acid hydroxyl and carbonyl stretch vibrations with increasing acid concentration due to an increase of the surface coverage by the acid molecules. At acid concentrations >5 M we observe first a saturation of these responses and then a decrease. For pure carboxylic acids we even observe a change of sign of the Im[χ⁽²⁾] response of the carbonyl vibration. The decrease of the response of the hydroxyl vibration and the decrease and sign change of the response of the carbonyl vibration indicate the formation of cyclic dimers, which only show a quadrupolar bulk response in the HD-VSFG spectrum because of their antiparallel conformation. We also find evidence for the presence of a quadrupolar response of the CH vibrations of the acid molecules.

Simulation of the Chiral Sum Frequency Generation Response of Supramolecular Structures Requires Vibrational Couplings

Chiral vibrational sum frequency generation (SFG) spectroscopy probes the structure of the solvation shell around chiral macromolecules. The dominant theoretical framework for understanding the origin of chiral SFG signals is based on the analysis of molecular symmetry, which assumes no interaction between molecules. However, water contains strong intermolecular interactions that significantly affect its properties. Here, the role of intermolecular vibrational coupling in the chiral SFG response of the O-H stretch of water surrounding an antiparallel β-sheet at the vacuum-water interface is investigated. Both intramolecular and intermolecular couplings between O-H groups are required to simulate the full lineshape of the chiral SFG signal. This dependence is also observed for a chiral water dimer, illustrating that this phenomenon is not specific to larger systems. We also find that a dimer of C3v molecules predicted to be chirally SFG-inactive by the symmetry-based theory can generate a chiral SFG signal when intermolecular couplings are considered, suggesting that even highly symmetric solvent molecules may produce chiral SFG signals when interacting with a chiral solute. The consideration of intermolecular couplings extends the prevailing theory of the chiral SFG response to structures larger than individual molecules and provides guidelines for future modeling.

Structure and sum-frequency generation spectra of water on neutral hydroxylated silica surfaces

Structural organization and vibrational sum-frequency generation (VSFG) spectra of water on crystalline and amorphous neutral silica surfaces were investigated by classical molecular dynamics simulations. The liquid phase represented with neat water and 1 M NaCl solution was analysed in terms of bonded interfacial layer (BIL), diffuse layer (DL) and bulk region. The simulations show that the structure of BIL depends on the surface morphology and density of surface OH groups. The water-silanol H-bond network and BIL structure are mainly insensitive to the presence of ions in the liquid phase. Molecules in DL of SiO₂/neat water interfaces preferentially orient their OH bonds towards the surfaces. This effect is directly related to an effective negative charge of formally neutral surfaces. Ions of the electrolyte solution affect the intermolecular structure in DL by screening the surface electric field and by the chaotropic effect. Calculated phase-sensitive VSFG (Im[χ⁽²⁾]) spectrum of BIL features low-frequency negative and high-frequency positive bands. Characteristics of the positive band reflect the strength of water-surface interactions and surface crystallinity, while the position and shape of the negative band are common to all interfaces. The Im[χ⁽²⁾] spectrum of DL is dominated by a contribution from the third-order χ⁽³⁾ susceptibility with the sign of the contribution directly related to the sign of electrostatic potential in the interfacial region. The DL spectrum is strongly affected by the presence of solvated ions. The computed intensity and Im[χ⁽²⁾] spectra of the amorphous silica/NaCl solution interface are in a good agreement with the conventional and phase-sensitive experimental VSFG spectra of fused SiO₂/water system at low pH, in contrast to the spectra of the amorphous silica/neat water interface. Origins of the discrepancy are discussed.

"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.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

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

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

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

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

Water Structure in the Submembrane Region of a Floating Lipid Bilayer: The Effect of an Ion Channel Formation and the Channel Blocker

The structure of water in the submembrane region of the bilayer of DPhPC floating (fBLM) on a monolayer of 1-thio-β-d-glucose (β-Tg)-modified gold nanoparticle film was studied by the surface-enhanced infrared absorption spectroscopy (SEIRAS). SEIRAS employs surface enhancement of the mean square electric field of the photon, which is acting on a few molecular layers above the film of gold nanoparticles. Therefore, it is uniquely suited to probe water molecules in the submembrane region and provides unique information concerning the structure of the hydrogen bond network of water surrounding the lipid bilayer. The IR spectra indicated that water with a strong hydrogen network is separating the membrane from the gold surface. This water is more ordered than the water in the bulk. When alamethicin, a peptide forming ion channels, is inserted into the membrane, the network is only slightly loosened. The addition of amiloride, an ion channel blocker, results in a significant decrease in the amount of water in the submembrane region. The remaining water has a significantly distorted hydrogen bond network. This study provides unique information about the effect of the ion channel on water transport across the bilayer. The electrode potential has a relatively small effect on water structure in the submembrane region. However, the IR studies demonstrated that water is less ordered at positive transmembrane potentials. The present results provide significant insight into the nature of hydration of a floating lipid bilayer on the gold electrode surface.

Unveiling Heterogeneity of Interfacial Water through the Water Bending Mode

The water bending mode provides a powerful probe of the microscopic structure of bulk aqueous systems because its frequency and spectral line shape are responsive to the intermolecular interactions. Furthermore, interpreting the bending mode response is straightforward, as the intramolecular vibrational coupling is absent. Nevertheless, bending mode has not been used for probing the interfacial water structure, as it has been yet argued that the signal is dominated by bulk effects. Here, through the sum-frequency generation measurement of the water bending mode at the water/air and water/charged lipid interfaces, we demonstrate that the bending mode signal is dominated not by the bulk but by the interface. Subsequently, we disentangle the hydrogen-bonding of water at the water/air interface using the bending mode frequency distribution and find distinct interfacial hydrogen-bonded structures, which can be directly related to the interfacial organization of water. The bending mode thus provides an excellent probe of aqueous interfacial structure.

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.

Nuclear Quantum Effect on the χ(2) Band Shape of Vibrational Sum Frequency Generation Spectra of Normal and Deuterated Water Surfaces

In the vibrational sum frequency generation (VSFG) spectrum of the air/water interface, there is an open question whether the imaginary spectrum of nonlinear susceptibility Im[χ⁽²⁾] has a positive band in the low-frequency tail of the OH stretching band. We previously elucidated by molecular dynamics (MD) analysis of the VSFG spectrum that the positive band could arise from particularly strong hydrogen-bonded water pairs at the water surface. This mechanism should be emphasized in OD stretching in comparison to OH, because OD forms stronger hydrogen bonds. Therefore, we calculated the Im[χ⁽²⁾] spectra of normal and deuterated water surfaces by MD simulation including the nuclear quantum effect and demonstrated that the low-frequency positive feature could arise in the tail of the OD stretching band. This positive feature is sensitive to the oxygen-oxygen distance of hydrogen-bonding pairs at the surface and hence the temperature and disappears with increasing temperature.

Hydrogen bonding and vibrational energy relaxation of interfacial water: A full DFT molecular dynamics simulation

The air-water interface has been a subject of extensive theoretical and experimental studies due to its ubiquity in nature and its importance as a model system for aqueous hydrophobic interfaces. We report on the structure and vibrational energy transfer dynamics of this interfacial water system studied with equilibrium and non-equilibrium molecular dynamics simulations employing a density functional theory -based description of the system and the kinetic energy spectral density analysis. The interfacial water molecules are found to make fewer and weaker hydrogen (H)-bonds on average compared to those in the bulk. We also find that (i) the H-bonded OH groups conjugate to the free OH exhibit rather low vibrational frequencies (3000-3500 cm^-1); (ii) the presence of a significant fraction (>10%) of free and randomly oriented water molecules at the interface ("labile water"), neither of whose OH groups are strong H-bond donors; (iii) the inertial rotation of free OH groups, especially from the labile water, contribute to the population decay of excited free OH groups with comparable rate and magnitude as intramolecular energy transfer between the OH groups. These results suggest that the labile water, which might not be easily detectable by the conventional vibrational sum frequency generation method, plays an important role in the surface water dynamics.

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.

Interactions of Sulfobetaine Zwitterionic Surfactants with Water on Water Surface

We carried out a combined study using surface tension, phase-sensitive frequency generation (SFG) vibrational spectroscopy, and MD simulations to investigate the industrially relevant zwitterionic surfactant N-dodecyl-N, N-dimethyl-3-ammonio-1-propanesulfonate (DDAPS) on water surface. The SFG Im(χ⁽²⁾) spectra showed that the interaction between DDAPS and water was different from those between biologically relevant zwitterionic phospholipids and water. While zwitterionic phospholipids were found to be anionic-like and flipped water molecules with their OHs pointing toward the air, DDAPS oriented water molecules with their OHs mostly pointing toward the liquid water. We built a new force field for the MD simulation which produced the correct surface tension of water with various DDPAS coverage. The MD simulation showed that the head groups of DDPAS were nearly parallel to the water surface. When the surface coverage of DDPAS was increased, the averaged tilting angle of DDPAS's tails decreased, but it had little effect on the orientation of the headgroup. The sulfobetaine zwitterionic surfactant was found to be more cationic-like because the positively charged group was more capable of orienting interfacial water.

Water's dual nature and its continuously changing hydrogen bonds

A model is proposed for liquid water that is a continuum between the ordered state with predominantly tetrahedral coordination, linear hydrogen bonds and activated dynamics and a disordered state with a continuous distribution of multiple coordinations, multiple types of hydrogen bond, and diffusive dynamics, similar to that of normal liquids. Central to water's heterogeneous structure is the ability of hydrogen to donate to either one acceptor in a conventional linear hydrogen bond or to multiple acceptors as a furcated hydrogen. Linear hydrogen bonds are marked by slow, activated kinetics for hydrogen-bond switching to more crowded acceptors and sharp first peaks in the hydrogen-oxygen radial distribution function. Furcated hydrogens, equivalent to free, broken, dangling or distorted hydrogens, have barrierless, rapid kinetics and poorly defined first peaks in their hydrogen-oxygen radial distribution function. They involve the weakest donor in a local excess of donors, such that barrierless whole-molecule vibration rapidly swaps them between the linear and furcated forms. Despite the low number of furcated hydrogens and their transient existence, they are readily created in a single hydrogen-bond switch and free up the dynamics of numerous surrounding molecules, bringing about the disordered state. Hydrogens in the ordered state switch with activated dynamics to make the non-tetrahedral coordinations of the disordered state, which can also combine to make the ordered state. Consequently, the ordered and disordered states are both connected by diffusive dynamics and differentiated by activated dynamics, bringing about water's continuous heterogeneity.