Low-temperature kinetics for the N + NO reaction: experiment guides the way.

The reaction N(4S) + NO(X2Π) → O(3P) + N2(X1Σ+g) plays a pivotal role in the conversion of atomic to molecular nitrogen in dense interstellar clouds and in the atmosphere. Here we report a joint experimental and computational investigation of the N + NO reaction with the aim of providing improved constraints on its low temperature reactivity. Thermal rates were measured over the 50 to 296 K range in a continuous supersonic flow reactor coupled with pulsed laser photolysis and laser induced fluorescence for the production and detection of N(4S) atoms, respectively. With decreasing temperature, the experimentally measured reaction rate was found to monotonously increase up to a value of (6.6 ± 1.3) × 10-11 cm3 s-1 at 50 K. To confirm this finding, quasi-classical trajectory simulations were carried out on a previously validated, full-dimensional potential energy surface (PES). However, around 50 K the computed rates decreased which required re-evaluation of the reactive PES in the long-range part due to a small spurious barrier with a height of ∼40 K in the entrance channel. By exploring different correction schemes the measured thermal rates can be adequately reproduced, displaying a clear negative temperature dependence over the entire temperature range. The possible astrochemical implications of an increased reaction rate at low temperature are also discussed.

Experimental and Theoretical Studies of Hyperthermal N + O2 Collisions.

The dynamics of hyperthermal N(4S) + O2 collisions were investigated both experimentally and theoretically. Crossed molecular beams experiments were performed at an average center-of-mass (c.m.) collision energy of ⟨Ecoll⟩ = 77.5 kcal mol-1, with velocity- and angle-resolved product detection by a rotatable mass spectrometer detector. Nonreactive (N + O2) and reactive (NO + O) product channels were identified. In the c.m. reference frame, the nonreactively scattered N atoms and reactively scattered NO molecules were both directed into the forward direction with respect to the initial direction of the reagent N atoms. On average, more than 90% of the available energy (⟨Eavl⟩ = 77.5 kcal mol-1) was retained in translation of the nonreactive products (N + O2), whereas a much smaller fraction of the available energy for the reactive pathway (⟨Eavl⟩ = 109.5 kcal mol-1) went into translation of the NO + O products, and the distribution of translational energies for this channel was broad, indicating extensive internal excitation in the nascent NO molecules. The experimentally derived c.m. translational energy and angular distributions of the reactive products suggested at least two dynamical pathways to the formation of NO + O. Quasiclassical trajectory (QCT) calculations were performed with a collision energy of Ecoll = 77 kcal mol-1 using two sets of potential energy surfaces, denoted as PES-I and PES-II, and these theoretical results were compared to each other and to the experimental results. PES-I is a reproducing kernel Hilbert space (RKHS) representation of multireference configurational interaction (MRCI) energies, while PES-II is a many-body permutation invariant polynomial (MB-PIP) fit of complete active space second order perturbation (CASPT2) points. The theoretical investigations were both consistent with the experimental suggestion of two dynamical pathways to produce NO + O, where reactive collisions may proceed on the doublet (12A') and quartet (14A') surfaces. When analyzed with this theoretical insight, the experimental c.m. translational energy and angular distributions were in reasonably good agreement with those predicted by the QCT calculations, although minor differences were observed which are discussed. Theoretical translational energy and angular distributions for the nonreactive N + O2 products matched the experimental translational energy and angular distributions almost quantitatively. Finally, relative yields for the nonreactive and reactive scattering channels were determined from the experiment and from both theoretical methods, and all results are in reasonable agreement.

Quantum and quasi-classical dynamics of the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction on its electronic ground state.

