Oxidation product characterization from ozonolysis of the diterpene ent -kaurene

. Diterpenes (C 20 H 32 ) are biogenically emitted volatile compounds that only recently have been observed in ambient air. They are expected to be highly reactive, and their oxidation is likely to form condens-able vapors. However, until now, no studies have investigated gas-phase diterpene oxidation. In this paper, we explored the ozonolysis of a diterpene, ent -kaurene, in a simulation chamber. Using state-of-the-art mass spectrometry, we characterized diterpene oxidation products for the ﬁrst time, and we identiﬁed several products with varying oxidation levels, including highly oxygenated organic molecules (HOM), monomers, and dimers. The most abundant monomers measured using a nitrate chemical ionization mass spectrometer were C 19 H 28 O 8 and C 20 H 30 O 5 , and the most abundant dimers were C 38 H 60 O 6 and C 39 H 62 O 6 . The exact molar yield of HOM from kaurene ozonolysis was hard to quantify due to uncertainties in both the kaurene and HOM concentrations, but our best estimate was a few percent, which is similar to values reported earlier for many monoterpenes. We also monitored the decrease in the gas-phase oxidation products in response to an increased condensation sink in the chamber to deduce their afﬁnity to condense. The oxygen content was a critical parameter affecting the volatility of products, with four to ﬁve O atoms needed for the main monomeric species to condense onto 80 nm particles. Finally, we report on the observed fragmentation and clustering patterns of kaurene in a Vocus proton-transfer-reaction time-of-ﬂight mass spectrometer. Our ﬁndings highlight similarities and differences between diterpenes and smaller terpenes during their atmospheric oxidation, but more studies on different diterpenes are needed for a broader view of their role in atmospheric chemistry.

BVOCs oxidation products can contribute to atmospheric aerosol loadings if their vapor pressures are low enough. Over the past decades, extensive studies have been conducted to investigate the oxidation mechanism and the role of terpenes in secondary organic aerosol (SOA) formation around the world (Mohr et al., 2017;Zhang et al., 2018;Li et al., 2021;Kontkanen et al., 2016;Ehn et al., 2012;Yassaa et al., 2012), in particular from isoprene and monoterpenes. Isoprene has a relatively low yield of SOA formation, while monoterpenes play a more vital role in the SOA formation with typical SOA yields of ∼ 5 %-10 % (Mutzel et al., 2016;Griffin et al., 1999;Mcfiggans et al., 2019;Zhang et al., 2018;Ehn et al., 2014). Also, sesquiterpenes have recently received growing interest and have been reported to be even more efficient SOA contributors, with reported SOA yields from 6 % up to 100 % (Jaoui et al., 2013;Lee et al., 2006;Boy et al., 2007;Li et al., 2011). However, with the yearly emissions of isoprene, monoterpenes, and sesquiterpenes being around 500, 160, and 30 Tg yr −1 , respectively (Guenther et al., 2012), the relative contributions of these different compound groups to SOA are not easy to determine. This is even truer for other terpene groups, for which very few studies exist to date.
Diterpenes have previously been identified in some extractions of trees such as leaf oils, pine needles, and resins (Kato, 2005;Keeling and Bohlmann, 2006;Lee et al., 2009;Tumen et al., 2010); however, they were for a long time thought to be not emitted into the atmosphere due to their fairly low volatility (Guenther, 2002). Matsunaga et al. (2012) were the first to report diterpene kaurene emissions from branch enclosure experiments of coniferous trees in Japan. After that, several studies reported the observation of diterpenes using similar branch enclosure systems with large variation in the estimated emission rates in different studies Lin et al., 2015;Yáñez-Serrano et al., 2018;Helin et al., 2020). For instance, the estimated emission rates of diterpenes were 1-3 orders of magnitude lower than that of monoterpenes and sesquiterpenes in the boreal forests (Helin et al., 2020). However, in a study of Mediterranean shrubs and a study in Japan during warm seasons (temperature ∼ 30 • C), the emission rates for diterpenes were comparable to those of monoterpenes (Matsunaga et al., 2012;Yáñez-Serrano et al., 2018). Additionally, diterpene emission rates could be further increased under stress conditions like heat and mechanical wounding (Lin et al., 2015).
Very recently, direct ambient observations of diterpenes have become feasible with the latest developments of mass spectrometric techniques. To our knowledge, three studies have reported diterpenes in the atmosphere, but diterpene concentrations were much lower compared to monoterpenes and sesquiterpenes (Chan et al., 2016;Li et al., 2020;. For example, Li et al. (2020) directly measured around 1.7 ppt (parts per trillion) of diterpenes in the French Landes forest during summertime, which is 10-1000 times lower than monoterpenes and sesquiterpenes. Moreover, four diterpenes (the most abundant one was kaurene) with concentrations ranging from 10 to 86 ppq (parts per quadrillion) were detected at a rural site in the Amazon region by Yee et al. (2018), while most detected sesquiterpenes were above 100 ppq. The low concentrations are likely primarily linked to the lower emission rates of the less volatile diterpenes compared to smaller terpenes, and their reactivity is also likely to be high once they are in the atmosphere. Although no diterpene oxidation products have been reported, it can be assumed that oxidation is a main sink of diterpenes, with unknown impacts on atmospheric chemistry and SOA formation.
In this study, we characterize the oxidation products formed via ozonolysis of the diterpene ent-kaurene (hereafter referred to as kaurene). The precursor and oxidation products in the gas phase were measured using state-of-the-art mass spectrometers (Jokinen et al., 2012;Krechmer et al., 2018). We present the first comprehensive summary of the identified oxidation products with different oxidation levels, including highly oxygenated organic molecules (HOM; Bianchi et al., 2019). We also probe the potential of the oxidation products to condense onto existing particles to form SOA. Finally, we discuss fragmentation and clustering patterns of kaurene using the Vocus proton-transfer-reaction time-of-flight (PTR-TOF) mass spectrometer.

Chamber setup
Kaurene ozonolysis experiments were conducted in a 2 m 3 Teflon (FEP, supplied by Vector Foiltec, Germany) chamber, the COALA 1 chamber, at the University of Helsinki, Finland. Further details of the chamber facility can be found elsewhere Peräkylä et al., 2020). During the experiments, the chamber was operated under steadystate conditions with a continuous inflow of ∼ 36 L min −1 , and the residence time was around 56 min. The injected flow consisted of purified air generated by a clean air system (AADCO, series 737-14, Ohio, USA), with additions of varying amounts of O 3 (Dasibi 1008-PC ozone generator) and kaurene, as well as 80 nm ammonium sulfate (AS) seed particles. Instruments sampled the majority of the outflow for chemical composition measurement both in the gas and the particle phases (Sect. 2.2), while the rest was flushed into an exhaust line. O 3 concentration in the chamber was monitored by a UV photometric analyzer (model 49P, Thermo-Environmental) during the experiments. Temperature, relative humidity (RH), and pressure in the chamber were monitored using a Vaisala temperature and humidity probe (INTERCAP ® HMP60) and a differential pressure sensor (Sensirion SDP1000-L025). All experiments were carried out at room temperature (26 ± 1 • C) under a slight overpressure condition to minimize leaks into the chamber and with the RH consistently below 1 %.
All instruments sampled the chamber air continuously, and the input into the chamber was varied (Fig. 1a) to achieve different concentrations and ratios of kaurene and O 3 . The O 3 concentrations ranged from ∼ 70 to ∼ 110 ppb for different experiments, and the kaurene concentration was at most ∼ 20 ppb, though most of the time well below 1 ppb (Fig. 1a). Seed aerosols were added at different times (the gray shading in Fig. 1) in order to initiate condensation onto aerosols rather than chamber walls, which is the dominant sink for low-volatility vapors. The total particle mass reached its maximum value of more than 200 µg m −3 in the experiment on the 31 January. The injection of kaurene was one of the most challenging parts of the experiments due to its low volatility. Kaurene is a solid at room temperature, and the kaurene in this study was purchased from OIChemIm s.r.o. Gas-phase kaurene was introduced into the chamber by flushing nitrogen (N 2 ) through a vial containing the solid kaurene, and the flow was subsequently directed into the chamber. In addition, a heater was placed under the vial, heating the air around it to promote the evaporation of kaurene. During the experiments, the vial's bottom temperature stayed below 45 • C, except for the last experiment day when the vial reached a temperature around 60 • C. Already at low heating, the solid kaurene had melted to liquid form, but we did 1 Comprehensive molecular characterization of secondary Organic AerosoL formation in the Atmosphere not see indications of decomposition of the kaurene from the heating in any of our mass spectra. Variable levels of kaurene in the chamber were achieved, either by turning the heater on and off or changing the flow rate of the N 2 carrier gas, which was controlled with a mass flow controller (MKS, G series, Andover, MA, USA). Due to the low volatility of kaurene, residual amounts were introduced into the chamber also when no active addition was performed, presumably due to evaporation from surfaces in the tubing and chamber.

