2D heterodyne-detected sum frequency generation study on the ultrafast vibrational dynamics of H2O and HOD water at charged interfaces

Double-Layer Distribution of Hydronium and Hydroxide Ions in the Air-Water Interface

The acid-base nature of the aqueous interface has long been controversial. Most macroscopic experiments suggest that the air-water interface is basic based on the detection of negative charges at the interface that indicates the enrichment of hydroxides (OH-), whereas microscopic studies mostly support the acidic air-water interface with the observation of hydronium (H3O+) accumulation in the top layer of the interface. It is crucial to clarify the interfacial preference of OH- and H3O+ ions for rationalizing the debate. In this work, we perform deep potential molecular dynamics simulations to investigate the preferential distribution of OH- and H3O+ ions at the aqueous interfaces. The neural network potential energy surface is trained based on density functional theory calculations with the SCAN functional, which can accurately describe the diffusion of these two ions both in the interface and in the bulk water. In contrast to the previously reported single ion enrichment, we show that both OH- and H3O+ surprisingly prefer to accumulate in interfaces but at different interfacial depths, rendering a double-layer ionic distribution within ∼1 nm near the Gibbs dividing surface. The H3O+ preferentially resides in the topmost layer of the interface, but the OH-, which is enriched in the deeper interfacial layer, has a higher equilibrium concentration due to the more negative free energy of interfacial stabilization [-0.90 (OH-) vs -0.56 (H3O+) kcal/mol]. The present finding of the ionic double-layer distribution may qualitatively offer a self-consistent explanation for the long-term controversy about the acid-base nature of the air-water interface.

Conical Intersections at Interfaces Revealed by Phase-Cycling Interface-Specific Two-Dimensional Electronic Spectroscopy (i2D-ES)

Conical intersections (CIs) hold significant stake in manipulating and controlling photochemical reaction pathways of molecules at interfaces and surfaces by affecting molecular dynamics therein. Currently, there is no tool for characterizing CIs at interfaces and surfaces. To this end, we have developed phase-cycling interface-specific two-dimensional electronic spectroscopy (i2D-ES) and combined it with advanced computational modeling to explore nonadiabatic CI dynamics of molecules at the air/water interface. Specifically, we integrated the phase locked pump pulse pair with an interface-specific electronic probe to obtain the two-dimensional interface-specific responses. We demonstrate that the nonadiabatic transitions of an interface-active azo dye molecule that occur through the CIs at the interface have different kinetic pathways from those in the bulk water. Upon photoexcitation, two CIs are present: one from an intersection of an optically active S2 state with a dark S1 state and the other from the intersection of the progressed S1 with the ground state S0. We find that the molecular conformations in the ground state are different for interfacial molecules. The interfacial molecules are intimately correlated with the locally populated excited state S2 being farther away from the CI region. This leads to slower nonadiabatic dynamics at the interface than in bulk water. Moreover, we show that the nonadiabatic transition from the S1 dark state to the ground state is significantly longer at the interface than that in the bulk, which is likely due to the orientationally restricted configuration of the excited state at the interface. Our findings suggest that orientational configurations of molecules manipulate reaction pathways at interfaces and surfaces.

Lipid Landscapes: Vibrational Spectroscopy for Decoding Membrane Complexity

Cell membranes are incredibly complex environments containing hundreds of components. Despite substantial advances in the past decade, fundamental questions related to lipid-lipid interactions and heterogeneity persist. This review explores the complexity of lipid membranes, showcasing recent advances in vibrational spectroscopy to characterize the structure, dynamics, and interactions at the membrane interface. We include an overview of modern techniques such as surface-enhanced infrared spectroscopy as a steady-state technique with single-bilayer sensitivity, two-dimensional sum-frequency generation spectroscopy, and two-dimensional infrared spectroscopy to measure time-evolving structures and dynamics with femtosecond time resolution. Furthermore, we discuss the potential of multiscale molecular dynamics (MD) simulations, focusing on recently developed simulation algorithms, which have emerged as a powerful approach to interpret complex spectra. We highlight the ongoing challenges in studying heterogeneous environments in multicomponent membranes via current vibrational spectroscopic techniques and MD simulations. Overall, this review provides an up-to-date comprehensive overview of the powerful combination of vibrational spectroscopy and simulations, which has great potential to illuminate lipid-lipid, lipid-protein, and lipid-water interactions in the intricate conformational landscape of cell membranes.

