Reply to "Comment on: 'Isolating Vibrational Polariton 2D-IR Transmission Spectra'"

In a Comment on our recent Letter, the authors take issue with our method of refining 2D-IR transmission spectra to remove a background contribution that arises from nonpolaritonic molecules in the cavity. In our response to their Comment, we describe how our approach was motivated by the previous work of the authors, and we present a spatially dependent molecule-cavity Tavis-Cummings model that can account for the significant response from localized molecules with nonzero oscillator strengths. The telltale signature of the localized molecule response is the spectral diffusion dynamics of the bare W(CO)₆ molecules in the polar butyl acetate solvent. Inhomogeneous broadening is absent from polaritonic states due to the extreme degree of exchange narrowing in coupling very large numbers of molecules to a cavity mode.

Isolating Polaritonic 2D-IR Transmission Spectra

Strong coupling between vibrational transitions in molecules within a resonant optical microcavity leads to the formation of collective, delocalized vibrational polaritons. There are many potential applications of "polaritonic chemistry", ranging from modified chemical reactivity to quantum information processing. One challenge in obtaining the polaritonic response is removing a background contribution due to the uncoupled molecules that generate an ordinary 2D-IR spectrum whose amplitude is filtered by the polariton transmission spectrum. We show that most features in 2D-IR spectra of vibrational polaritons can be explained by a linear superposition of this background signal and the true polariton response. Through a straightforward correction procedure, in which the filtered bare-molecule 2D-IR spectrum is subtracted from the measured cavity response, we recover the polaritonic spectrum.

Ultrafast vibrational dynamics of a solute correlates with dynamics of the solvent

Two-dimensional infrared (2D-IR) spectroscopy is used to measure the spectral dynamics of the metal carbonyl complex cyclopentadienyl manganese tricarbonyl (CMT) in a series of linear alkyl nitriles. 2D-IR spectroscopy provides direct readout of solvation dynamics through spectral diffusion, probing the decay of frequency correlation induced by fluctuations of the solvent environment. 2D-IR simultaneously monitors intramolecular vibrational energy redistribution (IVR) among excited vibrations, which can also be influenced by the solvent through the spectral density rather than the dynamical friction underlying solvation. Here, we report that the CMT vibrational probe reveals solvent dependences in both the spectral diffusion and the IVR time scales, where each slows with increased alkyl chain length. In order to assess the degree to which solute-solvent interactions can be correlated with bulk solvent properties, we compared our results with low-frequency dynamics obtained from optical Kerr effect (OKE) spectroscopy-performed by others-on the same nitrile solvent series. We find excellent correlation between our spectral diffusion results and the orientational dynamics time scales from OKE. We also find a correlation between our IVR time scales and the amplitudes of the low-frequency spectral densities evaluated at the 90-cm^-1 energy difference, corresponding to the gap between the two strong vibrational modes of the carbonyl probe. 2D-IR and OKE provide complementary perspectives on condensed phase dynamics, and these findings provide experimental evidence that at least at the level of dynamical correlations, some aspects of a solute vibrational dynamics can be inferred from properties of the solvent.

Direct comparison of amplitude and geometric measures of spectral inhomogeneity using phase-cycled 2D-IR spectroscopy

