HO2 + NO2: Kinetics, Thermochemistry, and Evidence for a Bimolecular Product Channel

A master equation (ME) analysis of available experimental data has been carried out on the reaction HO₂ + NO₂ + M ⇋ HO₂NO₂ + M (1a)/(-1a). The analysis, based on the ME code MESMER, uses both the association and dissociation kinetic data from the literature, and provides improved thermochemistry on reaction 1a. Our preferred model assigns two low-frequency vibrations of HO₂NO₂ as hindered rotors and couples these to the external rotations. This model gives Δ_rH°₀(1a) = -93.9 ± 1.0 kJ mol^-1, which implies that Δ_fH°₀ HO₂NO₂ = -42.0 ± 1.0 kJ mol^-1 (uncertainties are 2σ). A significant contributor to the uncertainty derives from modeling the interaction between the internal and external rotors. Using this improved kinetics for reaction 1a/-1a, data at elevated temperatures, 353-423 K, which show no evidence of the expected equilibration, have been reanalyzed, indicating that an additional reaction is occurring that masks the equilibration. Based on a published ab initio study, this additional channel is assigned to the bimolecular reaction HO₂ + NO₂ → H-NO₂ + O₂ (1b); H-NO₂ is nitryl hydride and has not previously been directly observed in experiments. The output of the master equation analysis has been parametrized and Troe expressions are provided for an improved description of k_1a(p,T) and k_-1a(p,T).

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.

Adaptive Fitting of Potential Energy Surfaces of Small to Medium-Sized Molecules in Sum-of-Product Form: Application to Vibrational Spectroscopy

A semi-automatic sampling and fitting procedure for generating sum-of-product (Born-Oppenheimer) potential energy surfaces based on a high-dimensional model representation is presented. The adaptive sampling procedure and subsequent fitting relies on energies only and can be used for re-fitting existing analytic potential energy surfaces in sum-of-product form or for direct fits from ab initio computa- tions. The method is tested by fitting ground electronic state potential energy surfaces for small to medium sized semi-rigid molecules, i.e., HFCO, HONO, and HCOOH, based upon ab initio computations at the CCSD(T)-F12/cc-pVTZ-F12 or MP2/aug-cc-pVTZ levels of theory. Vibrational eigenstates are computed using block improved relaxation in the Heidelberg MCTDH package and compared to available experimental and theoretical data. The new potential energy surfaces are compared to the best ones currently available for these molecules, in terms of accuracy, including of resulting vibrational states, required numbers of sampling points, and number of fitting parameters. The present procedure leads to compact expansions and scales well with the number of dimensions for simple potentials such as single or double wells.

Prediction of Reactive Nitrous Acid Formation in Rare-Earth MOFs via ab initio Molecular Dynamics

Reactive gas formation in pores of metal-organic frameworks (MOFs) is a known mechanism of framework destruction; understanding those mechanisms for future durability design is key to next generation adsorbents. Herein, an extensive set of ab initio molecular dynamics (AIMD) simulations are used for the first time to predict competitive adsorption of mixed acid gases (NO2 and H2 O) and the in-pore reaction mechanisms for a series of rare earth (RE)-DOBDC MOFs. Spontaneous formation of nitrous acid (HONO) is identified as a result of deprotonation of the MOF organic linker, DOBDC. The unique DOBDC coordination to the metal clusters allows for proton transfer from the linker to the NO2 without the presence of H2 O and may be a factor in DOBDC MOF durability. This is a previously unreported mechanisms of HONO formation in MOFs. With the presented methodology, prediction of future gas interactions in new nanoporous materials can be achieved.

Threshold photoionization shows no sign of nitryl hydride in methane oxidation with nitric oxide

No nitryl hydride was detected in partial oxidation of nitric oxide doped methane, despite recent theoretical reaction rates suggesting otherwise.

Hidden complexities in the reaction of H2O2 and HNO revealed by ab initio quantum chemical investigations

