Selective deuteration as a tool for resolving autoxidation mechanisms in α -pinene ozonolysis

. Highly oxygenated organic molecules (HOMs) from α -pinene ozonolysis have been shown to be signiﬁcant contributors to secondary organic aerosol (SOA), yet our mechanistic understanding of how the peroxy-radical-driven autoxidation leads to their formation in this system is still limited. The involved isomerisation reactions such as H-atom abstractions followed by O 2 additions can take place on sub-second timescales in short-lived intermediates, making the process challenging to study. Similarly, while the end-products and some-times radical intermediates can be observed using mass spectrometry, their structures remain elusive. Therefore,


Introduction
Highly oxygenated organic molecules (HOMs) are atmospheric compounds with low to extremely low volatilities formed from volatile organic compounds (VOCs) in the atmosphere via autoxidation (Crounse et al., 2013;Ehn et al., 2014;Bianchi et al., 2019).HOMs formed from α-pinene ozonolysis have been shown to contribute significantly to secondary organic aerosol (SOA) (Ehn et al., 2014); however, a mechanistic understanding of this system is limited, despite it having been studied computationally and experimentally (Kurtén et al., 2015;Iyer et al., 2021).Autoxidation involving peroxy radicals (RO 2 ) has been shown to be the source of the high oxygen content (Berndt et al., 2018), but the exact route or routes of isomerisation cannot be derived purely by the mass spectrometric methods (Jokinen et al., 2012;Riva et al., 2019a) typically used to detect the HOMs.In addition, the isomerisation reactions are very fast (Bianchi et al., 2019), which makes following them in real time typically not possible.One way to circumvent these limitations is to perturb the autoxidation process in a way that leads to different HOM products depending on the reaction mechanism.
Deuteration is a commonly used method for multiple applications in chemistry, including atmospheric studies (Rissanen et al., 2014;Zhang et al., 2017;Ye et al., 2018;Michelotti and Roche, 2019;Wang et al., 2020;Li et al., 2021).In deuteration, one or more hydrogen atoms in a molecule are exchanged to deuterium atoms (D).Deuteration can be a full deuteration where all hydrogen atoms are exchanged to deuterium atoms as used by Rissanen et al. (2014), or it can be done selectively, where specific hydrogen atoms are exchanged with deuterium atoms as used by Zhang et al. (2017) and Ye et al. (2018).The deuteration can have several impacts on the reactions taking place, including a decreased rate of reaction due to the higher mass of a D atom compared to a H atom, the so-called kinetic isotope effect (Laidler, 1987).Rissanen et al. (2014) used fully deuterated cyclohexene and found that the ozonolysis formed HOMs at greatly reduced yields compared to the non-deuterated cyclohexene.They also showed that D atoms were exchanged to H atoms in contact with water vapour in cases where the C−D bond was broken, and the D became attached to an oxygen through an O−D bond.In other words, the autoxidation was perturbed in predictable ways by the deuteration.Nevertheless, we expect that selective deuteration can provide substantially more information still.To our knowledge, selective deuteration has not been utilised in studying autoxidation pathways from α-pinene ozonolysis before.
In this study, we inspected α-pinene ozonolysis by simulating atmospheric conditions in a reaction chamber.We used non-deuterated as well as three selectively deuterated precursors to test how the resulting HOM spectrum might change and to gain insight into the autoxidation processes.We used a chemical ionisation Orbitrap mass spectrometer (Riva et al., 2019a) which has a high-enough resolving power to unambiguously distinguish between the different isotopes in the spectra.We compared the spectral distribution of different HOMs as well as the overall HOM yields for the different precursors.In particular, we examined the numbers of D atoms lost for the differently deuterated precursors for different HOM molecules, with the aim of assessing from which C atom H/D abstractions take place during autoxidation.We also compare our findings to the few suggested mechanisms from earlier publications.

Selective deuteration and autoxidation
In this section, we first describe the selectively deuterated precursors used in this study (Sect.2.1), after which we describe the main reaction pathways during α-pinene ozonolysis that are of relevance for this study (Sect.2.2).Lastly, we describe how selective deuteration can be used to study the oxidation processes (Sect.2.3).

Selectively deuterated precursors
In this study, we used four precursors, and their chemical structures are shown in Fig. 1.The three selectively deuterated precursors were synthesised at Northwestern University, Illinois, USA.The precursors are named based on the carbon that was deuterated and the number of deuterium atoms they have: D 0 is the standard α-pinene molecule without any isotopic alteration (C 10 H 16 ), 3 D 1 is α-pinene with one deuterium atom on C-3 (C 10 H 15 D, α-pinene-3-d1; see Fig. 1), 7 D 2 is α-pinene with two deuterium atoms on C-7 (C 10 H 14 D 2 , α-pinene-7,7-d2), and 10 D 3 is α-pinene with three deuterium atoms on C-10 (C 10 H 13 D 3 , α-pinene-10,10,10-d3).The vinylic C-3 carbon and the allylic C-10 methyl group carbon are expected to undergo changes during the ozonolysis process (Jenkin et al., 1997;Kurtén et al., 2015), while the cyclobutyl-ring C-7 carbon has not to our knowledge been shown to take part in any reactions in the initial ozonolysis reaction.
Synthesis of the 10 D 3 and 7 D 2 α-pinene species was carried out according to the published literature procedures (Upshur et al., 2016(Upshur et al., , 2019)), while access to the 3 D 1 α-pinene was accomplished through a new three-step route utilising nopinone as the starting material (see the Supplement for complete synthetic details).The purity of each compound was > 95 % as determined by 1 HNMR spectroscopy, and we were unable to observe residual proton resonances associated with the deuterated carbon positions for any sample, indicating complete deuterium incorporation within the detection limits of the spectrometer (see the Supplement for 1 HNMR spectral data).

