Experimental and Theoretical Investigations of the Radical-Radical Reaction: N2H3 + NO2

The N₂H₃ + NO₂ reaction plays a key role during the early stages of hypergolic ignition between N₂H₄ and N₂O₄. Here for the first time, the reaction kinetics of N₂H₃ in excess NO₂ was studied in 2.0 Torr of N₂ and in the narrow temperature range 298-348 K in a pulsed photolysis flow-tube reactor coupled to a mass spectrometer. The temporal profile of the product, HONO, was determined by direct detection of the m/z +47 amu ion signal. For each chosen [NO₂], the observed [HONO] trace was fitted to a biexponential kinetics expression, which yielded a value for the pseudo-first-order rate coefficient, k', for the reaction of N₂H₃ with NO₂. The slope of the plot of k' versus [NO₂] yielded a value for the observed bimolecular rate coefficient, k_obs, which could be fitted to an Arrhenius expression of (2.36 ± 0.47) × 10^-12 exp((520 ± 350)/T) cm³ molecule^-1 s^-1. The errors are 1σ and include estimated uncertainties in the NO₂ concentration. The potential energy surface of N₂H₃ + NO₂ was investigated by advanced ab initio quantum chemistry theories. It was found that the reaction occurs via a complex reaction mechanism, and all of the reaction channels have transition state energies below that of the entrance asymptote. The radical-radical addition forms the N₂H₃NO₂ adducts, while roaming-mediated isomerization reactions yield the N₂H₃ONO isomers, which undergo rapid dissociation reactions to several sets of distinct products. The RRKM multiwell master equation simulations revealed that the major product channel involves the formation of trans-HONO and trans-N₂H₂ below 500 K and the formation of NO + NH₂NHO above 500 K, which is nearly pressure independent. The pressure-dependent rate coefficients of the product channels were computed over a wide pressure-temperature range, which encompassed the experimental data.

First-Principles Molecular Dynamics of Monomethylhydrazine and Nitrogen Dioxide

The exploration of chemical reactions preceding ignition is essential for the development of ideal hypergolic propellants. Unexpected reaction pathways of a hypergolic mixture composed of monomethylhydrazine and nitrogen dioxide are predicted through a cooperative combination of (i) spin-unrestricted ab initio molecular dynamics (AIMD) and (ii) wave packet dynamics of protons. Ensembles of AIMD trajectories reveal a sequence of reaction steps for proton transfer and rupture of the C-N bond. The details of proton transfer are explored by wave packet dynamics on the basis of ab initio potential energy surfaces from AIMD trajectories. The possibility of spontaneous ignition of this hypergolic mixture at room temperature is predicted as a quantized feature of proton-transfer dynamics.

A Database of Formation Enthalpies of Nitrogen Species by Compound Methods (CBS-QB3, CBS-APNO, G3, G4)

Accurate thermochemical data for compounds containing C/H/N/O are required to underpin kinetics simulation and modeling of the reactions of these species in different environments. There is a dearth of experimental data so computational quantum chemistry has stepped in to fill this breach and to verify whether particular experiments are in need of revision. A number of composite model chemistries (CBS-QB3, CBS-APNO, G3, and G4) are used to compute theoretical atomization energies and hence enthalpies of formation at 0 and 298.15 K, and these are benchmarked against the best available compendium of values, the Active Thermochemical Tables or ATcT. In general the agreement is very good for some 28 species with the only discrepancy being for hydrazine. It is shown that, although individually the methods do not perform that well, collectively the mean unsigned error is <1.7 kJ mol(-1); hence, this approach provides a useful tool to screen published values and validate new experimental results. Using multiple model chemistries does have some drawbacks but can produce good results even for challenging molecules like HOON and CN2O2. The results for these smaller validated molecules are then used as anchors for determining the formation enthalpies of larger species such as methylated hydrazines and diazenes, five- and six-membered heterocyclics via carefully chosen isodesmic working reactions with the aim of resolving some discrepancies in the literature and establishing a properly validated database. This expanded database could be useful in testing the performance of computationally less-demanding density function methods with newer functionals that have the capacity to treat much larger systems than those tested here.

Computational Study of the Rate Coefficients for the Reactions of NO2 with CH3NHNH, CH3NNH2, and CH2NHNH2

The reactions of NO2 with cis-/trans-CH3NHNH, CH3NNH2 and CH2NHNH2 have been studied theoretically by quantum chemical calculations and steady-state unimolecular master equation analysis based on RRKM theory. The barrier heights for the roaming transition states between nitro (RNO2) and nitrite (RONO) isomerization reactions and those for the concerted HONO and HNO2 elimination reactions from RNO2 and RONO, affect the pressure dependences of the product-specific rate coefficients. At ambient temperature and pressure, the dominant product of the reactions of NO2 with cis-/trans-CH3NHNH and CH2NHNH2 would be expected to be HONO with trans-CH3NNH and CH2NNH2, respectively, whereas it is CH3N(NH2)NO2 for CH3NNH2 + NO2. The product-specific rate coefficients for the titled and related reactions on the same potential energy surfaces were proposed for kinetics modeling.