Quantum study of Eley-Rideal reaction and collision induced desorption of hydrogen atoms on a graphite surface. I. H-chemisorbed case

Collision induced (CI) processes involving hydrogen atoms on a graphite surface are studied quantum mechanically within the rigid, flat surface approximation, using a time-dependent wave packet method. The Eley-Rideal (ER) reaction and collision induced desorption (CID) cross sections are obtained with the help of two propagations which use different sets of coordinates, a "product" and a "reagent" set. Several adsorbate-substrate initial states of the target H atom in the chemisorption well are considered, and CI processes are studied over a wide range of projectile energy. Results show that (i) the Eley-Rideal reaction is the major reactive outcome and (ii) CID cross sections do not exceed 4 A2 and present dynamic thresholds for low values of the target vibrational quantum number. ER cross sections show oscillations at high energies which cannot be reproduced by classical and quasiclassical trajectory calculations. They are related to the vibrational excitation of the reaction products, which is a rather steep decreasing function of the collision energy. This behavior causes a selective population of the low-lying vibrational states and allows the quantization of the product molecular states to manifest itself in a collisional observable. A peak structure in the CID cross section is also observed and is assigned to the selective population of metastable states of the transient molecular hydrogen.     

I+HI'(v)→IH(v')+I'反应几率振荡行为的量子散射理论研究

本文基于振动绝热性分析, 用量子散射理论研究I+HI'(v)→IH(v')+I'反应几率的振荡行为, 其中双原子间相互作用势的振动本征态是Morse振子的严格解。在计算反应几率时采用分布高斯基(Distributed Gaussian Basis Sets; DGBS)展于方法, 既保持了反应几率的解析表达式, 又有效地简化了计算。计算结果明显地揭示了反应几率在低碰能区域的振荡行为。对反应机理也作了探讨。     

Quantum mechanical calculation of energy dependence of OCl/OH product branching ratio and product quantum state distributions for the O(1D) + HCl reaction on all three contributing electronic state potential energy surfaces

OCl/OH product branching ratios are calculated as a function of total energy for the O( (1) D) + HCl reaction using quantum wavepacket methods. The calculations take account of reaction on all the three electronic state potential energy surfaces which correlate with both reactants and products. Our results show that reaction on the excited electronic state surfaces has a large effect on the branching ratio at higher energies and that these surfaces must therefore be fully taken into account. The calculations use the potential energy surfaces of Nanbu and co-workers. Product vibrational and rotational quantum state distributions are also calculated as a function of energy for both product channels. Inclusion of the excited electronic state potential energy surfaces improves the agreement of the predicted product vibrational quantum state distributions with experiment for the OH product channel. For OCl agreement between theory and experiment is retained for the vibrational quantum state distributions when the excited electronic state potential energy surfaces are included in the analysis. For the rotational state distributions good agreement between theory and experiment is maintained for energies at which experimental results are available. At higher energies, above 0.7 eV of total energy, the OCl rotational state distributions predicted using all three electronic state potential energy surfaces shift to markedly smaller rotational quantum numbers.     

准经典轨线方法对Ca+CD3Ⅰ→CaⅠ+CD3同位素效应的动力学研究

在扩展Lond-Eyring-Polanyi-Sato(LEPS)势能面上,采用准经典轨线方法对反应Ca+CD3I→CaI+CD3进行了动力学计算,并讨论了该反应的同位素效应.在同位素效应作用下,产物CaI的振动态分布向低振动态转移,反应体系的散射截面在低碰撞能和高碰撞能处有较小的变化.同时,受到反应物的质量因子变化的影响,产物转动取向值减少,产物转动取向增强.仅有产物的角分布受同位素效应的影响不明显.     

Competition between dissociation and exchange processes in a collinear A + BC collision: Comparison of quantum and classical results

Previous calculations on a reactive/dissociative H + HD model system have been extended to higher collision energies. Exact quantum dissociation probabilities are now available in the 1–10 eV total energy range for initial vibrational quantum numbers v = 0,1 and 2. Comparison with quasiclassical results shows the absence of quantum tails at dissociation threshold, but large quantum effects at higher energies.     

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.     

Accurate time dependent wave packet calculations for the N + OH reaction

We present accurate quantum calculations of state-to-state cross sections for the N + OH → NO + H reaction performed on the ground (3)A' global adiabatic potential energy surface of Guadagnini et al. [J. Chem. Phys. 102, 774 (1995)]. The OH reagent is initially considered in the rovibrational state ν = 0, j = 0 and wave packet calculations have been performed for selected total angular momentum, J = 0, 10, 20, 30, 40,...,120. Converged integral state-to-state cross sections are obtained up to a collision energy of 0.5 eV, considering a maximum number of eight helicity components, Ω = 0,...,7. Reaction probabilities for J = 0 obtained as a function of collision energy, using the wave packet method, are compared with the recently published time-independent quantum mechanical one. Total reaction cross sections, state-specific rate constants, opacity functions, and product state-resolved integral cross-sections have been obtained by means of the wave packet method for several collision energies and compared with recent quasi-classical trajectory results obtained with the same potential energy surface. The rate constant for OH(ν = 0, j = 0) is in good agreement with the previous theoretical values, but in disagreement with the experimental data, except at 300 K.     