Autoxidation in α-pinene ozonolysis
As discussed above, autoxidation is a major source of condensable low-volatile vapours from α-pinene ozonolysis.The reaction of α-pinene with ozone starts by the attachment of ozone to the double bond, followed by bond scission and rapid isomerisations via a Criegee intermediate and vinyl hydroperoxide and OH loss, potential isomerisation, and O 2 addition to form a primary peroxy radical (C 10 H 15 O 4 ), as shown in Fig. 2a (Kurtén et al., 2015;Iyer et al., 2021).From these peroxy radicals, the process can continue with autoxidation, which is successive H shifts and oxygen additions.For pathways I to III, Kurtén et al. (2015) calculated that subsequent H shifts starting from these primary RO 2 were too slow to explain observed HOM formation rates.This was primarily due to the intact cyclobutyl ring causing strain for the abstraction of hydrogen atoms across the ring.However, Iyer et al. (2021) found that an isomerisation reaction was available for the vinoxy intermediate in pathway I due to excess energy remaining from the initial ozonolysis reaction.The excess energy allows the breaking of the cyclobutyl ring (pathway IV) and subsequent rapid H shifts leading to RO 2 with up to eight oxygen atoms (Fig. 2b).The O 8 −RO 2 product remains the only experimentally and computationally supported HOM-forming pathway in this system (Rissanen et al., 2014;Iyer et al., 2021) and is thus the most promising candidate to explain many of the observed closed-shell HOMs following termination reactions.
The autoxidation process can be terminated via unimolecular or bimolecular reactions, and the most relevant ones are introduced below.A common unimolecular Reaction (R1) of an RO 2 that has undergone autoxidation is the H abstraction from a C atom with a hydroperoxide functionality.This rapidly leads to the loss of a hydroxyl radical and the for-mation of a closed-shell HOM with a carbonyl functionality at said C atom (Crounse et al., 2013).The RO 2 can also undergo different bimolecular termination reactions (Atkinson and Arey, 2003).In the absence of NO x , as in this study, there are two relevant reaction partners, namely HO 2 and other RO 2 , and their main reaction products are listed below Reactions (R2)-(R6).In brief, RO 2 cross-reactions can form closed-shell monomers (Reaction R2), ROOR accretion products which are commonly referred to as dimers (Reaction R3), or alkoxy radicals RO (Reaction R4), with the latter able to undergo further isomerisation or termination reactions.Reactions of RO 2 with HO 2 typically form hydroperoxides ROOH (Reaction R5) but can also form RO (Reaction R6), especially for complex-enough RO 2 (Hasson et al., 2005;Orlando and Tyndall, 2012;Iyer et al., 2018).

Studying autoxidation with selective deuteration
During the oxidation and autoxidation processes, hydrogen atoms can be lost from the precursor molecule in multiple ways.If a carbon loses a hydrogen atom during an H shift, the carbon has actively participated in the autoxidation process, and we refer to this as the carbon being "active".More specifically, when discussing the experimental results, we only consider the carbon to be active when we observe the D abstractions taking place at fractions above our uncertainty margins.While we are primarily interested in H shifts by RO 2 , they can also occur by RO radicals or in the Criegee intermediate step, a reaction that takes place in all cases.When a hydrogen atom is shifted away from a carbon atom, the hydrogen atom ends up in a hydroperoxide group −OOH, except in the case of H shifts in RO, where the result is a hydroxyl group −OH.In both cases, the H atoms are now labile enough that they can be exchanged in collisions with water vapour molecules.In the special case where a deuterated carbon has been active and a D shift has taken place, the D can eventually be lost from the molecule entirely (example shown in Fig. 2c).With deuteration, we can thus observe the change in the molecular mass when the D/H exchange has taken place and thus specify which carbon in the precursor has been active.Additionally, if a deuterated precursor produces a mass spectral signature for a given HOM molecule where part of the signal has lost a D while another has not, it is an indication of different isomers for this composition.In other words, several pathways lead to the same elemental composition but with different C atoms being active.If the D atom was internally abstracted by the Criegee intermediate, it will be lost by ejection of the OD, as would be the case in pathway II (Fig. 2a, blue H) for the 10 D 3 precursor.Furthermore, deuterium atoms can be lost if the selectively deuterated carbon itself is lost through fragmentation during the HOM formation reaction.
Only two of the four pathways discussed earlier result in the precursor losing a deuterium atom.In the four pathways shown in Fig. 2a-b, the carbon and hydrogen atoms corresponding to those of the selectively deuterated precursors are highlighted with colours matching the colouring from Fig. 1.The atoms coloured with cyan correspond to the 3 D 1 precursor, the green ones correspond to the 7 D 2 precursor, and the blue ones correspond to the 10 D 3 precursor.In pathway II, the 10 D 3 precursor would lose one deuterium when OH is fragmented from the molecule, and in pathway IV with the continued steps in Fig. 2b, the deuterium in the 3 D 1 precursor would be shifted to a hydroperoxy group and hence readily exchanged to a hydrogen atom with water.
In addition to the deuteration impacting the molecular mass observed in the mass spectra, the H shifts themselves will become much slower since the deuterium atom being abstracted has twice the mass of a hydrogen atom.Rissanen et al. (2014) found that the HOM yields were roughly 100 times lower for ozonolysis of fully deuterated cyclohexene C 6 D 10 compared to that of standard cyclohexene C 6 H 10 .Some of the same effects should play a part in this oxidation system, namely the kinetic isotope effect and the slowing down of quantum tunnelling.The kinetic isotope effect is caused by the lower vibrational frequencies for the C−D bond compared to the C−H bond, and it decreases the rate of the abstraction (Laidler, 1987).Furthermore, quantum mechanical tunnelling is important for the H/D abstractions, and this causes an even larger decrease to the reaction rate (Laidler, 1987).
The exact total decrease in the reaction rate coefficients depends on the radical structure, but a typical range is a factor of 20-100 slower for D shifts versus H shifts (Crounse et al., 2011;Praske et al., 2018).However, it is important to keep in mind that the effect on the likelihood of the D shift taking place will still depend on the possible competing reactions.For example, if the H shift is 1000 times faster than any competing reaction, substituting for a D shift will barely impact the branching of the reaction, whereas with similar H-shift rates with any competing reaction, deuteration will largely shut down this pathway.

Methods
In this section, we first describe the experimental set-up and the used instruments (Sect.3.1).Then we describe how the data were analysed (Sect.3.2).

