Time resolved infrared absorption studies of geminate recombination and vibrational relaxation in OClO photochemistry

Ultrafast time-resolved infrared absorption studies of aqueous chlorine dioxide (OClO) photochemistry are reported. Following photoexcitation at 401 nm, the evolution in optical density at frequencies between 1000 to 1100 cm(-1) is monitored to investigate vibrational energy deposition and relaxation along the asymmetric-stretch coordinate following the reformation of ground-state OClO via geminate recombination of the primary photofragments. The measured kinetics are compared to two proposed models for the vibrational-relaxation dynamics along the asymmetric-stretch coordinate. This comparison demonstrates that the perturbation model derived from molecular dynamics studies is capable of qualitatively reproducing the observed kinetics, where the collisional model employed in previous UV-pump, visible probe experiments demonstrates poor agreement with experiment. The ability of the perturbation model to reproduce the optical-density evolution observed in these studies demonstrates that for aqueous OClO, frequency dependence of the solvent-solute coupling is important in defining the level-dependent vibrational relaxation rates along the asymmetric-stretch coordinate. The absence of optical-density evolution corresponding to the population of higher vibrational levels (n>8) along the asymmetric-stretch coordinate suggests that following geminate recombination, energy is initially deposited into a local Cl-O stretch, with the relaxation of vibrational energy from this coordinate providing for delayed vibrational excitation of the asymmetric- and symmetric-stretch coordinates relative to geminate recombination, as previously observed.     

Intermolecular hydrogen bonding in chlorine dioxide photochemistry: A time-resolved resonance Raman study

The geminate-recombination and vibrational-relaxation dynamics of chlorine dioxide (OClO) dissolved in ethanol and 2,2,2-trifluoroethanol (TFE) are investigated using time-resolved resonance Raman spectroscopy. Stokes spectra are measured as a function of time following photoexcitation using degenerate pump and probe wavelengths of 398 nm. For OClO dissolved in ethanol, subpicosecond geminate recombination occurs resulting in the reformation of ground-state OClO with a quantum yield of 0.5±0.1. Following recombination, intermolecular-vibrational relaxation of OClO occurs with a time constant of 31±10 ps. For OClO dissolved in TFE, recombination occurs with a time constant of 1.8±0.8 ps and a quantum yield of only 0.3±0.1. The intermolecular-vibrational-relaxation time constant of OClO in TFE is 79±27 ps. The reduced geminate-recombination quantum yield, delayed recombination, and slower vibrational relaxation for OClO in TFE is interpreted in terms of greater self-association of the solvent. Degenerate pump–probe experiments are also presented that demonstrate decay of the Cl-solvent charge-transfer complex on the ∼1-ns time scale in ethanol and TFE. This time is significantly longer than the abstraction times observed for other systems demonstrating that Cl hydrogen abstraction from alcohols occurs in the presence of a significant energy barrier.     

Ultrafast anisotropy dynamics of water molecules dissolved in acetonitrile

Infrared pump-probe experiments are performed on isolated H(2)O molecules diluted in acetonitrile in the spectral region of the OH stretching vibration. The large separation between water molecules excludes intermolecular interactions, while acetonitrile as a solvent provides substantial hydrogen bonding. Intramolecular coupling between symmetric and asymmetric modes results in the anisotropy decay to the frequency-dependent values of approximately 0-0.2 with a 0.2 ps time constant. The experimental data are consistent with a theoretical model that includes intramolecular coupling, anharmonicity, and environmental fluctuations. Our results demonstrate that intramolecular processes are essential for the H(2)O stretching mode relaxation and therefore can compete with the intermolecular energy transfer in bulk water.     

Vibrational coherence transfer characterized with Fourier-transform 2D IR spectroscopy

Two-dimensional infrared (2D IR) spectroscopy of the symmetric and asymmetric C[Triple Bond]O stretching vibrations of Rh(CO)(2)acac in hexane has been used to investigate vibrational coherence transfer, dephasing, and population relaxation in a multilevel vibrational system. The transfer of coherence between close-lying vibrational frequencies results in extra relaxation-induced peaks in the 2D IR spectrum, whose amplitude depends on the coherence transfer rate. Coherence transfer arises from the mutual interaction of the bright CO stretches with dark states, which in this case reflects the mutual d-pi(*) back bonding of the Rh center to both the terminal carbonyls and the acetylacenonate ligand. For 2D IR relaxation experiments with variable waiting times, coherent dynamics lead to the modulation of peak amplitudes, while incoherent population relaxation and exchange results in the growth of the relaxation-induced peaks. We have modeled the data by propagating the density matrix with the Redfield equation, incorporating all vibrational relaxation processes during all three experimental time periods and including excitation reorientation effects arising from relaxation. Coherence and population transfer time scales from the symmetric to the asymmetric stretch were found to be 350 fs and 3 ps, respectively. We also discuss a diagrammatic approach to incorporating all vibrational relaxation processes into the nonlinear response function, and show how coherence transfer influences the analysis of structural variables from 2D IR spectroscopy.     

