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

Molecular dynamics simulation study of water structure and dynamics on the gold electrode surface with adsorbed 4-mercaptobenzonitrile

Understanding water dynamics at charged interfaces is of great importance in various fields, such as catalysis, biomedical processes, and solar cell materials. In this study, we implemented molecular dynamics simulations of a system of pure water interfaced with Au electrodes, on one side of which 4-mercaptobenzonitrile (4-MBN) molecules are adsorbed. We calculated time correlation functions of various dynamic quantities, such as the hydrogen bond status of the N atom of the adsorbed 4-MBN molecules, the rotational motion of the water OH bond, hydrogen bonds between 4-MBN and water, and hydrogen bonds between water molecules in the interface region. Using the Luzar-Chandler model, we analyzed the hydrogen bond dynamics between a 4-MBN and a water molecule. The dynamic quantities we calculated can be divided into two categories: those related to the collective behavior of interfacial water molecules and the H-bond interaction between a water molecule and the CN group of 4-MBN. We found that these two categories of dynamic quantities exhibit opposite trends in response to applied potentials on the Au electrode. We anticipate that the present work will help improve our understanding of the interfacial dynamics of water in various electrolyte systems.

The hydrogen-bonding dynamics of water to a nitrile-functionalized electrode is modulated by voltage according to ultrafast 2D IR spectroscopy

We report the hydrogen-bonding dynamics of water to a nitrile-functionalized and plasmonic electrode surface as a function of applied voltage. The surface-enhanced two-dimensional infrared spectra exhibit hydrogen-bonded and non-hydrogen-bonded nitrile features in similar proportions, plus cross peaks between the two. Isotopic dilution experiments show that the cross peaks arise predominantly from chemical exchange between hydrogen-bonded and non-hydrogen-bonded nitriles. The chemical exchange rate depends upon voltage, with the hydrogen bond of the water to the nitriles breaking 2 to 3 times slower (>63 vs. 25 ps) under a positive as compared to a negative potential. Spectral diffusion created by hydrogen-bond fluctuations occurs on a ~1 ps timescale and is moderately potential-dependent. Timescales from molecular dynamics simulations agree qualitatively with the experiment and show that a negative voltage causes a small net displacement of water away from the surface. These results show that the voltage applied to an electrode can alter the timescales of solvent motion at its interface, which has implications for electrochemically driven reactions.

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.

Electrostatic switch mechanisms of membrane protein trafficking and regulation

Lipid-protein interactions are normally classified as either specific or general. Specific interactions refer to lipid binding to specific binding sites within a membrane protein, thereby modulating the protein's thermal stability or kinetics. General interactions refer to indirect effects whereby lipids affect membrane proteins by modulating the membrane's physical properties, e.g., its fluidity, thickness, or dipole potential. It is not widely recognized that there is a third distinct type of lipid-protein interaction. Intrinsically disordered N- or C-termini of membrane proteins can interact directly but nonspecifically with the surrounding membrane. Many peripheral membrane proteins are held to the cytoplasmic surface of the plasma membrane via a cooperative combination of two forces: hydrophobic anchoring and electrostatic attraction. An acyl chain, e.g., myristoyl, added post-translationally to one of the protein's termini inserts itself into the lipid matrix and helps hold peripheral membrane proteins onto the membrane. Electrostatic attraction occurs between positively charged basic amino acid residues (lysine and arginine) on one of the protein's terminal tails and negatively charged phospholipid head groups, such as phosphatidylserine. Phosphorylation of either serine or tyrosine residues on the terminal tails via regulatory protein kinases allows for an electrostatic switch mechanism to control trafficking of the protein. Kinase action reduces the positive charge on the protein's tail, weakening the electrostatic attraction and releasing the protein from the membrane. A similar mechanism regulates many integral membrane proteins, but here only electrostatic interactions are involved, and the electrostatic switch modulates protein activity by altering the stabilities of different protein conformational states.

Relationships between Molecular Structural Order Parameters and Equilibrium Water Dynamics in Aqueous Mixtures