Experimental set-up
A schematic of the experimental set-up is shown in Fig. 3.The experiments were conducted in the 2 m 3 "COALA" chamber (Teflon FEP chamber, Vector Foiltec, Germany) (Riva et al., 2019b;Peräkylä et al., 2020), a continuous stirred-tank reactor (CSTR) at the University of Helsinki, Finland.We monitored relative humidity (RH) and temperature with a temperature and humidity probe (INTERCAP ® HMP60, Vaisala).RH stayed below the detection limit of the probe (under 1 %) and temperature was 26.0 • C within the ±1.5 • C range.We monitored NO x (NO and NO 2 ) concentrations with a NO-NO 2 analyser (model T200UP, Teledyne), and the concentrations were under 0.3 ppb.We generated ozone by running 5 L min −1 clean air through an ozone generator (Dasibi 1008-PC ozone generator).We ran the chamber at a stable O 3 concentration (∼ 43 ± 2 ppb), monitoring the ozone concentration with a photometric O 3 analyser (model 400, Teledyne instruments; ozone monitor in Fig. 3).The total flow through the COALA chamber was kept at 50 L min −1 by adding 45 L min −1 flow of clean air generated in a zero-air generator (AADCO, Series 737-14, Ohio, USA).As a result, we achieved an average residence time of τ R ≈ 40 min inside the chamber.With these conditions, we estimate the α-pinene loss rate [α − pinene] loss to be 23 % using the equation where k 1 = 9.06 × 10 −17 cm 3 s −1 is the reaction rate coefficient for α-pinene ozonolysis at T = 26 • C (Atkinson et al., 2006).
During an experiment, we aimed to keep the precursor concentration stable when injecting the precursor into the chamber.We ran the chamber as a "steady-state chamber", meaning we continuously fed the precursor and O 3 into the chamber throughout the experiments, aiming for stable concentrations resulting from a balance between sources (inflow) and sinks (outflow, chemical reactions, and wall losses).As such, the absolute precursor consumption varied with the amount of precursor that was fed, which often varied even across one experiment.We injected precursor into the chamber using one of two methods: an overflow set-up for small amounts of precursor and a syringe pump set-up when there was enough precursor to use a 5 µL syringe to take a sample.In both cases, the evaporated α-pinene was injected into the chamber with a small 25-150 mL min −1 N 2 flow.We only had enough of one deuterated precursor ( 3 D 1 ) to be injected with the syringe pump set-up.With the other deuterated precursors ( 7 D 2 and 10 D 3 ), there was barely any visible liquid in the glass bottle, and thus we could not draw the precursors into the syringe to use with the syringe pump.Nevertheless, flushing N 2 over the bottle did release enough precursor to provide a spectrum of oxidation products, though the overflow set-up was harder to control and produced a less stable concentration in the chamber than the syringe pump set-up.The α-pinene injection set-ups are described in more detail in the Appendix (Sect.A1).
We monitored the α-pinene ozonolysis products with a NO − 3 chemical ionisation Orbitrap (CI-Orbitrap).The Orbitrap (Q Exactive Plus Orbitrap, Thermo Scientific) has an ultra-high mass resolution (280 000 Th/Th), meaning it is possible to separate between two hydrogen atoms (m H 2 = 2.0157) and a deuterium atom (m D = 2.0141) signal in the mass spectrum, hence enabling the use of selectively deuterated precursors.Coupled with the nitrate chemical ionisation inlet (NO − 3 CI inlet) which is selective towards highly oxidised products (Hyttinen et al., 2015), the instrument has been shown to be effective at detecting HOMs (Riva et al., 2019a;Bianchi et al., 2019).
The precursor concentration was monitored using the proton-transfer-reaction time-of-flight (PTR-ToF) mass spectrometer (Ionicon Analytik GmbH, Austria; Yuan et al., 2017).However, due to instrument failure, we conducted the last experiment with the 3 D 1 precursor using the Vocus PTR-ToF (Tofwerk AG/Aerodyne Research, Inc.; Krechmer et al., 2018) instead of the PTR-ToF.Nonetheless, both instruments are PTR-ToF mass spectrometers that can be used to monitor VOCs such as monoterpenes (Yuan et al., 2017), and they were calibrated with the same method and hence could fulfil the experimental purposes.
We used a 10 L min −1 sample flow rate with a 60 cmlong sampling line for the CI-Orbitrap, a 1 L min −1 sampling flow rate with a roughly 2.5 m sampling line for the ozone monitor, and a 5 L min −1 sampling flow rate with a roughly 1 m sampling line that was also split to house exhaust for the PTR-ToF mass spectrometers.As a result, we had a 16 L min −1 flow for sampling, and the excess flow from the chamber was fed to the house exhaust, keeping the chamber at a slight over-pressure to ensure contaminations from the room air could not leak into the chamber.

Data analysis
We preprocessed the raw data generated by the CI-Orbitrap using a new software called Orbitool (Orbitool version 1.40) developed specifically for the data analysis of time series data from Orbitrap (Cai et al., 2021).We denoised, mass calibrated, and time-averaged the data into 30 min steps using Orbitool.We then identified peaks corresponding to HOMs found in the atmosphere (Ehn et al., 2012) (Bianchi et al., 2019), but in this study we aimed to keep the data set more compact for ease of analysis by only inspecting the most prominent and atmospherically relevant signals.
We preprocessed the PTR-ToF and Vocus-PTR raw data using tofTools, a MATLAB-based software package (Junninen et al., 2010).Additionally, we conducted sensitivity calibrations to the PTR instruments to obtain calibration factors, which we used to determine the precursor concentrations from the ion signal.The calibration factors were 15 cps ppb −1 for the PTR-ToF used for the D 0 , 7 D 2 , and 10 D 3 precursors and 180 cps ppb −1 for the Vocus-PTR used for the 3 D 1 precursor.
Because the Orbitrap's sensitivity as a function of signal intensity is non-linear, we determined our instrument's sensitivity behaviour by analysing the α-pinene oxidation products and their natural isotopes following the method of Riva et al. (2020).When determining the instrument sensitivity, we excluded all isotopes that contained deuterium, because using selectively deuterated precursors distorts the signal distributions of the natural hydrogen isotopes.Based on our instrument's sensitivity behaviour, we created a correction function to correct the time series data for the nonlinearity.We also used the sensitivity behaviour to determine CI-Orbitrap's detection limit, which is the signal intensity normalised with the total reagent ion signal intensity above which the signal is ≥ 20 % of the true intensity.The detection limit is 10 −6 cps cps −1 , and it can be converted to a concentration of 10 4 molecules cm −3 by multiplying the signal intensity by the calibration coefficient c; however, the conversion introduces large uncertainties.As calibrations for HOMs remain limited due to a lack of suitable standards, we assumed that the calibration coefficient has the value c = 10 10 molecules cm −3 as an average of typically reported values (Jokinen et al., 2012;Ehn et al., 2014;Riva et al., 2019c).In order to be detected by the nitrate chemical ionisation mass spectrometer (CIMS), the molecule needs to form a very strong cluster with NO − 3 (Hyttinen et al., 2015), and the majority of HOMs are thus expected to be charged at the collision limit, motivating the use of one single c value for all HOMs.We estimate a large uncertainty of at least −67 %/ + 200 % arising from the listed assumptions.However, this uncertainty is of less importance in our study, as absolute HOM concentrations are not used for any of the main conclusions of our work.
Additionally, Riva et al. (2020) found that the relative ion transmission of the Orbitrap mass spectrometer is dependent on the mass of the ion.As we did not separately determine the ion transmission of our instrument during these measurements, we corrected our data assuming the same transmission curve as Riva et al. (2020) determined for their instrument.This assumption generates some uncertainty; however, as we are mainly focused on comparing adjacent ions in the mass spectrum, the transmission is unlikely to make a large difference for the purpose of this study.In contrast, in order to assess the impact of deuteration on the total HOM formation, we did also estimate total HOM yields which could have considerable errors if an erroneous correction is used.Nevertheless, as we will show, using the correction resulted in HOM yields that better matched those of previous studies, while using no correction caused unreasonably high values.Additionally, the HOM yields already have large uncertainties because of the assumptions about the calibration coefficient.
We calculated the total HOM concentration according to where {HOM signal} is the sum of found HOM signals, {total reagent ion signal} is the sum of the measured reagent ion monomer (NO − 3 ) and dimer (HNO 3 • NO − 3 ) signals, and c is the calibration coefficient.
We calculated the HOM yield γ from the production and loss terms of HOMs when the concentration of HOMs [HOM] did not change significantly over time (Ehn et al., 2014): from which we can solve for the yield γ : where k loss = 1 300 s ≈ 3.3 × 10 −3 s is the estimated loss rate for HOMs in the chamber dominated by the loss to the chamber walls (Peräkylä et al., 2020).

