Water Dynamics: Vibrational Echo Correlation Spectroscopy and Comparison to Molecular Dynamics Simulations

Interplay between anion-receptor and anion-solvent interactions in halide receptor complexes characterized with ultrafast infrared spectroscopies

The competition between host-guest binding and solvent interactions is a crucial factor in determining the binding affinities and selectivity of molecular receptor species. The interplay between these competing interactions, however, have been difficult to disentangle. In particular, the development of molecular-level descriptions of solute-solvent interactions remains a grand experimental challenge. Herein, we investigate the prototypical halide receptor meso-octamethylcalix[4]pyrrole (OMCP) complexed with either chloride or bromide anions in both dichloromethane (DCM) and chloroform (trichloromethane, TCM) solvent using ultrafast infrared transient absorption and 2D IR spectroscopies. OMCP·Br- complexes in both solvents display slower vibrational relaxation dynamics of the OMCP pyrrole NH stretches, consistent with weaker H-bonding interactions with OMCP compared to chloride and less efficient intermolecular relaxation to the solvent. Further, OMCP·Br- complexes show nearly static spectral diffusion dynamics compared to OMCP·Cl-, indicating larger structural fluctuations occur within chloride complexes. Importantly, distinct differences in the vibrational spectra and dynamics are observed between DCM and TCM solutions. The data are consistent with stronger and more perturbative solvent effects in TCM compared to DCM, despite DCM's larger dielectric constant and smaller reported OMCP·X- binding affinities. These differences are attributed to the presence of weak H-bond interactions between halides and TCM, in addition to competing interactions from the bulky tetrabutylammonium countercation. The data provide important experimental benchmarks for quantifying the role of solvent and countercation interactions in anion host-guest complexes.

Effects of Nanoconfinement on Dynamics in Concentrated Aqueous Magnesium Chloride Solutions

Water behavior in various natural and manufactured settings is influenced by confinement in organic or inorganic frameworks and the presence of solutes. Here, the effects on dynamics from both confinement and the addition of solutes are examined. Specifically, water and ion dynamics in concentrated (2.5-4.2 m) aqueous magnesium chloride solutions confined in mesoporous silica (2.8 nm pore diameter) were investigated using polarization selective pump-probe and 2D infrared spectroscopies. Fitting the rotational and spectral diffusion dynamics measured by the vibrational probe, selenocyanate, with a previously developed two-state model revealed distinct behaviors at the interior of the silica pores (core state) and near the wall of the confining framework (shell state). The shell dynamics are noticeably slower than the bulk, or core, dynamics. The concentration-dependent slowing of the dynamics aligns with behavior in the bulk solutions, but the spectrally separated water-associated and Mg2+-associated forms of the selenocyanate probe exhibit different responses to confinement. The disparity in the complete reorientation times is larger upon confinement, but the spectral diffusion dynamics become more similar near the silica surface. The length scales that characterize the transition from surface-influenced to bulk-like behavior for the salt solutions in the pores are discussed and compared to those of pure water and an organic solvent confined in the same pores. These comparisons offer insights into how confinement modulates the properties of different liquids.

Fast Structural Dynamics in Concentrated HCl Solutions: From Proton Hopping to the Bulk Viscosity

Concentrated acid solutions, particularly HCl, have been studied extensively to examine the proton hopping and infrared spectral signatures of hydronium ions. Much less attention has been given to the structural dynamics of concentrated HCl solutions. Here, we apply optical heterodyne detected-optical Kerr effect (OHD-OKE) measurements to examine HCl concentration-dependent dynamics from moderate (0.8 m) to very high (15.5 m) concentrations and compare the results to the dynamics of NaCl solutions, as Na+ is similar in size to the hydronium cation. Both HCl and NaCl OHD-OKE signals decay as triexponentials at all concentrations, in contrast to pure water, which decays as a biexponential. Two remarkable features of the HCl dynamics are the following: (1) the bulk viscosity is linearly related to the slowest decay constant, t3, and (2) the concentration-dependent proton hopping times, determined by ab initio MD simulations and 2D IR chemical exchange experiments, both obtained from the literature, fall on the same line as the slowest structural dynamics relaxation time, t3, within experimental error. The structural dynamics of hydronium/chloride/water clusters, with relaxation times t3, are responsible for the concentration dependence of microscopic property of proton hopping and the macroscopic bulk viscosity. The slowest time constant (t3), which does not have a counterpart in pure water, is 3 ps at 0.8 m and increases by a factor of ∼2 by 15.5 m. The two fastest HCl decay constants, t1 and t2, are similar to those of pure water and increase mildly with the concentration.

