Ultrafast vibrational energy redistribution within C–H and O–H stretching modes of liquid methanol

Probing the adsorption configuration of methanol at a charged air/aqueous interface using nonlinear spectroscopy

Investigating the effects of electrolyte ions on the adsorption configuration of methanol at a charged interface is important for studying the interface structure of electrolyte solutions and the oxidation mechanism of methanol in fuel cells. This study uses sum frequency generation (SFG) and heterodyne-detected second harmonic generation (HD-SHG) to investigate the adsorption configuration of methanol at the air/aqueous interface of 0.1 M NaClO4 solution, 0.1 M HClO4 solution and pure water. The results elucidate that the ion effect in the electrolyte solution affects the interface's charged state and the methanol's adsorption conformation at the interface. The negatively charged surface of the 0.1 M NaClO4 solution and the positively charged surface of the 0.1 M HClO4 solution arise from the corresponding specific ionic effects of the electrolyte solution. The orientation angle of methyl with respect to the surface normal is 43.4° ± 0.1° at the 0.1 M NaClO4 solution surface and 21.5° ± 0.2° at the 0.1 M HClO4 solution surface. Examining these adsorption configurations in detail, we find that at the negatively charged surface the inclined orientation angle (43.4°) of methanol favors the hydroxymethyl production by breaking the C-H bond, while at the positively charged surface the upright orientation angle (21.5°) of methanol promotes the methoxy formation by breaking the O-H bond. These findings not only illuminate the intricate ion effects on small organic molecules but also contribute to a molecular-level comprehension of the oxidation mechanism of methanol at electrode interfaces.

Intermolecular Interactions at the Silica-Liquid Interface Modulate the Fermi Resonance Coupling in Surface Methanol

The buried solid/liquid interface between hydrophilic fused silica and binary solvent mixtures of acetonitrile (MeCN) and methanol (MeOH) was studied with vibrational sum-frequency generation (vSFG) spectroscopy. Our data showed that at high relative concentrations of methanol, the Fermi resonance peak in the vSFG spectrum is greatly suppressed, and it progressively gains intensity as methanol is diluted with perdeuterated acetonitrile. This phenomenon is quantified by the Fermi resonance coupling coefficient, W, extracted using a two-level model, as well as the experimental intensity ratio, R, of the methyl Fermi resonance band to that of the symmetric stretch. At a 1.0 MeOH mole fraction, W and R values were 10 ± 10 cm^-1 and 0.01 ± 0.02, respectively, whereas at a 0.1 mole fraction, W and R increased to 46 ± 4 cm^-1 and 0.43 ± 0.16, respectively. This indicates that solvation with acetonitrile effectively tunes the Fermi coupling of methanol vibrations at the silica/liquid interface.

Study on hydrogen bonding network in aqueous methanol solution by Raman spectroscopy

The Raman spectra of aqueous methanol solution with various concentrations were measured at room temperature and atmospheric pressure. We found that the CO stretching vibration mode of methanol showed a significant blue shift at V_m (V_m represents the volume fraction of methanol) >0.4, while the CH symmetric and asymmetric stretching vibration modes exhibited red shift under the same conditions. These results illustrate that the variation of hydrogen bond between methanol and water molecules lead to a phase transition of the methanol-water complex at V_m = 0.4. Furthermore, the red shift of the CH vibration mode indicates that there is no hydrogen bond formed between the CH₃ group of methanol and water molecules. In addition, we found that the frequency shift of C-H is affected by the hydrogen bond C-O…H-O formed between methanol and water molecules, and the corresponding theoretical discussion is given. Finally, the phase transition process of methanol-water complex in methanol-water binary solution was given by theoretical analysis.

A "Universal" Spectroscopic Map for the OH Stretching Mode in Alcohols

Empirical maps are presented for the OH stretching vibrations in neat alcohols in which the relevant spectroscopic quantities are expressed in terms of the electric field exerted on the hydrogen atom by the surrounding liquid. It is found, by examination of the four lowest linear alcohols, methanol, ethanol, n-propanol, and n-butanol, that a single map can be used for alcohols with different alkyl groups. This "universal" map is in very good agreement with maps optimized for the individual alcohols but differs from those previously developed for water. This suggests that one map can be used for all alcohols, perhaps even those not examined in the present study. The universal map gives IR lineshapes in good agreement with measured spectra for isotopically dilute methanol and ethanol, while the two-dimensional IR photon echo spectra give results that differ from experiments. The role of non-Condon effects, reorientation dynamics, hydrogen bonding, and spectral diffusion is discussed.

