Crossed Molecular Beams and Theoretical Studies of the O(3P)+1,2-Butadiene Reaction: Dominant Formation of Propene+CO and Ethylidene+Ketene Molecular Channels

Detailed understanding of the mechanism of the combustion relevant multichannel reactions of O(³P) with unsaturated hydrocarbons (UHs) requires the identification of all primary reaction products, the determination of their branching ratios and assessment of intersystem crossing (ISC) between triplet and singlet potential energy surfaces (PESs). This can be best achieved combining crossed-molecular-beam (CMB) experiments with universal, soft ionization, mass-spectrometric detection and time-of-flight analysis to high-level ab initio electronic structure calculations of triplet/singlet PESs and RRKM/Master Equation computations of branching ratios (BRs) including ISC. This approach has been recently demonstrated to be successful for O(³P) reactions with the simplest UHs (alkynes, alkenes, dienes) containing two or three carbon atoms. Here, we extend the combined CMB/theoretical approach to the next member in the diene series containing four C atoms, namely 1,2-butadiene (methylallene) to explore how product distributions, branching ratios and ISC vary with increasing molecular complexity going from O(³P))+propadiene to O(³P)+1,2-butadiene. In particular, we focus on the most important, dominant molecular channels, those forming propene+CO (with branching ratio ∽0.5) and ethylidene+ketene (with branching ratio ∽0.15), that lead to chain termination, to be contrasted to radical forming channels (branching ratio ∽0.35) which lead to chain propagation in combustion systems.     

Prediction of product branching ratios in the C(3P)+C2H2-->l-C3H+H/c-C3H+H/C3+H2 reaction using ab initio coupled clusters calculations extrapolated to the complete basis set combined with Rice-Ramsperger-Kassel-Marcus and radiationless transition theories

Ab initio CCSD(T) calculations of intermediates and transition states on the singlet and triplet C3H2 potential energy surfaces extrapolated to the complete basis set limit are combined with statistical computations of energy-dependent rate constants of the C(3P)+C2H2 reaction under crossed molecular beam conditions. Rice-Ramsperger-Kassel-Marcus theory is applied for isomerization and dissociation steps within the same multiplicity and radiationless transition and nonadiabatic transition state theories are used for singlet-triplet intersystem crossing rates. The calculated rate constants are utilized to predict product branching ratios. The results demonstrate that, in qualitative agreement with available experimental data, c-C3H+H and C3+H2 are the most probable products at low collision energies, whereas l-C3H+H becomes dominant at higher Ec above approximately 25 kJ/mol.     

Crossed beam studies of elementary reactions of N and C atoms and CN radicals of importance in combustion

The dynamics of some elementary reactions of N(2D), C(3P,1D) and CN(X2 sigma +) of importance in combustion have been investigated by using the crossed molecular beam scattering method with mass spectrometric detection. The novel capability of producing intense, continuous beams of the radical reagents by a radio-frequency discharge beam source was exploited. From angular and velocity distribution measurements obtained in the laboratory frame, primary reaction products have been identified and their angular and translational energy distributions in the center-of-mass system, as well as branching ratios, have been derived. The dominant N/H exchange channel has been examined in the reaction N(2D) + CH4, which is found to lead to H + CH2NH (methylenimine) and H + CH3N (methylnitrene); no H2 elimination is observed. In the reaction N(2D) + H2O the N/H exchange channel has been found to occur via two competing pathways leading to HNO + H and HON + H, while formation of NO + H2 is negligible. Formation of H + H2CCCH (propargyl) is the dominant pathway, at low collision energy (Ec), of the C(3P) + C2H4 reaction, while at high Ec formation of the less stable C3H3 isomers (cyclopropenyl and/or propyn-1-yl) also occurs; the H2 elimination channel is negligible. The H elimination channel has also been found to be the dominant pathway in the C(3P,1D) + CH3CCH reaction leading to C4H3 isomers and, again, no H2 elimination has been observed to occur. In contrast, both H and H2 elimination, leading in comparable ratio to C3H + H and C3(X1 sigma g+) + H2(X1 sigma g+), respectively, have been observed in the reaction C(3P) + C2H2(X1 sigma g+). The occurrence of the spin-forbidden molecular pathway in this reaction, never detected before, has been rationalized by invoking the occurrence of intersystem crossing between triplet and singlet manifolds of the C3H2 potential energy surfaces. The reaction CN(X2 sigma +) + C2H2 has been found to lead to internally excited HCCCN (cyanoacetylene) + H. For all the reactions the dynamics have been discussed in the light of recent theoretical calculations on the relevant potential energy surfaces. Previous, lower resolution studies on C and CN reactions carried out using pulsed beams are noted. Finally, throughout the paper the relevance of these results to combustion chemistry is considered.     