Results and discussion
In the following, we first describe the experimental setup and how the oxidation measurements were conducted (Sect.4.1).We then describe how to interpret the data in the form of mass spectra (Sect.4.2).In Sect.4.3, we describe how the results were interpreted and what conclusions could be drawn concerning the involvement of different C atoms.Finally, in Sect.4.4, we evaluate the effects of the deuteration on the overall HOM yields.

Overview of the experiments
We conducted four experiments with three selectively deuterated α-pinenes and standard α-pinene, i.e. one experiment per precursor, Fig. 4 showing their time series.As can be seen from Fig. 4, optimal steady-state conditions were not achieved in all cases, in particular when the injection was performed using the overflow set-up (Fig. A1a), as was the case for 7 D 2 and 10 D 3 .However, we achieved a rate of change in precursor concentration of the order of 10 % h −1 or less in most cases.Thus, the conditions can be considered stable enough from the perspective of the oxidation chemistry, where HOM formation happens on timescales of seconds and HOM loss on timescales of minutes.For each experiment, we calculated average mass spectra for time ranges where the α-pinene reaction rate was comparable between the four precursors.The selected reaction rates k 1 [α-pinene][O 3 ] were 0.5 and 0.015 ppt s −1 , which we will refer to as high-concentration C high and low-concentration C low samples, respectively.The C low sample corresponds to precursor concentration ranges of 1.3-1.7,0.6-1.0, and 1.4-1.8ppb for the D 0 , 3 D 1 , and 7 D 2 precursors, respectively.The C high sample corresponds to 7.6-8.0,4.6-6.0,and 5.0-5.4ppb for the 3 D 1 , 7 D 2 , 10 D 3 precursors.
Due to limitations in the experiments, we could only run one experiment per precursor, which led to us not achieving a low-concentration sample with the 10 D 3 precursor.Similarly, only after the experiment did we conclude that there was no successful C high measurement with the standard αpinene.Thus, we performed the analyses with C high samples when comparing the changes in isotopic distributions due to deuteration, whereas we performed the yield comparisons against the standard α-pinene for the C low for all but the 10 D 3 precursor.
We have to consider the purity of the deuterated precursors in the chamber when analysing the data, as this is another source of uncertainty.The standard α-pinene has a natural isotope distribution where the D/H ratio is around 0.01 %, making D-containing isotopes often fall below the detection limit of our instrument.For the selectively deuterated samples, we can inspect their purity using PTR massspectrometer data.As PTR instruments are known to induce isomerisation and fragmentation of the protonated molecules (Tani et al., 2003;Li et al., 2022), we expect that the observed level of deuteration in the PTR will mainly provide a lower limit for the purity.The PTR measurements suggested that over 88 % of the measured α-pinene in our chamber contained exactly the number of D atoms specified, verifying that the samples are very pure and in line with the observed purity determined using 1 HNMR spectroscopy (see Sect. 2.1).
The uncertainties arising from the limited number of experiments we could conduct, the estimation of the precursor purity from PTR data, and the estimation of the calibration coefficient in addition to other uncertainty sources discussed later limit the extent to which we can draw conclusions from our data.For this reason, we opt for only broadly discussing the implications of our results and focusing on detailed mechanistic insights only when they appear unambiguous.

Interpreting the mass spectra
Interpreting the mass spectra requires some thought and is detailed below for easier reading.Each precursor can form HOM monomers with deuterium atoms up to the number the precursor initially contained.Depending on the formation reactions, the HOMs can have lost deuterium atoms during their formation.As a result, there can be HOM isotopes with the number of deuterium atoms in the molecule anywhere between the number of deuterium atoms in the precursor and zero.Similarly, HOM dimers can contain between twice the number of deuterium atoms in the precursor and zero deuterium atoms.
We show an example of the resulting mass spectra for the HOM C 10 H 14 O 11 when using deuterated precursors in Fig. 5.We have gathered the measured spectra from using each deuterated precursor into one figure.In the mass spectra, the HOMs are clustered with the nitrate ion, and hence the nominal mass for C 10 H 14 O 11 when using normal α-pinene (D 0 ) is 372 Th (310 + 62).We omit the reagent ion from all labels for clarity.
The HOM monomers formed from the 10 D 3 precursor can have zero to three deuterium atoms.As one can see from Fig. 5, where the 10 D 3 signal is shown with blue bars, we found isotopes containing three, two, and one https://doi.org/10.5194/acp-23-4373-2023Atmos.Chem.Phys., 23, 4373-4390, 2023  5).The highest signal by a large margin is C 10 H 12 D 2 O 11 , but in contrast to the 10 D 3 precursor, this HOM has lost zero deuterium atoms during its formation process.Assuming the deuteration does not affect the (preference of) reaction pathways significantly, these findings suggest that in the main pathway forming C 10 H 14 O 11 the C-7 carbon is inactive, whereas the C-10 carbon is active.We emphasise that we only consider the carbon to be active when we can observe the D abstractions taking place at fractions above our uncertainty margins.

