1 Introduction
Atmospheric submicron aerosols have an adverse effect on air quality, human health, and climate (Jimenez et al., 2009; Pöschl, 2005). Secondary organic aerosol (SOA), which is formed from oxidation of volatile organic compounds (VOCs) followed by gas-to-particle partitioning, comprises a large fraction (20 %–90 %) of the submicron aerosol mass (Jimenez et al., 2009; Zhang et al., 2007). It is confirmed that a significant proportion of SOA arises from biogenic VOC (BVOC) oxidation (Hallquist et al., 2009; Spracklen et al., 2011).
Isoprene is globally the most abundant non-methane volatile organic compound originating from vegetation, with emissions estimated to be 440–660 Tg yr (Guenther et al., 2012). Due to its high abundance, as well as its high reactivity with atmospheric oxidants, isoprene plays a significant role in tropospheric chemistry, and its chemistry affects the global aerosol burden and distribution (Carlton et al., 2009; Fry et al., 2018; Ng et al., 2008, 2017; Surratt et al., 2010), although its SOA yield is much lower than those of monoterpenes and sesquiterpenes (Friedman and Farmer, 2018; Kim et al., 2015; Marais et al., 2016; McFiggans, et al., 2019; Mutzel et al., 2016; Ng et al., 2007, 2008; Surratt et al., 2010; Thornton et al., 2020). Recent model simulations suggested the isoprene-derived SOA production is 56.7 Tg C yr, contributing up to 41 % of global SOA (Stadtler et al., 2018). Observations in the southeastern United States suggested that isoprene-derived SOA makes up 17 %–48 % of total organic aerosol (Hu et al., 2015; Kim et al., 2015; Marais et al., 2016). As a consequence, it is essential to fully characterize the potential of isoprene to form condensable products and its contribution to SOA formation (Carlton et al., 2009).
Although the majority of isoprene emissions is emitted by plants and is light-dependent, isoprene emitted in the day can persist in the boundary layer after sunset, and its mixing ratio can remain as high as several parts per billion (ppb) (Brown et al., 2009; Starn et al., 1998; Stroud et al., 2002; Warneke et al., 2004). During the daytime, isoprene is primarily oxidized by the hydroxyl radical (OH) and somewhat by ozone (O), but its main oxidizers shift to nitrate radical (NO) and O in the nighttime (Wennberg et al., 2018). Due to the higher reactivity of NO with isoprene ( 6.5 10 cm molecules s and 1.28 10 cm molecules s at 298 K, respectively, IUPAC), a considerable fraction of the residual isoprene would be oxidized by NO at night, and therefore nocturnal nitrate radical chemistry is typically thought to be of significant importance for isoprene, especially in regions where sufficient nitrogen oxides are available (Brown et al., 2009; Fry et al., 2018; Ng et al., 2017; Wennberg et al., 2018). Although reaction with NO only represents 5 %–6 % of isoprene loss, it accounts for a large proportion of organic nitrates derived from isoprene oxidation ( 40 %–50 %) (Wennberg et al., 2018). Therefore, reaction of isoprene with NO is a potential source of SOA. In addition, it is found from both laboratory and chamber experiments that the SOA yield of isoprene from NO oxidation is higher than that from OH or O oxidation, which is typically less than 5 % (Carlton et al., 2009; Dommen et al., 2009; Kleindienst et al., 2007; Kroll et al., 2006). For example, Ng et al. (2008) concluded the isoprene SOA yield from NO was in the range of 4.3 %–23.8 %, depending on RO fate (higher SOA yield when the experiments were dominated by RO RO rather than RO NO reaction). Rollins et al. (2009) also observed a high SOA yield from isoprene (14 %) when both of its double bonds were oxidized by NO. In an aircraft study in the southeastern United States, Fry et al. (2018) derived an isoprene-NO SOA yield as large as 27 % on average under high-NO conditions, although their mass yield estimation was indirect, and based on a molar yield determination of 9 5 %. In light of the relatively high SOA yield from NO oxidation, even though only a minor fraction of isoprene is oxidized by NO, the SOA formed at nighttime would still probably be comparable to that produced at daytime (Brown et al., 2009; Fry et al., 2018).
However, isoprene-NO chemistry (Wennberg et al., 2018; Vereecken et al., 2021) has received less attention than the extensively studied OH- or O-initiated oxidation (Barber et al., 2018; Novelli et al., 2020; Peeters et al., 2014; Wang et al., 2018; Wennberg et al., 2018; Whalley et al., 2012). It has been recognized that later-generation oxidation of isoprene by NO makes a more significant contribution to SOA formation (Carlton et al., 2009; Fry et al., 2018; Rollins et al., 2009). Nevertheless, although the importance of multi-generation NO oxidation of isoprene to SOA formation has been recognized, few studies extended the investigation beyond the first-generation oxidation, and details of isoprene-NO multi-generation chemistry are still lacking.
Organic compounds, especially highly oxygenated organic molecules (HOMs) that have low or extremely low volatility, contribute significantly to SOA formation by condensation, or even form new particles (Bianchi et al., 2019; Ehn et al., 2014; Kirkby et al., 2016, Tröstl et al., 2016). Previous studies have confirmed that low-volatility products from isoprene-NO reactions are the major precursors to SOA (Ng et al., 2008; Rollins et al., 2009; Schwantes et al., 2019). Here the low-volatility compounds refer to gas-phase products that allow fractions to exist in the particle phase and may include the groups of organic compounds with intermediate volatility (IVOC, 300 3 10 g m), semi-volatility (SVOC, 0.3 300 g m), low volatility (LVOC, 3 10 0.3 g m), and extremely low volatility (ELVOC, 3 10 g m) as proposed by Donahue et al. (2012). In general, SVOC, LVOC, and ELVOC can contribute to the SOA formation (Jimenez et al., 2009). In order to evaluate the potential of oxygenated products to form SOA, information about their vapor pressures is essential. However, due to the high degree of functionalization, low or extremely low volatility, and uncertainties in quantification and molecular structures, it is challenging to determine the exact vapor pressure of highly oxidized products. Detailed information on the volatilities of different-generation products is lacking, which impedes the assessment of their contribution to SOA formation.
In this work, we present the results of chamber experiments on isoprene oxidation by NO under near-atmospheric conditions, where NO was produced in situ by O reaction with NO. Subsequent characteristics of multi-generation chemistry of isoprene with NO are investigated. By examining the time evolution of various gas-phase products, we propose possible reaction mechanisms that help to get the possible functionalization of the products. Saturation vapor pressures of the major gas-phase products observed by HR-ToF-CIMS are predicted by using different parameterization methods that are widely used or state of the art in the literature. In addition, we estimate the vapor pressure derived from the equilibrium partitioning coefficient according to the condensation behavior of different products in experiments with and without seed aerosols. Based on these results, the volatility of the major oxidation products stemming from isoprene-NO reaction and their potential to form SOA are evaluated.
2 Experimental and methods
2.1 Atmospheric simulation chamber SAPHIR
All the data presented here were measured in the atmospheric simulation chamber SAPHIR (Simulation of Atmospheric PHotochemical In a large Reaction Chamber) at Forschungszentrum Jülich, Germany, which is designed to investigate the oxidation processes of both biogenic and anthropogenic trace gases and formation of secondary particles and pollutants under near-atmospheric conditions. The SAPHIR chamber is a double-walled Teflon (FEP) cylinder with a volume of 270 m (5 m in diameter and 18 m in length). The large volume-to-surface ratio (1 m) allows experiments to be conducted under natural conditions and reduces interference from the chamber walls. The chamber is equipped with a shutter system which can be opened to admit sunlight for photochemical experiments or closed to mimic nighttime conditions. There are two fans inside the chamber to ensure good mixing of trace gases (within 2 min). The chamber is filled with synthetic air made from mixing ultrapure nitrogen and oxygen (Linde, purity 99.99990 %) and is slightly over-pressured ( 35 Pa) to prevent intrusion of outside air into the chamber. Due to a small amount of leakage ( 7 m h) and gas consumption by instrument sampling, a replenishment flow is provided by a flow control, which leads to a dilution rate of 4 %–7 % h. A more detailed description of the chamber setup and its characterization can be found elsewhere (Rohrer et al., 2005).
2.2 Experiment description
A series of experiments investigating the oxidation of isoprene by NO were conducted in the SAPHIR chamber in August 2018 (ISOPNO campaign) under different chemical conditions. In this work, we primarily focus on an experiment conducted on 8 August 2018 that examined the fast oxidation of isoprene by NO (up to 130 pptv) without seed aerosols. The experiment was performed under dry (RH 5 %) and dark conditions and employed injections of O and NO as sources of NO, where O was generated by a silent discharge ozonizer (O3onia), and high-purity NO was introduced from a gas bottle (Linde, purity 99 %).
Before the experiment, the chamber was flushed overnight with a total amount of 1800 m synthetic air to minimize any remaining contamination. At the beginning of the experiment, the chamber air was slightly humidified (RH 0.1 %) by flushing water vapor from boiling Milli-Q® water into the chamber. Thereafter, O and NO were added to the chamber in succession, and their concentrations in the chamber after injection were approximately 100 and 25 ppbv, respectively, as shown in Fig. 1. After that, 10 ppbv of isoprene was injected using a GC syringe, initiating the reaction with NO. The period between the first and second injection is defined as “step I”, and so on for the other three periods. The second injection was done when isoprene from the first injection was almost completely consumed, to reach concentrations of O, NO, and isoprene in the chamber of 100, 30, and 10 ppbv, respectively. About 1.5 h later, the chemistry was further accelerated by a third injection of precursors, and accordingly the concentrations of O, NO, and isoprene reached 100, 25, and 10 ppbv, respectively. Two hours later, the fourth addition was made and the concentrations of O and NO increased to approximately 115 and 30 ppbv, respectively, aiming to promote further oxidation of early-generation products. In total the system was kept running for about 7.5 h. Calculation from measurements of isoprene, O, OH, NO, and dilution indicates that NO contributed for more than 90 % of the chemical losses of isoprene, as shown in Fig. S1, with reaction with O being a minor pathway in our system. The reaction of isoprene with OH was not considered as OH concentration was below the detection limit of the instrument in this study (Fig. S2). Thus, losses due to reaction with OH could not be quantified from the measurement but have been determined to contribute about 10 % of the isoprene losses according to a recently published modeling work based on the same campaign, with the contribution of the NO reaction accounting for up to 80 % accordingly (Vereecken et al., 2021).
Figure 1
Measurements of (a) O and NO, (b) NO and NO, (c) isoprene, and (d) temperature and relative humidity in the chamber during the experiment on 8 August 2018.
[Figure omitted. See PDF]
A complementary experiment was conducted on 14 August 2018 under similar conditions but with seed aerosols. Approximately 60 g m of ammonium sulfate aerosol was added at the beginning of the experiment. Thereafter, approximately 100 and 20 ppbv of O and NO was introduced to the chamber to produce NO, followed by the addition of 10 ppbv of isoprene about 30 min later (see Fig. S3). Another 6 ppbv of NO and 10 ppbv of isoprene were added about 1 h later to accelerate the reaction. At the last injection, only O ( 50 ppbv) and NO ( 7 ppbv) were added, similar to for the experiment without seeds. The experiment lasted for about 8 h. The results were used to investigate the condensation behavior of various gas-phase products from isoprene oxidation, aiming to estimate equilibrium partitioning coefficients and vapor pressures.
2.3 InstrumentationA high-resolution time-of-flight chemical ionization mass spectrometer (HR-ToF-CIMS, Aerodyne Research Inc., hereafter CIMS) was used to continuously measure the gas-phase products from isoprene oxidation by NO. The ToF-MS was operated in “” mode with the mass resolution power between 3000–4000 . In order to reduce the losses of HO radicals and HOM on the tubing, a customized inlet (Albrecht et al., 2019) was directly connected to the chamber. The CIMS was operated in negative ion mode using Br as the reagent ion, which is selective to polar species such as acids, hydroxy, or nitrooxy carbonyls, as well as HO radicals (Albrecht et al., 2019; Ng et al., 2008; Rissanen et al., 2019; Riva et al., 2019).
Bromide ions were generated by passing a mixture of 10 standard cubic centimeters per minute of 0.4 % CFBr in nitrogen and 2 standard liters of nitrogen per minute through a 370 MBq Po source (Type P-2021-5000, NDR Static Control LLC, USA), resulting in 10 ion counts per second (Albrecht et al., 2019). In our system, about 190 ions were identified for each mass spectrum on average, most of which were detected as adducts with Br, while some acidic compounds ( 7 % of the total) like nitric acid (HNO), glycolic acid (CHO), and malonic acid (CHO) were also detected as deprotonated ions. In addition, there were some ions ( 10 % of the total) identified as adducts with NO. The isotope distribution of Br and Br is approximately , and therefore two signals appear at MW 79 and MW 81 with MW being the molecular mass of the molecule that is detected as cluster with Br. In this work, we will use Thomson (Th) as the unit for mass-to-charge (), and the of molecules discussed in following include the mass contribution from Br 79) if there is no other annotation.
In order to have an indicator of the CIMS performance, perfluoropentanoic acid (PFPA, CFHO) was used as an internal standard. For calibration, five isolated peaks were used, including Br ( 79), HOBr ( 97), HNOBr ( 142), CFO ( 263), and CFHOBr ( 343), covering the mass range of dominant products. The averaged accuracies of all five calibrated masses were below 5 ppm over the whole measurement period. However, due to the low signal intensity, the PFPA cluster (CFOH, 527) was not suited for mass calibration, and there were no other suitable masses with sufficient intensity and high accuracy that could be used to calibrate the higher mass range. Therefore, peak fitting in the mass range between 300 and 500 Th might have higher uncertainties. The CIMS was optimized to gain a maximum signal of [HOBr] isotopes, which are weakly bounded clusters. This was achieved by adjusting step by step the electrostatic field in the transfer stage to minimize fragmentation. During the campaign, the settings of CIMS were kept unchanged to keep a similar performance. However, the signal of reagent ion Br decreased by about 65 % (from 100 000 to 34 000 counts s) over the campaign duration of 4 weeks. In order to minimize the effect of drift in performance, we used the normalized (by the sum of the total ion counts, norm. counts from here on) signals for analysis. The sensitivity for total carbon was calculated by determining the slope of wall-loss-corrected total carbon signals detected by CIMS (only the identified peaks were considered) versus isoprene consumed. As illustrated in Fig. S4a, the CIMS sensitivities were roughly identical in two experiments (0.026 0.002 norm. counts s ppbv on 8 August, and 0.022 0.001 norm. counts s ppbv on 14 August), indicating that different experimental conditions over two days had an insignificant impact on CIMS sensitivity for total carbon and thus the data from these days are comparable. In addition, an inter-comparison of measurements by Br CIMS and I CIMS were made. As shown in Fig. S4b, the measurements of CHNO from the two instruments are well linearly correlated with each other at the early oxidation stages. However, the correlation coefficients of measurements from two instruments deviated from experiment to experiment. This is probably related to different experimental conditions, which might lead to different chemical processes and thus formation of isomers. Since CIMS with different reagent ions might have different sensitivities to isomers and may be selective for different compounds, the correlation coefficients of measurements from Br, and I CIMS may differ from day to day. Moreover, the Br CIMS was not tuned during the campaign while the I CIMS was optimized from time to time. In general, the performance of Br CIMS was stable and the data taken by it are reliable.
The mass spectra data were processed using the software “Tofware” embedded in Igor as provided by Aerodyne Research Inc. (
Isoprene was measured by a Vocus proton transfer reaction time-of-flight mass spectrometer (Aerodyne Research Inc., hereafter Vocus), which has a higher mass-resolving power (nominal 10 000 ) and less inlet wall losses and sampling delays compared to traditional PTR-MS (Krechmer et al., 2018). The mixing ratio of O was monitored by an UV absorption instrument, and that of NO was monitored by a chemiluminescence instrument and a custom-built cavity ring-down spectrometer (CRDS). The concentrations of NO and NO were detected by two custom-built CRDS instruments (Dubé et al., 2006; Sobanski et al., 2016). In addition, temperature and pressure inside the chamber were monitored by an ultra-sonic anemometer and a pressure sensor, respectively. The relative humidity was primarily detected as water mixing ratio by a Picarro CRDS instrument (Crosson, 2008).
The particle number concentrations and their size distributions were measured by a condensation particle counter (TSI 3783, hereafter CPC) and a scan mobility particle sizer (TSI 3081 electrostatic classifier combined with TSI 3025 CPC, hereafter SMPS). The aerosol chemical composition was identified by a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS, Aerodyne Research Inc., hereafter AMS). The ionization efficiency of AMS was determined by using the monodisperse aerosol generated from NHNO and (NH)SO solutions. The collection efficiency (CE) could be estimated based on the particle mass concentration yielded from AMS and that derived from SMPS. In this study, the average CE value of 0.5 is used for correction.
2.4 Methods to estimate saturation vapor pressure
The pure liquid saturation vapor pressure is a thermodynamic metric relevant for the partitioning equilibrium of organic molecules, which determines their propensity to form SOA (Compernolle et al., 2011; O'Meara et al., 2014; Pankow and Asher, 2008). Due to their complex functionalities and low or extremely low volatility, it is challenging to determine the vapor pressures of highly oxidized molecules. As a result, theoretical and semiempirical methods are usually used for vapor pressure estimation. Commonly used semiempirical methods include composition–activity (CA), group-contribution (GC), and structure–activity (SA) methods. The CA methods are the easiest to use, as they only require information on molecular composition for estimation. They are widely applied in the context of the two-dimensional volatility basis set (2D-VBS) (Donahue et al., 2011). For GC methods, the exact functional groups are required to calculate the saturation vapor pressure. The SIMPOL.1 (Pankow and Asher, 2008), the parameterization as described by Nannoolal et al. (2008), and EVAPORATION (Compernolle et al., 2011) are three widely used GC methods. Structure–activity methods can provide more accurate estimates with sophisticated treatments of intramolecular interactions like intramolecular hydrogen bonding (Bilde et al., 2015). However, detailed molecular properties such as boiling point and evaporation enthalpy are required for estimation, which are generally obtained by complex and time-consuming quantum chemical calculations. Therefore, SA methods are not applied for vapor pressure estimation in this study.
Saturation concentration (, mass based) is related to saturation vapor pressure and can be calculated following Eq. (1) (Donahue et al., 2006). The log() is a metric used in the 2D-VBS method to evaluate the volatility of organic molecules.
1 where (8.206 10 m atm K mol) is the gas constant, (K) is the temperature, (g mol) is the molecular weight of compound , is the activity coefficient of compound and here is assumed to be 1 (Donahue et al., 2006), and (atm) is the pure liquid saturation vapor pressure at temperature (298 K).
In this study, different CA methods are applied to calculate the saturation vapor pressures of various oxidation products from isoprene reaction with NO. These include parameterizations that were constrained by chamber measurements as proposed by Donahue et al. (2011), Mohr et al. (2019), and Peräkylä et al. (2020). All of these three parameterization methods have included the effect of the presence of nitrate groups on vapor-pressure estimation. Further we test the GC methods proposed by Nannoolal et al. (2008), Pankow and Asher (2008, SIMPOL.1), and Compernolle et al. (2011, EVAPORATION). All the methods used in this study are summarized in Table 1. The calculations of EVAPORATION and the Nannool method were done via the online molecular and multiphase property prediction facility UManSysProp (
Summary of estimation methods of saturation vapor pressure used in this study.
Estimation method | Methodology | Input information | Reference | ||
---|---|---|---|---|---|
Molecular | Functional | Others | |||
formula | groups | ||||
Donahue et al. | CA | Donahue et al. (2011) | |||
Mohr et al. | ICA | Mohr et al. (2019) | |||
Peräkylä et al. | ICA | Peräkylä et al. (2020) | |||
Nannoolal et al. | GC | Nannoolal et al. (2008) | |||
SIMPOL.1 | GC | Pankow and Asher (2008) | |||
EVAPORATION | GC | Compernolle et al. (2011) | |||
This study | EXP |
Abbreviation of composition-activity method. Abbreviation of improved composition-activity method, which modified the parameterization based on chamber measurements to better fit HOMs. Abbreviation of group-contribution method. Boiling point parameterization method is also required to be defined. Abbreviation of experimental method.
In addition, we take advantage of the measurements in this study to calculate the gas-particle equilibrium partitioning coefficient () by comparing experiments with and without seed aerosols. The partitioning coefficient can be converted to saturation concentration by Eq. (2). 2 where and are the gas- and particle-phase concentrations (g m) of species , respectively, and is the organic aerosol concentration (g m). In this study, is signal of species from CIMS in the experiment with seeds, and is the difference of signals between experiments without and with seeds (under the same isoprene consumption condition). The in the experiment with seeds is in a range of 1–4 g m.
2.5 Pathways to the multifunctional oxidation products2.5.1 Basic peroxy and alkoxy radical chemistry
As mentioned before, information about molecular structures (at least functional groups) is required to calculate vapor pressures by using GC methods. Although the high-resolution ToF-CIMS allows for determining the chemical composition of the detected ions, it is unable to provide information about molecular structures, so that the constitutional or configurational isomers with the same mass cannot be distinguished without additional information. Fortunately, knowledge of detailed chemical formation mechanisms can help to infer the molecular structure information. However, the development of a comprehensive, multi-generational kinetic mechanism for NO-initiated oxidation of isoprene is outside the scope of the current paper. Instead, in order to link the observed mass peaks to representative molecular structures, we developed a framework tracing the chemical oxidation mechanisms by taking well-known oxidation steps to predict the most likely isomeric forms of the functionalized products formed in the isoprene oxidation. For this purpose, we rely on the extensive literature on isoprene, alkylperoxy radical, and alkoxy radical chemistry (Atkinson, 2007; Atkinson and Arey, 2003; Bianchi et al., 2019; Crounse et al., 2013; Ehn et al., 2014; Jenkin et al., 2015, 2019; Kwan et al., 2012; Mentel et al., 2015; Ng et al., 2008; Novelli et al., 2021; Orlando et al., 2003; Orlando and Tyndall, 2012; Rollins et al., 2009; Schwantes et al., 2015; Vereecken and Francisco, 2012; Vereecken and Peeters, 2009, 2010; Vereecken et al., 2021; Wennberg et al., 2018; Ziemann and Atkinson, 2012). This framework is depicted in the supporting information and will be discussed in more detail in Sect. 2.5.2 and 2.5.3. They are based on the following main reactivity trends.
Generally, RO radicals can react with other RO and HO radicals. There are three major channels for the reaction between two RO radicals, leading to alkoxy radicals (RO) (Reaction R1a), as well as termination products like alcohols, aldehydes, or ketones (Reaction R1b) and accretion products (Reaction R1c). These reactions should take place with the first-generation peroxy radicals, as well as with the higher-generation RO radicals formed in the later oxidation steps. Hydroperoxides can be formed from the reaction of RO with HO radicals (Reaction R2a). This reaction can also yield alkoxy radicals (Reaction R2b).
In the presence of NO, RO radicals can also react with NO and NO, leading to the formation of alkoxy radicals (Reaction R3a), organic nitrates (Reaction R3b), and peroxynitrates (Reaction R4) (including peroxyacyl nitrates, PANs, if R RC(O)–). The channel that results in RO radicals is the major pathway for the reaction of RO radicals with NO (Ziemann and Atkinson, 2012). However, reactions of RO radicals with NO (Reactions R3a and R3b) can be neglected in this study due to the high O concentration, which results in rapid conversion of NO to NO. The peroxynitrates formed from the reaction of RO with NO will undergo rapid thermal decomposition under our experimental conditions, with the exception of PANs. The reaction of RO with NO radicals mainly forms NO and alkoxy radicals (Reaction R5), which will continue the radical chains (Reaction R7).
In addition to bimolecular reactions, intramolecular rearrangement (H migration) is a competitive reaction pathway for RO radicals. RO radicals can undergo H migration to form a hydroperoxy functionality (–OOH) and a radical site that can subsequently recombine with an O molecule, leading to the formation of a new, more oxidized substituted RO (Reaction R6). This process is the so-called “autoxidation” path and has been confirmed as significantly important for SOA formation (Crounse et al., 2013; Ehn et al., 2014; Mentel et al., 2015; Praske et al., 2018; Rissanen et al., 2014). The rates of RO H migration are strongly dependent on the structure of RO radicals, and the most likely routes can be derived based on the structure–activity relationship proposed by Vereecken and Nozière (2020).
R6
The RO radicals formed in the reaction of RO RO typically have three accessible pathways, including isomerization by H migration (Reaction R7a), fragmentation (Reaction R7b), and (less important here) reaction with O (Reaction R7c). Like H migration in RO, rearrangement by H shift in RO radicals leads to the formation of more oxidized RO radicals. Fragmentation leads to smaller carbon chains, and this becomes more important for alkoxy radicals with a higher number of (oxygen-bearing) substituents (Vereecken and Peeters, 2009, 2010).
In addition to the above general reaction pathways, we include a number of other reactions in the framework, such as fragmentation of peroxy radicals, epoxidation of -OOH alkyl radicals, and unimolecular termination of nitrooxy or hydroperoxyl peroxy radicals. Details can be found in the supporting information.
2.5.2 Formation of first-generation productsHere “first-generation products” refers to the closed-shell compounds from the first attack of NO at the isoprene double bonds, while “second-generation products” follow an addition of NO to the remaining double bond (or any other oxidation reaction) of a first-generation product. Addition of a NO radical to one of isoprene double bonds and subsequent addition of O to the resulting (delocalized) radical sites leads to the formation of nitrooxy alkylperoxy radicals (INO, CHNO). Since isoprene contains two double bonds, NO can attack any of the four positions on the conjugated carbon bonds, resulting in eight possible INO isomers (including six constitutional and two conformational isomers), as shown in Scheme S1. However, both theoretical and experimental studies suggest that the addition occurs preferably at the primary and terminal carbons, wherein C1 addition seems to be preferred over C addition (Schwantes et al., 2015; Suh et al., 2001; Wennberg et al., 2018). As the GC methods have limited or no ability to distinguish between positional isomers (Kurten et al., 2016), we take exemplarily the products following the C addition for the vapor pressure analysis in this study.
The initial peroxy radicals (CHNO) can undergo rearrangement by H shift from C–H bonds with subsequent O addition, yielding new –OOH functionalized peroxy radicals (Reaction R6). Repeating this process can lead to the formation of a series of peroxy radicals and termination products with a stepwise increasing number of oxygen atoms by 2, as shown in the conceptual Scheme S2 in the Supplement. This is the RO autoxidation channel, and the molecular formula of peroxy radicals formed via consecutive O additions can be represented as CHNO ( 1, number of autoxidation steps). The autoxidation chain can be terminated when the H shift occurs at a carbon with an –OOH or –ONO group attached, leading to carbonyl formation with OH or NO loss (Anglada et al., 2016; Bianchi et al., 2019; Vereecken, 2008; Vereecken et al., 2004). The closed-shell products formed in these termination steps have the general molecular formula CHNO (OH loss channel) or CHO (NO loss channel).
The CHNO peroxy radicals can also react with HO radicals to form –OOH functionalized termination products with the general molecular formula CHNO (Reaction R2a), or yielding the alkoxy radicals CHNO (Reaction R2b). In addition, the CHNO peroxy radicals can react with other RO radicals (Reactions R1a–R1c). The reaction R1a leads to the formation of alkoxy radicals (CHNO) while Reaction R1b forms closed-shell products either with a carbonyl group (CHNO) or a hydroxyl group (CHNO). Alternatively, dimers can be formed following Reaction R1c, which have then two –ONO groups and at least eight oxygen atoms depending on the formula of RO radicals involved, as shown in Table S1.
The alkoxy radicals from Reactions R1a and R2b can undergo unimolecular rearrangement by H shift with subsequent O addition, similar to the RO radicals, forming new RO radicals with a –OH group (Reaction R7a). As mentioned above, when the H shift occurs at a carbon with an –OOH or –ONO group attached, the resulting intermediates tend to lose an OH group or NO (Bianchi et al., 2019), yielding the closed-shell carbonyl products with general formulas CHNO or CHO respectively, as shown in the conceptual scheme Scheme S3. The newly formed RO radicals from the alkoxy H shift channel can follow the peroxy pathways (Reactions R1–R6) like other RO radicals, leading to a range of compounds like hydroperoxides (Reaction R2a, CHNO), alcohols (Reaction R1b, CHNO), aldehydes (Reaction R1b, CHNO), and accretion products (Reaction R1c, CHNO), as depicted in Scheme S3. Alternatively, they can also yield alkoxy radicals again following Reactions R1a and R2b and continue in that manner. Furthermore, the alkoxy radicals can break apart into two fragments according to Reaction R7b.
In general, the alkoxy reaction pathways diversify the parity of the oxygen number of the products from the reaction of isoprene with NO, and the compounds formed via these reactions generally have one less or one more oxygen atom compared to those formed from straight peroxy reaction pathways. With the help of the mechanistic framework described above, we can infer the functionality of first-generation products. This is exemplified in Scheme S5 and S6 for the major first-generation C products. In addition, the reaction pathways and their corresponding structures of the first-generation C dimers (CHNO) are summarized in Scheme S13.
2.5.3 Formation of second-generation products
Nitrate radicals can oxidize the first-generation products once again at the double bond remaining ( (298 K) 3–11 10 cm molecules s, Wennberg et al., 2018). This leads eventually to “second-generation” products that contain at least two nitrogen atoms. Addition of a NO radical to the remaining double bond of the first-generation products results in the formation of dinitrooxy peroxy radicals. We assume that dinitrooxy peroxy radicals can undergo unimolecular and bimolecular reactions (Reactions R1–R6) in analogy to nitrooxy peroxy radicals, which lead to secondary products containing two or more nitrogen atoms, as summarized in the conceptual scheme Scheme S4.
The reaction of first-generation nitrooxy peroxy radicals with NO can also yield 2N compounds (Reaction R4); however, these 2N compounds ought to be under first-generation products by definition. Such species are not discussed in detail here but will be covered to catch the diversity of the functionalities for the vapor pressure estimation. With the help of this secondary reaction framework, we can propose functional groups for the major second-generation products. Schemes S8–S10 depict the detailed (possible) reaction pathways that lead to the formation of detected C dinitrates, as well as their possible structures. Furthermore, the proposed formation mechanism, and their structures for C trinitrates are shown in Scheme S12, while those for the second-generation C dimers (CHNO and CHNO) are depicted in Scheme S13.
2.5.4 Formation of fragmentation products
In addition to the multigenerational C and C products, fragmentation products can be formed from the reaction of isoprene with NO. As mentioned above, the alkoxy radicals can undergo C–C bond scission, producing a carbonyl compound and an alkyl fragment (Reaction R7b). As shown in Scheme S7, when the secondary nitrooxy alkoxy radicals from the further oxidation of C carbonyl compounds (CHO and CHO here) undergo unimolecular decomposition, C carbonyl products (CHNO and CHNO, respectively) are formed as well as formyl radicals. Since the bond fission can occur at different positions, the generation of more reactive C and C carbonyl compounds are possible. In addition, the C carbonyl compounds are possibly generated through peroxy radical arrangement by 1,4 H shift and subsequent acyl radical bond scission reactions (see Scheme S7). The C dinitrates can be formed following similar chemistry, as depicted in Scheme S11.
2.5.5 Candidate structures for vapor pressure estimation
Among all gas-phase products detected by CIMS, we selected 32 major representative organonitrates formed from isoprene oxidation by NO radicals. Their structures are rationalized by the corresponding molecular formulas and proposed formation mechanisms in the reaction framework. Table S2 summarizes all the exemplified structures used for vapor pressure estimation. The functional groups covered in the selected structures include nitrate, hydroxyl, ketone, aldehyde, carboxylic acid, peroxide, hydroperoxide, hydroperoxy acid, peroxynitrate, peroxyacyl nitrate, and epoxide. The structural information allows calculation of the saturation vapor pressure by GC methods.
3 Results and discussion
3.1 Chemical composition of oxidation products
Figure 2 illustrates the average mass spectra of the whole experiment measured by BrCIMS for isoprene-NO reaction. Chemical sum formulas were attributed to most of the detected ions. The gas-phase products were separated into two major groups according to their chemical composition, including monomers comprising C compounds and dimers containing C compounds. There were also products from decomposition reactions with C, which were merged into monomers. The monomers and dimers were further classified into five subgroups as follows. Monomers consisting of compounds with one nitrogen atom (hereafter 1N-monomers) and two or three N atoms (2N- or 3N-monomers) mainly accumulate in 220–280 Th, 300–340 Th, and 350–390 Th, respectively, while dimers containing compounds with two N atoms (2N-dimers) and three N atoms (3N-dimers) appear in 370–440 Th and 450–520 Th, respectively. As shown in Fig. 2, the signal intensities decrease from 1N-monomers, 2N-monomers, and 2N-dimers to 3N-monomers and 3N-dimers. Many of the compounds detected in this work were also observed in previous isoprene-NO systems (Kwan et al., 2012; Ng et al., 2008; Schwantes et al., 2015). In this work, only closed-shell products are considered for analysis.
Figure 2
Averaged mass spectra for isoprene-NO experiment on 8 August 2018. Molecular formulas were determined according to the accurate mass data provided by HR-ToF-CIMS.
[Figure omitted. See PDF]
The 1N-monomer CHNO at 242 is the dominant product formed from the NO-induced isoprene oxidation in our experiment, followed by the 1N-decomposition product CHNO at 228. In addition to CHNO, several analogues with molecular formulas CHNO and CHNO are in relatively high abundance. CHNO and CHNO are the major 2N- and 3N-monomers. Their signal intensities are 1 to 2 orders of magnitude lower than those of 1N-monomers. According to the chemical composition, the 1N-monomers are likely to be the first-generation products from NO oxidation of isoprene, while the 2N- and 3N-monomers probably arise from the further oxidation of 1N-monomers by NO, which therefore should be second- or later-generation products. As mentioned before, the reaction of nitrooxy alkylperoxy radicals with NO can lead to the formation of peroxynitrates (for the special case peroxyacyl nitrates, PAN-like) containing two N atoms. The peroxynitrates will decompose rapidly under experimental conditions, whereas the PAN-like compounds are more stable (with lifetimes ranging from minutes to weeks at 298 K and ambient temperature). Such C PAN-like compounds are isomers of aforementioned 2N-monomers but ought to be first-generation products. In addition to C-2N-monomers, we observe some C-2N-monomers with relatively high intensity, such as CHNO at 273 and CHNO at 291. It is proposed that such C dinitrates originate from the further oxidation of C carbonyl compounds followed by unimolecular decomposition (Schwantes et al., 2015; Wennberg et al., 2018), as shown in Scheme S11.
2N-Dimers are C compounds with 8–12 oxygen atoms (CHNO), and their signal intensities are relatively low compared to that of monomers, approximately 3 orders of magnitude lower. They might be ROOR products from the self- or cross-reaction of two nitrooxy peroxy radicals (Berndt et al., 2018). 3N-Dimers are molecules consisting of 12–16 oxygen atoms (CHNO). They are probably formed from further oxidation of 2N-dimers or from the cross-reaction of a nitrooxy peroxy radical with a dinitrooxy peroxy radicals.
3.2 Multi-generation chemistry3.2.1 Molecular composition for each step
As mentioned in Sect. 2.2, there were four injections during the experiment on 8 August (denoted as step I, II, III, IV in Fig. 3), wherein in the first three injections all components, O, NO, and isoprene, were added, while in the last step only O and NO were injected to promote the further oxidation of early-generation products. The extended oxidation time with reinjection of oxidants provides the opportunity to investigate the multi-generation oxidation chemistry of the isoprene-NO system. The mass spectra show only slow changes in the concentrations during the last period of each step, indicating weak chemical evolution. Therefore, we use integrated mass spectra over the last 10 min of each step for further analysis. Due to the similarity of the integrated mass spectra for step II and step III, the latter is omitted in Fig. 3.
Figure 3
Comparison of the chemical composition of each oxidation step. (a) Averaged mass spectra for step I, II, and IV, with the omitted spectrum of step III being very similar to that of step II. (b) Relative contribution of different chemical groups for each oxidation step. Only organic products were counted for analysis. “Others” refers to CHO compounds that do not contain nitrogen atoms (e.g., CHO and CHO).
[Figure omitted. See PDF]
As shown in Fig. 3a, large amounts of 1N-monomers were formed from NO oxidation of isoprene in step I, wherein CHNO, CHNO, and CHNO are the most abundant compounds in signal. The 2N-monomers, which are expected from further oxidation of 1N-monomers, are much less abundant compared to 1N-monomers, accounting for 5.0 % of the total organic signals. The 3N-monomers are even less abundant (0.04 %). The low contributions of second-generation products probably results from the relatively high concentration of isoprene in step I, reducing the possibility for further oxidation of first-generation products. These results indicate that the system is dominated by first-generation chemistry at the early stage and therefore the oxidation state of products is low. In addition to monomers, some 2N- and 3N-dimers are observed. They contribute 1.7 % and 0.2 %, respectively, to the total organic signals, as shown in Fig. 3b. The low signal intensity of dimers probably results from their small yield under our experimental conditions. In this case their contribution to SOA formation might be small. However, a part of the dimers condense onto the chamber wall due to their low volatility, so only a smaller portion exists in the gas phase (compare Table S3 and Fig. S5).
In step II, the secondary chemistry was accelerated by the further addition of O and NO, but the primary chemistry was also maintained by isoprene injection. As a result, more 1N-monomers (e.g., CHNO) were formed compared to step I, as well as dimers (e.g., CHNO and CHNO), as shown in Fig. 3a. The signals of 2N-monomers almost double in this period compared to those in step I, and their relative contribution increases from 5.0 % to 7.4 %. This is attributed to the further oxidation of first-generation products formed in step I. The relative contributions of different chemical groups exhibited in Fig. 3b clearly show that, although NO produced from the second addition of NO and O still primarily reacted with newly injected isoprene, reaction of NO with the first-generation oxidation products retaining a double bond was inevitable, leading to more second-generation 2N or 3N products compared to step I. The visibly increasing fraction of 2N-monomers indicates that the second-generation chemistry started to play a more important role than that in the early stage. In step III, the chemical process proceeded similarly and thus is not further discussed here.
Due to the favorable conditions for further oxidation, the signals of 1N-monomers (such as CHNO, CHNO, and CHNO), as well as 2N- and 3N-dimers, dropped dramatically in step IV, with their relative contributions decreasing to 58.1 %, 0.5 %, and 0.15 %, respectively. The decrease in signals of dimers is primarily ascribed to lack of isoprene, as there were fewer peroxy radicals under this condition, and hence fewer dimers were formed. In addition, their condensation on the wall and dilution also contributed to the decreasing signals. Furthermore, dimers with 2 or 3 nitrogen atoms possess at least one double bond in their molecular structures and can thus be further oxidized under high-NO conditions to form 4N- or 5N-dimers. However, only few 4N-dimers and no 5N-dimers were detected by CIMS, suggesting that the 4N- and 5N-dimers were either not formed, or if present, had lower absolute concentrations below the detection limit (approximately 5 10 and 5 10 molecules cm for salicylic acid and acetic acid, for an integration time of 60 s). condensed on the wall due to their low volatilities. In contrast, 2N- and 3N-monomers increase significantly, with their relative contributions ascending to 20.0 % and 0.29 %, respectively. This indicates that 2N- and 3N-monomers might be second- or later-generation products that are formed from the further oxidation of first-generation products. Additionally, unlike the C monomers, the signal of CHNO increased in step IV, indicating that there is a new formation pathway for CHNO under excess-NO conditions. No double bond can remain in such products, as otherwise they would be oxidized and their signal should decay instead.
In summary, the above findings confirm that multi-generation chemistry happened during the NO-initiated isoprene oxidation, and that the later generation oxidation was promoted by “excess” NO radicals.
3.2.2Carbon oxidation state ()
The oxidation state of carbon () is defined as the charge a carbon atom takes with assumption that it loses completely all electrons in bonds to more electronegative atoms and vice versa (Kroll et al., 2011). This quantity is a metric for the degree of oxidation and will increase with oxidation. Moreover, together with carbon number can be used to constrain the composition of organic mixtures and provide insights into their evolutions. The carbon oxidation state of a species is determined by the relative abundances and oxidation states of non-carbon atoms in the compound. Since we observed nitrate groups in the products, is defined by Eq. (3). In this study, the group-averaged is the signal-weighted mean average carbon oxidation state of compounds with the same carbon number, and the bulk-averaged is the signal-weighted mean average carbon oxidation state of all detected compounds in the system. 3 wherein , , and are the number of the respective atoms in the molecular formula.
Figure 4 shows the distribution of gas-phase products from the isoprene-NO system in the oxidation state versus carbon number (OS vs. ) space. The bulk-averaged is 0.35 in step I, wherein the smaller molecules (C) have higher oxidation states than the larger molecules. The group-averaged oxidation state of C compounds is relatively low (), indicating that both the oxidation and autoxidation degrees of isoprene are quite low during this period. This is consistent with the conclusion made previously from mass spectra results that at the early stage isoprene-NO oxidation was dominated by first-generation chemistry.
The system increases to 0.26 in step II, confirming that first-generation products were further oxidized after the second injection. During this step, the of most compound groups increases only weakly, except for that of the C compounds. The group-averaged of C compounds increases to 0.60 in step II, which is the major contributor to the increase in of the whole system. The increase in of C compounds is largely attributed to the formation of 2N-monomers expected from further oxidation of existing 1N-products formed in step I. This is confirmed by the detectable increase in 2N- and 3N-monomers in the mass spectra and their higher relative contributions to total signals (see Fig. 3). In addition to C compounds, the of C and C products increase significantly in step II.
In step IV, the secondary oxidation was largely accelerated by reinjection of O and NO, and hence the system oxidation degree increases, with the bulk-averaged growing substantially to 0.09. Similarly, the significant increase in system is mainly attributed to the C compounds, with their group-averaged increasing to 0.31. In addition, the of C compounds increased evidently despite their decreasing signals, suggesting C dimers were further oxidized as well in step IV. It is worth noting that the average carbon number decreases step by step with increasing . This is the case because of fewer C products, but more fragments were formed with the reaction proceeding, as shown in Fig. 4 by the decreasing peak areas of larger molecules but a converse trend for smaller molecules. One conceivable explanation for the decreasing dimers but increasing fragments with the increasing is that, with more highly oxidized RO formed under high-NO conditions, the prevailing fate of RO changes from dimerization to forming alkoxy radicals, which would undergo unimolecular decomposition rapidly, especially when there is a neighboring oxygen-containing functional group (Molteni et al., 2019).
Figure 4
Distribution of gas-phase products from isoprene oxidation by NO in the carbon oxidation state (OS) versus carbon number () space. Markers are colored by oxygen-to-carbon molar ratio and sized by the logarithm of peak areas. The group-averaged and bulk-averaged are the signal-weighted mean average carbon oxidation state of compounds with the same carbon number and of all detected compounds, respectively.
[Figure omitted. See PDF]
In the oxidation system, the increase in is attributed to the formation of bonds between carbon and oxygen as well as other electronegative atoms, and/or the breaking of bonds between carbon and hydrogen and other electropositive atoms (Kroll et al., 2011). The –ONO group has an oxidation state of 1, which means that addition of a –ONO group to isoprene will increase its by 0.2. According to our estimates, the values of system increased by 1.25 (step I), 0.09 (step II), and 0.35 (step IV), indicating that the increases in are not only due to addition of –ONO group(s) but also to other oxygen-containing functionalities. In addition to functionalization, it is possible that other reactions such as fragmentation and oligomerization which can increase or reduce the oxidation state were involved during the reaction.
As mentioned above, the average carbon oxidation state of a mixture of molecules largely depends on its chemical composition. Therefore, for different oxidation systems, their may differ due to different precursors and oxidation conditions. In our study, the of a NO-initiated isoprene oxidation system increased from 0.35 to 0.09 with further oxidation. For OH- and O-initiated systems, the average oxidation state of laboratory-generated isoprene SOA is reported to range from 1.3 to 0.2, as listed in Table S4. It seems that the SOA generated from chloride-initiated oxidation of isoprene is more oxidized compared to other isoprene oxidation systems, for which the can be as high as 1.8 according to limited studies (Wang and Ruiz, 2017). With regard to ambient measurements, the calculated values of organic aerosol and aerosol fractions fell into a wider range between 2 and 2, depending on the site position and the corresponding oxidation environment of that site (Table S4).
In summary, isoprene and its products undergo further oxidation by NO, leading to an increase in the degree of oxidation of products as the reaction proceeds. The increasing bulk is largely governed by the highly oxidized C compounds. In addition, more fragments but fewer dimers are formed as the increases, which can probably be explained by the change of RO fate from prevailing dimerization to fragmentation through the alkoxy radical channel.
3.2.3 Characteristics of different-generation products1N-monomers
To illustrate the multi-generation chemistry involved in the isoprene-NO reaction system, Fig. 5 shows the time evolution of the major gas-phase products. The signal of the most abundant compounds, CHNO, increases rapidly as soon as the reaction is initiated, reaching a maximum when its chemical production rate matches its loss rate (including chemical destruction, wall loss, dilution, etc.), and decreases slowly thereafter. Its time behavior in the first three steps is similar. In step IV, however, the injection of O and NO resulted in a strong decay of CHNO, owing to the occurrence of further oxidation by NO. The time behavior suggests that the CHNO signal is dominated by first-generation oxidation products, and the same conclusion can be made for CHNO and CHNO. According to the mechanistic framework developed above, the CHNO, CHNO, and CHNO compounds most likely correspond to hydroxyl nitrates, nitrooxy hydroperoxides, and hydroxy hydroperoxy nitrates, respectively, but other constitutional isomers are possible. They were already observed in previous studies and were proposed to form through reactions of INO radicals with RO, HO, and unimolecular rearrangement, as shown in Scheme S5 (Ng et al., 2008; Kwan et al., 2012; Schwantes et al., 2015; Wennberg et al., 2018).
Figure 5
Time evolution of selected gas-phase compounds measured during the isoprene – NO experiment on 8 August 2018. (a) Time series of O, NO, NO, and isoprene. (b)–(f) Time evolution of major 1N-monomers (CHNO and CHNO), 2N-monomers (CHNO, CHNO, and CHNO), 3N-monomers (CHNO), 2N-dimers (CHNO), and 3N-dimers (CHNO).
[Figure omitted. See PDF]
As shown in Fig. 5b, the temporal evolution of CHNO ( 274) is different to CHNO compounds, suggesting that it has a completely different formation pathway. Specifically, the formation rate of CHNO is initially much slower than that of CHNO but accelerates to become comparable to them later as the experiment proceeds, i.e., when a multitude of first-generation products are accumulated. This implies that CHNO is produced from the further oxidation of first-generation products, and its signal is dominated by second-generation products. Based on its molecular composition, CHNO could be the dihydroperoxy nitrate as shown in Scheme S5, but its formation through the reaction of HO with nitrooxy hydroperoxy radical from INO autoxidation suggests it should be first-generation products, not in accordance with the time behavior we actually observe. Consequently, we can conclude that it is not the major formation pathway that contributed to CHNO observed in this study. As shown in Scheme S7, the first-generation C hydroxy carbonyl (CHO, 179) can be further oxidized by NO and the resulting alkyl radical would rapidly recombine with O, producing a new peroxy radical, which then reacts with HO radicals to form CHNO (hydroxy hydroperoxy carbonyl nitrate). Similarly, the C hydroperoxy carbonyl (CHO, 195) can also lead to the formation of such CHNO (isomer of that formed through CHO channel) through further oxidation (see Scheme S7). According to the above two mechanisms, CHNO formed following such reaction pathways should be second-generation products, which is consistent with its time behavior.
Considering its similar time behavior to CHNO, the observed CHNO ( 228) signal is likewise thought to be dominated by second-generation products. Schwantes et al. (2015) proposed such a C product based on OH-initiated chemistry, but as the OH concentration in our system was below the detection limit during the experiment (see Fig. S2), this formation pathway cannot apply in our situation. Instead, we suggest that CHNO is formed through the unimolecular decomposition of the C alkoxy or acyl radicals, which result from further oxidation of the C hydroxy carbonyl (CHO, 179), as shown in Scheme S7. It should be pointed out here that there may be reaction pathways forming CHNO as first-generation products that are not considered here, whereas there is no doubt that the second-generation chemistry played a dominant role in CHNO formation according to its time evolution measured by CIMS.
Although CHNO and CHNO show similar time behaviors in the first three steps, it seems that they followed fairly different reaction pathways when the concentration of NO in the chamber increased dramatically in step IV. As shown in Fig. 5b, the signal of CHNO drops immediately after the injection of O and NO, while that of CHNO continues to increase, although its formation rate becomes slightly lower with increasing NO concentration. The decay of the CHNO signal can be explained by more chemical destruction or less production under high-NO conditions, wherein the latter seems more sensible in terms of its structure (no double bond remaining). As shown in Scheme S7, the second-generation CHNO and CHNO compounds share the same precursor in the CHO channel. Consequently, the production of CHNO through this pathway would be interrupted immediately after the injection of O and NO like CHNO. In reality, its signal might decay even faster due to the larger reaction rate of RO H shift (leading to the formation of CHNO) than that of RO reacting with HO (leading to the formation of CHNO). As presented by Vereecken and Nozière (2020), the rate coefficient of aldehydic H shift is 0.5 s (298 K), while the pseudo first-order rate coefficient of RO reacting with HO is 10 s ( (298 K) 5 10 cm molecules s (Atkinson, 2007), and [HO] 4 10 molecules cm), about 2 orders of magnitude smaller. This result implies that the increasing CHNO observed is contributed to by other formation pathways. As mentioned before, CHNO can also be produced by CHO oxidation. We find that the signal of CHNO ( 244), which results from CHO oxidation as well, continues increasing after the injection of O and NO. This tentatively confirms that the production of CHNO in step IV is mainly from the CHO oxidation channel. More experimental or theoretical studies are needed to provide insights into these differences.
2N- and 3N-monomers
As shown in Fig. 5c, 2N-monomers formed much slower than 1N-monomers in the early stage, but their formation rates were accelerated in step II and step III, probably due to the accumulation of first-generation products. According to our mechanistic framework, 2N-monomers are second-generation products resulting from the further oxidation of 1N-monomers by NO, which is consistent with their time behaviors detected by CIMS.
Like CHNO and CHNO, different 2N-monomers have similar behavior in the first three steps, but they are obviously different in step IV when the concentration of NO increased drastically in the chamber. For instance, the signals of CHNO, CHNO, and CHNO continue to increase after the injection of O and NO, while that of CHNO drops immediately. This is related to their detailed formation mechanisms, which are outside the scope of this study. Furthermore, CHNO and CHNO decay a little bit faster than CHNO and CHNO, which might be related to their volatility and will be further discussed in next section.
Different from other 2N-monomers, the signals of CHNO ( 301) increase continuously under high-NO conditions, although its net formation rate is almost zero at the end of step IV. The characteristics of CHNO under high-NO conditions reflects its different formation pathways from other dinitrates, and without having a comprehensive knowledge of its chemical mechanism, we are unable to tell what exactly leads to the differences. In the Master Chemical Mechanism (MCM v3.3.1), CHNO is proposed to be a PAN-like compound stemming from the C nitrooxy carbonyl (CHNO) (
In addition to C-2N-monomers, we observe some C dinitrates such as CHNO ( 273) and CHNO ( 291), and the signal intensity of CHNO is comparable to the major C-2N-monomers. C dinitrates have rarely been mentioned in previous isoprene-NO studies. As shown in Fig. 5c, CHNO has similar time behavior to C-2N-monomers and hence is thought to be a second-generation product. Wennberg et al. (2018) proposed that such a C dinitrate was generated from OH-initiated further oxidation of CHNO. However, this is not applicable here due to a lack of OH radicals in our system. Instead, we propose that the CHNO observed in this study is a dinitrooxy carbonyl compound resulting from NO oxidation of CHNO with subsequent unimolecular decomposition (see Scheme S11 for details).
As shown in Fig. 5d, 3N-monomers are generated more slowly than 1N-monomers, but their signals grow gradually as the experiment proceeds, with a significant increase especially for CHNO in the last step. Furthermore, we can see from Fig. 5c and d that the signals of C trinitrates in step IV appear anticorrelated to that of CHNO and CHNO. The gas-phase 3N-monomers have rarely been reported in the previous literature. Ng et al. (2008) observed a CHNO compound in the particle phase and assumed that it was produced from NO oxidation of the C hydroxy nitrate (CHNO). Similarly, CHNO and CHNO can be formed through NO reacting with dinitrooxy peroxy radicals, which result from corresponding first-generation nitrooxy compound (C hydroperoxy nitrate, CHNO, or C hydroxy hydroperoxy nitrate, CHNO) oxidation by NO radicals, as shown in Scheme S12. 3N-Monomers formed following such pathways are second-generation products by definition. Regarding the rising signals of 3N-monomers in step IV, one explanation is that although the reaction of dinitrooxy peroxy radicals with NO is not an oxidation process, their formation can be significantly facilitated by increasing NO concentration. It is also possible that 3N-monomers are formed through H abstraction of 2N-monomers. NO radicals can abstract the hydrogen of dihydroxy dinitrate (CHNO) or hydroxyl hydroperoxy dinitrate (CHNO) from the carbon with an –OH, –OOH or –ONO group attached, leading to alkyl radicals that can subsequently recombine with O and then react with NO or NO, yielding trinitrates or peroxynitrates containing three nitrogen atoms. 3N-Monomers stemming from such reactions ought to be third-generation products. However, we should point out that 3N-monomers formed following H-abstraction pathways are less likely because abstracting hydrogen from the hydroxyl, hydroperoxy, or nitrooxy carbon would lead to fragmentation in most cases (Bianchi et al., 2019).
In addition, it is interesting to note that the signal of CHNO increases continuously throughout step IV, whereas that of CHNOand CHNO drop after a short period of growth. Meanwhile, the production of CHNO is facilitated by the increasing NO concentration compared to that of CHNO and CHNO. Currently, we cannot explain what exactly causes these differences, but we suspect that there may be different chemical pathways forming different 3N-monomers that are not covered here and may also be related to their different physical properties, such as vapor pressures.
2N- and 3N-dimers
As shown in Fig. 5e, 2N-dimers (except for CHNO) display very similar time behavior to 1N-monomers, which form rapidly after each injection, indicating that the signals of 2N-dimers are dominated by first-generation products like most 1N-monomers. It is noted that the time behavior of CHNO ( 419) is completely different from that of other 2N-dimers. As illustrated in Fig. 5e, the production rate of CHNO is initially much slower compared to other dimers. Besides, its signal increases monotonically in the first two oxidation stages, whereas that of the others always increases first, approaching the maximum as its chemical production competes against the losses, and decreases gradually thereafter. The special time behavior of CHNO suggests that it has a different formation pathway from other 2N-dimers, and its signal is most likely dominated by secondary products. In addition, we find that the signal of CHNO always starts to decay earlier than that of CHNO and CHNO. If we assume that their production rates have the same order of magnitude (confirming by their formation rates after each injection), then it can be concluded that CHNO had additional chemical destruction, or its volatility is much lower than CHNO and CHNO and hence has more rapid loss on the wall. It seems the second hypothesis is more likely when comparing its signal with and without dilution and wall-loss corrections (see Fig. S5). More detailed discussion about volatilities of different isoprene organonitrates will be provided in the next section.
It is proposed that dimers (ROOR) are likely formed through the self- or cross-reaction of two peroxy radicals (Berndt et al., 2018). Consequently, the generation number of dimers depends only on how the involved peroxy radicals are formed. Table S1 summarizes the possible permutation scheme of 2N-dimers from RO RO reactions, and their structural information can be found in Scheme S13. For example, self-reaction of two C nitrooxy peroxy radicals (CHNO) leads to the formation of a CHNO compound, while recombination of two C nitrooxy hydroxyl peroxy radicals (CHNO) or a C nitrooxy peroxy radical (CHNO) with a C nitrooxy hydroperoxy peroxy radical (CHNO) results in CHNO compound. According to their time behavior, 2N-dimers (except for CHNO) are thought to be first-generation products, and from this fact we can infer that the peroxy radicals contributing to dimer formation are dominated by first-generation intermediates. With regard to CHNO, we conclude that it is most likely a secondary product considering its typical second-generation behavior. In other words, at least one of the two C nitrooxy peroxy radicals involved in formation of CHNO must be a secondary intermediate. As listed in Table S1, CHNO can be formed through CHNO CHNO or CHNO CHNO reactions, wherein CHNO and CHNO would be secondary peroxy radicals if they are formed through further NO oxidation of the C hydroxy carbonyl compounds (CHO or CHO), as shown in Scheme S7. In addition, it is possible that CHNO is formed from a C hydroxy peroxy radical CHO reacting with a C dinitrooxy hydroxy carbonyl peroxy radical CHNO (from CHNO oxidation by NO), as we observe high abundant CHO during the experiment, although CHO is assumed to be the major product of the OH-initiated chemistry.
Apart from 2N-dimers, we observe detectable signals at 450, 466, 482, 498, and 514, which are identified as 3N-dimers with molecular formulas CHNO. CHNO and CHNO were detected in the particle phase in a previous study, suggesting that they have low volatility and can contribute to SOA formation (Ng et al., 2008). As shown in Fig. 5f, 3N-dimers form much slower than 2N-dimers, but their productions are accelerated as the experiment proceeds. This is similar to the characteristics of second-generation 2N- and 3N-monomers to some degree, suggesting that the signals of 3N-dimers we observed are most likely dominated by secondary or even later-generation compounds.
It is worth noting that CHNO and CHNO have two completely different types of time behavior. The signals of CHNO, CHNO, and CHNO more or less increase in the first three oxidation steps and start to decline late in step III with increasing NO concentration. As depicted in Scheme S13, 3N-dimers can result from further oxidation of 2N-dimers or the cross-reaction of a first-generation nitrooxy peroxy radical with a secondary dinitrooxy peroxy radical. Accordingly, such 3N-dimers are thought to be second-generation products, and they would further react with NO due to the remaining double bond in their molecular structure, leading to severe chemical destruction of these compounds under high-NO conditions. This is consistent with the time behavior of CHNO, CHNO, and CHNO. In contrast, CHNO and CHNO are formed even more slowly, and their production in the first 4 h is close to zero. However, their signals start to climb late in step III, during which those of CHNO, CHNO, and CHNO decline. This suggests that CHNO and CHNO formed under high-NO conditions probably result from further reactions of CHNO. However, this assumption is highly uncertain, and more experimental and theoretical studies are needed to substantiate it. In terms of their time behavior, CHNO and CHNO are thought to be third- or even later-generation products.
3.3 Volatility distribution of isoprene nitrates3.3.1
estimated by experimental methods
Detailed information about the volatility of organic molecules is essential to evaluate their potential to form SOA. In order to investigate the potential contribution of various isoprene oxidation products to SOA formation, we use our (limited) experimental data to estimate the vapor pressure of different isoprene organonitrates on the basis of their condensation behavior. Figure 6 shows how the signals of gas-phase products change in experiments with and without seed aerosols (ammonium sulfate). Please note that while the two experiments were conducted under similar conditions, the procedures could not be kept fully identical as aerosol seeding required specific measures and the oxidation chemistry might be slightly altered (e.g., due to initiation of heterogeneous reactions).
Figure 6
Time evolution of selected major gas-phase products during experiments with (red) and without (blue) seed aerosols (ammonium sulfate). Signals have been corrected for dilution.
[Figure omitted. See PDF]
As shown in Fig. 6, the signals of most of the selected compounds decline when there are seed aerosols in the chamber, indicating that part of the condensable vapors is partitioned to the particle phase due to the introduction of condensation sinks. The decrease in signal differs for different products, mostly depending on their vapor pressures. As expected, the lower volatility of a compound, the higher the fraction that condenses. For instance, the signal of CHNO decreases by more than 70 % in an experiment with seed aerosols, compared to less than 40 % on average for other less-oxidized 1N-monomers. In some cases (e.g., CHNO and CHNO), however, the product signals in an experiment with seed aerosols are higher than that without seeds after the consumed isoprene exceeding a certain level. In addition, the signal of CHNO in the experiment with seeds is always higher compared to that without seeds. One explanation for this phenomenon is the effect of heterogeneous reactions. It is likely that some condensed compound (denoted as A) can react on the particle surface to form new products with the molecular composition of compound B, or alternatively forming a precursor of B. When they evaporate back to the gas phase, it can result in an increase in the signal of compound B. That is why a higher signal was observed for such compounds in an experiment with seeds than that without seeds, as observed for CHNO in this case.
Based on the observed condensation behavior of different products, we can derive their vapor pressures from the gas-particle equilibrium partitioning coefficients by Eq. (2). As depicted in Fig. 7, the saturation concentrations of different organonitrates show a decreasing tendency from 1N-, 2N-monomer and 3N-monomers to 2N- and 3N-dimers, suggesting that dimers have a higher propensity of condensation and contribute to SOA formation. This is partly related to their molecular weight, as larger molecules generally have lower vapor pressures. However, it cannot explain all the features of the volatility distribution. For example, CHNO (corresponding to no. 8 in Fig. 7) has higher mass than CHNO (corresponding to no. 7 in Fig. 7) but is predicted to have higher vapor pressure. In general, chemical composition and functionalities have significant effects on vapor pressure. For instance, the 2D-VBS composition-activity relationship suggests that each carbon and oxygen decreases by 0.475 and 1.75 decades, respectively (Donahue et al., 2011). Different functional groups also have a very different effect on volatility. For example, each hydroxyl group (–OH) or hydroperoxy group (–OOH) typically reduces the volatility by 2.4 to 2.5 decades, while the less polar carbonyl group ( O) reduces the volatility by 1 decade (Pankow and Asher, 2008; Donahue et al., 2011). The nitrooxy group (–ONO) has a similar reductive effect on vapor pressure, which typically reduces by 2.5 orders of magnitude (Pankow and Asher, 2008). Here, the irregularly high vapor pressure of CHNO is most likely attributed to the functional groups it contains. As listed in Table S2, CHNO is proposed to be a nitrooxy hydroxy hydroperoxyl compound, which consists of two highly polar functional groups –OH and –OOH, contributing to formation of intramolecular H bonding that can significantly increase the vapor pressure (Bilde et al., 2015; Kurten et al., 2016), while CHNO only contains a –OOH group and hence cannot form intramolecular H bonding. This explanation is also valid for CHNO and CHNO. In summary, these findings underline that the constitutional and configurational information of a molecule is critical for vapor-pressure estimation.
Figure 7
Saturation concentrations (in g m, at 298.15 K) of isoprene organonitrates estimated by using experimental and parameterization methods. The numbers correspond with the compound numbers of given in Table S2 (nos. 1–9, 10–18, 19–21, 22–27, and 28–32 corresponding to 1N-monomers, 2N-monomers, 3N-monomers, 2N-dimers, and 3N-dimers, respectively). Marker shapes indicate different isomers, with their size scaled by carbon oxidation state (OS).
[Figure omitted. See PDF]
3.3.2estimated by different parametrization methods
For comparison, we also adopt different parameterization methods to estimate the saturation vapor pressures of isoprene oxidation products based on their molecular composition and the proposed structures, with the results depicted in Fig. 7. In general, the saturation concentrations calculated by different parameterization methods show a similar volatility distribution to that calculated by the experimental method, with of 1N-, 2N-, and 3N-monomers and 2N- and 3N-dimers decreasing in turn. However, different parameterization methods lead to the predicted vapor concentrations with a variability of several orders of magnitude for the same compound, and the discrepancies become larger and larger with more complicated molecules. In addition, of structural isomers calculated by the same method could span several decades.
As shown in Fig. 7, the Donahue et al. parameterization mostly provides lower compared to the three GC methods, with a maximum discrepancy up to 12 orders of magnitude for dimers. With regard to smaller and less oxidized 1N-monomers, predicted values from different methods are in relatively good agreement with each other, whereas the disagreement increases to 11 orders of magnitude for 2N- and 3N-monomers. This is mainly the case because the organic molecules were regarded as a mixture of carbonyl ( O) and –OH functional groups in the Donahue et al. parameterization, and their relative abundance was assumed to be (Donahue et al., 2011). In consequence, the –OOH functional group in peroxides is treated as two –OH groups when adapting the method proposed by Donahue et al. (2011). However, it is demonstrated that the extra oxygen in peroxy moieties has little contribution to reduce vapor pressure (Pankow and Asher et al., 2008); hence treating –OOH equivalently to two –OH functional groups would underestimate the vapor pressures of hydroperoxyl compounds. Furthermore, organic compounds consisting of multiple polar functional groups (such as hydroperoxy, peroxy acid, and peroxide functional groups) tend to form intramolecular H bonding, which would increase the vapor pressure (Bilde et al., 2015; Kurten et al., 2016). All these issues contribute to an underestimation of the vapor pressures of multifunctional products when using the Donahue et al. parameterization. Mohr et al. (2019) improved the parameterization for vapor-pressure estimation by taking the presence of –OOH functional groups in HOM explicitly into consideration and revising the parameters to reduce the effect of –OOH on depressing . Consequently, the Mohr et al. parameterization effectively reduces the discrepancy between its estimates and those predicted by the GC methods, with the differences within 6 orders of magnitude. Nevertheless, there is a slight tendency to underestimate the vapor pressures of 3N-monomers and both 2N- and 3N-dimers. The Peräkylä et al. parameterization method, which was derived from measurements of the condensation behavior of HOM produced from -pinene ozonolysis, predicts similar to the Donahue et al. method for 1N-monomers, but higher for 2N- and 3N-monomers like the Mohr et al. method. As for dimers, especially for the 3N-dimers containing more multifunctional groups, the Peräkylä et al. method even predicts higher than the GC methods in most cases.
Three GC methods predict similar saturation vapor pressures for different isoprene nitrates in this work, with the differences within 5 orders of magnitudes. Generally, the SIMPOL.1 method always provides higher compared to another two methods, and the disagreement between methods becomes larger for molecules containing multifunctional groups. For instance, the vapor-pressure discrepancy between SIMPOL.1 and another two GC methods are both 2 orders of magnitude for CHNO and CHNO, but it increased up to 4 and 5 orders of magnitude, respectively, for CHNO and CHNO.
It is worth noting that the Nannoolal et al. method is able to distinguish between positional isomers (e.g., the estimated for two CHNO isomers are 0.858 and 0.333 g m, respectively), whereas such capacity of EVAPORATION method is limited (e.g., it is able to distinguish between the position isomers of CHNO, but it predicts identical for CHNO isomers). In this respect, the SIMPOL.1 method cannot distinguish between positional isomers at all. Moreover, the SIMPOL.1 method predicts smaller differences between functional group isomers for 1N-monomers and 3N-dimers compared to the Nannoolal et al. method and EVAPORATION, but there is no such regular pattern for 2N-monomers and 2N-dimers.
By comparing the results calculated by the experimental method with those by different parameterization methods, we can see that the GC methods predict higher-saturation concentrations for 1N-monomers than the experimental method, while the Donahue et al. and Peräkylä et al. method provide similar values. With regard to 2N-monomers, the GC methods predict higher vapor pressures compared to the experimental method, but the discrepancy decreases with decreasing saturation concentration. The disagreement of for 2N-monomers estimated by the experimental method and the Mohr et al. or Peräkylä et al. method are within 2 orders of magnitude. In terms of low-volatility dimers, however, the vapor pressures calculated by the experimental method were 1–3 orders of magnitude larger than that predicted by the parameterization methods except for the Peräkylä et al. method. The Peräkylä et al. method provides the most similar predictions to the experimental method for isoprene oxidation products in the full volatility range, with the disagreement within 1 order of magnitude.
In general, the vapor pressures estimated experimentally in this study are very close to those calculated by the Peräkylä et al. method for which the estimation parameters were also derived experimentally. The discrepancy between the experimental and the GC methods span several orders of magnitude depending on different compounds, with the GC methods predicting higher for less-functionalized 1N-monomers and approximate for 2N-monomers, but lower for highly functionalized dimers. It is difficult to tell which method is more reliable without any measured saturation vapor pressure data on such multifunctional organic nitrates. However, considering the fact that the existing GC methods tend to underestimate saturation vapor pressures of the highly functionalized organic molecules due to their limited capability to deal with intramolecular interactions (e.g., the intramolecular hydrogen bonding formed among polar functional groups), and the consistent results of two experimentally derived methods, we suggest that the experimental method might be a good choice to determine the volatility of highly oxidized compounds accurately.
3.3.3 Volatility distribution of isoprene nitrates and expected SOA yieldsAlthough the vapor pressures calculated by different methods show a variability of several orders of magnitude, the predicted volatility distributions of different organic groups are consistent. To eliminate the discrepancy caused by methods and get an average trend of the volatility distribution of various isoprene nitrates, we use the median value of calculated by different methods as the estimator of the vapor pressure for each nitrate compound.
The average carbon oxidation state is plotted against in Fig. 8 to describe the volatility distribution of organic nitrates formed from isoprene oxidation by NO. Generally, the volatility of measured gas-phase products spans a wide range from IVOCs to ELVOCs, wherein all of the 1N-monomers fall in the IVOC or SVOC range, suggesting that 1N-monomers have low potential to form SOA by simple condensation as long as the organic aerosol load is less than 200 g m. The addition of a second or third –NO functional group decreases of most 2N- and 3N-monomers by 2–3 decades compared with 1N-monomers, and most of them belong to SVOCs. They will start to condense in significant fractions if the organic aerosol load is in a range of 1–10 g m, which means 2N- and 3N-monomers with OS 0.8 may contribute to SOA formation under atmospheric conditions. With regard to dimers, all 3N-dimers and 2N-dimers (except for CHNO) are in LVOC or even ELVOC range, indicating isoprene dimers had high propensity to form SOA even at organic aerosol loads 1 g m. However, we would like to emphasize here that the signals of 2N- and 3N-dimers only account for less than 2 % on average of the total assigned signals, as shown in Fig. S6. This suggests that the SOA yield of isoprene from NO oxidation by condensation should be low under atmospheric conditions.
Figure 8
Volatility distribution of different organonitrates formed from NO-initiated isoprene oxidation. The volatility classes are indicated along the top with corresponding colors in the plot. The position of potential SOA contributors is determined depending on the exact functionalities of molecules adapted from Bianchi et al. (2019).
[Figure omitted. See PDF]
The fate of RO determines the product distribution directly and hence could substantially affect SOA yields and aerosol physicochemical properties (Boyd et al., 2015; Fry et al., 2018; Ng et al., 2008; Schwantes et al., 2015; Ziemann and Atkinson, 2012). Consequently, it would be helpful to provide SOA yields together with the fate of RO. In our experiment, reactions with HO and NO are the dominant loss channels for the initially formed RO from isoprene oxidation by NO, contributing for 53 % and 30 % of overall RO loss; RO RO reactions contributed a minor fraction ( 13 %) followed by unimolecular reactions with a contribution of 5 %, according to modeling results (Brownwood et al., 2021). More details about the modeling and the results can be found elsewhere (Brownwood et al., 2021; Vereecken et al., 2021).
In polluted urban regions, the fate of RO is typically dominated by RO NO, while in the more pristine environment, the RO HO reaction will dominate RO fate (Bianchi et al., 2019; Boyd et al., 2015; Brown and Stutz, 2012). RO HO was more important in the chamber than that in ambient conditions, and enhanced RO HO would potentially lead to less dimer formation by RO RO reactions and hence reduced SOA yields. However, a recent work from Brownwood et al. (2021) based on the same campaign as this study pointed out that the bulk aerosol composition and SOA yields were largely independent of RO fate. Similarly, Boyd et al. (2015) found for the -pinene-NO system that RO fate (“RO NO dominant” vs. “RO HO dominant”) had only few effects on SOA formation. Therefore, the SOA yield estimated in this study is expected to be comparable to that in the atmosphere.
Assuming that the dimers in the LVOC or ELVOC range will condense onto particles, we estimated a SOA mass yield for condensation of isoprene organic nitrates of about 5 % 2 %. This value is based on an averaged bulk organonitrate sensitivity of 0.019 norm. counts s ppbv and has been corrected for wall loss and dilution (see Fig. S7, with uncorrected SOA mass yield of about 2 %). The estimated SOA mass yield is within the range of those reported in the literature, but at the lower end (4.3 % to 23.8 % depending on RO fate, Ng et al., 2008; 0.7 % for first generation oxidation and 14 % after oxidation of both double bonds, Rollins et al., 2009; 27 % on average for ambient measurements, Fry et al., 2018). The SOA yield will probably become somewhat higher if taking the contribution of SVOCs (including CHNO, CHNO, and some other monomers, as shown in Fig. 8) into consideration. Our finding is commensurable with the SOA yield for isoprene organic nitrates of 2 %–6 % derived from HR-AMS measurements in the same campaign (Brownwood et al., 2021).
In addition, Br adduct ionization CIMS is selective for HO and less oxidized organic compounds (Albrecht et al., 2019; Rissanen et al., 2019), so it is reasonable to assume that there were more highly oxidized products that were not detected by Br CIMS. This assumption is confirmed by measurements with a NO CIMS performed in another isoprene-NO experiment in SAPHIR (Zhao et al., 2021). Zhao et al. (2021) observed a higher fraction of dimers and more highly oxidized monomers and dimers, as well as trimers (C compounds). As a consequence, the SOA yields derived from NO CIMS measurements is slightly higher.
From these points of view our yield is more a lower limit. However, even if we assume an error of a factor of 2, the SOA yield of isoprene organic nitrates by condensation is more likely in a range of about 10 % or less than in the higher range of 20 %–30 % published in the literature. Of course, by our method we cannot cover any liquid phase processes that would lead to additional SOA beyond the condensation of the target compounds.
4 Conclusions and implicationIn this work, a gas-phase experiment conducted in the SAPHIR chamber under near-atmospheric conditions in the dark was analyzed to primarily investigate the multi-generation chemistry of isoprene-NO system. The characteristics of a diversity of isoprene nitrates were measured by the CIMS using Br as the reagent ion. Isoprene 1N-,2N-, and 3N-monomers and 2N- and 3N-dimers have different time behaviors, indicating the occurrence of multi-generation oxidation during this process. Based on their specific time behaviors as well as the general knowledge of isoprene and radical chemistry, the possible formation mechanisms of these compounds are proposed.
In order to evaluate the potential contribution of various isoprene nitrates to SOA formation, different composition-activity and group-contribution methods were used to estimate their saturation vapor pressures. We also calculated the vapor pressures of isoprene oxidation products based on the gas-particle equilibrium coefficients derived from condensation measurements. The vapor pressures estimated by different methods spans several orders of magnitude, and the discrepancies increase as the compounds become highly functionalized. It shows that existing composition-activity methods (especially the Donahue et al. method) seriously underestimate the saturation vapor pressure of multifunctional low-volatility molecules compared to the experimental methods. The group-contribution methods seem to have a better performance than the CA methods in this aspect, but they still have a tendency to underestimate the vapor pressures of multifunctional molecules. We suggest that experimental methods are a good choice to estimate the volatility of highly oxidized compounds accurately.
According to our results, 1N-monomers and most 2N- and 3N-nitrates fall in the IVOC or SVOC range. Therefore, they have, with a few exceptions, low potential to form SOA at atmospheric organic aerosol loads. In contrast, 2N- and 3N-dimers are estimated to have low or extremely low volatility, indicating that they are significant contributors to SOA formation, although dimers constitute less than 2 % of the total explained signals. In this study, no new particle formation events were observed. Assuming that the dimers in the LVOC or ELVOC range will condense onto particles completely, we estimate a SOA mass yield of about 5 % 2 %, which is a lower limit if one takes a possible contribution of the minor dimer products as well as SVOC species into consideration. Both the volatility distribution and calculated SOA yields indicate that isoprene dimers formed from NO oxidation are the major contributors to SOA formation.
Data availability
All data given in figures can be displayed in tables or in digital form.
This includes the data given in the Supplement. Please send all requests for
data to [email protected] and [email protected]. The data used in this work are available on the EUROCHAMP database
(
The supplement related to this article is available online at:
Author contributions
HF, JNC, JLF, SSB, AW, and AKS designed the study. Instrument deployment and data analysis were carried out by RW, ET, SK, SRA, LH, AN, HF, RT, TH, PTMC, JS, FB, BB, and JAT. RW, LV, ET, DZ, JAT, MH, and TFM interpreted the compiled data set. RW, TFM, and LV wrote the paper. All co-authors discussed the results and commented on the paper.
Competing interests
The authors declare that they have no conflict of interest.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Special issue statement
This article is part of the special issue “Simulation chambers as tools in atmospheric research (AMT/ACP/GMD inter-journal SI)”. It is not associated with a conference.
Acknowledgements
Rongrong Wu gratefully acknowledges the fellowship from Helmholtz-OCPC (Office of China Postdoc Council) Postdoc Program for research support.
Financial support
This research has been supported by the European Research Council (ERC), (SARLEP grant agreement no. 681529), European Commission (EC) under the European Union's Horizon 2020 research and innovation program, (Eurochamp 2020 grant agreement no. 730997 and FORCeS grant agreement no. 821205), Vetenskapsrådet (VR, grant nos. 2014-05332 and 2018-04430), and Svenska Forskningsrådet Formas (grant nos. 2015-1537 and 2019-586). The article processing charges for this open-access publication were covered by the Forschungszentrum Jülich.
Review statement
This paper was edited by Eliza Harris and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2021. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Isoprene oxidation by nitrate radical (NO
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details















1 Institute of Energy and Climate Research, Troposphere (IEK-8), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany; College of Environmental Sciences and Engineering, State Key Joint Laboratory of Environmental Simulation and Pollution Control, Peking University, 100871, Beijing, China
2 Institute of Energy and Climate Research, Troposphere (IEK-8), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
3 Department of Chemistry and Molecular Biology, University of Gothenburg, 41296, Gothenburg, Sweden
4 Institute of Energy and Climate Research, Troposphere (IEK-8), Forschungszentrum Jülich GmbH, 52428 Jülich, Germany; present address: SOLIDpower GmbH, 52525 Heinsberg, Germany
5 Department of Atmospheric and Oceanic Sciences & Institute of Atmospheric Sciences, Fudan University, 200438, Shanghai, China
6 Atmospheric Chemistry Department, Max Planck Institut für Chemie, 55128 Mainz, Germany
7 Institut de Combustion, Aérothermique, Réactivité et Environnement (ICARE), UPR CNRS, 45071 Orléans, France
8 Department of Chemistry, Reed College, Portland, OR 97202, USA
9 Department of Atmospheric Sciences, University of Washington, Seattle, WA 98195, USA
10 NOAA Chemical Sciences Laboratory, Boulder, CO 80305, USA; Department of Chemistry, University of Colorado, Boulder, CO 80309, USA