A crossed-beam study of the F+HD-->DF+H reaction: the direct scattering channel

State-to-state differential cross sections of the title reaction are presented at four collision energies, ranging from 1.18 to 4.0 kcal /mol. Product angular distributions are predominantly backscattered at low energies and shift toward sideways (peaking near 150 degrees ) at higher energies. Experimental evidence for contributions from migratory trajectories was found in the more detailed angle-specific internal state distributions. The dynamics of this reaction is mostly governed by classical mechanics, and several major findings can qualitatively be rationalized. These "classical" behaviors serve as "references" and are to be contrasted to the attributes observed for the other isotopic product channel, HF+D, in a forthcoming paper.     

准经典轨线研究Li+HF(ν=0, j=0)→LiF+H反应动力学

基于Aguado等人拟合的APW势能面(PES), 运用准经典轨线(QCT)方法, 对反应Li+HF(ν=0, j=0)→LiF+H的动力学性质进行了计算. 主要研究了不同碰撞能条件下的反应截面、转动取向、产物散射角分布和竞争反应模式等. 结果表明, 该反应存在直接提取型和间接插入型两种反应模式, 在低能量下反应以间接插入反应模式为主, 能量大于200 meV时则以直接提取反应为主.     

Quantum Dynamics of Li+HF/DF Reaction Investigated by a State-to-State Time-dependent Wave Packet Approach

Using the reactant coordinate based time-dependent wave packet method, on the APW potential energy surface, the differential and integral cross sections of the Li+DF/HF(v=0, j=0, 1) reactions were calculated over the collision energy range from the threshold to 0.25 eV. The initial state-specified reaction rate constants of the title reaction were also calculated. The results indicate that, compared with the Li+DF reaction, the product LiF of Li+HF reaction is a little more rotationally excited but essentially similar. The initial rotational excitation from j=0 to 1 has little effect on the Li+DF reaction. However, the rotational excitation of DF does result in a little more rotationally excited product LiF. The different cross section of both reactions is forward biased in the studied collision energy range, especially at relatively high collision energy. The resonances in the Li+HF reaction may be identifiable as the oscillations in the product ro-vibrational state-resolved integral cross sections and backward scattering as a function of collusion energy. For the Li+HF reaction, the rate constant is not sensitive to the temperature and almost has no change in the temperature range considered. For the Li+DF reaction, the rate constant increase by a factor of about 10 in the temperature range of 100?300 K. Brief comparison for the total reaction probabilities and integral cross section of the Li+HF reaction has been carried out between ours and the values reported previously. The agreement is good, and the difference should come from the better convergence of our present calculations.     

Angular distributions for the F+H(2)-->HF+H reaction: the role of the F spin-orbit excited state and comparison with molecular beam experiments

We report quantum mechanical calculations of center-of-mass differential cross sections (DCS) for the F+H(2)-->HF+H reaction performed on the multistate [Alexander-Stark-Werner (ASW)] potential energy surfaces (PES) that describe the open-shell character of this reaction. For comparison, we repeat single-state calculations with the Stark-Werner (SW) and Hartke-Stark-Werner (HSW) PESs. The ASW DCSs differ from those predicted for the SW and HSW PES in the backward direction. These differences arise from nonadiabatic coupling between several electronic states. The DCSs are then used in forward simulations of the laboratory-frame angular distributions (ADs) measured by Lee, Neumark, and co-workers [J. Chem. Phys. 82, 3045 (1985)]. The simulations are scaled to match experiment over the range 12 degrees 相似文献   

Stereodynamics of the F + HD(v = 0, j = 1) reaction: direct vs. resonant mechanisms

The stereodynamics and mechanism of the F + HD(v = 0, j = 1) → HF (DF) + D (H) reactions have been thoroughly analysed at collision energies in the 0-160 meV range. Specifically, this study is focused on (i) the comparison between the stereodynamics of the collisions leading to HF and DF formation, and (ii) the stereodynamical fingerprints of the resonance that occurs at low collision energies in the HF channel and whose manifestation in the total cross section is greatly diminished for initial j > 0. While previous studies were limited to the analysis of integral cross sections (ICS), differential cross sections (DCS) and reaction probabilities, in the present work we have included the analysis of vectorial quantities such as the direction of the initial rotational angular momentum and internuclear axis, and their effect on reactivity. In particular, polarisation parameters (PP) and polarisation dependent differential cross sections (PDDCS), quantities that describe how the intrinsic HD rotational angular momentum and molecular axis polarisations contribute to reaction, are calculated and examined. The evolution of the PPs with the collision energy differs markedly between the two reaction channels. For the DF channel, the PP values are small and change very little in the energy range in which DF formation is appreciable. In contrast, rapid fluctuations in the magnitude and sign of the PPs are observed in the HF channel at low collision energies in and around the resonance. As the collision energy increases, direct (non-resonant) scattering prevails, and the various quantities are reasonably well accounted for by the QCT calculations, as in the case of the DF channel. The intrinsic directional information has been used to access the extent of control that can be achieved through polarisation of the HD molecule prior to collision. It was found that the same extrinsic preparation leads to very different outcomes on the HF channel DCS when the collision energy is close to the resonance. It is also shown that polarisation of the HD internuclear axis along the initial relative velocity enhances the effect of the resonance and allows its clear identification. Finally, the effect of different extrinsic preparations on the angle-velocity DCS is found to be strong, thus allowing considerable control of product angular distributions.     