Isotopic effect on the vibrational lifetime of the carbon-deuterium stretch excitation on graphene

The relaxation of vibrational energy in the H and D stretch modes has been studied on the graphene surface using ab initio calculations. The dissipation of the vibrational energy stored in the stretching modes proceeds through vibration-phonon coupling, while the dissipation through electronic excitations makes only minor contributions. Recently, we reported the fast relaxation of the H stretch energy on graphene [S. Sakong and P. Kratzer, J. Chem. Phys. 133, 054505 (2010)]. Interestingly, we predict the lifetime of the D stretch to be markedly longer compared to the relaxation of the H stretch. This is unexpected since the vibrational amplitudes at carbon atoms in the joint C-D vibrational modes are larger than in the joint C-H modes, due to the mass ratio m(D)/m(C) > m(H)/m(C). However, the vibrational relaxation rate for the D stretch is smaller than for the H stretch, because the energy is dissipated to an acoustic phonon of graphene in the case of C-D rather than an optical phonon as is the case in C-H, and hence, the corresponding phonon density of states is lower in the C-D case. To rationalize our findings, we propose a general scheme for estimating vibrational lifetimes of adsorbates based on four factors: the density of states of the phonons that mediates the transitions, the vibration-phonon coupling strength, the anharmonic coupling between local modes, and the number of quanta involved in the transitions. Mainly the first two of these factors are responsible for the differences in the lifetimes of the C-H and C-D stretches. The possible role of the other factors is illustrated in the context of vibrational lifetimes in other recently studied systems.     

An H/D isotopic substitution study of the H5O2+.Ar vibrational predissociation spectra: exploring the putative role of Fermi resonances in the bridging proton fundamentals

To clarify the nature of the motions contributing to the observed multiplet structures in the low-energy (900-1800 cm-1) vibrational spectrum of the H5O2+ "Zundel" ion, we report the evolution of its vibrational fingerprint with sequential H/D isotopic substitution in a predissociation study of the Ar complexes. Of particular interest is the D4HO2+ complex, which displays a single intense band in the vicinity of the asymmetric OHO stretch of the bridging proton, in contrast to the more complex multiplet observed for both H5O2+ and D5O2+ isotopologues. These intensity patterns are consistent with the recent assignment of the bridging proton band's doublet in the H5O2+.Ne spectrum to a 2 x 2 Fermi resonance interaction between the shared proton stretch and a complex background level primarily derived from the O-O stretch together with two quanta of the wagging vibration involving the pyramidal deformations of the flanking H2O groups (Vendrell, O.; Gatti, F.; Meyer, H.-D. Angew. Chem., Int. Ed. 2007, 46, 6918). In addition, the observed trends rule out assignment of the approximately 1800 cm-1 feature in H5O2+ to a combination band of the bridging proton vibration with the O-O stretch, providing a secure foundation for the previously reported scheme that attributes this band to the out-of-phase intramolecular bending fundamental. The observed feature occurs at an unusually high energy for typical HOH bends, however, and we explore the participation of the bridging proton in these eigenstates by following how the calculated harmonic spectrum evolves when artificially large masses are assigned to the proton. The empirical assignments are supported by anharmonic estimates of the isotope shifts evaluated by the diffusion Monte Carlo method.     

High-resolution IR spectroscopy of dimers of HDO with H2O in helium nanodroplets

