Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO)5.

Excitation of iron pentacarbonyl [Fe(CO)5], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO)5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental-theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M2,3-edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO)5. We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO)4 along the lowest singlet and triplet surfaces to form Fe(CO)3. Calculations of transient spectra in both singlet and triplet states as well as spin-orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO)4. Thus, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO)5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm.

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S-S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S-S and C-S bond cleavage, which involve the lowest excited electronic states S1 and S2. Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S1-S6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S-S and C-S bond dissociation coordinates. Finally, we present calculations of spin-orbit coupling constants at the Franck-Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S0-S6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck-Condon region and multiconfigurational character toward the bond-dissociation limits. We compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Generation and Implementation of Continuum Infrared Pulses for Broadband Detection in 2D IR Spectroscopy.

In recent years, new methods of generating continuum mid-infrared pulses through filamentation in gases have been developed for ultrafast time-resolved infrared vibrational spectroscopy. The generated infrared pulses can have thousands of wavenumbers of bandwidth, spanning the entire mid-IR region while retaining pulse length below 100 fs. This technology has had a significant impact on problems involving ultrafast structural dynamics in congested spectra with broad features, such as those found in aqueous solutions and molecules with strong intermolecular interactions. This study describes the recent advances in generating and characterizing these pulses and the practical aspects of implementing these sources for broadband detection in transient absorption and 2D IR spectroscopy.

Comparing the structures and photophysical properties of two charge transfer co-crystals.

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S0 state and ionic character in the S1 state. The high degree of charge transfer in the S1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S1 transitions.

Excited-State Dynamics during Primary C-I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy.

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to nπC═O* states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed nπC═O* and nσC-I* character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin-orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C-I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N4,5 edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C-I homolysis. By examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C-I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental-theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin-orbit coupled system.

Dramatic Conformer-Dependent Reactivity of the Acetaldehyde Oxide Criegee Intermediate with Dimethylamine Via a 1,2-Insertion Mechanism.

The reactivity of carbonyl oxides has previously been shown to exhibit strong conformer and substituent dependencies. Through a combination of synchrotron-multiplexed photoionization mass spectrometry experiments (298 K and 4 Torr) and high-level theory [CCSD(T)-F12/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ with an added CCSDT(Q) correction], we explore the conformer dependence of the reaction of acetaldehyde oxide (CH3CHOO) with dimethylamine (DMA). The experimental data support the theoretically predicted 1,2-insertion mechanism and the formation of an amine-functionalized hydroperoxide reaction product. Tunable-vacuum ultraviolet photoionization probing of anti- or anti- + syn-CH3CHOO reveals a strong conformer dependence of the title reaction. The rate coefficient of DMA with anti-CH3CHOO is predicted to exceed that for the reaction with syn-CH3CHOO by a factor of ∼34,000, which is attributed to submerged barrier (syn) versus barrierless (anti) mechanisms for energetically downhill reactions.

Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO)4 photodissociation at 261 nm.

We employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin-orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

Core spectroscopy of oxazole.

We have measured, analyzed, and simulated the ground state valence photoelectron spectrum, x-ray absorption (XA) spectrum, x-ray photoelectron (XP) spectrum as well as normal and resonant Auger-Meitner electron (AE) spectrum of oxazole at the carbon, oxygen, and nitrogen K-edge in order to understand its electronic structure. Experimental data are compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second-order and time-dependent density functional levels of theory. We demonstrate (1) that both N and O K-edge XA spectra are sensitive to the amount of dynamical electron correlation included in the theoretical description and (2) that for a complete description of XP spectra, additional orbital correlation and orbital relaxation effects need to be considered. The normal AE spectra are dominated by a singlet excitation channel and well described by theory. The resonant AE spectra, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.

A New Pathway for Intersystem Crossing: Unexpected Products in the O(3P) + Cyclopentene Reaction.