HOM isotopic distributions between deuterated precursors
Both radicals and closed-shell products can be found from the α-pinene ozonolysis HOM spectrum.Unfortunately, many of the inspected radical signals were covered under large background and oxidation product signals, making it impossible to say whether the radical signal was completely missing or just hidden under the dominant adjacent peak.For example, the radicals C 10 H 15 O 8 and C 10 H 14 DO 8 were typically not distinguishable due to very large signals of the closed-shell species C 10 H 13 DO 8 and C 10 H 16 O 8 , respectively.Hence, despite the high resolution of the Orbitrap, the significant difference in signal intensities prevented us from using the radical data for further analysis in this study, instead focusing our analysis solely on the closedshell molecules.With lower concentrations and/or shorter residence times, the radical signals would be easier to distinguish as they would account for a larger portion of all the signals (Jokinen et al., 2014;Molteni et al., 2019).We emphasise that the complication arises from the deuteration, as in normal experiments the radicals and closed-shell species would be at different integer masses and thus easily distinguishable.
We focused on compounds that Ehn et al. (2012) found to be the main HOMs in monoterpene-dominated regions like the Hyytiälä forestry field station located in the boreal forest.Due to the sheer number of different signals to look at as there are numerous H/D isotopes for every molecule, we limited the most detailed inspection to eight compounds that correspond to the most prominent 3 ) is at 372 Th.The colours show which selectively deuterated α-pinene is the precursor, and the values above the bars in matching colour tell how many deuterium atoms the HOMs have lost compared to the precursor.For example, the α-pinene with three deuterium atoms ( 10 D 3 , shown in blue) can form up to four different C 10 H 14 O 11 isotopes in the spectrum.However, in the 10 D 3 precursor's case there is no signal for C 10 H 14 O 11 without any deuterium atoms (−3D).(Ehn et al., 2014;Molteni et al., 2019).
By limiting our inspection to these compounds and focusing our main analysis on known ozonolysis products, we limit the effect OH reaction products have on our analysis.Firstly, earlier studies have concluded that HOM formation from OH-initiated reactions is a minor channel in α-pinene ozonolysis, only of the order of 10 % of the ozone-initiated oxidation (Ehn et al., 2014;Jokinen et al., 2014Jokinen et al., , 2015)).As such, the impact of OH reactions on the data would be barely detectable, especially within our high uncertainty range.Secondly, the limited inspection compounds are formed specifically from ozone reactions (Ehn et al., 2014;Molteni et al., 2019), as the OH-derived HOMs would normally have higher hydrogen atom contents and odd numbers of oxygen atoms.Thus, our main analysis specifically focuses on known ozonolysis products, and the effect of OH reactions is further negated.Thirdly, the likelihood of multigeneration oxidation, i.e. first-generation non-HOM oxidation products reacting with OH to form HOM, is quite low in our set-up.This is because the flush-out time of the chamber is only 40 min and the timescale of the wall loss of condensable products is shorter yet (of the order of a few minutes), thus making OH reactions a relatively small sink for the oxidation products in the chamber (Peräkylä et al., 2020;Bianchi et al., 2019).
Isotopic distributions of the four HOM monomers and four HOM dimers for the selectively deuterated precursors are shown in Fig. 6. Figure A1 shows the same distributions for a much wider range of HOMs, but for these compounds each isotope was not checked manually, and thus the results are less reliable than the results in Fig. 6 (see the Appendix for further discussion).The isotopic distributions are shown for C high data.In general, for the 3 D 1 and 10 D 3 precursors, the majority of the signal is generated by isotopes that have lost one or more deuterium atoms, whereas for the 7 D 2 precursor the signals are predominantly from isotopes that have not lost any deuterium atoms.In fact, the signal of HOMs that have lost one or more deuterium atoms for the 7 D 2 precursor is only at most 10 % of the total signal.Due to the large uncertainties in the data, we cannot give an exact value, but the fraction is clearly so small that the pathways involving abstraction from this carbon are likely not significant.In other words, the vinylic carbon C-3 and allylic methyl group carbon C-10 are more active in the formation process than the cyclobutyl-ring carbon C-7.
For the more active C-3 and C-10 carbons, we can also note a clear trend for the C 10 H 14 O x monomers and the C 20 H 30 O x dimers, where deuterium loss increases with increasing oxygen content.This makes sense qualitatively, because in order to reach a higher O-atom content, more H shifts need to take place, and almost all potentially abhttps://doi.org/10.Furthermore, the data for 3 D 1 and 10 D 3 precursors implicitly show that there are multiple structural isomers contributing to the different elemental compositions.If there was only one isomer, the number of lost D atoms would be identical for all molecules with a given composition, resulting in a signal of only one isotope.However, as we can see from Fig. 6 (and Fig. A3), there are signals for multiple isotopes.
We can now compare our results to those of Iyer et al. (2021), where they suggest a mechanism leading to the C 10 H 15 O 8 RO 2 radical shown in Fig. 2b.More specifically, we can look at the closed-shell product C 10 H 14 O 7 formed from the O 8 −RO 2 through OH loss and compare its expected behaviour based on the proposed mechanism to the behaviour observed when using selectively deuterated samples.According to the mechanism proposed by Iyer et al. (2021) (Fig. 2b), we see that, out of the deuterated positions used in this work (C-3, C-7, and C-10), only C-3 has undergone an H shift (cyan H).Hence, we would expect that the closed-shell C 10 H 14 O 7 HOM would always have lost D through exchange with water (Fig. 2c) when using 3 D 1 , while the HOM would never have lost a D atom when using 10 D 3 and 7 D 2 precursors.However, according to our observations, for the 3 D 1 signal of C 10 H 14 O 7 , roughly half still has the D attached (Fig. 6), meaning that the Iyer et al. (2021) pathway can at most explain half of the observed signals in our experiment.For the 10 D 3 C 10 H 14 O 7 , about 20 % has lost a D atom, while the other 80 % behaved according to expectations from Iyer et al. (2021).However, the Iyer et al. (2021) mechanism is not necessarily the only mechanism through which all of the 80 % of the signal is formed, as there can be other mechanisms that could explain parts of the signal.
Lastly, the C 19 H 28 O 11 dimer is often the most abundant dimer measured by the nitrate CIMS from α-pinene ozonolysis (Ehn et al., 2014;Rissanen et al., 2015) and also a dominant signal in night-time atmospheric measurements (Yan et al., 2016), but no suggestions have been made concerning which C atom is lost.In Fig. 6, the signal of this dimer for 10 D 3 precursor has always lost at least three D atoms, providing clear evidence for C-10 being the one that is lost.This is also true more generally for C 9 and C 19 HOMs according to Fig. A3, where we can see a clear difference in the isotopic distributions for the 10 D 3 precursor; i.e. the majority of the C 9 and C 19 signals lost all three D atoms.
To summarise the findings from our approach of studying selectively deuterated α-pinene, the observed roles of different C atoms in the autoxidation process are highlighted in Fig. 7.The carbons C-3 and C-10 (green stars) are active and often lose deuterium atoms during the HOM-formation process.Additionally, the fraction of signal that had lost deu- showing carbons that are active (carbons 3 and 10; green eight-point star) and inactive (carbon 7; light-blue diamond) in the formation process, carbons that were not inspected in this study (carbons 1, 4, 5, 8, and 9; magenta circle), and carbons that cannot be studied using selective deuteration (carbons 2 and 6; grey hexagon).
terium atoms increased with the increasing oxygen number of the HOM.In contrast, the carbon C-7 (cyan diamond) rarely lost deuterium atoms, meaning that it was mostly inactive in the HOM-formation process.As for the carbons C-2 and C-6 (grey hexagons), they can neither be active as defined in this work nor studied with selective deuteration because they are not attached to any hydrogen atoms.The carbons C-1, C-4, C-5, C-8, and C-9 were not inspected in this study, but selectively deuterated studies of these would be extremely interesting and provide crucial information for understanding the detailed autoxidation and HOM-formation mechanisms of this system.