The dynamics of the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (Ec = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO2 complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v' ≥ 15) and for the lowest product rotational states (j' ≤ 10). Opening of the product vibrational level CO(v' = 17) requires ∼0.2 eV from QCT and TDWP simulations with O2(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with Ec > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs.

Combining Machine Learning and Spectroscopy to Model Reactive Atom + Diatom Collisions.

The prediction of product translational, vibrational, and rotational energy distributions for arbitrary initial conditions for reactive atom + diatom collisions is of considerable practical interest in atmospheric re-entry. Because of the large number of accessible states, determination of the necessary information from explicit (quasi-classical or quantum) dynamics studies is impractical. Here, a machine-learned (ML) model based on translational energy and product vibrational states assigned from a spectroscopic, ro-vibrational coupled energy expression based on the Dunham expansion is developed and tested quantitatively. All models considered in this work reproduce final state distributions determined from quasi-classical trajectory (QCT) simulations with R2 ∼ 0.98. As a further validation, thermal rates determined from the machine-learned models agree with those from explicit QCT simulations and demonstrate that the atomistic details are retained by the machine learning which makes them suitable for applications in more coarse-grained simulations. More generally, it is found that ML is suitable for designing robust and accurate models from mixed computational/experimental data which may also be of interest in other areas of the physical sciences.

Machine learning product state distributions from initial reactant states for a reactive atom-diatom collision system.

A machine-learned model for predicting product state distributions from specific initial states (state-to-distribution or STD) for reactive atom-diatom collisions is presented and quantitatively tested for the N(4S) + O2(X3Σg -) → NO(X2Π) + O(3P) reaction. The reference dataset for training the neural network consists of final state distributions determined from quasi-classical trajectory (QCT) simulations for ∼2000 initial conditions. Overall, the prediction accuracy as quantified by the root-mean-squared difference (∼0.003) and the R2 (∼0.99) between the reference QCT and predictions of the STD model is high for the test set, for off-grid state-specific initial conditions, and for initial conditions drawn from reactant state distributions characterized by translational, rotational, and vibrational temperatures. Compared with a more coarse grained distribution-to-distribution (DTD) model evaluated on the same initial state distributions, the STD model shows comparable performance with the additional benefit of the state resolution in the reactant preparation. Starting from specific initial states also leads to a more diverse range of final state distributions, which requires a more expressive neural network compared with DTD. A direct comparison between QCT simulations, the STD model, and the widely used Larsen-Borgnakke (LB) model shows that the STD model is quantitative, whereas the LB model is qualitative at best for rotational distributions P(j') and fails for vibrational distributions P(v'). As such, the STD model can be well-suited for simulating nonequilibrium high-speed flows, e.g., using the direct simulation Monte Carlo method.

Photodissociation dynamics of N3.

The photodissociation dynamics of N3 + excited from its (linear) 3Σg -/(bent) 3A″ ground to the first excited singlet and triplet states is investigated. Three-dimensional potential energy surfaces for the 1A', 1A″, and 3A' electronic states, correlating with the 1Δg and 3Πu states in linear geometry, for N3 + are constructed using high-level electronic structure calculations and represented as reproducing kernels. The reference ab initio energies are calculated at the MRCI+Q/aug-cc-pVTZ level of theory. For following the photodissociation dynamics in the excited states, rotational and vibrational distributions P(v') and P(j') for the N2 product are determined from vertically excited ground state distributions. Due to the different shapes of the ground state 3A″ potential energy surface and the excited states, appreciable angular momentum j' ∼ 60 is generated in diatomic fragments. The lifetimes in the excited states extend to at least 50 ps. Notably, results from sampling initial conditions from a thermal ensemble and from the Wigner distribution of the ground state wavefunction are comparable.

The C(3P) + O2(3Σg-) → CO2 ↔ CO(1Σ+) + O(1D)/O(3P) reaction: thermal and vibrational relaxation rates from 15 K to 20 000 K.

Thermal rates for the C(3P) + O2(3Σg-) ↔ CO(1Σ+)+ O(1D)/O(3P) reaction are investigated over a wide temperature range based on quasi classical trajectory (QCT) simulations on 3-dimensional, reactive potential energy surfaces (PESs) for the 1A', (2)1A', 1A'', 3A' and 3A'' states. These five states are the energetically low-lying states of CO2 and their PESs are computed at the MRCISD+Q/aug-cc-pVTZ level of theory using a state-average CASSCF reference wave function. Analysis of the different electronic states for the CO2 → CO + O dissociation channel rationalizes the topography of this region of the PESs. The forward rates from QCT simulations match measurements between 15 K and 295 K whereas the equilibrium constant determined from the forward and reverse rates is consistent with that derived from statistical mechanics at high temperature. Vibrational relaxation, O + CO(ν = 1,2) → O + CO(ν = 0), is found to involve both, non-reactive and reactive processes. The contact time required for vibrational relaxation to take place is τ ≥ 150 fs for non-reacting and τ ≥ 330 fs for reacting (oxygen atom exchange) trajectories and the two processes are shown to probe different parts of the global potential energy surface. In agreement with experiments, low collision energy reactions for the C(3P) + O2(3Σg-, ν = 0) → CO(1Σ+) + O(1D) lead to CO(1Σ+, ν' = 17) with an onset at Ec ∼ 0.15 eV, dominated by the 1A' surface with contributions from the 3A' surface. Finally, the barrier for the COA(1Σ+) + OB(3P) → COB(1Σ+) + OA(3P) atom exchange reaction on the 3A' PES yields a barrier of ∼7 kcal mol-1 (0.300 eV), consistent with an experimentally reported value of 6.9 kcal mol-1 (0.299 eV).

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.

