Thermal Dissociation and Roaming Isomerization of Nitromethane: Experiment and Theory

Area-Selective Deposition by Cyclic Adsorption and Removal of 1-Nitropropane

The ever-greater complexity of modern electronic devices requires a larger chemical toolbox to support their fabrication. Here, we explore the use of 1-nitropropane as a small molecule inhibitor (SMI) for selective atomic layer deposition (ALD) on a combination of SiO2, Cu, CuOx, and Ru substrates. Results using water contact angle goniometry, Auger electron spectroscopy, and infrared spectroscopy show that 1-nitropropane selectively chemisorbs to form a high-quality inhibition layer on Cu and CuOx at an optimized temperature of 100 °C, but not on SiO2 and Ru. When tested against Al2O3 ALD, however, a single pulse of 1-nitropropane is insufficient to block deposition on the Cu surface. Thus, a new multistep process is developed for low-temperature Al2O3 ALD that cycles through exposures of 1-nitropropane, an aluminum metalorganic precursor, and coreactants H2O and O3, allowing the SMI to be sequentially reapplied and etched. Four different Al ALD precursors were investigated: trimethylaluminum (TMA), triethylaluminum (TEA), tris(dimethylamido)aluminum (TDMAA), and dimethylaluminum isopropoxide (DMAI). The resulting area-selective ALD process enables up to 50 cycles of Al2O3 ALD on Ru but not Cu, with 98.7% selectivity using TEA, and up to 70 cycles at 97.4% selectivity using DMAI. This work introduces a new class of SMI for selective ALD at lower temperatures, which could expand selective growth schemes to biological or organic substrates where temperature instability may be a concern.

A Combined Experimental and Modeling Study on Isopropyl Nitrate Pyrolysis

Alkyl nitrates thermally decompose by homolytic cleavage of the weak nitrate bond at very low temperatures (e.g., around 500 K at reaction times of a few seconds). This provides the opportunity to study the subsequent chemistry of the initially formed radical (or its subsequent pyrolysis products, if unstable) and nitrogen dioxide at such mild conditions. In this work this idea is applied to isopropyl nitrate (iPN) pyrolysis, which is studied in a tubular reactor at atmospheric pressure, temperatures ranging from 373 to 773 K, and residence times of around 2 s. At the experimental conditions, iPN decomposition starts at 473 K with O-N bond fission producing isopropoxy radical (i-C₃H₇O) and NO₂. i-C₃H₇O is rapidly converted to acetaldehyde (CH₃CHO), which is the most abundant product detected, and methyl radicals. Other major products detected are formaldehyde (CH₂O), methanol (CH₃OH), nitromethane (CH₃NO₂), NO, methane, formamide (CHONH₂), and methyl nitrite (CH₃ONO). Four literature nitrogen chemistry models─three of those augmented with iPN specific reactions─have been tested for their ability to predict the iPN decomposition and product profiles. The mechanism by the Curran group performs best, but it still underpredicts the observed high formaldehyde and methanol yields. A rate analysis indicates that the branching ratio of the reaction between methyl radicals and nitrogen dioxide is of significant importance. Based on recent theoretical and experimental data, new rate expressions for the two reactions CH₃ + NO₂ → CH₃O + NO and CH₃ + NO₂ + He → CH₃ONO₂ + He are calculated and incorporated in the kinetic models. It is shown that this change clearly improves the predictions, although additional work is needed to achieve good agreement between calculated and measured species profiles.

Unusual In-plane Aromaticity Facilitates Intramolecular Hydrogen Transfer in Long-Bonded cis-Isonitrosyl Methoxide