High resolution imaging with differential infrared absorption micro-spectroscopy

Although confocal infrared (IR) absorption micro-spectroscopy is well established for far-field chemical imaging, its scope remains restricted since diffraction limits the spatial resolution to values a little above half the radiation wavelength. Yet, the successful implementations of below-the-diffraction limit far-field fluorescence microscopies using saturated irradiation patterns for example for stimulated-emission depletion and saturated structured-illumination suggest the possibility of using a similar optical patterning strategy for infrared absorption mapping at high resolution. Simulations are used to show that the simple mapping of the difference in transmitted/reflected IR energy between a saturated vortex-shaped beam and a Gaussian reference with a confocal microscope affords the generation of high-resolution vibrational absorption images. On the basis of experimentally relevant parameters, the simulations of the differential absorption scheme reveal a spatial resolution better than a tenth of the wavelength for incident energies about a decade above the saturation threshold. The saturated structured illumination concepts are thus expected to be compatible with the establishment of point-like point-spread functions for measuring the absorbance of samples with a scanning confocal microscope recording the differential transmission/reflection.

A framework for far-field infrared absorption microscopy beyond the diffraction limit

A framework is proposed for infrared (IR) absorption microscopy in the far-field with a spatial resolution below the diffraction limit. The sub-diffraction resolution is achieved by pumping a transient contrast in the population of a selected vibrational mode with IR pulses that exhibit alternating central minima and maxima, and by probing the corresponding absorbance at the same wavelength with adequately delayed Gaussian pulses. Simulations have been carried out on the basis of empirical parameters emulating patterned thin films of octadecyltrichlorosilane and a resolution of 250 nm was found when probing the CH₂ stretches at 3.5 μm with pump energies less than ten times the vibrational saturation threshold.

Mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon-tetrachloride mixtures II: excited state hydrogen bonding structure and dynamics, infrared emission spectrum, and excited state lifetime

We present a mixed quantum-classical molecular dynamics study of the hydrogen-bonding structure and dynamics of a vibrationally excited hydroxyl stretch in methanol/carbon-tetrachloride mixtures. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the ground and first-excited adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are determined by Hellmann-Feynman forces obtained by taking the expectation value of the force with respect to the ground or excited vibrational wave functions. Polarizable force fields are used which were previously shown to reproduce the experimental infrared absorption spectrum rather well, for different isotopomers and over a wide composition range [Kwac, K.; Geva, E. J. Phys. Chem. B 2011, 115, 9184]. We show that the agreement of the absorption spectra with experiment can be further improved by accounting for the dependence of the dipole moment derivatives on the configuration of the classical degrees of freedom. We find that the propensity of a methanol molecule to form hydrogen bonds increases upon photoexcitation of its hydroxyl stretch, thereby leading to a sizable red-shift of the corresponding emission spectrum relative to the absorption spectrum. Treating the relaxation from the first excited to the ground state as a nonadiabatic process, and calculating its rate within the framework of Fermi's golden rule and the harmonic-Schofield quantum correction factor, we were able to predict a lifetime which is of the same order of magnitude as the experimental value. The experimental dependence of the lifetime on the transition frequency is also reproduced. Nonlinear mapping relations between the hydroxyl transition frequency and bond length in the excited state and the electric field along the hydroxyl bond axis are established. These mapping relations make it possible to reduce the computational cost of the mixed quantum-classical treatment to that of a fully classical treatment.

A mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon tetrachloride mixtures: equilibrium hydrogen-bond structure and dynamics at the ground state and the infrared absorption spectrum