Dynamics of Formamide-Water Mixtures Investigated by Linear and Nonlinear Infrared Spectroscopy

Formamide (FA) exhibits complete miscibility with water, offering a simplified model for exploring the solvation dynamics of peptide linkages in biophysical processes. Its liquid state demonstrates a three-dimensional hydrogen bonding network akin to water, reflecting solvent-like behavior. Analyzing the microscopic structure and dynamics of FA-water mixtures is expected to provide crucial insights into hydrogen bonding dynamics─a key aspect of various biophysical phenomena. This study is focused on the dynamics of FA-water mixtures using linear and femtosecond infrared spectroscopies. By using the intrinsic OD stretch and extrinsic probe SCN-, the local vibrational behaviors across various FA-water compositions were systematically investigated. The vibrational relaxation of OD stretch revealed a negligible impact of FA addition on the vibrational lifetime of water molecules, underscoring the mixture's water-like behavior. However, the reorientational dynamics of OD stretch slowed with increasing FA mole fraction (XFA), plateauing beyond XFA > 0.5. This suggests a correlation between OD's reorientational time and the strength of the hydrogen bond network, likely tied to the solution's changing dielectric constant. Conversely, the vibrational relaxation dynamics of SCN- was strongly correlated with XFA, highlighting a competition between water and FA molecules in solvating SCN-. Moreover, a linear relationship between rising viscosity and the prolonged correlation time of SCN-'s slow dynamics indicates that the solution's macroscopic viscosity is dictated by the extended structures formed between FA and water molecules. The relation between the reorientation dynamics of the SCN- and the macroscopic viscosity in aqueous FA-water mixture solutions was analyzed by using the Stokes-Einstein-Debye equations. The direct viscosity-diffusion coupling is observed, which can be attributed to the homogeneous dynamics feature in FA-water mixture solutions. The inclusion of these intrinsic and extrinsic probes not only enhances the comprehensiveness of our analysis but also provides valuable insights into various aspects of the dynamics within the FA-water system. This investigation sheds light on the fundamental dynamics of FA-water mixtures, emphasizing their molecular-level homogeneity in this binary mixture solution.

Water-in-Salt: Fast Dynamics, Structure, Thermodynamics, and Bulk Properties

We present concentration-dependent dynamics of highly concentrated LiBr solutions and LiCl temperature-dependent dynamics for two high concentrations and compare the results to those of prior LiCl concentration-dependent data. The dynamical data are obtained using ultrafast optical heterodyne-detected optical Kerr effect (OHD-OKE). The OHD-OKE decays are composed of two pairs of biexponentials, i.e., tetra-exponentials. The fastest decay (t1) is the same as pure water's at all concentrations within error, while the second component (t2) slows slightly with concentration. The slower components (t3 and t4), not present in pure water, slow substantially, and their contributions to the decays increase significantly with increasing concentration, similar to LiCl solutions. Simulations of LiCl solutions from the literature show that the slow components arise from large ion/water clusters, while the fast components are from ion/water structures that are not part of large clusters. Temperature-dependent studies (15-95 °C) of two high LiCl concentrations show that decreasing the temperature is equivalent to increasing the room temperature concentration. The LiBr and LiCl concentration dependences and the two LiCl concentrations' temperature dependences all have bulk viscosities that are linearly dependent on τcslow, the correlation time of the slow dynamics (weighted averages of t3 and t4). Remarkably, all four viscosity vs 1/τCslow plots fall on the same line. Application of transition state theory to the temperature-dependent data yields the activation enthalpies and entropies for the dynamics of the large ion/water clusters, which underpin the bulk viscosity.

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.

Cutting through the Noise: Extracting Dynamics from Ultrafast Spectra Using Dynamic Mode Decomposition

Coherent multidimensional spectroscopy provides experimental access to molecular structure and subpicosecond dynamics in solution. Dynamics are typically inferred from the evolution of lineshapes over a function of waiting time. Numerous spectral analysis methods, such as center/nodal line slope, have been developed to extract these dynamics. However, the extracted dynamics can depend heavily on subjective choices, such as the region selected for CLS analysis or the chosen models. In this study, we introduce a novel approach to extracting dynamics from ultrafast two-dimensional infrared (2D IR) spectra by using dynamic mode decomposition (DMD). As a data-driven method, DMD directly extracts spatiotemporal structures from the complex 2D IR spectra. We evaluated the performance of DMD in simulated and experimental spectra containing overlapped peaks. We show that DMD can retrieve the dynamics of overlapped transitions and cross peaks that are typically challenging to extract with traditional methods. In addition, we demonstrate that combining conditional generative adversarial neural networks with DMD can recover dynamics even at low signal-to-noise ratios. DMD methods do not require preliminary assumptions and can be readily extended to other multidimensional spectroscopies.