We investigated the reaction of O(3P) with cyclopentene at 4 Torr and 298 K using time-resolved multiplexed photoionization mass spectrometry, where O(3P) radicals were generated by 351 nm photolysis of NO2 and reacted with excess cyclopentene in He under pseudo-first-order conditions. The resulting products were sampled, ionized, and detected by tunable synchrotron vacuum ultraviolet radiation and an orthogonal acceleration time-of-flight mass spectrometer. This technique enabled measurement of both mass spectra and photoionization spectra as functions of time following the initiation of the reaction. We observe propylketene (41%), acrolein + ethene (37%), 1-butene + CO (19%), and cyclopentene oxide (3%), of which the propylketene pathway was previously unidentified experimentally and theoretically. The automatically explored reactive potential energy landscape at the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level and the related master equation calculations predict that cyclopentene oxide is formed on the singlet potential energy surface, whereas propylketene is first formed on the triplet surface. These calculations provide evidence that significant intersystem crossing can happen in this reaction not only around the geometry of the initial triplet adduct but also around that of triplet propylketene. The formation of 1-butene + CO is initiated on the triplet surface, with bond cleavage and hydrogen transfer occurring during intersystem crossing to the singlet surface. At present, we are unable to explain the mechanistic origins of the acrolein + ethene channel, and we thus refrain from assigning singlet or triplet reactivity to this channel. Overall, at least 60% of the products result from triplet reactivity. We propose that the reactivity of cyclic alkenes with O(3P) is influenced by their greater effective degree of unsaturation compared with acyclic alkenes. This work also suggests that searches for minimum-energy crossing points that connect triplet surfaces to singlet surfaces should extend beyond the initial adducts.

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films.

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. Herein, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N-N and axial NO2 stretching modes with the long-lived, excited phonon bath.

Ultraviolet photodissociation of gas-phase iron pentacarbonyl probed with ultrafast infrared spectroscopy.

It is well known that ultraviolet photoexcitation of iron pentacarbonyl results in rapid loss of carbonyl ligands leading to the formation of coordinatively unsaturated iron carbonyl compounds. We employ ultrafast mid-infrared transient absorption spectroscopy to probe the photodissociation dynamics of gas-phase iron pentacarbonyl following ultraviolet excitation at 265 and 199 nm. After photoexcitation at 265 nm, our results show evidence for sequential dissociation of iron pentacarbonyl to form iron tricarbonyl via a short-lived iron tetracarbonyl intermediate. Photodissociation at 199 nm results in the prompt production of Fe(CO)3 within 0.25 ps via several energetically accessible pathways. An additional 15 ps time constant extracted from the data is tentatively assigned to intersystem crossing to the triplet manifold of iron tricarbonyl or iron dicarbonyl. Mechanisms for formation of iron tetracarbonyl, iron tricarbonyl, and iron dicarbonyl are proposed and theoretically validated with one-dimensional cuts through the potential energy surface as well as bond dissociation energies. Ground state calculations are computed at the CCSD(T) level of theory and excited states are computed with EOM-EE-CCSD(dT).

Experimental and computational studies of Criegee intermediate reactions with NH3 and CH3NH2.