Unified picture of vibrational relaxation of OH stretch at the air/water interface

The elucidation of the energy dissipation process is crucial for understanding various phenomena occurring in nature. Yet, the vibrational relaxation and its timescale at the water interface, where the hydrogen-bonding network is truncated, are not well understood and are still under debate. In the present study, we focus on the OH stretch of interfacial water at the air/water interface and investigate its vibrational relaxation by femtosecond time-resolved, heterodyne-detected vibrational sum-frequency generation (TR-HD-VSFG) spectroscopy. The temporal change of the vibrationally excited hydrogen-bonded (HB) OH stretch band (ν=1→2 transition) is measured, enabling us to determine reliable vibrational relaxation (T1) time. The T1 times obtained with direct excitations of HB OH stretch are 0.2-0.4 ps, which are similar to the T1 time in bulk water and do not noticeably change with the excitation frequency. It suggests that vibrational relaxation of the interfacial HB OH proceeds predominantly with the intramolecular relaxation mechanism as in the case of bulk water. The delayed rise and following decay of the excited-state HB OH band are observed with excitation of free OH stretch, indicating conversion from excited free OH to excited HB OH (~0.9 ps) followed by relaxation to low-frequency vibrations (~0.3 ps). This study provides a complete set of the T1 time of the interfacial OH stretch and presents a unified picture of its vibrational relaxation at the air/water interface.

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.

Adsorption of SARS-CoV-2 Spike (N501Y) RBD to Human Angiotensin-Converting Enzyme 2 at a Lipid/Water Interface

The receptor binding domain (RBD) of spike proteins plays a crucial role in the process of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) attachment to the human angiotensin-converting enzyme 2 (ACE2). The N501Y mutation and later mutations introduced extra positive charges on the spike RBD and resulted in higher transmissibility, likely due to stronger binding with the highly negatively charged ACE2. Consequently, many studies have been devoted to understanding the molecular mechanism of spike protein binding with the ACE2 receptor. Most of the theoretical studies, however, have been done on isolated proteins. ACE2 is a transmembrane protein; thus, it is important to understand the interaction of spike proteins with ACE2 in a lipid matrix. In this study, the adsorption of ACE2 and spike (N501Y) RBD at a lipid/water interface was studied using the heterodyne-detected vibrational sum frequency generation (HD-VSFG) technique. The technique is a non-linear optical spectroscopy which measures vibrational spectra of molecules at an interface and provides information on their structure and orientation. It is found that ACE2 is effectively adsorbed at the positively charged 1,2-dipalmitoyl-3-trimethylammonium-propane (DPTAP) lipid monolayer via electrostatic interactions. The adsorption of ACE2 at the DPTAP monolayer causes a reorganization of interfacial water (D2O) from the D-down to the D-up orientation, indicating that the originally positively charged DPTAP interface becomes negatively charged due to ACE2 adsorption. The negatively charged interface (DPTAP/ACE2) allows further adsorption of positively charged spike RBD. HD-VSFG spectra in the amide I region show differences for spike (N501Y) RBD adsorbed at D2O, DPTAP, and DPTAP/ACE2 interfaces. A red shift observed for the spectra of spike RBD/DPTAP suggests that spike RBD oligomers are formed upon contact with DPTAP lipids.

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na2SO4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Herein, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

Character of the OH Bend-Stretch Combination Band in the Vibrational Spectra of the "Magic" Number H3O+(H2O)20 and D3O+(D2O)20 Cluster Ions

