Quantum real wave-packet dynamics of the N(S4)+NO(X̃Π2)→N2(X̃Σg+1)+O(P3) reaction on the ground and first excited triplet potential energy surfaces: Rate constants, cross sections, and product distributions

Accurate reproducing kernel-based potential energy surfaces for the triplet ground states of N2O and dynamics for the N + NO ↔ O + N2 and N2 + O → 2N + O reactions

Reproducing kernel-based potential energy surface based on MRCI+Q/aug-cc-pVTZ energies for the triplet states of N₂O and quasiclassical dynamical study for the reaction, dissociation and vibrational relaxation.

Reactive, Inelastic, and Dissociation Processes in Collisions of Atomic Oxygen with Molecular Nitrogen

We report the results of detailed calculations of reactive, inelastic, and dissociative processes in collisions of atomic oxygen with molecular nitrogen in their respective electronic ground states. Cross sections are calculated as a function of collision energy in the range 0.001-10 eV, considering the whole rovibrational ladder. Some problems related to the vibrational energy levels of the asymptotes of ³A″ and ³A' potential energy surfaces used in this work are solved by an appropriate scaling at the level of cross sections. The results are compared with data in the literature, obtaining excellent agreement with experimental thermal data for reactive processes on a very large temperature range, and reasonable agreement with indirect dissociative data. Significant discrepancies are observed with previous reactive state-to-state results calculated on less detailed potential energy surfaces. Inelastic results are compatible with extrapolation of experimental thermal rate coefficient for temperatures higher than 4500 K but completely fail to reproduce experimental data at room temperature. The issue is discussed, indicating the reasons and possible solutions to the problem, and a resonable rate coefficient is obtained combining experimental and theoretical results in the range 300-20000 K. Complete, accurate fits are provided for both reactive and dissociative state-to-state rate coefficients to use them in applicative numerical codes concerning air kinetics.

Reactive collisions for NO(2Π) + N(4S) at temperatures relevant to the hypersonic flight regime

Rate coefficients for the NO(²Π) + N(⁴S) reaction at high temperatures from quasiclassical trajectories using MRCI+Q PESs of the lowest triplet states.

Global triplet potential energy surfaces for the N2(X(1)Σ) + O((3)P) → NO(X(2)Π) + N((4)S) reaction

This work presents two global triplet potential energy surfaces (PESs) for the high-energy reaction N2(X(1)Σ) + O((3)P) → NO(X(2)Π) + N((4)S)-in particular, for the lowest energy (3)A' and (3)A″ PESs. In order to obtain the energies needed for fitting analytic surfaces, we carried out multireference configuration interaction (MRCI) calculations based on wave functions obtained from state-averaged complete active space self-consistent field calculations for 2280 geometries for the three lowest (3)A″ states and for 2298 geometries for the three lowest (3)A' states. The lowest-energy (3)A' and (3)A″ states at each of these geometries were then improved by applying the dynamically scaled external correlation (DSEC) method to all MRCI points, and the resulting DSEC energies were used for construction of the ground-state PES for each symmetry. The many-body component of the DSEC energies for the three-dimensional (3)A' and (3)A″ PESs was then least-squares fitted in terms of permutationally invariant polynomials in mixed exponential-Gaussian bond order variables. The global and local minima as well as the transition structures of both the (3)A' and the (3)A″ analytic PES were explored. In agreement with previous work, we find that the reverse reaction is barrierless on the (3)A″ surface along the minimum energy pathway. However, we have explored several new local minima and transition structures on the (3)A' PES. Furthermore, based on the newly found minima and transition structures, two independent reaction mechanisms have been illustrated for the reaction path on the (3)A' PES. The analytic surfaces may be used for dynamics calculations of electronically adiabatic reactive scattering and energy transfer.

State-to-state quantum dynamics of the N(4S) + CH(X2Π) → CN(X2Σ+,A2Π) + H(2S) reactions

The reactions between N((4)S) and CH(X(2)Π) lead to H((2)S) plus CN in its two lowest electronic states (X(2)Σ(+) and A(2)Π), which are responsible for the interstellar CN formation. Accurate quantum dynamics of these reactions are investigated on new global potential energy surfaces of the two lowest-lying triplet states of HCN (1(3)A' and 1(3)A") fitted to more than 37,000 points at the internally contracted multi-reference configuration interaction level with the Davidson correction. The pathways for these highly exothermic and barrierless reactions feature both the HCN and HNC wells. Long-lived resonances supported by these wells manifest in reaction probabilities as numerous oscillations, particularly for low J partial waves. The 1(3)A" state is found to be more reactive than the 1(3)A' state, due apparently to its more attractive nature in the entrance channel. The CN products in both electronic states are highly excited in both vibrational and rotational degrees of freedom. The near forward-backward symmetric differential cross sections are consistent with a complex-forming mechanism.

Coriolis coupling effects in the calculation of state-to-state integral and differential cross sections for the H+D2 reaction

The quantum wavepacket parallel computational code DIFFREALWAVE is used to calculate state-to-state integral and differential cross sections for the title reaction on the BKMP2 surface in the total energy range of 0.4-1.2 eV with D2 initially in its ground vibrational-rotational state. The role of Coriolis couplings in the state-to-state quantum calculations is examined in detail. Comparison of the results from calculations including the full Coriolis coupling and those using the centrifugal sudden approximation demonstrates that both the energy dependence and the angular dependence of the calculated cross sections are extremely sensitive to the Coriolis coupling, thus emphasizing the importance of including it correctly in an accurate state-to-state calculation.

State-to-state reactive differential cross sections for the H+H2→H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.

Time dependent quantum dynamics study of the O++H2(v=0,j=0)→OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.