On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.     

Sticking and desorption of hydrogen on graphite: a comparative study of different models

We study the physisorption of atomic hydrogen on graphitic surfaces with four different quantum mechanical methods: perturbation and effective Hamiltonian theories, close coupling wavepacket, and reduced density matrix propagation methods. Corrugation is included in the modeling of the surface. Sticking is a fast process which is well described by all methods. Sticking probabilities are of the order of a few percent in the collision energy range 0-25 meV, but are enhanced for collision energies close to those of diffraction resonances. Sticking also increases with surface temperature. Desorption is a slow process which involves multiphonon processes. We show, however, how to correct the close coupling wavepacket method to account for such phenomena and obtain correct time constants for initial state decay. Desorption time constants are in the range of 20-50 ps for a surface temperature of 300 K.     

Cross sections and low temperature rate coefficients for the H + CH+ reaction: a quasiclassical trajectory study

The H + CH(+) reaction is studied by quasiclassical trajectory (QCT) calculations, along with phase space theory (PST) and quantum rigid rotor calculations, employing a global single-valued potential energy surface recently derived by our group. We report QCT total cross sections for each of the three channels, for low collision energies and different reactant rotational quantum numbers. At the lowest collision energies, all cross sections exhibit a capture-like behaviour, as expected from a barrierless reaction. At higher energies, there are important dynamical effects coming from the opening of new channels in the inelastic and reactive exchange collisions. The inelastic cross sections turn out to largely increase, while the reactive abstraction cross sections are declining faster than predicted by the capture theory. A large value of the reactant rotational quantum number tends to suppress these dynamical effects. The QCT rate coefficients are reported for a temperature range from 1-700 K. Below 20 K, the abstraction and exchange QCT rate coefficients are almost constant, as predicted by the capture theory. Above this temperature, the abstraction rate coefficient declines, while the exchange and inelastic rate coefficients are increasing, due to the opening of new channels. A good agreement is observed between the experimental abstraction rate coefficient and the QCT and PST ones. The QCT inelastic results are also compared with those obtained from rigid rotor close coupling (CCRR) calculations in order to check the ability of this approach to provide a reliable estimate of the inelastic rate coefficients for a reactive system without a barrier. The laws of variation as a function of temperature are found to be very similar and the curves are parallel above 20 K. However, reaction is not allowed in the rigid rotor approximation, therefore the CCRR results are about twice as large as their QCT counterparts.     

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.     

Quasiclassical trajectory calculations to evaluate a kinematic constraint on internal energy in suprathreshold collision energy abstraction reactions

Experimentally observed product quantum state distributions across a wide range of abstraction reactions at suprathreshold collision energies have shown a strong bias against product internal energy. Only a fraction, sometimes quite a small fraction, of the energetically accessible product quantum states are populated. Picconatto et al. [J. Chem. Phys. 114, 1663 (2001)] noted a simple mathematical relationship between the highest-energy rovibrational states observed and the kinematics of the reaction system. They proposed a reaction model based on reaction kinematics that quantitatively explains this behavior. The model is in excellent agreement with measured quantum state distributions. The assumptions of the model invoke detailed characteristics of reactive trajectories at suprathreshold collision energies. Here we test those assumptions using quasiclassical trajectory calculations for the abstraction reactions H+HCl-->H2+Cl, D+HCl-->HD+Cl, and H+DCl-->HD+Cl. Trajectories were run on a potential-energy surface calculated with a London-Eyring-Polyani-Sato function with a localized 3-center term (LEPS-3C) previously shown to accurately reproduce experimentally observed product state distributions for the H+HCl abstraction reaction. The trajectories sample collision energies near threshold and also substantially above it. Although the trajectories demonstrate some aspects of the model, they show that it is not valid. However, the inadequacy of the proposed model does not invalidate the apparent kinematic basis of the observed energy constraint. The present results show that there must be some other molecular behavior rooted in the reaction kinematics that is the explanation and the source of the constraint.     

Quantum dynamics study of Li + HF reaction

We report rigorous quantum dynamics studies of the Li + HF reaction using the time-dependent wavepacket approach. The dynamics study is carried out on a recent ab initio potential energy surface, and state-selected reaction probabilities and cross sections are calculated up to 0.4 eV of collision energy. Many long-lived resonances (as long as 10 ps) at low collision energies (below 0.1 eV) are uncovered from the dynamics calculation. These long-lived resonances play a dominant role in the title reaction at low collision energies (below 0.1 eV). At higher energies, the direct reaction process becomes very important. The reaction probabilities from even rotational states exhibit a different energy dependence than those from odd rotational states. Our calculated integral cross section exhibits a broad maximum near the collision energy of 0.26 eV with small oscillations superimposed on the broad envelope which is reminiscent of the underlying resonance structures in reaction probabilities. The energy dependence of the present CS cross section is qualitatively different from the simple J-shifting approximation, in which a monotonic increase of cross section with collision energy was obtained. Received: 8 January 1997 / Accepted: 14 January 1997     

Methane dissociation on Ni(111): the effects of lattice motion and relaxation on reactivity

