Vibrationally Controlled Chemistry: Mode- and Bond-Selected Reaction of CH3D with Cl

Ultrafast dynamics of CN radical reactions with chloroform solvent under vibrational strong coupling

Polariton chemistry may provide a new means to control molecular reactivity, permitting remote, reversible modification of reaction energetics, kinetics, and product yields. A considerable body of experimental and theoretical work has already demonstrated that strong coupling between a molecular vibrational mode and the confined electromagnetic field of an optical cavity can alter chemical reactivity without external illumination. However, the mechanisms underlying cavity-altered chemistry remain unclear in large part because the experimental systems examined previously are too complex for detailed analysis of their reaction dynamics. Here, we experimentally investigate photolysis-induced reactions of cyanide radicals with strongly-coupled chloroform (CHCl3) solvent molecules and examine the intracavity rates of photofragment recombination, solvent complexation, and hydrogen abstraction. We use a microfluidic optical cavity fitted with dichroic mirrors to facilitate vibrational strong coupling (VSC) of the C-H stretching mode of CHCl3 while simultaneously permitting optical access at visible wavelengths. Ultrafast transient absorption experiments performed with cavities tuned on- and off-resonance reveal that VSC of the CHCl3 C-H stretching transition does not significantly modify any measured rate constants, including those associated with the hydrogen abstraction reaction. This work represents, to the best of our knowledge, the first experimental study of an elementary bimolecular reaction under VSC. We discuss how the conspicuous absence of cavity-altered effects in this system may provide insights into the mechanisms of modified ground state reactivity under VSC and help bridge the divide between experimental results and theoretical predictions in vibrational polariton chemistry.

Stretching-mode specificity in the Cl + CH3D(v1-I, v1-II, and v4 = 1; |jK〉) reactions: dependency on the initial |jK〉 selectivity

The title reactions were studied at a collisional energy of 5.4 kcal mol^-1 in a crossed-beam product-imaging experiment. The dynamics attributes of the dominant ground-state CH₂D(0₀) and the accompanied C-D bend-excited CH₂D(6₁) products were imaged in reactions with totally 16 ro-vibrationally selected states of the CH₃D(v_i, |jK〉) reagents. We found that all three vibrational excitations yielded marked |jK〉-dependent rate-enhancements in forming the (0₀, 0/1)_s product pairs. Furthermore, for a given rotational |jK〉-mode, a vibrational-mode propensity of v₄ > v₁-I > v₁-II in rate promotion and a clear manifestation of the Fermi-phase-induced interference effect of the latter two were observed. Compared to the reactivity of the rotationless state |jK〉 = |00〉, a minute rotational-excitation of all three stretch-excited CH₃D(v_i = 1) reagents could yield significantly higher reaction rates for the product pair (0₀, 0)_s, but not so for (0₀, 1)_s. The signals in forming the (6₁, 0)_s pair were clearly notable but smaller than that of the ground-state reaction product pair, (0₀, 0)_g. An opposite propensity of v₁-II ≈ v₁-I > v₄ with a milder dependency on the initial |jK〉-states was observed. The angular distributions of the (0₀, 0)_s pairs were nearly identical for all ro-vibrationally excited reagents, displaying the typical trait for a direct abstraction of the rebound mechanism. Similar distributions were found for the (6₁, 0)_s pairs; yet, both pairs deviated substantially from the peripheral feature of the ground-state reaction pair of (0₀, 0)_g. Those of the (0₀, 1)_s pairs in reactions with v₄-excitation featured a prominent forward-peaking distribution-suggestive of a time-delayed, resonance-mediated pathway, again with little dependency on the initial |jK〉-states. As for the reactions with the two Fermi-dyads, v₁-I and v₁-II, albeit showing globally similar distributions to that for v₄, a substantial variation with the initial rotational-mode excitation could be discerned in the forward-peaking features. To unravel the impact of the Fermi-phase on the |jK〉-dependent attributes, we adopted a comparative approach by contrasting the observations in reactions with the Fermi-dyad reagents (the superposition states) to those with the pure-state reagents. Remarkable distinctions are unveiled and elucidated with the unexplained results explicitly pointed out, which call for future theoretical investigations for deeper understanding.

Imaging the Mode-Specificity in Cl + CH3D(v1-I, v1-II, v4 = 1; |jK⟩ = |10⟩) → CH2D(41) + HCl(v)