Dynamics of the O(3P) + CH4 → OH + CH3 reaction is similar to that of a triatomic reaction

The O((3)P) + CH(4) reaction has been investigated using the quasi-classical trajectory (QCT) method and an ab initio pseudotriatomic potential energy surface (PES). This has been mainly motivated by very recent experiments which support the reliability of the triatomic modeling even at high collision energy ( = 64 kcal mol(-1)). The QCT results agree rather well with the experiments (translational and angular distributions of products); i.e., the ab initio pseudotriatomic modeling "captures" the essence of the reaction dynamics, although the PES was not optimized for high E(col). Furthermore, similar experiments on the O((3)P) + CD(4) reaction at moderate E(col) (12.49 kcal mol(-1)) have also been of a large interest here and, under these softer reaction conditions, the QCT method leads to results which are almost in quantitative agreement with experiments. The utility of the ab initio pseudotriatomic modeling has also been recognized for other analogous systems (X + CH(4)) but with very different PESs.     

Unimolecular decomposition of chemically activated pentatetraene (H2CCCCCH2) intermediates: A crossed beams study of dicarbon molecule reactions with allene

The reactions dynamics of the dicarbon molecule C2 in the 1Sigma (g)+ singlet ground state and 3Pi(u) first excited triplet state with allene, H2CCCH2(X1A1), was investigated under single collision conditions using the crossed molecular beam approach at four collision energies between 13.6 and 49.4 kJ mol(-1). The experiments were combined with ab initio electronic structure calculations of the relevant stationary points on the singlet and triplet potential energy surfaces. Our investigations imply that the reactions are barrier-less and indirect on both the singlet and the triplet surfaces and proceed through bound C5H4 intermediates via addition of the dicarbon molecule to the carbon-carbon double bond (singlet surface) and to the terminal as well as central carbon atoms of the allene molecule (triplet surface). The initial collision complexes isomerize to form triplet and singlet pentatetraene intermediates (H2CCCCCH2) that decompose via atomic hydrogen loss to yield the 2,4-pentadiynyl-1 radical, HCCCCCH2(X2B1). These channels result in symmetric center-of-mass angular distributions. On the triplet surface, a second channel involves the existence of a nonsymmetric reaction intermediate (HCCCH2CCH) that fragments through atomic hydrogen emission to the 1,4-pentadiynyl-3 radical [C5H3(X2B1)HCCCHCCH]; this pathway was found to account for the backward scattered center-of-mass angular distributions at higher collision energies. The identification of two resonance-stabilized free C5H3 radicals (i.e., 2,4-pentadiynyl-1 and 1,4-pentadiynyl-3) suggests that these molecules can be important transient species in combustion flames and in the chemical evolution of the interstellar medium.     

Dynamics of the O(3P) + C2H4 reaction: identification of five primary product channels (vinoxy, acetyl, methyl, methylene, and ketene) and branching ratios by the crossed molecular beam technique with soft electron ionization

The crossed molecular beam scattering technique with soft electron ionization (EI) is used to disentangle the complex dynamics of the polyatomic O(3P) + C2H4 reaction, which is of great relevance in combustion and atmospheric chemistry. Exploiting the newly developed capability of attaining universal product detection by using soft EI, at a collision energy of 54.0 kJ mol(-1), five different primary products have been identified, which correspond to the five exoergic competing channels leading to CH2CHO(vinoxy) + H, CH3CO(acetyl) + H, CH3(methyl) + HCO(formyl), CH2(methylene) + HCHO(formaldehyde), and CH2CO(ketene) + H2. From laboratory product angular and velocity distributions, center-of-mass product angular and translational energy distributions and the relative branching ratios for each channel have been obtained, affording an unprecedented characterization of this important reaction.     