Water's unique thermophysical properties and how it mediates aqueous interactions between solutes have long been interpreted in terms of its collective molecular structure. The seminal work of Errington and Debenedetti [Nature 2001, 409, 318-321] revealed a striking hierarchy of relationships among the thermodynamic, dynamic, and structural properties of water, motivating many efforts to understand (1) what measures of water structure are connected to different experimentally accessible macroscopic responses and (2) how many such structural metrics are adequate to describe the collective structural behavior of water. Diffusivity constitutes a particularly interesting experimentally accessible equilibrium property to investigate such relationships because advanced NMR techniques allow the measurement of bulk and local water dynamics in nanometer proximity to molecules and interfaces, suggesting the enticing possibility of measuring local diffusivities that report on water structure. Here, we apply statistical learning methods to discover persistent structure-dynamic correlations across a variety of simulated aqueous mixtures, from alcohol-water to polypeptoid-water systems. We investigate a variety of molecular water structure metrics and find that an unsupervised statistical learning algorithm (namely, sequential feature selection) identifies only two or three independent structural metrics that are sufficient to predict water self-diffusivity accurately. Surprisingly, the translational diffusivity of water across all mixed systems studied here is strongly correlated with a measure of tetrahedral order given by water's triplet angle distribution. We also identify a separate small number of structural metrics that well predict an important thermodynamic property, the excess chemical potential of an idealized methane-sized hydrophobe in water. Ultimately, we offer a Bayesian method of inferring water structure by using only structure-dynamics linear regression models with experimental Overhauser dynamic nuclear polarization (ODNP) measurements of water self-diffusivity. This study thus quantifies the relationships among several distinct structural order parameters in water and, through statistical learning, reveals the potential to leverage molecular structure to predict fundamental thermophysical properties. In turn, these findings suggest a framework for solving the inverse problem of inferring water's molecular structure using experimental measurements such as ODNP studies that probe local water properties.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

Elucidating the Molecular Structure of Hydrophobically Modified Polyethylenimine Nanoparticles and Its Potential Implications for DNA Binding

The structural properties of the polyethylenimine (PEI) polymer are generally tuned and selectively modified to reinforce its potential in a broad spectrum of applied domains of medicine, healthcare, material design, sensing, and electronic optimization. The selective modification of the polymer brings about changes in its interfacial characteristics and behavior. The current work involves the synthesis of naphthalimide conjugated polyethylenimine organic nanoparticles (NPEI-ONPs). The interfacial molecular structure of NPEI-ONPs is explored in an aqueous medium at pH 7.4 using surface tensiometry and sum-frequency generation vibrational spectroscopy (SFG-VS). The hydrophobic functionalization rendered a concentration-dependent surface coverage of NPEI-ONPs, where the SFG-VS analysis exhibited the molecular rearrangement of its hydrophobic groups at the interface. The interaction of NPEI-ONPs with double-stranded DNA (dsDNA) is carried out to observe the relevance of the synthesized nanocomposites in the biomedical domain. The bulk-specific studies (i.e., thermal denaturation, viscometry, zeta (ζ) potential, and ATR-FTIR) reveal the condensation of dsDNA in the presence of NPEI-ONPs, making its structure more compact. The interface-sensitive SFG-VS showcased the impact of the dsDNA and NPEI-ONP interaction on the interfacial molecular behavior of NPEI-ONPs at the air-aqueous interface. Our results exhibit the potential of such hydrophobically functionalized ONPs as promising candidates for developing biomedical sealants, substrate coatings, and other biomedical domains.

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.

Importance of Hydrogen Bonding in Crowded Environments: A Physical Chemistry Perspective

Cells are heterogeneous on every length and time scale; cytosol contains thousands of proteins, lipids, nucleic acids, and small molecules, and molecular interactions within this crowded environment determine the structure, dynamics, and stability of biomolecules. For decades, the effects of crowding at the atomistic scale have been overlooked in favor of more tractable models largely based on thermodynamics. Crowding can affect the conformations and stability of biomolecules by modulating water structure and dynamics within the cell, and these effects are nonlocal and environment dependent. Thus, characterizing water's hydrogen-bond (H-bond) networks is a critical step toward a complete microscopic crowding model. This perspective provides an overview of molecular crowding and describes recent time-resolved spectroscopy approaches investigating H-bond networks and dynamics in crowded or otherwise complex aqueous environments. Ultrafast spectroscopy combined with atomistic simulations has emerged as a powerful combination for studying H-bond structure and dynamics in heterogeneous multicomponent systems. We discuss the ongoing challenges toward developing a complete atomistic description of macromolecular crowding from an experimental as well as a theoretical perspective.

Choose your own adventure: Picosecond or broadband vibrational sum-frequency generation spectroscopy