Ammonia and amines are emitted into the troposphere by various natural and anthropogenic sources, where they have a significant role in aerosol formation. Here, we explore the significance of their removal by reaction with Criegee intermediates, which are produced in the troposphere by ozonolysis of alkenes. Rate coefficients for the reactions of two representative Criegee intermediates, formaldehyde oxide (CH2OO) and acetone oxide ((CH3)2COO) with NH3 and CH3NH2 were measured using cavity ring-down spectroscopy. Temperature-dependent rate coefficients, k(CH2OO + NH3) = (3.1 ± 0.5) × 10-20T2 exp(1011 ± 48/T) cm3 s-1 and k(CH2OO + CH3NH2) = (5 ± 2) × 10-19T2 exp(1384 ± 96/T) cm3 s-1 were obtained in the 240 to 320 K range. Both the reactions of CH2OO were found to be independent of pressure in the 10 to 100 Torr (N2) range, and average rate coefficients k(CH2OO + NH3) = (8.4 ± 1.2) × 10-14 cm3 s-1 and k(CH2OO + CH3NH2) = (5.6 ± 0.4) × 10-12 cm3 s-1 were deduced at 293 K. An upper limit of ≤2.7 × 10-15 cm3 s-1 was estimated for the rate coefficient of the (CH3)2COO + NH3 reaction. Complementary measurements were performed with mass spectrometry using synchrotron radiation photoionization giving k(CH2OO + CH3NH2) = (4.3 ± 0.5) × 10-12 cm3 s-1 at 298 K and 4 Torr (He). Photoionization mass spectra indicated production of NH2CH2OOH and CH3N(H)CH2OOH functionalized organic hydroperoxide adducts from CH2OO + NH3 and CH2OO + CH3NH2 reactions, respectively. Ab initio calculations performed at the CCSD(T)(F12*)/cc-pVQZ-F12//CCSD(T)(F12*)/cc-pVDZ-F12 level of theory predicted pre-reactive complex formation, consistent with previous studies. Master equation simulations of the experimental data using the ab initio computed structures identified submerged barrier heights of -2.1 ± 0.1 kJ mol-1 and -22.4 ± 0.2 kJ mol-1 for the CH2OO + NH3 and CH2OO + CH3NH2 reactions, respectively. The reactions of NH3 and CH3NH2 with CH2OO are not expected to compete with its removal by reaction with (H2O)2 in the troposphere. Similarly, losses of NH3 and CH3NH2 by reaction with Criegee intermediates will be insignificant compared with reactions with OH radicals.

Real-Time Probing of Electron Dynamics Using Attosecond Time-Resolved Spectroscopy.

Attosecond science has paved the way for direct probing of electron dynamics in gases and solids. This review provides an overview of recent attosecond measurements, focusing on the wealth of knowledge obtained by the application of isolated attosecond pulses in studying dynamics in gases and solid-state systems. Attosecond photoelectron and photoion measurements in atoms reveal strong-field tunneling ionization and a delay in the photoemission from different electronic states. These measurements applied to molecules have shed light on ultrafast intramolecular charge migration. Similar approaches are used to understand photoemission processes from core and delocalized electronic states in metal surfaces. Attosecond transient absorption spectroscopy is used to follow the real-time motion of valence electrons and to measure the lifetimes of autoionizing channels in atoms. In solids, it provides the first measurements of bulk electron dynamics, revealing important phenomena such as the timescales governing the switching from an insulator to a metallic state and carrier-carrier interactions.

Collective vibrations of water-solvated hydroxide ions investigated with broadband 2DIR spectroscopy.

The infrared spectra of aqueous solutions of NaOH and other strong bases exhibit a broad continuum absorption for frequencies between 800 and 3500 cm(-1), which is attributed to the strong interactions of the OH(-) ion with its solvating water molecules. To provide molecular insight into the origin of the broad continuum absorption feature, we have performed ultrafast transient absorption and 2DIR experiments on aqueous NaOH by exciting the O-H stretch vibrations and probing the response from 1350 to 3800 cm(-1) using a newly developed sub-70 fs broadband mid-infrared source. These experiments, in conjunction with harmonic vibrational analysis of OH(-)(H2O)n (n = 17) clusters, reveal that O-H stretch vibrations of aqueous hydroxides arise from coupled vibrations of multiple water molecules solvating the ion. We classify the vibrations of the hydroxide complex by symmetry defined by the relative phase of vibrations of the O-H bonds hydrogen bonded to the ion. Although broad and overlapping spectral features are observed for 3- and 4-coordinate ion complexes, we find a resolvable splitting between asymmetric and symmetric stretch vibrations, and assign the 2850 cm(-1) peak infrared spectra of aqueous hydroxides to asymmetric stretch vibrations.

A phenomenological approach to modeling chemical dynamics in nonlinear and two-dimensional spectroscopy.

