An intercomparison study of four different techniques for measuring the chemical composition of nanoparticles

Abstract


Introduction
So far there is no well-established instrument and technique to measure the complete chemical composition of ultrafine (< 100 nm) Secondary Organic Aerosol (SOA) particles.However, several analytical techniques have recently been developed in order to better advance our understanding on their chemistry.Techniques that are capable of measuring sub-30 nm particles include the Volatile Aerosol Component Analyzer (VACA) (Curtius et al., 1998), the Thermal Desorption Chemical Ionization Mass Spectrometry (TDCIMS) (Voisin et al., 2003), the Nano Aerosol Mass Spectrometer (NAMS) (Wang et al., 2006), the Aerosol time-of-flight mass spectrometer (Laitinen et al., 2009), the inlet for the size-resolved collection of aerosols (Phares and Collier, 2010), the Chemical Analyzer for Charged Ultrafine Particles (CAChUP) (Gonser and Held, 2013), the Electrostatic Precipitation-Electrospray Ionization Mass Spectrometry (EP-ESI-MS) (He et al., 2015), the Droplet Assisted Inlet Ionization (DAII) (Horan et al., 2017), and the Online Aerosol Chemical Characterization by Extractive Electrospray Ionization−Ultrahigh-Resolution Mass Spectrometer (EESI-Orbitrap) (Lee et al., 2020).Singleparticle analysis by mass spectrometry methods based on aerodynamics, light scattering and laser desorption ionization are suitable for particles with larger sizes.These methods include for example, the Aerosol Mass Spectrometer (AMS) described in Jayne et al. (2000), and the suite of single particle methods described in the review by Bzdek et al. (2012).The detection of particles with d < 100 nm using these techniques becomes difficult, because the scattering efficiency decreases when the particle diameter becomes smaller.
Using the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN), we used simultaneously four different techniques for measuring the chemical composition of ultrafine particles and inter-compare the results.
None of the techniques presented in this work represents the perfect instrument.In fact, a perfect instrument would be the one that is able to measure quantitatively all the hundreds of organic compounds that are present not only in the newly formed particles in the lab, but also in aerosol particles present in the ambient, i.e., with larger particles being present as well and at low mass concentration of the ultrafine particles.A perfect instrument should also be able to identify the molecular structures (including their isomeric and spatial configuration) at high time resolution and in real-time.Such an ideal instrument does not exist; and at present, a combination of techniques is required for a more complete characterization of SOA (Hallquist et al., 2009).
In order to compare the techniques mentioned above and to gain insights into their limitations (e.g., due to decomposition during evaporation, different ionization techniques, etc.), we performed α-pinene ozonolysis experiments at -50 ºC, -30 ºC and -10 ºC.For the experiments at -50 °C and -30 °C TD-DMA, FIGAERO and EESI-TOF were intercompared while for the experiment at -10 ºC FILTER, FIGAERO and, EESI-TOF were inter-compared.We carried out the in-depth inter-comparison by a) comparing the observed composition as a function of oxidation state and carbon number, b) estimating the volatility and comparing the fraction of ultralow-volatility (ULVOC), extremely low-volatility (ELVOC), low-volatility (LVOC), semi-volatile (SVOC), and intermediate-volatility (IVOC) organic compounds, and c) comparing the thermograms (for the thermal desorption techniques: TD-DMA and FIGAERO), and by performing positive matrix factorization (PMF) analysis to the thermograms.
Because the four techniques provide chemical composition, and more specifically the carbon, hydrogen and oxygen content (CHO), we determined the carbon oxidation state (OSc), which is a metric for the degree of the oxidation of organic species in the atmosphere (Kroll et al., 2011).It is calculated based on the ratios O:C and H:C and is useful to describe organic mixtures upon oxidation processes.In addition, we estimated the volatility (as introduced by Donahue et al. (2011) and modified by Stolzenburg et al. (2018)), and determined to which volatility classes the detected compounds belong.
Regarding the thermal desorption methods (TD-DMA and FIGAERO) we investigated the thermal behavior of the detected species.Both instruments first collect particles and subsequently the particles are evaporated in order to transfer their constituents to the gas-phase.When a temperature ramp is applied, the species that are adsorbed on the surface gradually desorb (as represented on a thermogram).In order to evaluate whether the thermal desorption methods lead to significant decomposition during evaporation, we applied a method called positive matrix factorization (PMF) (Paatero and Tapper, 1994;Buchholz et al., 2020), in which a dataset matrix is expressed in terms of the sum of factors matrices and a residual matrix.Thermal decomposition on FIGAERO particle phase data has been reported previously.D' Ambro et al. (2019), observed that some of the major components of IEPOX (isoprene-derivated epoxydiol) in SOA such as C5H10O3 and C5H12O4 are likely artifacts of thermal decomposition.Lopez-Hilfiker et al. (2015), Stark et al. (2017), Wang and Hildebrandt Ruiz (2018) have addressed the importance of considering thermal decomposition for assessing the chemical composition and volatility properties of SOA in techniques in which the aerosol is heated before or during the analysis.In this study, we evaluate thermal decomposition of TD-DMA samples for the first time and intercompare the results with FIGAERO.
Lastly, we present an overview on the advantages and disadvantages for the different methods.All methods presented here agreed on the most dominant compounds that are found in the nanoparticles.Nevertheless, they do see different parts of the organic spectrum.Therefore, the techniques are complementary.The four techniques described in this work are feasible for measuring sub-100 nm SOA particles and represent an important analytical development.

