Quasiclassical Trajectory Study of Atom-Exchange and Vibrational Relaxation Processes in Collisions of Atomic and Molecular Nitrogen

Formation Mechanisms of Electronically Excited Nitrogen Molecules from N + N2 and N + N + N Collisions Revealed by Full-Dimensional Potential Energy Surfaces

This work reports six new full-dimensional adiabatic potential energy surfaces (PESs) of the N3 system (four 4A″ states and two 2A″ states) at the MRCI + Q/AVQZ level of theory that correlated to N2(X1Σg+) + N(4S), N2(X1Σg+) + N(2D), N2(A3Σu+) + N(4S), N2(B3Πg) + N(4S), N2(W3Δu) + N(4S), and N(4S) + N(4S) + N(4S) channels. The neural networks with a proper account of the nuclear permutation invariant symmetry of N3 were employed to fit the PESs based on about 4000 ab initio points. The accuracy of the PESs was validated by excellent agreement on the equilibrium bond length, vertical excitation energy, and dissociation energy with experimental values. Two possible mechanisms of the formation of N2(A) were found. One is that the collision occurs between N2(X) and N(4S) in the 14A″ state, followed by a nonadiabatic transition through the conical intersection with the 24A″ PES, resulting in the formation of the N2(A) + N(4S) product. The other takes place in the collision among three N(4S) atoms in the adiabatic 24A″ state, and then, N2(A) + N(4S) is formed. This is the first systematical research of the N3 system focusing on the formation of the excited states of N2 via both adiabatic and nonadiabatic pathways.

High-Energy Quantum Dynamics of the 15N + o-14N14N Rovibrational Activation and Isotope Exchange Processes

We report full quantum reaction probabilities, computed within the framework of time-independent quantum mechanics using hyperspherical coordinates, for the 15N + 14N14N inelastic and reactive collision processes, restricted to total angular momentum J = 0, for kinetic energies up to 4.5 eV. We take advantage of the nonzero (i = 1) nuclear spin of 14N, leading to the existence of two nuclear spin isomers of 14N14N, namely, ortho- and para-14N14N, to restrict the study to the ortho molecular nitrogen species, with even rotational quantum number j = 0, 2, ... states. Specifically, we start with diatomic reagents ortho-14N14N in the initial rotational state j = 0. A comparison with similar works previously published by other groups using time-dependent wave packet and quasi-classical trajectory methods for the 14N + 14N14N fully symmetric collision is given. We find that reactive processes 15N + 14N14N involving atom exchange do not happen for collision energies less than 2.2 eV. Collisions at energies of around 2.0 eV are most effective for populating reactants' rovibrational states, that is, for inelastic scattering, whereas those at energies close to 5.0 eV yield a newly formed 14N15N isotopologue in a wide variety of excited vibrational levels.

Quantum and semiclassical studies of nonadiabatic electronic transitions between N(4S) and N(2D) by collisions with N2

The dynamics and kinetics of spin-forbidden transitions between N(2D) and N(4S) via collisions with N2 molecules are investigated using a quantum wave packet (WP) method and the semi-classical coherent switches with decay of mixing (CSDM) method. These electronic transition processes are competing with exchange reaction channels on both the doublet and quartet potential energy surfaces. The WP and CSDM quenching rate coefficients are found in reasonable agreement with each other, and both reproduce the previous theoretical results. For the excitation process, the agreement between the two approaches is dependent on the treatment of the zero-point energy (ZPE) in the product, because the high endoergicity of this process leads to severe violation of the vibrational ZPE. The Gaussian-binning (GB) method is found to improve the agreement with the quantum result. The excitation rate coefficients are found to be two orders of magnitude smaller than that of the adiabatic exchange reaction, underscoring the inefficient intersystem crossing due to the weak spin-orbit coupling between the two spin manifolds of the N3 system.

Quasiclassical Trajectory Study of the Si + SH Reaction on an Accurate Double Many-Body Expansion Potential Energy Surface