Simulating the isotropic Raman spectra of O-H stretching mode in liquid H2O based on a machine learning potential: the influence of vibrational couplings

In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.

Structural evolution of water-in-propylene carbonate mixtures revealed by polarized Raman spectroscopy and molecular dynamics

The liquid structure of systems wherein water is limited in concentration or through geometry is of great interest in various fields such as biology, materials science, and electrochemistry. Here, we present a combined polarized Raman and molecular dynamics investigation of the structural changes that occur as water is added incrementally to propylene carbonate (PC), a polar, aprotic solvent that is important in lithium-ion batteries. Polarized Raman spectra of PC solutions were collected for water mole fractions 0.003 ≤ χwater ≤ 0.296, which encompasses the solubility range of water in PC. The novel approach taken herein provides additional hydrogen bond and solvation characterization of this system that has not been achievable in previous studies. Analysis of the polarized carbonyl Raman band in conjunction with simulations demonstrated that the bulk structure of the solvent remained unperturbed upon the addition of water. Experimental spectra in the O-H stretching region were decomposed through Gaussian fitting into sub-bands and comparison to studies of dilute HOD in D2O. With the aid of simulations, we identified these different bands as water arrangements having different degrees of hydrogen bonding. The observed water structure within PC indicates that water tends to self-aggregate, forming a hydrogen bond network that is distinctly different from the bulk and dependent on concentration. For example, at moderate concentrations, the most likely aggregate structures are chains of water molecules, each with two hydrogen bonds.

Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces

The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.

Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy

The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.

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.

Dynamics of Concentrated Aqueous Lithium Chloride Solutions Investigated with Optical Kerr Effect Experiments

We report the dynamics of concentrated lithium chloride aqueous solutions over a range of moderate to high concentrations. Concentrations (1-29 to 1-3.3 LiCl-water) were studied in which, at the highest concentrations, there are far too few water molecules to solvate the ions. The measurements were made with optically heterodyne-detected optical Kerr effect experiments, a non-resonant technique able to observe dynamics over a wide range of time scales and signal amplitudes. While the pure water decay is a biexponential, the LiCl-water decays are tetra-exponentials at all concentrations. The faster two decays arise from water dynamics, while the slower two decays reflect the dynamics of the ion-water network. The fastest decay (t₁) is the same as pure water at all concentrations. The second decay (t₂) is also the same as that of pure water at the lower concentrations, and then, it slows with increasing concentration. The slower dynamics (t₃ and t₄), which do not have counterparts in pure water, arise from ion-water complexes and, at the highest concentrations, an extended ion-water network. Comparisons are made between the concentration dependence of the observed dynamics and simulations of structural changes from the literature, which enable the assignment of dynamics to specific ion-water structures. The concentration dependences of the bulk viscosity and the ion-water network dynamics are directly correlated. The correlation provides an atomistic-level understanding of the viscosity.

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.

Dynamics of Acrylamide Hydrogels, Polymers, and Monomers in Water Measured with Optical Heterodyne-Detected Optical Kerr Effect Spectroscopy

The ultrafast dynamics of acrylamide monomers (AAm), polyacrylamide (PAAm), and polyacrylamide hydrogels (PAAm-HG) in water were studied using optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy. Previous ultrafast infrared (IR) measurements of the water dynamics showed that at the same concentration of the acrylamide moiety, AAm, PAAm, and PAAm-HG exhibited identical water dynamics and that these dynamics slowed with increasing concentration. In contrast to the IR measurements, OHD-OKE experiments measure the dynamics of both the water and the acrylamide species, which occur on different time scales. In this study, the dynamics of all the acrylamide systems slowed with increasing concentration. We found that AAm exhibits tetraexponential decays, the longest component of which followed Debye-Stokes-Einstein behavior except for the highest concentration, 40% (w/v). Low concentrations of PAAm followed a single power law decay, while high concentrations of PAAm and all concentrations of PAAm-HG decayed with two power laws. The highest concentrations, 25% and 40%, of PAAm and PAAm-HG showed nearly identical dynamics. We interpreted this result as reflecting a similar extent of chain-chain interactions. At low concentrations, PAAm displays non-Markovian, single-chain dynamics (single power law), but PAAm displays entangled chain-chain interactions at high concentrations (two power laws). PAAm-HG has chain-chain interactions at all concentrations that arise from the cross-linking. At high concentrations, the dynamics of the entangled of PAAm become identical within error as those of the cross-linked PAAm-HG.