Instrumentation
We deployed a suite of online instrumentation measuring gas-and particle-phase species. The Vocus PTR-TOF (Tofwerk AG/Aerodyne Research, Inc.) (Krechmer et al., 2018) was applied to determine the reactant kaurene concentration and some of its oxidation products. The instrument allows for fast and continuous measurements of VOCs with subppt detection limits (Krechmer et al., 2018). Based on the previous evaluation, the Vocus can also detect various oxidation products of monoterpenes containing up to around six oxygen atoms . In this study, the instrument was run with the axial and radial voltages being 350 and 400 V, respectively. The pressure in the ionization regions (drift tube) was set to 1.4 mbar, giving an electric field strength / number density of gas E/N of 120 Td (townsend). However, as discussed in more detail in Sect. 2.3, the pressure drifted during the experiments, resulting in significant changes in the absolute sensitivity of the sampled molecules. Around 4.5 L min −1 of the chamber air was directed into the Vocus through 1 m long PTFE tubing (6 mm o.d., 4 mm i.d.) for further analysis. A water flow (HPLC grade or 18 M Milli-Q water) at the rate of 15 sccm (standard cubic centimeters per minute) was utilized to produce reagent ions in a discharged ion source. The mass-resolving power of the long TOF mass analyzer was 12 000-13 000, and data were recorded with a time resolution of 5 s.
For measuring HOM and other more oxygenated products from kaurene ozonolysis, a chemical-ionization atmospheric-pressure-interface time-of-flight mass spectrometer was used (CI-APi-TOF, Tofwerk AG/Aerodyne Research, Inc.) (Jokinen et al., 2012). We used nitric acid (HNO 3 ) as the reagent and an X-ray source to produce nitrate reagent ions (NO − 3 ). The sample molecules could be charged by collision with nitrate ion clusters (NO − 3 , HNO 3 NO − 3 , and (HNO 3 ) 2 NO − 3 ). After collisions, the sample molecules were either ionized via a direct proton transfer and detected as a deprotonated ion or, more typically, via clustering with reagent ions and detected as adducts (Jokinen et al., 2012). During the experiments, the instrument, equipped with a long TOF mass analyzer (mass-resolving power 13 000-14 000), was configured to measure ions up to 960 Th with a time resolution of 10 s. The sampling flow rate was around 10 L min −1 for all experiments.  . HOM monomers and dimers of kaurene ozonolysis detected by nitrate CI-APi-TOF are shown in panels (c)-(e). The light gray shading indicates periods when there were seed particles in the chamber. The black dashed lines are included in each panel to depict the kaurene ozonolysis rate ("Kaurene*O 3 " in the right y axis represents kaurene concentration times O 3 concentration, in units ppb 2 ). The labels "S1 and "S2" mark two steady-state periods used to determine the condensation behavior of oxidation products (Sect. 3.3).
A custom-built scanning mobility particle sizer (SMPS) and a long time-of-flight aerosol mass spectrometer (L-TOF-AMS, hereafter AMS, Aerodyne Research Inc.) were deployed to measure the particle phase. The L-TOF-AMS is similar to the high-resolution TOF-AMS (HR-TOF-AMS) described in DeCarlo et al. (2006) but has a longer TOF chamber for increased mass-resolving power.
Data of both the Vocus and the nitrate CI-APi-TOF were analyzed with a MATLAB tofTools package (version 607) (Junninen, 2014), while the TOF-AMS analysis software packages SQUIRREL (version 1.63H) and PIKA (version 1.23H)  within Igor Pro (version 6.37, WaveMetrics Inc.) were used for the AMS data. The SMPS data were analyzed with MATLAB and Igor Pro.
A system for total O 3 reactivity measurement (TORM) was also deployed during the experiments. The TORM consists of a reactor (three 2 L borosilicate glass bottles) in which sampled air reacts with a set concentration of O 3 (150 ppb in this study) produced by an O 3 generator. A modified O 3 analyzer (model 49i, Thermo Scientific, Waltham, MA, USA) operating in differential mode recorded directly the difference in O 3 before and after the reactor, which is used to calculate the total O 3 reactivity in the COALA chamber. The kaurene concentration could be estimated based on the calculated O 3 reactivity data at a given rate coefficient of kaurene-O 3 reaction, with the assumption that kaurene was the dominant sink for O 3 . A detailed explanation of the instrument and calculation can be found elsewhere (Helmig et al., 2021). During the experiments, the total flow through the reactor was 6.0 L min −1 with 5.0 L min −1 sampling flow from the COALA chamber and 1.0 L min −1 coming from the O 3 generator (dilution factor f = 0.83). The residence time in the reactor was estimated to be 60 s.
Multiphase adsorbent tubes (Tenax TA/Carbopack B, 60-80 mesh) (o.d. 1/4 in. × 3 1/2 in.; PerkinElmer Inc., Waltham, MA, USA) were utilized as one way to quantify diterpene concentration in the chamber due to the lack of any authentic diterpene standards for Vocus. We collected four Tenax tubes during the experiments, and details about the sampling parameters can be found in Table S1 in the Supplement. After sampling, a thermal desorption gas chromatograph (TD-GC, TurboMatrix 350 automatic TD unit, Clarus 680 GC, PerkinElmer Inc.) coupled with a quadrupole mass spectrometer (MS, Clarus SQ 8 T, PerkinElmer Inc.) was applied to determine and quantify kaurene in the Tenax samples. More details about the methodology and Tenax tubes have been reported by Helin et al. (2020).