We report a crossed-beam imaging experiment on the title reactions at two collisional energies (E_c) of 5.3 and 10 kcal mol^-1. Both the integral cross sections relative to the ground-state reactivity and the differential cross sections were measured and compared. We found that one-quantum excitations of the CH₃-stretching vibrations of the CH₃D reagent exerted profound mode-specificity in forming the umbrella-mode-excited CH₂D(4₁) products with the vibrational efficacy of v₄ > v₁-I > v₁-II at both E_c values. The concomitantly formed HCl(v) coproducts were vibrationally cold. Interestingly, the branching ratios of (v = 1)/(v = 0) appeared invariant to the initial stretch-modes of excitation at E_c = 5.3 kcal mol^-1, yet exhibited a pronounced mode-specific dependency in the order of v₁-II > v₁-I > v₄ at E_c = 10.3 kcal mol^-1. This large and E_c-dependent disparity between the two Fermi-coupled reagents, v₁-I and v₁-II, is particularly significant and could be another facet─in addition to that in the recently reported vibrational enhancement factors─of the Fermi-phase-induced interference effect manifested in the product vibrational branching ratio. The pair-correlated angular distributions (v_CH2D, v_HCl)_s = (4₁, 0)_s in the three stretch-excited reactions were globally alike and resembled that of the ground-state reaction pair (0₀, 0)_g, suggestive of a direct abstraction mechanism of the peripheral type. This is in sharp contrast to all other vibrationally excited pairs of (1₁, 0)_s, (3₁, 0)_s, and (6₁, 0)_s previously reported in the CH₂D + HCl isotopic channel, for which both the direct abstraction and a time-delayed resonance pathway partake.

Toward a Comprehensive Understanding of Mode-Specific Dynamics of Polyatomic Reactions: A Full-Dimensional Quantum Dynamics Study of the H + NH3 Reaction

Mode specificity not only sheds light on reaction dynamics but also opens the door for chemical reaction control. This work reports a state-of-the-art full-dimensional quantum dynamics study on the prototypical hydrogen abstraction reaction of hydrogen with ammonia, which serves as a benchmark for advancing our fundamental understanding of polyatomic reaction dynamics. By taking advantage of the (3 + 1) Radau-Jacobi coordinates, the bond-specific probabilities are resolved with the reactant NH₃ initiated from either a non-degenerate or degenerate stretching vibrational state. The observed different atom-specific abstraction probabilities from individual states of the degenerate pair are rationalized in the local mode representation according to the different vibrational energy deposited in each N-H bond. It is verified that the three H atoms in NH₃ have the same abstraction probability as that from the degenerate pair and the linear combination of the degenerate pair gives the same reaction probability. In addition, the symmetric and asymmetric stretching modes of the reactant NH₃ have comparable efficacies on driving the reaction.

Pair-Correlated Imaging of Cl + CH3D(v4, v1-I, v1-II = 1, |jK⟩) → CH2D(vi) + HCl(v)

The title reactions were studied at a collisional energy of 10.0 kcal mol^-1 in a crossed-beam, product-imaging experiment. In terms of integral cross sections, all three CH₃-stretching excited CH₃D(v_CH3 = 1) reagents promote the reactivity in forming the predominant product pair of (v_CH2D, v_HCl)_s = (0₀, 0/1)_s with a prominent mode-propensity of v₄ > v₁-I > v₁-II, where v₄ denotes the degenerate mode of CH₃ asymmetric stretch and v₁-I and v₁-II are a pair of Fermi-coupled, symmetric-stretch states. The vibrationally excited CH₂D product pairs of (6₁, 0)_s, (1₁, 0)_s, and (3₁, 0)_s appear to be minor channels and display a reverse propensity of v₄ < v₁-I ≈ v₁-II for (6₁, 0)_s, while v₄ > v₁-I for (1₁, 0)_s. Based on the observed angular distributions, we conjecture that, irrespective of the initial mode of excitation, the (0₀, 0)_s product pair proceeds by a direct abstraction of the peripheral type, whereas the (0₀,1)_s pair is mediated via a resonance pathway. Intriguingly, the angular distributions of the excited product pairs-(6₁, 0)_s, (1₁, 0)_s, and (3₁, 0)_s-are remarkably similar and comprise the traits of both the peripheral mechanism and resonance pathway. Possible interpretation and implication are suggested. In addition, due to the spectral overlap of the REMPI bands and heavily congested image features, a robust data analysis method is developed, which enables us to extract the dynamics attributes of a weak feature buried in the proximate, more intense ones with high fidelity.

Imaging pair-correlated reaction cross sections in F + CH3D(νb = 0, 1) → CH2D(ν4 = 1) + HF(ν)

