Using the thiocyanate anion as a vibrational probe chromophore in conjunction with infrared and NMR spectroscopy, we find that SCN– strongly associates with the cationic head group of dodecyltrimethylammonium bromide (DTAB) micelles, both in normal-phase and reverse micelles. In competition with chloride and iodide ions, we find no evidence for displacement of thiocyanate, in accord with the chaotropicity of the Hofmeister ordering, while lending support to a direct interaction picture of its origin. Ultrafast 2D-IR spectroscopy of the SCN– probe in a range of DTAB micelle sizes (w0 = 4 to w0 = 12) shows little if any size dependence on the time scale for spectral diffusion, which is found to be ∼3.5 times slower than in bulk water (both D2O and H2O). Normal-phase micelles studied with 2D-IR exhibit essentially the same spectral dynamics as do reverse micelles, indicating a lack of sensitivity to interfacial curvature. Combined with 1H NMR chemical shift perturbations, we conclude that the SCN– ions tightly associate with the head groups and are partially buried. The 3–4-fold slowdown in spectral diffusion is consistent with the excluded volume model for interfacial perturbation to hydrogen bond reorientation dynamics. On the basis of these observations and comparisons to previous studies of zwitterionic interfaces probed with phosphate transitions, we conclude that the SCN– spectral dynamics in both reverse- and normal-phase micelles is largely dominated by hydration contributions, and offers a promising probe of interfacial hydration at cationic interfaces. Addition of competitive anions alters neither the IR spectra nor the ultrafast dynamics, indicating that SCN– is robustly associated with the head groups.

Ultrafast picosecond time scale two-dimensional infrared (2D-IR) spectroscopy of a new water-soluble transition metal complex acting as a vibrational probe shows that over a range of concentration and poly(ethylene glycol) (PEG) molecular mass (2000, 8000, and 20000 Da) the time scale of the sensed hydration dynamics differs negligibly from bulk water (D2O). PEG is well-known to establish a highly stable hydration shell because the spacing between adjacent ethereal oxygens nearly matches water’s hydrogen-bonding network. Although these first-shell water molecules are likely significantly retarded, they present an interface to subsequent hydration shells and thus diminish the largely entropic perturbation to water’s orientational dynamics. In addition to the longer PEGs, a series of concentration-dependent 2D-IR measurements using aqueous PEG-400 show a pronounced hydration slowdown in the vicinity of the critical overlap concentration (c*). Comparison between these dynamical results and previously reported steady-state infrared spectroscopy of aqueous PEG-1000 solutions reveals a strikingly identical dependence on number of water molecules per ethylene oxide monomer, scaled according to the critical overlap concentration.

Hydrogenase enzymes enable organisms to use H2 as an energy source, having evolved extremely efficient biological catalysts for the reversible oxidation of molecular hydrogen. Small-molecule mimics of these enzymes provide both simplified models of the catalysis reactions and potential artificial catalysts that might be used to facilitate a hydrogen economy. We have studied two diiron hydrogenase mimics, μ-pdt-[Fe(CO)3]2 and μ-edt-[Fe(CO)3]2 (pdt = propanedithiolate, edt = ethanedithiolate), in a series of alkane solvents and have observed significant ultrafast spectral dynamics using two-dimensional infrared (2D-IR) spectroscopy. Since solvent fluctuations in nonpolar alkanes do not lead to substantial electrostatic modulations in a solute’s vibrational mode frequencies, we attribute the spectral diffusion dynamics to intramolecular flexibility. The intramolecular origin is supported by the absence of any measurable solvent viscosity dependence, indicating that the frequency fluctuations are not coupled to the solvent motional dynamics. Quantum chemical calculations reveal a pronounced coupling between the low-frequency torsional rotation of the carbonyl ligands and the terminal CO stretching vibrations. The flexibility of the CO ligands has been proposed to play a central role in the catalytic reaction mechanism, and our results highlight that the CO ligands are highly flexible on a picosecond time scale.