The fundamental transitions that contribute to the diffuse OH stretching spectrum of water are known to increase in width and intensity with increasing red shift from the free OH frequency. In contrast, the profile of the higher-energy combination band involving the OH stretching and the intramolecular HOH bending modes displays a qualitatively different spectral shape with a much faster falloff on the lower-energy side. We elucidate the molecular origin of this difference by analyzing the shapes of the combination bands in the IR spectra of cryogenically cooled H₃O⁺(H₂O)₂₀ and D₃O⁺(D₂O)₂₀ clusters. The difference in the shapes of the bands is traced to differences in the dependence of their transition dipole matrix elements on the hydrogen-bonding environment. The fact that individual transitions across the combination band envelope have similar intensities makes it a useful way to determine the participation of various sites in extended H-bonding networks.

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.

A liquid surface height controller for surface spectroscopy

We present a simple and inexpensive liquid surface height controller that can monitor and maintain the height of a liquid surface in a surface-sensitive experiment. The system is based on a commercial laser pointer, universal serial bus webcam, syringe pump, and homemade control software. The system can sense changes in the height of the surface of ±1 µm, and the maximum range of the device without readjustment is around 2.5 mm. The intended use of the device is to maintain the height of a sample at the air-water interface in a sum-frequency generation spectroscopy measurement, which constantly changes due to water evaporation. A demonstration of the system maintaining the height of a water surface to a tolerance of ±5 µm over a period of 8 h is shown to illustrate the stability of a system controlled by this device.

Progress in phase-sensitive sum frequency generation spectroscopy

Sum frequency generation (SFG) spectroscopy is a unique and powerful tool for investigating surfaces and interfaces at the molecular level. Phase-sensitive SFG (PS-SFG) is an upgraded technique that can overcome the inherent drawbacks of conventional SFG. Here we review several methods of PS-SFG developed and reported in 1990-2020. We introduce how and by which group each PS-SFG method was designed and built in terms of interferometer implementation for optical heterodyne detection, with one exception of a recent numerical method that does not rely on interferometry. We also discuss how PS-SFG solved some typical problems for aqueous interfaces that were once left open by conventional SFG. These problems and their solutions are good examples to demonstrate why PS-SFG is essential. In addition, we briefly note a few terminology issues related with PS-SFG to avoid confusion.

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

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

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

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

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

Mapping the temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, infrared-infrared photodissociation mass spectrometry to the D₃O⁺·(HDO)(D₂O)₁₉ isotopologue of the "magic number" protonated water cluster, H⁺·(H₂O)₂₁ The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90 to 120 K). Starting at 100 K, large "jumps" are observed between two OH frequencies separated by ∼300 cm^-1, indicating migration of the OH group from the bound OH site at 3,350 cm^-1 to the free position at 3,686 cm^-1 Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Spatially dependent H-bond dynamics at interfaces of water/biomimetic self-assembled lattice materials

Understanding hydrogen-bond interactions in self-assembled lattice materials is crucial for preparing such materials, but the role of hydrogen bonds (H bonds) remains unclear. To gain insight into H-bond interactions at the materials' intrinsic spatial scale, we investigated ultrafast H-bond dynamics between water and biomimetic self-assembled lattice materials (composed of sodium dodecyl sulfate and β-cyclodextrin) in a spatially resolved manner. To accomplish this, we developed an infrared pump, vibrational sum-frequency generation (VSFG) probe hyperspectral microscope. With this hyperspectral imaging method, we were able to observe that the primary and secondary OH groups of β-cyclodextrin exhibit markedly different dynamics, suggesting distinct H-bond environments, despite being separated by only a few angstroms. We also observed another ultrafast dynamic reflecting a weakening and restoring of H bonds between bound water and the secondary OH of β-cyclodextrin, which exhibited spatial uniformity within self-assembled domains, but heterogeneity between domains. The restoration dynamics further suggest heterogeneous hydration among the self-assembly domains. The ultrafast nature and meso- and microscopic ordering of H-bond dynamics could contribute to the flexibility and crystallinity of the material--two critically important factors for crystalline lattice self-assemblies--shedding light on engineering intermolecular interactions for self-assembled lattice materials.

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.

