The F + HD(v = 0, 1; j = 0, 1) reactions: stereodynamical properties of orbiting resonances.

The excitation functions (reaction cross-section as a function of collision energy) of the F + HD(v = 0, 1; j = 0, 1) benchmark system have been calculated in the 0.01-6 meV collision energy interval using a time-independent hyperspherical quantum dynamics methodology. Special attention has been paid to orbiting resonances, which bring about detailed information on the three-atom interaction during the reactive encounter. The location of the resonances depends on the rovibrational state of the reactants HD(v,j), but is the same for the two product channels HF + D and DF + H, as expected for these resonances that are linked to the van der Waals well at the entrance. The resonance intensities depend both on the entrance and on the exit channels. The peak intensities for the HF + D channel are systematically larger than those for DF + H. Vibrational excitation leads to an increase of the peak intensity by more than an order of magnitude, but rotational excitation has a less drastic effect. It deceases the resonance intensity of the F + HD(v = 1) reaction, but increases somewhat that of F + HD(v = 0). Polarization of the rotational angular momentum with respect to the initial velocity reveals intrinsic directional preferences in the F + HD(v = 0, 1; j = 1) reactions that are manifested in the resonance patterns. The helicities (Ω = 0, Ω = ±1) possible for j = 1 contribute to the resonances, but that from Ω± 1 is, in general, dominant and in some cases exclusive. It corresponds to a preferential alignment of the HD internuclear axis perpendicular to the initial direction of approach and, thus, to side-on collisions. This work also shows that external preparation of the reactants, following the intrinsic preferences, would allow the enhancement or reduction of specific resonance features, and would be of great help for their eventual experimental detection.

Competing Dynamical Mechanisms in Inelastic Collisions of H + HF.

The dynamics of inelastic collisions between HF and H has been investigated in detail by means of time-independent quantum mechanical calculations on the LWA-78 potential energy surface ( Li , G. ; et al. J. Chem. Phys. 2007 , 127 , 174302 ). Reaction probabilities, differential cross sections, and three-vector correlations have been calculated and analyzed. Our results show that there are two competing collision mechanisms that correlate with low and high impact parameters and show very different stereodynamical preferences. The mechanism promoted by high impact parameters is the only one present at low collision energies. We also observe the presence of an apparent threshold in the inelastic cross section for relatively high initial HF rotational quantum numbers, which is associated with the larger energy difference between adjacent rotational quantum states with increasing rotation.

Orbiting resonances in the F + HD (v = 0, 1) reaction at very low collision energies. A quantum dynamical study.

Time-independent, fully converged, quantum dynamical calculations have been performed for the F + HD (v = 0, j = 0) and F + HD (v = 1, j = 0) reactions on an accurate potential energy surface down to collision energies of 0.01 meV. The two isotopic exit channels, HF + D and DF + H, have been investigated. The calculations reproduce satisfactorily the Feshbach resonance structures for collision energies between 10 and 40 meV, previously reported in the literature for the HF + D channel. Contrary to the results of a former literature work, vibrational excitation of HD is found to enhance reactivity in all cases down to the lowest collision energy investigated. Shape-type orbiting resonances are found for collision energies lower than 2 meV. The resonances appear as peaks in the reaction cross sections that are associated to specific values of the total angular momentum, J. In contrast with the Feshbach resonances at higher energies, the orbiting resonance structure, which is caused by the van der Waals well of the entrance channel, is identical for the HF + D and DF + H exit channels. The orbiting resonance peaks for F + HD (v = 0) are very small, but those for F + HD (v = 1) could be observed, in principle, with a combination of Raman pumping and merged beams methods.

Influence of vibration in the reactive scattering of D + MuH: the effect of dynamical bonding.