Nitroxyl (HNO) and hydrogen peroxide have both been implicated in a variety of reactions relevant to environmental and physiological processes and may contribute to a unique, unexplored, pathway for the production of nitrous acid (HONO) in soil. To investigate the potential for this reaction, we report an in-depth investigation of the reaction pathway of H₂O₂ and HNO forming HONO and water. We find the breaking of the peroxide bond and a coupled proton transfer in the first step leads to hydrogen nitryl (HNO₂) and an endogenous water, with an extrapolated NEVPT2 (multireference perturbation theory) barrier of 29.3 kcal mol^-1. The first transition state is shown to possess diradical character linking the far peroxide oxygen to the bridging, reacting, peroxide oxygen. The energy of this first step, when calculated using hybrid density functional theory, is shown to depend heavily on the amount of Hartree-Fock exchange in the functional, with higher amounts leading to a higher barrier and more diradical character. Additionally, high amounts of spin contamination cause CCSD(T) to significantly overestimate the TS1 barrier with a value of 36.2 kcal mol^-1 when using the stable UHF wavefunction as the reference wavefunction. However, when using the restricted Hartree-Fock reference wavefunction, the TS1 CCSD(T) energy is lowered to yield a barrier of 31.2 kcal mol^-1, in much better agreement with the NEVPT2 result. The second step in the reaction is the isomerization of HNO₂ to trans-HONO through a Grotthuss-like mechanism accepting a proton from and donating a proton to the endogenous water. This new mechanism for the isomerization of HNO₂ is shown to have an NEVPT2 barrier of 23.3 kcal mol^-1, much lower than previous unimolecular estimates not including an explicit water. Finally, inclusion of an additional explicit water is shown to lower the HNO₂ isomerization barrier even further.

Bimolecular Reactions between Dimethylnitramine and Its Radical Decomposition Products

Bimolecular reactions between intact nitramines and their radical decomposition products can accelerate thermal decomposition, yet the detailed mechanisms of such reactions are not well understood. We have used density functional theory at the M06/6-311++G(3df,3pd) level to locate transition structures and compute 0 K activation barriers for various gas-phase reactions that may contribute to radical-assisted decomposition of dimethylnitramine (DMNA, (CH₃)₂NNO₂). Our calculations indicate that H abstraction from DMNA is the lowest-barrier mechanism for most radicals and a subsequent N-N β-scission in the alkyl radical 3 leads to an imine intermediate and NO₂. H abstraction is thus responsible for conversion of most radicals to NO₂. Also, among the nine radicals considered, NO is found to be least reactive and its reactions with DMNA yield dimethylnitrosoamine (DMNSA, (CH₃)₂NNO), a known product of DMNA decomposition.

Modeling Nitrogen Species as Pollutants: Thermochemical Influences

To simulate emissions of nitrogen-containing compounds in practical combustion environments, it is necessary to have accurate values for their thermochemical parameters, as well as accurate kinetic values to describe the rates of their formation and decomposition. Significant disparity is observed in the literature for the former, and we therefore present herein high-accuracy ab initio gas-phase thermochemistry for 60 nitrogenous compounds, many of which are important in the formation and consumption chemistry of NOx species. Several quantum-chemical composite methods (CBS-APNO, G3, and G4) were utilized to derive enthalpies of formation via the atomization method. Entropies and heat capacities were calculated from traditional statistical thermodynamics, with oscillators treated as anharmonic based on ro-vibrational property analyses carried out at the B3LYP/cc-pVTZ level of theory. The use of quantum chemical methods, along with the treatments of anharmonicities and hindered rotors, ensures accurate enthalpy of formation, entropy, and heat capacity values across the temperature range 298.15-3000 K. The implications of these results for atmospheric and combustion modeling are discussed.

Photoelectron Spectroscopy of cis-Nitrous Acid Anion (cis-HONO(-))

We report photoelectron spectra of cis-HONO(-) formed from an association reaction of OH(-) and NO in a pulsed, plasma-entrainment ion source. The experimental data are assigned to the cis-HONO(-) isomer, which is predicted to be the global minimum on the anion potential energy surface. We do not find evidence for a significant contribution from trans-HONO(-). Electron photodetachment of cis-HONO(-) with 1613, 1064, 532, 355, and 301 nm photons accesses the ground X̃ (1)A' (S0) and excited ã (3)A″ (T1) states of neutral HONO. The photoelectron spectrum resulting from detachment forming cis-HONO (S0) exhibits a long vibrational progression, dominated by overtones and combination bands involving the central O-N stretching and ONO bending vibrations. This indicates that there is a significant change in the central O-N bond length between cis-HONO(-) and cis-HONO (S0). The electron affinity (EA) of cis-HONO is determined to be 0.356(8) eV. We also report the dissociation energy (D0) of cis-HONO(-), forming OH(-) + NO, as 0.594(9) eV, which is a factor of 4 decrease in the central O-N bond strength compared to neutral cis-HONO. The T1 state of cis-HONO is shown to be ∼2.3 eV higher in energy than cis-HONO (S0). Electron photodetachment to form cis-HONO (T1) accesses a transition state along the HO-NO bond dissociation coordinate. The resulting photoelectron spectrum exhibits broad peaks spaced by the terminal N═O stretching frequency. Electronic structure calculations and photoelectron spectrum simulations reported here show very good agreement with the experimental data.