In Situ Monitoring of the Unsaturated Phospholipid Monolayer Oxidation in Ambient Air by HD-SFG Spectroscopy

The pulmonary surfactant monolayer is indispensable for the respiratory system. Recently, it was reported that some unsaturated lipids of the pulmonary surfactants are oxidized by low-level ozone in ambient air. However, the molecular-level understanding of the reaction mechanism is still limited due to technical difficulties. We applied heterodyne-detected sum frequency generation (HD-SFG) spectroscopy to probe the reaction process of an unsaturated phospholipid monolayer (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine, POPC), which is one of the major lipids in the pulmonary surfactant, under low-level ozone (30 ± 5 ppb). The HD-SFG spectroscopy realized the accurate peak assignments of the spectra and the identification of molecular species with high sensitivity, which were impossible with previous measurements. The time-resolved spectra indicated that the C═C moiety in the unsaturated alkyl chain is selectively oxidized by ozone with a time constant of 22 ± 3 min by first-order reaction kinetics. Furthermore, it was revealed for the first time that the reaction product of the POPC monolayer under low-level ozone is not the carboxylic form but the aldehyde form based on the vibrational spectroscopy results. The present study has deepened our molecular-level understanding of the oxidation mechanism of unsaturated lipids that are widely found in many biological systems.

Hidden Isolated OH at the Charged Hydrophobic Interface Revealed by Two-Dimensional Heterodyne-Detected VSFG Spectroscopy

Water around hydrophobic groups mediates hydrophobic interactions that play key roles in many chemical and biological processes. Thus, the molecular-level elucidation of the properties of water in the vicinity of hydrophobic groups is important. We report on the structure and dynamics of water at two oppositely charged hydrophobic ion/water interfaces, that is, the tetraphenylborate-ion (TPB^- )/water and tetraphenylarsonium-ion (TPA⁺ )/water interfaces, which are clarified by two-dimensional heterodyne-detected vibrational sum-frequency generation (2D HD-VSFG) spectroscopy. The obtained 2D HD-VSFG spectra of the anionic TPB^- interface reveal the existence of distinct π-hydrogen bonded OH groups in addition to the usual hydrogen-bonded OH groups, which are hidden in the steady-state spectrum. In contrast, 2D HD-VSFG spectra of the cationic TPA⁺ interface only show the presence of usual hydrogen-bonded OH groups. The present study demonstrates that the sign of the interfacial charge governs the structure and dynamics of water molecules that face the hydrophobic region.

Orientation independent vibrational dynamics of lipid-bound interfacial water

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

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.

Heterogenized Molecular Catalysts: Vibrational Sum-Frequency Spectroscopic, Electrochemical, and Theoretical Investigations