Vibrational sum-frequency generation (VSFG) spectroscopy is a method capable of measuring chemical structure and dynamics within the interfacial region between two bulk phases. At the core of every experimental system is a laser source that influences the experimental capabilities of the VSFG spectrometer. In this article, we discuss the differences between VSFG spectrometers built with picosecond and broadband laser sources as it will impact everything from material costs, experimental build time, experimental capabilities, and more. A focus is placed on the accessibility of the two different SFG systems to newcomers in the SFG field and provides a resource for laboratories considering incorporating VSFG spectroscopy into their research programs. This Tutorial provides a model decision tree to aid newcomers when determining whether the picosecond or femtosecond laser system is sufficient for their research program and navigates through it for a few specific scenarios.

How Far Is "Bulk Water" from Interfaces? Depends on the Nature of the Surface and What We Measure

Using systematic molecular dynamics (MD) simulations, we revisit the question: At what distance from an interface do the properties of "bulk water" get recovered? We have considered three different kinds of interfaces: nonpolar (hydrophobic; isooctane-water interface), charged (negative; AOT bilayer), and polar (zwitterionic; POPC bilayer). In order to interrogate the extent of perturbation of the interfacial water molecules as a function of the distance from the interface, we utilize a diverse range of structural and dynamical parameters. To capture the structural perturbations, we look into local density (translational order), local tetrahedral order parameter, and dipolar orientation of the water molecules. We also explore the anisotropic diffusion of the water molecules in the direction perpendicular to the interface as well as the planar diffusion parallel to the interface in a distance dependent manner. In addition, the orientational time correlation functions have been computed to understand the extent of slowdown in the rotational dynamics. As expected, the electrostatic field emanating from the charged AOT interface seems to have the highest long-range effect on the orientational order and dynamics of the water molecules, whereas specific interactions like hydrogen bonding and electrostatic interaction lead to significant trapping and kinetic slowdown for both AOT and POPC (zwitterionic) very close to the interface. Our analysis highlights that not only the length-scale of perturbation depends on the nature of the interfaces and specific interactions but also the type of water property that we measure/calculate. Different water properties seem to have widely different length-scale of perturbation. Orientational order parameters seem to be perturbed to a much longer length-scale as compared to translational order parameters. The global orientational order of water can be perturbed even up to ∼4-5 nm near the negatively charged AOT surface in the absence of any extra electrolyte. This observation has significant implication toward the interpretation of experimental measurements as well since different spectroscopic techniques would probe different parameters or water properties with possible mutual disagreement and inconsistency between different types of measurements. Thus, our study provides a broader and unifying perspective toward the aspect of "context dependent" structural and dynamical perturbation of "interfacial water".

Direct Evidence for a Surface and Bulk Specific Response in the Sum-Frequency Generation Spectrum of the Water Bend Vibration

We study the bending mode of pure water and charged aqueous surfaces using heterodyne-detected vibrational sum-frequency generation spectroscopy. We observe a low (1626  cm^{-1}) and a high (1656  cm^{-1}) frequency component that can be unambiguously assigned to an interfacial dipole and a bulk quadrupolar response, respectively. We thus demonstrate that probing the bending mode provides structural and quantitative information on both the surface and the bulk.

Water at charged interfaces

The ubiquity of aqueous solutions in contact with charged surfaces and the realization that the molecular-level details of water-surface interactions often determine interfacial functions and properties relevant in many natural processes have led to intensive research. Even so, many open questions remain regarding the molecular picture of the interfacial organization and preferential alignment of water molecules, as well as the structure of water molecules and ion distributions at different charged interfaces. While water, solutes and charge are present in each of these systems, the substrate can range from living tissues to metals. This diversity in substrates has led to different communities considering each of these types of aqueous interface. In this Review, by considering water in contact with metals, oxides and biomembranes, we show the essential similarity of these disparate systems. While in each case the classical mean-field theories can explain many macroscopic and mesoscopic observations, it soon becomes apparent that such theories fail to explain phenomena for which molecular properties are relevant, such as interfacial chemical conversion. We highlight the current knowledge and limitations in our understanding and end with a view towards future opportunities in the field.

Vibrational Sum-Frequency Generation Hyperspectral Microscopy for Molecular Self-Assembled Systems