Kaurene quantification and correction of Vocus time series
In Vocus, the sampled air is charged in a focusing ionmolecule reactor (FIMR) by combining with the reagent ions produced from a low-pressure discharge reagent ion source.
The FIMR is a glass tube with four quadrupole rods installed outside in a radial pattern. The FIMR pressure is controlled by a valve between the reactor and a mechanical pump. The FIMR pressure typically remains constant due to the control where additional pumping is applied to it to keep a set pressure. However, in our study, the pinhole through which the sample air entered the FIMR became partially blocked by seed aerosol particles. When sufficiently blocked, the pressure control could no longer be maintained since even without additional FIMR pumping, the pressure dropped below the desired 1.4 mbar, as the flow through the pinhole had decreased. As the pinhole clogging was not noticed during the experiments, the FIMR pressures kept decreasing in a stepwise pattern from the initial ∼ 1.4 mbar at the beginning of the experiments to ∼ 0.9 mbar in the end. The changes in FIMR pressure are expected to impact both the sensitivity and the fragmentation behavior of the Vocus by affecting the frequency and energy of ion-molecule collisions. Hence, the raw signal intensity of kaurene could not directly depict the actual changes of kaurene in the chamber over the course of the experiments. Therefore, we needed to correct the sensitivity changes as far as possible to simulate the actual variation of kaurene using a simple model. For this purpose, we utilized the observed changes in chamber O 3 during kaurene additions. The details are described in the following.
The O 3 concentration was monitored online for all experiments. When we injected kaurene into an O 3 -stable chamber, the O 3 loss rate increased due to reactions with kaurene (in addition to ever-present flush-out). For any time point, we can write the change in O 3 concentration in the chamber as where Q O 3, in is the injection rate of O 3 ; [kaurene] and [O 3 ] are the current kaurene and O 3 concentrations, respectively; k flush out is the rate at which air is flushed out of the chamber; and k 1 is the kaurene-O 3 reaction rate coefficient. For every experiment starting with a stable O 3 condition, we simulated a kaurene time series. We solve Eq. (1) for [kaurene], and take it as the expected kaurene concentration at time i (kau_exp i ). The rate d [O 3 ] dt was estimated based on the average change of O 3 from time point (i − 1) to the time point (i + 1) (time resolution: 15 min = 900 s), giving the expression finally as The average residence time in the COALA chamber concerning flush out was 56 min; thus, the reciprocal k flush out = 3 × 10 −4 s −1 . Q O 3, in (unit cm −3 s −1 ) was estimated for each experiment from the steady-state periods before kaurene injection, when d[O 3 ] dt = 0. The kaurene-O 3 reaction rate coefficient k 1 was set to be 5 × 10 −16 cm 3 s −1 based on the model results of the rate coefficient tests of kaurene and O 3 performed in the beginning of our experiments (see the Supplement). This value is fairly close to the only value (1.4 ± 8.2 ×10 −15 cm 3 s −1 ) so far reported in the literature (Helin et al., 2020).
The simulated kaurene concentrations were only reliable when kaurene concentrations were high enough (> ∼ 1 ppb) to make a noticeable impact on the O 3 concentrations. We used these periods to scale the measured Vocus raw kaurene signal (kau_sig i ) to the simulated values (kau_exp i ) by a correction factor, cor_fac i : We further assumed that this factor stayed constant as long as the pressure in the FIMR did not change; thus, we also extrapolated the correction factors to periods with low kaurene concentrations. In this way, a complete and consecutive time series of correction factors for kaurene was achieved (see the blue line in Fig. 2a).
All the above corrections bring with them considerable uncertainty and potential error sources. In order to assess the reliability of the correction method, we compared our acquired kaurene time series with those measured by TORM and the four collected Tenax samples. First, quantified kaurene concentrations by Tenax tubes are shown in Table S1. The average Vocus signal intensity of kaurene during the sampling time for each Tenax tube was calculated, and the corresponding correction factors required by the Vocus to match the Tenax concentrations can be found in Fig. 2a. The correction factors we derived above were in the same order of magnitude as those deduced from the Tenax samples, although the Tenax samples consistently suggested lower kaurene concentrations. However, the third Tenax sample (Table S1) on 28 January was 2 orders of magnitude lower, which leads us to expect that this sample was unsuccessful for unknown reasons. For comparison with TORM, we scaled the raw Vocus signal intensity data of kaurene to match the estimated kaurene based on TORM. The required scaling factor time series is shown in Fig. 2a for periods when the kaurene was above the detection limit (∼ 1 ppb) of TORM. Overall, the TORM calibration factors compared well with those from the O 3 depletion method above. The good agreement becomes evident when plotting the final kaurene time series from TORM and the corrected Vocus data (Fig. 2b), lending confidence to the Vocus correction method. The figure also shows the importance of utilizing the Vocus data, as the TORM was not able to quantify kaurene concentrations below ∼ 1 ppb due to background. Despite this, the correlation between the instruments was quite good (R 2 = 0.87, Fig. 2c). Based on the comparison of our simulated kaurene results with these two independent methods, we conclude that our correction approach was successful. However, the uncertainties remain large due to the needed corrections, the lack of authentic standards, and the dependence on needing to know the kaurene reaction rate coefficient with O 3 . We estimate that the absolute concentrations have an uncertainty of at least a factor of 3 (+200 %/ − 67 %), while the uncertainty in the relative changes over the course of our study is closer to a factor of 2 (+100 %/ − 50 %).