We present a mixed quantum-classical molecular dynamics study of the structure and dynamics of the hydroxyl stretch in methanol/carbon tetrachloride mixtures. One of the methanol molecules is tagged, and its hydroxyl stretch is treated quantum-mechanically, while the remaining degrees of freedom are treated classically. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the corresponding adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are in turn affected by the quantum-mechanical state of the tagged hydroxyl stretch via the corresponding Hellmann-Feynman forces. The ability of five different force-field combinations to reproduce the experimental absorption infrared spectrum of the hydroxyl stretch is examined for different isotopomers and on a wide range of compositions. It is found that, in addition to accounting for the anharmonic nature of the hydroxyl stretch, one also has to employ polarizable force fields and account for the damping of the polarizability at short distances. The equilibrium ground-state hydrogen-bonding structure and dynamics is analyzed, and its signature on the absorption infrared spectrum of the hydroxyl stretch is investigated in detail. Five different hydroxyl stretch subpopulations are identified and spectrally assigned: monomers (α), hydrogen-bond acceptors (β), hydrogen-bond donors (γ), simultaneous hydrogen-bond donors and acceptors (δ), and simultaneous hydrogen-bond donors and double-acceptors (ε). The fundamental transition frequencies of the α and β subpopulations are found to be narrowly distributed and to overlap, thereby giving rise to a single narrow band whose intensity is significantly diminished by rotational relaxation. The fundamental transition frequency distributions of the γ, δ, and ε subpopulations are found to be broader and to partially overlap, thereby giving rise to a single broad band which is red-shifted relative to the αβ band. The γδε band is also found to be inhomogeneously broadened and unaffected by rotational relaxation. The exchange rates between the different subpopulations and corresponding branching ratios are reported and explained. Finally, nonlinear mapping relations between the hydroxyl transition frequency and bond length and the electric field along the hydroxyl bond axis are established, which can be used to reduce the computational cost of the mixed quantum-classical treatment to that of a purely classical molecular dynamics simulation.

Efforts toward developing probes of protein dynamics: vibrational dephasing and relaxation of carbon-deuterium stretching modes in deuterated leucine

The spectral position of C-D stretching absorptions in the so-called "transparent window" of protein absorption (1800-2300 cm(-1)) makes them well suited as probes of protein dynamics with high temporal and structural resolution. We have previously incorporated single deuterated amino acids into proteins to site-selectively follow protein folding and ligand binding by steady-state FT IR spectroscopy. Ultimately, our goal is to use C-D bonds as probes in time-resolved IR spectroscopy to study dynamics and intramolecular vibrational energy redistribution (IVR) in proteins. As a step toward this goal, we now present the first time-resolved experiments characterizing the population and dephasing dynamics of selectively excited C-D bonds in a deuterated amino acid. Three differently deuterated, Boc-protected leucines were selected to systematically alter the number of additional C-D bonds that may mediate IVR out of the initially populated bright C-D stretching mode. Three-pulse photon echo experiments show that the steady-state C-D absorption linewidths are broadened by both homogeneous and inhomogeneous effects, and transient grating experiments reveal that IVR occurs on a subpicosecond time scale and is nonstatistical. The results have important implications for the interpretation of steady-state C-D spectra and demonstrate the potential utility of C-D bonds as probes of dynamics and IVR within a protein.

Vibrational Relaxation of Normal and Deuterated Liquid Nitromethane

Anti-Stokes Raman scattering is used to monitor vibrational energy redistribution in the ambient temperature liquids nitromethane (NM-h3) and perdeuterated nitromethane (NM-d3) after ultrafast IR excitation of either the symmetric or asymmetric CH- or CD-stretch transitions. The instantaneous populations of most of the fifteen NM vibrations are determined with good accuracy, and a global fitting procedure with a master equation is used to fit all the data. The pump pulses excite not only CH- or CD-stretches but also certain combinations of bending and nitro stretching fundamentals. The coupled vibrations that comprise the initial state are revealed via the instantaneous rise of the anti-Stokes transients associated with each vibrational fundamental. In contrast to many other polyatomic liquids studied previously, there is little energy exchange among the CH-stretch (or CD-stretch) excitations, which is attributed to the nearly free rotation of the methyl group in NM. The vibrational cooling process, which is the multistep return to a thermalized state, occurs in three stages in both NM-h3 and NM-d3. In the first stage, the parent CH- or CD-stretch decays in a few picoseconds, exciting all lower-energy vibrations. In the second stage, the midrange vibrations decay in 10-15 ps, exciting the lower-energy vibrations. In the third stage, these lower-energy vibrations decay into the bath in tens of picoseconds. The initial excitations are thermalized in approximately 150 ps in NM-h3 and there is little dependence on which CH-stretch is excited. VC is somewhat faster in NM-d3 with more dependence on the initial CD-stretch, taking approximately 100 ps with symmetric CD-stretch excitation and approximately 120 ps with asymmetric CD-stretch excitation. Comparison is made with earlier nonequilibrium molecular dynamics simulations of VC [Kabadi, V. N.; Rice, B. M. Molecular dynamics simulations of normal mode vibrational energy transfer in liquid nitromethane. J. Phys. Chem. A 2004, 108, 532-540]. The simulations do a good job of reproducing the observed VC process and in addition they predicted the slow interconversion among CH-stretch excitations and the slower relaxation of the asymmetric CH-stretch now observed here.