We present an approach for calculating nonlinear spectroscopic observables, which overcomes the approximations inherent to current phenomenological models without requiring the computational cost of performing molecular dynamics simulations. The trajectory mapping method uses the semi-classical approximation to linear and nonlinear response functions, and calculates spectra from trajectories of the system's transition frequencies and transition dipole moments. It rests on identifying dynamical variables important to the problem, treating the dynamics of these variables stochastically, and then generating correlated trajectories of spectroscopic quantities by mapping from the dynamical variables. This approach allows one to describe non-Gaussian dynamics, correlated dynamics between variables of the system, and nonlinear relationships between spectroscopic variables of the system and the bath such as non-Condon effects. We illustrate the approach by applying it to three examples that are often not adequately treated by existing analytical models--the non-Condon effect in the nonlinear infrared spectra of water, non-Gaussian dynamics inherent to strongly hydrogen bonded systems, and chemical exchange processes in barrier crossing reactions. The methods described are generally applicable to nonlinear spectroscopy throughout the optical, infrared and terahertz regions.

Observation of a Zundel-like transition state during proton transfer in aqueous hydroxide solutions.

It is generally accepted that the anomalous diffusion of the aqueous hydroxide ion results from its ability to accept a proton from a neighboring water molecule; yet, many questions exist concerning the mechanism for this process. What is the solvation structure of the hydroxide ion? In what way do water hydrogen bond dynamics influence the transfer of a proton to the ion? We present the results of femtosecond pump-probe and 2D infrared experiments that probe the O-H stretching vibration of a solution of dilute HOD dissolved in NaOD/D(2)O. Upon the addition of NaOD, measured pump-probe transients and 2D IR spectra show a new feature that decays with a 110-fs time scale. The calculation of 2D IR spectra from an empirical valence bond molecular dynamics simulation of a single NaOH molecule in a bath of H(2)O indicates that this fast feature is due to an overtone transition of Zundel-like H(3)O(2)(-) states, wherein a proton is significantly shared between a water molecule and the hydroxide ion. Given the frequency of vibration of shared protons, the observations indicate the shared proton state persists for 2-3 vibrational periods before the proton localizes on a hydroxide. Calculations based on the EVB-MD model argue that the collective electric field in the proton transfer direction is the appropriate coordinate to describe the creation and relaxation of these Zundel-like transition states.

Proton transfer in concentrated aqueous hydroxide visualized using ultrafast infrared spectroscopy.

While it is generally recognized that the hydroxide ion can rapidly diffuse through aqueous solution due to its ability to accept a proton from a neighboring water molecule, a description of the OH(-) solvation structure and mechanism of proton transfer to the ion remains controversial. In this report, we present the results of femtosecond infrared spectroscopy measurements of the O-H stretching transition of dilute HOD dissolved in NaOD/D(2)O. Pump-probe, photon echo peak shift, and two-dimensional infrared spectroscopy experiments performed as a function of deuteroxide concentration are used to assign spectral signatures that arise from the OH(-) ion and its solvation shell. A spectral feature that decays on a ∼110 fs time scale is assigned to the relaxation of transiently formed configurations wherein a proton is equally shared between a HOD molecule and an OD(-) ion. Over picosecond waiting times, features appear in 2D IR spectra that are indicative of the exchange of population between OH(-) ions and HOD molecules due to deuteron transfer. The construction of a spectral model that includes spectral relaxation, chemical exchange, and thermalization processes, and self-consistently treats all of our data, allows us to qualitatively explain the results of our experiments and gives a lower bound of 3 ps for the deuteron transfer kinetics.

Ultrafast 2D IR anisotropy of water reveals reorientation during hydrogen-bond switching.