The dynamics of the D + MuH(v = 1) reaction has been investigated using time-independent quantum mechanical calculations. The total reaction cross sections and rate coefficients have been calculated for the two exit channels of the reaction leading, respectively, to DMu + H and DH + Mu. Over the 100-1000 K temperature range investigated the rate coefficients for the DMu + H channel are of the order of 10(-10) cm(3) s(-1) and those for the DH + Mu channel vary between 1 × 10(-12) and 8 × 10(-11) cm(3) s(-1). These results point to a virtually barrierless reaction for the DMu + H channel and to the presence of a comparatively small barrier for the DH + Mu channel and are consistent with the profiles of their respective collinear vibrationally adiabatic potentials (VAPs). The effective barrier in the VAP of the DH + Mu channel is located in the reactant valley and, consequently, translation is found to be more efficient than vibration for the promotion of the reaction over a large energy interval in the post threshold region. Below this barrier, the DH + Mu channel can be accessible through an indirect mechanism implying crossing from the DMu + H pathway. The most salient feature found in the present study is revealed in the total reaction cross section for the DMu + H channel, which shows a sharp resonance caused by the presence of a deep well in the vibrationally adiabatic potential. This well has a dynamical origin, reminiscent of that found recently in the vibrationally bonded BrMuBr complex [Fleming, et al., Angew. Chem., Int. Ed., 2014, 53, 1], and is due to the stabilizing effect of the light Mu atom oscillating between the heavier H and D isotopes and to the bond softening associated with vibrational excitation of MuH.

A semiclassical treatment of the ℓ-j correlation in atom-diatom collisions.

The explicit consideration of the vector correlations is an essential step when it comes to determining the mechanism of chemical reactions. Usual vector correlations involve initial and final relative velocity vectors and rotational angular momenta. However, the correlation between the orbital, ℓ, and rotational, j, angular momenta has seldom received any attention. In this article, we present a semiclassical methodology capable of describing the ℓ-j correlation that may serve as a connection between the quantum and quasiclassical treatments. Using the scattering matrix in the orbital angular momentum representation, the ℓ-j correlation is expressed as a probability density function of the angle formed by both vectors. This technique is exemplified through the H + D2 reaction and its accuracy is appraised by comparing with results derived from calculations based on quasiclassical trajectories.

Comparative dynamics of the two channels of the reaction of D + MuH.

The dynamics of the asymmetric D + MuH (Mu = Muonium) reaction leading to Mu exchange, DMu + H, and H abstraction, DH + Mu, channels has been investigated using time-independent quantum mechanical (QM) calculations. Reaction probabilities, cross sections, cumulative reaction probabilities, and rate coefficients were determined for the two exit channels of the reaction. Quasiclassical trajectory (QCT) calculations were also performed in order to check the reliability of the method for this reaction and to discern the genuine quantum effects. Overall, the Mu exchange channel exhibits more structured reaction probabilities and cross sections with much larger rate coefficients than the H abstraction counterpart. Over the 100-1000 K temperature interval considered in this study, the QM rate coefficients for the Mu exchange vary between ≈5 × 10(-15) and 2 × 10(-11) cm(3) s(-1) and those for the generation of DH + Mu between 2 × 10(-18) and 3.5 × 10(-12) cm(3) s(-1). In common with the rest of the isotopologues of the H + H2 system, the height of the respective barriers in the collinear (symmetric stretch) vibrationally adiabatic potential energy curves matches the classical total energy threshold very accurately. Indeed, the lower and narrower vibrationally adiabatic collinear barrier as compared with that for the DH + Mu formation determines the preponderance of the DMu + H channel. Comparison of QM and QCT results and their analysis show that tunneling accounts for the reactivity at energies below the height of these barriers and that its effect on the rate coefficients becomes appreciable below 300 K. As expected, with growing temperature the contribution of tunneling to the global reactivity decreases markedly, but the rate coefficients are still much higher for the Mu exchange channel due to the effect of MuH rotational excitation that boosts the formation of DMu while diminishing the H abstraction channel that leads to DH formation. The analysis of the thermal cumulative reaction probabilities of the two channels indicates that at the lowest energies/temperatures the reaction into the DH + Mu channel takes place via'leakage' from collisions proceeding along the DMu + H reaction path.