An accurate potential energy surface (PES) for the HSiS system based on MRCI+Q calculations extrapolated to the complete basis set limit is presented. Modeled with the double many-body expansion (DMBE) method, the PES provides an accurate description of the long-range interactions, including electrostatic and dispersion terms decaying as R^-4, R^-5, R^-6, R^-8, R^-10 that are predicted from dipole moments, quadrupole moments, and dipolar polarizabilities, which are also calculated at the MRCI+Q level. The novel PES is then used in quasiclassical trajectory calculations to predict the rate coefficients of the Si + SH → SiS + H reaction, which has been shown to be a major source of the SiS in certain regions of the interstellar medium. An account of the zero-point energy leakage based on various nonactive models is also given. It is shown that the reaction is dominated by long-range forces, with the mechanism Si + SH → SiSH → SSiH → SiS + H being the most important one for all temperatures studied. Although SSiH corresponds to the global minimum of the PES, the contribution from the direct reaction Si + SH → SSiH → SiS + H is less than 0.5% for temperatures higher than 500 K. The rovibrational distributions of the products are calculated using the momentum Gaussian binning method and show that as the temperature is increased the average vibrational quantum number decreases while the rotational distribution spreads up to larger values.

Nonequilibrium internal energy distributions during dissociation

In this work, we propose a model for nonequilibrium vibrational and rotational energy distributions in nitrogen using surprisal analysis. The model is constructed by using data from direct molecular simulations (DMSs) of rapidly heated nitrogen gas using an ab initio potential energy surface (PES). The surprisal-based model is able to capture the overpopulation of high internal energy levels during the excitation phase and also the depletion of high internal energy levels during the quasi-steady-state (QSS) dissociation phase. Due to strong coupling between internal energy and dissociation chemistry, such non-Boltzmann effects can influence the overall dissociation rate in the gas. Conditions representative of the flow behind strong shockwaves, relevant to hypersonic flight, are analyzed. The surprisal-based model captures important molecular-level nonequilibrium physics, yet the simple functional form leads to a continuum-level expression that now accounts for the underlying energy distributions and their coupling to dissociation.

Modeling cusps in adiabatic potential energy surfaces

A method for modeling cusps on adiabatic potential energy surfaces without the need for any adiabatic-to-diabatic transformation is presented and shown to be successfully applied to the (2)A″ state of NO2. The more complicated case of a system with permutationally equivalent crossing seams is also examined and illustrated by considering the two first (2)A' states of the nitrogen trimer.

Electronic Quenching in N((2)D) + N2 Collisions: A State-Specific Analysis via Surface Hopping Dynamics

The electronic quenching reaction N((2)D) + N2 → N((4)S) + N2 is studied using the trajectory surface hopping method and employing two doublet and one quartet accurate potential energy surfaces. State-specific properties are analyzed, such as the dependence of the cross section on the initial quantum state of the reactants, vibrational energy transfer, and rovibrational distribution of the product N2 molecule in thermalized conditions. It is found that rotational energy on the reactant N2 molecule is effective in promoting the reaction, whereas vibrational excitation tends to reduce the reaction probability. For initial states and collision energy thermalized in an initial bath, it is found that the products are "hotter", both vibration and rotation wise.

A high-level ab initio study of the N2 + N2 reaction channel

A new six-dimensional (6D) global potential energy surface (PES) is proposed for the full range description of the interaction of the N2(1Σg+)+N2(1Σg+) system governing collisional processes, including N atom exchange. The related potential energy values were determined using high-level ab initio methods. The calculations were performed at a coupled-cluster with single and double and perturbative triple excitations level of theory in order to have a first full range picture of the PES. Subsequently, in order to accurately describe the stretching of the bonds of the two interacting N2 molecules by releasing the constraints of being considered as rigid rotors, for the same molecular geometries higher level of theory multi reference calculations were performed. Out of the calculated values a 6D 4-atoms global PES was produced for use in dynamical calculations. The ab initio calculations were made possible by the combined use of High Throughput Computing and High Performance Computing techniques within the frame of a computing grid empowered molecular simulator.