HOM yields
Until now, we have only looked at the relative isotope distributions of specific elemental compositions from separate precursors.However, as described in Sect.2.3, the deuteration may also slow down the autoxidation, which could result in an overall decrease in HOM yields.In this section, we inspect how selective deuterations affected the yields of different HOMs in comparison with HOM yields from standard α-pinene ozonolysis.
Following the definition of HOM in Sect.3.2, the total HOM yields (Fig. 8a) are calculated from the sum of the signals of all found HOMs for the C low samples, except for the 10 D 3 precursor for which the C high sample was used.Not all HOM peaks were necessarily found, especially in the lower masses and with lower signal intensities due to the high background and instrumental limitations.Nevertheless, after correcting the data for the CI-Orbitrap's relative ion transmission and non-linearity of the sensitivity as a func- tion of signal intensity, the HOM yield of the D 0 precursor is roughly 5 %, which is similar to values reported earlier for α-pinene ozonolysis experiments which ranged from 1.2 (+1.2/ − 0.72) to 7.0 ± 3.5 % (Bianchi et al., 2019).Regardless, we again emphasise that there are large uncertainties involved in determining this absolute value, and we are mainly interested in the possible change in the HOM yields with different precursors.This comparison should be much more robust than the determination of absolute HOM yields.
We report the HOM yields grouped by their chemical compositions in Fig. 8a due to their differing behaviour.C 5 H 6 O 7 is separated due to the signal being much higher than the rest of the HOM signals.The yields of HOM monomers with seven O atoms and HOM monomers with eight or more O atoms were separated due to the potentially differing volatilities between the two groups of molecules, the O 7 monomers having a higher volatility (Peräkylä et al., 2020).Lastly, we give one yield value to the rest of the HOMs, namely the HOM dimers.
Surprisingly, the selective deuterations used in this study did not decrease the HOM yields significantly.We can see in Fig. 8a that the obtained total HOM yields for precursors 3 D 1 , 7 D 2 , and 10 D 3 were 3.7 %, 6.3 %, and 3.5 %, respectively.Considering the instrumental uncertainties and the differences in precursor concentrations between the experiments, we conclude that the ±20 %-40 % variability is within the uncertainty limits of our experiment.However, this change can be compared to the change of 2 orders of magnitude observed by Rissanen et al. (2014) for fully deuterated vs. non-deuterated cyclohexene, suggesting that deuteration had a minimal role, if any, in the HOM yields of our precursors.
Furthermore, the detailed spectra from the deuterated precursors were similar to that of standard α-pinene.The compound-wise comparisons of HOMs between D 0 and deuterated precursors are shown in Fig. 8b-d.We calculated the compound-wise yields by using the sum of all found isotope signals corresponding to a specific HOM composition.As an example, to calculate the compound-wise yield for C 10 H 14 O 11 for the 10 D 3 precursor, we take all the found isotope signals corresponding to C 10 H 14 O 11 (blue bars in Fig. 5), sum them, and use that sum to calculate the yield for https://doi.org/10.5194/acp-23-4373-2023Atmos.Chem.Phys., 23, 4373-4390, 2023 the specific compound.What is more, we limited the inspection in Fig. 8b-d to the compounds that were found in both the deuterated precursor's spectrum and the D 0 precursor's spectrum.As can be seen from the figures, the data points fall very close to the 1 : 1 line in all cases, meaning that the total concentrations of HOMs in the deuterated precursors' spectra do not differ significantly from the concentrations of HOMs formed from the D 0 precursor.Additionally, there is no clear trend that more-oxygenated HOMs would be more impacted by the deuterations, although we showed that they typically had lost more deuterium atoms.This may be an indication that the D shifts were still fast enough to outcompete other reaction pathways despite the deuteration.On the other hand, it is possible that the autoxidation could proceed through the next most competitive pathway not shut down by the deuteration and end up losing deuterium atoms later in the process, especially in the case of the more-oxidised products.
The compound-wise HOM yields for 10 D 3 are in general slightly lower than the yields for the D 0 precursor (Fig. 8d), which is mainly caused by us comparing the highconcentration C high 10 D 3 to the low-concentration C low D 0 sample.We saw a similar decrease in HOM yields for the other two selectively deuterated precursors when using C high data, which may be explained by higher precursor loadings that decrease the RO 2 lifetimes, thus hampering the HOM formation.In other words, comparing the C high 10 D 3 to the C low D 0 sample would only exaggerate the HOM yield decrease we expect to see from the deuteration.As the difference is still small enough, we can conclude that the decrease caused by selective deuteration is not significant for the 10 D 3 HOMs either, matching what was found for the other two selectively deuterated precursors.
Notably, the O 7 yield in Fig. 8a is smaller for the 3 D 1 sample.A difference can also be seen in the mass spectra (Fig. A2); i.e. the masses 308 (C 10 H 14 O 7 ) and 310 (C 10 H 16 O 7 ) Th have markedly lower signals in the spectrum for 3 D 1 compared to the D 0 spectrum.Figure A3 shows that, for C 10 H 14,16 O 7 , there is a clear difference in the fraction of D lost between 3 D 1 and 10 D 3 , but at C 10 H 14,16 O 9 , the difference is gone.Likewise, the yields of HOMs with eight or more O atoms (grey bars in Fig. 8a) are very close between the different labelled compounds.From a RO 2 radical perspective, it is clear that, to form C 10 H 15 O 10 , one necessarily has gone through C 10 H 15 O 8 , and as such it is surprising that the products of the less-oxidised radicals can decrease while the more-oxidised products do not.It is possible that different structures of C 10 H 15 O 8 can form, and the ones that can undergo further H shifts are not hampered by deuteration, while the isomers that cannot undergo H shifts and thus will terminate to e.g.C 10 H 14 O 7 are hampered.Based on the data we have, the latter option is supported; however, without data on the radicals, this finding remains highly speculative.