Direct Correlation between Short-Range Vibrational Spectral Diffusion and Localized Ion-Cage Dynamics of Water-in-Salt Electrolytes

The molecular dynamics simulations of a "water-in-salt" electrolyte, lithium bis(trifluoromethyl sulfonyl) imide (LiNTf₂), with a varying concentration range of 3 to 20 m were performed to establish a direct connection between a dynamic property like the ion-cage lifetime with the short-range vibrational stretching frequency shift of the used probe, HOD. The properties reported here are compared to that obtained from experiments performed at the same concentrations. The time-series wavelet transform was adopted as a preferable mathematical tool for calculating the instantaneous fluctuating frequencies of the probe O-D stretch mode and the concentration-dependent vibrational stretch spectral signature based on the variable functions associated with a particular chemical bond derived from classical molecular dynamics trajectories. The decay time constants of frequency fluctuations and the lifetime of the ion cage (τ_IC) were estimated as a function of salt concentration. Herein, we emphasize the correlation between the slowest time constant (τ₃) of the decay of O-D stretch frequency fluctuations and the timescales associated with the lifetime of ion cages (τ_IC). The results exhibit that the existing relationships were also concentration-dependent. Therefore, this study highlights the connection between the ionic motions that regulate the overall system dynamics with the short-range vibrational frequency shift of the used probe, which was used similar to experiments. It also provides an understanding of the interionic interactions and the dynamical and spectral properties of the electrolytic mixtures. We establish a direct correlation between short-range frequency profile and localized ion-cage lifetime, which can fill the gap of understanding between viscosity, vibrational frequency, and ion-cage dynamics of electrolytes.

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy (2DIR) was applied to phylloquinone (PhQ), an important biological cofactor, to elucidate the impact of hydrogen bonding on the ultrafast dynamics and energetics of the carbonyl stretching modes. 2DIR measurements were performed on PhQ dissolved in hexanol, which served as the hydrogen bonding solvent, and hexane, which served as a non-hydrogen bonding control. Molecular dynamics simulations and quantum chemical calculations were performed to aid in spectral assignment and interpretation. From the position of the peaks in the 2DIR spectra, we extracted the transition frequencies for the fundamental, overtone, and combination bands of hydrogen bonded and non-hydrogen bonded carbonyl groups of PhQ in the 1635-1680 cm^-1 region. We find that hydrogen bonding to a single carbonyl group acts to decouple the two carbonyl units of PhQ. Through analysis of the time-resolved 2DIR data, we find that hydrogen bonding leads to faster vibrational relaxation as well as an increase in the inhomogeneous broadening of the carbonyl groups. Overall, this work demonstrates how hydrogen bonding to the carbonyl groups of PhQ presents in the 2DIR spectra, laying the groundwork to use PhQ as a 2DIR probe to characterize the ultrafast fluctuations in the local environment of natural photosynthetic complexes.

Distinct Hydrogen Bonding Dynamics Underlies the Microheterogeneity in DMF-Water Mixtures

The hydrogen bonding interaction between the amide functional group and water is fundamental to understanding the liquid-liquid heterogeneity in biological systems. Herein, the structure and dynamics of the N,N-dimethylformamide (DMF)-water mixtures have been investigated by linear and nonlinear IR spectroscopies, using the hydroxyl stretch and extrinsic probe of thiocyanate as local vibrational reporters. According to vibrational relaxation dynamics measurements, the orientational dynamics of water is not directly tied to those of DMF molecules. Wobbling-in-a-cone analysis demonstrates that the water molecules have varying degrees of angular restriction depending on their composition due to the formation of specific water-DMF networks. Because of the preferential solvation by DMF molecules, the rotational dynamics of the extrinsic probe is slowed significantly, and its rotational time constants are correlated to the change of solution viscosity. The unique structural dynamics observed in the DMF-water mixtures is expected to provide important insights into the underlying mechanism of microscopic heterogeneity in binary mixtures.

Dephasing Processes in the Molecular Dye Lumogen-F Orange Characterized by Two-Dimensional Electronic Spectroscopy