The hydrogen-atom transfer from methoxy radical to nitric oxide, leading to the formation of formaldehyde and nitroxyl, represents a secondary reaction of photodissociation of methyl nitrite, which is used as rocket fuel. In this study, we explored the potential energy profile of the hydrogen-atom transfer using the electronic structure calculations at the DLPNO-CCSD(T)/aug-cc-pVTZ level of theory for two isomeric forms (cis and trans) of the pre-reaction complex. The cis-oriented pre-reaction complex has a weak elongated O─O bond, which gets further elongated in the hydrogen transfer transition state. This O─O bond stabilizes the pre-reaction complex by 32.9 kJ/mol. The O─O-induced stabilization is even greater for the transition state (48.2 kJ/mol), which was unexpected because of the larger O─O distance in the transition state structure. To address this paradox, we performed the electronic structure analysis of the reaction participants using the valence bond (VB) theory, natural resonance theory, topological analysis of the electron density and its derivatives, and analysis of the electron localization function distribution. This combined analysis led to the conclusion that the cis-transition state for hydrogen transfer, instead of being directly stabilized by the O─O interaction, gained substantial stabilization from the in-plane five-center six-electron aromaticity.

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

We explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

Photodissociation Dynamics of CO-Forming Channels on the Ground-State Surface of Methyl Formate at 248 nm: Direct Dynamics Study and Assessment of Generalized Multicenter Impulsive Models

The photodissociation dynamics of methyl formate in the electronic ground state S₀, initiated by a 248 nm-wavelength laser, is studied by direct dynamics simulations. We analyze five channels, where four of them have as products CH₃OH + CO, one leading to the formation of three fragments, H₂CO + H₂ + CO, and a channel characterized by a roaming transition state. The analysis of energy distribution among the degrees of freedom of the product and the comparison with experimental results previously published by other groups provide the ingredients to distinguish the examined dissociation pathways. The interpretation of the results proves that the characterization of dissociation mechanisms must rely on a dynamics approach involving multiple electronic states, including considerations on the features of the S₁/S₀ conical intersection. Here, we also assess the generalized multicenter impulsive model, GMCIM, that has been designed for dissociation processes with exit barriers, and the energy distribution in the products is predicted on the basis of information from the saddle points and the intrinsic reaction coordinates. Main features, advantages, limits, and future perspectives of the method are reported and discussed.

The impact of NOx addition on the ignition behaviour of n-pentane

Modern engine concepts present several opportunities for nitrogen combustion chemistry, particularly the interaction of NOx (NO + NO2) with fuel fragments and products of partial combustion.

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.

Theoretical Investigations on the Possibility of Prebiotic HCN Formation via O-Addition Reactions

Until now, reactions between methane photolysis products (CH₃^•, CH₂) and active N atom or reactive NO radical are proposed as routes of HCN formation in the prebiotic Earth. Scientists think that the reducing atmosphere of primitive Earth was made of H₂, He, N₂, NO, CH₄, H₂O, CO₂, etc., and there was no molecular oxygen. However, it has been evident from experiments that the vacuum ultraviolet (VUV) photolysis of CO₂ can produce atomic oxygen. Therefore, it can be presumed that atomic oxygen was likely present in early Earth's atmosphere. Was there any impact of atomic oxygen in production of early atmospheric HCN for the emergence of life? To hunt for the answer, we have employed computational methods to study the mechanism and kinetics of CH₃NO + O(¹D) and CH₂NO^• + O(³P) addition reactions. Current study suggests that the addition of O(¹D) into nitrosomethane (CH₃NO) and the addition of O(³P) into nitrosomethylene radical (CH₂NO^•) can efficiently produce HCN through an effectively barrierless pathway. At STP, Bartis-Widom phenomenological loss rate coefficients of O(¹D) and O(³P) are obtained as 2.47 × 10^-12 and 4.67 × 10^-11 cm³ molecule^-1 s^-1, respectively. We propose that addition reactions of atomic oxygen with CH₃NO and CH₂NO^• might act as a potential source for early atmospheric HCN.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

Auto-ignition behaviors of nitromethane in diluted oxygen in a rapid compression machine: Critical conditions for ignition, ignition delay times measurements, and kinetic modeling interpretation