The title reactions were studied in a crossed-beam experiment at collisional energies (Ec) from 0.5 to 4.7 kcal mol-1. The νb (ν4) vibrational mode denotes the bending (umbrella) motion of the CH3D reactant (CH2D product). Using a time-sliced, velocity-map imaging technique, we extracted the state-specific, pair-correlated integral and differential cross sections. As with other isotopically analogous ground-state reactions, an inverted vibrational population of the HF coproduct was observed. Both the step-like excitation function near the threshold and the oscillatory forward-backward peakings in the Ec-evolution of the two dominant pair-correlated angular distributions at lower Ec suggest a resonance-mediated, time-delay mechanism. As Ec increases, the angular distribution of the HF(ν = 2) product evolves into a smooth and broad swath in the backward hemisphere, indicative of a direct rebound mechanism. One quantum excitation of the bending modes of CH3D(νb = 1) promotes the reaction rate by two- to three-fold up to Ec = 2.1 kcal mol-1. Broadly speaking, all major findings are qualitatively in line with previous results in the reactions of the F atom with other isotopologues. However, the rainbow feature recently observed in the CH2D(00) + HF(ν = 3) product channel is entirely absent. A possible rationale is put forward, which reinforces the previous reactive rainbow conjecture and calls for future theoretical investigations.

Mode specific dynamics in the H2 + SH → H + H2S reaction

Full-dimensional quantum dynamics and quasi-classical trajectory studies indicate strong mode selectivity in the H₂ + SH reaction.

Quantum state resolved gas–surface reaction dynamics experiments: a tutorial review

This tutorial review provides an introduction and overview of quantum state resolved reactivity measurements for methane chemisorption on transition metal surfaces.

Vibrational excitation induced proton transfer in hydrated Nafion membranes

We study the energy relaxation and structural relaxation dynamics of hydrated protons in Nafion membranes at different hydration levels using femtosecond infrared transient absorption spectroscopy. At low hydration levels we observe that the excitation of the proton vibration of an Eigen-like proton hydration structure leads to a structural relaxation process in which the Eigen-like structure evolves to a Zundel-like proton hydration structure. This reorganization leads to a transfer of the proton charge and closely follows the mechanism of infrared-induced adiabatic proton transfer that has been proposed by S. Hammes-Schiffer, J. T. Hynes, and others. At high hydration levels, the spectral dynamics are dominated by vibrational energy relaxation and subsequent cooling of the proton hydration structures and the surrounding water molecules. Using a kinetic analysis of the transient spectral data, we determine the rates of proton transfer, vibrational energy relaxation, and cooling as a function of hydration level. We find that infrared-induced proton transfer occurs at all hydration levels but becomes less observable at high hydration levels due to the increasingly dominant influence of the vibrational energy relaxation.

Infrared-activated proton transfer in aqueous nafion proton-exchange-membrane nanochannels

We report on the observation of a strong reorganization of the proton hydration structure in hydrated Nafion membranes following single-quantum excitation of a proton vibration with ∼4 μm light pulses. The reorganization takes place with a time constant of 170 ± 20  fs and leads to a strong red shift of the excited proton vibration and the rise of new waterlike O-H stretch absorption bands. These observations can be explained from a vibrational-excitation-induced change of the proton-hydration structure that involves transfer of the proton charge. The results are consistent with recent quantum molecular dynamics simulations of proton transfer in Nafion membranes.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.

Mode Selectivity for a "Central" Barrier Reaction: Eight-Dimensional Quantum Studies of the O((3)P) + CH4 → OH + CH3 Reaction on an Ab Initio Potential Energy Surface

The dynamics of a combustion reaction, namely, O((3)P) + CH4 → OH + CH3, is investigated with an eight-dimensional quantum model that includes representatives of all vibrational modes of CH4 and with a full-dimensional quasi-classical trajectory (QCT) method. The calculated excitation functions for the ground vibrational state CH4 agree well with experiment. Both quantum and QCT results suggest that excitation of the stretching modes of CH4 enhances the reaction, while the bending and umbrella modes have a smaller impact on reactivity, again consistent with experimental findings. However, none of the vibrational excitations has comparable efficiency in promoting the reaction as translational energy.

Theoretical Study of the Validity of the Polanyi Rules for the Late-Barrier Cl + CHD3 Reaction

The Polanyi rules, which state that vibrational energy is more efficient in promoting a late-barrier reaction than translational energy, were questioned recently by an experimental unexpected finding that the CH stretch excitation is no more effective in promoting the late-barrier Cl + CHD3 reaction than the translational energy. However, the present quantum dynamics study on the best-available potential energy surface for the title reaction reveals that the CH stretch excitation does promote the reaction significantly, except at low collision energies. Further studies should be carried out to solve the disagreements between theory and experiment on the reaction.