Molecular dyes are finding more and more applications in photonics and quantum technologies, such as polaritonic optical microcavities, organic quantum batteries and single-photon emitters for quantum sensing and metrology. For all these applications, it is of crucial importance to characterize the dephasing mechanisms. In this work we use two-dimensional electronic spectroscopy (2DES) to study the temperature dependent dephasing processes in the prototypical organic dye Lumogen-F orange. We model the 2DES maps using the Bloch equations for a two-level system and obtain a dephasing time T₂ = 53 fs at room temperature, which increases to T₂ = 94 fs at 86 K. Furthermore, spectral diffusion processes are observed and modeled by a combination of underdamped and overdamped Brownian oscillators. Our results provide useful design parameters for advanced optoelectronic and photonic devices incorporating dye molecules.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

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.

Cationic Effects on the Structural Dynamics of the Metal Ion-Crown Ether Complexes Investigated by Ultrafast Infrared Spectroscopy

It is usually believed that the binding affinity and selectivity of an alkali metal ion with crown ether are defined by the size matching model. However, the underlying mechanism of the specific host-guest interactions and the structural dynamics of the metal ions confined in the cavity of the crown ethers in the solutions are still not clear. In this report, a series of alkali thiocyanate salts (XSCN; X = Li, Na, K, and Cs) complexed with 18-crown-6 (a typical crown ether) in the chloroform solutions were studied by the polarization-selective infrared pump-probe spectroscopy and the ultrafast two-dimensional infrared (2D IR) spectroscopy. The SCN^- counteranions were employed as the local vibrational probe to reveal the specific host-guest interactions in the crown ether complexes. The rotational dynamics and spectral diffusion of SCN^- vibration were both measured by ultrafast IR spectroscopy, and it was found that the metal cations hosted by the crown ethers can have a pronounced effect on the rotational dynamics of the counteranions. The reorientational time constants of the SCN^- vibration in the complexation follow the order Li⁺ > Na⁺ > K⁺ ≃ Cs⁺. More importantly, the spectral diffusion dynamics of SCN^-, which quantifies the decay of the correlation of the frequency fluctuations in the complexation, was also affected by the metal ions but showed a different order of cationic effect. A detailed analysis of the 2D IR data showed that the spectral diffusion of SCN^- counteranion clearly decayed with two different time scales in the complex of 18-crown-6 with K⁺. The 3-4-fold slowdown in spectral diffusion indicated that the fluctuation of SCN^- vibrational transition frequency was strongly affected by the K⁺ cation due to the geometric constraint imposed by the crown ether. The results should help the researchers to unravel the specific host-guest interactions and further reveal the origination of the binding selectivity of crown ether for metal cations in the condensed phases from the perspective of structural dynamics.

Least-Squares Fitting of Multidimensional Spectra to Kubo Line-Shape Models

We report a comprehensive study of the efficacy of least-squares fitting of multidimensional spectra to generalized Kubo line-shape models and introduce a novel least-squares fitting metric, termed the scale invariant gradient norm (SIGN), that enables a highly reliable and versatile algorithm. The precision of dephasing parameters is between 8× and 50× better for nonlinear model fitting compared to that for the centerline-slope (CLS) method, which effectively increases data acquisition efficiency by 1–2 orders of magnitude. Whereas the CLS method requires sequential fitting of both the nonlinear and linear spectra, our model fitting algorithm only requires nonlinear spectra but accurately predicts the linear spectrum. We show an experimental example in which the CLS time constants differ by 60% for independent measurements of the same system, while the Kubo time constants differ by only 10% for model fitting. This suggests that model fitting is a far more robust method of measuring spectral diffusion than the CLS method, which is more susceptible to structured residual signals that are not removable by pure solvent subtraction. Statistical analysis of the CLS method reveals a fundamental oversight in accounting for the propagation of uncertainty by Kubo time constants in the process of fitting to the linear absorption spectrum. A standalone desktop app and source code for the least-squares fitting algorithm are freely available, with example line-shape models and data. We have written the MATLAB source code in a generic framework where users may supply custom line-shape models. Using this application, a standard desktop fits a 12-parameter generalized Kubo model to a 10⁶ data-point spectrum in a few minutes.

Proton Transfer from a Photoacid to a Water Wire: First Principles Simulations and Fast Fluorescence Spectroscopy

