Michelle
Färber
a,
Luc
Vereecken
a,
Hendrik
Fuchs
ab,
Georgios I.
Gkatzelis
a,
Franz
Rohrer
a,
Sergej
Wedel
a,
Andreas
Wahner
a and
Anna
Novelli
*a
aInstitute for Energy and Climate Research, IEK-8: Troposphere, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany. E-mail: m.faerber@fz-juelich.de; a.novelli@fz-juelich.de
bDepartment of Physics, University of Cologne, 50932 Cologne, Germany
First published on 17th January 2024
The formation of peroxynitrates (RO2NO2) from the reaction of peroxy radicals (RO2) and nitrogen dioxide (NO2) and their subsequent redissociation are typically not included in chemical mechanisms. This is often done to save computational time as the assumption is that the equilibrium is strongly towards the RO2 + NO2 reaction for most conditions. Exceptions are the reactions of the methyl peroxy radical due to its abundance in the atmosphere and of acyl-RO2 radicals due to the long lifetime of peroxyacyl nitrates RO2NO2 (PANs). In this study, the nighttime oxidation of cis-2-butene and trans-2-hexene in the presence of NO2 is investigated in the atmospheric simulation chamber SAPHIR, Forschungszentrum Jülich, Germany, at atmospherically-relevant conditions at different temperatures (≈276 K, ≈293 K, ≈305 K). Measured concentrations of peroxy and hydroperoxy radicals as well as other trace gases (ozone, NO2, volatile organic compounds) are compared to state-of-the-art zero-dimensional box model calculations. Good model-measurement agreement can only be achieved when reversible RO2 + NO2 reactions are included for all RO2 species using literature values available from the latest SAR by [Jenkin et al., Atmos. Chem. Phys., 2019, 19, 7691]. The good agreement observed gives confidence that the SAR, derived originally for aliphatic RO2, can be applied to a large range of substituted RO2 radicals, simplifying generalised implementation in chemical models. RO2NO2 concentrations from non-acyl RO2 radicals of up to 2 × 10 cm−3 are predicted at 276 K, impacting effectively the kinetics of RO2 radicals. Under these conditions, peroxy radicals are slowly regenerated downwind of the pollution source and may be lost in the atmosphere through deposition of RO2NO2. Based on this study, 60% of RO2 radicals would be stored as RO2NO2 at a temperature of 10 °C and in the presence of a few ppbv of NO2. The fraction increases further at colder temperatures and/or higher NO2 mixing ratios. This does not only affect the expected concentrations of RO2 radicals but, as the peroxynitrates can react with OH radicals or photolyse, they could comprise a net sink for RO2 radicals as well as increase the production of NOx (= NO + NO2) in different locations depending on their lifetime. Omitting this chemistry from the kinetic model can lead to misinterpreted product formation and may prevent reconciling observations and model predictions.
NO2 + O3 → NO3+ O2, | (1) |
In recent field studies utilising the laser-induced-fluorescence (LIF) technique for radical measurements, discrepancies between measured and modelled RO2 radical concentrations were found at night at moderate O3 (∼(1–70) ppbv), NO2 (∼(3–50) ppbv), and low NO (∼(0–0.3) ppbv) mixing ratios in the vicinity of Beijing,5,6 China, Wangdu,7 China, and London,8 UK. Just recently, two studies focussing on investigating the oxidation of isoprene9 and anthropogenic alkenes10 by NO3 resulted in new detailed chemical mechanisms for these compounds as intricate as their oxidation by OH. These studies highlighted the inability of detecting a large fraction of NO3 -containing isoprene-RO29 as well as NO3–RO2 from short-chain alkenes.10 Both the isomerisation (isoprene)9 and the decomposition (short-chain alkenes)10 of the NO3 -alkoxy radical do not lead to the formation of HO2 or OH, which is required for the RO2 measurement by the LIF technique:
![]() | (2) |
However, for long-chain alkenes such as for example trans-2-hexene, all RO2 are detectable.10 The discrepancies found at nighttime in the field5–8 may be partly explainable by the missing detection of the LIF instrument of single short-chain alkenes, contained in the measured air mixture. Still large discrepancies remain between measured and modelled RO2 radicals for high NO (up to ∼100 ppbv) and NO2 (up to ∼70 ppbv).5–8,11
Despite improvements in the agreement between measured and modelled RO2 radicals observed in the experiments in the SAPHIR chamber,10 discrepancies are still observed in particular right after the injection of the VOC when the RO2 radical production is the highest, as well as in the time dependence of the RO2 radical concentrations. The conditions in the chamber were chosen to facilitate the formation of NO3 radicals, and were characterised by zero NO and medium NO2 (17–40 ppbv) and ozone (7–30 ppbv). For these conditions, RO2 radicals formed in the oxidation by either O3 or NO3 are assumed to be mainly lost by their reaction with HO2 or RO2. NO2 will be the dominant reaction partner for these RO2 radicals, producing a short-lived peroxynitrate (RO2NO2) which is in a very fast equilibrium (τRO2NO2 ≈ 0.2 s at 298 K) with its decomposition in the lower troposphere:
![]() | (3) |
The impact of eqn (3) on the RO2 radical and the formed RO2NO2 is thought to be negligible due to the instability of RO2NO2 in lower tropospheric oxidation processes, therefore, the reaction of RO2 and NO2 forming RO2NO2 is omitted in most atmospheric models for non-acyl RO2 radicals except for CH3O2.12 For the latter, the formation of alkyl peroxynitrates has been thought to be mostly relevant at cold temperatures in polar regions,13,14 in the upper troposphere, where Browne et al.14 found that including the methyl peroxynitrate (CH3O2NO2) chemistry is relevant at temperatures below 240 K, or in biomass burning plumes that are lofted to high altitudes.14 In contrast, the reactions of acyl-RO2 radicals with NO2 are typically implemented in all chemical models due to the long lifetime (∼40–45 min at 298 K) of the corresponding peroxyacyl nitrate RO2NO2 (PANs) and, consequently, their relevance on regional and even global scales.13
A recent study by Khan et al.15 investigated the global effect of the reversible formation of RO2NO2 on NOx, i.e. NO + NO2, OH, and O3, using the generic rate coefficients from Jenkin et al.12 for more than 40 non-acyl RO2 radicals. Loss rates for the formed RO2NO2 species via reaction of OH and photolysis were also included. In their study, it was shown that even on the ground, up to 25% more NOx would be expected at the equator as photolysis and reaction with OH accelerate the consumption of RO2NO2. This indicates that their role and importance might need to be reevaluated.
A measure for the ratio of RO2 and RO2NO2 concentrations at equilibrium is given by the equilibrium constant K, the ratio of the forward and backward reaction rate constants kRO2+NO2 and kRO2NO2, respectively:
K = kRO2+NO2/kRO2NO2. | (4) |
For the alkyl peroxynitrates CH3O2NO2 and ethyl peroxynitrate (C2H5O2NO2), the equilibrium constant is 7 × 10−11 cm3 at 276 K, respectively (Table 1). In contrast, peroxyacetyl nitrate (CH3C(O)O2) has an equilibrium constant of 1 × 10−8 cm3 due to its longer lifetime. While at high temperatures formation of alkyl peroxynitrates is not important, the equilibrium shifts more towards the RO2NO2 products at lower temperatures, making it worthwhile to evaluate the impact of RO2NO2 formation during winter conditions.
RO2 or RO2NO2 | k 0 | k ∞ | F c |
k
276![]() |
k
292![]() ![]() |
k
305![]() ![]() |
---|---|---|---|---|---|---|
a Recommended by IUPAC.16 b Determined from isomeric mixtures formed from the reaction of chlorine with butane17,18 and recommended by Jenkin et al.12 c Determined from isomeric mixtures formed from the reaction of chlorine with hexane18 and recommended by Jenkin et al.12 d Pressure-independent generic rate coefficient recommended by Jenkin et al.12 | ||||||
Forward reaction, kRO2+NO2 [cm3 s−1] | ||||||
CH3O2a | 1.2 × 10−30 (T/300)−6.9[M] | 1.8 × 10−11 | 0.36 | 7.1 × 10−12 | 6.2 × 10−12 | 5.5 × 10−12 |
C2H5O2a | 1.3 × 10−29 (T/300)−6.2 [M] | 8.8 × 10−12 | 0.31 | 5.6 × 10−12 | 5.3 × 10−12 | 5.0 × 10−12 |
n- and sec- C4H9O2b | 9.6 × 10−12 | 9.6 × 10−12 | 9.6 × 10−12 | 9.6 × 10−12 | ||
RO2d | 9.0 × 10−12 (=kfPN) | 9.0 × 10−12 | 9.0 × 10−12 | 9.0 × 10−12 | ||
Reverse reaction, kRO2NO2 [s−1] | ||||||
CH3O2NO2a | 9.0 × 10−5exp (−9690/T)[M] | 1.1 × 1016 exp (−10![]() |
0.36 | 0.1 | 0.7 | 3.2 |
C2H5O2NO2a | 4.8 × 10−4exp (−9285/T)[M] | 8.8 × 1015 exp (−10![]() |
0.31 | 0.2 | 1.7 | 7.5 |
n- and sec- C4H9O2NO2b | 8.3 × 1015 exp (−10![]() |
0.4 | 3.2 | 14.3 | ||
C 6H13O2NO2 isomersc | 7.5 × 1015 exp (−10![]() |
0.4 | 2.8 | 12.9 | ||
RO2NO2d | 7.6 × 1015 exp (−10![]() |
0.3 | 2.6 | 11.8 |
In this work, the impact of the RO2 + NO2 reaction on modelled RO2 radical concentrations is studied for RO2 radicals from the oxidation of cis-2-butene and trans-2-hexene by O3 and NO3 at different temperatures (≈276 K, ≈293 K, ≈305 K). Experiments were performed in the atmospheric simulation chamber SAPHIR at Forschungszentrum Jülich, Germany. Nighttime conditions with high NO2 (>20 ppbv) and moderate O3 (∼10 ppbv) mixing ratios were tested which facilitate the formation of NO3 and mimic conditions often found at night. Concentrations of HO2 and RO2 radicals, as well as of O3, NO2, and VOCs were measured and are compared to zero-dimensional box model calculations, incorporating state-of-the-art oxidation schemes. The reaction of non-acyl RO2 radicals with NO2 forming RO2NO2 and its backward reaction are implemented in the chemical mechanisms of cis-2-butene and trans-2-hexene by using literature values,12,16–18 and the implications of the findings on the nighttime atmospheric RO2 chemistry are discussed. The newly introduced ozonolysis scheme and RO2 isomerisation reactions, suggested by Novelli et al.10 for cis-2-butene and trans-2-hexene, as well as their temperature dependence are tested and model-measurement comparisons are used to improve the chemical mechanisms for these species.
Nighttime experiments for cis-2-butene and trans-2-hexene were performed at different temperatures. The nighttime oxidation of cis-2-butene was studied at two different temperatures (T ≈ 276 K (cold), 295 K (medium)). The nighttime chemistry of trans-2-hexene was investigated at three temperatures (T ≈ 276 K (cold), 292 K (medium), 305 K (hot)). At low temperatures (T ≈ 280 K), a pure ozonolysis experiment in the absence of NO2 was also performed for trans-2-hexene.
The experimental procedure of the experiments is shown in Fig. 1–3. Before each experiment, the chamber was flushed to remove trace gases from the previous experiment. No detectable OH reactivity, equivalent to the inverse lifetime of OH, was observed in the dark, clean, and dry chamber. Cis-2-Butene (Air Liquide, 1% in N2, purity 99.4%) or trans-2-hexene (Sigma Aldrich, purity 97%) were first injected in the clean, dry, and dark chamber, reaching mixing ratios close to atmospheric conditions (<6 ppbv). Afterwards, 170–210 ppmv of carbon monoxide (CO, Air Liquide, purity 99.997%) was added to the chamber acting as OH scavenger by converting OH into HO2 to avoid its reaction with cis-2-butene or trans-2-hexene. In one oxidation experiment of trans-2-hexene, 142 ppmv of methane (CH4, Air Liquide, purity 99.5%) were injected instead to aim for lower HO2 radical concentrations. Since more reactions are involved to form HO2, a slower production of HO2 can be achieved, and higher RO2 radical concentrations are expected due to the formation of CH3O2. Nitrogen dioxide (Linde, 500 ppmv in N2, purity 99.991%) of up to 35 ppbv was then injected in all experiments with the exception of the ozonolysis study. Afterwards, ozone generated by a silent discharge ozoniser (O3onia) was injected to reach mixing ratios between 12 ppbv and 33 ppbv (96 ppbv in the ozonolysis experiment), initiating the production of NO3 in the reaction with NO2 (eqn (1)). Due to the fast reaction of cis-2-butene and trans-2-hexene with O3 and NO3, they are consumed on a timescale of ∼2–5 hours, and the alkene was injected a second time and third time (experiments with cis-2-butene at 295 K and trans-2-hexene at 305 K). To boost the oxidation by NO3, NO2 was injected a second time as well. In all experiments, no NO was present in the chamber and OH radical concentrations were below the detection limit, leading to a negligible contribution of OH chemistry. All experiments were conducted in dry conditions (H2O < 0.07 ppmv) and in complete darkness.
The OH reactivity (kOH) was measured by a pump-and-probe technique utilising LIF to detect the amount of OH reacting with the sampled air in a flow tube.30,31 Time series of cis-2-butene and trans-2-hexene and calibrated time series of acetaldehyde were measured by proton-transfer-reaction time-of-flight mass spectrometry (PTR-ToF-MS, Ionicon).32,33 Normalised counts of cis-2-butene and trans-2-hexene are converted to parts per billion (ppb) by using the observed OH reactivity at the point in time of the VOC injection. Ozone was detected by UV absorption (Ansyco), and CO, CH4, formaldehyde, and water vapour were monitored utilising cavity ring-down spectroscopy (CRDS, Picarro). Furthermore, NO and NO2 were measured by a chemilumiscence instrument (EcoPhysics). More detailed information about the performance of the instruments can be found in the ESI‡ Section D. NO3 and N2O5 measurements were not available for any experiment, therefore corresponding wall loss rates were introduced and adjusted to match the observed VOC decay.
• Ozonolysis scheme for cis-2-butene and trans-2-hexene including the formation of OCHCH2OO˙ (ethanal-2-peroxy, MCM notation: HCOCH2O2) and O
CHCHO(OO˙)CH2CH3 (butanal-2-peroxy, MCM notation: BUTALAO2). The subsequent bimolecular chemistry for these radicals follows the MCM.
• Isomerisation reactions for ethanal-2-peroxy and butanal-2-peroxy based on quantum chemical calculations, leading to the formation of OC(OOH)CH2OO˙ (ethyl-1-peracid-2-peroxy) and O
C(OOH)CH(OO˙)C2H5 (butyl-1-peracid-2-peroxy), respectively (Fig. S1, ESI‡).
• Updated alkoxy decomposition rates for RO radicals formed in the oxidation of cis-2-butene and trans-2-hexene by NO3 based on quantum chemical calculations.
• Optimised yield of the RO and OH radical from the reaction of RO2 with HO2 for the first generation NO3–RO2 (MCM notation: C62NO33O2, C63NO32O2) formed in the oxidation of trans-2-hexene by NO3.
The follow-up chemistry for additionally included RO2 and RO radicals was implemented following the SARs in Jenkin et al.,12 Vereecken et al.,10,36–38 and Novelli et al.10 For the model-measurement comparison when NO3 chemistry is contributing, the fraction of RO2 which is not measurable by the radical instrument10 (Section 2.2) must be taken into account, i.e. distinguishing between “total RO2” and “detectable RO2”.
The effect of the isomerisation reactions of the RO2 radicals from the ozonolysis reaction, ethanal-2-peroxy and butanal-2-peroxy, on the predicted RO2 speciation is investigated by comparing model runs using the FZJ mechanisms with and without including RO2 isomerisation reactions. In another sensitivity model run, not only reactions of methyl peroxy (CH3O2) and PAN-like RO2 (RC(O)OO˙) with NO2 are included, as implemented in the FZJ mechanism, but also the reaction of NO2 with all formed non-acyl RO2 species. In this sensitivity run, RO2 + NO2 reactions forming RO2NO2 as well as the corresponding backward reactions (eqn (3)) are introduced for all formed RO2 species. Rate constants for the forward and backward reactions for CH3O2, ethyl peroxy (C2H5O2), and acetonyl peroxy (CH3 C(
O)CH2O2) are taken from recommendations by IUPAC.16 Recommendations in Jenkin et al.,12 partially based on Zabel et al.,18 are used for the forward and backward reaction rate constants for C4 peroxy radicals, and for the backward reaction rate constant for C6–RO2NO2 in this work. For all RO2 for which rate constants cannot be taken from literature, such as for n-C3H7O2, a pressure-independent generic forward reaction rate constant, kfPN, and a rounded average of reported backward reaction rate constants for C2 to C8 alkyl RO2, kbPN, recommended by Jenkin et al.,12 are used. However, the contribution of n-C3H7O2 to the total RO2 amounts to 2% and is thus negligible. It is important to mention that all the rate coefficients used from literature and SAR are for non-NO3 substituted RO2 radicals but are used here due to the lack of specific values for nitrate RO2 radicals. An overview of temperature-dependent rate coefficients for the formation and decomposition of non-acyl RO2NO2 is shown in Table 1, together with the corresponding rate coefficients at T = 276 K, 292 K, and 305 K. We also determined the kbPN and kfPN values by a fitting procedure against the experimental data, retrieving values that are within 30% of the literature data above. This supports applying literature data for alkyl-RO2 to our NO3-RO2 radicals, but the fitted values carry a large uncertainty and are not used in the results and discussion below
For cold conditions, there is a large difference in the total modelled RO2 between the FZJ mechanisms with and without including the additional formation of RO2NO2, with the latter predicting a much lower concentration (by a factor of 7 on average) in the first two hours of oxidation and a much slower increase. When the detectable modelled RO2 are compared to the measured RO2 radical concentrations at low temperatures, the mechanism including non-acyl RO2NO2 formation (FZJ + RO2NO2 model) reproduces the data very well with a model-to-measurement ratio of 0.84, and excellent reproduction of the time dependence. The mechanism without including RO2 + NO2 reactions (FZJ model) can reasonably reproduce the detectable RO2 after one hour after the cis-2-butene injections (average model-to-measurement ratio = 0.8). However, their concentration is largely overestimated by up to a factor of 4.5 right after the injection, with a distinctly earlier concentration peak. The smaller difference between the two model results observed for the second injection of cis-2-butene is due to a larger contribution of CH3O2 (from close to 0% after the first cis-2-butene injection to ∼ 20% after the second). As the reaction of CH3O2 with NO2 is included in both the FZJ and MCM mechanisms, the effect of the correction is smaller.
Differences of total RO2 radical concentrations between the two FZJ mechanisms are also observed in the experiment at medium temperatures. However, the difference is much less pronounced and the average ratio of concentrations obtained between the two FZJ mechanisms is 1.1. A small underestimation of the measured RO2 (by a factor of 1.2) by both mechanisms can be observed by comparing it to the detectable RO2. The FZJ mechanism without alkyl-RO2 + NO2 reactions would predict a small peak of the detectable RO2, while concentrations are better described if these reactions are considered (Fig. 1). From the second injection on, RO2 radical concentrations predicted by both mechanisms converge to the same value.
The measured HO2 radical concentrations are well described by both FZJ mechanisms and agree within 36% (FZJ model) and 21% (FZJ + RO2NO2 model) after the first VOC injection and within 5% (both FZJ models) after the second VOC injection in the experiment at cold conditions and within 7% (both FZJ mechanisms) in the experiment at medium temperatures. A comparison of the concentration of acetaldehyde, a major product of the ozonolysis of cis-2-butene, obtained by the different models is shown in Fig. S11 (ESI‡) for the experiment at cold conditions. The different mechanisms predict similar product concentrations with a model-to-measurement ratio of 0.7.
The addition of RO2 + NO2 reactions leads to a less steep increase of the RO2 radical concentrations which is in much better agreement with the measurements after the first injection of trans-2-hexene (within 2% compared to 20% for the FZJ model). Similar to cis-2-butene, at 276 K the FZJ mechanism expects a fast increase of RO2 radicals, which is not observed in the measurements.
As ozone contributes up to 40% to the oxidation of trans-2-hexene, and discrepancies were observed earlier when ozone dominated the oxidation at 292 K in the study by Novelli et al.,10 an ozonolysis experiment was performed for trans-2-hexene. Fig. S2 (ESI‡) shows the comparison between measured and modelled HO2 and RO2 radicals for the MCM and the FZJ mechanisms. An agreement within 20% and 30% between measured and modelled HO2 is found for the MCM and the FZJ mechanism, respectively. Although the FZJ mechanism improves the agreement between measured and modelled RO2 radicals, still a ratio of modelled to measured of 1.7 remains at the point of injection of trans-2-hexene. The impact of this discrepancy, in the presence of NO2, was estimated by assuming an overall overestimation of a factor of 1.7 for RO2 formed from ozonolysis and correcting their concentrations modelled by the FZJ mechanism including RO2NO2 formation by this factor. As a consequence, the total modelled RO2 radical concentrations would be 25% lower (cold temperatures) at the peak RO2 after the first injection (Fig. 2) with an average ratio of modelled to measured RO2 of ≈0.7. Although the agreement between measured and modelled RO2 is decreased when an overestimation of RO2 of 1.7 is assumed (from 12% to 37% for the first 2.5 hours of oxidation), the modelled RO2 shows the same behaviour in reaching the peak RO2 radical concentrations and a worse agreement is found for the mechanism without additional NO2 reactions (model-to-measurement ratio of 1.9 for the first 2.5 hours of oxidation). More details about the interpretation of the ozonolysis experiment can be found in the ESI‡ Section A.
As compared to cis-2-butene (Fig. 1) and to the ozonolysis experiment (Fig. S2, ESI‡), the HO2 radicals are largely underestimated by both models. HO2 radicals are mainly formed from the reaction of CO with OH (≈61% of the total production rate) which is formed directly from the ozonolysis of trans-2-hexene, among other production pathways such as isomerisation reactions of second-generation RO2 (30% of the total production rate). A similar discrepancy between measured and modelled HO2 radicals was observed when performing the same experiment but using CH4 instead of CO as an OH scavenger (Fig. 3).
The observed model-measurement discrepancy could be due to an artefact in the LIF instrument detecting HO2, if there was a mechanism artificially producing HO2 in the HOx detection cell. An instrumental interference by NO2 and peroxynitric acid (HNO4, formed in the reaction of HO2 with NO2) can be excluded because this would also impact the HO2 radical concentration measured in the experiment with cis-2-butene, for which no model-measurement discrepancy for HO2 was observed (Fig. 1). Furthermore, an interference from the contemporary decomposition of RO2 radicals is minimised by running the HOx detection cell at reduced NO (Section 2.2) and no indication of an interference was observed when changing the NO concentrations in the HOx cell. A clear correlation of the model-measurement difference in HO2 radical concentrations is only found with the modelled non-acyl RO2NO2 concentrations (Fig. 4), indicating that the model-measurement discrepancy of HO2 is likely due to an instrumental interference in the HO2 measurements from RO2NO2. This would be consistent with the good agreement observed for the ozonolysis experiment (Fig. S2, ESI‡). An interference of ≈2% of the non-acyl RO2NO2 in the HOx cell would be enough to explain the measured HO2 radical concentration (Fig. 2 and 3).
![]() | ||
Fig. 4 Correlation between observed model-measurement deviations of HO2 with modelled non-acyl RO2NO2 concentrations. Displayed are NO3 oxidation experiments of trans-2-hexene. |
It is currently not clear how the non-acyl RO2NO2 would be detected as HO2 in the HOx cell of the LIF instrument, especially as no discrepancy was observed for cis-2-butene. RO2 radicals are obtained by subtracting the measured HO2 and OH radicals from the sum of all three radical species (OH + HO2 + RO2), as detected in the ROx LIF system. Therefore, an interference in the HOx measurement would result in a lower measured RO2 radical concentration as compared to the “real” value. Using the modelled HO2 radical instead of the measured value would increase the RO2 radical concentrations by 20% and 30% in the experiments at low temperatures with CO (Fig. S15, ESI‡) and CH4 (Fig. 3) as OH scavenger, respectively. As the HO2 radical concentrations for the conditions of the experiments in this study are low, their impact is limited and no change in the time profile of the RO2 radical is observed. More details can be found in Section G of the ESI.‡
Two more experiments at 292 K and 305 K were performed with trans-2-hexene (Fig. S12, ESI‡). Although a worse agreement between measured and modelled RO2 was observed (average model-measurement ratio of 1.5 and 1.2 for 292 K and 305 K, respectively) as compared with the experiment at 276 K, a shift in the peak for measured RO2 radical concentration could be observed which does not show any “delay” as found in the experiment at cold temperatures (Fig. 2). This is consistent with the faster decomposition rate of the formed RO2NO2 at those temperatures and is also reproduced by both models where negligible differences are seen already maximum 30 minutes after the VOC injection.
As compared to the experiment at 276 K (cold), a better agreement between modelled and measured HO2 radical is obtained. Model calculations tend to overestimate the HO2 radicals at 276 K. Just recently, a study by McKee et al.39 showed that at high temperatures production of nitryl hydride (HNO2) from the reaction of HO2 and NO2 can be expected. By including this reaction in the model calculations of the experiment at hot conditions (T = 305 K), the comparison of the model with the measurement improves leading to a model-to-measurement ratio of 1.3 on average (compared to 1.6 if the HNO2 formation is not considered, Fig. S12, ESI‡). More discussion about possible reasons for the observed discrepancies found at 292 K and 305 K for trans-2-hexene can be found in the ESI‡ Section G.
The comparison of acetaldehyde concentrations predicted by the mechanisms is shown in Fig. S14 (ESI‡), showing similar concentrations for the different temperatures. The mechanisms overall underestimate the measured acetaldehyde concentrations by a factor of 1.7 and 3 in the experiments at medium and hot temperatures, respectively, but in the experiment at cold conditions a good agreement is achieved (within 15%). The reason for this temperature-dependent discrepancy cannot be easily explained and needs further investigations.
The RO2NO2 concentrations expected from the FZJ mechanism are systematically lower than the RO2NO2 concentrations expected from the FZJ + RO2NO2 mechanism, including non-acyl RO2NO2 formation, as only the formation of CH3O2NO2 is implemented (Section 2.3). A strong increase of peak RO2NO2 concentrations from 140 pptv to up to 480 pptv is observed over the temperature range of 295–276 K. While the main loss path for RO2 radicals in the FZJ model is the loss in the reaction with HO2 (k∼ (2 − 10) × 10−3 s−1), the reaction with NO2 (k∼ (5 − 9) s−1, Table 1) is the predominant, temporary loss of RO2 in the FZJ + RO2NO2 mechanism, including formation of non-acyl RO2NO2, due to the large NO2 mixing ratios (up to 40 ppbv).
The fraction of RO2 stored in the reservoir RO2NO2 species can be derived from the ratio of corresponding concentrations: RO2NO2/(RO2 + RO2NO2) (Table 2, Fig. 5 and Fig. S13, ESI‡). In the experiment at cold conditions, the conversion of RO2via reaction with NO2 leads to more than 90% of non-acyl RO2 that is present in form of the reservoir species RO2NO2 throughout the experiment. In comparison, 80% to 40% of the non-acyl RO2 are stored as RO2NO2 throughout the experiment at medium temperatures (Fig. 5), in which the temperature increased from 290 K to 300 K, leading to a shift of equilibrium towards RO2. When one of the VOCs is injected a large fraction of the RO2 radicals formed is converted into the reservoir species RO2NO2, leading to the slower increase of the free, measurable RO2 radical concentrations in the FZJ + RO2NO2 mechanism, including additional RO2NO2 formation, compared to the FZJ mechanism (Fig. 1 and 2). This emphasises the importance of including the formation of peroxynitrates in chemical mechanisms, especially at lower temperatures where the RO2NO2/RO2 ratio can reach a factor of 20.
Temperature [K] | K|eqn (3) [cm3] |
---|---|
305 | 7.0 × 10−13 |
292 | 3.2 × 10−12 |
276 | 2.5 × 10−11 |
At temperatures typically reached in mid latitudes, RO2NO2 concentrations are expected to be on the order of 1 × 1010 cm−3 at the conditions tested in this work and can therefore impact the fate of ROx species. In cities as Beijing,5 China, temperatures around 255 K can be reached in winter, leading to a significant fraction of non-acyl RO2 stored as reservoir RO2NO2 species, and thus to different RO2 radical concentrations. Low temperatures are also present at high altitudes, where biomass burning plumes can introduce polluted conditions, allowing alkyl-RO2NO2 other than CH3O2NO2 to be formed and impacting the effective RO2 radical concentration.
One important finding of this study is that the rate coefficients used for forward and backward reactions (Table 1 and 2) do an excellent job in bringing measurements and model results in agreement despite not being derived for the NO3-substituted RO2 radicals as investigated here. This indicates that the substituent on the RO2 radical does not seem to have an impact on the reaction rate. This allows to generalise the finding beyond the molecules and conditions explored in the SAPHIR chamber.
Fig. 6 shows the fraction of RO2NO2/(RO2 +RO2NO2) for different temperatures and NO2 mixing ratios. Conditions observed in different cities5,6,8,40 are highlighted and are found to favour the formation of non-acyl RO2NO2. Especially at wintertime in Beijing, a significant amount of RO2 is expected to be stored as RO2NO2. In contrast, NO2 and temperature observed in September in a boreal forest2 do not facilitate the production of RO2NO2. Though, for relatively low values of NO2 of ∼1 ppbv and ∼10 °C (280 K) that can also be found in forested environments, more than 20% of RO2 radical is stored as RO2NO2. For a mixing ratio of NO2 up to 10 ppbv, the fraction can increase up to 80%. Although a large part of the RO2NO2 will decompose back to RO2 and NO2, depending on temperature and lifetime, RO2NO2 can be transported and contribute to the NOx levels further away from their emission sources. In addition, nothing is known about the additional loss rate of RO2NO2 by reaction with OH and photolysis. Khan et al.15 introduced these loss reactions based on similarity with other molecules and show they could have an impact in particular at the equator.
![]() | ||
Fig. 6 Fraction of peroxy radicals stored as RO2NO2 as a function of temperature and NO2 mixing ratios. Contour lines mark certain RO2NO2/(RO2 + RO2NO2) ratios. Values were determined from a steady-state calculation, assuming the generic forward and backward reaction rate coefficients kfPN and kbPN (Table 1) for reaction (3), respectively. White stars mark regimes that were observed in a boreal forest,2 in London,8 Pasadena,40 and Beijing at cold5 and medium temperatures.6 |
The measured time profiles of the RO2 radical concentrations are reproduced best when accounting for the reaction of RO2 with NO2 for all RO2, not solely those for CH3O2 and acyl-RO2 radicals as typically included in atmospheric models. Especially for cold conditions, including the formation of reservoir RO2NO2 species significantly improves the modelled RO2 and an average agreement with the measured RO2 within 16% is achieved for cis-2-butene. This reaction allows reversible formation of alkyl peroxynitrate reservoir species, leading to a reduction of the effective concentration of free RO2 radicals that is especially pronounced at lower temperatures. In addition, the RO2 radical concentrations observed are equally (within 20% on average) described by both FZJ mechanisms (FZJ model with and without including additional RO2NO2 formation) above 290 K. With decreasing temperatures, including the formation and subsequent decomposition of RO2NO2 in the model reduces the predicted total free RO2 radical concentration and delays peak RO2 radical concentrations as observed in the experiments, thereby significantly improving the shape of modelled RO2 time profiles.
A discrepancy between modelled and measured HO2 radical concentration, which for the low-temperature conditions reached a factor of 3 on average, was observed for trans-2-hexene. Sensitivity analyses indicate that this discrepancy correlates best with a measurement interference by non-acyl RO2NO2, where an interference of ≈2% would be sufficient to reproduce the measured HO2. However, as this discrepancy does not influence the time-dependence of the RO2 radical concentration profiles in the chamber experiments, it has no impact on the conclusions.
Under the conditions of the experiments in this work, up to 2 × 10 cm−3 of non-acyl RO2NO2 are expected at 276 K, comprising about 95% of RO2. Under winter conditions in polluted cities such as Beijing, China, where temperatures of 255 K are reached,5 even higher fractions of the RO2 present as RO2NO2 reservoir species are expected. This will have consequences for radical chemistry, the spatial distribution of RO2 and NO2, and the loss of radicals through deposition of RO2NO2.13
The general forward and backward reaction rates from SAR12 which refer to non-substituted RO2 radicals do an excellent job in bringing measured and modelled RO2 radicals in agreement for the conditions of this study (NO3–RO2). This suggests that substituents on the RO2 radical have a small effect on the reaction rates making the implementation in global models easier. Although the largest impact of these reactions can be expected for cold temperatures (higher altitudes), already for NO2 mixing ratios as low as few ppbv and a temperature of 280 K, a fraction of RO2NO2/(RO2 + RO2NO2) of 20% can be expected (Fig. 6). Given the large uncertainties of the follow-up chemistry of RO2NO2, it is not easy to assess, a priori, their effect in different environments and more studies are needed to shed light on their chemistry. Generally, the impact of the formation of RO2NO2 manifests likely as a time delay in the response of the chemical mixture to changes in VOC emission and oxidation, NOx levels, temperature, and other environmental factors.
Discrepancies remain between modelled and measured RO2 radicals in the ozonolysis experiment of trans-2-hexene and for NO2 + O3 experiments of trans-2-hexene at higher temperatures (292 K and 305 K). Uncertainties in the ozonolysis may be related to the chemistry of peracid-substituted RO2 radicals formed in RO2H-migration reactions. Insufficient data are available to resolve this issue at this time, but it is shown that the discrepancy is too small to significantly affect the O3 + NO3 nighttime experiments in this work.
Footnotes |
† The data of the experiments in the SAPHIR chamber used in this work are available on the EUROCHAMP data home page (https://data.eurochamp.org/). More information can be found in the Data availability statement. |
‡ Electronic supplementary information (ESI) available: The supplement related to this article is available online, and contains detailed information about the ozonolysis experiment of trans-2-hexene, the NO3 interference of the radical instrument, the modified alkoxy decomposition rate used for cis-2-butene-RO, instrumentation details, contributions of NO3 and O3 to the VOC oxidation, the time series of measured acetaldehyde concentrations, and the discussion of the nighttime oxidation of trans-2-hexene at medium and hot temperatures. See DOI: https://doi.org/10.1039/d3cp04163h |
This journal is © the Owner Societies 2024 |