We report on rotationally resolved IR spectra of dimers of HDO as a deuterium (d) donor with H(2)O, HDO, and D(2)O embedded in superfluid Helium nanodroplets in the 2650-2660 and 2725-2740 cm(-1) regions of the O-D donor stretch and symmetric acceptor stretch vibrations, respectively. By comparing spectra at different levels of deuteration we were able to unambiguously assign the donor stretch signals of H(2)O···DOH, HDO···DOH, and D(2)O···DOH. For H(2)O···DOH, three ΔK(a) = 0 sub-bands were found that were assigned to transitions from the lower and upper acceptor switching states of K(a) = 0 and the lower acceptor switching state of K(a) = 1. In addition, b- and c-type transitions in the acceptor stretch region of HDO···DOH were observed that allowed us to determine the acceptor switching splitting of Δv? = 5.68 cm(-1) in the HDO···DOH vibrational ground state. We suggest that the dominating broadening mechanism is intervibrational relaxation due to coupling of the rovibrational levels of the chromophore via internal droplet excitations.     

Vibrational relaxation of CN stretch of pseudo-halide anions (OCN-, SCN-, and SeCN-) in polar solvents

The vibrational relaxation dynamics of pseudo-halide anions XCN- (X = O, S, Se) in polar solvents were studied to understand the effect of charge on solute-to-solvent intermolecular energy transfer (IET) and solvent assisted intramolecular vibrational relaxation (IVR) pathways. The T1 relaxation times of the CN stretch in these anions were measured by IR pump/IR probe spectroscopy, in which the 0-1 transition was excited, and the 0-1 and 1-2 transitions were monitored to follow the recovery of the ground state and decay of the excited state. For these anions in five solvents, H2O, D2O, CH3OH, CH3CN, and (CH3)2SO, relaxation rates followed the trend of OCN- > SCN- > SeCN-. For these anions and isotopes of SCN-, the relaxation rate was a factor of a few (2.5-10) higher in H2O than in D2O. To further probe the solvent isotope effect, the relaxation rates of S12C14N-, S13C14N-, and S12C15N- in deuterated methanols (CH3OH, CH3OD, CH3OH, CD3OD) were compared. Relaxation rate was found to be affected by the change of solvent vibrational band at the CN- stretching mode (CD3 symmetric stretch) and lower frequency regions, suggesting the presence of both direct IET and solvent assisted IVR relaxation pathways. The possible relaxation pathways and mechanisms for the observed trends in solute and solvent dependence were discussed.     

Dynamic structures of aqueous oxalate and the effects of counterions seen by 2D IR

Two dimensional vibrational echo spectra of oxalate in the carboxylate asymmetric stretch region in D(2)O show two transitions having anomalously slow spectral diffusion and a third transition having relaxation properties typical of the free carboxylate ion. Quantitative analysis of the frequency shifts of the carboxylate asymmetric stretch modes caused by a singly charged cation in the oxalate hydration shell supports that ion pairs can be responsible for these new transitions. Experimental evidence and DFT calculations are consistent with oxalate forming a mixture of "side-on" and "end on" contact ion pairs wherein the carboxylate groups are protected from mobile heavy water molecules.     

Vibrational dynamics of deoxyguanosine 5'-monophosphate following UV excitation

The relaxation dynamics of the DNA nucleotide deoxyguanosine 5'-monophosphate (dGMP) following 266 nm photoexcitation has been studied by transient IR spectroscopy with femtosecond time resolution. The induced dynamics of the amide I (carbonyl) stretch, the asymmetric guanine ring stretch and the phosphate asymmetric stretch are monitored in the region 1000-1800 cm(-1). Excitation and subsequent rapid internal conversion to a "hot" ground state is reflected by depletion of the vibrational ground states of the amide I stretch and guanine ring stretch. However, the vibrational ground state of the phosphate is left unperturbed, indicating the absence of vibrational coupling between the guanine ring system and the phosphate group. The vibrational ground state of the amide I is repopulated in 2.5 ps (±0.2 ps) while it takes 3.7 ps (±0.5 ps) to repopulate the guanine ring vibration. This article discusses two possible relaxation pathways of dGMP, as well as the implications of the weak phosphate dynamics.     

On the role of nonbonded interactions in vibrational energy relaxation of cyanide in water