Vibrational relaxation of C-D stretching vibrations in CDCl3, CDBr3, and CDI3

We present time-resolved transient grating measurements of the vibrational relaxation rates of the C-D stretching vibrations of deuterated haloforms in benzene and acetone. We compare our results with previous measurements of excited C-H stretches in the same solvents to obtain insight into the solvent effect on the vibrational relaxation. In deuterated molecules, there are more low-order-coupled states and the states are closer in energy to the C-D stretch than in the unlabeled isotopologs. Therefore, the relaxation is faster for the deuterated molecules. The relaxation also shows a significant solvent dependence. Bromoform and iodoform form charge-transfer complexes with both benzene and acetone which enhance the relaxation rate. For chloroform, hydrogen bonding to acetone is expected to be a more favorable interaction. Surprisingly, however, the vibrational relaxation of CDCl(3) is slower in acetone than in benzene.

Classical vs Quantum Vibrational Energy Relaxation Pathways in Solvated Polyatomic Molecules

Vibrational energy relaxation (VER) of solvated polyatomic molecules can occur via different pathways. In this paper, we address the question of whether treating VER classically or quantum-mechanically can lead to different predictions with regard to the preferred pathway. To this end, we consider the relaxation of the singly excited asymmetric stretch of a rigid, symmetrical, and linear triatomic molecule (A-B-A) in a monatomic liquid. In this case, VER can occur either directly to the ground state or indirectly via intramolecular vibrational relaxation (IVR) to the symmetric stretch. We have calculated the rates of these two different VER pathways via classical mechanics and the linearized semiclassical (LSC) method. When the mass of the terminal A atoms is significantly larger than that of the central B atom, we find that LSC points to intermolecular VER as the preferred pathway, whereas the classical treatment points to IVR. The origin of this trend reversal appears to be purely quantum-mechanical and can be traced back to the significantly weaker quantum enhancement of solvent-assisted IVR in comparison to that of intermolecular VER.

Multidimensional infrared spectroscopy of the N-H bond motions in formamide

The heterodyned two-dimensional (2D) IR spectra and equilibrium dynamics of the N-H stretching motion of DCONHD in deuterated formamide, DCOND(2), were studied with 80 fs pulses at 3 microm. The time evolution of the heterodyned 2D IR spectra, pump-probe spectra, and photon echo peak shift demonstrate that interstate dynamics is occurring by relaxation of the original N-H excitation. The N-H vibrational frequency correlation function can be expressed as a sum of three exponentials with correlation times 0.24 ps, 0.8 ps, and 11 ps. The intermediate component is attributed to motions of the N-Hcdots, three dots, centeredO unit involving only slight angular variations of the N-H bond. The slow component is attributed to the structure breaking and making. The anisotropy decay confirmed that the significant angular N-H bond motion occurs on the 11 ps time scale. The fast component, which is the least well determined, might correspond to the modulation of the H-bond distance without angular motion. The correlation coefficient between the pumped and relaxed state distributions was +0.51, implying that the excited state phase memory is only slightly diminished by the relaxation of the N-H excitation. The relaxed modes are concluded to be local to the driven N-H mode.

Vibrational Energy Relaxation of Polyatomic Molecules in Liquid Solution via the Linearized Semiclassical Method

