Lanczos Subspace Time-Independent Wave Packet Calculations of S (1D) + H2 Reactive Scattering

Spin-Dependent Stereochemistry: A Nonadiabatic Quantum Dynamics Case Study of S + H2 → SH + H Reaction

We study the spin-dependent stereodynamics of the S + H2 → SH + H reaction by using full-dimensional quantum dynamics calculations with zero total nuclear angular momentum along the triplet 3A″ states and singlet 1A' states. We find that the interplay between the electronic spin direction and the molecular geometry has a measurable influence on the singlet-triplet intersystem crossing reaction probabilities. Our results show that for some incident scattering angles in the body-fixed frame, the relative difference in intersystem crossing reaction probabilities (as determined between spin up and spin down initial states) can be as large as 15%. Our findings are an ab initio demonstration of spin-dependent nonadiabatic dynamics, which we hope will shine light as far as understanding the chiral-induced spin selectivity effect.

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.

Reduced-Order Modeling, Error Estimation, and the Role of the Start-Vector: The Recursive Residue Generation Method Revisited

The recursive residue generation method (RRGM) [Wyatt, R. E. Adv. Chem. Phys. 1989, 73, 231] is re-derived using the formalism of reduced-order modeling [Bai, Z. Appl. Numerical Math. 2002, 43, 9]. A stopping criteria for the RRGM recursions is proposed, on the basis of an expression for an upper bound to the absolute error [Bai, Z.; Ye, Q. Electron. Trans. Numerical Anal. 1998, 7, 1]. It is further pointed out that, in general, the start-vector has a negligible effect on the convergence of the RRGM.

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.

Determination of the internal state distribution of the SD product from the S(1D)+D2 reaction

The S(1D)+D2-->SD+D reaction has been studied through a photolysis-probe experiment in a cell. S(1D) reagent was prepared by 193 nm photolysis of CS2, and the SD(X 2Pi) product was detected by laser fluorescence excitation. The nascent rotational/fine-structure state distribution of the SD(X 2Pi) product was determined. This reaction, previously studied theoretically and in a crossed molecular beam experiment, is known to proceed through formation and decay of a long-lived collision complex involving the deep well in the H2S ground electronic state. The determined SD rotational state distribution in the v=0 vibrational level was found to be approximately statistical, with a small preference for formation of the F1 (Omega=3/2) fine-structure manifold over F2 (Omega=1/2). The branching into the Lambda doublet levels was also investigated, and essentially equal populations of levels of A' and A" symmetry were found. The present results are compared with previous investigations of this reaction and the analogous O(1D)+D2 reaction.