Quantized dynamical bottlenecks and transition state control of the reaction of D with H2: Effect of varying the total angular momentum

Rate coefficients for the O + H2 and O + D2 reactions: how well ring polymer molecular dynamics accounts for tunelling

We present here extensive calculations of the O(3P) + H2 and O(3P) + D2 reaction dynamics spanning the temperature range from 200 K to 2500 K. The calculations have been carried out using fully converged time-independent quantum mechanics (TI QM), quasiclassical trajectories (QCT) and ring polymer molecular dynamics (RPMD) on the two lowest lying adiabatic potential energy surfaces (PESs), 13A' and 13A'', calculated by Zanchet et al. [J. Chem. Phys., 2019, 151, 094307]. TI QM rate coefficients were determined using the cumulative reaction probability formalism on each PES including all of the total angular momenta and the Coriolis coupling and can be considered to be essentially exact within the Born-Oppenheimer approximation. The agreement between the rate coefficients calculated by using QM and RPMD is excellent for the reaction with D2 in almost the whole temperature range. For the reaction with H2, although the agreement is very good above 500 K, the deviations are significant at lower temperatures. In contrast, the QCT calculations largely underestimate the rate coefficients for the two isotopic variants due to their inability to account for tunelling. The differences found in the disagreements between RPMD and QM rate coefficients for the reactions for both the isotopologues are indicative of the ability of the RPMD method to accurately describe systems where tunelling plays a relevant role. Considering that both reactions are dominated by tunelling below 500 K, the present results show that RPMD is a very powerful tool for determining rate coefficients. The present QM rate coefficients calculated on adiabatic PESs slightly underestimate the best global fits of the experimental measurements, which we attribute to the intersystem crossing with the singlet 11A' PES.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Inelastic Scattering of H Atoms from Surfaces

We have developed an instrument that uses photolysis of hydrogen halides to produce nearly monoenergetic hydrogen atom beams and Rydberg atom tagging to obtain accurate angle-resolved time-of-flight distributions of atoms scattered from surfaces. The surfaces are prepared under strict ultrahigh vacuum conditions. Data from these experiments can provide excellent benchmarks for theory, from which it is possible to obtain an atomic scale understanding of the underlying dynamical processes governing H atom adsorption. In this way, the mechanism of adsorption on metals is revealed, showing a penetration-resurfacing mechanism that relies on electronic excitation of the metal by the H atom to succeed. Contrasting this, when H atoms collide at graphene surfaces, the dynamics of bond formation involving at least four carbon atoms govern adsorption. Future perspectives of H atom scattering from surfaces are also outlined.