Rearrangements of the hydrogen bond network of liquid water are believed to involve rapid and concerted hydrogen bond switching events, during which a hydrogen bond donor molecule undergoes large angle molecular reorientation as it exchanges hydrogen bonding partners. To test this picture of hydrogen bond dynamics, we have performed ultrafast 2D IR spectral anisotropy measurements on the OH stretching vibration of HOD in D(2)O to directly track the reorientation of water molecules as they change hydrogen bonding environments. Interpretation of the experimental data is assisted by modeling drawn from molecular dynamics simulations, and we quantify the degree of molecular rotation on changing local hydrogen bonding environment using restricted rotation models. From the inertial 2D anisotropy decay, we find that water molecules initiating from a strained configuration and relaxing to a stable configuration are characterized by a distribution of angles, with an average reorientation half-angle of 10°, implying an average reorientation for a full switch of ≥20°. These results provide evidence that water hydrogen bond network connectivity switches through concerted motions involving large angle molecular reorientation.

Structural rearrangements in water viewed through two-dimensional infrared spectroscopy.

Compared with other molecular liquids, water is highly structured because of its ability to form up to four hydrogen bonds, resulting in a tetrahedral network of molecules. However, this underlying intermolecular structure is constantly in motion, exhibiting large fluctuations and reorganizations on time scales from femtoseconds to picoseconds. These motions allow water to play a key role in a number of chemical and biological processes. By exploiting the fact that the OH stretching frequency of dilute HOD in liquid D(2)O is highly dependent upon the configuration of the neighbor nearest to the proton, researchers have been able to track water's time-dependent structure using two-dimensional infrared (2D IR) spectroscopy, which tags molecules at an initial frequency and then watches as that frequency evolves with respect to time. Recent advances in molecular dynamics simulation techniques allow for the calculation of 2D IR spectra, providing an atomistic interpretation tool of 2D IR spectra in terms of the underlying dynamics of the liquid. In this Account, we review recent ultrafast 2D IR studies at MIT that provide new information on the mechanism of hydrogen-bond rearrangements in liquid water. The 2D IR spectra of the OH stretching vibration of HOD in D(2)O appear highly asymmetric. In the frequency range indicative of hydrogen-bonded molecules (<3300 cm(-1)), the 2D spectra remain fairly compact. By contrast, in the frequency range in which molecules having weak or broken hydrogen bonds absorb (>3500 cm(-1)), the 2D spectra broaden over a time scale of approximately 60 fs, consistent with librations (hindered rotations) of water molecules. This broadening indicates that molecules forming weak or broken hydrogen bonds are unstable and reorient rapidly to return to a hydrogen-bonded configuration. These conclusions are supported by the results of molecular dynamics simulations, which suggest that water molecules undergo a large-angle reorientation during the course of hydrogen-bond exchange. The transition state for hydrogen-bond rearrangements is found to resemble a bifurcated hydrogen bond. Roughly half of the hydrogen-bond exchange events in the simulation are found to involve the insertion of a water molecule across a hydrogen bond, suggesting that hydrogen-bond exchange in water involves the correlated motion of water molecules as far away as the second solvation shell. The combination of ultrafast 2D IR spectroscopy with simulation-based modeling is leading to self-consistent descriptions of the underlying dynamics in liquid water. Moreover, these results also demonstrate a more general, unique characteristic of the spectroscopy: if a spectral signature of the transition state exists, then 2D IR can effectively serve as a transition-state spectroscopy.

Sub-picosecond to Sub-nanosecond Vibrational Energy Transfer Dynamics in Pentaerythritol Tetranitrate.

The time scale associated with shock-induced detonation is a key property of energetic materials that remains poorly understood. Herein, we test aspects of one potential mechanism, the phonon up-pumping mechanism, where shock compression excites lattice phonon modes, transferring energy to intramolecular vibrations leading to chemical bond cleavage and reaction. Using ultrafast infrared pump-probe spectroscopy on pentaerythritol tetranitrate (PETN), we reveal sub-picosecond vibrational energy transfer (VET) from the photoexcited band at 1660 cm-1 into every other infrared-active mode in the probed frequency range 800-1800 cm-1. Energy transfer processes remain incomplete at 150 ps. Computational predictions from density functional theory are used in tandem to elucidate VET pathways in PETN.