Proton transfer reactions are ubiquitous in chemistry, especially in aqueous solutions. We investigate photoinduced proton transfer between the photoacid 8-hydroxypyrene-1,3,6-trisulfonate (HPTS) and water using fast fluorescence spectroscopy and ab initio molecular dynamics simulations. Photoexcitation causes rapid proton release from the HPTS hydroxyl. Previous experiments on HPTS/water described the progress from photoexcitation to proton diffusion using kinetic equations with two time constants. The shortest time constant has been interpreted as protonated and photoexcited HPTS evolving into an "associated" state, where the proton is "shared" between the HPTS hydroxyl and an originally hydrogen bonded water. The longer time constant has been interpreted as indicating evolution to a "solvent separated" state where the shared proton undergoes long distance diffusion. In this work, we refine the previous experimental results using very pure HPTS. We then use excited state ab initio molecular dynamics to elucidate the detailed molecular mechanism of aqueous excited state proton transfer in HPTS. We find that the initial excitation results in rapid rearrangement of water, forming a strong hydrogen bonded network (a "water wire") around HPTS. HPTS then deprotonates in ≤3 ps, resulting in a proton that migrates back and forth along the wire before localizing on a single water molecule. We find a near linear relationship between the emission wavelength and proton-HPTS distance over the simulated time scale, suggesting that the emission wavelength can be used as a ruler for the proton distance. Our simulations reveal that the "associated" state corresponds to a water wire with a mobile proton and that the diffusion of the proton away from this water wire (to a generalized "solvent-separated" state) corresponds to the longest experimental time constant.

Machine Learning Approach for Describing Water OH Stretch Vibrations

A machine learning approach employing neural networks is developed to calculate the vibrational frequency shifts and transition dipole moments of the symmetric and antisymmetric OH stretch vibrations of a water molecule surrounded by water molecules. We employed the atom-centered symmetry functions (ACSFs), polynomial functions, and Gaussian-type orbital-based density vectors as descriptor functions and compared their performances in predicting vibrational frequency shifts using the trained neural networks. The ACSFs perform best in modeling the frequency shifts of the OH stretch vibration of water among the types of descriptor functions considered in this paper. However, the differences in performance among these three descriptors are not significant. We also tried a feature selection method called CUR matrix decomposition to assess the importance and leverage of the individual functions in the set of selected descriptor functions. We found that a significant number of those functions included in the set of descriptor functions give redundant information in describing the configuration of the water system. We here show that the predicted vibrational frequency shifts by trained neural networks successfully describe the solvent-solute interaction-induced fluctuations of OH stretch frequencies.

Coherent Acoustic Interferometry during the Photodriven Oxygen Evolution Reaction Associates Strain Fields with the Reactive Oxygen Intermediate (Ti-OH*)

The oxygen evolution reaction (OER) from water requires the formation of metastable, reactive oxygen intermediates to enable oxygen-oxygen bond formation. Conversely, such reactive intermediates could also structurally modify the catalyst. A descriptor for the overall catalytic activity, the first electron and proton transfer OER intermediate from water, (M-OH*), has been associated with significant distortions of the metal-oxygen bonds upon charge-trapping. Time-resolved spectroscopy of in situ, photodriven OER on transition metal oxide surfaces has characterized M-OH* for the charge trapping and the symmetry of the lattice distortions by optical and vibrational transitions, respectively, but had yet to detect an interfacial strain field arising from a surface coverage M-OH*. Here, we utilize picosecond, coherent acoustic interferometry to detect the uniaxial strain normal to the SrTiO₃/aqueous interface directly caused by Ti-OH*. The spectral analysis applies a fairly general methodology for detecting a combination of the spatial extent, magnitude, and generation time of the interfacial strain through the coherent oscillations' phase. For lightly n-doped SrTiO₃, we identify the strain generation time (1.31 ps), which occurs simultaneously with Ti-OH* formation, and a tensile strain of 0.06% (upper limit 0.6%). In addition to fully characterizing this intermediate across visible, mid-infrared, and now GHz-THz probes on SrTiO₃, we show that strain fields occur with the creation of some M-OH*, which modifies design strategies for tuning catalytic activity and provides insight into photo-induced degradation so prevalent for OER. To that end, the work put forth here provides a unique methodology to characterize intermediate-induced interfacial strain across OER catalysts.

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.

Ultrafast Dynamics and Liquid Structure in Mesoporous Silica: Propagation of Surface Effects in a Polar Aprotic Solvent

