Secondary Organic Aerosol Mass Yields from NO3 Oxidation of α-Pinene and Δ-Carene: Effect of RO2 Radical Fate.

Dark chamber experiments were conducted to study the SOA formed from the oxidation of α-pinene and Δ-carene under different peroxy radical (RO2) fate regimes: RO2 + NO3, RO2 + RO2, and RO2 + HO2. SOA mass yields from α-pinene oxidation were <1 to ∼25% and strongly dependent on available OA mass up to ∼100 µg m-3. The strong yield dependence of α-pinene oxidation is driven by absorptive partitioning to OA and not by available surface area for condensation. Yields from Δ-carene + NO3 were consistently higher, ranging from ∼10-50% with some dependence on OA for <25 µg m-3. Explicit kinetic modeling including vapor wall losses was conducted to enable comparisons across VOC precursors and RO2 fate regimes and to determine atmospherically relevant yields. Furthermore, SOA yields were similar for each monoterpene across the nominal RO2 + NO3, RO2 + RO2, or RO2 + HO2 regimes; thus, the volatility basis sets (VBS) constructed were independent of the chemical regime. Elemental O/C ratios of ∼0.4-0.6 and nitrate/organic mass ratios of ∼0.15 were observed in the particle phase for both monoterpenes in all regimes, using aerosol mass spectrometer (AMS) measurements. An empirical relationship for estimating particle density using AMS-derived elemental ratios, previously reported in the literature for non-nitrate containing OA, was successfully adapted to organic nitrate-rich SOA. Observations from an NO3- chemical ionization mass spectrometer (NO3-CIMS) suggest that Δ-carene more readily forms low-volatility gas-phase highly oxygenated molecules (HOMs) than α-pinene, which primarily forms volatile and semivolatile species, when reacted with NO3, regardless of RO2 regime. The similar Δ-carene SOA yields across regimes, high O/C ratios, and presence of HOMs, suggest that unimolecular and multistep processes such as alkoxy radical isomerization and decomposition may play a role in the formation of SOA from Δ-carene + NO3. The scarcity of peroxide functional groups (on average, 14% of C10 groups carried a peroxide functional group in one test experiment in the RO2 + RO2 regime) appears to rule out a major role for autoxidation and organic peroxide (ROOH, ROOR) formation. The consistently substantially lower SOA yields observed for α-pinene + NO3 suggest such pathways are less available for this precursor. The marked and robust regime-independent difference in SOA yield from two different precursor monoterpenes suggests that in order to accurately model SOA production in forested regions the chemical mechanism must feature some distinction among different monoterpenes.

A Four Carbon Organonitrate as a Significant Product of Secondary Isoprene Chemistry.

Oxidation of isoprene by nitrate radicals (NO3) or by hydroxyl radicals (OH) under high NOx conditions forms a substantial amount of organonitrates (ONs). ONs impact NOx concentrations and consequently ozone formation while also contributing to secondary organic aerosol. Here we show that the ONs with the chemical formula C4H7NO5 are a significant fraction of isoprene-derived ONs, based on chamber experiments and ambient measurements from different sites around the globe. From chamber experiments we found that C4H7NO5 isomers contribute 5%-17% of all measured ONs formed during nighttime and constitute more than 40% of the measured ONs after further daytime oxidation. In ambient measurements C4H7NO5 isomers usually dominate both nighttime and daytime, implying a long residence time compared to C5 ONs which are removed more rapidly. We propose potential nighttime sources and secondary formation pathways, and test them using a box model with an updated isoprene oxidation scheme.

Heterogeneous Nucleation Drives Particle Size Segregation in Sequential Ozone and Nitrate Radical Oxidation of Catechol.

