On the role of hydrogen-bond exchanges in the spectral diffusion of water

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

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

A structure-dynamics relationship enables prediction of the water hydrogen bond exchange activation energy from experimental data

It has long been understood that the structural features of water are determined by hydrogen bonding (H-bonding) and that the exchange of, or "jumps" between, H-bond partners underlies many of the dynamical processes in water. Despite the importance of H-bond exchanges there is, as yet, no direct method for experimentally measuring the timescale of the process or its associated activation energy. Here, we identify and exploit relationships between water's structural and dynamical properties that provide an indirect route for determining the H-bond exchange activation energy from experimental data. Specifically, we show that the enthalpy and entropy determining the radial distribution function in liquid water are linearly correlated with the activation energies for H-bond jumps, OH reorientation, and diffusion. Using temperature-dependent measurements of the radial distribution function from the literature, we demonstrate how these correlations allow us to infer the value of the jump activation energy, Ea,0, from experimental results. This analysis gives Ea,0 = 3.43 kcal mol-1, which is in good agreement with that predicted by the TIP4P/2005 water model. We also illustrate other approaches for estimating this activation energy consistent with these estimates.

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.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.

Effects of polarizability and charge transfer on water dynamics and the underlying activation energies

A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy-entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.

Theoretical Two Dimensional Infrared Spectroscopy of Aqueous Solutions of tert-Butyl Alcohol: Variation of the Dynamics of Spectral Diffusion along the Percolation Transition

Binary mixtures of water and tert-butyl alcohol (TBA) are known to exhibit the so-called percolation transition where small clusters of TBA molecules span into large aggregates beyond a threshold concentration of the alcohol. In the present study, we have investigated the linear and two-dimensional infrared spectral features of aqueous solutions of TBA for varying concentration of the alcohol along the percolation transition. The percolation transition is characterized through calculations of intermolecular radial distribution functions and average size of the largest cluster of TBA molecules. It is found that, with variation of alcohol concentration, the radial distribution functions of the central carbon atoms of TBA molecules show a nonmonotonic change in the height of the first peak and also the size of the largest cluster of TBA molecules show a jump in the increase of its size for TBA mole fraction between 0.04 and 0.06 corresponding to a transition from smaller clusters to larger spanning aggregates. However, it is found that the linear infrared spectrum of water does not exhibit any noticeable changes on variation of TBA concentration along the percolation transition. Subsequently, two-dimensional infrared (2DIR) spectra and vibrational frequency time correlation function of water are calculated for all the TBA-water solutions considered in this study. The spectral diffusion of water calculated from 2DIR is found to slow down with increase of the TBA concentration. The time scales of spectral diffusion of water, as characterized by the relaxation of frequency time correlation function, 2DIR metric of central line slope, and also the hydrogen bond time correlation functions, are found to exhibit a noticeable jump along the percolation transition. The hydrophilic group of TBA is found to retard the water dynamics more effectively than the hydrophobic groups. Also, the jump in the dynamical slowdown along the percolation transition is found to be more significant for water molecules at the hydrophilic sites.

Temperature Dependence of Non-Condon Effects in Two-Dimensional Vibrational Spectroscopy of Water

Non-Condon effects in vibrational spectroscopy refers to the dependence of a molecule's vibrational transition dipole and polarizability on the coordinates of the surrounding environment. Earlier studies have shown that such effects can be pronounced for hydrogen-bonded systems like liquid water. Here, we present a theoretical study of two-dimensional vibrational spectroscopy under the non-Condon and Condon approximations at varying temperatures. We have performed calculations of both two-dimensional infrared and two-dimensional vibrational Raman spectra to gain insights into the temperature dependence of non-Condon effects in nonlinear vibrational spectroscopy. The two-dimensional spectra are calculated for the OH vibration of interest in the isotopic dilution limit where the coupling between the oscillators is ignored. Generally, both the infrared and Raman line shapes undergo red shifts with decrease in temperature due to strengthening of hydrogen bonds and decrease in the fraction of OH modes with weaker or no hydrogen bonds. The infrared line shape is further red-shifted under the non-Condon effects at a given temperature, while the Raman line shape does not show any such red shift due to non-Condon effects. The spectral dynamics becomes slower on decrease of temperature due to slower hydrogen bond relaxation and, for a given temperature, the spectral diffusion occurs at a faster rate upon inclusion of non-Condon effects. The time scales of spectral diffusion extracted from different metrics agree well with each other and also with experiments. The changes in the spectrum due to non-Condon effects are found to be more significant at lower temperatures.

Direct calculation of the temperature dependence of 2D-IR spectra: Urea in water

A method for directly calculating the temperature derivative of two-dimensional infrared (2D-IR) spectra from simulations at a single temperature is presented. The approach is demonstrated by application to the OD stretching spectrum of isotopically dilute aqueous (HOD in H₂O) solutions of urea as a function of concentration. Urea is an important osmolyte because of its ability to denature proteins, which has motivated significant interest in its effect on the structure and dynamics of water. The present results show that the temperature dependence of both the linear IR and 2D-IR spectra, which report on the underlying energetic driving forces, is more sensitive to urea concentration than the spectra themselves. Additional physical insight is provided by calculation of the contributions to the temperature derivative from different interactions, e.g., water-water, water-urea, and urea-urea, present in the system. Finally, it is demonstrated how 2D-IR spectra at other temperatures can be obtained from only room temperature simulations.

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH [Formula: see text]. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Water Dynamics and Structure of Highly Concentrated LiCl Solutions Investigated Using Ultrafast Infrared Spectroscopy