The effects of lattice motion and relaxation on the dissociative adsorption of methane on a Ni(111) surface are explored. Electronic structure methods based on the density functional theory are used to compute the potential energy surface for this reaction. It is found that, in the transition state and product regions, there are forces causing the Ni atom over which the molecule dissociates to move out of the surface. In order to examine the extent to which the lattice might pucker during this reaction, high dimensional fully quantum scattering calculations are carried out. It is found that a significant amount of lattice puckering can occur, even at large collision energies, lowering the barrier to reaction and increasing the dissociative sticking probability. This is shown to be in contrast to the predictions of the surface oscillator model. While we observe similar puckering forces for this reaction on Pt(111), our calculations suggest that the puckering on this surface will be considerably less due to the larger metal atom mass. The "laser off" reactivities of CD(3)H on Ni(111) are computed, and it is demonstrated that there can be significant contributions to the reactivity from vibrationally excited molecules, particularly at lower collision energies, or when a large nozzle temperature is required to attain the necessary collision energy for reaction. Comparisons are made with recent experiments with regard to the variation of reactivity with collision energy, vibrational state, and surface temperature.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.     

Quantum and quasiclassical state-to-state dynamics of the NH + H reaction: competition between abstraction and exchange channels

Quantum and quasiclassical state-to-state dynamics for the NH + H' reaction at high collision energies up to 1.6 eV was studied on an accurate ab initio potential energy surface. Both of the endothermic abstraction (NH + H' → N + HH') and thermoneutral exchange (NH + H' → H + NH') channels were investigated from the same set of wave packets using an efficient coordinate transformation method. It is found that the abstraction represents a minor reaction channel in the energy range studied, primarily due to endothermicity. The cross section for the abstraction reaction increases monotonically with the collision energy, while that for the exchange reaction is relatively energy insensitive. As a result, the thermal rate constant for the abstraction reaction follows the Arrhenius law, where that for the exchange reaction is nearly temperature independent. Finally, it is shown that the quantum mechanical results can be reasonably reproduced by the Gaussian-binning quasiclassical trajectory method and to a lesser extent by a quantum statistical model.     

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.     

用含时波包方法将SVRT模型应用到F+CD4→CD3+DF反应中

In order to study the quantum reaction dynamics of large molecular systems, the time dependent quantum wave packet approach was used to study the F+CD4→CD3+DF reaction systems. The semirigid vibrating rotor model proposed by J.Z.H. Zhang was used on the MJ1 potential energy surface. The barrier height of the MJ1 PES was about 66 meV. In the semirigid vibrating rotor model, the fragment CD3 was fixed to in the geometry, its transition state value, because from the reactant to the transition state the C-D bond in the CD3 group almost remains constant, which can be treated as a spectator bond. The numerical calculation showed that there were oscillatory structures in the energy dependence of the calculated integral cross section. Those structures are generally associated with dynamic resonances. Cross section and rate constant were calculated based on the MJ1 PES of the ground state. These results are comparable to the results of previous calculations and reaction dynamic experiment results. At low temperature and collision energies, the tunneling effect works most remarkably in the reaction process to make the D abstraction easier. At high temperatures and collision energies, the rate constant is higher than the experimental results.     

Bimolecular reactions of vibrationally excited molecules. Roaming atom mechanism at low kinetic energies

Quasiclassical trajectory calculations have been performed for the H + H'X(v) → X + HH' abstraction and H + H'X(v) → XH + H' (X = Cl, F) exchange reactions of the vibrationally excited diatomic reactant at a wide collision energy range extending to ultracold temperatures. Vibrational excitation of the reactant increases the abstraction cross sections significantly. If the vibrational excitation is larger than the height of the potential barrier for reaction, the reactive cross sections diverge at very low collision energies, similarly to capture reactions. The divergence is quenched by rotational excitation but returns if the reactant rotates fast. The thermal rate coefficients for vibrationally excited reactants are very large, approach or exceed the gas kinetic limit because of the capture-type divergence at low collision energies. The Arrhenius activation energies assume small negative values at and below room temperature, if the vibrational quantum number is larger than 1 for HCl and larger than 3 for HF. The exchange reaction also exhibits capture-type divergence, but the rate coefficients are larger. Comparisons are presented between classical and quantum mechanical results at low collision energies. At low collision energies the importance of the exchange reaction is enhanced by a roaming atom mechanism, namely, collisions leading to H atom exchange but bypassing the exchange barrier. Such collisions probably have a large role under ultracold conditions. The calculations indicate that for roaming to occur, long-range attractive interaction and small relative kinetic energy in the chemical reaction at the first encounter are necessary, which ensures that the partners can not leave the attractive well. Large orbital angular momentum of the primary products (equivalent to large rotational excitation in a unimolecular reaction) is favorable for roaming.     

Quantum dynamics of Br + HD reaction