Secondary organic aerosol formation via condensation of organic vapors onto existing aerosol transforms the chemical composition and size distribution of ambient aerosol, with implications for air quality and Earth's radiative balance. Gas-to-particle conversion is generally thought to occur on a continuum between equilibrium-driven partitioning of semivolatile molecules to the pre-existing mass size distribution and kinetic-driven condensation of low volatility molecules to the pre-existing surface area size distribution. However, we offer experimental evidence in contrast to this framework. When catechol is sequentially oxidized by O3 and NO3 in the presence of (NH4)2SO4 seed particles with a single size mode, we observe a bimodal organic aerosol mass size distribution with two size modes of distinct chemical composition with nitrocatechol from NO3 oxidation preferentially condensing onto the large end of the pre-existing size distribution (∼750 nm). A size-resolved chemistry and microphysics model reproduces the evolution of the two distinct organic aerosol size modes─heterogeneous nucleation to an independent, nitrocatechol-rich aerosol phase.

Gas-Particle Partitioning and SOA Yields of Organonitrate Products from NO3-Initiated Oxidation of Isoprene under Varied Chemical Regimes.

Alkyl nitrate (AN) and secondary organic aerosol (SOA) from the reaction of nitrate radicals (NO3) with isoprene were observed in the Simulation of Atmospheric PHotochemistry In a large Reaction (SAPHIR) chamber during the NO3Isop campaign in August 2018. Based on 15 day-long experiments under various reaction conditions, we conclude that the reaction has a nominally unity molar AN yield (observed range 90 ± 40%) and an SOA mass yield of OA + organic nitrate aerosol of 13-15% (with ∼50 µg m-3 inorganic seed aerosol and 2-5 µg m-3 total organic aerosol). Isoprene (5-25 ppb) and oxidant (typically ∼100 ppb O3 and 5-25 ppb NO2) concentrations and aerosol composition (inorganic and organic coating) were varied while remaining close to ambient conditions, producing similar AN and SOA yields under all regimes. We observe the formation of dinitrates upon oxidation of the second double bond only once the isoprene precursor is fully consumed. We determine the bulk partitioning coefficient for ANs (K p ∼ 10-3 m3 µg-1), indicating an average volatility corresponding to a C5 hydroxy hydroperoxy nitrate.

Coupling of organic and inorganic aerosol systems and the effect on gas-particle partitioning in the southeastern US.

Several models were used to describe the partitioning of ammonia, water, and organic compounds between the gas and particle phases for conditions in the southeastern US during summer 2013. Existing equilibrium models and frameworks were found to be sufficient, although additional improvements in terms of estimating pure-species vapor pressures are needed. Thermodynamic model predictions were consistent, to first order, with a molar ratio of ammonium to sulfate of approximately 1.6 to 1.8 (ratio of ammonium to 2× sulfate, RN/2S ≈ 0.8 to 0.9) with approximately 70% of total ammonia and ammonium (NH x ) in the particle. Southeastern Aerosol Research and Characterization Network (SEARCH) gas and aerosol and Southern Oxidant and Aerosol Study (SOAS) Monitor for AeRosols and Gases in Ambient air (MARGA) aerosol measurements were consistent with these conditions. CMAQv5.2 regional chemical transport model predictions did not reflect these conditions due to a factor of 3 overestimate of the nonvolatile cations. In addition, gas-phase ammonia was overestimated in the CMAQ model leading to an even lower fraction of total ammonia in the particle. Chemical Speciation Network (CSN) and aerosol mass spectrometer (AMS) measurements indicated less ammonium per sulfate than SEARCH and MARGA measurements and were inconsistent with thermodynamic model predictions. Organic compounds were predicted to be present to some extent in the same phase as inorganic constituents, modifying their activity and resulting in a decrease in [H+]air (H+ in µgm-3 air), increase in ammonia partitioning to the gas phase, and increase in pH compared to complete organic vs. inorganic liquid-liquid phase separation. In addition, accounting for nonideal mixing modified the pH such that a fully interactive inorganic-organic system had a pH roughly 0.7 units higher than predicted using traditional methods (pH = 1.5 vs. 0.7). Particle-phase interactions of organic and inorganic compounds were found to increase partitioning towards the particle phase (vs. gas phase) for highly oxygenated (O : C≥0.6) compounds including several isoprene-derived tracers as well as levoglu-cosan but decrease particle-phase partitioning for low O: C, monoterpene-derived species.

