Narrow Subthreshold Quantum Mechanical Resonances in the Li + HF → H + LiF Reaction

Quantum Dynamics Calculations of Na (32S, 32P) + HF → NaF + H Reactions

The dynamics of the reactive scatterings of the ground (Na(3S)) and first excited (Na(3P)) state sodium atoms from hydrogen fluoride (HF) molecules is investigated by performing comprehensive Coriolis-coupled quantum wave packet calculations on the recent ab initio potential energy surface (PES). In the Na(3S) + HF reaction, the nonadiabatic effect is found to be negligibly small, and the reactivity is significantly promoted by the initial vibrational excitations being in line with the available experimental result. Excellent quantitative agreement between theory and experiment is also achieved for the initial state specified integral cross sections and rate constants for v ≥ 3 vibrational states. However, the calculated rate constant for v = 2 significantly underestimates the experimental result. The possible cause for the disagreement is discussed in detail. For the Na(3P) + HF scattering, which can lead to the formation of either the ground state NaF + H product or Na(3S) + HF reactant via the harpooning process, the calculated result for the integral cross section shows excellent agreement with the available experimental result indicating the reasonable accuracy of the interstate coupling term of the employed PES.

Partial Potential Energy Surfaces and Their Application to Reaction Resonances

Feshbach resonances are not restricted to small reactive systems such as F + H2 and I + HI but can be found in many reactive systems. In this paper, the concept of the partial potential energy surface (PPES) is introduced. It is shown that the dynamic “Lake Eyring” explains very well the existence of reactive resonances in elementary chemical reactions. In particular, the PPESs of the Cl + CH3CH2Br and Cl + CH3CH2CH2Br systems, including the minimum energy path and the vibrational potential energy curves, were constructed using quantum chemistry methods. Based on the constructed PPESs, the scattering resonance states of these reactions could be examined and the resonance state lifetimes were estimated.