Global ab initio potential energy surfaces for the O2(3Σg-)+O2(3Σg-) interaction

Completely ab initio global potential energy surfaces (PESs) for the singlet and triplet spin multiplicities of rigid O(2)((3)Σ(g)(-))+O(2)((3)Σ(g)(-)) are reported for the first time. They have been obtained by combining an accurate restricted coupled cluster theory with singles, doubles, and perturbative triple excitations [RCCSD(T)] quintet potential [Bartolomei et al., J. Chem. Phys. 128, 214304 (2008)] with complete active space second order perturbation theory (CASPT2) or, alternatively, multireference configuration interaction (MRCI) calculations of the singlet-quintet and triplet-quintet splittings. Spherical harmonic expansions, containing a large number of terms due to the high anisotropy of the interaction, have been built from the ab initio data. The radial coefficients of these expansions are matched at long range distances with analytical functions based on recent ab initio calculations of the electric properties of the monomers [M. Bartolomei, E. Carmona-Novillo, M. I. Hernández, J. Campos-Martínez, and R. Hernández-Lamoneda, J. Comput. Chem. (2010) (in press)]. The singlet and triplet PESs obtained from either RCCSD(T)-CASPT2 or RCCSD(T)-MRCI calculations are quite similar, although quantitative differences appear in specific terms of the expansion. CASPT2 calculations are the ones giving rise to larger splittings and more attractive interactions, particularly in the region of the absolute minima (in the rectangular D(2h) geometry). The new singlet, triplet, and quintet PESs are tested against second virial coefficient B(T) data and, their spherically averaged components, against integral cross sections measured with rotationally hot effusive beams. Both types of multiconfigurational approaches provide quite similar results, which, in turn, are in good agreement with the measurements. It is found that discrepancies with the experiments could be removed if the PESs were slightly more attractive. In this regard, the most attractive RCCSD(T)-CASPT2 PESs perform slightly better than the RCCSD(T)-MRCI counterpart.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Crossed beams study of the reaction 1CH2+C2H2-->C3H3+H

The reaction of electronically excited singlet methylene (1CH2) with acetylene (C2H2) was studied using the method of crossed molecular beams at a mean collision energy of 3.0 kcal/mol. The angular and velocity distributions of the propargyl radical (C3H3) products were measured using single photon ionization (9.6 eV) at the advanced light source. The measured distributions indicate that the mechanism involves formation of a long-lived C3H4 complex followed by simple C-H bond fission producing C3H3+H. This work, which is the first crossed beams study of a reaction involving an electronically excited polyatomic molecule, demonstrates the feasibility of crossed molecular beam studies of reactions involving 1CH2.     

Quasiclassical trajectory calculations of the reaction C+C2H2-->l-C3H, c-C3H+H, C3+H2 using full-dimensional triplet and singlet potential energy surfaces

Full-dimensional, density functional theory (B3LYP/6-311g(d,p))-based potential energy surfaces (PESs) are reported and used in quasi-classical calculations of the reaction of C with C(2)H(2). For the triplet case, the PES spans the region of the reactants, the complex region (with numerous minima and saddle points) and the products, linear(l)-C(3)H+H, cyclic(c)-C(3)H+H and c-(3)C(3)+H(2). For the singlet case, the PES describes the complex region and products l-C(3)H+H, c-C(3)H+H and l-(1)C(3)+H(2). The PESs are invariant under permutation of like nuclei and are fit to tens of thousands of electronic energies. Energies and harmonic frequencies of the PESs agree well the DFT ones for all stationary points and for the reactant and the products. Dynamics calculations on the triplet PES find both l-C(3)H and c-C(3)H products, with l-C(3)H being dominant at the energies considered. Limited unimolecular reaction dynamics on the singlet PES find both products in comparable amounts as well as the C(3)+H(2) product.     

Electron ionization of acetylene