Two-dimensional infrared (2D-IR) spectroscopy provides access to equilibrium dynamics with the extraction of the frequency-fluctuation correlation function (FFCF) from the measured spectra. Several different methods of obtaining the FFCF from experimental spectra, such as the center line slope (CLS), ellipticity, phase slope, and nodal line slope, all depend on the geometrical nature of the 2D line shape and necessarily require spectral extent in order to achieve a measure of the FFCF. Amplitude measures, on the other hand, such as the inhomogeneity index, rely only on signal amplitudes and can, in principle, be computed using just a single point in a 2D spectrum. With a pulse shaper-based 2D-IR spectrometer, in conjunction with phase cycling, we separate the rephasing and nonrephasing signals used to determine the inhomogeneity index. The same measured data provide the absorptive spectrum, needed for the CLS. Both methods are applied to two model molecular systems: tungsten hexacarbonyl (WCO₆) and methylcyclopentadienyl manganese tricarbonyl [Cp'Mn(CO)₃, MCMT]. The three degenerate IR modes of W(CO)₆ lack coherent modulation or noticeable intramolecular vibrational redistribution (IVR) and are used to establish a baseline comparison. The two bands of the MCMT tripod complex include intraband coherences and IVR as well as likely internal torsional motion on a few-picosecond time scale. We find essentially identical spectral diffusion, but faster, non-equilibrium dynamics lead to differences in the FFCFs extracted with the two methods. The inhomogeneity index offers an advantage in cases where spectra are complex and energy transfer can mimic line shape changes due to frequency fluctuations.

Transmission Mode 2D-IR Spectroelectrochemistry of In Situ Electrocatalytic Intermediates

Unraveling electrocatalytic mechanisms, as well as fundamental structural dynamics of intermediates, requires spectroscopy with high time and frequency resolution that can account for nonequilibrium in situ concentration changes inherent to electrochemistry. Two-dimensional infrared (2D-IR) spectroscopy is an ideal candidate, but several technical challenges have hindered development of this powerful tool for spectroelectrochemistry (SEC). We demonstrate a transmission-mode, optically transparent thin-layer electrochemical (OTTLE) cell adapted to 2D-IR-SEC to monitor the important Re(bpy)(CO)₃Cl CO₂-reduction electrocatalyst. 2D-IR-SEC reveals pronounced differences in both spectral diffusion time scales and spectral inhomogeneity in the singly reduced catalyst, [Re(bpy)(CO)₃Cl]^•-, relative to the starting Re(bpy)(CO)₃Cl. Cross-peaks between well-resolved symmetric vibrations and congested low-frequency bands enable direct assignment of all distinct species during the electrochemical reaction. With this information, 2D-IR-SEC provides new mechanistic insights regarding unproductive, catalyst-degrading dimerization. 2D-IR-SEC opens new experimental windows into the electrocatalysis foundation of future energy conversion and greenhouse gas reduction.

Mechanistic Study of Charge Separation in a Nonfullerene Organic Donor-Acceptor Blend Using Multispectral Multidimensional Spectroscopy

Organic photovoltaics (OPVs) based on nonfullerene acceptors are now approaching commercially viable efficiencies. One key to their success is efficient charge separation with low potential loss at the donor-acceptor heterojunction. Due to the lack of spectroscopic probes, open questions remain about the mechanisms of charge separation. Here, we study charge separation of a model system composed of the donor, poly[(2,6-(4,8-bis(5-(2-ethylhexyl)thiophen-2-yl)-benzo[1,2-b:4,5-b']dithiophene))-alt-(5,5-(1',3'-di-2-thienyl-5',7'-bis(2-ethylhexyl)benzo[1',2'-c:4',5'-c']dithiophene-4,8-dione) (PBDB-T), and the nonfullerene acceptor, 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-indanone))-5,5,11,11-tetrakis(4-hexylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (ITIC), using multidimensional spectroscopy spanning the visible to the mid-infrared. We find that bound polaron pairs (BPPs) generated within ITIC domains play a dominant role in efficient hole transfer, transitioning to delocalized polarons within 100 fs. The weak electron-hole binding within the BPPs and the resulting polaron delocalization are key factors for efficient charge separation at nearly zero driving force. Our work provides useful insight into how to further improve the power conversion efficiency in OPVs.

The Photoactive Excited State of the B12-Based Photoreceptor CarH

We have used transient absorption spectroscopy in the UV-visible and X-ray regions to characterize the excited state of CarH, a protein photoreceptor that uses a form of B₁₂, adenosylcobalamin (AdoCbl), to sense light. With visible excitation, a nanosecond-lifetime photoactive excited state is formed with unit quantum yield. The time-resolved X-ray absorption near edge structure difference spectrum of this state demonstrates that the excited state of AdoCbl in CarH undergoes only modest structural expansion around the central cobalt, a behavior similar to that observed for methylcobalamin rather than for AdoCbl free in solution. We propose a new mechanism for CarH photoreactivity involving formation of a triplet excited state. This allows the sensor to operate with high quantum efficiency and without formation of potentially dangerous side products. By stabilizing the excited electronic state, CarH controls reactivity of AdoCbl and enables slow reactions that yield nonreactive products and bypass bond homolysis and reactive radical species formation.

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 XANES Monitors Femtosecond Sequential Structural Evolution in Photoexcited Coenzyme B12

Polarized X-ray absorption near-edge structure (XANES) at the Co K-edge and broadband UV-vis transient absorption are used to monitor the sequential evolution of the excited-state structure of coenzyme B₁₂ (adenosylcobalamin) over the first picosecond following excitation. The initial state is characterized by sub-100 fs sequential changes around the central cobalt. These are polarized first in the y-direction orthogonal to the transition dipole and 50 fs later in the x-direction along the transition dipole. Expansion of the axial bonds follows on a ca. 200 fs time scale as the molecule moves out of the Franck-Condon active region of the potential energy surface. On the same 200 fs time scale there are electronic changes that result in the loss of stimulated emission and the appearance of a strong absorption at 340 nm. These measurements provide a cobalt-centered movie of the excited molecule as it evolves to the local excited-state minimum.

Antivitamins B12 in a Microdrop: The Excited-State Structure of a Precious Sample Using Transient Polarized X-ray Absorption Near-Edge Structure

Polarized transient X-ray absorption near-edge structure (XANES) was used to probe the excited-state structure of a photostable B₁₂ antivitamin (Coβ-2-(2,4-difluorophenyl)-ethynylcobalamin, F₂PhEtyCbl). A drop-on-demand delivery system synchronized to the LCLS X-ray free electron laser pulses was implemented and used to measure the XANES difference spectrum 12 ps following excitation, exposing only ∼45 μL of sample. Unlike cyanocobalamin (CNCbl), where the Co-C bond expands 15-20%, the excited state of F₂PhEtyCbl is characterized by little change in the Co-C bond, suggesting that the acetylide linkage raises the barrier for expansion of the Co-C bond. In contrast, the lower axial Co-N_DMB bond is elongated in the excited state of F₂PhEtyCbl by ca. 10% or more, comparable to the 10% elongation observed for Co-N_DMB in CNCbl.

Probing the Excited State of Methylcobalamin Using Polarized Time-Resolved X-ray Absorption Spectroscopy

We use picosecond time-resolved polarized X-ray absorption near-edge structure (XANES) measurements to probe the structure of the long-lived photoexcited state of methylcobalamin (MeCbl) and the cob(II)alamin photoproduct formed following photoexcitation of adenosylcobalamin (AdoCbl, coenzyme B₁₂). For MeCbl, we used 520 nm excitation and a time delay of 100 ps to avoid the formation of cob(II)alamin. We find only small spectral changes in the equatorial and axial directions, which we interpret as arising from small (<∼0.05 Å) changes in both the equatorial and axial distances. This confirms expectations based on prior UV-visible transient absorption measurements and theoretical simulations. We do not find evidence for the significant elongation of the Co-C bond reported by Subramanian [ J. Phys. Chem. Lett. 2018 , 9 , 1542 - 1546 ] following 400 nm excitation. For AdoCbl, we resolve the difference XANES contributions along three unique molecular axes by exciting with both 540 and 365 nm light, demonstrating that the spectral changes are predominantly polarized along the axial direction, consistent with the loss of axial ligation. These data suggest that the microsecond "recombination product" identified by Subramanian et al. is actually the cob(II)alamin photoproduct that is produced following bond homolysis of MeCbl with 400 nm excitation. Our results highlight the pronounced advantage of using polarization-selective transient X-ray absorption for isolating structural dynamics in systems undergoing atomic displacements that are strongly correlated to the exciting optical polarization.

Multispectral multidimensional spectrometer spanning the ultraviolet to the mid-infrared

Multidimensional spectroscopy is the optical analog to nuclear magnetic resonance, probing dynamical processes with ultrafast time resolution. At optical frequencies, the technical challenges of multidimensional spectroscopy have hindered its progress until recently, where advances in laser sources and pulse-shaping have removed many obstacles to its implementation. Multidimensional spectroscopy in the visible and infrared (IR) regimes has already enabled respective advances in our understanding of photosynthesis and the structural rearrangements of liquid water. A frontier of ultrafast spectroscopy is to extend and combine multidimensional techniques and frequency ranges, which have been largely restricted to operating in the distinct visible or IR regimes. By employing two independent amplifiers seeded by a single oscillator, it is straightforward to span a wide range of time scales (femtoseconds to seconds), all of which are often relevant to the most important energy conversion and catalysis problems in chemistry, physics, and materials science. Complex condensed phase systems have optical transitions spanning the ultraviolet (UV) to the IR and exhibit dynamics relevant to function on time scales of femtoseconds to seconds and beyond. We describe the development of the Multispectral Multidimensional Nonlinear Spectrometer (MMDS) to enable studies of dynamical processes in atomic, molecular, and material systems spanning femtoseconds to seconds, from the UV to the IR regimes. The MMDS employs pulse-shaping methods to provide an easy-to-use instrument with an unprecedented spectral range that enables unique combination spectroscopies. We demonstrate the multispectral capabilities of the MMDS on several model systems.

Solvent Quality Controls Macromolecular Structural Dynamics of a Dendrimeric Hydrogenase Model

We report a spectroscopic investigation of the ultrafast dynamics of the second-generation poly(aryl ether) dendritic hydrogenase model using two-dimensional infrared (2D-IR) spectroscopy to probe the metal carbonyl vibrations of the dendrimer and a reference small molecule, [Fe(μ-S)(CO)₃]₂. We find that the structural dynamics of the dendrimer are reflected in a slow phase of the spectral diffusion, which is absent from [Fe(μ-S)(CO)₃]_2, and we relate the slow phase to the quality of the solvent for poly(aryl ether) dendrimers. We observe a solvent-dependent modulation of the initial phase of vibrational relaxation of the carbonyl groups, which we attribute to an inhibition of solvent assistance in the intramolecular vibrational redistribution process for the dendrimer. There is also a clear solvent dependence of the vibrational frequencies of both the dendrimer and [Fe(μ-S)(CO)₃]₂. Our data represent the first 2D-IR study of a dendritic complex and provide insight into the solvent dependence of molecular conformation in solution and the ultrafast dynamics of moderately sized, conformationally mobile compounds.

Ultrafast X-ray Absorption Near Edge Structure Reveals Ballistic Excited State Structural Dynamics

Polarized ultrafast time-resolved X-ray absorption near edge structure (XANES) allows characterization of excited state dynamics following excitation. Excitation of vitamin B₁₂, cyanocobalamin (CNCbl), in the αβ-band at 550 nm and the γ-band at 365 nm was used to uniquely resolve axial and equatorial contributions to the excited state dynamics. The structural evolution of the excited molecule is best described by a coherent ballistic trajectory on the excited state potential energy surface. Prompt expansion of the Co cavity by ca. 0.03 Å is followed by significant elongation of the axial bonds (>0.25 Å) over the first 190 fs. Subsequent contraction of the Co cavity in both axial and equatorial directions results in the relaxed S₁ excited state structure within 500 fs of excitation.

Solvent exchange in preformed photocatalyst-donor precursor complexes determines efficiency

In homogeneous photocatalytic reduction of CO2, it is widely assumed that the primary electron transfer from the sacrificial donor to the catalyst is diffusion controlled, thus little attention has been paid to optimizing this step. We present spectroscopic evidence that the precursor complex is preformed, driven by preferential solvation, and two-dimensional infrared spectroscopy reveals triethanolamine (donor)/tetrahydrofuran (solvent) exchange in the photocatalyst's solvation shell, reaching greatest magnitude at the known optimal concentration (∼20% v/v TEOA in THF) for catalytically reducing CO2 to CO. Transient infrared absorption shows the appearance of the singly reduced catalyst on an ultrafast (<70 ps) time scale, consistent with non-diffusion controlled electron transfer within the preformed precursor complex. Identification of preferential catalyst-cosolvent interactions suggests a revised paradigm for the primary electron transfer, while illuminating the pivotal importance of solvent exchange in determining the overall efficiency of the photocycle.

Dynamic Flexibility of Hydrogenase Active Site Models Studied with 2D-IR Spectroscopy

Hydrogenase enzymes enable organisms to use H₂ 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)₃]₂ and μ-edt-[Fe(CO)₃]₂ (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.

NOESY-Like 2D-IR Spectroscopy Reveals Non-Gaussian Dynamics

We have identified an unexpected signature of non-Gaussian dynamics in a conventional 2D IR measurement on a system with rapid intermolecular vibrational energy transfer. In a ternary mixture of the CO₂ reduction photocatalyst, ReCl(bpy)(CO)₃, NaSCN, and THF solvent, preferential association between the metal carbonyl catalyst and the NaSCN ion pairs facilitates intermolecular energy transfer on a few picoseconds time scale. Monitoring the cross peak between the highest frequency metal carbonyl band and the CN bands of NaSCN contact ion pairs, we find a striking time evolution of the cross-peak position on the detection axis. This frequency shift, which is due to spectral diffusion following intermolecular energy transfer, occurs with a time scale that is distinct from either the donor or acceptor spectral diffusion measured simultaneously. We argue that the energy transfer, a second-order Förster process, effectively increases the dimensionality of the 2D-IR spectroscopy and thus enables sensitivity to non-Gaussian dynamics.

Histidine Orientation Modulates the Structure and Dynamics of a de Novo Metalloenzyme Active Site

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.

Monitoring equilibrium reaction dynamics of a nearly barrierless molecular rotor using ultrafast vibrational echoes

Using rapidly acquired spectral diffusion, a recently developed variation of heterodyne detected infrared photon echo spectroscopy, we observe ∼3 ps solvent independent spectral diffusion of benzene chromium tricarbonyl (C6H6Cr(CO)3, BCT) in a series of nonpolar linear alkane solvents. The spectral dynamics is attributed to low-barrier internal torsional motion. This tripod complex has two stable minima corresponding to staggered and eclipsed conformations, which differ in energy by roughly half of kBT. The solvent independence is due to the relative size of the rotor compared with the solvent molecules, which create a solvent cage in which torsional motion occurs largely free from solvent damping. Since the one-dimensional transition state is computed to be only 0.03 kBT above the higher energy eclipsed conformation, this model system offers an unusual, nearly barrierless reaction, which nevertheless is characterized by torsional coordinate dependent vibrational frequencies. Hence, by studying the spectral diffusion of the tripod carbonyls, it is possible to gain insight into the fundamental dynamics of internal rotational motion, and we find some evidence for the importance of non-diffusive ballistic motion even in the room-temperature liquid environment. Using several different approaches to describe equilibrium kinetics, as well as the influence of reactive dynamics on spectroscopic observables, we provide evidence that the low-barrier torsional motion of BCT provides an excellent test case for detailed studies of the links between chemical exchange and linear and nonlinear vibrational spectroscopy.

Dynamics of rhenium photocatalysts revealed through ultrafast multidimensional spectroscopy

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.

Heterogeneous preferential solvation of water and trifluoroethanol in homologous lysozymes

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.

Crowding induced collective hydration of biological macromolecules over extended distances

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.

Site-specific measurements of lipid membrane interfacial water dynamics with multidimensional infrared spectroscopy

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.

Site-specific coupling of hydration water and protein flexibility studied in solution with ultrafast 2D-IR 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 D(2)O/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 D(2)O, 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.

Site-specific hydration dynamics of globular proteins and the role of constrained water in solvent exchange with amphiphilic cosolvents

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.

Ultrafast α-like relaxation of a fragile glass-forming liquid measured using two-dimensional infrared spectroscopy

Ultrafast two-dimensional infrared (2D-IR) spectroscopy is used to study the picosecond dynamics of a vibrational probe molecule dissolved in a fragile glass former. The spectral dynamics are observed as the system is cooled to within a few degrees of the glass transition temperature (T(g)). We observe nonexponential relaxation of the frequency-frequency correlation function, similar to what has been reported for other dynamical correlation functions. In addition, we see evidence for α-like relaxation, typically associated with long-time, cooperative molecular motion, on the ultrafast time scale. The data suggests that the spectral dynamics are sensitive to cooperative motion occurring on time scales that are necessarily longer than the observation time.

Local-mode approach to modeling multidimensional infrared spectra of metal carbonyls

We present a general approach for modeling multidimensional infrared spectra based on a combination of phenomenological fitting and ab initio electronic structure calculations. The vibrational Hamiltonian is written in terms of bilinearly coupled Morse oscillators that represent local carbonyl stretches. This should be contrasted with the previous approach, where the anharmonic Hamiltonian was given in terms of normal-mode coordinates ( Baiz et al. J. Phys. Chem. A 2009 , 113 , 9617 ). The bilinearly coupled Morse oscillator Hamiltonian is parametrized such that the frequencies and couplings are consistent with experiment, and the anharmonicities are computed by density functional theory. The advantages of the local-mode versus normal-mode approaches are discussed, as well as the ability of different density functionals to provide accurate estimates of the model parameters. The applicability and usefulness of the new approach are demonstrated in the context of the recently measured multidimensional infrared spectra of dimanganese decacarbonyl. The shifts in local site frequencies, couplings, and anharmonicities due to hydrogen bonding to the individual carbonyls are explored. It is found that, even though the effect of hydrogen bonding is nonlocal, it is additive.

Molecular theory and simulation of coherence transfer in metal carbonyls and its signature on multidimensional infrared spectra

We present a general and comprehensive theoretical and computational framework for modeling ultrafast multidimensional infrared spectra of a vibrational excitonic system in liquid solution. Within this framework, we describe the dynamics of the system in terms of a quantum master equation that can account for population relaxation, dephasing, coherence-to-coherence transfer, and coherence-to-population transfer. A unique feature of our approach is that, in principle, it does not rely on any adjustable fitting parameters. More specifically, the anharmonic vibrational Hamiltonian is derived from ab initio electronic structure theory, and the system-bath coupling is expressed explicitly in terms of liquid degrees of freedom whose dynamics can be obtained via molecular dynamics simulations. The applicability of the new approach is demonstrated by employing it to model the recently observed signatures of coherence transfer in the two-dimensional spectra of dimanganese decacarbonyl in liquid cyclohexane. The results agree well with experiment and shed new light on the nature of the molecular interactions and dynamics underlying the spectra and the interplay between dark and bright states, their level of degeneracy, and the nature of their interactions with the solvent.

Solvent-hindered intramolecular vibrational redistribution

Ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations of Mn(2)(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 T(1) relaxation rate unchanged.

Solvent-dependent spectral diffusion in a hydrogen bonded "vibrational aggregate"

Two-dimensional infrared spectroscopy (2DIR) is used to measure the viscosity-dependent spectral diffusion of a model vibrational probe, Mn(2)(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.

Ultrabroadband detection of a mid-IR continuum by chirped-pulse upconversion

Ultrabroadband mid-IR continuum pulses can be detected on a single-shot basis using chirped-pulse upconversion. Converting the mid-IR pulse to the visible reduces the fractional bandwidth and enables use of a silicon CCD camera. Removing the cross-phase modulation of the chirped pulse results in 1 cm(-1) resolution over a 600 cm(-1) detected bandwidth.