Rhenium and manganese bipyridyl tricarbonyl complexes have attracted intense interest for their promising applications in photocatalytic and electrocatalytic CO₂ reduction in both homogeneous and heterogenized systems. To date, there have been extensive studies on immobilizing Re catalysts on solid surfaces for higher catalytic efficiency, reduced catalyst loading, and convenient product separation. However, in order for the heterogenized molecular catalysts to achieve the combination of the best aspects of homogeneous and heterogeneous catalysts, it is essential to understand the fundamental physicochemical properties of such heterogeneous systems, such as surface-bound structures of Re/Mn catalysts, substrate-adsorbate interactions, and photoinduced or electric-field-induced effects on Re/Mn catalysts. For example, the surface may act to (un)block substrates, (un)trap charges, (de)stabilize particular intermediates (and thus affect scaling relations), and shift potentials in different directions, just as protein environments do. The close collaboration between the Lian, Batista, and Kubiak groups has resulted in an integrated approach to investigate how the semiconductor or metal surface affects the properties of the attached catalyst. Synthetic strategies to achieve stable and controlled attachment of Re/Mn molecular catalysts have been developed. Steady-state, time-resolved, and electrochemical vibrational sum-frequency generation (SFG) spectroscopic studies have provided insight into the effects of interfacial structures, ultrafast vibrational energy relaxation, and electric field on the Re/Mn catalysts, respectively. Various computational methods utilizing density functional theory (DFT) have been developed and applied to determine the molecular orientation by direct comparison to spectroscopy, unravel vibrational energy relaxation mechanisms, and quantify the interfacial electric field strength of the Re/Mn catalyst systems. This Account starts with a discussion of the recent progress in determining the surface-bound structures of Re catalysts on semiconductor and Au surfaces by a combined vibrational SFG and DFT study. The effects of crystal facet, length of anchoring ligands, and doping of the semiconductor on the bound structures of Re catalysts and of the substrate itself are discussed. This is followed by a summary of the progress in understanding the vibrational relaxation (VR) dynamics of Re catalysts covalently adsorbed on semiconductor and metal surfaces. The VR processes of Re catalysts on TiO₂ films and TiO₂ single crystals and a Re catalyst tethered on Au, particularly the role of electron-hole pair (EHP)-induced coupling on the VR of the Re catalyst bound on Au, are discussed. The Account also summarizes recent studies in quantifying the electric field strength experienced by the catalytically active site of the Re/Mn catalyst bound on a Au electrode based on a combined electrochemical SFG and DFT study of the Stark tuning of the CO stretching modes of these catalysts. Finally, future research directions on surface-immobilized molecular catalyst systems are discussed.

Polyatomic Iodine Species at the Air-Water Interface and Its Relevance to Atmospheric Iodine Chemistry: An HD-VSFG and Raman-MCR Study

Iodine plays a key role in tropospheric ozone destruction, atmospheric new particle formation, as well as growth. Air-water interface happens to be an important reaction site pertaining to such phenomena. However, except iodide (I^-), the behavior of other iodine species, for example, triiodide (I₃^-) and iodate (IO₃^-, the most abundant iodine species in seawater) at the aqueous interface and their effect on the interfacial water are largely unknown. Using interface-specific vibrational spectroscopy (heterodyne-detected vibrational sum frequency generation), we recorded the imaginary-χ⁽²⁾ spectra (Imχ⁽²⁾; χ⁽²⁾ is the second-order electric susceptibility in OH stretch region) of the air-water interface in the presence of IO₃^-, I₃^-, and I^- (≤0.3 M) in the aqueous subphase. The Imχ⁽²⁾ spectra reveal that the chaotropic I₃^- is the most surface-active anion among the iodine species studied and decreases the vibrational coupling and hydrogen-bonding of interfacial water. Interestingly, the IO₃^-, even being a kosmotrope, is quite prevalent in the interfacial region and preferentially orients the interfacial water as "H-down" (i.e., water dipole moment is pointed toward the bulk water). Mapping of the OH stretch response of ion-affected water at interface (i.e., ΔImχ⁽²⁾ = Imχ⁽²⁾_{air-water-iodine salt} - Imχ⁽²⁾_air-water) with that in the hydration shell of the respective ion (hydration shell water response is obtained by Raman multivariate curve resolution spectroscopy) reveals a correlative link between the ion's influence on the interfacial water and their hydration shell structure. The distinct water structure of stronger as well as weaker H-bonding in the hydration shell of the polyatomic IO₃^- anion promotes the anion to stay at the interfacial region. Thus, the surface prevalence of the iodine species and their effect on the interfacial water are perceived to be crucial for the transfer of iodine from seawater to the atmosphere across the marine boundary layer and the chemistry of iodine at aqueous aerosol surface.

Phase-Sensitive Sum-Frequency Generation Measurements Using a Femtosecond Nonlinear Interferometer