Vibrational energy relaxation (VER) of polyatomic, as opposed to diatomic, molecules can occur via different, often solvent assisted, intramolecular and/or intermolecular pathways. In this paper, we apply the linearized semiclassical (LSC) method for calculating VER rates in the prototypical case of a rigid, symmetrical and linear triatomic molecule (A-B-A) in a monatomic liquid. Starting at the first excited state of either the symmetric or asymmetric stretches, VER can occur either directly to the ground state or indirectly via intramolecular vibrational relaxation (IVR). The VER rate constants for the various pathways are calculated within the framework of the Landau-Teller formalism, where they are expressed in terms of two-time quantum-mechanical correlation functions. The latter are calculated by the LHA-LSC method, which puts them in a "Wignerized" form, and employs a local harmonic approximation (LHA) in order to compute the necessary multidimensional Wigner integrals. Results are reported for the LHL/Ar model of Deng and Stratt [J. Chem. Phys. 2002, 117, 1735], as well as for CO(2) in liquid argon and in liquid neon. The LHA-LSC method is shown to give rise to significantly faster VER and IVR rates in comparison to the classical treatment, particularly at lower temperatures. We also find that the type and extent of the quantum rate enhancement is strongly dependent on the particular VER pathway. Finally, we find that the classical and semiclassical treatments can give rise to opposite trends when it comes to the dependence of the VER rates on the solvent.

Quasiclassical Trajectory Simulations of OH(v) + NO2 → HONO2* → OH(v‘) + NO2: Capture and Vibrational Deactivation Rate Constants

Quasiclassical trajectory calculations are used to investigate the dynamics of the OH(v) + NO(2) --> HONO(2) --> OH(v') + NO(2) recombination/dissociation reaction on an analytic potential energy surface (PES) that gives good agreement with the known structure and vibrational frequencies of nitric acid. The calculated recombination rate constants depend only weakly on temperature and on the initial vibrational energy level of OH(v). The magnitude of the recombination rate constant is sensitive to the potential function describing the newly formed bond and to the switching functions in the PES that attenuate inter-mode interactions at long range. The lifetime of the nascent excited HONO(2) depends strongly not only on its internal energy but also on the identity of the initial state, in disagreement with statistical theory. This disagreement is probably due to the effects of slow intramolecular vibrational energy redistribution (IVR) from the initially excited OH stretching mode. The vibrational energy distribution of product OH(v') radicals is different from statistical distributions, a result consistent with the effects of slow IVR. Nonetheless, the trajectory results predict that vibrational deactivation of OH(v) via the HONO(2) transient complex is approximately 90% efficient, almost independent of initial OH(v) vibrational level, in qualitative agreement with recent experiments. Tests are also carried out using the HONO(2) PES, but assuming the weaker O-O bond strength found in HOONO (peroxynitrous acid). In this case, the predicted vibrational deactivation efficiencies are significantly lower and depend strongly on the initial vibrational state of OH(v), in disagreement with experiments. This disagreement suggests that the actual HOONO PES may contain more inter-mode coupling than found in the present model PES, which is based on HONO(2). For nitric acid, the measured vibrational deactivation rate constant is a useful proxy for the recombination rate, but IVR randomization of energy is not complete, suggesting that the efficacy of the proxy method must be evaluated on a case-by-case basis.

Quasi-Classical Trajectory Simulations of Intramolecular Vibrational Energy Redistribution in HONO2 and DONO2

By use of an analytic potential energy surface developed in this work for nitric acid, the quasi-classical trajectory method was used to simulate intramolecular vibrational energy redistribution (IVR). A method was developed for monitoring the average vibrational energy in the OH (or OD) mode that uses the mean-square displacement of the bond length calculated during the trajectories. This method is effective for both rotating and nonrotating molecules. The calculated IVR time constant for HONO(2) decreases exponentially with increasing excitation energy, is almost independent of rotational temperature, and is in excellent agreement with the experimental determination (Bingemann, D.; Gorman, M. P.; King, A. M.; Crim, F. F. J. Chem.Phys. 1997, 107, 661). In DONO(2), the IVR time constants show more complicated behavior with increasing excitation energy, apparently due to 2:1 Fermi-resonance coupling with lower frequency modes. This effect should be measurable in experiments.

Control of Raman lasing in the nonimpulsive regime

We explore coherent control of stimulated Raman scattering in the nonimpulsive regime. Optical pulse shaping of the coherent pump field leads to control over the stimulated Raman output. A model of the control mechanism is investigated.

Watching vibrational energy transfer in liquids with atomic spatial resolution

Ultrafast spectroscopy was used to study vibrational energy transfer between vibrational reporter groups on different parts of a molecule in a liquid. When OH stretching vibrations of different alcohols were excited by mid-infrared laser pulses, vibrational energy was observed to move through intervening CH2 or CH groups, taking steps up and down in energy, ending up at terminal CH3 groups. For each additional CH2 group in the path between OH and CH3, the time for vibrational energy transfer increased by about 0.4 picosecond.