In highly concentrated salt solutions, the water hydrogen bond (H-bond) network is completely disrupted by the presence of ions. Water is forced to restructure as dictated by the water-ion and ion-ion interactions. Using ultrafast polarization-selective pump-probe (PSPP) spectroscopy measurements of the OD stretch of dilute HOD, we demonstrate that the limited water-water H-bonding present in concentrated lithium chloride solutions (up to four waters per ion pair) is, on average, stronger than that occurring in bulk water. Furthermore, information on the orientational dynamics and the angular restriction of water H-bonded to both water oxygens and chloride anions was obtained through analysis of the frequency-dependent anisotropy decays. It was found that, when the salt concentration increased, the water showed increasing restriction and slowing at frequencies correlated with strong H-bonding. The angular restriction of the water molecules and strengthening of water-water H-bonds are due to the formation of a water-ion network not present in bulk water and dilute salt solutions. The structural evolution of the ionic medium was also observed through spectral diffusion of the OD stretch using 2D IR spectroscopy. Compared to bulk water, there is significant slowing of the biexponential spectral diffusion dynamics. The slowest component of the spectral diffusion (13 ps) is virtually identical to the time for complete reorientation of HOD measured with the PSPP experiments. This result suggests that the slowest component of the spectral diffusion reflects rearrangement of water molecules in the water-ion network.

Using Activation Energies to Elucidate Mechanisms of Water Dynamics

Recent advances in the calculation of activation energies are shedding new light on the dynamical time scales of liquid water. In this Perspective, we examine how activation energies elucidate the central, but not singular, role of the exchange of hydrogen-bond (H-bond) partners that rearrange the H-bond network of water. The contributions of other motions to dynamical time scales and their associated activation energies are discussed along with one case, vibrational spectral diffusion, where H-bond exchanges are not mechanistically significant. Nascent progress on outstanding challenges, including descriptions of non-Arrhenius effects and activation volumes, are detailed along with some directions for future investigations.

Altered relaxation dynamics of excited state reactions by confinement in reverse micelles probed by ultrafast fluorescence up-conversion

Chemical reactions in confined environments are important in areas as diverse as heterogenous catalysis, environmental chemistry and biochemistry, yet they are much less well understood than the equivalent reactions in either the gas phase or in free solution. The understanding of chemical reactions in solution was greatly enhanced by real time studies of model reactions, through ultrafast spectroscopy (especially when supported by molecular dynamics simulation). Here we review some of the efforts that have been made to adapt this approach to the investigation of reactions in confined media. Specifically, we review the application of ultrafast fluorescence spectroscopy to measure reaction dynamics in the nanoconfined water phase of reverse micelles, as a function of the droplet radius and the charge on the interface. Methods of measurement and modelling of the reactions are outlined. In all of the cases studied (which are focused on ultrafast intramolecular reactions) the effect of confinement was to suppress the reaction. Even in the largest micelles the result in the bulk aqueous phase was not usually recovered, suggesting an important role for specific interactions between reactant and environment, for example at the interface. There was no simple one-to-one correspondence with direct measures of the dynamics of the confined phase. Thus, understanding the effect of confinement on reaction rate appears to require not only knowledge of the dynamics of the reaction in solutions and the effect of confinement on the medium, but also of the interaction between reactant and confining medium.

Identical Water Dynamics in Acrylamide Hydrogels, Polymers, and Monomers in Solution: Ultrafast IR Spectroscopy and Molecular Dynamics Simulations

The dynamics and structure of water in polyacrylamide hydrogels (PAAm-HG), polyacrylamide, and acrylamide solutions are investigated using ultrafast infrared experiments on the OD stretch of dilute HOD/H₂O and molecular dynamics simulations. The amide moiety of the monomer/polymers interacts strongly with water through hydrogen bonding (H-bonding). The FT-IR spectra of the three systems indicate that the range of H-bond strengths is relatively unchanged from bulk water. Vibrational population relaxation measurements show that the amide/water H-bonds are somewhat weaker but fall within the range of water/water H-bond strengths. A previous study of water dynamics in PAAm-HG suggested that the slowing observed was due to increasing confinement with concentration. Here, for the same concentrations of the amide moiety, the experimental results demonstrate that the reorientational dynamics (infrared pump-probe experiments) and structural dynamics (two-dimensional infrared spectroscopy) are identical in the three acrylamide systems studied. Molecular dynamics simulations of the water orientational relaxation in aqueous solutions of the acrylamide monomer, trimer, and pentamer are in good agreement with the experimental results and are essentially chain length independent. The simulations show that there is a slower, low-amplitude (<7%) decay component not accessible by the experiments. The simulations examine the dynamics and structure of water H-bonded to acrylamide, in the first solvent shell, and beyond for acrylamide monomers and short chains. The experiments and simulations show that the slowing of water dynamics in PAAm-HG is not caused by confinement in the polymer network but rather by interactions with individual acrylamide moieties.

Molecular Simulations of Phase Equilibria and Transport Properties in a Model CO2-Expanded Lithium Perchlorate Electrolyte

Carbon-dioxide (CO₂)-expanded liquids, in which a significant mole fraction of CO₂ is dissolved into an organic solvent, have been of significant interest, especially as catalytic support media. Because the CO₂ mole fraction and density can be controlled over a significant range by changing the CO₂ partial pressure, the transport properties of these solvents are highly tunable. Recently, these liquids have garnered interest as potential electrolyte solutions for catalytic electrochemistry; however, little is currently known about the influence of the electrolyte on CO₂ expansion. In the present work, we use molecular-dynamics simulations to study diffusion and viscosity in a model lithium perchlorate electrolyte in CO₂-expanded acetonitrile and demonstrate that these properties are highly dependent on the concentration of the electrolyte. Our present results indicate that the electrolyte slows down diffusion of both CO₂ and MeCN, and that the slowed diffusion in the former is driven by changes in the activation entropy, whereas slowed diffusion in the latter is driven by changes in the activation energy.