In this article we report the results of three-dimensional time-dependent quantum wavepacket calculations carried out for the Br + HD( v = 0, j = 0) reaction in the collision energy range 0.0-1.2 eV. An accurate potential energy surface computed by Kurosaki was used for the dynamical calculations. Both reactive channels, BrH + D and BrD + H, show vibrational enhancement of the reaction cross sections. For the three initial vibrational states considered, the production of BrD channel dominates over that of BrH for the considered collision energy range. The two arrangement channels exhibit different initial rotational state dependence. The cross section for the formation of BrD is almost independent of j whereas the same for the formation of BrH increases with increase in j. A comparison with the results on an e-LEPS surface shows that the two surfaces behave very differently with respect to the cross section for the initial rotational states.     

Rotational excitation cross sections of para-H2 + para-H2 collisions. A full-dimensional wave-packet propagation study using an exact form of the kinetic energy

A full-dimensional quantum dynamical study of the rotational excitation in para-para H2 + H2 collisions using the potential-energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] is reported. The multiconfiguration time-dependent Hartree algorithm is utilized to propagate wave packets and the cross sections for collision energies up to 1.2 eV are determined by a flux analysis through the interaction of the wave packet with a complex absorbing potential. Calculations for a collection of total angular momenta up to J = 70 are performed; the missing channels are obtained with a J-interpolation algorithm.     

QCT and QM calculations of the Cl(2P) + NH3 reaction: influence of the reactant well on the dynamics

A detailed dynamics study, using both quasi-classical trajectory (QCT) and reduced-dimensional quantum mechanical (QM) calculations, was carried out to understand the reactivity and mechanism of the Cl((2)P) + NH(3)→ HCl + NH(2) gas-phase reaction, which evolves through deep wells in the entry and exit channels. The calculations were performed on an analytical potential energy surface recently developed by our group, PES-2010 [M. Monge-Palacios, C. Rangel, J. C. Corchado and J. Espinosa-Garcia, Int. J. Quantum. Chem., 2011], together with a simplified model surface, mod-PES, in which the reactant well is removed to analyze its influence. The main finding was that the QCT and QM methods show a change of the reaction probability with collision energy, suggesting a change of the atomic-level mechanism of reaction with energy. This change disappeared when the mod-PES was used, showing that the behaviour at low energies is a direct consequence of the existence of the reactant well. Analysis of the trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) > 5 kcal mol(-1)) practically all trajectories are direct, at low energies (E(coll) < 3 kcal mol(-1)) the trajectories are indirect, i.e., with the mediation of a trapping complex in the entry and/or the exit wells. The reactant complex allows repeated encounters between the reactants, increasing the reaction probability at low energies. The differential cross section results reinforce this change of mechanism, showing also the influence of the reactant well on this reaction. Thus, the PES-2010 surface yields a forward-backward symmetry in the scattering, while when the reactant well is removed with the mod-PES the shape is more isotropic.     

Classical theory of collisional entropy production

A classical dynamical theory of elementary collision processes is formulated in analogy to the quantum theory of the dynamical scattering matrix, which can be defined for a pure quantum stationary scattering state. The elements of this matrix are probability amplitudes for transitions between internal states defined for given values of a reaction coordinate. The squared magnitudes of these amplitudes, modeled in the proposed classical theory, define normalized internal state population distributions suitable for information theoretical analysis. Statistical entropy and surprisal are defined as dynamical functions of a reaction coordinate. This formalism differs fundamentally from concepts based on the classical Liouville equation.     

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.     

State-to-state quantum dynamics of the H + HBr reaction: competition between the abstraction and exchange reactions

Quantum state-to-state dynamics for the H + HBr(υ(i) = 0, j(i) =0) reaction was studied on an accurate ab intio potential energy surface for the electronic ground state of BrH(2). Both the H + HBr → H(2) + Br abstraction reaction and the H' + HBr → H'Br + H exchange reaction were investigated up to a collision energy of 2.0 eV. It was found that the abstraction channel is dominant at lower collision energies, while the exchange channel becomes dominant at higher collision energies. The total integral cross section of the abstraction reaction at a collision energy of 1.6 eV was found to be 1.37 A?(2), which is larger than a recent quantum mechanical result (1.06 A?(2)) and still significantly smaller than the experimental value (3 ± 1 A?(2)). Meanwhile, similar to the previous theoretical study, our calculations also predicted much hotter product rotational state distributions than those from the experimental study. This suggests that further experimental investigations are highly desirable to elucidate the dynamic properties of the title reactions.     

Quantum effects in an exoergic, barrierless reaction at high collision energies

The exoergic Eley-Rideal hydrogen recombination on graphite surfaces is known to produce vibrationally hot product molecules. Realistic quantum scattering calculations at normal incidence over a wide range of collision energies show that the degree of vibrational excitation of the reaction product is a steep decreasing function of the collision energy. This results in selective population of the lower-lying vibrational levels and gives rise to an oscillating structure in the total reaction cross-sections at high energies, which cannot be reproduced by classical and quasi-classical trajectory calculations. An analogous quantum structure appears in the total collision-induced desorption cross-sections and is assigned to selective population of the metastable states of the transient molecular hydrogen.     

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

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