Relative partial ionization cross sections and precursor specific relative partial ionization cross sections for fragment ions formed by electron ionization of C2H2 have been measured using time-of-flight mass spectrometry coupled with a 2D ion-ion coincidence technique. We report data for the formation of H+, H+2, C2+, C+/C2+ 2, CH+/C2H+2, CH+2, C+2, and C2H+ relative to the formation of C2H+2, as a function of ionizing electron energy from 30-200 eV. While excellent agreement is found between our data and one set of previously published absolute partial ionization cross sections, some discrepancies exist between the results presented here and two other recent determinations of these absolute partial ionization cross sections. We attribute these differences to the loss of some translationally energetic fragment ions in these earlier studies. Our relative precursor-specific partial ionization cross sections enable us, for the first time, to quantify the contribution to the yield of each fragment ion from single, double, and triple ionization. Analysis shows that at 50 eV double ionization contributes 2% to the total ion yield, increasing to over 10% at an ionizing energy of 100 eV. From our ion-ion coincidence data, we have derived branching ratios for charge separating dissociations of the acetylene dication. Comparison of our data to recent ab initio/RRKM calculations suggest that close to the double ionization potential C2H2+2 dissociates predominantly on the ground triplet potential energy surface (3Sigma*g) with a much smaller contribution from dissociation via the lowest singlet potential energy surface (1Delta g). Measurements of the kinetic energy released in the fragmentation reactions of C2H2+2 have been used to obtain precursor state energies for the formation of product ion pairs, and are shown to be in good agreement with available experimental data and with theory.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

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.     

Probing anisotropic interaction potentials of unsaturated hydrocarbons with He*(2 3S) metastable atom: attractive-site preference of sigma-direction in C2H2 and pi-direction in C2H4

State-resolved collision energy dependence of Penning ionization cross sections of acetylene (C2H2) and ethylene (C2H4) with He*(2 3S) metastable atoms was observed in a wide collision energy range from 20 to 350 meV. A recently developed discharge nozzle source with a liquid N2 circulator was employed for the measurements in the low-energy range from 20 to 80 meV. Based on classical trajectory calculations for the energy dependence of the partial ionization cross sections, anisotropic potential energy surfaces for the present systems were obtained by optimizing ab initio model potentials for the chemically related systems Li+C2H2 and C2H4. In the case of C2H2, the global minimum was found to be located around the H atom along the molecular axis with a well depth of 48 meV (ca. 1.1 kcal/mol). On the other hand, a dominant attractive well with a depth of 62 meV (ca. 1.4 kcal/mol) was found in the piCC electron region of C2H4. These findings were discussed in connection with orbital interactions between molecular orbitals of the target molecules and atomic orbitals of the metastable atom. It is concluded that sigma-type unoccupied molecular orbitals of C2H2 and a piCC-type highest occupied molecular orbital of C2H4 play a significant role for the attractive-site preference of sigma direction in C2H2 and pi direction in C2H4, respectively.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

Ab initio energies and product branching ratios for the O+C3H6 reaction

Intermediate and transition-state energies have been calculated for the O+C3H6 (propene) reaction using the compound ab initio CBS-QB3 and G3 methods in combination with density functional theory. The lowest-lying triplet and singlet potential energy surfaces of the O-C3H6 system were investigated. RRKM statistical theory was used to predict product branching fractions over the 300-3000 K temperature and 0.001-760 Torr pressure ranges. The oxygen atom adds to the C3H6 terminal olefinic carbon in the primary step to generate a nascent triplet biradical, CH3CHCH2O. On the triplet surface, unimolecular dissociation of CH3CHCH2O to yield H+CH3CHCHO is favored over the entire temperature range, although the competing H2CO+CH3CH product channel becomes significant at high temperature. Rearrangement of triplet CH3CHCH2O to CH3CH2CHO (propanal) via a 1,2 H-atom shift has a barrier of 122.3 kJ mol(-1), largely blocking this reaction channel and any subsequent dissociation products. Intersystem crossing of triplet CH3CHCH2O to the singlet surface, however, leads to facile rearrangement to singlet CH3CH2CHO, which dissociates via numerous product channels. Pressure was found to have little influence over the branching ratios under most conditions, suggesting that the vibrational self-relaxation rates for p相似文献   