Conclusions
In this study, we tested the method of selective deuteration for assessing autoxidation mechanisms that are otherwise challenging to elucidate.We monitored the highly oxygenated oxidation products of three selectively deuterated αpinenes to calculate the fractions of HOM that had lost D atoms.We found that, in the cases where carbons C-3 and C-10 had been deuterated, D atoms were frequently lost, indicating that HOM formation happens through H abstractions from these C atoms.Contrastingly, no significant D-atom loss was observed in the case where C-7 had been deuterated, suggesting that this C atom was rarely involved in any isomerisation reactions during autoxidation and HOM formation.Our findings not only support the existence of the channel Iyer et al. (2021) postulated, but also indicate that there are additional, so far unknown oxidation channels for HOM formation from α-pinene ozonolysis.Furthermore, we found that the used deuterations did not significantly affect the HOM yields, contrasting with our expectations based on kinetic isotope effects.
Our experiments demonstrate that selective deuteration can be a very powerful tool for identifying autoxidation mechanisms.What is more, we can identify many improvements to these initial investigations.Firstly, the interpretation of the results would be simpler if we were able to observe the radicals themselves.Even better still would be to monitor less-oxygenated radicals, which has so far been shown to be possible only using few different instruments (e.g.Berndt et al., 2018).The radical signals are more straightforward to interpret because these compounds are at an earlier stage of the oxidation processes, removing different types of branching possibilities that can take place when the inspected closed-shell products are formed.Unfortunately, due to the large background and overlapping isotope signals that possibly obfuscated some of the signals at the smaller masses, we were not able to inspect the radical signals with the experimental design we used in this study.
Secondly, for a more complete picture of the autoxidation pathways, a more comprehensive set of selectively deuterated compounds would be beneficial.Lastly, it might be possible to investigate the most prominent signals with tandem mass spectrometry (MS 2 ) and gain further insight into their structures (Tomaz et al., 2021).
We have shown that using selectively deuterated precursors for studying the ozonolysis and oxidation of α-pinene is a viable method for mechanism development to gain insight into the autoxidation process still veiled in limited understanding of exact reaction pathways.We expect that utilising this approach and possibly other related isotope-based approaches will be the most promising route forward for understanding autoxidation propagation on a mechanistic level, both for α-pinene and other complex atmospheric VOCs of importance.As such, we are planning a continuation study with a more mechanistic orientation.The schematics of the α-pinene injection set-ups used in the experiments are shown in Fig. A1.The overflow injection set-up shown in Fig. A1a) consists of a gas bubbler and a small glass vial containing a very small amount (few tens of microlitres) of the precursor.The glass vial is covered with parafilm with a small hole in the film to control the rate at which the precursor is flushed out.A small flow of N 2 (25-150 mL min −1 ) is flushed over the glass vial through the gas bubbler to carry the evaporated precursor into the chamber.Due to the very small amounts of precursor samples we had, we wanted to use all the precursor sample from the vial.Thus, we also added the cap of the glass vial inside the bubbler in case any of the precursor was on it.This set-up was used when there was not enough precursor to take a sample in a syringe, and the resulting precursor concentration in the chamber was less stable than with the syringe pump injection set-up.
The syringe pump injection set-up shown in Fig. A1b was used when there was enough precursor to take a sample with a 5 µL syringe.We used a syringe pump to inject the precursor into a stream of 0.5 L min −1 N 2 flow carrying the precursor into the chamber.The precursor is injected into the stream using a small T piece, with one of its openings sealed with a septa plug through which the syringe is put and the two others connected to the N 2 flow.This set-up was much more controllable than the overflow set-up; however, there was a communication problem between the pump and the controlling software which caused the pump to inject at an order of magnitude faster rate within a certain range of set injection rates.In some cases, this resulted in far-too-high precursor concentrations in the chamber, increasing the large background and prohibiting us from using those data for analysis.

A2 HOM mass spectra
Figure A2 shows the mass spectra for all the precursors used, highlighting the closely inspected compounds and their isotopes used in Fig. 6

A3 HOM isotopic distributions for all found HOMs
Figure A3 shows the isotopic distributions for all found HOMs, including the compounds shown already in Fig. 6.Notably, unlike the signals for the compounds in Fig. 6, most of these signals have not been checked for overlapping signals, and hence we cannot say for sure whether there was no signal at all for some isotopes or that some of the signals were obfuscated by larger nearby signals.Additionally, since many of these signals are clearly smaller than those in Fig. 6, it is possible that only one of the isotopes' signals is above the detection limit of our instrument, which would lead to a 100 % contribution from that isotope, although others would have been almost equally high.Therefore, caution should be taken when interpreting the presented distributions.
We can see that C 9 and C 19 compounds have usually lost at least three deuterium atoms when using the 10 D 3 precursor.In other words, it is likely that the carbon lost in the formation of the C 9 and C 19 compounds is C-10, fragmenting with all the deuterium atoms connected to it.The behaviour differs significantly from that of C 10 and C 20 compounds, which have at most lost two deuterium atoms.
Additionally, when using the 3 D 1 precursor, the C 9 compounds' signals consist mostly of the isotopes that have lost one deuterium atom, and C 19 compounds mostly lost one https://doi.org/10.5194/acp-23-4373-2023Atmos.Chem.Phys., 23, 4373-4390, 2023  The colour shows how many deuterium atoms have been lost from the precursor during the formation of the HOM isotope that corresponds to the main HOM on the x axis.For example, in the case of the 3 D 1 precursor and the C10H14O7 compound, -0D corresponds to C10H13DO7 and -1D corresponds to C10H14O7.The entry for -4D is excluded from the upper legend as there were no corresponding compounds in the inspected monomers.deuterium and sometimes two deuterium atoms, meaning that the C-3 carbon is active in the formation process of these compounds.In contrast, these compounds' signals consist mostly of the isotopes that have lost zero deuterium atoms when using the 7 D 2 precursor, meaning that the C-7 carbon is rarely active.Figure A3 also contains the isotopic distributions for C 5 H 6 O 7 , for which the carbons C-3 and C-10 have in practice always lost all of the deuterium atoms and C-7 has never lost any.This likely means that, of the five carbons lost from α-pinene when forming this compound, two of them are the C-3 and C-10 carbons.
Furthermore, we can also see the general trends depicted in Fig. 6 in Fig. A3, where C-3 and C-10 are active and C-7 is inactive in the formation process.

A4 HOM yield behaviour as a function of the average number of D atoms lost
Figure A4 shows the comparison of HOM yields compound per compound between the D0 precursor and the deuterated precursors.It is otherwise the same as Fig. 8b-d, except that the colour shows the average number of deuterium atoms lost per compound.This modification could high-light the effect of the kinetic isotope effects better; however, despite the modification, we do not see a clear effect from the kinetic isotope effects.Nonetheless, Fig. A4 highlights the differing behaviours of the precursors, i.e. how 7 D 2 HOMs mostly did not lose deuterium atoms, whereas 3 D 1 and 10 D 3 did (the colour scales have a larger range for the latter two precursors).Additionally, we can see that, for the 3 D 1 precursor, most of the HOMs have lost on average some deuterium atoms (most points are in the green-coloured region), whereas for 10 D 3 precursor HOMs, there are clearly some HOMs that have lost more deuterium atoms than others (some points are yellow, while others are in the blue-to green-coloured region).

Figure 1 .
Figure 1.α-pinene with numbered carbons following IUPAC numbering (a) and the molecular structures of the selectively deuterated precursors used in this study (b-e).The precursors are named based on the number of the carbon that is deuterated and the number of deuterium atoms added to the molecule.The carbons with the deuterium atoms are highlighted with colours specific to the precursors.(b) Normal α-pinene is D 0 , (c) α-pinene with one deuterium atom is 3 D 1 (cyan), (d) α-pinene with two deuterium atoms is 7 D 2 (green), and (e) α-pinene with three deuterium atoms is 10 D 3 (blue).

