Nonadiabatic energy transfer studies of O(1D)+N2(X 1Σg+)→O(3P)+N2(X 1Σg+) by time-dependent wave packet

Coupled 3D (J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

We had calculated adiabatic potential energy surfaces (PESs), nonadiabatic, and spin-orbit (SO) coupling terms among the lowest three electronic states (12A', 22A', and 12A″) of the F + H2 system using the multireference configuration interaction (MRCI) level of theory, and the adiabatic-to-diabatic transformation equations were solved to formulate the diabatic Hamiltonian matrix [J. Chem. Phys. 2020, 153, 174301] for the entire region of the nuclear configuration space. The accuracy of such diabatic PESs is explored by performing scattering calculations to evaluate integral cross sections (ICSs) and rate constants. The nonadiabatic and SO effects are studied by utilizing coupled 3D time-dependent wave packet formalism with zero and nonzero total angular momentum on multiple adiabatic/diabatic surfaces calculation. We depict the convergence profiles of reaction probabilities for the reactive as well as nonreactive processes on various electronic states at different collision energies with respect to total angular momentum including all helicity quantum numbers. Finally, total ICSs are calculated as functions of collision energies for the initial rovibrational state (v = 0, j = 0) of the H2 molecule along with the temperature-dependent rate coefficient, where those quantities are compared with previous theoretical and experimental results.

Accurate Calculation of Rate Constant and Isotope Effect for the F + H2 Reaction by the Coupled 3D Time-Dependent Wave Packet Method on the Newly Constructed Ab Initio Ground Potential Energy Surface

We employ coupled three-dimensional (3D) time dependent wave packet formalism in hyperspherical coordinates for reactive scattering problem on the newly constructed ab initio calculated ground adiabatic potential energy surface for the F + H₂/D₂ reaction. The convergence profiles for various reactive channels are depicted at low collision energy regimes with respect to the total angular momentum (J) quantum numbers. For two different reactant diatomic molecules (H₂ and D₂) initially at their respective ground roto-vibrational state (v = 0, j = 0), calculated state-to-state as well as total integral cross sections as a function of collision energy, temperature dependent rate constants, and the kinetic isotope effect for various reactivity profiles of F + H₂ and F + D₂ reactions are presented along with previous theoretical and experimental results.

Dynamical calculations of O(3P) + OH(2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates

We have carried out quantum dynamics calculations for the O + OH → H + O₂ reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO₂ using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O₂ → O + OH) and backward (O + OH → H + O₂) reactions of this reactive system, the equilibrium constant of the reversible process [H + O₂ ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.

Charge Transfer Processes for H + H2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces

The dynamics of the H + H₂⁺ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (2¹A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H₂⁺ (²Σ_g⁺) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.

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.

Fully coupled (J > 0) time-dependent wave-packet calculations using hyperspherical coordinates for the H + O2 reaction on the CHIPR potential energy surface

ICS calculation by time dependent wavepacket approach for H + O₂ reaction using non-zero J values.

First-principles dynamics of collisional intersystem crossing: resonance enhanced quenching of C(1D) by N2

Intersystem crossing is a common and important nonadiabatic process in molecular systems. Revealed with a first-principles investigation, the quenching of C(¹D) by N₂ is efficient due to multiple passages via long-lived collisional resonances, despite relatively small spin–orbit couplings.

3D time-dependent wave-packet approach in hyperspherical coordinates for the H + O2 reaction on the CHIPR and DMBE IV potential energy surfaces

3D wavepacket quantum dynamics methodology ICS calculation of H + O₂ reaction on the CHIPR and DMBE IV PESs by J-shifting scheme.

The photodissociation dynamics of OCS at 248nm: The S(PJ3) atomic angular momentum polarization

The dissociation of OCS has been investigated subsequent to excitation at 248 nm using velocity map ion imaging. Speed distributions, speed dependent translational anisotropy parameters, and the atomic angular momentum orientation and alignment are reported for the channel leading to S((3)P(J)). The speed distributions and beta parameters are in broad agreement with previous work and show behavior that is highly sensitive to the S-atom spin-orbit state. The data are shown to be consistent with the operation of at least two triplet production mechanisms. Interpretation of the angular momentum polarization data in terms of an adiabatic picture has been used to help identify a likely dissociation pathway for the majority of the S((3)P(J)) products, which strongly favors production of J=2 fragment atoms, correlated, it is proposed, with rotationally hot and vibrationally cold CO cofragments. For these fragments, optical excitation to the 2 (1)A(') surface is thought to constitute the first step, as for the singlet dissociation channel. This is followed by crossing, via a conical intersection, to the ground 1 (1)A(') state, from where intersystem crossing occurs, populating the 1 (3)A(')1 (3)A(")((3)Pi) states. The proposed mechanism provides a qualitative rationale for the observed spin-orbit populations, as well as the S((3)P(J)) quantum yield and angular momentum polarization. At least one other production mechanism, leading to a more statistical S-atom spin-orbit state distribution and rotationally cold, vibrationally hot CO cofragments, is thought to involve direct excitation to either the (3)Sigma(-) or (3)Pi states.

Spin-orbit effect in the energy pooling reaction O2(aΔ1)+O2(aΔ1)→O2(bΣ1)+O2(XΣ3)

Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.

Time-Dependent Wave Packet Studies on the Cl + HCl Hydrogen Exchange Reaction

The initiation of the hydrogen exchange reaction Cl((2)P)+HCl --> ClH+Cl((2)P) by excitation of the HCl molecular stretch to v=2 is studied for total angular momentum quantum number J=(1)/(2) and both even and odd parity. The calculations were performed using a time-dependent propagation from an initial quasi-bound state and employed all three relevant potential energy surfaces and the nonadiabatic couplings between them. Coriolis and spin-orbit coupling were also taken into account. The electronic and HCl rotational distributions of the products in both dissociation channels are analyzed, and the results are interpreted using features of the potential energy surfaces.