Rhenium catalysts have shown promise to promote carbon neutrality by reducing a prominent greenhouse gas, CO2, to CO and other starting materials. Much research has focused on identifying intermediates in the photocatalysis mechanism as well as time scales of relevant ultrafast processes. Recent studies have implemented multidimensional spectroscopies to characterize the catalyst’s ultrafast dynamics as it undergoes the many steps of its photocycle.Two-dimensional infrared (2D-IR) spectroscopy is a powerful method to obtain molecular structure information while extracting time scales of dynamical processes with ultrafast resolution. Many observables result from 2D-IR experiments including vibrational lifetimes, intramolecular redistribution time scales, and, unique to 2D-IR, spectral diffusion, which is highly sensitive to solute–solvent interactions and motional dynamics.Spectral diffusion, a measure of how long a vibrational mode takes to sample its frequency space due to multiple solvent configurations, has various contributing factors. Properties of the solvent, the solute’s structural flexibility, and electronic properties, as well as interactions between the solvent and solute, complicate identifying the origin of the spectral diffusion. With carefully chosen experiments, however, the source of the spectral diffusion can be unveiled.Within the context of a considerable body of previous work, here we discuss the spectral diffusion of several rhenium catalysts at multiple stages in the catalysis. These studies were performed in multiple polar liquids to aid in discovering the contributions of the solvent. We also performed electronic ground state 2D-IR and electronic excited state transient-2D-IR experiments to observe how spectral diffusion changes upon electronic excitation. Our results indicate that with the original Lehn catalyst in THF, relative to the ground state, the spectral diffusion slows by a factor of 3 in the equilibrated triplet metal-to-ligand charge transfer state. We attribute this slowdown to a decrease in dielectric friction as well as an increase in molecular flexibility. It is possible to partially simulate the charge transfer by altering the electron density moderately by adding electron donating or withdrawing substituents symmetrically to the bipyridine ligand. We find that unlike the significant electronic structure change induced by MLCT, such small substituent effects do not influence the spectral diffusion. A solvent study in THF, DMSO, and CH3CN found there to be an explicit solvent dependence that we can correlate to the solvent donicity, which is a measure of its nucleophilicity. Future studies focused on the solvent effects on spectral diffusion in the crucial photoinitiated state can illuminate the role the solvent plays in the catalysis.

The ultrafast dynamics of a de novo metalloenzyme active site is monitored using two-dimensional infrared spectroscopy. The homotrimer of parallel, coiled coil α-helices contains a His3-Cu(I) metal site where CO is bound and serves as a vibrational probe of the hydrophobic interior of the self-assembled complex. The ultrafast spectral dynamics of Cu-CO reveals unprecedented ultrafast (2 ps) nonequilibrium structural rearrangements launched by vibrational excitation of CO. This initial rapid phase is followed by much slower ∼40 ps vibrational relaxation typical of metal-CO vibrations in natural proteins. To identify the hidden coupled coordinate, small molecule analogues and the full peptide were studied by QM and QM/MM calculations, respectively. The calculations show that variation of the histidines’ dihedral angles in coordinating Cu controls the coupling between the CO stretch and the Cu–C–O bending coordinates. Analysis of different optimized structures with significantly different electrostatic field magnitudes at the CO ligand site indicates that the origin of the stretch–bend coupling is not directly due to through-space electrostatics. Instead, the large, ∼3.6 D dipole moments of the histidine side chains effectively transduce the electrostatic environment to the local metal coordination orientation. The sensitivity of the first coordination sphere to the protein electrostatics and its role in altering the potential energy surface of the bound ligands suggests that long-range electrostatics can be leveraged to fine-tune function through enzyme design.

Although it is well known that water is essential for biological function, it has been a challenge to determine how water behaves near biomacromolecular interfaces, and what role water plays in influencing the dynamics of the biochemical machinery. By adopting a vibrational labeling strategy coupled with ultrafast two-dimensional infrared (2D-IR) spectroscopy, it has recently become possible to study hydration dynamics, site specifically at the surface of proteins and model membranes. We review our recent progress in measuring hydration dynamics in contexts ranging from small-molecule solutes to biomacromolecules in dilute, viscous, and crowded environments.The combination of site-specific chemical labeling, ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulation has enabled detailed investigations of biomolecule hydration dynamics.