The CLOUD chamber experiment
We conducted experiments in the CLOUD chamber at CERN to study pure biogenic new particle formation (NPF), without the presence of nitrogen oxides (NOx).The CLOUD chamber is a stainless-steel cylinder with a volume of 26.1 m 3 which has been extensively described by Kirkby et al. (2011) and Duplissy et al. (2016).To create the particles, NPF was induced by continuously adding α-pinene and ozone into the chamber.The monoterpene concentration was regulated by an evaporation source, in which dry nitrogen (N2) passes through an evaporator containing liquid α-pinene at a precisely controlled temperature.Ozone was produced by exposing cryogenic O 2 to UV light and was introduced directly into the chamber via a separate line.The relative humidity was adjusted with a temperature-controlled Nafion humidifier using ultrapure Millipore water.All the precursor gases were homogenously mixed in the chamber by two magnetically driven Teflon fans placed at the top and at the bottom of the chamber.The temperature was kept constant by an insulated thermal housing, which surrounds the chamber.The α-pinene mixing ratio was measured by a proton transfer reaction time-of-flight (PTR-TOF) mass spectrometer (Graus et al., 2010;Breitenlechner et al., 2017), whereas ozone was measured by a TEI 49C ozone analyzer (Thermo Environmental Instruments).The experiments relevant for this work were performed at -50 ºC, -30 ºC, and -10 °C.The α-pinene mixing ratio ranged between 1 and 8 ppbv and ozone was approximately 100 ppbv.
Table 1 presents the most important features for the instruments used in this work.We categorize the techniques according to certain criteria: continuous or discontinuous operation mode, evaporation method, phase measured, ionization technique, reagent ion, target substances, occurrence of significant thermal decomposition and, whether the technique allows to perform size-resolved analysis of aerosol particles.

Thermal Desorption-Differential Mobility Analyzer (TD-DMA) coupled to a NO3 -chemical ionizationatmospheric-pressure-interface-time-of-flight (CI-APi-TOF) mass spectrometer
The TD-DMA coupled to a NO3 -CI-APi-TOF analyzes the chemical composition of nanoparticles in a semi-continuous mode of operation.The design and characterization have been reported by Wagner et al. (2018).This method allows gas-and particle-phase measurements using the same ionization technique.Individual results of gas-and particle-phase comparison of the same chemical system as in this study were reported in Caudillo et al. (2021).While the gas-phase measurement is taking place with the CI-APi-TOF mass spectrometer, the TD-DMA samples particles from the chamber.The TD-DMA can perform size-resolved and non-size-resolved measurements.In any case, the particles are first charged by an X-Ray source and then transferred to the differential mobility analyzer (DMA) unit.When a size-resolved measurement is desired, a specific voltage is applied to the central electrode inside the DMA unit.Subsequently, a sheath flow will carry only particles with specific electrical mobility and will conduct them through the DMA.In contrast, during a non-size-resolved measurement, no voltage to the central electrode and no sheath-flow are applied, and a fraction of the particles charged by the X-ray source will pass through the DMA unit.For the experiments reported in the present study, we performed non-sizeresolved measurements in order to maximize the mass collected and to be comparable to the other methods that do not perform size selection.
The particle collection takes place by electrostatic precipitation on a platinum/rhodium (90:10) filament placed inside the central electrode.After a certain collection time (~3 hours for the experiment at -30 ºC, and ~5 hours for the experiment at -50 ºC), an electric current is applied to the filament, which causes its direct heating.We estimate based on the filament resistance, that the temperature gradually increased up to approximately 600 ºC in a period of ~1 minute (details of the heating curve will be discussed in Section 3.2).
The vapors that evaporate from the heated particles are carried by an N2 flow to the nitrate CI-APi-TOF for chemical composition analysis.For chemical ionization of the vapors, nitrate reagent ions (HNO3)n NO3 -with n = 0-2 are created by a corona discharge (Kürten et al., 2011;Kürten et al., 2014).Some of the vapor molecules are ionized and subsequently detected by the APi-TOF mass spectrometer.A second heating cycle of the particle collecting filament is performed afterwards (without particle collection) in order to estimate the instrumental background due to the heating of the inlet line; this enables a more accurate determination of particle constituents.The particle constituents are estimated by subtracting the second heating from the first heating.Besides the particle and background estimations, a second heating up to ~600 ºC ensures that the filament is clean and avoids memory effects for the next measurement.With the nitrate ionization technique, sulfuric acid (Jokinen et al., 2012), iodic acid (He et al., 2021), methane sulfonic acid (Shen et al., 2022), and highly oxygenated molecules (Kirkby et al., 2016;Simon et al., 2020) can be detected.

Filter Inlet for Gases and AEROsols (FIGAERO) coupled to an I -high-resolution time-of-flight chemicalionization mass spectrometer (HRToF-CIMS)
The Filter Inlet for Gases and AEROsols (FIGAERO) coupled to an I -high-resolution time-of-flight chemical-ionization mass spectrometer (HRToF-CIMS) was first described by Lopez-Hilfiker et al. (2014) and optimized for CLOUD operation conditions by Wang et al. (2020).FIGAERO uses a multi-port to measure in alternation both gas-and particle-phase following the same general procedure as the TD-DMA/CI-APi-TOF.While the gas-phase is analyzed, particle collection takes place on a polytetrafluoroethylene (PTFE) filter, and after a certain collection time (in this study 15 minutes), the filter is automatically moved into the ion-molecule reactor and exposed to a pure N2 gas flow.The N2 flow is gradually heated to evaporate the particles by thermal desorption using a temperature-programmed heating curve.For the measurements reported in this study, the temperature was slowly ramped from room temperature up to 150 ºC in approximately 15 minutes (an example of the heating curve is discussed in Section 3.2).When the heating cycle ends a new collection starts, and the process repeats.Thus, 15-minute collection period is followed by 15-minute desorption-period, which implies two collections per hour (resulting in a resolution of 30 minutes).The detection technique is based on iodide-molecule adduct ionization.Iodine ions I -, I3 -and (H2O)I -are generated from a solution of methyl iodide (CH3I) and a Po-210 ion source (Lee et al., 2014).With this soft ionization technique, the FIGAERO HRToF-CIMS can detect intermediate oxygenated organic molecules (Wang et al., 2020;Stolzenburg et al., 2018;Lopez-Hilfiker et al., 2016), organosulfates, and inorganic acids such as sulfuric acid and nitric acid.