Alkoxy Radical Bond Scissions Explain the Anomalously Low Secondary Organic Aerosol and Organonitrate Yields From α-Pinene + NO3.

Oxidation of monoterpenes (C10H16) by nitrate radicals (NO3) constitutes an important source of atmospheric secondary organic aerosol (SOA) and organonitrates. However, knowledge of the mechanisms of their formation is incomplete and differences in yields between similar monoterpenes are poorly understood. In particular, yields of SOA and organonitrates from α-pinene + NO3 are low, while those from Δ3-carene + NO3 are high. Using computational methods, we suggest that bond scission of the nitrooxy alkoxy radicals from Δ3-carene lead to the formation of reactive keto-nitrooxy-alkyl radicals, which retain the nitrooxy moiety and can undergo further reactions to form SOA. By contrast, bond scissions of the nitrooxy alkoxy radicals from α-pinene lead almost exclusively to the formation of the relatively unreactive and volatile product pinonaldehyde (C10H16O2), thereby limiting organonitrate and SOA formation. This hypothesis is supported by laboratory experiments that quantify products of the reaction of α-pinene + NO3 under atmospherically relevant conditions.

Nitrate radicals and biogenic volatile organic compounds: oxidation, mechanisms, and organic aerosol.

Oxidation of biogenic volatile organic compounds (BVOC) by the nitrate radical (NO3) represents one of the important interactions between anthropogenic emissions related to combustion and natural emissions from the biosphere. This interaction has been recognized for more than 3 decades, during which time a large body of research has emerged from laboratory, field, and modeling studies. NO3-BVOC reactions influence air quality, climate and visibility through regional and global budgets for reactive nitrogen (particularly organic nitrates), ozone, and organic aerosol. Despite its long history of research and the significance of this topic in atmospheric chemistry, a number of important uncertainties remain. These include an incomplete understanding of the rates, mechanisms, and organic aerosol yields for NO3-BVOC reactions, lack of constraints on the role of heterogeneous oxidative processes associated with the NO3 radical, the difficulty of characterizing the spatial distributions of BVOC and NO3 within the poorly mixed nocturnal atmosphere, and the challenge of constructing appropriate boundary layer schemes and non-photochemical mechanisms for use in state-of-the-art chemical transport and chemistry-climate models. This review is the result of a workshop of the same title held at the Georgia Institute of Technology in June 2015. The first half of the review summarizes the current literature on NO3-BVOC chemistry, with a particular focus on recent advances in instrumentation and models, and in organic nitrate and secondary organic aerosol (SOA) formation chemistry. Building on this current understanding, the second half of the review outlines impacts of NO3-BVOC chemistry on air quality and climate, and suggests critical research needs to better constrain this interaction to improve the predictive capabilities of atmospheric models.

Modeling the Current and Future Roles of Particulate Organic Nitrates in the Southeastern United States.

Organic nitrates are an important aerosol constituent in locations where biogenic hydrocarbon emissions mix with anthropogenic NOx sources. While regional and global chemical transport models may include a representation of organic aerosol from monoterpene reactions with nitrate radicals (the primary source of particle-phase organic nitrates in the Southeast United States), secondary organic aerosol (SOA) models can underestimate yields. Furthermore, SOA parametrizations do not explicitly take into account organic nitrate compounds produced in the gas phase. In this work, we developed a coupled gas and aerosol system to describe the formation and subsequent aerosol-phase partitioning of organic nitrates from isoprene and monoterpenes with a focus on the Southeast United States. The concentrations of organic aerosol and gas-phase organic nitrates were improved when particulate organic nitrates were assumed to undergo rapid (τ = 3 h) pseudohydrolysis resulting in nitric acid and nonvolatile secondary organic aerosol. In addition, up to 60% of less oxidized-oxygenated organic aerosol (LO-OOA) could be accounted for via organic nitrate mediated chemistry during the Southern Oxidants and Aerosol Study (SOAS). A 25% reduction in nitrogen oxide (NO + NO2) emissions was predicted to cause a 9% reduction in organic aerosol for June 2013 SOAS conditions at Centreville, Alabama.