Enhancement of processes ranging from gas sorption to ion conduction in a liquid can be substantial upon nanoconfinement. Here, the dynamics of a polar aprotic solvent, 1-methylimidazole (MeIm), in mesoporous silica (2.8, 5.4, and 8.3 nm pore diameters) were examined using femtosecond infrared vibrational spectroscopy and molecular dynamics simulations of a dilute probe, the selenocyanate (SeCN^-) anion. The long vibrational lifetime and sensitivity of the CN stretch enabled a comprehensive investigation of the relatively slow time scales and subnanometer distance dependences of the confined dynamics. Because MeIm does not readily donate hydrogen bonds, its interactions in the hydrophilic silanol pores differ more from the bulk than those of water confined in the same mesopores, resulting in greater structural order and more dramatic slowing of dynamics. The extent of surface effects was quantified by modified two-state models used to fit three spatially averaged experimental observables: vibrational lifetime, orientational relaxation, and spectral diffusion. The length scales and the models (smoothed step, exponential decay, and simple step) describing the transitions between the distinctive shell behavior at the surface and the bulk-like behavior at the pore interior were compared to those of water. The highly nonuniform distributions of the SeCN^- probe and antiparallel layering of MeIm revealed by the simulations guided the interpretation of the results and development of the analytical models. The results illustrate the importance of electrostatic effects and H-bonding interactions in the behavior of confined liquids.

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

Proteins function as ensembles of interconverting structures. The motions span from picosecond bond rotations to millisecond and longer subunit displacements. Characterization of functional dynamics on all spatial and temporal scales remains challenging experimentally. Two-dimensional infrared spectroscopy (2D IR) is maturing as a powerful approach for investigating proteins and their dynamics. We outline the advantages of IR spectroscopy, describe 2D IR and the information it provides, and introduce vibrational groups for protein analysis. We highlight example studies that illustrate the power and versatility of 2D IR for characterizing protein dynamics and conclude with a brief discussion of the outlook for biomolecular 2D IR.

Vibrational exciton delocalization precludes the use of infrared intensities as proxies for surfactant accumulation on aqueous surfaces

Surface-sensitive vibrational spectroscopy is a common tool for measuring molecular organization and intermolecular interactions at interfaces. Peak intensity ratios are typically used to extract molecular information from one-dimensional spectra but vibrational coupling between surfactant molecules can manifest as signal depletion in one-dimensional spectra. Through a combination of experiment and theory, we demonstrate the emergence of vibrational exciton delocalization in infrared reflection–absorption spectra of soluble and insoluble surfactants at the air/water interface. Vibrational coupling causes a significant decrease in peak intensities corresponding to C–F vibrational modes of perfluorooctanoic acid molecules. Vibrational excitons also form between arachidic acid surfactants within a compressed monolayer, manifesting as signal reduction of C–H stretching modes. Ionic composition of the aqueous phase impacts surfactant intermolecular distance, thereby modulating vibrational coupling strength between surfactants. Our results serve as a cautionary tale against employing alkyl and fluoroalkyl vibrational peak intensities as proxies for concentration, although such analysis is ubiquitous in interface science.

Coupling between surfactant molecules at the air/water interface bleeds intensity into a diffuse background, such that single-wavelength vibrational intensity is effectively depleted at high surface coverage.

Examining the Role of Different Molecular Interactions on Activation Energies and Activation Volumes in Liquid Water

There are a large number of force fields available to model water in molecular dynamics simulations, which each have their own strengths and weaknesses in describing the behavior of the liquid. One particular weakness in many of these models is their description of dynamics away from ambient conditions, where their ability to reproduce measurements is mixed. To investigate this issue, we use the recently developed fluctuation theory for dynamics to directly evaluate measures of the local temperature and pressure dependence: the activation energy and the activation volume. We examine these activation parameters for hydrogen-bond jump exchange times, OH reorientation times, and diffusion coefficients calculated from the SPC/E, SPC/Fw, TIP3P-PME, TIP3P-PME/Fw, OPC3, TIP4P/2005, TIP4P/Ew, E3B2, and E3B3 water models. Activation energy decompositions available through the fluctuation theory approach provide mechanistic insight into the origins of different temperature dependences between the various models, as well as the influence of three-body effects and flexibility.

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