Vibrational mode and collision energy effects on reaction of H2CO+ with C2D4

We report the effects of collision energy (Ecol) and five different H2CO+ vibrational modes on the reaction of H2CO+ with C2D4 over the center-of-mass E(col) range from 0.1 to 2.1 eV. Properties of various complexes and transition states were also examined computationally. Seven product channels are observed. Charge transfer (CT) has the largest cross section over the entire energy range, substantially exceeding the hard sphere cross section at high energies. Competing with CT are six channels involving transfer of one or more hydrogen atoms or protons and one involving formation of propanal, followed by hydrogen elimination. Despite the existence of multiple deep wells on the potential surface, all reactions go by direct mechanisms, except at the lowest collision energies, where short-lived complexes appear to be important. Statistical complex decay appears adequate to account for the product branching at low collision energies, however, even at the lowest energies, the vibrational effects are counter to statistical expectations. The pattern of Ecol and vibrational mode effects provide insight into factors that control reaction and interchannel competition.     

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.     

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.     

H+ versus D+) transfer from HOD+ to N2: mode- and bond-selective effects

Reactions of HOD(+) with N(2) have been studied for HOD(+) in its ground state and with one quantum of excitation in each of its vibrational modes: (001)--predominately OH stretch, 0.396 eV, (010)--bend, 0.153 eV, and (100)--predominately OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 4 eV. The cross sections for both H(+) and D(+) transfer rise slowly from threshold with increasing collision energy; however, all three vibrational modes enhance reaction much more strongly than equivalent amounts of collision energy and the enhancements remain large even at high collision energy, where the vibration contributes less than 10% of the total energy. Excitation of the OH stretch enhances H(+) transfer by a factor of ～5, but the effect on D(+) transfer is only slightly larger than that from an equivalent increase in collision energy, and smaller than the effect from the much lower energy bend excitation. Similarly, OD stretch excitation strongly enhances D(+) transfer, but has essentially no effect beyond that of the additional energy on H(+) transfer. The effects of the two stretch vibrations are consistent with the expectation that stretching the bond that is broken in the reaction puts momentum in the correct coordinate to drive the system into the exit channel. From this perspective it is quite surprising that bend excitation also results in large (factor of 2) enhancements of both H(+) and D(+) transfer channels, such that its effect on the total cross section at collision energies below ～2 eV is comparable to those from the two stretch modes, even though the bend excitation energy is much smaller. For collision energies above ～2 eV, the vibrational effects become approximately proportional to the vibrational energy, though still much larger than the effects of equivalent addition of collision energy. Measurements of the product recoil velocity distributions show that reaction is direct at all collision energies, with roughly half the products in a sharp peak corresponding to stripping dynamics and half with a broad and approximately isotropic recoil velocity distribution. Despite the large effects of vibrational excitation on reactivity, the effects on recoil dynamics are small, indicating that vibrational excitation does not cause qualitative changes in the reaction mechanism or in the distribution of reactive impact parameters.     

Internal Energy State Distribution of Product CaBr in the Collision Reactions Ca+C2H5Br和Ca+nC3H7Br

The internal energy distributions of product CaBr in the collision reactions Ca+C2H5Br and Ca+nC3H7Br are studied by using the quasiclassical trajectory method. The average vibrational, rotational and translational energies and total available energies of the product CaBr molecules are calculated. The results indicate that when the collision energy is equal to 7.54 kJ/mol the energy of product CaBr is mainly the vibrational energy. As the reactant collision energy increases, the average translational and rotational energies of the product CaBr increase, the average vibrational energy decreases slightly, and the most probable vibrational state shifts to lower vibrational energy levels. The internal states of reagents have little influence on the internal energy distribution of the product. The bigger the radical group is, the higher ratio of the vibrational energy to the available energy of the product is. There exist two competitive reaction paths for the collision reactions Ca+C2H5Br and Ca+nC3H7Br, the migratory encounter and direct reaction paths. The former produces high vibrational excited state product CaBr and the latter causes C-Br bond to break. When the collision energy increases, the reactions tend to the latter path.     

Comparison of experimental time-of-flight spectra of the HF products from the F+H2 reaction with exact quantum mechanical calculations

High resolution HF product time-of-flight spectra measured for the reactive scattering of F atoms from n-H2(p-H2) molecules at collision energies between 69 and 81 meV are compared with exact coupled-channel quantum mechanical calculations based on the Stark-Werner ab initio ground state potential energy surface. Excellent agreement between the experimental and computed rotational distributions is found for the HF product vibrational states v'=1 and v'=2. For the v'=3 vibrational state the agreement, however, is less satisfactory, especially for the reaction with p-H2. The results for v'=1 and v'=2 confirm that the reaction dynamics for these product states is accurately described by the ground electronic state 1 (2)A' potential energy surface. The deviations for HF(v'=3, j' > or =2) are attributed to an enhancement of the reaction resulting from the 25% fraction of excited ((2)P(12)) fluorine atoms in the reactant beam.     

