Quantum mechanical calculation of product state distributions for the O(1D)+H2→OH+H reaction on the ground electronic state surface

Hydrogen atom quantum diffusion in solid parahydrogen: The H + N2O → cis-HNNO → trans-HNNO reaction

The diffusion and reactivity of hydrogen atoms in solid parahydrogen at temperatures between 1.5 K and 4.3 K are investigated by high-resolution infrared spectroscopy. Hydrogen atoms are produced within solid parahydrogen as the by-products of the 193 nm in situ photolysis of N₂O, which induces a two-step tunneling reaction, H + N₂O → cis-HNNO → trans-HNNO. The second-order rate constant for the first step to form cis-HNNO is found to be inversely proportional to the N₂O concentration after photolysis, indicating that the hydrogen atoms move through solid parahydrogen via quantum diffusion. This reaction only readily occurs at temperatures below 2.8 K, not due to an increased rate constant for the first reaction step at low temperatures but rather due to an increased selectivity to the reaction. The rate constant for the second step of the reaction mechanism involving unimolecular isomerization is shown to be independent of the N₂O concentration as expected. The inverse concentration dependence of the rate constant for the reaction step that involves the hydrogen atom demonstrates clearly that quantum diffusion influences the reactivity of the hydrogen atoms in solid parahydrogen, which does not have an analogy in classical reaction kinetics.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Dynamics study on the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction: insertion-abstraction mechanism

Time-dependent wave packet calculations are carried out for two reaction channels of the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction. The potential well on the excited state potential energy surface makes the reaction preferable to proceed through the insertion reaction path. The dominance of the NaD + H reaction channel and product rotational state distributions are found to be in agreement with the characteristics of typical adiabatic insertion reactions. However, significant forward scattering peaks in the differential cross sections (DCS) are found to be inconsistent with the forward-backward symmetric scattering characteristic of typical adiabatic insertion reactions, which indicate that the Na(3p) + HD reaction is dominated by a direct reaction mechanism. The comparison between adiabatic and non-adiabatic calculated DCSs reveals that the non-adiabatic couplings in the reaction could reduce the lifetime of the intermediate complex. Finally, the insertion-abstraction mechanism is put forward for the non-adiabatic Na(3p) + HD reaction.

A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

The O(¹D) + H₂ reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(¹D) + D₂ reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.

A review of dynamical resonances in A + BC chemical reactions

The concept of the transition state has played an important role in the field of chemical kinetics and reaction dynamics. Reactive resonances in the transition-state region can dramatically enhance the reaction probability; thus investigation of the reactive resonances has attracted great attention from chemical physicists for many decades. In this review, we mainly focus on the recent progress made in probing the elusive resonance phenomenon in the simple A  +  BC reaction and understanding its nature, especially in the benchmark F/Cl  +  H₂ and their isotopic variants. The signatures of reactive resonances in the integral cross section, differential cross section (DCS), forward- and backward-scattered DCS, and anion photodetachment spectroscopy are comprehensively presented in individual prototype reactions. The dynamical origins of reactive resonances are also discussed in this review, based on information on the wave function in the transition-state region obtained by time-dependent quantum wave-packet calculations.

CHEMICAL REACTION DYNAMICS WITH REAL WAVE PACKETS

An approach to carrying out accurate quantum dynamics simulations of chemical reactions, termed the real wave packet (RWP) method, is outlined. The method focuses on propagation of just the real part of a complex-valued wave packet, halving computational memory and effort requirements in comparison with comparable high accuracy quantum propagation methods. Applications to 3-atom and 4-atom chemical reactions are reviewed. Potential future directions are indicated.

LANCZOS ALGORITHMS AND CROSS CORRELATION FUNCTIONS Cif(E)

Cross correlation (CC) functions Cif(E) play an important role in chemical physics. They appear in the description of reactive scattering, photo-dissociation, photo-electron spectroscopy and electron transfer to mention a few. In this paper, we discuss two methods based on the Lanczos algorithm to compute the CC function for several initial and final states at the same time, without diagonalization. The property of the coupled two-term recursions variant of the Lanczos algorithm that yields a decomposition [Formula: see text] of the tridiagonal Lanczos matrix is crucial. The first method, the quasi minimal-recursive residue generation method (QM-RRGM) is based on solving a set of linear equations whereas the second method is based on a band-Lanczos method. The computational cost is of the same order of magnitude for both methods and is given by the number of matrix-vector multiplications in the underlying Lanczos method. Only a small set of scalars needs to be updated each recursion. The methods are compared for a model problem, the continuum resonance Raman cross section for a collinear model of CH2IBr . Both methods shows similar convergence properties. By adding a pre-conditioner, the rate of convergence can be increased dramatically.

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.

State-to-state reaction probabilities for the H+O2(v,j)→O+OH(v′,j′) reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.

Coriolis coupling effects in the calculation of state-to-state integral and differential cross sections for the H+D2 reaction

The quantum wavepacket parallel computational code DIFFREALWAVE is used to calculate state-to-state integral and differential cross sections for the title reaction on the BKMP2 surface in the total energy range of 0.4-1.2 eV with D2 initially in its ground vibrational-rotational state. The role of Coriolis couplings in the state-to-state quantum calculations is examined in detail. Comparison of the results from calculations including the full Coriolis coupling and those using the centrifugal sudden approximation demonstrates that both the energy dependence and the angular dependence of the calculated cross sections are extremely sensitive to the Coriolis coupling, thus emphasizing the importance of including it correctly in an accurate state-to-state calculation.

State-to-state reactive differential cross sections for the H+H2→H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.

Time dependent quantum dynamics study of the O++H2(v=0,j=0)→OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.

Reactions of C(1D) with H2 and its deuterated isotopomers, a wave packet study

Using a Chebyshev wave packet method, total and state-resolved reaction probabilities (J=0) were calculated for the reactions of C(1D) with various hydrogen isotopomers (H2, D2, and HD, nu i=0, j i=0) on a recent ab initio potential energy surface. For all the isotopic variants, it was found that the initial state specified reaction probabilities have no energy threshold and are strongly oscillatory, indicative of the involvement of long-lived resonances in this barrierless reaction. The J=0 product vibrational and rotational distributions for all three isotopic reactions, and the CH/CD branching ratio for the C+HD reaction, show strong dependence on the collision energy, further underscoring the important role played by the resonances. The generally decaying vibrational distributions and highly excited rotational distributions, which corroborate an insertion mechanism, and the dominance of the CD+H channel in the C+HD reaction are consistent with existing experimental observations. Initial state specified integral cross sections and rate constants were estimated using a capture model. The estimated rate constants were found to be close and in the order kHD>kH2>kD2. Finally, a method to calculate branching ratio in the C+HD reaction is proposed.

A wave packet based statistical approach to complex-forming reactions

A wave packet based statistical model is suggested for complex-forming reactions. This model assumes statistical formation and decay of the long-lived reaction complex and computes reaction cross sections and their energy dependence from capture probabilities. This model is very efficient and reasonably accurate for reactions dominated by long-lived resonances, as confirmed by its application to the C((1)D)+H(2) reaction.

Quantum scattering calculations on chemical reactions

This review discusses recent quantum scattering calculations on bimolecular chemical reactions in the gas phase. This theory provides detailed and accurate predictions on the dynamics and kinetics of reactions containing three atoms. In addition, the method can now be applied to reactions involving polyatomic molecules. Results obtained with both time-independent and time-dependent quantum dynamical methods are described. The review emphasises the recent development in time-dependent wave packet theories and the applications of reduced dimensionality approaches for treating polyatomic reactions. Calculations on over 40 different reactions are described.