Quantum dynamics study on the product branching for the C(3P) + C2H2 reaction: cyclic-C3H versus linear-C3H

Reduced-dimensionality quantum reactive scattering calculations for the C(3P) + C2H2 reaction have been carried out in order to understand the product branching dynamics of cyclic-C3H + H and linear-C3H + H. Our model treats only two degrees of freedom but can explicitly describe both of the C3H isomer product channels. The lowest triplet potential energy surface has been obtained by the hybrid density-functional method at the B3LYP/6-31G(d,p) level of theory. The calculated reaction probabilities were found to be dominated by resonance consistent with the complex-formation potential, and the results show that cyclic-C3H is preferentially formed via the cyclic-C3H2 intermediate produced by insertion of C(3P) into the CC bond. We have found that the isomerization from the cyclic-C3H2 to linear-C3H2 intermediate is suppressed by a barrier separating potential wells corresponding to these two intermediates. It has also been found that the energy dependence of the calculated total reaction cross section is in good agreement with the result of crossed molecular beam experiments.     

On the role of dynamical barriers in barrierless reactions at low energies: S(1D) + H2

Reaction probabilities as a function of total angular momentum (opacity functions) and the resulting reaction cross sections for the collision of open shell S((1)D) atoms with para-hydrogen have been calculated in the kinetic energy range 0.09-10 meV (1-120 K). The quantum mechanical hyperspherical reactive scattering method and quasi-classical trajectory and statistical quasi-classical trajectory approaches were used. Two different ab initio potential energy surfaces (PESs) have been considered. The widely used reproducing kernel Hilbert space (RKHS) PES by Ho et al. [T.-S. Ho, T. Hollebeek, H. Rabitz, S. D. Chao, R. T. Skodje, A. S. Zyubin, and A. M. Mebel, J. Chem. Phys 116, 4124 (2002)] and the recently published accurate double many-body expansion (DMBE)/complete basis set (CBS) PES by Song and Varandas [Y. Z. Song and A. J. C. Varandas, J. Chem. Phys. 130, 134317 (2009)]. The calculations at low collision energies reveal very different dynamical behaviors on the two PESs. The reactivity on the RKHS PES is found to be considerably larger than that on the DMBE/CBS PES as a result of larger reaction probabilities at low total (here also orbital) angular momentum values and to opacity functions which extend to significantly larger total angular momentum values. The observed differences have their origin in two major distinct topographic features. Although both PESs are essentially barrierless for equilibrium H-H distances, when the H-H bond is compressed the DMBE/CBS PES gives rise to a dynamical barrier which limits the reactivity of the system. This barrier is completely absent in the RHKS PES. In addition, the latter PES exhibits a van der Walls well in the entrance channel which reduces the height of the centrifugal barrier and is able to support resonances. As a result, a significant larger cross section is found on this PES, with marked oscillations attributable to shape resonances and/or to the opening of partial wave contributions. The comparison of the results on both PESs is illustrative of the wealth of the dynamics at low collision energy. It is also illuminating about the difficulties encountered in modeling an all-purpose global potential energy surface.     

Exploring the accuracy level of new potential energy surfaces for the F + HD reactions: from exact quantum rate constants to the state-to-state reaction dynamics

Exact quantum reactive scattering calculations in the collision energy range 1-250 meV have been carried out for both the isotopic product channels of the title system. The dynamical studies compares an ab initio potential energy surface (PES) recently appeared in the literature (J. Chem. Phys., 2008, 129, 011103) with other phenomenological PESs. Vibrational branching ratios, cross sections and rate constants are presented and compared with molecular beam scattering experiments as well as with chemical kinetics data. In particular, the agreement of the vibrational branching ratios with experimental measurements is improved with respect to previous studies on other PESs, mainly because of the presence of a broad peak in the HF(v' = 3) integral cross section completely absent in the previous simulations. This feature, observed by molecular beam experiments, is the fingerprint of a new reaction mechanism operative in the dynamics described by the new PES. A conjecture for its origin, able to explain many of its characteristic aspects, is analyzed and discussed.