In this review, we discuss the recent developments and applications of vibrational sum-frequency generation (VSFG) microscopy. This hyperspectral imaging technique can resolve systems without inversion symmetry, such as surfaces, interfaces and noncentrosymmetric self-assembled materials, in the spatial, temporal, and spectral domains. We discuss two common VSFG microscopy geometries: wide-field and confocal point-scanning. We then introduce the principle of VSFG and the relationships between hyperspectral imaging with traditional spectroscopy, microscopy, and time-resolved measurements. We further highlight crucial applications of VSFG microscopy in self-assembled monolayers, cellulose in plants, collagen fibers, and lattice self-assembled biomimetic materials. In these systems, VSFG microscopy reveals relationships between physical properties that would otherwise be hidden without being spectrally, spatially, and temporally resolved. Lastly, we discuss the recent development of ultrafast transient VSFG microscopy, which can spatially measure the ultrafast vibrational dynamics of self-assembled materials. The review ends with an outlook on the technical challenges of and scientific potential for VSFG microscopy.

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.

Versatility of Reverse Micelles: From Biomimetic Models to Nano (Bio)Sensor Design

This paper presents an overview of the principal structural and dynamics characteristics of reverse micelles (RMs) in order to highlight their structural flexibility and versatility, along with the possibility to modulate their parameters in a controlled manner. The multifunctionality in a large range of different scientific fields is exemplified in two distinct directions: a theoretical model for mimicry of the biological microenvironment and practical application in the field of nanotechnology and nano-based sensors. RMs represent a convenient experimental approach that limits the drawbacks of the conventionally biological studies in vitro, while the particular structure confers them the status of simplified mimics of cells by reproducing a complex supramolecular organization in an artificial system. The biological relevance of RMs is discussed in some particular cases referring to confinement and a crowded environment, as well as the molecular dynamics of water and a cell membrane structure. The use of RMs in a range of applications seems to be more promising due to their structural and compositional flexibility, high efficiency, and selectivity. Advances in nanotechnology are based on developing new methods of nanomaterial synthesis and deposition. This review highlights the advantages of using RMs in the synthesis of nanoparticles with specific properties and in nano (bio)sensor design.

Hydrogen bond dynamics of interfacial water molecules revealed from two-dimensional vibrational sum-frequency generation spectroscopy

Vibrational sum-frequency generation (vSFG) spectroscopy allows the study of the structure and dynamics of interfacial systems. In the present work, we provide a simple recipe, based on a narrowband IR pump and broadband vSFG probe technique, to computationally obtain the two-dimensional vSFG spectrum of water molecules at the air-water interface. Using this technique, to study the time-dependent spectral evolution of hydrogen-bonded and free water molecules, we demonstrate that at the interface, the vibrational spectral dynamics of the free OH bond is faster than that of the bonded OH mode.

Water Layer at Hydrophobic Surface: Electrically Dead but Dynamically Alive?

The origin of the anomalous low value of the static dielectric constant (SDC) of confined water has been addressed and unearthed. While the low value is partly due to the different dielectric boundaries, a significant role is played by the "electrically dead layer" (EDL). As the observed dielectric constant is the harmonic mean of the grid-wise SDCs, the first layer, having the smallest SDC, makes a disproportionately large contribution. This enhanced contribution, in turn, arises from the orientationally ordered surface water molecules. They exhibit reduced fluctuations in collective dipole moment, as the molecules remain partly caged due to water-surface interactions. This phenomenon is found to be universal. We study the structure and dynamics of the water molecules which characterize the EDL. We demonstrate that while the EDL remains alive at a molecular level, with a finite residence time, it displays time scales not substantially different compared to the distant water layers.

Resolving the Controversy over Dipole versus Quadrupole Mechanism of Bend Vibration of Water in Vibrational Sum Frequency Generation Spectra

Recently, there has been controversy over whether the HOH bend signal of water in the vibrational sum frequency generation (VSFG) spectrum arises from the conventional dipole mechanism or the quadrupole mechanism. Here, we show that the Im χ⁽²⁾ (the imaginary part of the second-order nonlinear susceptibility) spectra of the HOH bend mode of water at oppositely charged monolayer/water interfaces all exhibit positive bands, irrespective of the difference in the sign of the charge at the interface. Furthermore, it is found that the peak frequency of the HOH bend band substantially changes depending on the chemical structure of the charged headgroup located at the interface. These results demonstrate that the VSFG signal of the HOH bend vibration is generated from interfacial water with the interfacial quadrupole mechanism that is associated with the large field gradient of incident lights localized in a very thin region at the interface.

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.

Interfacial Vibrational Dynamics of Ice Ih and Liquid Water