The dynamics of a vibrational frequency in a condensed phase environment, i.e., the spectral diffusion, has attracted considerable interest over the last two decades. A significant impetus has been the development of two-dimensional infrared (2D-IR) photon-echo spectroscopy that represents a direct experimental probe of spectral diffusion, as measured by the frequency-frequency time correlation function (FFCF). In isotopically dilute water, which is perhaps the most thoroughly studied system, the standard interpretation of the longest timescale observed in the FFCF is that it is associated with hydrogen-bond exchange dynamics. Here, we investigate this connection by detailed analysis of both the spectral diffusion timescales and their associated activation energies. The latter are obtained from the recently developed fluctuation theory for the dynamics approach. The results show that the longest timescale of spectral diffusion obtained by the typical analysis used cannot be directly associated with hydrogen-bond exchanges. The hydrogen-bond exchange time does appear in the decay of the water FFCF, but only as an additional, small-amplitude (<3%) timescale. The dominant contribution to the long-time spectral diffusion dynamics is considerably shorter than the hydrogen-bond exchange time and exhibits a significantly smaller activation energy. It thus arises from hydrogen-bond rearrangements, which occur in between successful hydrogen-bond partner exchanges, and particularly from hydrogen bonds that transiently break before returning to the same acceptor.

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.

Ultrafast Dynamics at the Lipid-Water Interface: DMSO Modulates H-Bond Lifetimes

Dimethyl sulfoxide (DMSO) is a common cosolvent and cryopreservation agent used to freeze cells and tissues. DMSO alters the H-bond structure of water, but its interactions with biomolecules and, specifically, with biological interfaces remain poorly understood. Here we investigate the effects of DMSO on the H-bond dynamics at the lipid-water interface using a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations. Ester carbonyl absorption spectra show that DMSO dehydrates the interface, and simulations show that the area per lipid is decreased. Ultrafast 2D IR spectra measure the time scales of frequency fluctuations at the ester carbonyl positions located precisely between the hydrophobic and hydrophilic regions of the membrane. 2D IR measurements show that low DMSO concentrations (<10 mol %) induce ∼40% faster H-bond dynamics compared with pure water, whereas increased concentrations (>10-20 mol %) once again slow down the dynamics. This slow-fast-slow trend is described in terms of two different solvation regimes. Below 10 mol %, DMSO weakens the interfacial H bond, leading to faster "bulk-like" dynamics, whereas above 10 mol %, water molecules become "relatively immobilized" as the H-bond networks becoming disrupted by the H-bond donor/acceptor imbalance at the interface. These studies are an important step toward characterizing the environments around lipid membranes, which are essential to numerous biological processes.

Two-dimensional infrared spectral explorations into bilayer and monolayer self-assemblies of amphiphilic polypeptides

Poly(2-(3-((2-hydroxyethyl)amino)-3-oxopropyl)ethyleneamido) (PHAOE) is an amphiphilic polypeptide. The self-assembly is significant, but the ultrafast dynamic analyses of the peptide self-assembly are exiguous and worth further exploring. In this investigation, the temporal dynamic characteristics of the aggregates and unaggregated PHAOEs are mined by the two-dimensional infrared (2D IR) spectroscopy. The homogeneous and inhomogeneous diffusion processes of the carbonyl stretching modes of the unaggregated PHAOEs are slower than those of the self-assemblies. The inhomogeneous spectral diffusion proportion of the biopolymer PHAOE in methanol is greater than that in dimethyl sulfoxide (DMSO). The solvation shells surround the aggregates and unaggregated PHAOEs in the protic solvent methanol, but there are not any solvation shells around the aggregates or unaggregated PHAOEs in the dipolar solvent DMSO. The massive hydrogen-bonded monolayer self-assembly has merely an aggregate of PHAOEs and no solvation shell in DMSO. But the hydrogen-bonded bilayer self-assembly has a self-assembled methanol shell and an interior aggregate of PHAOEs in methanol. The self-assemblies of PHAOEs motivate the methanols to self-assemble. The large delocalized amide structure results in the fast spectral diffusion of the carbonyl stretching mode.Communicated by Ramaswamy H. Sarma.

Molecular origins of bulk viscosity in liquid water

The shear viscous response of water is closely associated with changes in network connectivity on the sub ps timescale. The bulk viscous response is shown here to be associated with local density fluctuations and rotational motion around 1–3 ps.

Measuring Ultrafast Spectral Diffusion and Correlation Dynamics by Two-Dimensional Electronic Spectroscopy

The frequency fluctuation correlation function (FFCF) measures the spectral diffusion of a state's transition while the frequency fluctuation cross-correlation function (FXCF) measures the correlation dynamics between the transitions of two separate states. These quantities contain a wealth of information on how the chromophores or excitonic states interact and couple with its environment and with each other. We summarize the experimental implementations and theoretical considerations of using two-dimensional electronic spectroscopy to characterize FFCFs and FXCFs. Applications can be found in systems such as the chlorophyll pigment molecules in light-harvesting complexes and CdSe nanomaterials.