Secondary organic aerosol formation and organic nitrate yield from NO3 oxidation of biogenic hydrocarbons.

The secondary organic aerosol (SOA) mass yields from NO3 oxidation of a series of biogenic volatile organic compounds (BVOCs), consisting of five monoterpenes and one sesquiterpene (α-pinene, ß-pinene, Δ-3-carene, limonene, sabinene, and ß-caryophyllene), were investigated in a series of continuous flow experiments in a 10 m(3) indoor Teflon chamber. By making in situ measurements of the nitrate radical and employing a kinetics box model, we generate time-dependent yield curves as a function of reacted BVOC. SOA yields varied dramatically among the different BVOCs, from zero for α-pinene to 38-65% for Δ-3-carene and 86% for ß-caryophyllene at mass loading of 10 µg m(-3), suggesting that model mechanisms that treat all NO3 + monoterpene reactions equally will lead to errors in predicted SOA depending on each location's mix of BVOC emissions. In most cases, organonitrate is a dominant component of the aerosol produced, but in the case of α-pinene, little organonitrate and no aerosol is formed.

Vibrational overtone initiated unimolecular dissociation of HOCH2OOH and HOCD2OOH: evidence for mode selective behavior.

The vibrational overtone induced unimolecular dissociation of HMHP (HOCH(2)OOH) and HMHP-d(2) (HOCD(2)OOH) into OH and HOCH(2)O (HOCD(2)O) fragments is investigated in the region of the 4nu(OH) and 5nu(OH) bands. The unimolecular dissociation rates in the threshold region, corresponding to the 4nu(OH) band, exhibit measurable differences associated with excitation of the OH stretch of the alcohol versus the peroxide functional group, with the higher energy alcohol OH stretching state exhibiting a slower dissociation rate compared to the lower energy peroxide OH stretch in both HMHP and HMHP-d(2). Predictions using the Rice-Ramsperger-Kassel-Marcus theory give rates that are in reasonably good agreement with the measured dissociation rate for the alcohol OH stretch but considerably differ from the measured rates for the peroxide OH stretch in both isotopomers. The present results are interpreted as suggesting that the extent of intramolecular vibrational energy redistribution (IVR) is different for the two OH stretching states associated with the two functional groups in HMHP, with IVR being substantially less complete for the peroxide OH stretch. Analysis of the OH fragment product state distributions in conjunction with phase-space theory simulation gives a D(0) value of 38+/-0.7 kcal/mole for breaking the peroxide bond in HMHP.

OH-stretch vibrational spectroscopy of hydroxymethyl hydroperoxide.

We report measurement and analysis of the photodissociation spectrum of hydroxymethyl hydroperoxide (HOCH(2)OOH) and its partially deuterated analogue, HOCD(2)OOH, in the OH-stretching region. Spectra are obtained by Fourier transform infrared spectroscopy in the 1nu(OH) and 2nu(OH) regions, and by laser induced fluorescence detection of the OH fragment produced from dissociation of HOCH(2)OOH initiated by excitation of the 4nu(OH) and 5nu(OH) overtone regions (action spectroscopy). A one-dimensional local-mode model of each OH chromophore is used with ab initio calculated OH-stretching potential energy and dipole moment curves at the coupled-cluster level of theory. Major features in the observed absorption and photodissociation spectra are explained by our local-mode model. In the 4nu(OH) region, explanation of the photodissocation spectrum requires a nonuniform quantum yield, which is estimated by assuming statistical energy distribution in the excited state. Based on the estimated dissociation threshold, overtone photodissociation is not expected to significantly influence the atmospheric lifetime of hydroxymethyl hydroperoxide.

Rotational spectroscopy and dipole moment of cis-cis HOONO and DOONO.