Phase-sensitive sum-frequency spectroscopy is a unique tool to interrogate the vibrational structure of interfaces. A precise understanding of the interfacial structure often relies on accurately determining the phase of χ⁽²⁾, which has recently been demonstrated using a nonlinear interferometer in conjunction with a frequency-scanning picosecond laser system. Here, we implement nonlinear interferometry using a femtosecond laser system for broadband sum-frequency generation. The phase of the vibrational response from a self-assembled monolayer of octadecanethiol on gold is determined using the nonlinear femtosecond interferometer. The results are compared to those obtained using the more traditional heterodyne-detected phase measurements. Both methods give a similar phase spectrum and phase uncertainty. We also discuss the origin of the phase uncertainties and provide guidelines for further improvement.

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.

Interpreting Quasi-Thermal Effects in Ultrafast Spectroscopy of Hydrogen-Bonded Systems

Vibrational excitation of molecules in the condensed phase relaxes through vibrational modes of decreasing energy to ultimately generate an equilibrium state in which the energy is distributed among low-frequency modes. In ultrafast vibrational spectroscopy, changes in the vibrational features of hydrogen-bonded NH and OH stretch modes are typically observed to persist long after these high-frequency vibrations have relaxed. Due to the resemblance to the spectral changes caused by heating the sample, these features are typically described as arising from a hot ground state. However, these spectral features appear on ultrafast time scales that are much too fast to result from a true thermal state, and significant differences between the thermal difference spectrum and the induced quasi-thermal changes in ultrafast spectroscopy are often observed. Here, we examine and directly compare the thermal and quasi-thermal responses of the hydrogen-bonded homodimer of 7-azaindole with temperature-dependent FTIR spectroscopy and ultrafast mid-IR continuum spectroscopy. We find that the thermal difference spectra contain contributions from both dissociation of the hydrogen bonds and from frequency shifts due to changes in the thermal population of low-frequency modes. The transient spectra in ultrafast vibrational spectroscopy are also found to contain two contributions: initial frequency shifts over 2.3 ± 0.11 ps associated with equilibration of the initial excitation, and frequency shifts associated with the excitation of several fingerprint modes, which decay over 21.8 ± 0.11 ps, giving rise to a quasi-thermal response caused by a distribution of fingerprint modes being excited within the sample ensemble. This resembles the thermal frequency shifts due to population changes of low-frequency modes, but not the overall thermal spectrum, which is dominated by features caused by dimer dissociation. These findings provide insight into the changes in the vibrational spectrum from different origins and are important for assigning, analyzing, and comparing features in thermal and ultrafast vibrational spectroscopy of hydrogen-bonded complexes.

The structure of a lanthanide complex at an extractant/water interface studied using heterodyne-detected vibrational sum frequency generation

Eu³⁺ at an extractant/water interface is bound to extractants from the upper side and to water molecules from the lower side, and forms a unique interfacial complex.

Structure at the air/water interface in the presence of phenol: a study using heterodyne-detected vibrational sum frequency generation and molecular dynamics simulation

A simple, neutral organic molecule, phenol, forms a specific hydrogen-bonding structure with water at the air/water interface.

Ultrafast structural molecular dynamics investigated with 2D infrared spectroscopy methods