Insights into energy flow dynamics at ice surfaces are essential for understanding chemical dynamics relevant to atmospheric and geographical sciences. Here, employing ultrafast surface-specific spectroscopy, we report the interfacial vibrational dynamics of ice I_h. A comparison to liquid water surfaces reveals accelerated vibrational energy relaxation and dissipation at the ice surface for hydrogen-bonded OH groups. In contrast, free-OH groups sticking into the vapor phase exhibit substantially slower vibrational dynamics on ice. The acceleration and deceleration of vibrational dynamics of these different OH groups at the ice surface are attributed to enhanced intermolecular coupling and reduced rotational mobility, respectively. Our results highlight the unique properties of free-OH groups on ice, putatively linked to the high catalytic activities of ice surfaces.

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.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

Dynamic coupling of a hydration layer to a fluid phospholipid membrane: intermittency and multiple time-scale relaxations

Understanding the coupling of a hydration layer and a lipid membrane is crucial to gaining access to membrane dynamics and understanding its functionality towards various biological processes.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

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.

Hydrogen bonding and molecular orientations across thin water films on sapphire

HYPOTHESIS

Water vapor binding to metal oxide surfaces produces thin water films with properties controlled by interactions with surface hydroxo sites. Hydrogen bonding populations vary across films and induce different molecular orientations than at the surface of liquid water. Identifying these differences can open possibilities for tailoring film-mediated catalytic reactions by choice of the supporting metal oxide substrate.

EXPERIMENTS

The (0001) face of a single sapphire (α-Al₂O₃) sample exposed to water vapor and the surface of liquid water were probed by polarization dependent Sum Frequency Generation-Vibration Spectroscopy (SFG-VS). Molecular dynamics (MD) provided insight into the hydrogen bond populations and molecular orientations across films and liquid water.

FINDINGS

SFG-VS revealed a submonolayer film on sapphire exposed to 43% relative humidity (R.H.), and a multilayer film at 78% R.H. Polarization dependent SFG-VS spectra showed that median tilt angles of free OH bonds on the top of films are at ∼43° from the normal of the (0001) face but at 38° on neat liquid water. These values align with MD simulations, which also show that up to 36% of all OH bonds on films are free. This offers new means for understanding how interfacial reactions on sapphire-supported water films could contrast with those involving liquid water.

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.

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.

Polar Interactions Play an Important Role in the Energetics of the Main Phase Transition of Phosphatidylcholine Membranes

Conformational changes of membrane proteins are accompanied by deformation in the surrounding lipid bilayer. To gain insight into the energetics of membrane deformation, the phase behavior of dimyristoylphosphatidylcholine (DMPC) membranes in the presence of the dipole potential, ψ_d, modifiers was investigated by differential scanning calorimetry. 7-Ketocholesterol, which weakens ψ_d and reduces membrane-perpendicular dipole-dipole repulsion, causes a discrete second peak on the high-temperature side of the main transition, whereas 6-ketocholestanol, which strengthens ψ_d and increases membrane-perpendicular dipole-dipole repulsion, merely produces a shoulder. Measurements on pure DMPC vesicles showed that the observed temperature profile could not be explained by a single endothermic process, that is, breaking of van der Waals forces between hydrocarbon chains alone. Removal of NaCl from the buffer caused an increase in the main transition temperature and the appearance of an obvious shoulder, implicating polar interactions. Consideration of the phosphatidylcholine (PC) head group dipole moment indicates direct interactions between PC dipoles that are unlikely to account for the additional process. It seems more likely that the breaking of an in-plane hydrogen-bonded network consisting of hydrating water dipoles together with zwitterionic lipid head groups is responsible. The evidence presented supports the idea that the breaking of van der Waals forces between lipid tails required for the main phase transition of PC membranes is coupled to partial breaking of a hydrogen-bonded network at the membrane surface.

Simulations of Nanoseparated Charged Surfaces Reveal Charge-Induced Water Reorientation and Nonadditivity of Hydration and Mean-Field Electrostatic Repulsion

We perform atomistic simulations of nanometer-separated charged surfaces in the presence of monovalent counterions at fixed water chemical potential. The counterion density profiles are well described by a modified Poisson-Boltzmann (MPB) approach that accounts for nonelectrostatic ion-surface interactions, while the effects of smeared-out surface-charge distributions and dielectric profiles are found to be relatively unimportant. The simulated surface interactions are for weakly charged surfaces well described by the additive contributions of hydration and MPB repulsions, but already for a moderate surface charge density of σ = -0.77 e/nm² this additivity breaks down. This we rationalize by a combination of different effects, namely, counterion correlations as well as the surface charge-induced reorientation of hydration water, which modifies the effective water dielectric constant as well as the hydration repulsion.

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

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