The Low-Lying Electronic States of NO2: Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum

UV-spectrum and photodecomposition of peroxynitrous acid in the troposphere

The UV spectrum of peroxynitrous acid, HOONO, was computed at the B3LYP/AVTZ and MCSCF/AVTZ levels using the fewest switches surface hopping algorithm. Due to large-amplitude vibrational motions of this molecule, the maxima in the simulated spectra are displaced from the positions of vertical excitations. The three lowest excited electronic singlet states, which are all repulsive, can be reached by UV absorption. The photolysis products are determined, and the photolysis rate constant is provided for the first time. We found that near the tropopause the photolysis rate constant J ≈ 6 × 10-4 s-1, exceeds that for thermal decomposition by two orders of magnitude. The photolysis lifetime is about 30 minutes. Thus, photolysis is an important process and should be included in atmospheric models.

A white cell based broadband transient UV-vis absorption spectroscopy with pulsed laser photolysis reactors for chemical kinetics under variable temperatures and pressures

UV-vis spectroscopy is widely used for kinetic studies in physical chemistry, as species' absolute cross-sections are usually less sensitive to experimental conditions (i.e., temperature and pressure). Here, we present the design and characterization of a multipass UV-vis absorption spectroscopy white cell coupled to a pulsed-laser photolysis flow reactor. The glass reactor was designed to facilitate studies of gas phase chemical reactions over a range of conditions (239-293 K and 10-550 Torr). Purged windows mitigate contamination from chemical precursors and photolysis products. We report the measured impact of this purging on temperature uniformity and the absorption length and present some supporting flow calculations. The combined optical setup is unique and enables the photolysis laser to be coaligned with a well-defined absorption pathlength probe beam. This alignment leverages the use of one long-pass filter to increase the spectrum flatness and increase the light intensity vs other systems that use two dichroic mirrors. The probe beam is analyzed with a dual exit spectrograph, customized to split the light between an intensified CCD and photomultiplier tube, enabling simultaneous spectrum and single wavelength detection. This multipass system yields a pathlength of ∼450 cm and minimum observable concentrations of ∼3.7 × 1011 molecule cm-3 (assuming cross-sections ∼1.2 × 10-17 cm2). The temperature profile across the reaction region is ±2 K, defined by the worst-case temperature of 239 K, validated by measurements of the N2O4 equilibrium constant. Finally, the system is implemented to study the simplest Criegee intermediate, demonstrating the instrument performance and advantages of simultaneous spectrum and temporal profile measurements.

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.

Formation of N(2D) from Hyperthermal Collisions between O(3P) and NO(X2Π)

Hyperthermal collisions between O(3P) and NO(X2Π) could lead to the formation of the first electronically excited atomic nitrogen (N(2D)), which plays a key role in plasma formation in shock-heated air. This process is facilitated mainly by four doublet states, and to a much lesser extent by two quartet states. In this work, we report quasi-classical trajectory studies of this reactive process using the four analytical adiabatic potential energy surfaces for the doublet states developed previously from fitting high-level ab initio data. The reactions were found to be strongly enhanced by vibrational excitation of the NO reactant, consistent with the existence of potential energy barriers in the exit channel. Despite the large endothermicity of the reaction, the rate coefficient is significant at high temperatures, suggesting a possible role of this reaction in the hyperthermal kinetics in the shock layer of a hypersonic (re)entry vehicle.

Temperature Dependence of the Electronic Absorption Spectrum of NO2