Vibrational mode and collision energy effects on reaction of H2CO+ with CO2

The effects of collision energy (Ecol) and five different modes of H2CO+ vibration on the title reaction have been studied over the center-of-mass Ecol range from 0.1 to 3.2 eV, including measurements of product ion recoil velocity distributions. Electronic structure and Rice-Ramsperger-Kassel-Marcus calculations were used to examine properties of various complexes and transition states that might be important along the reaction coordinate. Two product channels are observed, corresponding to Hydrogen Transfer (HT) and Proton Transfer (PT). Both channels are endothermic with similar onset energies of approximately 0.9 eV; however, HT dominates over the entire Ecol range and accounts for 70-85% of the total reaction cross section. Both HT and PT occur by direct mechanisms over the entire Ecol range, and have similar dependence on reactant vibrational and collision energy. Despite these similarities, and the fact that the two channels are nearly isoenergetic and differ only in which product moiety carries the charge, their dynamics appear quite different. PT occurs primarily in large impact parameter stripping collisions, where most of the available energy is partitioned to product recoil. HT, in contrast, results in internally hot products with little recoil energy and a more forward-backward symmetric product velocity distribution. Vibration is found to affect the reaction differently in different collision energy regimes. The appearance thresholds are found to depend only on total energy, i.e., all modes of vibration are equivalent to Ecol. With increasing Ecol, vibrational energy becomes increasingly effective, relative to Ecol, at driving reaction. For HT, this transition occurs just above threshold, while for PT it begins at roughly twice the threshold energy.     

Quantum dynamical study of the O(1D)+HCl reaction employing three electronic state potential energy surfaces

Quantum dynamical calculations are reported for the title reaction, for both product arrangement channels and using potential energy surfaces corresponding to the three electronic states, 1 1A', 2 1A', and 1 1A", which correlate with both reactants and products. The calculations have been performed for J=0 using the time-dependent real wavepacket approach by Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. Reaction probabilities for both product arrangement channels on all three potential energy surfaces are presented for total energies between 0.1 and 1.1 eV. Product vibrational state distributions at two total energies, 0.522 and 0.722 eV, are also presented for both channels and all three electronic states. Product rotational quantum state distributions are presented for both product arrangement channels and all three electronic states for the first six product vibrational states.     

A transition state view on reactive scattering: initial state-selected reaction probabilities for the H + CH4 → H2 + CH3 reaction studied in full dimensionality

Initial state-selected reaction probabilities for the H+CH(4)→H(2)+CH(3) reaction are computed for vanishing total angular momentum by full-dimensional calculations employing the multiconfigurational time-dependent Hartree approach. An ensemble of wave packets completely describing reactivity for total energies up to 0.58 eV is constructed in the transition state region by diagonalization of the thermal flux operator. These wave packets are then propagated into the reactant asymptotic region to obtain the initial state-selected reaction probabilities. Reaction probabilities for reactants in all rotational states of the vibrational 1A(1), 1F(2), and 1E levels of methane are presented. Vibrational excitation is found to decrease reactivity when reaction probabilities at equivalent total energies are compared but to increase reaction probabilities when the comparison is done at the basis of equivalent collision energies. Only a fraction of the initial vibrational energy can be utilized to promote the reaction. The effect of rotational excitation on the reactivity differs depending on the initial vibrational state of methane. For the 1A(1) and 1F(2) vibrational states of methane, rotational excitation decreases the reaction probability even when comparing reaction probabilities at equivalent collision energies. In contrast, rotational energy is even more efficient than translational energy in increasing the reaction probability when the reaction starts from the 1E vibrational state of methane. All findings can be explained employing a transition state based interpretation of the reaction process.     

State to State Photodissociation Dynamics of Vibrationally Excited D2O in B Band

The state-to-state photodissociassion dynamics for the B band of D₂O have been explored from quantum dynamical calculations including the electronic ~X and ~B states. The calculations were carried out using a Chebyshev real wave packet method. The calculated absorption spectra, product state distributions, and branching ratios from different initial vibrational states show di?erent dynamic features, due to the different shapes of the vibrational wavefunctions. The initial bending mode (0,1,0) generates two lobes with a shallow minimum on the absorption spectrum and a slight inverted vibrational population of OD(~X )product at high total energies. The rotational state distributions of OD(~X , v=0) product are highly inverted and depend weakly on the initial state and total energy. On the other hand, the ro-vibrational distributions of OD(A~) product strongly oscillate with the total energy, which are dominated by the long-living resonances and depend sensitively on the potential surfaces. The antisymmetric stretching mode (0,0,1) has large OD( ~ A)/OD(~X ) branching ratios at high total energies, which indicates that the B band dissociation proceeds mainly via the adiabatic pathway in some cases.     

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.     

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.