Gaussian binning of the vibrational distributions for the Cl + CH4(v(4/2) = 0, 1) → H + CH3Cl(n(1)n(2)n(3)n(4)n(5)n(6)) reactions

We test several binning techniques to obtain mode-specific final-state distributions for polyatomic reactions. Normal mode analysis is done after an exact transformation to the Eckart frame. Standard histogram binning (HB) and three different variants of the energy-based Gaussian binning (1GB) are employed to obtain the probabilities of the vibrational states. We consider the two major issues of the polyatomic quasiclassical product analysis, i.e., (1) rounding the classical action to the nearest integer can result in unphysical states and (2) the normal-mode analysis can break down for highly distorted geometries. We show that 1GB can handle issue 1 when the total vibrational energy is evaluated in the normal mode space using the harmonic approximation and both issues 1 and 2 can be solved when the total vibrational energy is calculated exactly in the Cartesian space. We found that anharmonicity in the quantized energy levels does not have a significant effect on the final-state distributions. Quasiclassical trajectory calculations are performed for the reactant ground-state and bending-excited Cl((2)P(3/2)) + CH(4)(v(4/2) = 0, 1) → H + CH(3)Cl reactions using an ab initio potential energy surface. The product analysis techniques are successfully applied to the CH(3)Cl product molecules and some qualitative features of the results are discussed.

Dynamics of the O(3P) + CHD3(vCH = 0,1) reactions on an accurate ab initio potential energy surface

Recent experimental and theoretical studies on the dynamics of the reactions of methane with F and Cl atoms have modified our understanding of mode-selective chemical reactivity. The O + methane reaction is also an important candidate to extend our knowledge on the rules of reactivity. Here, we report a unique full-dimensional ab initio potential energy surface for the O((3)P) + methane reaction, which opens the door for accurate dynamics calculations using this surface. Quasiclassical trajectory calculations of the angular and vibrational distributions for the ground state and CH stretching excited O + CHD(3)(v(1) = 0,1) → OH + CD(3) reactions are in excellent agreement with the experiment. Our theory confirms what was proposed experimentally: The mechanistic origin of the vibrational enhancement is that the CH-stretching excitation enlarges the reactive cone of acceptance.

Quasiclassical trajectory calculations of correlated product distributions for the F + CHD3(v1 = 0, 1) reactions using an ab initio potential energy surface

We report quasiclassical trajectory (QCT) calculations of the correlated product distributions and branching ratios of the reactions F+CHD(3)(v(1)=0,1)-->HF(v)+CD(3)(v) and DF(v)+CHD(2)(v) using a recently published ab initio-based full-dimensional global potential energy surface [G. Czako et al., J. Chem. Phys. 130, 084301 (2009)]. Harmonic normal mode analysis is done for the methyl products to determine the classical actions of each normal mode and then standard histogram binning and Gaussian binning (GB) methods are employed to obtain quantum state-resolved probabilities of the products. QCT calculations have been performed for both the vibrationally ground state and the CH stretching excited F+CHD(3)(v(1)=0,1) reactions at eight different collision energies in the 0.5-7.0 kcal/mol range. HF and DF vibrationally state-resolved rotational and angular distributions, CD(3) and CHD(2) mode-specific vibrational distributions, and correlated vibrationally state-specific distributions for the product pairs have been obtained and the correlated results were compared to the experiment. We find that the use of GB can be advantageous especially in the threshold regions. The CH stretching excitation in the reactant does not change the CD(3) vibrational distributions significantly, whereas the HF molecules become vibrationally and rotationally hotter. On the other hand in the case of the DF+CHD(2) channel the initially excited CH stretch appears mainly "intact" in the CHD(2) product and the DF distributions are virtually the same as formed from the ground state CHD(3) reaction. The computed results qualitatively agree with recent crossed molecular beam experiment [W. Zhang et al., Science 325, 303 (2009)] that (a) CHD(2)(v(1)=1) is the most populated product state of the F+CHD(3)(v(1)=1) reaction and this reaction produces much less CHD(2)(v=0) compared to the reaction F+CHD(3)(v=0); (b) the cross section ratio of CHD(2)(v(1)=1):CHD(2)(v=0) formed from the reactions F+CHD(3)(v(1)=1):F+CHD(3)(v=0) is less than 1 and shows little collision energy dependency; (c) the reactant CH stretch excitation increases the DF:HF ratio at low collision energies; (d) the correlated vibrational and angular distributions for DF(v)+CHD(2)(v(1)=0,1) from the ground state and stretch-excited reactions, respectively, are almost identical.