New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O(3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results

In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(³P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the ³A' and ³A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the ³A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.

State-to-State Dynamics of Elementary Bimolecular Reactions

The study of state-to-state dynamics of elementary bimolecular reactions has provided remarkable insights into chemical reactivity at the most fundamental level. This review covers exciting developments in this important field in the past decade. I focus on recent studies of quantum-state-resolved molecular-beam reactive-scattering studies of elementary chemical reactions, from triatomic to polyatomic systems. Researchers have made great advances in the fundamental understanding of many elementary chemical reactions through state-to-state dynamics studies. The strong interaction between theory and experiment has significantly enhanced our understanding of the dynamics of these reactions. I hope this review provides a glimpse of this exciting research field to both experts and beginners.

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.

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.

Theories of reactive scattering

This paper is an overview of the theory of reactive scattering, with emphasis on fully quantum mechanical theories that have been developed to describe simple chemical reactions, especially atom-diatom reactions. We also describe related quasiclassical trajectory applications, and in all of this review the emphasis is on methods and applications concerned with state-resolved reaction dynamics. The review first provides an overview of the development of the theory, including a discussion of computational methods based on coupled channel calculations, variational methods, and wave packet methods. Choices of coordinates, including the use of hyperspherical coordinates are discussed, as are basis set and discrete variational representations. The review also summarizes a number of applications that have been performed, especially the two most comprehensively studied systems, H+H2 and F+H2, along with brief discussions of a large number of other systems, including other hydrogen atom transfer reactions, insertion reactions, electronically nonadiabatic reactions, and reactions involving four or more atoms. For each reaction we describe the method used and important new physical insight extracted from the results.

The state-to-state-to-state model for direct chemical reactions: Application to D+H2→HD+H

A simple theoretical model is developed to predict the state-to-state dynamics of direct chemical reactions. Motivated by traditional ideas from transition state theory, expressions are derived for the reactive S matrix that may be computed using the local transition state dynamics. The key approximation involves the use of quantum bottleneck states to represent the near separable dynamics taking place near the transition state. Explicit expressions for the S matrix are obtained using a Franck-Condon treatment for the inelastic coupling between internal states of the collision complex. It is demonstrated that the energetic thresholds for various initial reagent states of the D+H(2) reaction can be understood in terms of our theory. Specifically, the helicity of the reagent states are found to correlate directly to the symmetry of the quantum bottleneck states, which thus possess very different thresholds. Furthermore, the rotational product state distributions for D+H(2) are found to be associated with interfering pathways through the quantum bottleneck states.

Quasiclassical determination of reaction probabilities as a function of the total angular momentum

This article presents a quasiclassical trajectory (QCT) method to determine the reaction probability as a function of the total angular momentum J for any given value of the initial rotational angular momentum j. The proposed method is based on a discrete sampling of the total and orbital angular momenta for each trajectory and on the development of equations that have a clear counterpart in the quantum-mechanical (QM) case. The reliability of the method is illustrated by comparing QCT and time-dependent wave-packet QM results for the H+D(2)(upsilon=0,j=4,10) reaction. The small discrepancies between both sets of calculations, when they exist, indicate some genuine quantum effects. In addition, a procedure to extract the reaction probabilities as a function of J when trajectories are calculated in the usual way using a continuous distribution of impact parameters is also described.

Quantum mechanical reaction rate constants by vibrational configuration interaction: the OH + H2->H2O + H reaction as a function of temperature

The thermal rate constant of the 3D OH + H(2)-->H(2)O + H reaction was computed by using the flux autocorrelation function, with a time-independent square-integrable basis set. Two modes that actively participate in bond making and bond breaking were treated by using 2D distributed Gaussian functions, and the remaining (nonreactive) modes were treated by using harmonic oscillator functions. The finite-basis eigenvalues and eigenvectors of the Hamiltonian were obtained by solving the resulting generalized eigenvalue equation, and the flux autocorrelation function for a dividing surface optimized in reduced-dimensionality calculations was represented in the basis formed by the eigenvectors of the Hamiltonian. The rate constant was obtained by integrating the flux autocorrelation function. The choice of the final time to which the integration is carried was determined by a plateau criterion. The potential energy surface was from Wu, Schatz, Lendvay, Fang, and Harding (WSLFH). We also studied the collinear H + H(2) reaction by using the Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential energy surface. The calculated thermal rate constant results were compared with reported values on the same surfaces. The success of these calculations demonstrates that time-independent vibrational configuration interaction can be a very convenient way to calculate converged quantum mechanical rate constants, and it opens the possibility of calculating converged rate constants for much larger reactions than have been treated until now.

Ab initiorate constants from hyperspherical quantum scattering: Application to H+CH4→H2+CH3

A general and practical procedure is described for calculating rate constants for chemical reactions using a minimal number of ab initio calculations and quantum-dynamical computations. The method exploits a smooth interpolating functional developed in the hyperspherical representation. This functional is built from two Morse functions and depends on a relatively small number of parameters with respect to conventional functionals developed to date. Thus only a small number of ab initio points needs to be computed. The method is applied to the H + CH4 --> H2 + CH3 reaction. The quantum scattering calculations are performed treating explicitly the bonds being broken and formed. All the degrees of freedom except the breaking and forming bonds are optimized ab initio and harmonic vibrational frequencies and zero-point energies for them are calculated at the MP2(full) level with a cc-pVTZ basis set. Single point energies are calculated at a higher level of theory with the same basis set, namely CCSD(T, full). We report state-to-state cross sections and thermal rate constants for the title reaction and make comparisons with previous results. The calculated rate constants are in good agreement with experiments.

Interference of quantized transition-state pathways in the H + D2 -> D + HD chemical reaction

The collision-energy dependence of the state-resolved differential cross section at a specific backward-scattering angle for the reaction H + D2 --> D + HD is measured with the D-atom Rydberg "tagging" time-of-flight technique. The reaction was modeled theoretically with converged quantum scattering calculations that provided physical interpretation of the observations. Oscillations in the differential cross sections in the backward-scattering direction are clearly observed and are attributed to the transition-state structures that originate from the interferences of different quantized transition-state pathways.

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.

Forward scattering due to slow-down of the intermediate in the H + HD --> D + H(2) reaction

Quantum dynamical processes near the energy barrier that separates reactants from products influence the detailed mechanism by which elementary chemical reactions occur. In fact, these processes can change the product scattering behaviour from that expected from simple collision considerations, as seen in the two classical reactions F + H(2) --> HF + H and H + H(2) --> H(2) + H and their isotopic variants. In the case of the F + HD reaction, the role of a quantized trapped Feshbach resonance state had been directly determined, confirming previous conclusions that Feshbach resonances cause state-specific forward scattering of product molecules. Forward scattering has also been observed in the H + D(2) --> HD + D reaction and attributed to a time-delayed mechanism. But despite extensive experimental and theoretical investigations, the details of the mechanism remain unclear. Here we present crossed-beam scattering experiments and quantum calculations on the H + HD --> H(2) + D reaction. We find that the motion of the system along the reaction coordinate slows down as it approaches the top of the reaction barrier, thereby allowing vibrations perpendicular to the reaction coordinate and forward scattering. The reaction thus proceeds, as previously suggested, through a well-defined 'quantized bottleneck state' different from the trapped Feshbach resonance states observed before.

Scattering resonances in the simplest chemical reaction

Recent studies of state-resolved angular distributions show the participation of reactive scattering resonances in the simplest chemical reaction. This review is intended for those who wish to learn about the state-of-the-art in the study of the H + H2 reaction family that has made this breakthrough possible. This review is also intended for those who wish to gain insight into the nature of reactive scattering resonances. Following a tour across several fields of physics and chemistry where the concept of resonance has been crucial for the understanding of new phenomena, we offer an operational definition and taxonomy of reactive scattering resonances. We introduce simple intuitive models to illustrate each resonance type. We focus next on the last decade of H + H2 reaction dynamics. Emphasis is placed on the various experimental approaches that have been applied to the search for resonance behavior in the H + H2 reaction family. We conclude by sketching the road ahead in the study of H + H2 reactive scattering resonances.