Direct dynamics classical trajectory simulations of the O+ + CH4 reaction at hyperthermal energies

A Born-Oppenheimer direct dynamics simulation of the O(+) + CH(4) reaction dynamics at hyperthermal energies has been carried out with the PM3 (ground quartet state) Hamiltonian. Calculations were performed at various collision energies ranging from 0.5 to 10 eV with emphasis on high energy collisions where this reaction is relevant to materials erosion studies in low Earth orbit and geosynchronous Earth orbit. Charge transfer to give CH(4)(+) is the dominant channel arising from O(+) + CH(4) collisions in this energy range, but most of the emphasis in our study is on collisions that lead to reaction. All energetically accessible reaction channels were found, including products containing carbon-oxygen bonds, which is in agreement with the results of recent experiments. After correcting for compensating errors in competing reaction channels, our excitation functions show quantitative agreement with experiment (for which absolute magnitudes of cross sections are available) at high collision energies (several eV). More detailed properties, such as translational and angular distributions, show qualitative agreement. The opacity function reveals a high selectivity for producing OH(+) at high impact parameters, CH(3)(+)/CH(2)(+)/H(2)O(+) at intermediate impact parameters, and H(2)CO(+)/HCO(+)/CO(+) at small impact parameters. Angular distributions for CH(3)(+)/CH(2)(+)/OH(+) are forward scattered at high collision energies which implies the importance of direct reaction mechanisms, while reaction complexes play an important role at lower energies, especially for the H(2)O(+) product. Finally, we find that the nominally spin-forbidden product CH(3)(+) + OH can be produced by a spin-allowed pathway that involves the formation of the triplet excited product CH(3)(+)(?(3)E). This explains why CH(3)(+) can have a high cross section, even at very low collision energies. The results of this work suggest that the PM3 method may be applied directly to the study of O(+) reactions with small alkane molecules and polymer surfaces.     

碰撞能量对反应Sr＋HF转动取向的影响

在推广LEPS势能面上,用经典轨线方法,研究了反应碰撞能量对反应Sr＋HF的转动取向的影响.计算结果与产物轨道角动量模型进行比较.计算结果表明,随着碰撞能量的增加,产物转动取向越强烈.     

Quantum dynamics of the Li + HF --> H + LiF reaction at ultralow temperatures

Quantum-mechanical calculations are reported for the Li+HF(v=0,1,j=0)-->H+LiF(v',j') bimolecular scattering process at low and ultralow temperatures. Calculations have been performed for zero total angular momentum using a recent high-accuracy potential-energy surface for the X2A' electronic ground state. For Li+HF(v=0,j=0), the reaction is dominated by resonances due to the decay of metastable states of the Li cdots,...F-H van der Waals complex. Assignment of these resonances has been carried out by calculating the eigenenergies of the quasibound states. We also find that while chemical reactivity is greatly enhanced by vibrational excitation, the resonances get mostly washed out in the reaction of vibrationally excited HF with Li atoms. In addition, we find that at low energies, the reaction is significantly suppressed due to the less-efficient tunneling of the relatively heavy fluorine atom.     

The excitation function for Li + HF → LiF + H at collision energies below 80 meV

We have measured the dependence of the relative integral cross section of the reaction Li + HF → LiF + H on the collision energy (excitation function) using crossed molecular beams. By varying the intersection angle of the beams from 37° to 90° we covered the energy range 25 meV ≤ E(tr) ≤ 131 meV. We observe a monotonous rise of the excitation function with decreasing energy over the entire energy range indicating that a possible translational energy threshold to the reaction is significantly smaller than 25 meV. The steep rise is quantitatively recovered by a Langevin-type excitation function based on a vanishing threshold and a mean interaction potential energy ∝R(-2.5) where R is the distance between the reactants. To date all threshold energies deduced from ab initio potentials and zero-point vibrational energies are at variance with our results, however, our findings support recent quantum scattering calculations that predict significant product formation at collision energies far below these theoretical thresholds.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Rotationally resolved reactive scattering: imaging detailed Cl+C2H6 reaction dynamics

