Quantum dynamics of dissociative adsorption of H2 on Ni(100) using the dynamical Lie algebraic method

A dynamical Lie algebraic method has been applied to treating the quantum dynamics of dissociative adsorption of H₂ on a static flat metal surface. An LEPS potential energy surface has been used to describe the interaction of H₂ with Ni(100) surface. The dependence of the initial state-selected dissociation probability was obtained analytically on the initial kinetic energy and time. A comparison with other theoretical calculations and experiments is made. The results show that the method can be effectively used to describe the dynamics of reactive gas-sdace scattering. Project supported by the National Natural Science Foundation of China (Grant No. 19694033) and partially by the State Key Laboratory of Theoretical and Computational Chemistry of Jilin University (Grant No. 9801).     

Communication: state-to-state quantum dynamics study of the OH + CO → H + CO2 reaction in full dimensions (J = 0)

A full dimensional state-to-state quantum dynamics study is carried out for the prototypical complex-formation OH + CO → H + CO(2) reaction in the ground rovibrational initial state on the Lakin-Troya-Schatz-Harding potential energy surface by using the reactant-product decoupling method. With three heavy atoms and deep wells on the reaction path, the reaction represents a huge challenge for accurate quantum dynamics study. This state-to-state calculation is the first such a study on a four-atom reaction other than the H(2) + OH ? H(2)O + H and its isotope analogies. The product CO(2) vibrational and rotational state distributions, and product energy partitioning information are presented for ground initial rovibrational state with the total angular momentum J = 0.     

Wave packet dynamics of H2(v1=8-14)+H2(v2=0-2): the role of the potential energy surface on different reactive and dissociative processes

A time-dependent wave packet method has been used to study different competing products of H(2)+H(2) collisions: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. A three-degree-of-freedom reduced dimensionality model has been used for five different geometries of the colliding complex (parallel H, crossed X, collinear L, and two T-shaped geometries T(I) and T(II)), with reactants in selected vibrational states with one diatom vibrationally "hot" and the other one vibrationally "cold." Product probabilities have been calculated using two potential energy surfaces [J. Chem. Phys. 101, 4004 (1994); J. Chem. Phys. 116, 666 (2002)] in order to compare their performance in the dynamics. The regions of the potential energy surfaces responsible of the threshold behavior of the probabilities have been identified. Overall, we have found that the most recent potential energy surface is less anisotropic, provides a smaller propensity for insertion-type processes, and gives lower energy thresholds.     

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.     

D_2在Ni(100)表面离解吸附的量子动力学研究

用量子含时波包法研究了Ｄ２在镍表面顶－桥位上离解吸附量子动力学．计算了不同入射动能及初始振转态的离解几率．讨论了分子的同核对称性、转动取向和振动激发对离解几率的影响，并与其他理论计算结果做了比较．     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

Diffusion dynamics of the li atom on amorphous carbon: A direct molecular orbital-molecular dynamics study

Direct molecular orbital-molecular dynamics (MO-MD) calculation was applied to diffusion processes of the Li atom on a model surface of amorphous carbon and compared with the diffusion mechanism of Li+ ion. A carbon sheet composed of C96H24 was used as the model surface. The total energy and energy gradient on the full dimensional potential energy surface of the LiC96H24 system were calculated at each time step in the trajectory calculation. The optimized structure, where the Li atom is located at the center of mass of the model surface, was used as the initial structure at time zero. Simulation temperatures were chosen in the range of 200-1250 K. The dynamics calculations showed that the Li atom vibrates around the initial position below 250 K, and it moves above 300 K. At middle temperature, the Li atom translates freely on the surface. At higher temperature (1000 K), the Li atom moves from the center to edge region of the model surface and is trapped in the edge. The activation energy calculated for the Li atom is larger than that for the Li+ ion. This difference is due to the fact that the Li atom diffuses together with an unpaired electron on the carbon surface. The diffusion mechanism of the Li atom was discussed on the basis of the theoretical results.     

Adsorption dynamics of water on Pt{110}-(1 x 2)

The dynamics of H(2)O adsorption on Pt{110}-(1 x 2) is studied using supersonic molecular beam and temperature programed desorption techniques. The sticking probabilities are measured using the King and Wells method at a surface temperature of 165 K. The absolute initial sticking probability s(0) of H(2)O is 0.54+/-0.03 for an incident kinetic energy of 27 kJmol. However, an unusual molecular beam flux dependence on s(0) is also found. At low water coverage (theta<1), the sticking probability is independent of coverage due either to diffusion in an extrinsic precursor state formed above bilayer islands or to incorporation into the islands. We define theta=1 as the water coverage when the dissociative sticking probability of D(2) on a surface predosed with water has dropped to zero. The slow falling H(2)O sticking probability at theta>1 results from compression of the bilayer and the formation of multilayers. Temperature programed desorption of water shows fractional order kinetics consistent with hydrogen-bonded islands on the surface. A remarkable dependence of the initial sticking probability on the translational (1-27 kJ/mol) and internal energies of water is observed: s(0) is found to be essentially a step function of translational energy, increasing fivefold at a threshold energy of 5 kJ/mol. The threshold migrates to higher energies with increasing nozzle temperature (300-700 K). We conclude that both rotational state and rotational alignment of the water molecules in the seeded supersonic expansion are implicated in dictating the adsorption process.     

