Ab initio derived analytical fits of the two lowest triplet potential energy surfaces and theoretical rate constants for the N(4S)+NO(X 2Π) system

Rovibrational-Specific QCT and Master Equation Study on N2(X1Σg+) + O(3P) and NO(X2Π) + N(4S) Systems in High-Energy Collisions

This work presents a detailed investigation of the energy-transfer and dissociation mechanisms in N₂(X¹Σ_g⁺) + O(³P) and NO(X²Π) + N(⁴S) systems using rovibrational-specific quasiclassical trajectory (QCT) and master equation analyses. The complete set of state-to-state kinetic data, obtained via QCT, allows for an in-depth investigation of the Zel'dovich mechanism leading to the formation of NO molecules at microscopic and macroscopic scales. The master equation analysis demonstrates that the low-lying vibrational states of N₂ and NO have dominant contributions to the NO formation and the corresponding extinction of N₂ through the exchange process. For the considered temperature range, it is found that nearly 50% of the dissociation processes for N₂ and NO molecules occur in the quasi-steady-state (QSS) regime, while for the Zel'dovich reaction, the distribution of the reactants does not reach the QSS conditions. Furthermore, using the QSS approximation to model the Zel'dovich mechanism leads to overestimating NO production by more than a factor of 4 in the high-temperature range. The breakdown of this well-known approximation has profound consequences for the approaches that heavily rely on the validity of QSS assumption in hypersonic applications. Finally, the investigation of the rovibrational state population dynamics reveals substantial similarities among different chemical systems for the energy-transfer and the dissociation processes, providing promising physical foundations for the use of reduced-order strategies in other chemical systems without significant loss of accuracy.

An Experimental and Theoretical Investigation of the Gas-Phase C(3P) + N2O Reaction. Low Temperature Rate Constants and Astrochemical Implications

The reaction between atomic carbon in its ground electronic state, C(³P), and nitrous oxide, N₂O, has been studied below room temperature due to its potential importance for astrochemistry, with both species considered to be present at high abundance levels in a range of interstellar environments. On the experimental side, we measured rate constants for this reaction over the 50-296 K range using a continuous supersonic flow reactor. C(³P) atoms were generated by the pulsed photolysis of carbon tetrabromide at 266 nm and were detected by pulsed laser-induced fluorescence at 115.8 nm. Additional measurements allowing the major product channels to be elucidated were also performed. On the theoretical side, statistical rate theory was used to calculate low temperature rate constants. These calculations employed the results of new electronic structure calculations of the ³A″ potential energy surface of CNNO and provided a basis to extrapolate the measured rate constants to lower temperatures and pressures. The rate constant was found to increase monotonically as the temperature falls (k_C(³P)+N2O (296 K) = (3.4 ± 0.3) × 10^-11 cm³ s^-1), reaching a value of k_C(³P)+N2O (50 K) = (7.9 ± 0.8) × 10^-11 cm³ s^-1 at 50 K. As current astrochemical models do not include the C + N₂O reaction, we tested the influence of this process on interstellar N₂O and other related species using a gas-grain model of dense interstellar clouds. These simulations predict that N₂O abundances decrease significantly at intermediate times (10³ - 10⁵ years) when gas-phase C(³P) abundances are high.

Reconciling experimental and theoretical vibrational deactivation in low-energy O + N2 collisions

Molecular dynamics calculations of inelastic collisions of atomic oxygen with molecular nitrogen are known to show orders of magnitude discrepancies with experimental results in the range from room temperature to many thousands of degrees Kelvin. In this work, we have achieved an unprecedented quantitative agreement with experiments even at low temperature, by including a non-adiabatic treatment involving vibronic states on newly developed potential energy surfaces. This result paves the way for the calculation of accurate and detailed databases of vibrational energy exchange rates for this collisional system. This is bound to have an impact on air plasma simulations under a wide range of conditions and on the development of Very Low Earth Orbit (VLEO) satellites, operating in the low thermosphere, objects of great technological interest due to their potential at a competitive cost.

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 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.

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.

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