Extractive electrospray Na + ionization time-of-flight mass spectrometer (EESI-TOF)
The Extractive electrospray ionization time-of-flight mass spectrometer (Lopez-Hilfiker et al., 2019) is a technique used for online particle-phase measurements without batch collection.This technique aims to provide the chemical composition of organic particles in real time (1 second).It is also possible to measure the gas-phase by using the dual configuration (Lee et al., 2022).In the beginning of the sampling process, the aerosol sample passes through the inlet line where a carbon denuder is located to remove the gas-phase molecules.The particles then collide with electrospray droplets and the soluble components are extracted and ionized.Then, solvent evaporation occurs in a heated stainless-steel capillary, leading to Coulomb fissioning and/or direct ion desorption.Finally, ions enter a time-of-flight mass spectrometer.The EESI-TOF uses here an electrospray solution of pure water doped with 100 ppm NaI and is running in the positive ion mode.This enables the measurements of SOA species as adducts with Na + .With this ionization method, most organic compounds that are relevant for atmospheric SOA particles can be analyzed, such as intermediate oxygenated organic molecules with the exception of species that are not oxygenated and organosulfates (Lopez-Hilfiker et al., 2019).

Offline analysis of filters (FILTER) using Ultra-high-performance liquid chromatography (UHPLC) and heated electrospray ionization (HESI) coupled to an Orbitrap high-resolution mass spectrometer (HRMS)
This procedure was optimized and described in detail by Ungeheuer et al. (2021).The method enhances the separation of organic compounds with high-resolution, and enables the determination of the accurate mass.The analysis consists mainly of four steps: sampling, extraction, separation, and detection.
First, the particles were collected from a flow of 5 l min -1 on a 47 mm diameter Emfab™ Filter (Pall Life Science, USA) during approximately 17 hours for the experiment at -10 ºC (see Table 2 and Figure S1).After sampling, the filter was stored at -18 °C to avoid possible losses by evaporation.The filter was cut into small pieces (approximately 3 x 3 mm) and extracted two times in 0.2 ml solution (mixture of 90 % water and 10 % methanol) for 20 min.After each extraction step, the extract was filtered through a syringe filter (PTFE with a pore size of 0.2 µm).For chromatographic separation a gradient of ultrapure water (eluent A, Milli-Q Reference A+, Merck Millipore) and methanol (eluent B, Optima LC/MS Grade, Fisher Scientific) was applied.Both eluents were mixed with 0.1 % formic acid (v/v) for improved chromatographic performance.
The injection volume was 5 µl, the flow rate was 400 µl min -1 , and the temperature was 40 °C.The gradient started with 1 % eluent B (0-0.5 min), increased linearly to 99 % B (0.5-14 min), stayed at 99 % B (14-16 min), backflushed in 1 min, and equilibrated in 3 min, resulting in a total run time of 20 min.Negative ionization mode was used for the detection, in which the molecular ions [M-H] -are produced by deprotonation.The ion source setting used for this purpose were: -3.5 kV spray voltage, 40 psi sheath gas, 8 arbitrary units auxiliary gas, and 262.5 °C capillary temperature.The scan range in full MS was 50-750 m/z with a mass resolving power of about 70k at m/z 200.For data-dependent MS 2 (ddMS 2 ) the resolution was 17.5k.
Fragments were produced in a higher-energy collisional dissociation (HCD) cell with stepped collision energies of 15, 30, and 45 eV.

Data processing
TD-DMA, FIGAERO and EESI-TOF data were processed using IGOR Pro 7 (WaveMetrics, Inc., USA) and Tofware (Version 3.1.2,Aerodyne Inc., USA).The data from the offline method was processed with Compound Discoverer 3.2 (Thermo Fisher Scientific).The postprocessing was done using MATLAB R2022a (MathWorks, Inc., USA).TD-DMA data was corrected by the mass-dependent transmission efficiency in the mass classifier (Heinritzi et al., 2016) and normalized by the nitrate reagent ions.FIGAERO data was averaged to 1 minute and normalized by the reagent ions.EESI-TOF signals were averaged to 10 seconds and normalized by the most abundant electrospray ion (NaINa + ).In order to align the sampling times for the different techniques and perform a more direct particle-phase comparison, we selected as a reference the TD-DMA and FILTER collection times.Thus, for the comparison at -30 ºC and -50 ºC, EESI-TOF particle signals were averaged, and FIGAERO particle signals were integrated during the period where the TD-DMA collected particles (~3 and 5 hours, respectively).While for the comparison at -10 ºC, the EESI-average and FIGAEROintegration period corresponded to the time where the particles were collected with the FILTER for the offline analysis (~17 hours).Figure 1 provides an overview of a representative experiment at -30 ºC.The overview of the experiments at -10 ºC and -50 ºC are shown in Fig. S1 in the supplement.Table 2 summarizes the sampling conditions during the experiments reported in this study (for the purpose of the intercomparison), including the particle number concentration, mass concentration and median mass diameter (MMD) calculated from the SMPS.The MMD indicates the particle size measured by the SMPS for which 50 % of the aerosol mass is contained in smaller particles and 50 % is contained i n larger particles.
From Table 2 it can be seen that for the experiments at -30 °C and -50 °C MMDav < 100 nm, while for the experiment at -10 °C, MMDav ~106 nm.

Positive matrix factorization (PMF) analysis
One of the main questions we want to answer in this work is whether the thermal desorption methods (TD-DMA and FIGAERO) experience significant decomposition during the desorption process.To answer this question, we utilized positive matrix factorization (PMF) analysis.This method was originally described by Paatero and Tapper (1994) for analyzing time series of variable (e.g.mass spectra data) from ambient observations, and it was implemented by Buchholz et al. (2020) to thermal desorption data for identifying different processes during particle evaporation.We therefore applied the same procedure as Buchholz et al. (2020) to the TD-DMA and FIGAERO thermal desorption profiles (for the α-pinene oxidation experiment at -30 ºC and 20 % RH only).For the analysis, we processed separately 1-second TD-DMA and 1second FIGAERO thermograms.Since the FIGAERO measures in a semi-continuous mode, we chose a representative thermogram.Both TD-DMA and FIGAERO data sets were background subtracted.For the TD-DMA background, we used the second heating cycle that is performed immediately after the first heating (described in Section 2.1.2).For the FIGAERO background, we used a period where no significant particle load was present in the chamber (at the beginning of the experiment).We considered only the organic compounds and excluded the reagent ions for this analysis.We ran the PMF software using the CNerror scheme (based on the noise of each ion) and up to 10 different solutions.4-factor TD-DMA and 6-factor FIGAERO solutions (discussed later in Section 3.2.1)were chosen as the most interpretable results by a) comparing the residuals and by looking at which solution captured the total signal and certain species the best (e.g., C8H12O4), b) by finding an equilibrium between good reconstructed signal and physically interpretable results.This means, that for the solutions presented here, likely a higher number of factors improve the residuals, nevertheless, we chose the solution with the smallest number of factors that can still provide realistic information.