Chemical dynamics of vibrationally excited molecules: Controlling reactions in gases and on surfaces

Experimental studies of the chemical reaction dynamics of vibrationally excited molecules reveal the ability of different vibrations to control the course of a reaction. This Perspective describes those studies for the prototypical reaction of vibrationally excited methane and its isotopologues in gases and on surfaces and looks to the prospects of similar studies in liquids. The influences of vibrational excitation on the C-H bond cleavage in a single collision reaction with Cl and in dissociative adsorption on a Ni surface bear some striking similarities. Both reactions are bond-selective processes in which the initial preparation of a molecular eigenstate containing a large component of C-H stretching results in preferential cleavage of that bond. It is possible to cleave either the C-H bond or C-D bond in the reaction of Cl with CH3D, CH2D2, or CHD3 and, similarly, to use initial excitation of the C-H stretch to promote dissociation of CHD3 to CD3 and H on a Ni surface. Different vibrational modes, such as the symmetric and antisymmetric stretches in CH3D or CH4, lead to very different reactivities, and molecules with the symmetric stretching vibration excited can be as much as 10 times more reactive than ones with the antisymmetric stretch excited. The origin of this behavior lies in the change in the vibrational motion induced by the interaction with the atomic reaction partner or the surface.

Tracking the energy flow along the reaction path

We report a comprehensive study of the quantum-state correlation property of product pairs from reactions of chlorine atoms with both the ground-state and the CH stretch-excited CHD(3). In light of available ab initio theoretical results, this set of experimental data provides a conceptual framework to visualize the energy-flow pattern along the reaction path, to classify the activity of different vibrational modes in a reactive encounter, to gain deeper insight into the concept of vibrational adiabaticity, and to elucidate the intermode coupling in the transition-state region. This exploratory approach not only opens up an avenue to understand polyatomic reaction dynamics, even for motions at the molecular level in the fleeting transition-state region, but it also leads to a generalization of Polanyi's rules to reactions involving a polyatomic molecule.

Probing remnants of invariants to mediate energy exchange in highly chaotic many-dimensional systems

A technique is presented to scrutinize a piece of remnants of invariants [R. B. Shirts and W. P. Reinhardt, J. Chem. Phys. 77, 15 (1982)] buried in chaos in many degrees of freedom (DOF) dynamical systems in terms of canonical perturbation theory based on Lie transforms. The transformed canonical variables are often evaluated by the truncation of the coordinate transformation at a finite order in the original Hamiltonian system. However, the truncation of canonical variables gives rise to a loss of the symplectic property of the system. This results in apparent abrupt fluctuation of the action integrals, which yields a misinterpretation. We demonstrate, in a three-DOF Hamiltonian system of HCN isomerization reaction, that our technique can detect remnants of invariants buried in the potential well even at energies higher than the potential barrier, although the conventional truncation scheme fails to do so. This technique makes it possible to shed light on the physical insight into how the reactive mode exchanges its energy with the other modes and through which resonance the energy exchange takes place in reacting systems.

Ab initio design of picosecond infrared laser pulses for controlling vibrational-rotational excitation of CO molecules

Optimal control of rovibrational excitations of the CO molecule using picosecond infrared laser pulses is described in the framework of the electric-nuclear Born-Oppenheimer approximation [G. G. Balint-Kurti et al., J. Chem. Phys. 122, 084110 (2005)]. The potential energy surface of the CO molecule in the presence of an electric field is calculated using coupled cluster theory with a large orbital basis set. The quantum dynamics of the process is treated using a full three dimensional treatment of the molecule in the laser field. The detailed mechanisms leading to efficient control of the selected excitation processes are discussed.

Disentangling mode-specific reaction dynamics from overlapped images

The hydrogen abstraction reaction between atomic chlorine and C-H stretch-excited CHD(3) was studied under crossed-beam conditions. Prior to collisions, an infrared (IR) laser was used to pump up a fraction of CHD(3) to nu(1) = 1. A time-sliced velocity imaging technique was exploited to image the recoil velocity distribution of the state-selected product CD(3)(nu = 0). For energetic reasons, the IR-on image shows severely overlapped features arising from both the excited and the un-pumped ground-state reagents. A novel threshold method was then developed to directly determine the fraction of IR-excited CHD(3) reagents, which in turn enables us to disentangle the state-selected dynamics from the overlapped images. The results reveal significant differences from previous experimental reports.