The spectral dynamics of a series of rhenium photocatalysts, fac-Re(4,4′-R2-bpy)(CO)3Cl, where R = H, methyl, t-butyl, and carboxylic acid, as well as Re(1,10-phenanthroline)(CO)3Cl were observed in multiple aprotic solvents using two-dimensional infrared spectroscopy (2DIR). The carbonyl vibrational stretching frequencies showed slight variations due to the electron-donating or -withdrawing nature of the substituents on the bipyridine. The different substituents had minimal to no influence on the spectral diffusion time scales of the compounds within a particular solvent, but among the three different solvents investigated (DMSO, THF, and CH3CN), we find the spectral diffusion times to correlate with the solvent’s donor number (DN). Because the donicity is a measure the Lewis basicity of the solvent, these findings may help establish a more complete dynamical picture of the photocatalysis, where the first chemical step following optical excitation is electron transfer from a sacrificial donor to the rhenium complex.

Using a derivative of the vitamin biotin labeled with a transition-metal carbonyl vibrational probe in a series of aqueous N,N-dimethylformamide (DMF) solutions, we observe a striking slowdown in spectral diffusion dynamics with decreased DMF concentration. Equilibrium solvation dynamics, measured with the rapidly acquired spectral diffusion (RASD) technique, a variant of heterodyne-detected photon–echo peak shift experiments, range from 1 ps in neat DMF to ∼3 ps in 0.07 mole fraction DMF/water solution. Molecular dynamics simulations of the biotin–metal carbonyl solute in explicit aqueous DMF solutions show marked preferential solvation by DMF, which becomes more pronounced at lower DMF concentrations. The simulations and the experimental data are consistent with an interpretation where the slowdown in spectral diffusion is due to solvent exchange involving distinct cosolvent species. A simple two-component model reproduces the observed spectral dynamics as well as the DMF concentration dependence, enabling the extraction of the solvent exchange time scale of 8 ps. This time scale corresponds to the diffusive motion of a few Å, consistent with a solvent-exchange mechanism. Unlike most previous studies of solvation dynamics in binary mixtures of polar solvents, our work highlights the ability of vibrational probes to sense solvent exchange as a new, slow component in the spectral diffusion dynamics.

Cytoplasmic osmolytes can significantly alter the thermodynamic and kinetic properties of proteins relative to those under dilute solution conditions. Spectroscopic experiments of lysozymes in cosolvents indicate that such changes may arise from the heterogeneous, site-specific hydrophobic interactions between protein surface residues and individual solvent molecules. In pursuit of an accurate and predictive model for explaining biomolecular interactions, we study the averaged structural characteristics of mixed solvents with homologous lysozyme solutes using all-atom molecular dynamics. By observing the time-averaged densities of different aqueous solutions of trifluoroethanol, we deduce trends in the heterogeneous solvent interactions over each protein’s surface, and investigate how the homology of protein structure does not necessarily translate to similarities in solvent structure and composition—even when observing identical side chains.

A detailed understanding of photocatalyzed reaction dynamics requires a sensitive means of investigating the transient catalytically active species. Ideally, the method should be able to compare the electronically excited photocatalyst directly to the ground state species. We use equilibrium and transient two-dimensional infrared (2DIR and t-2DIR) spectroscopy to study the ground and excited state spectral dynamics of [Re(CO)3(bpy)Cl] in tetrahydrofuran (THF). We leverage the long-lived triplet excited state of the molecule to re-establish an equilibrated state relative to intersystem crossing dynamics and external solvent fluctuations, allowing access to the dynamics experienced by the excited state photocatalyst. The decay of frequency correlations within the excited triplet state species differs significantly from the ground state (slower by a factor of 3), indicating that the electronic excitation and subsequent metal-to-ligand charge transfer and associated structural changes are sufficient to perturb the spectral dynamics as sensed by the carbonyl ligands. In addition, we observe a 2-fold slowdown in ground state spectral dynamics around the in-phase symmetric vibrational mode compared to the two lower frequency, out-of-phase symmetric and asymmetric modes. Following electronic absorption and metal-to-ligand charge transfer the symmetry of the vibrational modes are disrupted, and all vibrational modes experience inhomogeneous broadening and spectral diffusion. The qualitative change in broadening mechanisms arises from the charge redistribution, indicating that direct comparisons of vibrational spectral dynamics on different electronic states—reported here for the first time—can be highly sensitive indicators of changes in electronic structure and in the concomitant solvation dynamics that underlie the microscopic details of charge transfer reactions.