Accurate potential energy surfaces (PESs) have been determined for the 3A' and 3A'' states of N2O using electronic structure calculations at the multireference configuration interaction level with Davidson correction (MRCI+Q) and the augmented Dunning-type correlation consistent polarized triple zeta (aug-cc-pVTZ) basis set. More than 20 000 MRCI+Q/aug-cc-pVTZ energies are represented using a reproducing kernel Hilbert space (RKHS) scheme. The RKHS PESs successfully describe all reactant channels with high accuracy and all minima and transition states connecting them are determined. Quasiclassical trajectory (QCT) simulations are then used to determine reaction rates for N + NO and O + N2 collisions. Vibrational relaxation N2(ν = 1) → N2(ν = 0) and dissociation of N2→ 2N for O + N2 collisions are also investigated using QCT. The agreement between results obtained from the QCT simulations and from available experiments is favourable for reaction and vibrational relaxation rates, which provides a test for the accuracy of the PESs. The PESs can be used to calculate more detailed state-to-state observables relevant for applications to hypersonic reentry.

The N(4S) + O2(X3Σ) ↔ O(3P) + NO(X2Π) reaction: thermal and vibrational relaxation rates for the 2A', 4A' and 2A'' states.

The kinetics and vibrational relaxation of the N(4S) + O2(X3Σ-g) ↔ O(3P) + NO(X2Π) reaction is investigated over a wide temperature range based on quasiclassical trajectory simulations on 3-dimensional potential energy surfaces (PESs) for the lowest three electronic states. Reference energies at the multi reference configuration interaction level are represented as a reproducing kernel and the topology of the PESs is rationalized by analyzing the CASSCF wavefunction of the relevant states. The forward rate matches one measurement at 1575 K and is somewhat lower than the high-temperature measurement at 2880 K whereas for the reverse rate the computations are in good agreement for temperatures between 3000 and 4100 K. The temperature-dependent equilibrium rates are consistent with results from JANAF and CEA results. Vibrational relaxation rates for O + NO(ν = 1) → O + NO(ν = 0) are consistent with a wide range of experiments. This process is dominated by the dynamics on the 2A' and 4A' surfaces which both contribute similarly up to temperatures T ∼ 3000 K, and it is found that vibrationally relaxing and non-relaxing trajectories probe different parts of the potential energy surface. The total cross section depending on the final vibrational state monotonically decreases which is consistent with early experiments and previous simulations but at variance with other recent experiments which reported an oscillatory cross section.

Near dissociation states for H2+-He on MRCI and FCI potential energy surfaces.

New potential energy surfaces (PES) have been constructed for H2+-He using a reproducing kernel Hilbert space (RKHS) representation from an extensive number of ab initio energies computed at the multi reference and full configuration interaction levels of theory. For the MRCI PES the long-range interaction region of the PES is described by analytical functions and is connected smoothly to the short range interaction region, represented as a RKHS. All bound ro-vibrational states for the ground electronic state of H2+-He are calculated using two different methods to determine quantum bound states. Comparing transition frequencies for the near-dissociation states for ortho- and para-H2+-He allows assignment of the 15.2 GHz line to a J = 2 e/f parity doublet of ortho-H2+-He whereas the experimentally determined 21.8 GHz line is only consistent with a (J = 0) → (J = 1) e/e transition in para-H2+-He.