The nitrogen dioxide (NO₂) radical is composed of the two most abundant elements in the atmosphere, where it can be formed in a variety of ways including combustion, detonation of energetic materials, and lightning. Relevant also to smog and ozone cycles, together these processes span a wide range of temperatures. Remarkably, high-resolution NO₂ electronic absorption spectra have only been reported in a narrow range below about 300 K. Previously, we reported [ J. Phys. Chem. A 2021, 125, 5519-5533] the construction of quasi-diabatic potential energy surfaces (PESs) for the lowest four electronic states (X̃, Ã, B̃, and C̃) of NO₂. In addition to three-dimensional PESs based on explicitly correlated MRCI(Q)-F12/VTZ-F12 ab initio data, the geometry dependence of each component of the dipoles and transition dipoles was also mapped into fitted surfaces. The multiconfigurational time-dependent Hartree (MCTDH) method was then used to compute the 0 K electronic absorption spectrum (from the ground rovibrational initial state) employing those energy and transition dipole surfaces. Here, in an extension of that work, we report an investigation into the effects of elevated temperature on the spectrum, considering the effects of the population of rotationally and vibrationally excited initial states. The calculations are complemented by new experimental measurements. Spectral contributions from hundreds of rotational states up to N = 20 and from 200 individually-characterized vibrational states were computed. A spectral simulation tool was developed that enables modeling the spectrum at various temperatures─by weighting individual spectral contributions via the partition function, or for pure excited initial states, which can be probed via transient absorption spectroscopy. We validate these results against experimental absorption spectroscopy data at high temperatures, as well as via a new measurement from the (1,0,1) initial vibrational state.

Franck-Condon spectra of unbound and imaginary-frequency vibrations via correlation functions: a branch-cut free, numerically stable derivation

Molecular electronic spectra can be represented in the time domain as auto-correlation functions of the initial vibrational wavepacket. We present a derivation of the harmonic vibrational auto-correlation function that is valid for both real and imaginary harmonic frequencies. The derivation rests on Lie algebra techniques that map otherwise complicated exponential operator arithmetic to simpler matrix formulae. The expressions for the zero- and finite-temperature harmonic auto-correlation functions have been carefully structured both to be free of branch-cut discontinuities and to remain numerically stable with finite-precision arithmetic. Simple extensions correct the harmonic Franck-Condon approximation for the lowest-order anharmonic and Herzberg-Teller effects. Quantitative simulations are shown for several examples, including the electronic absorption spectra of F2, HOCl, CH2NH, and NO2.

Reactive and Nonreactive Collisions between NO(X2Π) and O(3P) under Hyperthermal Conditions

Quasiclassical trajectory calculations are performed for hyperthermal collisions between NO(X²Π) and O(³P) on recently developed potential energy surfaces for the lowest doublet and quartet states of the NO₂ system. Three product channels are investigated, and their branching fractions are in reasonably good agreement with the recent crossed molecular beam study at 84 kcal/mol of collision energy. The dominant inelastic channel has a strong forward scattering bias and a high translational energy distribution with limited internal excitation in the scattered NO. The exchange channel has significantly higher NO internal excitation and is also forward biased. The abstraction channel producing internally excited O₂ has the smallest branching fraction and a broader angular distribution also with a forward peak. The angular and translational energy distributions in the three channels are consistent with experiment, but the agreement is not always quantitative. The sources of the differences are discussed. Finally, the impact of NO vibrational excitation on the reactive channels and the corresponding rate coefficients are reported.

Electronic Structure of Nitrobenzene: A Benchmark Example of the Accuracy of the Multi-State CASPT2 Theory

The electronic structure of nitrobenzene (C₆H₅NO₂) has been studied by means of the complete active space self-consistent field (CASSCF) and multi-state second-order perturbation (MS-CASPT2) methods. To this end, an active space of 20 electrons distributed in 17 orbitals has been selected to construct the reference wave function. In this work, we have calculated the vertical excitation energies and the energy barrier for the dissociation of the molecule on the ground state into phenyl and nitrogen dioxide. After applying the corresponding vibrational corrections to the electronic energies, it is demonstrated that the MS-CASPT2//CASSCF values obtained in this work yield an excellent agreement between calculated and experimental data. In addition, other active spaces of lower size have been applied to the system in order to check the active space dependence in the results.