Ultrafast, multi-dimensional infrared (IR) spectroscopy has been advanced in recent years to a versatile analytical tool with a broad range of applications to elucidate molecular structure on ultrafast timescales, and it can be used for samples in a many different environments. Following a short and general introduction on the benefits of 2D IR spectroscopy, the first part of this chapter contains a brief discussion on basic descriptions and conceptual considerations of 2D IR spectroscopy. Outstanding classical applications of 2D IR are used afterwards to highlight the strengths and basic applicability of the method. This includes the identification of vibrational coupling in molecules, characterization of spectral diffusion dynamics, chemical exchange of chemical bond formation and breaking, as well as dynamics of intra- and intermolecular energy transfer for molecules in bulk solution and thin films. In the second part, several important, recently developed variants and new applications of 2D IR spectroscopy are introduced. These methods focus on (i) applications to molecules under two- and three-dimensional confinement, (ii) the combination of 2D IR with electrochemistry, (iii) ultrafast 2D IR in conjunction with diffraction-limited microscopy, (iv) several variants of non-equilibrium 2D IR spectroscopy such as transient 2D IR and 3D IR, and (v) extensions of the pump and probe spectral regions for multi-dimensional vibrational spectroscopy towards mixed vibrational-electronic spectroscopies. In light of these examples, the important open scientific and conceptual questions with regard to intra- and intermolecular dynamics are highlighted. Such questions can be tackled with the existing arsenal of experimental variants of 2D IR spectroscopy to promote the understanding of fundamentally new aspects in chemistry, biology and materials science. The final part of the chapter introduces several concepts of currently performed technical developments, which aim at exploiting 2D IR spectroscopy as an analytical tool. Such developments embrace the combination of 2D IR spectroscopy and plasmonic spectroscopy for ultrasensitive analytics, merging 2D IR spectroscopy with ultra-high-resolution microscopy (nanoscopy), future variants of transient 2D IR methods, or 2D IR in conjunction with microfluidics. It is expected that these techniques will allow for groundbreaking research in many new areas of natural sciences.

Cooperative Hydrogen-Bond Dynamics at a Zwitterionic Lipid/Water Interface Revealed by 2D HD-VSFG Spectroscopy

Molecular-level elucidation of hydration at biological membrane interfaces is of great importance for understanding biological processes. We studied ultrafast hydrogen-bond dynamics at a zwitterionic phosphatidylcholine/water interface by two-dimensional heterodyne-detected vibrational sum frequency generation (2D HD-VSFG) spectroscopy. The obtained 2D spectra confirm that the anionic phosphate and cationic choline sites are individually hydrated at the interface. Furthermore, the data show that the dynamics of water at the zwitterionic lipid interface is not a simple sum of the dynamics of the water species that hydrate to the separate phosphate and choline. The center line slope (CLS) analysis of the 2D spectra reveals that ultrafast hydrogen-bond fluctuation is not significantly suppressed around the phosphate at the zwitterionic lipid interface, which makes the hydrogen-bond dynamics look similar to that of the bulk water. The present study indicates that the hydrogen-bond dynamics at membrane interfaces is not determined only by the hydrogen bond to a specific site of the interface but is largely dependent on the water dynamics in the vicinity and other nearby moieties, through the hydrogen-bond network.

Phase Acrobatics: The Influence of Excitonic Resonance and Gold Nonresonant Background on Heterodyne-Detected Vibrational Sum Frequency Generation Emission

We show how heterodyne-detected vibrational sum frequency generation (HD-VSFG) spectroscopy can discriminate between the excitonic and monomeric properties of a helical, nanotube molecular aggregate by monitoring the phase of the VSFG emission associated with different polarization configurations. By keeping track of the "phase acrobatics" associated with the added phase of the nonresonant SFG emission of gold as well as that of the double-resonance conditions achieved when the SF frequency is resonant with an electronic exciton transition, we discover that for aggregates of tetra(sulfonatophenyl)porphyrin (TSPP) the PPP-polarized spectra exhibit double-resonance conditions while SSP-polarized spectra exhibit resonance only with the ground-state vibration. Along with observed shifts in the vibrational frequency, intensity differences, and sign flips in the imaginary second-order susceptibility, χ_s,Im⁽²⁾, we conclude that PPP-polarized HD-VSFG spectra reflect the delocalized, excitonic nature of the molecular aggregate, while the SSP-polarized HD-VSFG spectra measure the localized, monomeric nature of the molecular subunits. It is implied from this study that HD-VSFG spectroscopy can be uniquely utilized to measure the excitonic and monomeric properties associated with molecular assemblies for a single sample.

Vibrational Coupling at the Topmost Surface of Water Revealed by Heterodyne-Detected Sum Frequency Generation Spectroscopy