HOM formation
Nearly all atmospheric VOCs oxidation will produce peroxyl radical (RO 2 ) intermediates, and the final distribution of oxidation products is largely determined by the reaction pathways available to these RO 2 species (Atkinson and Arey, 2003). Autoxidation is one important propagation process for RO 2 recently shown to be of significance in the atmosphere (Crounse et al., 2013;Ehn et al., 2014;Bianchi et al., 2019). In the autoxidation process, intramolecular H-shift reactions can convert an RO 2 into a carbon-centered radical with a hydroperoxide functionality. A new RO 2 can then be formed through the subsequent addition of oxygen. In some precursor molecules, this process can repeat multiple times, and the rapid addition of oxygen in this process can eventually result in very high oxygen contents of the produced RO 2 .
The final formation of a closed-shell molecule from RO 2 takes place through different termination reactions. For instance, RO 2 can decompose in a unimolecular way, e.g., by OH loss to form a carbonyl product (Rissanen et al., 2014). In addition, a number of bimolecular reactions can efficiently terminate RO 2 and form a variety of organic species. The main terminating reaction partners in the atmosphere are NO, HO 2 , and other peroxyl radicals (Vereecken and Fran-cisco, 2012). In the specific reaction of an RO 2 with another R O 2 , ROOR dimers can be formed (Berndt et al., 2018;Finlayson-Pitts and Pitts Jr, 1999).
Kaurene has an exocyclic double bond, as shown in Fig. 3, and can react with both O 3 and OH. In reaction with O 3 , two separate molecules will form following the scission at the initial double bond. Out of the resulting C 1 and C 19 species, one will typically form a "primary" RO 2 and the other a closedshell carbonyl. The initial steps can be considered equivalent to the ozonolysis of the monoterpene β-pinene, which also has a lone double bond attached to a ring structure (Zhang and Zhang, 2005). If the carbonyl is formed on the C 1 fragment, we expected to produce a C 19 H 29 O 3 peroxyl radical (Fig. 3). In addition to the oxidation by O 3 , also OH is of importance, as it is formed during the experiments as a byproduct of kaurene-O 3 reactions. Kaurene oxidation by OH initiate via either addition of OH to the double bond or abstraction of an H atom. Through OH attachment, the primary RO 2 will be C 20 H 33 O 3 (Fig. 3), and H abstraction by OH will result in C 20 H 31 O 2 as the primary RO 2 (Fig. 3). These three radicals, initiated by ozonolysis and OH oxidation, are expected to be the starting points for the majority of observed large oxidation products, though no mechanistic studies exist on the radical reactions in kaurene oxidation. In this study, the potential of the radicals to undergo autoxidation and achieve higher oxygen contents was of particular interest. Oxidation products with five or more oxygen atoms are classified as HOM in this work.