The reaction N+NO-->N(2)+O was studied by means of the time-dependent real wave-packet (WP) method and the J-shifting approximation. We consider the ground 1 (3)A(") and first excited 1 (3)A(') triplet states, which correlate with both reactants and products, using analytical potential energy surfaces (PESs) recently developed in our group. This work extends our previous quantum dynamics study, and probabilities, cross sections, and rate constants were calculated and interpreted on the basis of the different shapes of the PESs (barrierless 1 (3)A(") and with barrier 1 (3)A(') surfaces, respectively). The WP rate constant (k(1)) shows a weak dependence on T(200-2500 K), as the dominant contribution to reactivity is provided by the barrierless ground PES. There is a good agreement of WP k(1) with the measurements and variational transition state theory (VTST) data, and also between the WP and VTST k(1)(1 (3)A(")) results. Nevertheless, there is a large discrepancy between the WP and VTST k(1)(1 (3)A(')) results. Product state distributions were also calculated for the much more reactive 1 (3)A(") PES. There is an excellent agreement with the experimental average fraction of vibrational energy in N(2)(25+/-3%), the only measured dynamics property of this reaction.

N(4S) formation following the 193.3-nm ArF laser irradiation of NO and NO2 and its application to kinetic studies of N(4S) reactions with NO and NO2

Formation of the ground-state nitrogen atom, N((4)S), following 193.3-nm ArF laser irradiation of NO and NO(2) was detected directly by a technique of laser-induced fluorescence (LIF) spectroscopy at 120.07 nm. Tunable vacuum ultraviolet (VUV) laser radiation around 120.07 nm was generated by two-photon resonance four-wave sum frequency mixing in Hg vapor. Photoexcitation processes of NO and NO(2) giving rise to the N((4)S) formation are discussed on the basis of the Doppler profiles of the nascent N((4)S) atoms produced from the photolysis of NO and NO(2) and the photolysis laser-power dependence of the N((4)S) signal intensities. Using laser flash photolysis and vacuum ultraviolet laser-induced fluorescence detection, the kinetics of the reactions of N((4)S) with NO and NO(2) have been investigated at 295 +/- 2 K. The rate constants for the reactions of N((4)S) with NO and NO(2) were determined to be (3.8 +/- 0.2) x 10(-11) and (7.3 +/- 0.9) x 10(-12) cm(3) molecule(-1) s(-1), respectively, where the quoted uncertainties are 2sigma statistical uncertainty including estimated systematic error.

Ab initiopotential energy surface, variational transition state theory, and quasiclassical trajectory studies of the F+CH4→HF+CH3 reaction

In this work we present a study of the F+CH(4)-->HF+CH(3) reaction (DeltaHdegrees(298 K)=-32.0 kcal mol(-1)) using different methods of the chemical reaction theory. The ground potential energy surface (PES) is characterized using several ab initio methods. Full-dimensional rate constants have been calculated employing the variational transition state theory and using directly ab initio data. A triatomic analytical representation of the ground PES was derived from ab initio points calculated at the second- and fourth-order Møller-Plesset levels with the 6-311+G(2df,2pd) basis set, assuming the CH(3) fragment to be a 15 a.m.u. pseudoatom in the fitting process. This is suggested from experiments that indicate that the methyl group is uncoupled to the reaction coordinate. A dynamics study by means of the quasiclassical trajectory (QCT) method and employing this analytical surface was also carried out. The experimental data available on the HF internal states distributions are reproduced by the QCT results. Very recent experimental information about the reaction stereodynamics is also borne out by our QCT calculations. Comparisons with the benchmark F+H(2) and analogous Cl+CH(4) reactions are established throughout.

Ab initio analytical potential energy surface and quasiclassical trajectory study of the O+(4S)+H2(X 1Σg+)→OH+(X 3Σ−)+H(2S) reaction and isotopic variants

An analytical potential energy surface (PES) representation of the O(+)((4)S)+H(2)(X (1)Sigma(g) (+)) system was developed by fitting around 600 CCSD(T)/cc-pVQZ ab initio points. Rate constant calculations for this reaction and its isotopic variants (D(2) and HD) were performed using the quasiclassical trajectory (QCT) method, obtaining a good agreement with experimental data. Calculations conducted to determine the cross section of the title reaction, considering collision energies (E(T)) below 0.3 eV, also led to good accord with experiments. This PES appears to be suitable for kinetics and dynamics studies. Moreover, the QCT results show that, although the hypotheses of a widely used capture model are not satisfied, the resulting expression for the cross section can be applied within a suitable E(T) interval, due to errors cancellation. This could be a general situation regarding the application of this simple model to ion-molecule processes.