Unraveling vibrational coupling is the key to consistently interpret vibrational spectra of complex molecular systems. The vibrational spectrum of the water surface heavily suffers from vibrational coupling, which hinders complete understanding of the molecular structure and dynamics of the water surface. Here we apply heterodyne-detected sum frequency generation spectroscopy to the water surface and accomplish the assignment of a weak vibrational band located at the lower energy side of the free OH stretch. We find that this band is due to a combination mode of the hydrogen-bonded OH stretch and a low-frequency intermolecular vibration, and this combination band appears in the surface vibrational spectrum through anharmonic vibrational coupling that takes place exclusively at the topmost surface.

Interferometric 2D Sum Frequency Generation Spectroscopy Reveals Structural Heterogeneity of Catalytic Monolayers on Transparent Materials

Molecular monolayers exhibit structural and dynamical properties that are different from their bulk counterparts due to their interaction with the substrate. Extracting these distinct properties is crucial for a better understanding of processes such as heterogeneous catalysis and interfacial charge transfer. Ultrafast nonlinear spectroscopic techniques such as 2D infrared (2D IR) spectroscopy are powerful tools for understanding molecular dynamics in complex bulk systems. Here, we build on technical advancements in 2D IR and heterodyne-detected sum frequency generation (SFG) spectroscopy to study a CO₂ reduction catalyst on nanostructured TiO₂ with interferometric 2D SFG spectroscopy. Our method combines phase-stable heterodyne detection employing an external local oscillator with a broad-band pump pulse pair to provide the first high spectral and temporal resolution 2D SFG spectra of a transparent material. We determine the overall molecular orientation of the catalyst and find that there is a static structural heterogeneity reflective of different local environments at the surface.

Change of the isoelectric point of hemoglobin at the air/water interface probed by the orientational flip-flop of water molecules

Nonlinear vibrational spectroscopy reveals that the isoelectric point of proteins can largely change when the proteins are adsorbed at the air/water interface.

Femtosecond Hydrogen Bond Dynamics of Bulk-like and Bound Water at Positively and Negatively Charged Lipid Interfaces Revealed by 2D HD-VSFG Spectroscopy

Interfacial water in the vicinity of lipids plays an important role in many biological processes, such as drug delivery, ion transportation, and lipid fusion. Hence, molecular-level elucidation of the properties of water at lipid interfaces is of the utmost importance. We report the two-dimensional heterodyne-detected vibrational sum frequency generation (2D HD-VSFG) study of the OH stretch of HOD at charged lipid interfaces, which shows that the hydrogen bond dynamics of interfacial water differ drastically, depending on the lipids. The data indicate that the spectral diffusion of the OH stretch at a positively charged lipid interface is dominated by the ultrafast (<∼100 fs) component, followed by the minor sub-picosecond slow dynamics, while the dynamics at a negatively charged lipid interface exhibit sub-picosecond dynamics almost exclusively, implying that fast hydrogen bond fluctuation is prohibited. These results reveal that the ultrafast hydrogen bond dynamics at the positively charged lipid-water interface are attributable to the bulk-like property of interfacial water, whereas the slow dynamics at the negatively charged lipid interface are due to bound water, which is hydrogen-bonded to the hydrophilic head group.

Bend Vibration of Surface Water Investigated by Heterodyne-Detected Sum Frequency Generation and Theoretical Study: Dominant Role of Quadrupole

Heterodyne-detected vibrational sum frequency generation spectroscopy was applied to the water surface for measuring the imaginary part of second-order nonlinear susceptibility (Im χ((2))) spectrum in the bend frequency region for the first time. The observed Im χ((2)) spectrum shows an overall positive band around 1650 cm(-1), contradicting former theoretical predictions. We further found that the Im χ((2)) spectrum of NaI aqueous solution exhibits an even larger positive band, which is apparently contrary to the flip-flop orientation of surface water. These unexpected observations are elucidated by calculating quadrupole contributions beyond the conventional dipole approximation. It is indicated that the Im χ((2)) spectrum in the bend region has a large quadrupole contribution from the bulk water.