Dynamics of the reactions of muonium and deuterium atoms with vibrationally excited hydrogen molecules: tunneling and vibrational adiabaticity.

Quantum mechanical (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the exchange reactions of D and Mu (Mu = muonium) with hydrogen molecules in their ground and first vibrational states. In all the cases considered, the QM rate coefficients, k(T), are in very good agreement with the available experimental results. In particular, QM calculations on the most accurate potential energy surfaces (PESs) predict a rate coefficient for the Mu + H(2) (ν = 1) reaction which is very close to the preliminary estimate of its experimental value at 300 K. In contrast to the D + H(2) (ν = 0,1) and the Mu + H(2) (ν = 0) reactions, the QCT calculations for Mu + H(2) (ν = 1) predict a much smaller k(T) than that obtained with the accurate QM method. This behaviour is indicative of tunneling. The QM reaction probabilities and total reactive cross sections show that the total energy thresholds for the reactions of Mu with H(2) in ν = 0 and ν = 1 are very similar, whereas for the corresponding reaction with D the ν = 0 total energy threshold is about 0.3 eV lower than that for ν = 1. The results just mentioned can be explained by considering the vibrational adiabatic potentials along the minimum energy path. The threshold for the reaction of Mu with H(2) in both ν = 0 and ν = 1 states is the same and is given by the height of the ground vibrational adiabatic collinear potential, whereas for the D + H(2) reaction the adiabaticity is preserved and the threshold for the reaction in ν = 1 is very close to the height of the ν = 1 adiabatic collinear barrier. For Mu + H(2) (ν = 1) the reaction takes place by crossing from the ν = 1 to the ν = 0 adiabat, since the exit channel leading to MuH (ν = 1) is not energetically accessible. At the lowest possible energies, the non-adiabatic vibrational crossing implies a strong tunneling effect through the ν = 1 adiabatic barrier. Absence of tunneling in the classical calculations results in a threshold that coincides with the height of the ν = 1 adiabatic barrier. Most interestingly, the expected tunneling effect in the reaction of Mu with hydrogen molecules occurs for H(2) (ν = 1) but not for H(2) (ν = 0) where zero-point-energy effects clearly dominate.

Three-vector correlation in statistical reactions: the role of the triatomic parity.

This article presents a methodology for the determination of the k-j-k' three-vector correlation assuming a statistical model for atom-diatom reactions; k and k' are the reagent-approach and product-recoil directions, respectively, and j is the rotational angular momentum of the reagent diatomic. Although the polarization of reagent angular momentum is in most cases negligible, conservation of the triatomic parity imposes a certain polarization for some combinations involving low reagent and product rotational states. Statistical and quantum-mechanical polarization-dependent differential cross sections were calculated for the barrierless D(+) + H(2)(v = 0,j) → HD(v' = 0,j') + H(+) reaction. The agreement between the two is in most cases excellent, confirming the statistical character of the reaction at low and moderate collision energies.

Dynamics of the D+ + H2 and H+ + D2 reactions: a detailed comparison between theory and experiment.

An extensive set of experimental measurements on the dynamics of the H(+) + D(2) and D(+) + H(2) ion-molecule reactions is compared with the results of quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical quasiclassical trajectory (SQCT) calculations. The dynamical observables considered include specific rate coefficients as a function of the translational energy, E(T), thermal rate coefficients in the 100-500 K temperature range. In addition, kinetic energy spectra (KES) of the D(+) ions reactively scattered in H(+) + D(2) collisions are also presented for translational energies between 0.4 eV and 2.0 eV. For the two reactions, the best global agreement between experiment and theory over the whole energy range corresponds to the QCT calculations using a gaussian binning (GB) procedure, which gives more weight to trajectories whose product vibrational action is closer to the actual integer QM values. The QM calculations also perform well, although somewhat worse over the more limited range of translational energies where they are available (E(T) < 0.6 eV and E(T) < 0.2 eV for the H(+) + D(2) and D(+) + H(2) reactions, respectively). The worst agreement is obtained with the SQCT method, which is only adequate for low translational energies. The comparison between theory and experiment also suggests that the most reliable rate coefficient measurements are those obtained with the merged beams technique. It is worth noting that none of the theoretical approaches can account satisfactorily for the experimental specific rate coefficients of H(+) + D(2) for E(T)≤ 0.2 eV although there is a considerable scatter in the existing measurements. On the whole, the best agreement with the experimental laboratory KES is obtained with the simulations carried out using the state resolved differential cross sections (DCSs) calculated with the QCT-GB method, which seems to account for most of the observed features. In contrast, the simulations with the SQCT data predict kinetic energy spectra (KES) considerably cooler than those experimentally determined.

Orientation effects in Cl + H2 inelastic collisions: characterization of the mechanisms.

Based on quantum mechanical scattering (QM) calculations, we have analyzed the polarization of the product hydrogen molecule in Cl + H(2) (v = 0, j = 0) inelastic collisions. The spatial arrangements adopted by the rotational angular momentum and internuclear axis of the departing molecule have been characterized and used to prove that two distinct mechanisms, corresponding to different dynamical regimes, are responsible for the inelastic collisions. Such mechanisms, named as low-b and high-b, correlate with well defined ranges of impact parameter values, add in an essentially incoherent way, and can be clearly differentiated through the quantum mechanical polarization moment that measures the orientation of the products rotational angular momentum with respect to the scattering plane. Other directional effects turn out to fail when it comes to distinguishing the mechanisms. Quasiclassical trajectories (QCT) calculations have been used as a supplement to the purely quantum mechanical analysis. By combining QM and QCT results, which are in very good agreement, we have succeeded in obtaining a clear and meaningful picture of how the two types of collisions take place.

Cumulative reaction probabilities and transition state properties: a study of the H+ + H2 and H+ + D2 proton exchange reactions.

Cumulative reaction probabilities (CRPs) have been calculated by accurate (converged, close coupling) quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical QCT (SQCT) methods for the H(+) + H(2) and H(+) + D(2) reactions at collision energies up to 1.2 eV and total angular momentum J = 0-4. A marked resonance structure is found in the QM CRP, most especially for the H(3)(+) system and J = 0. When the CRPs are resolved in their ortho and para contributions, a clear steplike structure is found associated with the opening of internal states of reactants and products. The comparison of the QCT results with those of the other methods evinces the occurrence of two transition states, one at the entrance and one at the exit. At low J values, except for the quantal resonance structure and the lack of quantization in the product channel, the agreement between QM and QCT is very good. The SQCT model, that reflects the steplike structure associated with the opening of initial and final states accurately, clearly tends to overestimate the value of the CRP as the collision energy increases. This effect seems more marked for the H(+) + D(2) isotopic variant. For sufficiently high J values, the growth of the centrifugal barrier leads to an increase in the threshold of the CRP. At these high J values the discrepancy between SQCT and QCT becomes larger and is magnified with growing collision energy. The total CRPs calculated with the QCT and SQCT methods allowed the determination of the rate constant for the H(+) + D(2) reaction. It was found that the rate, in agreement with experiment, decreases with temperature as expected for an endothermic reaction. In the range of temperatures between 200 and 500 K the differences between SQCT and QCT rate results are relatively minor. Although exact QM calculations are formidable for an exact determination of the k(T), it can be reliably expected that their value will lie between those given by the dynamical and statistical trajectory methods.

Vibrationally inelastic collisions of H+D2: a comparison of quantum mechanical, quasiclassical, and experimental results.

A detailed comparison of quantum mechanical (QM) and quasiclassical trajectory (QCT) integral and differential cross sections (DCSs) as well as opacity functions is presented in this work for the vibrationally inelastic collisions of H+D(2)(v=0,j=0)-->H+D(2)(v(')=3,j(')) at 1.72 eV collision energy. These results are also compared with the experimental differential cross sections by Greaves et al. [Nature (London) 454, 88 (2008)]. The agreement between QCT and QM results is fairly good but some differences are appreciable, and it is shown that the experimental results are in a somewhat better agreement with the calculated QM DCS. The present results and their analysis confirm that the vibrational excitation takes place by elongation of the D-D bond in a "tug-of-war" mechanism, where the incoming H atom and one of the D atoms compete for the formation of a bond with the other D atom, as proposed by Greaves et al. It is also found that these collisions may give rise to the formation of short-lived collision complexes (tau(coll)=35-50 fs) that can be traced back to the presence of relatively deep wells in the potential surface when the original D-D bond is stretched. The analysis of the trajectories into v(')=3 reveals that most of them cross at least twice the reaction barrier via a recrossing mechanism.

A comparison of quantum and quasiclassical statistical models for reactions of electronically excited atoms with molecular hydrogen.

A detailed comparison of statistical models based on the quasiclassical trajectory (SQCT) and quantum mechanical (SQM) methods is presented in this work for the C((1)D)+H(2), S((1)D)+H(2), O((1)D)+H(2) and N((2)D)+H(2) insertion reactions. Reaction probabilities, integral (ICS) and differential (DCS) cross sections at different levels of product's state resolution are shown and discussed for these reactions. The agreement is in most cases excellent and indicates that the effect of tunneling through the centrifugal barrier is negligible. However, if there exists a dynamical barrier, as in the case of the N((2)D)+H(2) reaction, some of the SQM results can be slightly different than those calculated with the SQCT model. The rationale of the observed similarities and discrepancies can be traced back to the specific topologies of the potential energy surfaces for each of the reactions examined. The SQCT model is sensitive enough to show the relatively small inaccuracies resulting from the decoupling inherent to the centrifugal sudden approximation when used in the SQM calculations. In addition, the effect of ignoring the parity conservation is also examined. This effect is in general minor except in particular cases such as the DCS from initial rotational state j=0, which requires, in order to reproduce the sharp forward and backward peaks, the explicit conservation of parity.

Cumulative reaction probabilities and transition state properties: a study of the F + H2 reaction and its deuterated isotopic variants.

A comparative quantum mechanical (QM) and quasiclassical trajectory (QCT) study of the cumulative reaction probabilities (CRPs) is presented in this work for the F + H(2) reaction and its isotopic variants for low values of the total angular momentum J. The agreement between the two sets of calculations is very good with the exception of some features whose origin is genuinely QM. The agreement also extends to the CRP resolved in the helicity quantum number k. The most remarkable feature is the steplike structure, which becomes clearly distinct when the CRPs are resolved in odd and even rotational states j. The analysis of these steps shows that each successive increment is due to the opening of the consecutive rovibrational states of the H(2) or D(2) molecule, which, in this case, nearly coincide with those of the transition state. Moreover, the height of each step reflects the number of helicity states compatible with a given J and j values, thus indicating that the various helicity states for a specific j have basically the same contribution to the CRPs at a given total energy. As a consequence, the dependence with k of the reactivity is practically negligible, suggesting very small steric restrictions for any possible orientation of the reactants. This behavior is in marked contrast to that found in the D + H(2) reaction, wherein a strong k dependence was found in the threshold and magnitude of the CRP. The advantages of a combined QCT and QM approaches to the study of CRPs are emphasized in this work.

Stringent test of the statistical quasiclassical trajectory model for the H+3 exchange reaction: a comparison with rigorous statistical quantum mechanical results.

A complete formulation of a statistical quasiclassical trajectory (SQCT) model is presented in this work along with a detailed comparison with results obtained with the statistical quantum mechanical (SQM) model for the H+ +D2 and H+ +H2 reactions. The basic difference between the SQCT and the SQM models lies in the fact that trajectories instead of wave functions are propagated in the entrance and exit channels. Other than this the two formulations are entirely similar and both comply with the principle of detailed balance and conservation of parity. Reaction probabilities, and integral and differential cross sections (DCS's) for these reactions at different levels of product's state resolution and from various initial states are shown and discussed. The agreement is in most cases excellent and indicates that the effect of tunneling through the centrifugal barrier is negligible. Some differences are found, however, between state resolved observables calculated by the SQCT and the SQM methods which makes use of the centrifugal sudden (coupled states) approximation (SQM-CS). When this approximation is removed and the full close coupling treatment is used in the SQM model (SQM-CC), an almost perfect agreement is achieved. This shows that the SQCT is sensitive enough to show the relatively small inaccuracies resulting from the decoupling inherent to the CS approximation. In addition, the effect of ignoring the parity conservation is thoroughly examined. This effect is in general minor except in particular cases such as the DCS from initial rotational state j=0. It is shown, however, that in order to reproduce the sharp forward and backward peaks the conservation of parity has to be taken into account.

Constraints at the transition state of the D + H2 reaction: quantum bottlenecks vs. stereodynamics.

This article presents a quasiclassical trajectory method for the calculation of cumulative reaction probabilities by sampling of the helicity quantum number of the reagents (k). The method is applied to the D + H(2) reaction at various total angular momentum (J) values, and the helicity-resolved quasiclassical cumulative reaction probabilities are compared to their quantum mechanical counterparts. The agreement between the two sets of results is fairly good. In particular, k-dependent, J-independent reaction thresholds found with quantum methods are reproduced by the quasiclassical calculations. The shift of these thresholds with increasing k, which has been previously attributed to the quantum bottleneck states taking part in the reaction, is revisited and discussed also in terms of the reaction stereodynamics.

A statistical quasiclassical trajectory model for atom-diatom insertion reactions.

A statistical model based on the quasiclassical trajectory method is presented in this work for atom-diatom insertion reactions. The basic difference between this and the corresponding statistical quantum model (SQM) lies in the fact that trajectories instead of wave functions are propagated in the entrance and exit channels. Other than this the two formulations are entirely similar. In particular, it is shown that conservation of parity can be taken into account in a natural and precise way in the statistical quasiclassical trajectory (SQCT) model. Additionally, the SQCT model complies with the principle of detailed balance and overcomes the problem of the zero point energy in the products. As a test, the model is applied to the H3+ and H+D2 exchange reactions. The excellent agreement between the SQCT and SQM results, especially in the case of the differential cross sections, indicates that the effect of tunneling through the centrifugal barrier is negligible. The effect of ignoring quantum mechanical parity conservation is also investigated.

Analysis of the H + D2 reaction mechanism through consideration of the intrinsic reactant polarisation.

The effect of reactant polarisation on the dynamics of the title reaction at collision energies up to 1.6 eV is analysed in depth. The analysis takes advantage of two novel theoretical concepts: intrinsic reaction properties and stereodynamical portraits. Exact quantum methods are used to determine the polarisation moments that quantify the intrinsic reactant polarisation at various levels of detail, including or not product state and/or scattering angle resolution. The data is then examined with the aid of stereodynamical portraits, which facilitate the rationalisation of the stereochemical effects that are relevant for the reaction dynamics. This allows for detailed characterisations of the so-called direct and delayed reaction mechanisms.

Mechanism and control of the F+H2 reaction at low and ultralow collision energies.

This article uses theoretical methods to study the dependence on stereodynamical factors of the mechanism and reactivity of the F+H2 reaction at low and ultralow collision energies. The impact of polarization of the H2 reactant on total and state-to-state integral and differential cross sections is analyzed. This leads to detailed pictures of the reaction mechanism in the cold and ultracold regimes, accounting, in particular, for distinctions associated with the various product states and scattering angles. The extent to which selection of reactant polarization allows for external control of the reactivity and reaction mechanism is assessed. This reveals that even the simplest of reactant polarization schemes allows for fine, product state-selective control of differential and (for reactions involving more than a single, zero orbital angular momentum partial wave) integral cross sections.

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results.

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.