Chemical composition comparison
Figure 2 shows the OSC, calculated as OSC = 2  O:C -H:C (an approximation stated by Kroll et al. (2011)) against the number of carbon atoms for α-pinene oxidation products in the particle-phase at -30 ºC and 20 % RH as measured by three different techniques: TD-DMA, FIGAERO, and EESI-TOF.For all techniques, the highest intensities correspond to compounds with 10 carbon atoms (C10) for which the oxidation state varies between 0.5 and -1.5 depending on the measurement technique.Compounds with more than 10 carbon atoms were also detectable by the TD-DMA, FIGAERO and EESI-TOF.The TD-DMA and FIGAERO detected compounds with less than 5 carbon atoms and OSC > 0 which, in contrast, are not detected by the EESI-TOF (this feature will be discussed in Section 3.2.1).In order to simplify the comparison, we calculated the fraction of species containing less than 10 carbon atoms (C<10), 10 carbon atoms (C10), and more than 10 carbon atoms (C>10), since this can provide an insight of the detected fraction of monomers and dimers for each technique (Figure 2d).Approximately 42 %, 32 %, and 23 % of the signals correspond to C<10, 47 %, 65 % and, 53 % to C10 and, around 11 %, 3 %, and 24 % to C>10 measured by the TD-DMA, FIGAERO and EESI-TOF, respectively.Figure S2 in the Supplement displays the results for the experiments at -50 ºC.In every case, C10 represents the highest fraction detected by all the techniques in this experiment.Nevertheless, we do see significant differences between the techniques for C<10, C10, and C>10.
The chemical ionization utilized by the TD-DMA (nitrate, NO3 -) is more sensitive towards highly oxygenated species, while the FIGAERO (iodide, I -) detects intermediately oxygenated species with higher sensitivity.In Figures S3 and   S4 in the Supplement, we present the results (number of oxygen atoms vs number of carbon atoms) at -30 ºC and -50 ºC, respectively.From the figures (S3 and S4) we observe that more oxygenated species contribute more to the total signal on TD-DMA than in FIGAERO, this observation is consistent with the sensitivity that one would expect according to the chemical ionization.The electrospray ionization (Na + ) for the EESI-TOF is usually more sensitive towards intermediate oxygenated species, even though in the results presented here, it seems to capture very well the whole spectrum.
Besides the reagent ion selectivity, several factors can explain the quantitative differences.For example, both the TD-DMA and FIGAERO detect a lower fraction of C>10 (11 % and 3 % compared to 24 % for the EESI-TOF).The TD-DMA and FIGAERO techniques are based on thermal desorption, which may cause decomposition of thermally unstable compounds during evaporation (as discussed in more detail in Section 3.2.1).On the other hand, the resulting fractions of C>10 can be influenced by the ionization method employed: chemical ionization and electrospray are soft ionization techniques, for which one can expect little fragmentation.Thus, we presume that thermal decomposition during evaporation could be the most significant factor that explain these differences.
Figure 3 presents the results for the experiment conducted at -10 ºC (α-pinene oxidation products at 80 % RH) for particles collected on a FILTER (Fig. 3a) and later analyzed by the UHPLC-HESI-HRMS method, and as measured by the FIGAERO and EESI-TOF (Fig. 3b and 3c, respectively).Overall, fewer compounds are detected by UHPLC-HESI-HRMS than by FIGAERO and EESI-TOF.The highest intensities in Fig. 3a (FILTER) correspond to C8H12O4, C9H14O4, C10H16O3-6, C17H26O8, and C19H28O7.Ions with the same formulas are also detected by the FIGAERO and EESI-TOF, but the contribution to the total signal differs.The results at -10 ºC (number of oxygen atoms vs number of carbon atoms) are shown in Figure S5 in the Supplement.By applying the UHPLC-HESI-HRMS method, it is possible to distinguish between compounds with identical chemical formula (isomers).For the experiment reported here, two isomers for C8H12O4, C10H16O3, and C10H16O5, as well as three isomers for C10H16O4 and C10H16O6 were detected.The detection of these isomers is enabled by the chromatographic separation (their interaction with a reversed-phase column results in different retention times and therefore makes the separation feasible).However, complementary experiments are needed to investigate the molecular structure.Furthermore, Figure 3d shows the contributions of the C<10, C10 and C>10 fractions to the total signal for the FILTER, FIGAERO and EESI-TOF.At -10 ºC, the fraction of compounds with more than 10 carbon atoms (C>10) has the smallest contribution to the total signal (14 %, 5 % and, 24 % for the FILTER, FIGAERO and EESI-TOF, respectively).The fractions C<10 and C10 do not seem to have a clear tendency, they both contribute substantially to the total signal in each technique.
For the methods utilizing electrospray ionization (EESI-TOF and UHPLC-HESI-HRMS), the polarity plays an important role.The EESI-TOF ran in the positive mode, allowing the SOA species to be detected as adducts with Na + .While HESI utilized in the offline method ran in the negative mode, in which molecular ions are produced by deprotonation.Surdu et al. (2021) used an aerosol growth model based on the condensation of organic vapors, and demonstrated that the chemical composition measured by the EESI-TOF (in the positive mode as in this study) is consistent with the expected condensed oxidation products with small differences.On the other hand, Ungeheuer et al. ( 2021), who utilized the UHPLC-HESI-HRMS method, reported that the detection of ester molecules was accomplish only in the positive ionization mode (1000 compounds were detected in the positive mode while only 16 were detected in the negative mode).Besides the polarity, several other factors can influence the electrospray ionization response, such as, analyte chargeability and surface activity (Cech and Enke, 2001), and therefore impact on the completeness of the chemical composition of aerosol particles.