The rotational spectrum of cis-cis HOONO has been studied over a broad range of frequencies, 13-840 GHz, using pulsed beam Fourier-transform microwave spectroscopy and room-temperature flow cell submillimeter spectroscopy. The rotational spectrum of the deuterated isotopomer, cis-cis DOONO, has been studied over a subset of this range, 84-640 GHz. Improved spectroscopic constants have been determined for HOONO, and the DOONO spectrum is analyzed for the first time. Weak-field Stark effect measurements in the region of 84-110 GHz have been employed to determine the molecular dipole moments of cis-cis HOONO [mu(a) = 0.542(8) D, mu(b) = 0.918(15) D, mu = 1.07(2) D] and DOONO [mu(a) = 0.517(9) D, mu(b) = 0.930(15) D, mu = 1.06(2) D]. The quadrupole coupling tensor in the principal inertial axis system for the 14N nucleus has been determined to be chi(aa) = 1.4907(25) MHz, chi(bb) = -4.5990(59) MHz, chi(ab) = 3.17(147) MHz, and chi(cc) = 3.1082(59) MHz. Coordinates of the H atom in the center-of-mass frame have been determined with use of the Kraitchman equations, /aH/ = 0.516 A and /bH/ = 1.171 A. The inertial defects of HOONO and DOONO are consistent with a planar equilibrium structure with significant out-of-plane H atom torsional motion. Comparisons of the present results are made to ab initio calculations.

Role of OH-stretch/torsion coupling and quantum yield effects in the first OH overtone spectrum of cis-cis HOONO.

