the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 28 Mar 2024
Ether and ester formation from peroxy radical recombination: A qualitative reaction channel analysis
Abstract. The least volatile organic compounds participating in atmospheric new-particle formation are very likely accretion products from self- and cross-reactions of peroxy radicals (RO2). It has long been assumed that the only possible accretion product channel in this reaction is that forming a peroxide (RO2 + RO2 → ROOR + O2), but it has recently been discovered that a rapid alkoxy radical (RO) decomposition may precede the accretion step of the mechanism, forming slightly fragmented but more stable ether (ROR) or ester (RC'(O)OR) accretion products. In this work, the atmospheric implications of this new reaction channel have been explored further using a modified version of the GECKO-A software to generate a large amount of representative RO2 + RO2 reactive pairs formed from the oxidation of typical primary hydrocarbons, and applying Structure-activity relationships (SAR) to predict the potential accretion products. This data is analysed in terms of formation of low-volatility products, and new discoveries are presented on what kind of RO2 are especially efficient (and which are surprisingly inefficient) at forming accretion products. These findings are discussed in terms of atmospheric relevance of these new RO2 + RO2 reaction channels. As the generation of this data rests on several simplifications and assumptions, many open questions worthy of later studies are also raised.
- Preprint
(3563 KB) - Metadata XML
-
Supplement
(2176 KB) - BibTeX
- EndNote
Status: open (until 09 May 2024)
-
RC1: 'Comment on egusphere-2024-920', William Carter, 08 Apr 2024
reply
This paper discusses an analysis of the formation of accretion products from the cross-reactions of peroxy (RO2) radicals formed in atmospheric reactions of several representative compounds using a modified version of the Gecko-A atmospheric chemical mechanism generation system. Previously it was thought that the major pathways for such cross-reactions was either dissociation forming alkoxy radicals or H-shifts (when possible) forming an alcohol and carbonyl, but addition forming peroxides + O2 has also been considered. More recently, experimental and theoretical studies indicated that accretion reactions forming more stable (yet lower molecular weight) ethers and esters may also occur in some cases. This is investigated in this work, where the Gecko-A system is modified to include estimates of these accretion reactions for peroxy radicals formed from representative compounds. A number of the estimates are uncertain and approximations had to be used in this analysis because it is not possible in practice to represent all the possible peroxy + peroxy combinations in explicit mechanisms, and because the levels and distributions of peroxy radicals formed, depend significantly on the environment where they react. Therefore, as the authors admit, the results are largely qualitative, and do not necessarily represent any particular environment. However, this appears to be the first attempt to estimate the distributions formation of accretion products in reactions of organics in atmospheric systems, which are potentially important sources of secondary organic aerosol in the atmosphere that are not represented in current models. For this reason, I consider this work to be of high scientific significance and appropriate for publication in this journal.
Although the methods they used to estimate the branching ratios for the peroxy + peroxy reactions have uncertainties, especially for the larger radicals where accretion reactions may be more significant, I believe their chemical mechanism estimation methods represent the current state of the science and are probably as good as can reasonably be done given our current state of knowledge. Unfortunately, the Gecko-A tool that they used was not ideally suited for this study because it did not consider H-shift auto-oxidation isomerization reactions of peroxy radicals, which will affect predictions of distributions of peroxy radicals formed and which last long enough to react with other peroxy radicals. The authors recognize this limitation and attempts are made to deal with it by removing at least one type of radical where this is known to be fast from their analysis (acyl peroxy radicals with ‑OOH groups). However, this isomerization is either fast or at least non-negligible for radicals with -CHO groups or with allylic hydrogens, and these are formed from many compounds, including some examined in this study. Fortunately, structure and reactivity (SAR) methods now exist for H-shift reactions of a wide variety of peroxy radicals (e.g., Vereecken and Nozière, 2020, https://doi.org/10.5194/acp-20-7429-2020), and those could be incorporated in mechanism generation systems to account for this type of reaction. It would have been better, and certainly within our current state of knowledge, had they also modified Gecko-A to include these auto-oxidation reactions as well as the accretion reactions. As it is, a study like this would need to be repeated once this modification is made.
The greatest problem I have with this work is that it did not provide any indication of the overall importance of these peroxy + peroxy reactions in atmospheric systems, or how the estimated accretion product yields and distributions might depend on environmental conditions. Environmental conditions are important in affecting both the importance of peroxy + peroxy reactions compared to competing reactions with NOx and HO2, and also the distribution of peroxy radicals present. Peroxy + peroxy reactions can be important in laboratory systems in experiments with high concentrations and no NOx, but are negligible at high or moderate NOx levels characteristic of urban daytime scenarios, and tend to be much less important than reactions with HO2 under lower NOx conditions in the atmosphere. Therefore, the significance of these reactions as a source of low-volatility compounds in the real atmosphere is uncertain. Unfortunately, this work provided no useful information in this regard, which yields given only relative to the total amounts of peroxy + peroxy reactions.
A better approach for this study would be to define a set of standard environments where levels of NOx, HOx, O3, peroxy radicals, light intensity and spectrum and reaction times are specified, so that total yields of accretion and all other products, relative to the amount of compound reacted, can be unambiguously determined. This gives information not only on product yields and distributions, but also how they depend on the environment if more than one standard environment is used. For the purpose of this study, it is also necessary to specify a distribution of peroxy radicals that serve as co-reactants for the peroxy radicals formed from the compound of interest, unless the scenario being simulated has only one organic reactant initially present. That may be a good representation of laboratory experiments carried out to determine product yields, but is a poor representation of the atmospheric environment where a wide variety of compounds are reacting, most being lower the molecular weight compounds than compounds of interest for SOA formation. However, if a peroxy radical distribution that is reasonably representative of the environment is specified as an input, with Gecko calculating only the peroxy radicals formed from individual compounds of interest, then the modified Gecko simulations could give quantitative predictions of types of accretion products formed from the compounds in various environments. It may not be practical to use Gecko to simulate complex atmospheric mixtures, but the MCM mechanism, which explicitly represents many of the most important types of radicals formed, should be sufficient for this purpose.
However, the use of scenarios with simple mixtures, as employed in two of the three sets of calculations presented in this study, is not a good representation of either the atmosphere or of the types of laboratory experiments generally used to study product formation from individual compounds. The other set of calculations used beta-caryophyllene alone, which may be a good representation of product yield experiments with that compound, but is a very poor representation of real atmospheres.
I found their discussion of how they estimated product yields in the approach they employed to be unclear and questionable. The discussion of how yields are determined for the purpose of determining which radical pairs can be neglected, and for the purpose of reporting relative yields in the various figures, is unclear and needs to be improved. On lines 110-115 they suggested that a maximum yield of 100% was used for competitive bimolecular reactions, while in Section 2.1.2 they discuss calculating relative yields based on rate constants and levels of co-reactants. I'm not sure the approach they discuss in Section 2.1.2 is valid for acyl peroxy radicals, since reaction with NO2 is not a final sink for these radicals, since the PAN compounds formed in the acyl peroxy + NO2 reaction can decompose and re-form these radicals. What is needed is presentation of simple examples in the SI for exactly what calculations are used to judge whether reactions of a peroxy pair should be neglected, and for how the relative yields that are presented in the figures are derived. However, if they used actual yields in calculations with standard environments, as suggested above, then all this complexity and questionable treatments regarding estimated yields can be eliminated.
I could not locate a discussion of how they handle cases where both of the reacting peroxy radicals form alkoxy radicals with rapid decompositions. If they both form alkoxys with rapid b-scission reactions that form alkyl radicals, could they react together and form accretion products with new C-C bonds, rather than ethers or esters?
There are several other concerns or issues I have with this paper. These are summarized below, in approximate order they appear in the manuscript and supplement.
Main Manuscript:
I do not understand why the fact that several reaction channels may form the same accretion products, mentioned in Section 2.2.3, is a problem. Certainly Gecko-A must have the code needed to tell whether a molecule formed in one reaction is the same as that formed in another. For the purpose of this paper, I think this should be sufficient for the code to compute the number of products formed without double counting. Or am misunderstanding something in this section?
In Table 3 in the "Statistics" section, it indicates that the first generation n-decane reactions could form an acyl peroxy radical in one generation. I don't understand how this could be the case. In addition, the column headers on Table 3 need to be improved. How are the numbers in the "yield" column computed and what are they relative to? It can't be relative to the amounts of compounds reacted because the yields of peroxy radicals on OH or NO3 reactions should be close to 100%, or higher in cases where consecutive reactions are important. If they are relative to the amount of reactant present, then they would depend on the reaction time and the levels OH, O3, or NO3 present.
I am assuming that the yields shown plots shown on Figures 2-5 are relative to the total amount of accretion products formed from a given compound, since it looks like they may sum up to 100%. If so, this should be stated explicitly. If not, the meaning of yields in these figures, and what they are relative to, need to be clarified. Note that as mentioned above these relative yields give no useful information on how important these products may actually be in the atmospheric environment, which might be very low or negligible for many conditions.
Figures 2b, 2d, 5b, 5d, and perhaps others seem to have more colors on the bar graphs than there are in the legends. If lighter and darker shades of the same color on the bar graphs have significance, then this needs to be mentioned.
In figure 3, it is very difficult to determine which structures are being plotted because the font of the group designations on the X axis is so small it is unreadable without greatly expanding the screen. This figure needs to be redone so that the main information, the type of group, can be more clearly seen. In addition, the plot labels should have "terpene and beta-caryophyllene peroxy radicals", not "terpene peroxy radicals".
Line 424 in section 3.3 states that "the code insures that peroxy radicals with an NO2 group in the alpha position never forms accretion products". This should be corrected to state that the NO2 is in the beta position. It may be in the alpha position of the alkyl radical formed in the decomposition, but it is in the beta position in the peroxy radical.
Why are the product yields shown on the volatility distribution plots on Figures 9b and 9d be almost an order of magnitude lower than on the analogous plots of carbon or oxygen number distributions on Figures 2-5? Are they relative to a different quantity? Note that the volatility distributions seem to include compounds of all volatilities, so the total amounts should be the same as in the atom number plots on the earlier figures.
Supplement:
As indicated above, there needs to be a section giving a more complete description of the algorithm(s) used to derive yields for various purposes, and it should include examples of such calculations. The examples would need to as simple as possible while incorporating all the types of situations that may need to be considered. With an adequate discussion of yields in the supplement, some of the complexity in sections 2.1.2 and 2.1.5 could well be moved to the supplement, so the main paper focuses more on the chemistry.
Section S1 is useful, and is sufficiently important that it may be appropriate to include some of this discussion in the main text. This includes possibly Figure S3, which I found to be interesting. However, I think Figure S3 could be improved by using only one plot with different symbols for the complexes with multiple stereoisomers, and also include the 35 kJ/mole cutoff. The line fit should only consider the points below or near that cutoff.
I found the bubble plots on Figure S19 sufficiently interesting that maybe this figure should also be promoted to the main text. Similar bubble plots comparing C and O numbers would also be interesting.
Citation: https://doi.org/10.5194/egusphere-2024-920-RC1
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 193 | 56 | 9 | 258 | 19 | 4 | 7 |
- HTML: 193
- PDF: 56
- XML: 9
- Total: 258
- Supplement: 19
- BibTeX: 4
- EndNote: 7
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1