Volatility classes
We estimated the volatilities of the detected particle-phase compounds and associated them with defined volatility classes.We used the parametrization introduced in Donahue et al. (2011) and modified by Stolzenburg et al. (2018).This approach has been also discussed in Simon et al. (2020), and Wang et al. (2020).The volatility was approximated from the number of carbon and oxygen atoms in the specific molecules, and it was first defined at 300 K.By using the Clausius-Clapeyron equation, the volatility was then shifted according to the corresponding experimental temperature.The evaporation enthalpy was approximated according to Donahue et al. (2011) and Epstein et al. (2009).Thereafter the volatility was associated with any of the following classes: ultralow-volatility (ULVOC), extremely low-volatility (ELVOC), low-volatility (LVOC), semivolatile (SVOC), and intermediate-volatility (IVOC) organic compounds.
Figures 4 and 5 show the contribution of each volatility class to the total particle signal for the corresponding experiment and for each technique.The results at -50 ºC are reported in the Figure S6 in the Supplement.For the experiment at -30 ºC (Fig. 4), LVOC constitute the higher fraction for the TD-DMA, FIGAERO and, EESI-TOF.EESI-TOF detects the highest fraction of ULVOC (~12 % compared to 6 %, and 2 % measured by TD-DMA and FIGAERO, respectively).An IVOC fraction (~8 %) is only detected by the EESI-TOF.At -10 ºC (Figure 5), semi-volatile organic compounds (SVOC) contribute the most to the total particle signal for the FILTER, FIGAERO and EESI-TOF.Very small fractions of ULVOC are also detected by all the techniques.The EESI-TOF detects a higher fraction of IVOC (20 %) than FILTER and FIGAERO (2 % and 4 %), respectively.Taking into account the particle load (~1-3 ug m -3 ) and size of the particles (diameter < 100 nm), it is possible that a significant fraction of IVOC measured by the EESI-TOF results from measurement artefacts, as seen in previous studies using EESI (Surdu et al., 2021).Several reasons (or a combination of them) can explain this feature: a) the Na + ionization technique may be more sensitive to lower oxygenated organic compounds than the I -or NO3 -techniques; b) likely, some amount of the gas-phase broke through the charcoal denuder (although its efficiency is > 99 %) and reached the detector.Lee et al. (2021a) reported that the EESI-TOF is more sensitive toward gas-phase analytes as compared to their particle-phase counterparts and, c) it can possibly occur some ion-induced fragmentation.
Overall, we observed that the contribution of the lowest volatility classes (ULVOC, ELVOC and LVOC) increases as the temperature decrease.This observation reflects two opposing temperature effects, as discussed in Ye et al. (2019) based on FIGAERO results: autoxidation and thus the extent of oxidation is reduced at low temperature, but any given compound is much less volatile at low temperature because of the strong dependency between saturation concentration and temperature.

Thermal desorption methods: TD-DMA and FIGAERO
Figure 6 shows the thermograms obtained by FIGAERO and TD-DMA for three different species detected in the α-pinene ozonolysis experiment at -30 ºC.Fig. 5a shows that the C8H12O4 signal measured by both FIGAERO and TD-DMA is broad and exhibits a multimodal behavior, two maxima are observed at approximately 50-60 ºC and at ~120-150 ºC.As described in Section 2, the FIGAERO temperature is slowly ramped up to 150 ºC in approximately 15 minutes while the TD-DMA temperature increased up to 600 ºC in approximately 1 minute.From Fig. 6a it can be seen at ~250 ºC the C8H12O4 signal measured by the TD-DMA reached background levels.This might suggest that higher than 150 ºC FIGAERO temperatures are needed for completely evaporating this SOA component collected at -30 ºC.In contrast, Fig. 6b and c show that the normalized intensity of C9H14O4 and C10H16O6 first increased reached a maximum at around 40-60 ºC, and gradually decreased (sharp peak).A similar trend is observed in both FIGAERO and TD-DMA.
Additionally, we display in Fig. 6 the second TD-DMA heating.A second heating is performed immediately after the first heating, without particle collection.This performance allows us to estimate the signal coming from the particles a nd the signal coming from the background due to the inlet line.From the background measurement, it can be seen that the particle constituents measured by the TD-DMA have been efficiently evaporated.Hyttinen et al. (2022) reported the Tmax values of several particle-phase compounds measured by FIGAERO during α-pinene ozonolysis experiments.The reported values for C8H12O4, C9H14O4, and C10H16O6 are ~80 ºC, ~73 ºC, and ~61 ºC.These values do not distant from the ones measured by both TD-DMA and FIGAERO during the experiment reported here, with the exception of the second maxima observed in C8H12O4. .