A new ab initio potential energy surface of LiClH (1A') system and quantum dynamics calculation for Li + HCl (v = 0, j = 0-2) → LiCl + H reaction

A new ab initio potential energy surface (PES) for the ground state of Li + HCl reactive system has been constructed by three-dimensional cubic spline interpolation of 36 654 ab initio points computed at the MRCI+Q/aug-cc-pV5Z level of theory. The title reaction is found to be exothermic by 5.63 kcal/mol (9 kcal/mol with zero point energy corrections), which is very close to the experimental data. The barrier height, which is 2.99 kcal/mol (0.93 kcal/mol for the vibrationally adiabatic barrier height), and the depth of van der Waals minimum located near the entrance channel are also in excellent agreement with the experimental findings. This study also identified two more van der Waals minima. The integral cross sections, rate constants, and their dependence on initial rotational states are calculated using an exact quantum wave packet method on the new PES. They are also in excellent agreement with the experimental measurements.

QUASICLASSICAL TRAJECTORY STUDY OF STEREODYNAMICS FOR THE REACTIONS Li+HF/DF/TF

Stereodynamics of the reaction Li + HF (v = 0,j = 0) → LiF + H and its isotopic variants on the ground electronic state (12A′) potential energy surface (PES) are studied by employing the quasiclassical trajectory (QCT) method. At a collision energy of 2.2 kcal/mol, product rotational angular momentum distributions, P(θr) and P(ϕr), are calculated in the center-of-mass (CM) frame. The results demonstrate that the product rotational angular momentum j′ is not only aligned along the direction perpendicular to the reagent relative velocity vector k, but also oriented along the negative y-axis. The four generalized polarization-dependent differential cross sections (PDDCSs) are also computed. The PDDCS00 distribution shows a sideways scattering for the reaction Li + HF and a strongly backward scattering for the reaction Li + DF . However, it displays both the forward and backward scatterings for the reaction Li + TF . These features demonstrate that the Li + HF and Li + DF reactions proceed predominantly through the direct reaction mechanism. However, the Li + TF reaction undergoes both the direct and indirect reaction mechanisms. The PDDCS21- distribution indicates that the product angular distributions are anisotropic.

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.

Reactivity enhancement of ultracold O(P3)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.

Computing resonance energies, widths, and wave functions using a Lanczos method in real arithmetic

We introduce new ideas for calculating resonance energies and widths. It is shown that a non-Hermitian-Lanczos approach can be used to compute eigenvalues of H+W, where H is the Hamiltonian and W is a complex absorbing potential (CAP), without evaluating complex matrix-vector products. This is done by exploiting the link between a CAP-modified Hamiltonian matrix and a real but nonsymmetric matrix U suggested by Mandelshtam and Neumaier [J. Theor. Comput. Chem. 1, 1 (2002)] and using a coupled-two-term Lanczos procedure. We use approximate resonance eigenvectors obtained from the non-Hermitian-Lanczos algorithm and a very good CAP to obtain very accurate energies and widths without solving eigenvalue problems for many values of the CAP strength parameter and searching for cusps. The method is applied to the resonances of HCO. We compare properties of the method with those of established approaches.

Heavy atom tunneling in chemical reactions: Study of H+LiF collisions

The H+LiF(X (1)sigma(+),upsilon=0-2,j=0)-->HF(X (1)sigma(+),upsilon',j')+Li(2S) bimolecular process is investigated by means of quantum scattering calculations on the chemically accurate X 2A' LiHF potential energy surface of Aguado et al. [A. Aguado, M. Paniagua, C. Sanz, and J. Roncero, J. Chem. Phys. 119, 10088 (2003)]. Calculations have been performed for zero total angular momentum for translational energies from 10(-7) to 10(-1) eV. Initial-state selected reaction probabilities and cross sections are characterized by resonances originating from the decay of metastable states of the H...F-Li and Li...F-H van der Waals complexes. Extensive assignment of the resonances has been carried out by performing quasibound states calculations in the entrance and exit channel wells. Chemical reactivity is found to be significantly enhanced by vibrational excitation at low temperatures, although reactivity appears much less favorable than nonreactive processes due to the inefficient tunneling of the relatively heavy fluorine atom strongly bound in van der Waals complexes.

Quantum dynamics of the Li+HF→H+LiF reaction at ultralow temperatures

Quantum-mechanical calculations are reported for the Li+HF(v=0,1,j=0)-->H+LiF(v',j') bimolecular scattering process at low and ultralow temperatures. Calculations have been performed for zero total angular momentum using a recent high-accuracy potential-energy surface for the X2A' electronic ground state. For Li+HF(v=0,j=0), the reaction is dominated by resonances due to the decay of metastable states of the Li cdots,...F-H van der Waals complex. Assignment of these resonances has been carried out by calculating the eigenenergies of the quasibound states. We also find that while chemical reactivity is greatly enhanced by vibrational excitation, the resonances get mostly washed out in the reaction of vibrationally excited HF with Li atoms. In addition, we find that at low energies, the reaction is significantly suppressed due to the less-efficient tunneling of the relatively heavy fluorine atom.

Lifetime of reactive scattering resonances: Q-matrix analysis and angular momentum dependence for the F+H2 reaction by the hyperquantization algorithm

We report a study on the behavior with total angular momentum J of several resonances occurring at collision energies below or slightly above the reaction barrier in the F+H2-->HF+H reaction. Resonance positions and widths are extracted from exact time-independent quantum mechanical calculations using the hyperquantization algorithm and Smith's Q-matrix formalism which exploits complete S-matrix information. The results confirm previous work but provide much greater insight. Identification of quasi-bound states responsible for the resonances based on adiabatic models for the long-range atom-molecule interactions both in the entrance and exit channels, is successful except for the feature occurring at the lowest energy, which is found to overlap with an exit-channel resonance for J approximately 7. The two features are analyzed as overlapping resonances and their excellent Lorentzian fits, with well-behaved J-dependences of positions and widths, support the interpretation of the low-energy feature as a resonance to be associated to the triatomic transition state of the reaction. Resonance role on the reactive observables (integral cross sections and angular distributions) is investigated. The mechanism leading to forward scattering in the reactive differential cross section is commented, while the effects on rate constants, as well as the sensitivity of the resonance pattern to modification of the potential energy surface, are fully discussed elsewhere.

A reduced dimensionality quasiclassical and quantum study of the proton transfer reaction H3O++H2O→H2O+H3O+

We report quantum and quasiclassical calculations of proton transfer in the reaction H(3)O(+)+H(2)O in three degrees of freedom, the two OH(+) bond lengths and the OH(+)O angle. The reduced dimensional potential energy surface is obtained from the full dimensional OSS3(p) energy function of H(5)O(2) (+) [L. Ojamae, I. Shavitt, and S. J. Singer, J. Chem. Phys. 109, 5547 (1998)], with an additional long-range correction to reproduce the correct ion-molecule interaction. This surface is used to perform both quasiclassical trajectory and quantum reactive scattering calculations of the zero total angular momentum cumulative reaction probability and cross sections for initial rotational states 0, 1, and 2. Comparison of these quantities are made to assess the importance of quantum effects in this reduced dimensional reaction. Additional quasiclassical cross sections are calculated to obtain the thermal rate constant for the reaction.