Figure 2 .
Figure 2. (a) Reaction pathways for α-pinene ozonolysis up until the primary peroxy radical (RO 2 ) from Kurtén et al. (2015) (I-III) and Iyer et al. (2021) (IV).(b) Continued reaction pathway for the IV pathway up to a peroxy radical containing eight oxygen atoms (O 8 −RO 2 )(Iyer et al., 2021).In panels (a) and (b) the carbon and hydrogen atoms corresponding to deuterated ones in the selectively deuterated precursors are coloured matching the colouring in Fig.1; i.e. cyan corresponds to the 3 D 1 precursor, green corresponds to the 7 D 2 precursor, and blue corresponds to the 10 D 3 precursor.(c) Deuterium atoms are readily exchanged between HOMs and water vapour if the deuterium is in an −OD or −OOD group, resulting in HOMs without deuterium and semi-heavy water, i.e. a water molecule with one hydrogen atom and one deuterium atom instead of two hydrogen atoms(Rissanen et al., 2014).

Figure 3 .
Figure 3.The experimental set-up consists of a continuous stirred-tank reactor (CSTR).Additionally, there are ozone and zero-air generators, an α-pinene injection set-up, an ozone monitor, a nitrate CI-Orbitrap, and a PTR or a Vocus.

Figure 4 .
Figure 4. Time series of the experiments showing total HOM concentration (green solid line), precursor concentration (black dashed line), α-pinene ozonolysis rate multiplied by 10 (blue dotted line), and ozone concentration multiplied by 0.1 (red dash-dotted line).Additionally, sample times for a low-concentration sample (C low ) and a high-concentration sample (C high ) are shown with a magenta "\" hatched area and an olive cross-hatched area, respectively.

Figure 5 .
Figure 5. Example spectra using C 10 H 14 O 11 to show resulting HOM signals when using different selectively deuterated precursors.The nitrate cluster of C 10 H 14 O 11 (C 10 H 14 O 11 • NO −3 ) is at 372 Th.The colours show which selectively deuterated α-pinene is the precursor, and the values above the bars in matching colour tell how many deuterium atoms the HOMs have lost compared to the precursor.For example, the α-pinene with three deuterium atoms ( 10 D 3 , shown in blue) can form up to four different C 10 H 14 O 11 isotopes in the spectrum.However, in the 10 D 3 precursor's case there is no signal for C 10 H 14 O 11 without any deuterium atoms (−3D).

Figure 6 .
Figure 6.Isotopic distributions of chosen HOMs as proportions of signals using high-concentration C high data.The colours show how many deuterium atoms are lost from the molecule during the formation process.For example, in the case of 3 D 1 precursor and C 10 H 14 O 7 compounds, −0D corresponds to C 10 H 13 DO 7 and −1D corresponds to C 10 H 14 O 7 .The entries for −3D and −4D are excluded from the upper legend as there were no corresponding compounds in the inspected monomers.

Figure 7 .
Figure 7. Molecular structure of α-pinene with numbered carbonsshowing carbons that are active (carbons 3 and 10; green eight-point star) and inactive (carbon 7; light-blue diamond) in the formation process, carbons that were not inspected in this study (carbons 1, 4, 5, 8, and 9; magenta circle), and carbons that cannot be studied using selective deuteration (carbons 2 and 6; grey hexagon).

Figure 8 .
Figure 8.(a) HOM yields for α-pinene samples, where yields for C 5 H 6 O 7 (mint green), O 7 monomers (orange), O 8+ monomers (grey), and dimers (blue) are plotted separately.The corresponding precursor names and concentrations are shown on the x axis.(b-d) Comparison of HOM yields compound per compound between the D 0 sample and 3 D 1 sample (b), 7 D 2 sample (c), and 10 D 3 sample (d), where the mint green cross corresponds to C 5 H 6 O 7 and is excluded from the data trend and the colour shows the number of oxygen atoms in the compound.C low and C high refer to the α-pinene concentration of the sample.Note that each point in panels (b)-(d) corresponds to the yield of the sum of all found isotopes per main compound.
Figure A1.α-pinene injection set-ups used in the experiments.(a) The overflow injection set-up was a gas bubbler into which the small glass vial holding the precursor sample (and its cap) were inserted.(b) The syringe pump injection set-up consists of a syringe pump, a syringe, a metal T piece, and a septa plug.

Figure A2 .
Figure A2.Mass spectra of HOMs from selectively deuterated α-pinene precursors at unit mass resolution.Masses corresponding to HOM signals used on closer inspection are marked with a colour, whereas other peaks in the mass spectra are shown in black.(a-i) and (a-ii): the closer-inspection compounds are 1.C 10 H 14 O 7 , 2. C 10 H 14 O 9 , 3. C 10 H 16 O 9 , 4. C 10 H 14 O 11 , 5. C 19 H 28 O 11 , 6. C 20 H 30 O 14 , 7. C 20 H 30 O 16 , and 8. C 20 H 30 O 18 .The y scales are different between the monomer and dimer plots.

Figure A3 .
Figure A3.Isotopic distributions of all inspected peaks.The colour shows how many deuterium atoms have been lost from the precursor during the formation of the HOM isotope that corresponds to the main HOM on the x axis.For example, in the case of the 3 D 1 precursor and the C10H14O7 compound, -0D corresponds to C10H13DO7 and -1D corresponds to C10H14O7.The entry for -4D is excluded from the upper legend as there were no corresponding compounds in the inspected monomers.

Figure A4 .
Figure A4.Comparison of HOM yields compound per compound with the number of deuterium atoms lost on average shown with the colour between the D 0 sample and the 3 D 1 sample (a), the 7 D 2 sample (b), and the 10 D 3 sample (c), where the mint green cross corresponds to C 5 H 6 O 7 and is excluded from the data trend.C low and C high refer to the α-pinene concentration of the sample.Note that each point corresponds to the yield of the sum of all found isotopes per main compound.
signals in α-pinene ozonolysis HOM mass spectra: four HOM monomers C 10 H 14 O 7 , C 10 H 14 O 9 , C 10 H 16 O 9 , and C 10 H 14 O 11 and four HOM dimers C 19 H 28 O 11 , C 20 H 30 O 14 , C 20 H 30 O 16 , and C 20 H 30 O 18 5194/acp-23-4373-2023 Atmos.Chem.Phys., 23, 4373-4390, 2023stractable H/D atoms need to have been abstracted for the most-oxygenated HOMs to form.For C-10, the monomers lost at most one D atom, while the C 20 dimers lost up to two D atoms, indicating that it is rare to have more than one H/D abstraction from the same C atom.The C 19 dimer is discussed later.