The vibrationally excited cyanide ion (CN(-)) in H2O or D2O relaxes back to the ground state within several tens of picoseconds. Pump-probe infrared spectroscopy has determined relaxation times of T1 = 28 ± 7 and 71 ± 3 ps in H2O and D2O, respectively. Atomistic simulations of this process using nonequilibrium molecular dynamics simulations allow determination of whether it is possible at all to describe such a process, what level of accuracy in the force fields is required, and whether the information can be used to understand the molecular mechanisms underlying vibrational relaxation. It is found that, by using the best electrostatic models investigated, absolute relaxation times can be described rather more qualitatively (T1(H2O) = 19 ps and T1(D2O) = 34 ps) whereas the relative change in going from water to deuterated water is more quantitatively captured (factor of 2 vs 2.5 from experiment). However, moderate adjustment of the van der Waals ranges by less than 20% (for NVT) and 7.5% (for NVE), respectively, leads to almost quantitative agreement with experiment. Analysis of the energy redistribution establishes that the major pathway for CN(-) relaxation in H2O or D2O proceeds through coupling to the water-bending plus libration mode.     

Ultrafast vibrational energy relaxation of the water bridge

We report the energy relaxation of the OH stretch vibration of HDO molecules contained in an HDO:D(2)O water bridge using femtosecond mid-infrared pump-probe spectroscopy. We found that the vibrational lifetime is shorter (~630 ± 50 fs) than for HDO molecules in bulk HDO:D(2)O (~740 ± 40 fs). In contrast, the thermalization dynamics following the vibrational relaxation are much slower (~1.5 ± 0.4 ps) than in bulk HDO:D(2)O (~250 ± 90 fs). These differences in energy relaxation dynamics strongly indicate that the water bridge and bulk water differ on a molecular scale.     

The vibrational structure of the oxygen K-shell spectra in acenaphthenequinones: an ab initio study

The vibrational structure of the K-shell O1s → π? of acenaphthenequinone C(12)H(6)O(2) and its halogenated compound C(12)H(2)Br(2)Cl(2)O(2) has been simulated using an entirely ab initio approach. For both molecules, analysis of the calculated Franck-Condon factors confirm without ambiguity that, contrary to initial claims, the C-H stretching modes are not modified in the core states and are not excited. For C(12)H(6)O(2), the vibrational fine structure appears to be mainly due to three modes, involving C=O? asymmetric stretch and in-plane ring deformation modes, due to the symmetry breaking of the core state. For C(12)H(2)Br(2)Cl(2)O(2), the vibrational excitation arises essentially from the C=O? asymmetric stretch, with numerous secondary peaks arising from hot and combination bands. For both molecules, these bands are probably responsible for the asymmetry deduced in the experimental fits using a unique Morse potential and initially assigned to anharmonic effects.     

Nonperturbative vibrational energy relaxation effects on vibrational line shapes

A general formulation of nonperturbative quantum dynamics of solutes in a condensed phase is proposed to calculate linear and nonlinear vibrational line shapes. In the weak solute-solvent interaction limit, the temporal absorption profile can be approximately factorized into the population relaxation profile from the off-diagonal coupling and the pure-dephasing profile from the diagonal coupling. The strength of dissipation and the anharmonicity-induced dephasing rate are derived in Appendix A. The vibrational energy relaxation (VER) rate is negligible for slow solvent fluctuations, yet it does not justify the Markovian treatment of off-diagonal contributions to vibrational line shapes. Non-Markovian VER effects are manifested as asymmetric envelops in the temporal absorption profile, or equivalently as side bands in the frequency domain absorption spectrum. The side bands are solvent-induced multiple-photon effects which are absent in the Markovian VER treatment. Exact path integral calculations yield non-Lorentzian central peaks in absorption spectrum resulting from couplings between population relaxations of different vibrational states. These predictions cannot be reproduced by the perturbative or the Markovian approximations. For anharmonic potentials, the absorption spectrum shows asymmetric central peaks and the asymmetry increases with anharmonicity. At large anharmonicities, all the approximation schemes break down and a full nonperturbative path integral calculation that explicitly accounts for the exact VER effects is needed. A numerical analysis of the O-H stretch of HOD in D(2)O solvent reveals that the non-Markovian VER effects generate a small recurrence of the echo peak shift around 200 fs, which cannot be reproduced with a Markovian VER rate. In general, the nonperturbative and non-Markovian VER contributions have a stronger effect on nonlinear vibrational line shapes than on linear absorption.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

The relaxation of OH (v=1) and OD (v=1) by H2O and D2O at temperatures from 251 to 390 K