Positive matrix factorization results
The results of the PMF analysis of the TD-DMA data are shown in Figure 7 which contains the factor mass spectra (Fig. 7ad), the factor thermograms (Fig. 7e), and the contribution of each factor to the total signal (Fig. 7f).We found that four factors are the best choice to reconstruct the TD-DMA data and to provide the most interpretable results (the residuals are shown in Fig. S7 in the Supplement).We numbered the factors according to their peak desorption temperatures (Fig. 7e).F1 TD-DMA which desorbs at the very first stages of the heating cycle, includes organic compounds with molecular masses between 150 and 250 Da (Mav = 211.6Da), with a carbon, hydrogen and oxygen average content (CHOac) of 9.2, 14.6 and 5.4, respectively.F1 TD-DMA contains mainly compounds in the monomer region (see Fig. S8 in the Supplement).F2 TD-DMA desorbs right after F1 TD-DMA.The mass average (Mav) is 230.7 Da and the CHOac is 10.0, 16.1 and 5.9.Compounds in the monomer region also contribute to this factor (see Fig. S8 in the Supplement).F3 TD-DMA shows a clear contribution of both monomers and dimers for the time when the time series shows a broadened peak.CHOac is 10.8, 16.7 and 5.5 and Mav = 234.5 Da.Lastly, F4 TD-DMA is dominated (∼ 60 % of the signal intensity) by a high signal with molecular mass of 172.18 Da, which corresponds to C8H12O4.This is reflected by lower values of CHOac and Mav compared to the other factors (CHOac = 8.6,12.8,and 4.5,Mav = 188.2 Da).By looking closer into F4TD-DMA (Fig. 7d), we observe that some compounds with mass < 200 Da also contribute to this factor.By integrating each factor thermogram (Fig. 7e), we calculated that F1 TD-DMA and F2 TD-DMA contribute to ~70 % of the total signal while F3 TD-DMA and F4 TD-DMA make up ~ 30 % of the total signal as shown in Figure 7f.
For the purpose of comparison, we present in Figure 8 the results from applying PMF to the FIGAERO thermal desorption data for a solution with 6 factors.Figure 8a-f contains the factor mass spectra, Fig. 8g the factor thermograms, and Fig. 8h shows the contribution of each factor to the total particle signal (the residuals and the factors expressed in terms of their oxygen content and mass are shown in Figures S9 and S10 in the Supplement, respectively).F1 FIGAERO, F2 FIGAERO and F3 FIGAERO (Fig. 8a-c) show a distinct contribution from monomers, and similar mass spectra, but display different thermal profiles (in Fig. 8g).Specifically, F2 FIGAERO and F3 FIGAERO exhibit well-defined thermal profiles (~15 ºC difference in Tmax).However, F1 FIGAERO shows a broader profile with no distinct maximum.We suspect that F1 FIGAERO can be related to some of the following causes to some extent: a) Limited resolution of the chosen PMF solution at T < 50 ºC, likely to the presence of two neighbouring factors that were not resolved completely; b) interference from volatile material already evaporating at the beginning of the thermogram and, c) adsorption of gaseous compounds.Reason b could be related to the procedure initiating the desorption, where the filter is flushed with N2 at ambient temperature before starting the heating ramp, which would likely affect the most volatile material.F4 FIGAERO and F5 FIGAERO (Fig. 8d and e As mentioned previously, we observed some factors with similar CHOac and Mav but different thermal behavior, possibly due to the presence of isomers.Molecules with the same composition but different structure and functional groups may exhibit different volatilities.In fact, functionality is one of the driving factors that determines volatility (Pankow and Asher, 2008;Wang et al., 2020).The mass spectrometry techniques reported here are not able to determine the molecular structure.We further note that the factors F4 TD-DMA and F6 FIGAERO desorb mainly at the very last stage of the heating curves, although they both have a clear contribution of compounds with mass < 200 Da, and one of lowest oxygen content of all factors (i.e., the lowest degree of oxidation).We suspect that, F4 TD-DMA (Fig. 7d) and F6 FIGAERO (Fig. 8g) are comprised primarily of products of thermal decomposition; the heat applied to desorb the particles instead cleaves certain chemical bonds in (larger) molecules before these could desorb.Small compounds are generally expected to desorb before the transmitted thermal energy (i.e., the desorption temperature) is high enough to cause such decomposition.However, F4 TD- DMA and F6 FIGAERO thermal profiles also exhibit a small peak at lower temperatures (observed more clearly in Fig. 8g at ~40-50 °C), which is likely direct desorption.This may suggest that the low temperature peak stems from monomers desorbing directly, while the broad high temperature peak represents decomposing dimers/oligomers which are then detected at the composition of the corresponding monomers.Previous studies applying PMF analysis to FIGAERO thermal desorption data also observed the presence of one or more factors dominated by thermal decomposition products for α-pinene and sesquiterpene derived SOA (Buchholz et al., 2020;Li et al., 2021).Those studies concluded that thermal decomposition was the main volatilization process at desorption temperatures above 100 °C with differences between the observed SOA types.
By applying PMF analysis to thermal desorption data we observed that often, several factors are needed to explain the behavior of a single ion.One example is shown in the thermal profile of C8H12O4 for both FIGAERO and TD-DMA (in the Supplement Fig. S11).Particularly, F4 TD-DMA and F6 FIGAERO explain the C8H12O4 signal at higher temperature.This is consistent with previous observations.Lopez -Hilfiker et al. (2015) reported a significant contribution of thermal decomposition to the detection of C8H12O4 in the α-pinene ozonolysis system and stated that small acids present in higherthan-expected concentrations in SOA are likely entirely due to thermal decomposition.In our previous work (Caudillo et al., 2021), we presented individual results of gas-and particle-phase of the same chemical system as in this study using the same ionization and detention scheme.We found that C8H12O4 contributed ~10 times more to the particles than to the gas-phase.
Furthermore, in the present study, the presence of the other PMF factors suggests that either there are at least 3 isomers with distinguishable volatility, or that there are different thermal decomposition processes occurring at different desorption temperatures which all form C 8 H 12 O 4 as a stable product.There are some studies that suggest possible decomposition pathways.For example, Hyttinen et al. (2022) investigated the two possible thermal decomposition reactions (dehydration and decarboxylation) proposed by Yang et al. (2021), and explored which reactants provide C8H12O4 as a product.They reported two isomers of C8H12O4 formed from C9H12O6 decarboxylation and one isomer formed from C8H14O5 dehydration.Certainly, the mechanisms that explain our observation remain uncertain and need to be further investigated.

Discussion of advantages and disadvantages of methods for measuring sub-100 nm SOA
Table 3 summarizes some advantages and disadvantages that should be considered when applying the methods presented here to the measurement of laboratory-generated sub-100 nm diameter SOA.When measuring particle chemical composition, the time needed for collecting enough particles (mass) should be carefully considered.This fact becomes a challenge when analyzing nanoparticles, since the small particles do not contribute significantly to the overall SOA mass.In that regard, EESI-TOF and FIGAERO provide a faster response (every 1 second and 30 minutes, respectively) than the other two methods, and allows a nearly real time monitoring.This is especially convenient when the chemical composition changes continuously (i.e., in complex environments or during oxidative flow reactor or chamber experiments).In fact, the EESI-TOF's total particle signal exhibited a good correlation with the mass concentration calculated from the Scanning Mobility Particle Sizer (SMPS) measurements (r 2 > 0.94, Fig. S12 in the Supplementary material).Despite the fact that there is a size-dependence on EESI-TOF sensitivity, EESI-TOF sensitivity decreases as the size of the particles increases (Lee et al., 2021b).The particle collection periods for TD-DMA and FILTER (offline analysis) were much longer and depended on the particle load and limit of detection.Besides the low time resolution, a main disadvantage of longer collection times is that aerosol aging may occur.This can potentially change the chemical composition and therefore lead to inaccurate aerosol speciation.Several studies have reported positive and negative artifacts caused by adsorption of gases on the collection surfaces, longer sampling periods, and volatilization of organic species either during collection or during storage (Turpin et al., 1994;Subramanian et al., 2004;Kristensen et al., 2016).
The mass spectrometers coupled to TD-DMA and FIGAERO (nitrate CI-APi-TOF and iodide HRToF-CIMS, respectively) can perform gas-phase measurements while the particle collection takes places.EESI-TOF in the dualconfiguration can measure both particle-and gas-phase quasi simultaneously.This allows for a direct comparison between gas-and particle-phase.Nevertheless, these mass spectrometry techniques can only identify chemical formulas but with some limitations.Thus, it is not possible for example to provide structural information or identification of isomers.In contrast, the UHPLC-HESI-HRMS offline method has the advantage of being able to distinguish between clusters, molecules, and isomers based on the chromatographic separation.Furthermore, the fragmentation pattern (via MS 2 -experiments) can provide hints to interpret the functional groups and can be used for unambiguous compound identification.Hence, UHPLC-HRMS can provide robust analytical insight of the stable compounds.
The thermal desorption methods (TD-DMA and FIGAERO) exhibit significant thermal decomposition of compounds with desorption temperatures above 100 °C.For the EESI this seems to be less of an issue in general but in some specific studies thermal decomposition was found to be relevant (Bell et al., 2021).PMF analysis of the thermal desorption data from the TD-DMA and FIGAERO could separate the contribution of products from thermal decomposition from those directly desorbing.However, even with this method, it is not possible to obtain information about the original compounds decomposing and their true volatility.The observed decomposition temperature can be used as an upper limit for volatility (i.e., their true volatility is lower than that associated with the apparent desorption temperature).
For the FILTER method, the compounds collected on the filter have to be extracted into a liquid phase for the UHPLC separation.The choice of solvents for this extraction will determine which fraction of organic compounds will be analyzed.The water-methanol mixture used in this study will extract polar, hydrophilic compounds similar to the watersoluble organic carbon category.Note that the exposure to water (or other solvents) may lead to chemical reactions, e.g., hydrolysis of (hydro) peroxides.The selective extraction and potential aqueous phase chemistry may explain the smaller number of compounds detected with the FILTER method.However, for the compounds that do get analyzed, a much deeper understanding can be achieved (e.g., separation of isomers).
The ionization technique also plays a role on the final detection.If the ionization technique is not soft enough, this can result in fragmentation and affects the final response in the detection.In principle, the ionization techniques utilized by the instruments reported here are soft, meaning that no significant fragmentation occurs during the detection.Nevertheless, the ionization efficiency is different between the techniques.For example, with the nitrate reagent ion, highly oxygenated species can be better detected, while the ionization techniques used for FIGAERO and EESI-TOF (I -, and Na + , respectively) are more sensitive to intermediately oxygenated organic compounds.The UHPLC-HESI-HRMS can be operated in both polarity modes, however, and therefore detect species are either able to donate protons (in the negative mode) or form clusters with protons or sodium (in the positive mode).
The complete characterization of species in the particle-phase in terms of chemical formula and structure represents an analytical challenge.In this sense, the full identification of organic compounds is only possible by combining different techniques.

Conclusions
In this study, we presented an overview of four different methods for measuring the chemical composition of ultrafine particles and we described their capabilities to detect organic compounds.Specifically, we reported the particle-phase composition from α-pinene ozonolysis at -50 ºC, -30 ºC and -10 ºC.In all the cases, the highest portion of detected compounds correspond to species with 10 and less than 10 carbon atoms (C10, C<10).The EESI-TOF generally detected a higher fraction (compared to the other techniques: TD-DMA, FIGAERO and FILTERS) of compounds with more than 10 carbon atoms (C>10).In terms of volatility classes, EESI-TOF detected the higher fraction of ULVOC in all the experiments reported here, especially for those at lower temperatures (-50 ºC and -30 ºC).We presume that several factors can explain these differences, i.e., thermal decomposition of large compounds (for the thermal desorption methods), for which we applied positive matrix factorization on the thermal profiles and suggested a 4-Factor solution for TD-DMA and a 6-Factor solution for the FIGAERO.Specifically, we suspect that F4 TD-DMA and F6 FIGAERO might be related to thermal decomposition to some extent.The PMF factors dominated by direct desorption can be interpreted as volatility classes, characterized by their Tmax values (the peak in the respective temperature desorption profiles).Nevertheless, further calibration experiments are needed to determine the relation between Tmax and saturation concentration.With the offline method UHPLC-HESI-HRMS, we were able to verify the presence of isomers (two isomers for C8H12O4, C10H16O3, C10H16O5, and three isomers for C10H16O4 and C10H16O6), which represents an important advantage over the online methods reported here.
While the methods generally agree on the most important compounds that are found in the nanoparticles, they all have their strengths and shortcomings.A major limit of these methods is that the measurements of the chemical compounds are not quantitative and only rough estimates of the exact contributions of a compound to the overall chemical composition can be made.However, knowing the limitations of each method and using combinations of the available methods can provide deeper insights into the chemical composition and volatility of nanoparticles.Table 1.Instruments for measuring particle-phase chemical composition used in the CLOUD chamber experiments.
a TD-DMA can measure both size-resolved and non-size-resolved.For this work it was chosen the non-size-resolved in order to maximize the mass collected and to intercompare with the particle-phase instruments, b in this work, only the particle-phase measurements are reported, c gas-and particle-phase can be measured by using the dual-EESI-TOF configuration.

TD
) show contribution of both monomers and dimers, with a very similar CHOac and Mav.However, they show different thermogram behavior (~20 ºC difference for Tmax, Fig 8g).F6 FIGAERO shows mainly contributions from compounds with low mass (< 200 Da) and desorbs mainly at the very end of the heating curve.The contribution of each factor to the total signal is shown in Fig 8h.