In this work, the weakest thermodynamic conditions for the auto-ignition of mixtures containing nitromethane were experimentally determined by using the rapid compression machine facility. Results show there is a narrow weak ignition region between ignition and non-ignition. The weak ignition region would disappear with the increase of the EOC (end of compression) pressure and nitromethane concentration. In addition, the ignition delay times for successful auto-ignition for different nitromethane concentrations and equivalence ratio mixtures were measured and compared. Results show that the dependence of nitromethane ignition on the equivalence ratio is weak. Subsequently, the measured ignition delay time data were employed to validate several kinetic models in literature and our previous model shows better agreement with experimental results, as well as other available literature data. Sensitivity analysis for the model reveals the importance of unimolecular decomposition and H-abstraction reactions for the ignition delay times in the temperature range studied herein. Finally, critical conditions for nitromethane ignition under extended conditions that are beyond the ability of the experimental facility were predicted.

Shock Tube Laser Schlieren Study of the Pyrolysis of Isopropyl Nitrate

The decomposition of isopropyl nitrate was measured behind incident shock waves using laser schlieren densitometry in a diaphragmless shock tube. Experiments were conducted over the temperature range of 700-1000 K and at pressures of 71, 126, and 240 Torr. Electronic structure theory and RRKM Master Equation methods were used to predict the decomposition kinetics. RRKM/ME parameters were optimized against the experimental data to provide an accurate prediction over a broader range of conditions. The initial decomposition i-C₃H₇ONO₂ ⇌ i-C₃H₇O + NO₂ has a high-pressure limit rate coefficient of 5.70 × 10²²T^-1.80 exp[-21287.5/T] s^-1. A new chemical kinetic mechanism was developed to model the chemistry after the initial dissociation. A new shock tube module was developed for Cantera, which allows for arbitrarily large mechanisms in the simulation of laser schlieren experiments. The present work is in good agreement with previous experimental studies.

A modular, multi-diagnostic, automated shock tube for gas-phase chemistry

A new shock tube has been constructed for investigations of high-temperature chemical kinetics with an emphasis on combustion chemistry. This instrument includes a diaphragmless driver and electrical control of valving. A diaphragmless design significantly improves repeatability of experimental conditions vs the use of diaphragms and leads to an approximate order of magnitude reduction in turnaround time between experiments. Electrical control of valves, combined with diaphragmless operation, also enables remote and automated operation of the shock tube. The design allows for both incident and reflected shock experiments with multiple diagnostics. The performance of the shock tube is demonstrated by reproducing previous literature measurements on the unimolecular decomposition of isobutyl nitrite and cyclohexene.

Toward reliable characterization of energetic materials: interplay of theory and thermal analysis in the study of the thermal stability of tetranitroacetimidic acid (TNAA)

The thermal stability of energetic materials, being of the utmost importance for safety issues, is often considered in terms of kinetics, e.g., the Arrhenius parameters of the decomposition rate constant. The latter, in turn, are commonly determined using conventional thermoanalytical procedures with the use of simple Kissinger or Ozawa methods for kinetic data processing. However, thermal decomposition of energetic materials typically occurs via numerous exo- and endothermal processes including fast parallel reactions, phase transitions, autocatalysis, etc. This leads to numerous drawbacks of simple approaches. In this paper, we proposed a new methodology for characterization of the thermochemistry and thermal stability of melt-cast energetic materials, which is comprised of a complementary set of experimental and theoretical techniques in conjunction with a suitable kinetic model. With the aid of the proposed methodology, we studied in detail a novel green oxidizer, tetranitroacetimidic acid (TNAA). The experimental mass loss kinetics in the melt was perfectly fitted with a model comprised of zero-order reaction (sublimation or evaporation) and first-order thermal decomposition of TNAA with the effective Arrhenius parameters Ea = 41.0 ± 0.2 kcal mol-1 and log(A/s-1) = 20.2 ± 0.1. We rationalized the experimental findings on the basis of highly accurate CCSD(T)-F12 quantum chemical calculations. Computations predict that thermolysis of TNAA involves an intricate interplay of multiple decomposition channels of the three tautomers, which are equilibrated via either monomolecular reactions or concerted double hydrogen atom transfer in the H-bonded dimers; the calculated Arrhenius parameters of the effective rate constant coincide well with experiment. Most importantly, calculations provide detailed mechanistic evidence missing in the thermoanalytical experiment and explain formation of the experimentally observed primary products N2O and NO2. Along with the kinetics and mechanism of decomposition, the proposed approach yields accurate thermochemistry and phase change data of TNAA.