Ultrafast two-dimensional infrared (2D-IR) spectroscopy reveals picosecond protein and hydration dynamics of crowded hen egg white lysozyme (HEWL) labeled with a metal–carbonyl vibrational probe covalently attached to a solvent accessible His residue. HEWL is systematically crowded alternatively with polyethylene glycol (PEG) or excess lysozyme in order to distinguish the chemically inert polymer from the complex electrostatic profile of the protein crowder. The results are threefold: (1) A sharp dynamical jamming-like transition is observed in the picosecond protein and hydration dynamics that is attributed to an independent-to-collective hydration transition induced by macromolecular crowding that slows the hydration dynamics up to an order of magnitude relative to bulk water. (2) The interprotein distance at which the transition occurs suggests collective hydration of proteins over distances of 30–40 Å. (3) Comparing the crowding effects of PEG400 to our previously reported experiments using glycerol exposes fundamental differences between small and macromolecular crowding agents.

Using an implementation of heterodyne-detected vibrational echo spectroscopy, we show that equilibrium spectral diffusion caused by solvation dynamics can be measured in a fraction of the time required using traditional two-dimensional infrared spectroscopy. Spectrally resolved, heterodyne-detected rephasing and nonrephasing signals, recorded at a single delay between the first two pulses in a photon echo sequence, can be used to measure the full waiting time dependent spectral dynamics that are typically extracted from a series of 2D-IR spectra. Hence, data acquisition is accelerated by more than 1 order of magnitude, while permitting extremely fine sampling of the spectral dynamics during the waiting time between the second and third pulses. Using cymantrene (cyclopentadienyl manganese tricarbonyl, CpMn(CO)3) in alcohol solutions, we compare this novel approach—denoted rapidly acquired spectral diffusion (RASD)—with a traditional method using full 2D-IR spectra, finding excellent agreement. Though this approach is largely limited to isolated vibrational bands, we also show how to remove interference from cross-peaks that can produce characteristic modulations of the spectral dynamics through vibrational quantum beats.

One route to accessing site-specific dynamical information available with ultrafast multidimensional infrared spectroscopy is the development of robust and versatile vibrational probes. Here we synthesize and characterize a vibrationally labeled cholesterol derivative, (cholesteryl benzoate) chromium tricarbonyl, to probe model lipid membranes, focusing specifically on the membrane–water interface. Utilizing FTIR and polarized-ATR spectroscopies, we determine the location of the chromium tricarbonyl motif to be situated at the water–membrane interface with an orientation of 46 ± 2° relative to the vector normal to the membrane surface. We test the dynamical sensitivity of the (cholesteryl benzoate) chromium tricarbonyl label with two different nonlinear infrared spectroscopy methods, both of which show that the probe is well-suited to the study of membrane dynamics as well as the dynamics of water at the membrane interface. The metal carbonyl vibrational probe located at the surface of a bicelle exhibits spectral diffusion dynamics induced by membrane hydration water that is roughly three times slower than observed using a nearly identical vibrational probe in bulk water.

A fundamental aspect of Fourier transform (FT) spectroscopy is the inverse relationship between frequency resolution and the maximum scanned time delay. In situations where essential chemical information is contained in spectral peak amplitudes rather than in their detailed shapes, it is possible to dramatically reduce the experimental acquisition time of time domain methods such as two-dimensional infrared (2D-IR) spectroscopy. By introducing compressed sensing to the analysis and experimental design of 2D-IR spectroscopy, we show that waiting-time-dependent 2D peak amplitudes reproduce conventional FT acquisition and analysis but can be recorded in a fraction of the time. Peak amplitude data are often sufficient for measuring intramolecular vibrational redistribution, vibrational coherence, chemical exchange, population, and orientational relaxation, as well as spectral diffusion.Keywords: compressed sampling; Fourier transform; matching pursuit; signal processing; spectral analysis; spectroscopy;