Theoretical Study on the Dynamics of the Reaction of HNO((1)A') with HO2((2)A″)

We used stochastic one-dimensional chemical master equation (CME) simulation to gain insight into the dynamics of the reaction of HNO((1)A') with HO2((2)A″). The reaction takes place over a multiwell, multichannel potential energy surface that is based on the computations at the CBS-QB3 level of theory. The calculated multipath potential energy surface consists of three potential wells and three van der Waals complexes. In solving the master equation, the Lennard-Jones potential is used to model the collision between the collider gases. The fractional population of different intermediates and products in the early stages of the reaction is examined to determine the role of the energized intermediates and van der Waals complexes on the kinetics of the title reaction. The major products of the title reaction at lower temperatures are OH, HNO2, HNOH, and O2(X(3)Σg(-)). The temperature- and pressure-dependence of the reaction over a wide range of temperature (300-3000 K) and pressure (0.1-2000 Torr) are studied. No sign of pressure dependence was being observed for the title reaction over the stated range of pressure. The calculated rate constants from the CME simulation are compared with those obtained from the RRKM-SSA method that is based on strong collision assumption. Our results indicate that the strong collision assumption increases the calculated rate constant for the formation of the main products (HNO2 + OH) by a factor of 2 at 300 K and 1 atm pressure, compared to the results of CME simulation, although the results are in good agreement at higher temperatures.

Nitrate formation from atmospheric nitrogen and oxygen photocatalysed by nano-sized titanium dioxide

The concentration of nitrate in aquatic systems is rising with the development of modern industry and agriculture, causing a cascade of environmental problems. Here we describe a previously unreported nitrate formation process. Both indoor and outdoor experiments are conducted to demonstrate that nitrate may be formed from abundant atmospheric nitrogen and oxygen on nano-sized titanium dioxide surfaces under UV or sunlight irradiation. We suggest that nitric oxide is an intermediate product in this process, and elucidate its formation mechanisms using first-principles density functional theory calculations. Given the expanding use of titanium dioxide worldwide, such a titanium dioxide-mediated photocatalysis process may reveal a potentially underestimated source of nitrate in the environment, which on one hand may lead to an increasing environmental pollution concern, and on the other hand may provide an alternative, gentle and cost-effective method for nitrate production.

Detection and structure of HOON: microwave spectroscopy reveals an O-O bond exceeding 1.9 Å

Nitric oxide (NO) reacts with hydroxyl radicals (OH) in the gas phase to produce nitrous acid, HONO, but essentially nothing is known about the isomeric nitrosyl-O-hydroxide (HOON), owing to its perceived instability. We report the detection of gas-phase HOON in a supersonic molecular beam by Fourier transform microwave spectroscopy and a precise determination of its molecular structure by further spectroscopic analysis of its (2)H, (15)N, and (18)O isotopologs. HOON contains the longest O-O bond in any known molecule (1.9149 ± 0.0005 Å) and appears surprisingly stable, with an abundance roughly 3% that of HONO in our experiments.

FONO: a difficult case for theory. The ELF and ELI-D topological studies on the chemical bonding using correlated wavefunctions

The complicated nature of the chemical bonding in cis and trans isomers of F-O-N=O is discussed based on the results obtained from the topological analysis of electron localization function (η) (ELF), electron localizability index (Y(D)(σ)), and electron density (ρ). The calculations have been performed for correlated wavefunctions using the CCSD and CASSCF methods. The F-O1 bond with non-bonding basins, V(F) and V(')(O1), belongs to the protocovalent type (η,Y(D)(σ)) and its total population ranges between 0.2 and 0.4e. The central N-O1 bond in the cis form is protocovalent (η, Y(D)(σ)) with two basins, V(N) and V(O1). The total population oscillates between 0.7 and 0.9e. In the trans isomer, topology of ELF depends on used method. At the CCSD level only one non-bonding basin, V(N), is observed (η). Its population is about 0.5e. According to the definition of a heteronuclear charge-shift (CS) bond, only N-O1 bond in trans-FONO belongs to the CS class. A relation between η- and ρ-topology and N-O1 bond length is discussed.