Diffusion dynamics of the Li+ ion on a model surface of amorphous carbon: a direct molecular orbital dynamics study

Diffusion processes of the Li+ ion on a model surface of amorphous carbon (Li+C96H24 system) have been investigated by means of the direct molecular orbital (MO) dynamics method at the semiempirical AM1 level. The total energy and energy gradient on the full-dimensional AM1 potential energy surface were calculated at each time step in the dynamics calculation. The optimized structure, where Li+ is located in the center of the cluster, was used as the initial structure at time zero. The dynamics calculation was carried out in the temperature range 100-1000 K. The calculations showed that the Li+ ion vibrates around the equilibrium point below 200 K, while the Li+ ion moves on the surface above 250 K. At intermediate temperatures (300 K < T < 400 K), the ion moves on the surface and falls in the edge regions of the cluster. At higher temperatures (600 K < T), the Li+ ion transfers freely on the surface and edge regions. The diffusion pathway of the Li+ ion was discussed on the basis of theoretical results.     

A new ab initio interaction energy surface and high-resolution spectra of the H2-CO van der Waals complex

A new four-dimensional intermolecular potential-energy surface for the H(2)-CO complex is presented. The ab initio points have been computed on a five-dimensional grid including the dependence on the H-H separation (the C-O separation was fixed). The surface has then been obtained by averaging over the intramolecular vibration of H(2). The coupled-cluster supermolecular method with single, double, and noniterative triple excitations has been used to calculate the interaction energy. The correlation part of the interaction energy has been obtained from extrapolations based on calculations in a series of basis sets. An analytical fit of the ab initio potential-energy surface has the global minimum of -93.049 cm(-1) at the intermolecular separation of 7.92 bohr for the linear geometry with the C atom pointing toward the H(2) molecule. For the other linear geometry, with the O atom pointing toward H(2), the local minimum of -72.741 cm(-1) has been found for the intermolecular separation of 7.17 bohr. The potential has been used to calculate the rovibrational energy levels of the para-H(2)-CO complex. The results agree very well with those observed by McKellar [A. R. W. McKellar J. Chem. Phys. 108, 1811 (1998)]: the discrepancies are smaller than 0.1 cm(-1). The calculated dissociation energy is equal to 19.527 cm(-1) and significantly smaller than the value of 22 cm(-1) estimated from the experiment. Predictions of rovibrational energy levels for ortho-H(2)-CO have also been done and can serve as a guidance to assign recorded experimental spectra. The interaction second virial coefficient has been calculated and compared with the experimental data.     

State-selective predissociation dynamics of methylamines: the vibronic and HD effects on the conical intersection dynamics

The photodissociation dynamics of methylamines (CH(3)NH(2) and CD(3)ND(2)) on the first electronically excited state has been investigated using the velocity map ion imaging technique probing the H or D fragment. Two distinct velocity components are found in the H(D) translational energy distribution, implying the existence of two different reaction pathways for the bond dissociation. The high H(D) velocity component with the small internal energy of the radical fragment is ascribed to the N-H(D) fragmentation via the coupling of S(1) to the upper-lying S(2) repulsive potential energy surface along the N-H(D) bond elongation axis. Dissociation on the ground S(0) state prepared via the nonadiabatic dynamics at the conical intersection should be responsible for the slow H(D) fragment. Several S(1) vibronic states of methylamines including the zero-point level and nnu(9) states (n=1, 2, or 3) are exclusively chosen in order to explore the effect of the initial quantum content on the chemical reaction dynamics. The branching ratio of the fast and slow components is found to be sensitive to the initial vibronic state for the N-H bond dissociation of CH(3)NH(2), whereas it is little affected in the N-D dissociation event of CD(3)ND(2). The fast component is found to be more dominant in the translational distribution of D from CD(3)ND(2) than it is in that of H from CH(3)NH(2). The experimental result is discussed with a plausible mechanism of the conical intersection dynamics.     

Competition between dissociation and exchange processes in a collinear A + BC collision. I. Exact quantum results

Exact quantum results for collision-induced dissociation on a reactive surface are presented. A modified LEPS potential-energy surface modeling the H + HD → H₂ + D system has been used. HD and H₂ bearing respectively 7 and 6 bound states. This system has been chosen because it displays significant reactive scattering for total energies above the dissociation threshold. Calculations have been performed using the time-dependent wavepacket method for two initial vibrational quantum numbers of the HD molecule (v = 0, 2). For each vibrational quantum number, two wavepackets with overlapping energy distributions have been run, covering a total energy range up to more than three times the dissociation energy. Comparison with previous collision-induced dissociation calculations shows that the dissociation is greatly enhanced by the presence of concomitant reactive scattering. A vibrational enhancement effect is also observed above the dissociation threshold; for higher energies the system exhibits a pronounced vibrational inhibition effect.     