There is considerable evidence for the slaving of biomolecular dynamics to the motions of the surrounding solvent environment, but to date there have been few direct experimental measurements capable of site-selectively probing both the dynamics of the water and the protein with ultrafast time resolution. Here, two-dimensional infrared spectroscopy (2D-IR) is used to study the ultrafast hydration and protein dynamics sensed by a metal carbonyl vibrational probe covalently attached to the surface of hen egg white lysozyme dissolved in D2O/glycerol solutions. Surface labeling provides direct access to the dynamics at the protein–water interface, where both the hydration and the protein dynamics can be observed simultaneously through the vibrational probe’s frequency–frequency correlation function. In pure D2O, the correlation function shows a fast initial 3 ps decay corresponding to fluctuations of the hydration water, followed by a significant static offset attributed to fluctuations of the protein that are not sampled within the <20 ps experimental window. Adding glycerol increases the bulk solvent viscosity while leaving the protein structurally intact and hydrated. The hydration dynamics exhibit a greater than 3-fold slowdown between 0 and 80% glycerol (v/v), and the contribution from the protein’s dynamics is found to slow in a nearly identical fashion. In addition, the magnitude of the dynamic slowdown associated with hydrophobic hydration is directly measured and shows quantitative agreement with predictions from molecular dynamics simulations.

Water is capable of assisting exceptionally rapid vibrational relaxation within dissolved solute species. Although ultrafast dynamics of metal carbonyl complexes have long served as models for vibrational relaxation, all reports to-date have investigated nonaqueous solutions due to the insolubility of the vast majority of metal carbonyl complexes in water. Using the water-soluble complex [RuCl2(CO)3]2, which belongs to a class known as “carbon monoxide (CO) releasing molecules” (CORM), we report the first ultrafast vibrational relaxation measurements of a metal carbonyl complex in water and compare this relaxation with polar organic solvents, namely, methanol. The vibrational relaxation, measured by two-dimensional IR (2D-IR) spectroscopy, is an order of magnitude faster in H2O (3.12 ± 0.29 ps) than in methanol (42.25 ± 3 ps). The accelerated relaxation times of the coupled CO units in H2O and D2O is interpreted as resulting from the enhancement of intramolecular relaxation pathways through additional coupling induced by the solvent. In addition, the vibrational lifetime shows a significant isotope dependence: in D2O the relaxation time is 4.27 ± 0.27 ps, a difference of roughly 30%. We interpret these measurements in terms of a nonresonant channel primarily arising from water’s reorientational dynamics, which occur primarily through large angular jumps, as well as a resonant transfer of vibrational energy from the carbonyl bands to the libration-bend combination band. These measurements indicate that metal carbonyls, which are among the strongest IR transitions, are exquisitely sensitive to the presence of water and hold promise as IR analogs of EPR spin labels.

The thermodynamic driving forces for protein folding, association, and function are often determined by protein–water interactions. With a novel covalently bound labeling approach, we have used sensitive vibrational probes, site-selectively conjugated to two lysozyme variants—in conjunction with ultrafast two-dimensional infrared (2D-IR) spectroscopy—to investigate directly the protein–water interface. By probing alternatively a topologically flat, rigid domain and a flexible domain, we find direct experimental evidence for spatially heterogeneous hydration dynamics. The hydration environment around globular proteins can vary from exhibiting bulk-like hydration dynamics to dynamically constrained water, which results from stifled hydrogen bond switching dynamics near extended hydrophobic surfaces. Furthermore, we leverage preferential solvation exchange to demonstrate that the liberation of dynamically constrained water is a sufficient driving force for protein–surface association reactions. These results provide an intuitive picture of the dynamic aspects of hydrophobic hydration of proteins, illustrating an essential function of water in biological processes.

Two-dimensional infrared (2DIR) spectroscopy is applied to both (Cp)2Fe2(CO)4 and its ruthenium analog (Cp)2Ru2(CO)4 in order to study the vibrational dynamics of these two systems. Combining the results of 2DIR spectroscopy and DFT calculations, the different structural forms of both the iron and the ruthenium complexes were characterized, furthering the previous assignment of the linear IR spectrum by determining the transition frequencies associated with the different isomeric forms. Monitoring the time-dependent amplitudes of the cross peaks enabled the observation of equilibrium energy transfer dynamics between different vibrational modes of the cis-B (Cp)2Fe2(CO)4 and the gauche-NB (Cp)2Ru2(CO)4 complexes. Treating the energy transfer as an equilibrium process, we extracted the rate constants associated with both the uphill and the downhill transfer of vibrational energy, finding that the difference in the rate constants of the two metal complexes maps to the difference in the energy gap between the two modes involved.

Ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations of Mn2(CO)10 in a series of linear alcohols reveal that the rate of intramolecular vibrational redistribution among the terminal carbonyl stretches is dictated by the average number of hydrogen bonds formed between the solute and solvent. The presence of hydrogen bonds was found to hinder vibrational redistribution between eigenstates, while leaving the overall T1 relaxation rate unchanged.