HOM quantification
Based on earlier observations, HOM formed from kaurene oxidation are expected to be efficiently detected by the ni- trate CI-APi-TOF (Bianchi et al., 2019). However, the approaches for valid and suitable calibrations still remain limited for these highly oxygenated and reactive molecules, severely hampering their reliable quantification. For estimation of HOM concentrations in this study, we adopted a similar approach as previously reported (Jokinen et al., 2012(Jokinen et al., , 2014Ehn et al., 2014) to convert the measured HOM ion signals to HOM concentrations using the following equation: Here, [HOM] is the estimated concentration of HOM in the chamber, C is a calibration factor determining the sensitivity, and HOM i · NO − 3 is the signal intensity of a molecule classified as HOM and measured by the nitrate CI-APi-TOF as a cluster with NO − 3 . This approach implicitly assumes an identical sensitivity to all HOM. The ionization of the majority of the HOM is expected to proceed at the collision limit (Bianchi et al., 2019;Ehn et al., 2014;Hyttinen et al., 2015), but factors like mass-dependent transmission Heinritzi et al., 2016) of the instrument may cause additional uncertainties. Still, the main challenge is to determine the calibration factor C, accounting for both charging efficiency and HOM losses in the sampling lines. In this study, we chose C = 1 × 10 10 cm −3 based on typical values reported earlier (Jokinen et al., 2012(Jokinen et al., , 2014. This approach comes with significant uncertainties, estimated to be at least a factor of 2 (+100 %/ − 50 %). A strongly mass-dependent ion transmission would be the most likely reason for a clearly larger deviation. Nevertheless, given the large uncertainty also in the kaurene concentration, we only use the calibration factors to approximate the HOM concentrations and yields.

Molar yield calculation
The change rate of HOM concentration in our chamber was determined by the production and loss rates according to Eq. (5): HOM can be formed from reactions with both O 3 and OH. Most of the time, the kaurene concentration was very low in the chamber and may not have been the major sink for the formed OH. Instead, OH was likely lost to reactions with O 3 and potential contaminant molecules in the chamber. In addition, it is also hard to estimate the relative importance of the two oxidants, because neither the O 3 nor OH reaction rate coefficients with kaurene are well known. Thus, for simplicity, the production of HOM in the equation above will be written as purely the ozonolysis reaction, i.e., k 1 γ [kaurene][O 3 ], where k 1 is the kaurene-O 3 reaction rate co-5626 Y. Luo et al.: Oxidation product characterization from ozonolysis of the diterpene ent-kaurene efficient, and γ is the HOM molar yield from the kaurene-O 3 reaction. Note that γ will be overestimated if the kaurene-OH reactions contribute substantially to the measured HOM. The kaurene-O 3 reaction rate coefficient k 1 was assumed to be 5 × 10 −16 cm 3 s −1 as shown in Sect. 2.3. The loss of HOM mainly consists of three parts: chemical reactions, flush out from the chamber, and the condensation onto the walls and potential aerosol particles; they can be written as k loss [HOM], where k loss is the total loss coefficient of HOM. Under steady-state conditions, the concentration of HOM is constant, and we can write it as d [HOM] dt Thus, All data used to estimate HOM molar yield were selected from the periods when there were no seed particles in the chamber. According to an earlier study in the COALA chamber by Peräkylä et al. (2020), the lifetime of HOM in the gas phase under a typical situation without seed particles could be assumed to be 200 s, primarily controlled by losses to chamber walls. Therefore, 0.005 s −1 was used for k loss here. The flush-out rate is over an order of magnitude slower, and considering the low oxidant concentrations in our chamber, it is unlikely that the chemical loss of HOM would be faster than the wall loss rate. The HOM quantification alone had a large uncertainty, which we can only estimate to be at least +100 %/ − 50 %. Also the kaurene quantification required several assumptions, as did k 1 . As for HOM quantification, we can only make estimations on the uncertainties of these parameters. If we assume that these three are the major error sources, and each of them have the same uncertainty of +100 %/−50 %, the final uncertainty (via error propagation) becomes +173 %/ − 87 %.

Results
The ozonolysis of kaurene in the gas phase was investigated in the COALA chamber under different oxidation conditions. In the following, the characteristics of oxidation products from kaurene ozonolysis, including chemical identification, are first described (Sect. 3.1), followed by HOM molar yield estimates (Sect. 3.2). Section 3.3 discusses the volatilities of the oxidation products, and Sect. 3.4 depicts fragmentation patterns of kaurene and the formation of NH + 4 clusters in the Vocus.

HOM and other oxidation products of kaurene ozonolysis
We observed the oxidation products of kaurene ozonolysis in the gas phase using Vocus and nitrate CI-APi-TOF. Prod-uct distributions are shown in Fig. 4 including both less oxidized species and HOM, and the temporal behavior of some important oxidation products are shown in Fig. 1. When discussing detected compounds (M) in the following sections, they were either detected after protonation (MH + ) in the Vocus or as clusters with nitrate ions (M·NO − 3 ) in the CI-APi-TOF. In other words, the identified species are presented after omitting the reagent ions, but if the specific ions need to be addressed, they will contain the charging ion as well as the label for the charge ( + or − ). When masses are discussed, they will refer to the mass-to-charge ratio in the spectra where the compound was identified, including the corresponding reagent ion.
We first compared the spectra of the reactant-free chamber (no kaurene or O 3 injection) with that after injecting ∼ 5 ppb kaurene. As shown in Fig. 4b, a predominant peak of kaurene, C 20 H 32 , along with a series of hydrocarbon ions was detected by the Vocus. These hydrocarbons with carbon numbers from 10 to 16 were observed between 130 and 230 Th. The time series of those hydrocarbons showed high similarities to the kaurene time series with very high correlation coefficients (R 2 >0.95) at different FIMR pressures, indicating that they most likely are fragments of kaurene formed in the Vocus. More details about the fragments of kaurene in the Vocus will be discussed in Sect. 3.4.
After introducing O 3 (∼ 70 ppb), several less oxidized kaurene oxidation products (C 18−20 H 26−34 O 1−3 , complete peak list in Table S3.) were observed between 270 and 330 Th in the Vocus spectra (Fig. 4b). The highest concentrations were observed from compounds identified as C 19 H 30 O, C 19 H 32 O 2 , and C 20 H 32 O 3 , of which the time series is shown in Fig. 1. Compounds with identical elemental compositions have been observed in the headspace samples of heated pine needles and spruce twigs (Helin et al., 2020), but to the best of our knowledge, these compounds have not been reported in the ambient air. The ions containing an N atom are believed to be adducts with NH + 4 , as will be discussed in Sect. 3.4.
In the nitrate CI-APi-TOF, various HOM monomers (C 18−20 H 24−34 O 5−13 ) and dimers (C 37−40 H 58−68 O 5−12 ) were detected in the mass-to-charge ratio (m/z) ranges 400-550 and 650-750 Th, respectively ( Fig. 4c and d). Several radicals with C 19 and C 20 carbon skeletons were observed as well (Table S4). Overall, the spectra contain a wide range of molecules with different compositions, and some speculation about possible formation pathways is given below. However, we stress that our data are not optimal for deducing detailed oxidation mechanisms, and in some cases, several pathways may lead to the same product, while for some other products, no viable formation pathways were identified. Among the closed-shell HOM monomers, C 19   . Background-subtracted (diff) spectra of kaurene and kaurene oxidation products. The spectrum in (a) shows the difference between the Vocus spectra before and after kaurene injection; the spectrum in (b) displays the change of Vocus spectrum before and after O 3 was injected into the chamber with kaurene. Comparison of nitrate CI-APi-TOF spectra before and after O 3 injection into the kaurene-existing chamber is shown in (c) and (d). The compounds in (a) and (b) are labeled as what they were actually detected by Vocus, whereas all peaks labeled in (c) and (d) are detected as a cluster with NO − 3 . molecular way (for the former) or via reaction with HO 2 (for the latter). Corresponding highly oxidized RO 2 species were detected (C 19 H 29 O 7,9 ,  The signal intensities of HOM dimers were about 1 to 2 orders of magnitude lower than those of monomers (Fig. 4d) species would most easily be explained by the involvement of these radicals. Another unanticipated finding was that we did not detect any dimers with oxygen numbers above 12, and the most abundant dimer signals were identified as compounds with 6-8 oxygen atoms. This is similar to the detected monomers, which is in contrast to earlier findings of monoterpenes and 5628 Y. Luo et al.: Oxidation product characterization from ozonolysis of the diterpene ent-kaurene sesquiterpene (Quéléver et al., 2019;Jokinen et al., 2016;Li et al., 2021;Barreira et al., 2021;Kirkby et al., 2016;Rissanen et al., 2014), where dimers typically contained clearly more oxygen than the monomers. We can only speculate on the reasons for this, and a possible explanation is that the most oxygenated RO 2 species have very short lifetimes and thus terminate through unimolecular channels before colliding with another RO 2 . It is also possible that the nitrate CI-APi-TOF might be more sensitive towards less oxygenated dimers if there are 40 C atoms rather than 20 in their skeleton, if a larger molecule increases the likelihood of forming a stable cluster with NO − 3 . Mass-dependent transmission may also play a role, as the detected dimers were around 700 Th and above, where sensitivity is most likely decreasing .
Interestingly, a group of HOM with 21-29 carbon atoms was observed in the m/z range from 500 to 560 Th (Fig. 4c), and they exhibited similar temporal behaviors to those of other HOM (Fig. 1). The dominant peaks were identified as C 24,26,27,28 H 36,38,40,42 O 6−8 species. We conclude that the most probable explanation for these HOM is that they were the dimers formed via the reactions of kaurene-derived RO 2 with smaller RO 2 derived from various C 2−9 organic contaminants in the chamber. Both Vocus and nitrate CI-APi-TOF observed such contaminants during the experiments (Fig. S2), which is the likely source being fairly high-loading monoterpene experiments performed in the chamber prior to the kaurene experiments. Low-volatility gases and particles deposited on the chamber walls can then lead to compounds slowly off-gassing from the walls. In reactions with OH, these C 2−9 precursors can then form RO 2 that terminate by reacting with the C 19−20 radicals, ending up forming the observed C 21 -C 29 HOM dimers. As most of our experiments were done with very low kaurene concentrations (<1 ppb), the contaminants could amount to an equally high (or even higher) OH sink than the kaurene itself. By comparing the most abundant observed contaminants and the C 24 -C 28 dimers, we also tried to assess which RO 2 species from kaurene were the most important. In most cases, C 19 H 29 O 3 and C 19 H 31 O 2 seemed to be the most plausible radicals to explain the most abundant observed dimers. These were also thought to be the main contributors to C 38 -C 39 dimers of kaurene oxidation. Unfortunately, due to their low oxygen content, we cannot directly detect these radicals with nitrate CI-APi-TOF.
During most of the experiments, all HOM detected by nitrate CI-APi-TOF exhibited similar temporal behaviors as the kaurene time series (Fig. 1), as expected due to the low variability in O 3 concentrations. The products measured by the nitrate CI-APi-TOF reacted rapidly to changes in the injected kaurene, while the less oxidized products detected by the Vocus showed a slower response to kaurene changes, especially when kaurene was removed.

HOM molar yield estimation
We calculated the molar yield for HOM from kaurene ozonolysis (without OH scavenger) as described in Sect. 2.5. Figure 5 shows the measured HOM concentration as a function of the kaurene ozonolysis rate (k 1 × [kaurene] × [O 3 ]). In the kaurene ozonolysis experiments at low oxidation rates (Fig. 5b), the total HOM concentrations had a near-linear dependence on the amount of kaurene reacting with O 3 . We can estimate a molar HOM yield of ∼ 2 % in this range, though keeping in mind the large uncertainties in the quantification of kaurene and HOM. At higher oxidation rates (Fig. 5a), the increase in total HOM concentration exhibited a logarithmic dependence. The slower growth of HOM concentration at higher oxidation rates might be explained by increasing sinks in the chamber, e.g., due to particle formation. The RO 2 lifetimes will also decrease, which may hamper the autoxidation when bimolecular reactions become faster. When the entire data range was included, the molar HOM yield ranged between 0.1 % and 10 %.
HOM yields from kaurene ozonolysis, as for any other system, will always depend on the atmospheric conditions, in particular the type and concentration of bimolecular reaction partners (mainly NO, RO 2 , HO 2 ). The reaction partner will determine the possible branching pathways, while their absolute concentrations will affect RO 2 lifetimes and thereby the potential to undergo autoxidation. According to previous studies, the reported kaurene concentrations in the atmosphere range from ppq to a few ppt (Li et al., 2020;. In our study, concentrations were higher (Fig. S3a  and c), but the ∼ 2 % HOM yield was estimated based on periods when kaurene concentrations were <50 ppt (Fig. S3b). As the atmosphere is much more complex, with various other species present, it is possible that the RO 2 concentrations in our experiments are, in fact, lower than in atmosphere, despite the kaurene concentrations being higher. Nevertheless, we think that the oxidation conditions in our experiments are quite close to atmospheric conditions and that the potential bias in HOM yields due to this is marginal when compared to the overall uncertainty of HOM quantification (Sect. 2.4.2). There are two main exceptions to this. Firstly, if NO concentrations would be very high, it could lead to very different oxidation pathways, with unknown effects on the HOM yields. Secondly, in the atmosphere it would be much less likely for two kaurene-derived RO 2 species to react, as the majority of the RO 2 would come from different precursors; thus, the C 37−40 dimers would be unlikely to form in the atmosphere.
Although large uncertainties are included in the derivation of our HOM yield from kaurene ozonolysis, yields of a few percent have been reported for a variety of other BVOCs earlier (Table 1). More studies are needed to better quantify the HOM formation in this system and get further chemical insights into the detailed oxidation pathways. Although the global emission budgets of diterpenes remain unquantified, when combining the limited emission rates reported previ-   ously (Matsunaga et al., 2012;Yáñez-Serrano et al., 2018) with the HOM yield estimated here, it is possible that HOM formation from diterpenes, as well as their influence on SOA formation, are of larger importance than previously thought.

Volatilities of the oxidation products
To deduce the ability of the oxidation products from kaurene ozonolysis to contribute to aerosol formation, we investigated their behavior as a function of the available condensation sink (CS) in the chamber. This method has been applied to study the volatilities of the α-pinene and cyclohexene oxidation products by Peräkylä et al. (2020) and Räty et al. (2021). In this methodology, we assume that the condensation onto seed particles is the main driver of the concentration changes of oxidation products during the experiment, although other changes (e.g., the changes of precursors and intermediates due to seed injection) can have a non-negligible influence (Peräkylä et al., 2020). This approach does not provide direct estimates of volatilities, but it does provide a good separation between readily condensing vapors (accommodation coefficient near unity) and the more volatile ones that are unaffected by the condensation sink (i.e., with very low accommodation coefficients). In this experiment, kaurene, O 3 , and seed particles were injected into the chamber until it reached a steady state (referred to as S1 in Figs. 1c and S4b). The injection of seed particles was then stopped, causing a decrease in the CS and a corresponding increase in the concentrations of several compounds, ultimately reaching a new steady-state condition (referred to as S2 in Figs. 1c and S4b). For each observed compound, we then calculated the concentration ratio at S1 and S2, henceforth termed the fraction remaining (FR) (Fig. 6).
From Fig. 6a, we found a clear transition of volatilities when the molecular mass increased, in line with previous findings for α-pinene and cyclohexene systems (Peräkylä et al., 2020;Räty et al., 2021). Interestingly, all species detected by the Vocus had an FR around 1, i.e., not affected by the CS, suggesting that the Vocus primarily detects VOCs and IVOCs (intermediate-volatility organic compounds). Kaurene (and all its fragments) were in this range, as expected. Signals with FR much above unity may either be semivolatile compounds that never reached a steady state during the experiments or molecules that were preferentially formed in the particle phase before evaporating into the gas phase. A sharp drop in the FR, from ∼ 1 to ∼ 0.1, was found in the mass range 300-350 Th, and the oxidation products with more oxygen tended to have lower values. This trend is consistent with that in the mass range above 350 Th (Fig. 6b), where 5630 Y. Luo et al.: Oxidation product characterization from ozonolysis of the diterpene ent-kaurene Figure 6. Seed injection behavior; the "fraction remaining" corresponds to the ratio of the signal with seed particles to the signal without seed. We emphasize that the molar mass indicated on the x axis is without the reagent ions, unlike many previous studies (Peräkylä et al., 2020;Räty et al., 2021). Species with less than 15 carbon numbers in the composition are colored with gray, while the color scale indicates the oxygen content for species with 15 or more carbon atoms. The areas of the circles (nitrate CI-APi-TOF) and diamonds (Vocus) are scaled linearly to the magnitude of each compound's signal when there were no seed particles in the chamber (different scaling for the two instruments). Compounds with the average signal intensities lower than 3 times the standard deviation are excluded from the plot. The text at the top shows roughly which compound groups are found in which parts of the mass range. The "Vocus fragments" refers to the kaurene fragments, and "Vocus OP" refers to C 18−20 oxidation products. "Nitrate monomers", "Nitrate C 21−29 dimers", and "Nitrate dimers" represent the C 18−20 HOM monomers, C 21−29 HOM dimers, and C 37−40 HOM dimers measured with nitrate CI-APi-TOF, respectively. A zoomed-in view of the mass range above 300 Th is shown in (b). most HOM were found. This indicates that these species condensed efficiently and were mainly LVOCs (low-volatility organic compounds) or ELVOCs (extremely low-volatility organic compounds). Similar results were reported for HOM formed from α-pinene ozonolysis by Peräkylä et al. (2020). The relatively large spread of the FR above 350 Th may indicate that the seed particles also impacted the formation pathways, in addition to providing a larger CS. For example, some longer-lived RO 2 might also have been lost to the particles, as discussed by Peräkylä et al. (2020).
In the mass range 200-300 Th, the FR of nitrate CI-APi-TOF species were also around one even though some have quite high oxygen content (Fig. S5a). These compounds were primarily contaminants from the chamber walls, typically with less than 10 carbon atoms (Fig. S5c). While the high Oatom content might suggest lower volatility, the FR around unity suggests that the CS was not a noticeable sink for these species. On the one hand, this is credible since their existence in the chamber during these experiments means that they evaporate efficiently from the chamber walls. On the other hand, it is also possible that there is large storage on the walls, and they have reached a balance between the gas-phase and wall concentrations. When these species condensed onto the seed particles, more came out from the chamber walls, keeping the measured gas-phase concentra-tions stable. Either way, although not of particular importance for this study, the mere existence of these compounds in the continuously flushed chamber strongly suggests that they are SVOCs (semivolatile organic compounds).
In summary, these results show that the oxidation products of kaurene ozonolysis with approximately five or more oxygen atoms were mainly LVOCs or ELVOCs, and will thus efficiently condense to form organic aerosols. In line with earlier studies (Räty et al., 2021), we also found that the oxygen content is the most important parameter affecting the volatility. Despite the longer carbon chain, thus higher mass, of the diterpenes compared to, for example, monoterpenes or cyclohexene, the transition from condensing to non-condensing was not dramatically different. For those compounds (C 6 and C 10 ), similar experiments indicated that approximately seven O atoms were needed, whereas for the C 20 kaurene the number was around five.

Fragmentation and NH + 4 clusters in Vocus PTR-TOF
Fragmentation is a ubiquitous process when measuring VOCs with PTR instruments, complicating the identification and interpretation of the mass spectra. Previous studies have shown that the settings of PTR-MS (in essence E/N ), the types of VOCs, and environmental conditions (e.g., RH)

Conclusions
This study presents the first characterization of gas-phase oxidation products from diterpene ozonolysis. We studied the reaction of the diterpene kaurene with O 3 in a simulation chamber, monitoring the products using a Vocus PTR-TOF and a nitrate ion-based CI-APi-TOF. A group of HOM with carbon numbers in the range C 21 -C 29 was surprisingly observed in this study, and their temporal behaviors were similar to those of other oxidation products. These C 21 -C 29 HOM were assumed to be accretion products formed via the reactions of kaurene-derived RO 2 with some smaller RO 2 formed from OH oxidation of evaporated contaminants from the chamber walls. The HOM molar yield from kaurene ozonolysis without using an OH scavenger was estimated to be ∼ 2 %, comparable to those reported for other terpenes, though the exact value comes with a large uncertainty due to the uncertain quantification of both the HOM and the precursor. The volatilities of the oxidation products formed in kaurene ozonolysis were also investigated by adding seed aerosol as an additional CS. None of the species detected by the Vocus seem to be affected by the increasing CS, indicating they were VOCs and IVOCs. However, most HOM, with five or more oxygen atoms, detected with the nitrate CI-APi-TOF were assumed to be LVOCs or ELVOCs, because they efficiently condensed onto the seed particles.
For detecting diterpenes in air, the exact behavior of the VOCs in a PTR instrument is of importance. Here, fragment distributions of kaurene in a Vocus PTR instrument were reported. We observed seven main kaurene fragments: C 10 H + 15 , C 11 H + 17 , C 12 H + 19 , C 13 H + 21 , C 14 H + 23 , C 15 H + 23 , and C 16 H + 25 . The lost neutral fragments are consistent with those observed earlier for sesquiterpenes. The fragmentation process of kaurene was found to be a function of the Vocus FIMR pressure, with the fragment signals increasing rapidly and even surpassing that of the parent ion when FIMR pressure decreased. In addition, the fragment distribution shifted towards smaller fragments as the FIMR pressure decreased. Interestingly, we also observed that kaurene and two oxidation products were detected as clusters with NH + 4 , in addition to the typical protonated molecular ions in the Vocus. The ratio between the NH + 4 cluster signal and the protonated signal was less than 2 % for C 20 H 32 , but the ratios ranged from 15 % to 70 % for two oxidation products C 19 H 30 O 1,2 . Future studies should evaluate the capability of different methods to detect diterpenes and the corresponding oxidation products, as it is possible that these larger molecules may be more efficiently detected as clusters.
Based on the estimated HOM yield and low volatility of kaurene oxidation products, the influence of diterpenes on SOA formation may be more important than previously thought. However, more studies are required to determine the global emission budget of diterpenes. In addition, fur-ther work should try to understand the oxidation pathways and HOM contributions of other diterpenes in the laboratory and the ambient air, in order to produce a better picture of the role of diterpenes in atmospheric SOA formation.
Data availability. Data are available upon request by contacting the corresponding authors.
Author contributions. ME and AMYS initiated the study. ME and OG designed the experiments. YL, OG, HL, FG, YZ, MM, and OP conducted the measurements and operated the chamber. YL performed the main data analysis of gas-phase and related model results and was supervised by ME. OG helped with the kaurene-O 3 rate coefficient estimate performed by YL. FG analyzed the particlephase data. APP and AL performed and analyzed the TORM. YL plotted the main figures and wrote the original draft. All the authors, including JP, discussed the results and commented on the paper.

Competing interests.
The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open-access funding was provided by the Helsinki University Library.
Review statement. This paper was edited by Ivan Kourtchev and reviewed by two anonymous referees.