Figure 1 .
Figure 1.Experimental overview of a representative biogenic new particle formation experiment (α-pinene ozonolysis at -30 °C and 20 % RH).(a) Mixing ratio in ppbv for the precursor gases, α-pinene and ozone.(b) Particle size distribution measured by the SMPS; the color scale represents the log 10 of the normalized particle concentration per cubic centimeter (cm -3 ).The median mass diameter (MMD) is shown with a black dashed line.(c) Particle number concentration in cm -3 measured by the CPC with a cut-off diameter of 2.5 nm and mass concentration in μg m -3 (obtained by integrating the normalized mass concentration from the SMPS).(d) Particle-phase signal measured continuously by the EESI-TOF, the gray shaded area refers to the period where the EESI-TOF was averaged for the intercomparison with TD-DMA and FIGAERO.(e) Particle-phase measured by FIGAERO, the gray shaded areas refer to the particle collection period and the yellow shaded areas to the desorption period.FIGAERO measured in a semicontinuous mode, namely 15-minute particle collections followed by 15-minute desorption periods.In order to intercompare with EESI-TOF and TD-DMA, FIGAERO signals were integrated during the period where the TD-DMA collected particles (~3 hours), seven FIGAERO particle samples were integrated during the 3 hours comparison period.(f) Particle-phase measured by the TD-DMA; The TD-DMA collection period was approximately 3 hours while the desorption period lasted around 1 minute followed by a second heating for estimating the background.