Influence of nuclear exchange on nonadiabatic electron processes in H(+)+H2 collisions

H(+)+H(2) collisions are studied by means of a semiclassical approach that explicitly accounts for nuclear rearrangement channels in nonadiabatic electron processes. A set of classical trajectories is used to describe the nuclear motion, while the electronic degrees of freedom are treated quantum mechanically in terms of a three-state expansion of the collision wavefunction. We describe electron capture and vibrational excitation, which can also involve nuclear exchange and dissociation, in the E = 2-1000 eV impact energy range. We compare dynamical results obtained with two parametrizations of the potential energy surface of H(3)(+) ground electronic state. Total cross sections for E > 10 eV agree with previous results using a vibronic close-coupling expansion, and with experimental data for E < 10 eV. Additionally, some prototypical features of both nuclear and electron dynamics at low E are discussed.     

Hyperthermal reactions of O and O2 with a hydrocarbon surface: direct C-C bond breakage by O and H-atom abstraction by O2

A C-C bond-breaking reaction has been observed when a beam containing hyperthermal oxygen was directed at a continuously refreshed saturated hydrocarbon liquid (squalane) surface. The dynamics of this C-C bond-breaking reaction have been investigated by monitoring time-of-flight and angular distributions of the volatile product, OCH3 or H2CO. The primary product is believed to be the methoxy radical, OCH3, but if this radical is highly internally excited, then it may undergo secondary dissociation to form formaldehyde, H2CO. Either the primary or the secondary product may scatter directly into the gas phase before thermal equilibrium with the surface is reached, or they may become trapped on the surface and desorb in thermal equilibrium with the surface. Direct, single-collision scattering events that produce a C-C bond-breaking product are described with a kinematic picture that allows the determination of the effective surface mass encountered by an incident O atom, the atom-surface collision energy in the center-of-mass frame, and the fraction of the center-of-mass collision energy that goes into translation of the scattered gaseous product and the recoiling surface fragment. The dynamical behavior of the C-C bond-breaking reaction is compared with that of the H-atom abstraction reaction, which was the subject of an earlier study. Another reaction, H-atom abstraction by O2 (which is present in the hyperthermal beam), has also been observed, and the dynamics of this reaction are compared with the inelastic scattering dynamics of O2 and the dynamics of H-atom abstraction by O. The dynamics involving direct inelastic and reactive scattering of O2 are also described in terms of a kinematic picture where the incident O2 molecule is viewed as interacting with a local region of the surface.     

Classical trajectory study of the dynamics of the reaction of Cl atoms with ethane

We present an electronic-structure and dynamics study of the Cl + C2H6 --> HCl + C2H5 reaction. The stationary points of the ground-state potential energy surface have been characterized using various electronic-structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis limit, using geometries and harmonic frequencies obtained at the MP2/aug-cc-pVTZ level, are in agreement with the experimental reaction energy. Ab initio information has been used to reparameterize a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in key regions of the potential energy surface. The improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories. Computed kinetic energy release and scattering angle distributions at a collision energy of approximately 5.5 kcal mol(-1) are in reasonable agreement with experiments, but no evidence was found for the low translational energy HCl products scattered in the backward hemisphere reported in recent experiments.     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

Photodissociation of hydrogen iodide on the surface of large argon clusters: the orientation of the librational wave function and the scattering from the cluster cage

A set of photodissociation experiments and simulations of hydrogen iodide (HI) on Arn clusters, with an average size n = 139, has been carried out for different laser polarizations. The doped clusters are prepared by a pick-up process. The HI molecule is then photodissociated by a UV laser pulse and the outgoing H fragment is ionized by resonance enhanced multiphoton ionization in a (2 + 1) excitation scheme within the same laser pulse at the wavelength of 243 nm. The measured time-of-flight spectra are transformed into hydrogen kinetic energy distributions. They exhibit a strong fraction of caged H atoms at zero-kinetic energy and peaks at the unperturbed cage exit for both spin-orbit channels nearly independent of the polarization. At this dissociation wavelength, the bare HI molecule exhibits a strict state separation, with a parallel transition to the spin-orbit excited state and perpendicular transitions to the ground state. The experimental results have been reproduced using molecular simulation techniques. Classical molecular dynamics was used to estimate the HI dopant distribution after the pick-up procedure. Subsequently, quasi-classical molecular dynamics (Wigner trajectories approach) has been applied for the photodissociation dynamics. The following main results have been obtained: (i) The HI dopant lands on the surface of the argon cluster during the pick-up process, (ii) zero-point energy plays a dominant role for the hydrogen orientation in the ground state of HI-Arn surface clusters, qualitatively changing the result of the photodissociation experiment upon increasing the number of argon atoms, and, finally, (iii) the scattering of hydrogen atoms from the cage which originate from different dissociation states seriously affects the experimentally measured kinetic energy distributions.