Theories and simulations of roaming

The phenomenon of roaming in chemical reactions has now become both commonly observed in experiment and extensively supported by theory and simulations. Roaming occurs in highly-excited molecules when the trajectories of atomic motion often bypass the minimum energy pathway and produce reaction in unexpected ways from unlikely geometries. The prototypical example is the unimolecular dissociation of formaldehyde (H₂CO), in which the "normal" reaction proceeds through a tight transition state to yield H₂ + CO but for which a high fraction of dissociations take place via a "roaming" mechanism in which one H atom moves far from the HCO, almost to dissociation, and then returns to abstract the second H atom. We review below the theories and simulations that have recently been developed to address and understand this new reaction phenomenon.

First kinetic study of the atmospherically important reactions BrHg˙ + NO2 and BrHg˙ + HOO

We use computational chemistry to determine the rate constants and product yields for the reactions of BrHg˙ with the atmospherically abundant radicals NO₂ and HOO. The reactants, products, and well-defined transition states are characterized using CCSD(T) with large basis sets. The potential energy profiles for the barrierless addition of HOO and NO₂ to BrHg˙ are characterized using CASPT2 and RHF-CCSDT, and the rate constants are computed as a function of temperature and pressure using variational transition state theory and master equation simulations. The calculated rate constant for the addition of NO₂ to BrHg˙ is larger than that for the addition of HOO by a factor of up to two under atmospheric conditions. For the reaction of HOO with BrHg˙ the addition reaction entirely dominates competing HOO + BrHg˙ reaction channels. The addition of NO₂ to BrHg˙ initially produces both BrHgNO₂ and BrHgONO, but after a few seconds under atmospheric conditions the sole product is syn-BrHgONO. A previously unsuspected reaction channel for BrHg˙ + NO₂ competes with the addition to yield Hg + BrNO₂. This reaction reduces the mercury oxidation state in BrHg˙ from Hg(i) to Hg(0) and slows the atmospheric oxidation of Hg(0). While the rate constant for this reduction channel is not well-constrained by the present calculations, it may be as much as 18% as large as the oxidation channel under some atmospheric conditions. As no experimental kinetic or product yield data are available for the reactions studied here, this work will provide guidance for atmospheric modelers and experimental kineticists.

Flash Pyrolysis of t-Butyl Hydroperoxide and Di-t-butyl Peroxide: Evidence of Roaming in the Decomposition of Organic Hydroperoxides

Thermal decomposition of t-butyl hydroperoxide and di-t-butyl peroxide was investigated using flash pyrolysis (in a short reaction time of <100 μs) and vacuum-ultraviolet (λ = 118.2 nm) single-photon ionization time-of-flight mass spectrometry (VUV-SPI-TOFMS) at temperatures up to 1120 K and quantum computational methods. Acetone and methyl radical were detected as the predominant products in the initial decomposition of di-t-butyl peroxide via O-O bond fission. In the initial dissociation of t-butyl hydroperoxide, acetone, methyl radical, isobutylene, and isobutylene oxide products were identified. The novel detection of the unimolecular formation of isobutylene oxide, as supported by the computational study, was found to proceed via a roaming hydroxyl radical facilitated by a hydrogen-bonded intermediate. This new pathway could provide a new class of reactions to consider in the modeling of the low temperature oxidation of alkanes.

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.