We report rate coefficients for the relaxation of OH(v=1) and OD(v=1) by H2O and D2O as a function of temperature between 251 and 390 K. All four rate coefficients exhibit a negative dependence on temperature. In Arrhenius form, the rate coefficients for relaxation (in units of 10(-12) cm3 molecule-1 s-1) can be expressed as: for OH(v=1)+H2O between 263 and 390 K: k=(2.4+/-0.9) exp((460+/-115)/T); for OH(v=1)+D2O between 256 and 371 K: k=(0.49+/-0.16) exp((610+/-90)/T); for OD(v=1)+H2O between 251 and 371 K: k=(0.92+/-0.16) exp((485+/-48)/T); for OD(v=1)+D2O between 253 and 366 K: k=(2.57+/-0.09) exp((342+/-10)/T). Rate coefficients at (297+/-1 K) are also reported for the relaxation of OH(v=2) by D2O and the relaxation of OD(v=2) by H2O and D2O. The results are discussed in terms of a mechanism involving the formation of hydrogen-bonded complexes in which intramolecular vibrational energy redistribution can occur at rates competitive with re-dissociation to the initial collision partners in their original vibrational states. New ab initio calculations on the H2O-HO system have been performed which, inter alia, yield vibrational frequencies for all four complexes: H2O-HO, D2O-HO, H2O-DO and D2O-DO. These data are then employed, adapting a formalism due to Troe (J. Troe, J. Chem. Phys., 1977, 66, 4758), in order to estimate the rates of intramolecular energy transfer from the OH (OD) vibration to other modes in the complexes in order to explain the measured relaxation rates-assuming that relaxation proceeds via the hydrogen-bonded complexes.     

Rotational dynamics of solvated carbon dioxide studied by infrared, Raman, and time-resolved infrared spectroscopies and a molecular dynamics simulation

Rotational dynamics of solvated carbon dioxide (CO(2)) has been studied. The infrared absorption band of the antisymmetric stretch mode in acetonitrile is found to show a non-Lorentzian band shape, suggesting a non-exponential decay of the vibrational and/or rotational correlation functions. A combined method of a molecular dynamics (MD) simulation and a quantum chemical calculation well reproduces the observed band shape. The analysis suggests that the band broadening is almost purely rotational, while the contribution from the vibrational dephasing is negligibly small. The non-exponential rotational correlation decay can be explained by a simple rotor model simulation, which can treat large angle rotations of a relatively small molecule. A polarized Raman study of the symmetric stretch mode in acetonitrile gives a rotational bandwidth consistent with that obtained from the infrared analysis. A sub-picosecond time-resolved infrared absorption anisotropy measurement of the antisymmetric stretch mode in ethanol also gives a decay rate that is consistent with the observed rotational bandwidths.     

Vibrational relaxation of the H2O bending mode in liquid water

We have studied the vibrational relaxation of the H(2)O bending mode in an H(2)O:HDO:D(2)O isotopic mixture using infrared pump-probe spectroscopy. The transient spectrum and its delay dependence reveal an anharmonic shift of 55+/-10 cm(-1) for the H(2)O bending mode, and a value of 400+/-30 fs for its vibrational lifetime.     

Structure of the aqueous solvated electron from resonance Raman spectroscopy: lessons from isotopic mixtures

The structure and thermodynamics of the hydrated electron are probed with resonance Raman spectroscopy of isotopic mixtures of H(2)O and D(2)O. The strongly enhanced intramolecular bends of e(-)(H(2)O) and e(-)(D(2)O) produce single downshifted bands, whereas the e(-)(HOD) bend consists of two components: one slightly upshifted from the 1,446 cm(-1) bulk frequency to 1,457 cm(-1) and the other strongly downshifted to approximately 1,396 cm(-1). This 60 cm(-1) split and the 200 (120) cm(-1) downshifts of the OH (OD) stretch frequencies relative to bulk water reveal that the water molecules that are Franck-Condon coupled to the electron are in an asymmetric environment, with one proton forming a strong hydrogen bond to the electron. The downshifted bend and librational frequencies also indicate significantly weakened torsional restoring forces on the water molecules of e(-)(aq), which suggests that the outlying proton is a poor hydrogen bond donor to the surrounding solvent. A 1.6-fold thermodynamic preference of the electron for H(2)O is observed based on the relative intensities of the e(-)(H(2)O) and e(-)(D(2)O) bands in a 50:50 isotopic mixture. This equilibrium isotope effect is consistent with the downshifted vibrational frequencies and a relative reduction of the zero-point energy of H(2)O bound to the electron. Our results enhance the cavity model of the solvated electron and support only those models that contain water monomers as opposed to other molecular species.