A joint theoretical and experimental investigation is undertaken to study the effects of OH-stretch/HOON torsion coupling and of quantum yield on the previously reported first overtone action spectrum of cis-cis HOONO (peroxynitrous acid). The minimum energy path along the HOON dihedral angle is computed at the coupled cluster singles and doubles with perturbative triples level with correlation consistent polarized quadruple zeta basis set, at the structure optimized using the triple zeta basis set (CCSD(T)/cc-pVQZ//CCSD(T)/cc-pVTZ). The two-dimensional ab initio potential energy and dipole moment surfaces for cis-cis HOONO are calculated as functions of the HOON torsion and OH bond length about the minimum energy path at the CCSD(T)/cc-pVTZ and QCISD/AUG-cc-pVTZ (QCISD-quadratic configuration interaction with single and double excitation and AUG-augmented with diffuse functions) level of theory/basis, respectively. The OH-stretch vibration depends strongly on the torsional angle, and the torsional potential possesses a broad shelf at approximately 90 degrees , the cis-perp conformation. The calculated electronic energies and dipoles are fit to simple functional forms and absorption spectra in the region of the OH fundamental and first overtone are calculated from these surfaces. While the experimental and calculated spectra of the OH fundamental band are in good agreement, significant differences in the intensity patterns are observed between the calculated absorption spectrum and the measured action spectrum in the 2nu(OH) region. These differences are attributed to the fact that several of the experimentally accessible states do not have sufficient energy to dissociate to OH+NO(2) and therefore are not detectable in an action spectrum. Scaling of the intensities of transitions to these states, assuming D(0)=82.0 kJ/mol, is shown to produce a spectrum that is in good agreement with the measured action spectrum. Based on this agreement, we assign two of the features in the spectrum to Deltan=0 transitions (where n is the HOON torsion quantum number) that are blue shifted relative to the origin band, while the large peak near 7000 cm(-1) is assigned to a series of Deltan=+1 transitions, with predominant contributions from torsionally excited states with substantial cis-perp character. The direct absorption spectrum of cis-cis HOONO (6300-6850 cm(-1)) is recorded by cavity ringdown spectroscopy in a discharge flow cell. A single band of HOONO is observed at 6370 cm(-1) and is assigned as the origin of the first OH overtone of cis-cis HOONO. These results imply that the origin band is suppressed by over an order of magnitude in the action spectrum, due to a reduced quantum yield. The striking differences between absorption and action spectra are correctly predicted by the calculations.

Rotational spectrum of cis-cis HOONO.

The pure rotational spectrum of cis-cis peroxynitrous acid, HOONO, has been observed. Over 220 transitions, sampling states up to J'=67 and Ka'=31, have been fitted with an rms uncertainty of 48.4 kHz. The experimentally determined rotational constants agree well with ab initio values for the cis-cis conformer, a five-membered ring formed by intramolecular hydrogen bonding. The small, positive inertial defect Delta=0.075667(60) amu A2 and lack of any observable torsional splittings in the spectrum indicate that cis-cis HOONO exists in a well-defined planar structure at room temperature.

Cis-cis and trans-perp HOONO: action spectroscopy and isomerization kinetics.

The weakly bound HOONO product of the OH+NO2+M reaction is studied using the vibrational predissociation that follows excitation of the first OH overtone (2nu1). We observe formation of both cis-cis and trans-perp conformers of HOONO. The trans-perp HOONO 2nu1 band is observed under thermal (223-238 K) conditions at 6971 cm(-1). We assign the previously published (warmer temperature) HOONO spectrum to the 2nu1 band at 6365 cm(-1) and 2nu1-containing combination bands of the cis-cis conformer of HOONO. The band shape of the trans-perp HOONO spectrum is in excellent agreement with the predicted rotational contour based on previous experimental and theoretical results, but the apparent origin of the cis-cis HOONO spectrum at 6365 cm(-1) is featureless and significantly broader, suggesting more rapid intramolecular vibrational redistribution or predissociation in the latter isomer. The thermally less stable trans-perp HOONO isomerizes rapidly to cis-cis HOONO with an experimentally determined lifetime of 39 ms at 233 K at 13 hPa (in a buffer gas of predominantly Ar). The temperature dependence of the trans-perp HOONO lifetime in the range 223-238 K yields an isomerization barrier of 33+/-12 kJ/mol. New ab initio calculations of the structure and vibrational mode frequencies of the transition state perp-perp HOONO are performed using the coupled cluster singles and doubles with perturbative triples [CCSD(T)] model, using a correlation consistent polarized triple zeta basis set (cc-pVTZ). The energetics of cis-cis, trans-perp, and perp-perp HOONO are also calculated at this level [CCSD(T)/cc-pVTZ] and with a quadruple zeta basis set using the structure determined at the triple zeta basis set [CCSD(T)/cc-pVQZ//CCSD(T)/cc-pVTZ]. These calculations predict that the anti form of perp-perp HOONO has an energy of DeltaE0=42.4 kJ/mol above trans-perp HOONO, corresponding to an activation enthalpy of DeltaH298 (double dagger 0)=41.1 kJ/mol. These results are in good agreement with statistical simulations based on a model developed by Golden, Barker, and Lohr. The simulated isomerization rates match the observed decay rates when modeled with a trans-perp to cis-cis HOONO isomerization barrier of 40.8 kJ/mol and a strong collision model. The quantum yield of cis-cis HOONO dissociation to OH and NO2 is also calculated as a function of photon excitation energy in the range 3500-7500 cm(-1), assuming D0=83 kJ/mol. The quantum yield is predicted to vary from 0.15 to 1 over the observed spectrum at 298 K, leading to band intensities in the action spectrum that are highly temperature dependent; however, the observed relative band strengths in the cis-cis HOONO spectrum do not change substantially with temperature over the range 193-273 K. Semiempirical calculations of the oscillator strengths for 2nu1(cis-cis HOONO) and 2nu1(trans-perp HOONO) are performed using (1) a one-dimensional anharmonic model and (2) a Morse oscillator model for the OH stretch, and ab initio dipole moment functions calculated using Becke, Lee, Yang, and Parr density functional theory (B3LYP), Møller-Plesset pertubation theory truncated at the second and third order (MP2 and MP3), and quadratic configuration interaction theory using single and double excitations (QCISD). The QCISD level calculated ratio of 2nu1 oscillator strengths of trans-perp to cis-cis HOONO is 3.7:1. The observed intensities indicate that the concentration of trans-perp HOONO early in the OH+NO2 reaction is significantly greater than predicted by a Boltzmann distribution, consistent with statistical predictions of high initial yields of trans-perp HOONO from the OH+NO2+M reaction. In the atmosphere, trans-perp HOONO will isomerize nearly instantaneously to cis-cis HOONO. Loss of HOONO via photodissociation in the near-IR limits the lifetime of cis-cis HOONO during daylight to less than 45 h, other loss mechanisms will reduce the lifetime further.