We present the first implementation of transient vibrational Stark-effect spectroscopy as an ultrafast probe of solvation dynamics. The method is applied to the phototriggered intramolecular charge-transfer reaction of Betaine-30, where the vibrational Stark shifts of the nearby solvent molecules—arising from the change in the electrostatic environment—are measured using a three-pulse photon echo probe. This new experiment provides a direct subpicosecond measure of the chromophore’s excited-state dynamics and back electron transfer as viewed from the solvent’s perspective. We develop a simple ab initio model that offers semiquantitative prediction of the experimental Stark shifts.

Two-dimensional infrared spectroscopy (2DIR) is used to measure the viscosity-dependent spectral diffusion of a model vibrational probe, Mn2(CO)10 (dimanganese decacarbonyl, DMDC), in a series of alcohols with time scales ranging from 2.67 ps in methanol to 5.33 ps in 1-hexanol. Alcohol−alkane solvent mixtures were found to produce indistinguishable linear IR spectra, while still demonstrating viscosity-dependent spectral diffusion. Using a vibrational exciton model to characterize the inhomogeneous energy landscape, several analogies emerge with multichromophoric electronic systems, such as J-aggregates and light-harvesting protein complexes. An excitonic, local vibrational mode Hamiltonian parametrized to reproduce the vibrational structure of DMDC serves as a starting point from which site energies (i.e., local carbonyl frequencies) are given Gaussian distributed disorder. The model gives excellent agreement with both the linear IR spectrum and the inhomogeneous widths extracted from 2DIR, indicating the system can be considered to be a “vibrational aggregate.” This model naturally leads to exchange narrowing due to disorder-induced exciton localization, producing line widths consistent with our 1D and 2D measurements. Further, the diagonal disorder alone effectively reduces the molecular symmetry, leading to the appearance of Raman bands in the IR spectrum in accord with the measurements. Here, we show that the static inhomogeneity of the excitonic model with disorder successfully captures the essential details of the 1D spectrum while predicting the degree of IR activity of forbidden modes as well as the inhomogeneous widths and relative magnitudes of the transition moments.

Cyclopentadienylmolybdenum(II) tricarbonyl dimer exists in two different equilibrium conformations: trans and gauche. Ultrafast photoexcitation in the ultraviolet cleaves the Mo−Mo bond, permitting observation of the subsequent geminate rebinding reaction (trebinding = 31.6 ps) by monitoring infrared bleach recoveries at frequencies corresponding to the CO stretches of the trans and gauche isomers, the time-resolved measurements revealed that the monomers rebind in the trans configuration only. Further insight into the rebinding reaction was obtained by mapping the full potential energy surface along the reaction coordinate using electronic-structure methods.

Transient two-dimensional infrared (2DIR) spectroscopy is applied to the photodissociation of Mn2(CO)10 to 2 Mn(CO)5 in cyclohexane solution. By varying both the time delay between the 400 nm phototrigger and the 2DIR probe as well as the waiting time in the 2DIR pulse sequence, we directly determine the orientational relaxation of the vibrationally hot photoproduct. The orientational relaxation slows as the photoproduct cools, providing a measure of the transient temperature decay time of 70 ± 16 ps. We compare the experimental results with molecular dynamics simulations and find near quantitative agreement for equilibrium orientational diffusion time constants but only qualitative agreement for nonequilibrium thermal relaxation. The simulation also shows that the experiment probes an unusual regime of thermal excitation, where the solute is heated while the solvent remains essentially at room temperature.

Ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations of Mn2(CO)10 in a series of linear alcohols reveal that the rate of intramolecular vibrational redistribution among the terminal carbonyl stretches is dictated by the average number of hydrogen bonds formed between the solute and solvent. The presence of hydrogen bonds was found to hinder vibrational redistribution between eigenstates, while leaving the overall T1 relaxation rate unchanged.