The hydrogen atom abstraction reaction of Cl (2P3/2) with ethane has been studied using the crossed molecular beam technique with dc slice imaging at collision energies from 3.2 to 10.4 kcal/mol. The products HCl (v,J) (v = 0, J = 0-5) were state-selectively detected using 2+1 resonance enhanced multiphoton ionization. The images were used to obtain the center-of-mass frame product angular distributions and translational energy release distributions. Two general features were found in all probed HCl quantum states at 6.7 kcal/mol collision energy, and these features have distinct translational energy release and angular distributions, as described for HCl (v = 0, J = 2) in a recent preliminary report [Li et al., J. Chem. Phys. 124, 011102 (2006)]. The results for HCl (v = 0, J = 2) at four collision energies were also compared to investigate the energy-dependent dynamics. We discuss the reaction in terms of a variety of models of polyatomic reaction dynamics. The dynamics of this well studied system are more complicated than can be accounted for by a single mechanism, and the results call for further theoretical and experimental investigations.     

Quasi-classical trajectory study of the S + OH → SO + H reaction: from reaction probability to thermal rate constant

First quasi-classical trajectory calculations have been carried out for the S((3)P) + OH(X?(2)Π) → SO(X?(3)Σ(-)) + H((2)S) reaction on an ab initio global potential energy surface for the ground electronic state, X?(2)A', of HSO. Cross sections, computed for collision energies up to 1 eV, show no energy threshold and decrease with the increasing collision energy. Rate constants have been calculated in the 5-500 K temperature range. The thermal rate constant is in good agreement with approximate quantum results, while a disagreement is found at 298 K with the experimental data. Product energy distributions have also been reported at four collision energies from 0.001 to 0.5 eV. The shapes of the rovibrational and angular distributions suggest the formation of an intermediate complex that is more and more long-lived as the collision energy increases.     

Monte Carlo simulation of RRKM unimolecular decomposition in molecular beam experiments. III. Application to angular distributions and branching fractions for the system K + CsF

The experiments on the reactive systems K + CsF and K + RbF by Stolte et al. are unique in that product angular distributions and branching fractions with variation of both translational and rotational energies have been obtained. The results for the reaction K + CsF with translational energy variation have now been re-evaluated using an algorithm recently published by us. A satisfactory fit to the experimental data is found by a long-lived collision complex model including strict angular momentum conservation in both the reactive and the non-reactive channels. No “osculating” behaviour has to be added to the model. Angular momentum limits have to be included in the model to reach agreement with the measured branching fractions and cross sections. This indicates that the internal dynamics of the complexes is important. The calculation shows that the reactants giving reactive complexes have large rotational and vibrational energies, while their translational energy is nearly average. The internal energy of reactive pairs of molecules decreases with increasing translational energy. The question how rapidly the branching fraction varies with translational energy is thus more complicated to answer than believed previously.     

Barrierless reactions between two closed-shell molecules. I. Dynamics of F(2)+CH(3)SCH(3) reaction

A detailed experimental and theoretical investigation of the first-reported barrierless reaction between two closed-shell molecules [J. Chem. Phys. 127, 101101 (2007)] is presented. The translational energy and angular distributions of two product channels, HF+CH(2)SFCH(3) and F+CH(3)SFCH(3), determined at several collision energies, have been analyzed to reveal the dynamics of the studied reaction. Detailed analysis of the experimental and computational results supports the proposed reaction mechanism involving a short-lived F-F-S(CH(3))(2) intermediate, which can be formed without any activation energy. Other possible reaction mechanisms have been discriminated. The decay of the intermediate and competition between the two product channels have been discussed.     

Low temperature kinetics, crossed beam dynamics and theoretical studies of the reaction S((1)D) + CH4 and low temperature kinetics of S((1)D) + C2H2