Figure 2 .
Figure 2. Carbon oxidation state OSc against the number of carbon atoms for α-pinene oxidation products in the particle-phase at -30 ºC and 20 % RH measured by three different techniques.(a) TD-DMA: Thermal Desorption-Differential Mobility Analyzer coupled to a NO3 -chemical ionization-atmospheric-pressure-interface-time-of-flight mass spectrometer.(b) FIGAERO: Filter Inlet for Gases and AEROsols coupled to an I -high-resolution time-of-flight chemical-ionization mass spectrometer, and (c) EESI-TOF: Extractive Electrospray Na + Ionization time-of-flight mass spectrometer.The level of α-pinene was between 1 and 6 ppbv while the ozone level was ~100 ppbv.The carbon oxidation state is calculated as follows: OSc = 2xO:C -H:C.The marker sizes in (a), (b), and (c) represent the intensities normalized by the total signal in each system.(d) Fraction of species in the particle-phase containing less than 10 carbon atoms (C<10), 10 carbon atoms (C10), and more than 10 carbon atoms (C>10).The fraction was calculated by normalizing the intensities by the total signal in each system.

Figure 3 .
Figure 3. Carbon oxidation state OSc against the number of carbon atoms for α-pinene oxidation products in the particle-phase at -10 ºC and 80 % RH measured by three different techniques.(a) FILTER: Offline analysis of filters using Ultra-high-performance liquid chromatography (UHPLC) and heated electrospray ionization (HESI) coupled to an Orbitrap high-resolution mass spectrometer (HRMS).(b) FIGAERO: Filter Inlet for Gases and AEROsols coupled to an I -high-resolution time-of-flight chemical-ionization mass spectrometer, and (c) EESI-TOF: Extractive Electrospray Na + Ionization time-of-flight mass spectrometer.The level of α-pinene was between 1 and 3 ppbv while the ozone level was ~100 ppbv.The carbon oxidation state is calculated as follows: OSc = 2*O:C -H:C.The symbol sizes in (a), (b), and (c) represent the intensities normalized by the total signal in each system.(d) Fraction of species in the particle-phase containing less than 10 carbon atoms (C<10), 10 carbon atoms (C10), and more than 10 carbon atoms (C>10).The fraction was calculated by normalizing the intensities by the total signal in each system.

Figure 4 .
Figure 4. Distribution of volatility classes for α-pinene oxidation products in the particle-phase at -30 ºC and 20 % RH measured by three different techniques: (TD-DMA) Thermal Desorption-Differential Mobility Analyzer coupled to a NO3 -chemical ionization-atmospheric-pressure-interface-time-of-flight mass spectrometer, (FIGAERO) Filter Inlet for Gases and AEROsols coupled to an I -high-resolution time-of-flight chemical-ionization mass spectrometer, and (EESI-TOF) Extractive Electrospray Na + Ionization time-of-flight mass spectrometer.The level of α-pinene was between 1 and 6 ppbv while the ozone level was ~100 ppbv.The volatility classes (ULVOC, ELVOC, LVOC, SVOC, IVOC) were defined as in Donahue et al. (2012) and in Schervish and Donahue (2020).

Figure 5 .
Figure 5. Distribution of volatility classes for α-pinene oxidation products in the particle-phase at -10 ºC and 80 % RH measured by three different techniques: (FILTER) Offline analysis of filters using Ultra-high-performance liquid chromatography and heated electrospray ionization coupled to an Orbitrap high-resolution mass spectrometer, (FIGAERO) Filter Inlet for Gases and AEROsols coupled to an I -high-resolution time-of-flight chemical-ionization mass spectrometer, and (EESI-TOF) Extractive Electrospray Na + Ionization time-of-flight mass spectrometer.The level of α-pinene was between 1 and 3 ppbv while the ozone level was ~100 ppbv.The volatility classes (ULVOC, ELVOC, LVOC, SVOC, IVOC) were defined as in Donahue et al. (2012) and in Schervish and Donahue (2020).

Figure 6 .
Figure 6.FIGAERO and TD-DMA thermal desorption profiles for three different compounds detected in α-pinene ozonolysis experiment at -30 ºC and 20 % RH.(a) C8H12O4, (b) C9H14O4, and, (c) C10H16O6.FIGAERO and TD-DMA intensities are normalized by the reagent ions and expressed in normalized counts per second (ncps), FIGAERO signals have been divided by 1E3 in (a) and (c) and by 5E4 in (b).FIGAERO temperature is slowly ramped up to 150 ºC in approximately 15 minutes while TD-DMA temperature increased up to 600 ºC in approximately 1 minute.The TD-DMA temperature is an estimate based on the resistance of the filament.For the TD-DMA two heating profiles are needed for determining the particle signal and the background due to the heating of the inlet line.

Table 2 .
Sample size conditions during the experiments performed at -30º C, -10 ºC and -50 ºC.For the purposes of this intercomparison study, the sampling times were aligned base on the TD-DMA and FILTER collection times.