The reaction between sulfur atoms in the first electronically excited state, S((1)D), and methane (CH(4)), has been investigated in a complementary fashion in (a) crossed-beam dynamics experiments with mass spectrometric detection and time-of-flight (TOF) analysis at two collision energies (30.4 and 33.6 kJ mol(-1)), (b) low temperature kinetics experiments ranging from 298 K down to 23 K, and (c) electronic structure calculations of stationary points and product energetics on the CH(4)S singlet potential energy surface. The rate coefficients for total loss of S((1)D) are found to be very large (ca. 2 × 10(-10) cm(3) molec(-1) s(-1)) down to very low temperatures indicating that the overall reaction is barrier-less. Similar measurements are also performed for S((1)D) + C(2)H(2), and also for this system the rate coefficients are found to be very large (ca. 3 × 10(-10) cm(3) molec(-1) s(-1)) down to very low temperatures. From laboratory angular and TOF distributions at different product masses for the reaction S((1)D) + CH(4), it is found that the only open reaction channel at the investigated collision energies is that leading to SH + CH(3). The product angular, T(θ), and translational energy, P(E'(T)), distributions in the center-of-mass frame are derived. The reaction dynamics are discussed in terms of two different micromechanisms: a dominant long-lived complex mechanism at small and intermediate impact parameters with a strongly polarized T(θ), and a direct pickup-type (stripping) mechanism occurring at large impact parameters with a strongly forward peaked T(θ). Interpretation of the experimental results on the S((1)D) + CH(4) reaction kinetics and dynamics is assisted by high-level theoretical calculations on the CH(4)S singlet potential energy surface. The dynamics of the SH + CH(3) forming channel are compared with those of the corresponding channel (leading to OH + CH(3)) in the related O((1)D) + CH(4) reaction, previously investigated in crossed-beams in other laboratories at comparable collision energies. The possible astrophysical relevance of S((1)D) reactions with hydrocarbons, especially in the chemistry of cometary comae, is discussed.     

State-correlation matrix of the product pair from F + CD(4)--> DF(nu') + CD(3)(0 v(2) 0 0)

The title reaction was investigated under crossed-beam conditions at three different collision energies, E(c) = 8.4, 2.76 and 1.46 kcal mol(-1). The combination of using a (2 + 1) resonance-enhanced multiphoton ionization for tagging state-specific CD(3) products and exploiting a time-sliced velocity imaging for ion detection allows us to reveal the coincident information of the two product pairs in a state-correlated manner. The pair-correlated results are reported for the two product vibrators -- (v(2) = 0, v'), (v(2) = 1, v'), (v(2) = 2, v') and (v(2) = 3, v')-and the dynamics attributes we examined include product state distributions, energy disposals and angular distributions. Together with our earlier communications, a rather complete picture of the correlated dynamics of the title reaction emerges. One of the major findings, the anti-correlated excitations of the two product vibrators at all four energies of this study, can qualitatively be understood by kinematics arguments.     

Vibrational and translational enhancement of endothermic reactions: K + HCl(υ = 0,1) and K + HF(υ = 0,1)

The relative integral cross section for the two endothermic reactions K + HCl(υ = 0 and 1) → KCl + H and K + HF(υ = 0 and 1) → KF + H has been measured as a function of the collision energy E using the crossed molecular beam technique. The vibrationally excited state (υ = 1) has been populated thermally by heating the beam source to temperatures around 2000 K. The variation of the collision energy from thermal up to around 2.1 eV was achieved by seeding the K-beam with various carrier gases. The molecular reaction product was detected by surface ionization in connection with a time-of-flight method. The total energy threshold of the reactions has been found to be equal to or only slightly above the corresponding endothermicities. This suggests a vanishing or very low barrier crest on the potential energy hypersurfaces which is contradictory to recent theoretical results. The inclusion of tunneling in case of K + HF leads to a negligible rise of the barrier heights. The efficacy of translational and vibrational energy in promoting the reactive process has been directly compared over a wide range of collision energies. For K + HCl the vibrational enhancement of the reactivity descends with increasing E from approximately a factor of 10 at E = 0.08 eV to around unity for E ? 0.5 eV. The good agreement of this experimental result with phase space calculations suggests that the marked enhancements are predominantly caused by the long-range attraction between reagents in connection with an “early” barrier on the potential energy surface. In case of K + HF vibrational energy is by a factor of up to 380 more favourable in promoting the reaction than the same amount of translational energy. Again, with rising collision energy its efficacy decreases but promotes the reaction still by a factor of 70 at E = 1.7 eV. Since phase space theory fails here the reaction is certainly non-statistical and we conclude that the observed large efficacy of vibrational energy is due to a “late” barrier. The proposed barrier positions for the two systems are in accordance with theoretical results.     

Collision dynamics of O(3P) + DMMP using a specific reaction parameters potential form

Starting from previous benchmark CBS-QB3 electronic structure calculations (Conforti, P. F.; Braunstein, M.; Dodd, J. A. J. Phys. Chem. A 2009, 113, 13752), we develop two global potential energy surfaces for O((3)P) + DMMP collisions, using the specific reaction parameters approach. Each surface is simultaneously fit along the three major reaction pathways: hydrogen abstraction, hydrogen elimination, and methyl elimination. We then use these surfaces in classical dynamics simulations and compute reactive cross sections from 4 to 10 km s(-1) collision velocity. We examine the energy disposal and angular distributions of the reactive and nonreactive products. We find that for reactive collisions, an unusually large amount of the initial collision energy is transformed into internal energy. We analyze the nonreactive and reactive product internal energy distributions, many of which fit Boltzmann temperatures up to ~2000 K.