Marcel
Müller
*,
Fabrice
Stefanetti
and
Ulrich K.
Krieger
*
Institute for Atmospheric and Climate Science, ETH Zurich, Universitätstrasse 16, 8092 Zurich, Switzerland. E-mail: marcel.mueller@env.ethz.ch; ulrich.krieger@env.ethz.ch
First published on 28th November 2022
Unsaturated organic compounds in aerosol particles undergo oxidative ageing via heterogeneous reactions with atmospheric oxidants. Ozonolysis of linoleic acid has been serving as a proxy for this process and the linoleic acid decay has been shown to deviate from a linear dependence on ozone concentration. We use electrodynamic balance–mass spectrometry to measure mass spectra from levitated droplets before and after exposure to specific ozone mixing ratios. We find almost identical reactive uptake coefficients for ozone mixing ratios from 1 to 10 ppm, suggesting that a limitation due to droplet surface coverage with ozone is absent in this regime. However, a strong increase of the uptake coefficient for ozone mixing ratios below approximately 0.2 ppm is observed. Based on measurements using an oxygen atmosphere without ozone, which also show a degradation of linoleic acid, we attribute the apparent increase of the uptake coefficient to the oxidation of linoleic acid with molecular oxygen. These findings are consistent with an autoxidation mechanism, as proposed previously. Our results highlight that strongly elevated oxidant concentrations can mask the impact of autoxidation, which is slow but may dominate under atmospheric conditions.
Environmental significanceMany laboratory setups investigating the ozonolysis of model systems mimic atmospheric processes but apply strongly elevated ozone concentrations. Electrodynamic balance–mass spectrometry allows to investigate single levitated droplets, which may be trapped for several days under a wide range of oxidant concentrations, and provides molecular-level insight into particle composition. In the oxidation of linoleic acid, the influence of the reaction conditions on the decay of linoleic acid and on the product spectra reveals the importance of oxidation with molecular oxygen. Therefore, in systems with competing reactions, care needs to be taken, that elevated concentrations of one oxidant do not dominate the fate of the reactant. |
A particular chemical compound, which has been investigated in several studies as a model compound is linoleic acid (Fig. 1). This doubly unsaturated acid was used to study the uptake of ozone onto organic phases,10,11 to investigate the influence of ozonolysis on particle sizes and their cloud condensation nucleation activity,12 and to measure products from linoleic acid ozonolysis using MS techniques,13 Raman-, and infrared spectroscopy.14–17 Besides the oxidation via ozonolysis, linoleic acid is known to undergo free radical oxidation with molecular oxygen (autoxidation). The autoxidation mechanism and products of fatty acids and fatty acid esters have been studied in the fields of air pollution toxicity and food science.18–23 In particular the bis-allylic hydrogen of linoleic acid and its esters lead to a high susceptibility to hydrogen abstraction.
In studies on the oxidative ageing of linoleic acid in deposited droplets and films, it has been reported that the consumption of linoleic acid does not scale linearly with the ozone mixing ratio in the gas phase. Two effects have been suggested to cause this non-linearity. (i) An effect of particle surface coverage by ozone, which could be described with a Langmuir–Hinshelwood mechanism.15,16 (ii) A competition between oxidation by ozone and autoxidation.14,17 Because the autoxidation is expected to proceed largely independent of the ozone concentration, its relative importance increases towards low ozone concentrations.
Recent experiments have mainly focused on detection methods based on Raman- and IR spectroscopy,15–17 which are sensitive to certain functional groups but not on a molecular level. Based on these studies, we apply electrodynamic balance–mass spectrometry (EDB–MS) to investigate the influence of the ozone concentration on the decay of linoleic acid in single levitated particles. Experiments based on particle levitation in an EDB allow to store particles for long durations of several days without contact between the particles and a support. Therefore, the investigation and the comparison of atmospheric aerosol processes that require a wide range of reaction durations (e.g. due to strongly differing oxidant concentrations) is rendered possible.24 By coupling an EDB with MS, molecular-level information about the droplet composition can be obtained during the ageing process.
In brief, single droplets of linoleic acid were produced from a linoleic acid solution (10 wt% linoleic acid (≥99%, Sigma-Aldrich) in methanol (UHPLC for MS, Sigma-Aldrich), freshly prepared or stored at <6 °C for less than one week) using a droplet on demand generator. Droplets were charged inductively and injected into a quadrupole type EDB. In the EDB, droplets are confined by the electrodynamic field from AC and DC electrodes, with the DC voltages compensating the gravitational force of the droplets. Once injected, the droplets are stored for at least 15 minutes in a nitrogen atmosphere before being exposed to oxidants. During this time, all methanol evaporates from the droplet, which is indicated by reaching a steady voltage required to counteract gravitation.
The two-dimensional angular scattering pattern (phase function) of particles illuminated by a laser is captured and the mean peak-to-peak angle in the phase function is used to estimate the droplet size by comparison to simulated phase functions.36,37 The mean droplet size in this study was 25.4 ± 2.5 μm. The obtained MS signal was normalised with respect to the initial droplet volume.
Droplet ageing experiments were carried out by flowing a dry gas stream of 300 sccm through the EDB (RH < 15%). For measurements at ozone mixing ratios ≤1 ppm, a calibration ozone generator (Photometric O3 Calibrator – Model 401, Advanced Pollution Instrumentation, USA) was used, otherwise an ambient ozone simulator (AOS 2 with bypassed pump, BMT Messtechnik, Germany) was used. Both ozone generators use UV-light for generating ozone and, hence, we expect no significant NOx production during ozone generation. However, while the generator used for low ozone mixing ratios uses a built-in zero air generator, the one used for higher mixing ratios uses the untreated pressurised air supply. Oxidation experiments without ozone were carried out with air from the same supply. For measurements of the ozone mixing ratio, an electrochemical ozone sensor was used. We measured the concentrations of NOx in the air before and after ozone production using a chemiluminescence NOx sensor (Model 200a, Advanced Pollution Instrumentation, USA). For all experiments with ozone mixing ratios ≤ 1 ppm, NOx levels were below detection limit (1 ppb), whereas at all other measurements the NOx levels were those of outside air (typically below 20 ppb). Ozone levels in the in-house supplied air were below detection limits (1 ppb, measured with Ozone Monitor BMT 932).
For the analysis of droplet compositions, droplets were ejected from the trap into a home-built evaporation unit, held at 190 °C. The total gas flow to the mass spectrometer was 1.5 l min−1, consisting of a gas flow through the evaporation unit (300 sccm N2), a control gas flow (5 sccm N2 guided through a gas washing bottle filled with a mixture of 1% toluene in ethanol) and complemented by a nitrogen flow. Subsequent to the droplet injection into the evaporation unit, the evaporation unit was flushed for 8 min, during which the MS signal was recorded. A cold plasma dielectric barrier discharge ion source (SICRIT SC-30X Ionisation Set, Plasmion GmbH, Germany) was used for ionisation of gas-phase species.38–40 The analysis by MS was carried out with a triple quadrupole mass spectrometer (QTRAP 4500, AB Sciex LLC, USA) at 2000 Da s−1 and unit resolution. Preliminary measurements of a broad range of mass-to-charge ratios (m/z) were carried out in negative Q3 mode from m/z 40 to 340 (total scan time of 0.155 s). Measurements focusing on the quantification of linoleic acid measured a limited set of m/z ranges to achieve a higher temporal resolution during the evaporation of droplets. The probed ranges were at m/z 84–86, 100–102, 120–122, 142–144, 156–160, 169–173, 184–188, 196–204, 208–212, and 276–283 with a pause between mass ranges of 5 ms (total scan time of 0.070 s).
Tandem MS experiments were carried out in negative mode with a scan rate of 2000 Da s−1 from m/z 5 to 212 (0.11 s total scan time). The collision cell was operated with N2, the collision energy was set to −20 V, and the collision cell exit potential to −11 V. The tandem MS signal was summed over a duration of 10 s after droplet evaporation started and assigned to bins of integer m/z.
Further data processing was carried out with MATLAB (Version R2019a, MathWorks, USA). For the preliminary comparison of entire mass spectra, the droplet signal was integrated over a duration of 240 s from the start of the evaporation, corrected for the background (measured during the last 30 s of this duration), and corrected for the signal obtained from the gas flows without droplet ejection. For the analysis of quantification experiments, an integration duration of 480 s was used.
Besides the peak at m/z 279, there are intense peaks at m/z 278 and 280. The signal at m/z 280 is assigned to the isotope signal of linoleic acid with one 13C atom, the peak at m/z 278 may stem from a deprotonated radical species (e.g. the pentadienyl radical), probably formed during ionisation. Although we did not see evidence for radical species in previous experiments, charged radical species are known to be involved in dielectric barrier discharge ionisation.40 A strong contribution to the spectrum at an even m/z is not expected for non-radical ions consisting only of carbon, oxygen and hydrogen atoms due to their valence and atomic mass.
Further minor peaks were detected at m/z 293, 295, 311, and 325. These stem potentially from impurities in the injected solution or from adducts formed in the ionisation process. As the intensity of these peaks goes down with increased exposure to ozone, it seems unlikely that they are due to the presence of linoleic acid oxidation products.
Pure, unreacted linoleic acid droplets of different sizes were used to estimate the process sensitivity of evaporation, ionisation and MS analysis, and to ensure a linear relationship between the volume and the MS signal of the droplets. In Fig. 7 in the appendix, the signal from 17 droplets collected on eight different days over the course of four months is plotted against the estimated droplet volume and reveals a linear trend within the uncertainties of the measurements. The signal variability from ionisation and mass spectrometry was estimated from averaged measurements of the control-gas species (ethanol and toluene), and nitrate generated in the ion source, that were taken on every day of MS measurements. The strongest variability of approximately 11% (standard deviation relative to the mean value) was obtained for deprotonated ethanol at m/z 45 (7.2 ± 0.8 × 105 cps). This variability was considerably lower than the estimated uncertainty in particle volumes of about 30% (from a conservative estimation of a 10% uncertainty in radius measurements).
(1) |
(2) |
With the method presented in this paper, it is possible to combine the measurements of several particles that were aged for different reaction times to reconstruct a decay curve for the amount of linoleic acid in the droplets. This was done for five different ozone mixing ratios (0.03, 0.2, 1, 5 and 10 ppm). As can be seen in Fig. 8 in the appendix, the consumption of linoleic acid can be described well with an exponential decay. Exponential curves were fitted to the data taking into account the uncertainty in particle size and calculating a common [LA]0 for all fits using Bayesian inference as described by Müller et al.35
Eqn (2) suggests that the decay of linoleic acid depends on the product of the ozone concentration in the droplet and time ([O3] × t). It is a common assumption in aerosol ageing experiments that the reaction progress and the outcome of aerosol ageing can be described with the ozone exposure (product of applied ozone concentration and reaction time). Note that the exposure metric is equivalent to [O3] × t, if [O3] scales with the gas-phase concentration of ozone. This allows to compensate short reaction times in laboratory experiments with high ozone mixing ratios to reach atmospherically relevant ozone exposures. There are, however, a few cases (including linoleic acid oxidation) where this assumption breaks down14,17,35,41–43 and therefore the extrapolation from a laboratory setting to atmospheric condition needs to be assessed carefully in each case.
A compilation of all data from experiments with different ozone mixing ratios is plotted against exposure in Fig. 3. For high ozone mixing ratios, the decay curves lie almost perfectly on top of each other. However, the curves seem to deviate from the common curve for ozone mixing ratios lower than 1 ppm. Apparently, the degradation of linoleic acid is more efficient at low ozone concentrations.
Fig. 3 Linoleic acid decay with ozone exposure for different ozone mixing ratios. Each data point corresponds to an individual droplet that has been exposed to a specific gas phase for a certain time. Uncertainties on the single data points are depicted in Fig. 8 in the appendix and are left out here for clarity. The axis on the right gives the corresponding number of molecules based on an initial droplet of 25 μm radius. Exponential curves are based on point estimates from exponential fits taking into account the uncertainty in droplet size. |
Fig. 4 Linoleic acid decay in droplets exposed to air. Error bars of aged droplets depict the uncertainty from droplet volume estimation. |
In contrast, the composition of droplets stored in a nitrogen atmosphere for comparable times did not show any measurable deviation in the composition in comparison to unreacted droplets, nor any measurable size change.
The spectrum of a droplet exposed to air for five days is shown in Fig. 2c. Because of the different experiment durations, a direct comparison of spectra from different oxidation conditions (Fig. 2) is complicated. It seems, however, that the spectra from the two different oxidants (Fig. 2b and c) are rather similar. Both spectra are dominated by a peak at m/z 171 and many of the intensive peaks are found in both spectra. An interesting feature is the peak at m/z 211 in the ozonolysis experiment, which is seemingly absent in the autoxidation experiment. In turn, the signal from m/z 213 is much more intensive after exposure to ozone. Ongoing decrease in mass and size (not shown) indicates evaporation of semi-volatile compounds even after all linoleic acid was consumed, potentially due to ongoing secondary reactions.
Also, there is evidence that the pentadiene substructure may undergo radical-induced hydrogen abstraction to produce the relatively stable pentadienyl radical, which can participate in further reactions (c).14,17,19–22 Note that hydroperoxyl radicals can act as hydrogen abstractors and therefore maintain an autocatalytic cycle. As the double bonds are preserved in the autoxidation, the resulting hydroperoxides may undergo ozonolysis as well, although reaction rates may drop for conjugated systems.50
When the decay of linoleic acid follows an exponential decay, the reactive uptake of ozone on linoleic acid (γ, ratio of ozone molecules hitting the surface divided by the number of linoleic acid molecules consumed) may be calculated from
(3) |
With the incoming flux of ozone molecules on the surface of a droplet Fin, the number concentration of molecules in the gas phase nO3, the mean thermal velocity and the surface to volume ratio A/V.15,16
For the initial reactive uptake coefficient on droplets, eqn (1) and (3) may be combined with A/V = 3/r for spherical geometries to give
(4) |
The initial reactive uptake coefficient was calculated for all ozone mixing ratios and the obtained coefficients are shown in dependence on the ozone mixing ratios in Fig. 6, together with literature values.
Fig. 6 Reactive uptake coefficients for different ozone mixing ratios. Open symbols denote experiments that were run under reduced pressure (the x axis on the top does not apply to these). Note that Moise and Rudich10 and Thornberry and Abbatt11 measured only the loss rate of ozone, and Broekhuizen12 investigated droplet size changes to derive an uptake coefficient, while all others derived the uptake coefficient from the loss rate of linoleic acid. Literature data were taken from Moise and Rudich,10 Thornberry and Abbatt,11 Broekhuizen,12 Hearn and Smith,13 Zeng et al.15 and He et al.16 Apparent reaction rates from Chu et al.17 were converted to reactive uptake coefficients with eqn (4) using r = 70 μm. |
Our derived reactive uptake coefficients show a relatively constant value for ozone mixing ratios ≥ 1 ppm and an increase towards lower mixing ratios. There is a considerable spread in the reported reactive uptake coefficients in the literature. Our data for ozone mixing ratios > 0.2 ppm are rather on the higher side but generally in agreement. For low ozone mixing ratios, the literature data are sparse but the general trend of a steep increase of reactive uptake coefficients towards mixing ratios below approximately 0.2 ppm ozone (logarithmic abscissa in Fig. 6) is emerging from other studies as well.15,17 In the studies by Moise and Rudich10 and Thornberry and Abbatt,11 only the loss rate of ozone was monitored, which may explain the low uptake coefficient, as autoxidation is not detected (see Sect. 4.2).
In the recent literature, a Langmuir–Hinshelwood surface saturation behaviour has been discussed for the surface adsorption of ozone on linoleic acid.15–17 Such a saturation behaviour would lead to a significant decrease of the reactive uptake coefficient with increasing ozone mixing ratios. However, our reactive uptake coefficient level off towards higher ozone mixing ratios and stay at a constant value around 1.5 × 10−3. Therefore, our data suggest that for levitated droplets with radii of approximately 25 μm there is no significant surface saturation effect at ozone mixing ratios below 10 ppm. Surface saturation was reported in studies of films or deposited droplets of markedly bigger sizes (50 μm and 60–80 μm radius). Those have lower surface-to-volume ratios, which may have enhanced a surface saturation effect.
Besides the discussion of surface saturation, there are indications of autoxidation processes from Raman spectra of levitated and deposited droplets.14,17 Both reactions, the ozonolysis and the autoxidation can lead to the consumption of linoleic acid. However, the ozonolysis is expected to scale linearly with the gas-phase ozone concentration (in the absence of surface saturation phenomena), whereas the autoxidation should only depend on the availability of oxygen and a radical chain initiation but proceed independently from the applied ozone concentration. This means that the relative importance of the autoxidation increases towards low-ozone conditions, which explains the seemingly higher efficiency to consume linoleic acid (Fig. 3) and the increase in the reactive uptake coefficient (Fig. 6) towards lower mixing ratios. Note that in our experiments the oxygen concentration in air is not significantly changing and can be assumed identical in all experiments.
We obtain here an estimate of the potential autoxidation rate from the linoleic acid decay in droplets exposed to air, in the absence of ozone. Because of the potential acceleration of the autoxidation by ozonolysis products, we use the high autoxidation rates from the radical chain reaction period for the comparison.
The e-folding lifetime of linoleic acid in the radical chain reaction period (t ≥ 40 h) is 29 ± 7 h. This is equal to a pseudo-first order rate constant kautox of 0.034 ± 0.007 h−1. We use this value as an estimate for the loss rate of linoleic acid due to autoxidation for droplets of 25 μm radius at an oxygen mixing ratio of 21%, as applied in all our experiments.
The loss of linoleic acid that is exclusively due to ozonolysis is derived from experiments at high ozone mixing ratios, where it is assumed that the decay of linoleic acid is dominated by ozonolysis because the reactive uptake coefficient for ozone remain largely constant (γ ≈ 1.5 × 10−3).
At low ozone mixing ratios, this reactive uptake coefficient may now be combined with the autoxidation rate kautox. We investigate here the conditions in our ozonolysis experiment with the lowest ozone mixing ratio (0.03 ppm). For two independent decay routes in this regime, the expected decay may be calculated in analogy to eqn (3) using
(5) |
For 0.03 ppm ozone, using the initial uptake coefficient from high ozone concentrations and the initial linoleic acid concentration, eqn (5) yields an exponential lifetime of 15.9 h, which is within the error of the experimentally inferred e-folding lifetime of 18 ± 5 h. This indicates that the apparent increase in reactive uptake of ozone on linoleic acid at low mixing ratios may be explained with autoxidative loss of linoleic acid. Also, it implies that the autoxidation induction period is shortened because of the ozone chemistry, and that the autoxidation rate reaches almost the same magnitude as in the radical chain reaction period of droplets aged in air.
It has been proposed in the literature that the induction period in autoxidation reactions may be influenced by byproducts from ozonolysis reactions.21,51 For example, Criegee intermediates have been found to lead to OH production in gas-phase reactions,53,54 and similar reactions have been proposed to occur in the particle phase.55
The uptake coefficient derived by Moise and Rudich10 from measurements of the ozone loss represent the uptake of ozone since their experiment did not detect the linoleic acid decay. Their obtained value should therefore be compared with the uptake coefficients for ozone-dominated linoleic acid decay. Indeed, it agrees well with our uptake coefficients for ozone mixing ratios ≥ 1 ppm.
The product spectrum from ozonolysis (Fig. 2b) is overall similar to the spectrum from autoxidation (Fig. 2c). The similarities in the spectra from air and ozone oxidation suggest that products with a similar thermal decay pattern dominate the droplet composition in both experiments. This indicates that autoxidation-caused isomerisation is not dominating but autoxidation results in products with a similar carbon backbone as ozonolysis. Therefore, when using thermal evaporation and unit mass resolution MS, the product evolution is only of limited use to distinguish between autoxidation and ozonolysis.
An interpretation of the product spectra and a discussion of the evolution of the product peaks with exposure at different mixing ratios can be found in the appendix.
When comparing the product spectra from ozonolysis (Fig. 2b) and autoxidation (Fig. 2c), it seems that in the experiment with air, peaks at m/z 187 and 213 are gaining slightly in strength (relative to m/z 171) whereas the signal at 115, 197 and 211 seem less pronounced. However, note that the relative intensities may be influenced by evaporation of semi-volatile compounds as the reaction times for these two experiments differ strongly. While under autoxidation conditions the linoleic acid peak decreases, the signal from m/z 295 remains approximately constant.
Analogous to the assessment of the decay of linoleic acid, the product formation can be compared at different ozone mixing ratios. A selection of typical product formation curves is shown in Fig. 9. Generally, there are MS peaks, for which the signal intensity increases rather monotonous with increasing ozone exposure (e.g. m/z 143, 157, 171, 173), and peaks, for which the signal intensity increases at first and later decreases again (e.g. m/z 197, 199, 209, 211). An increase followed by a decrease in concentration could be explained by evaporation of semi-volatile compounds or by further secondary reactions. However, a strong contribution of evaporation to the decrease can be ruled out since this should be much more pronounced in long experiments under low ozone concentrations. It is therefore most likely that the main contributor to the strong signal at these m/z is formed in the reaction of linoleic acid with ozone but reacts in a subsequent reaction and decreases in concentration. A compound, which might explain this behaviour is 12-oxodec-9-enoic acid, one of the two primary formed aldehydes with low vapour pressure (see above). Because it still contains one double bond, it can react further with ozone, which causes the decay. Note that the corresponding signal (at m/z 211) does not drop to zero, even at a very high ozone exposure of 30 ppm h. To exclude contributions to this signal from different species with the same nominal mass, tandem MS spectra were measured from droplets at exposures when the highest signal is expected and at the maximum exposure shown in Fig. 9. The resulting spectra (shown in Fig. 10) are almost identical, indicating that the dominant contribution to the peak at m/z 211 is the same in early and late stages of the experiment. Therefore, it seems plausible that a high particle viscosity and (partial) phase change inhibit the reaction through low ozone- and reactant diffusion. A phase change at high exposures was reported in the literature14 and is observed in our experiments for exposures of approximately 4–5 ppm h as significant distortions in the measured phase function.
There appears to be no significant dependence on ozone concentration when scaling the measured product signal with exposure. This can be attributed to the similarities of the product spectra from autoxidation and ozonolysis. Also, ozonolysis is expected to strongly dominate in all experiments with ozone mixing ratios ≥ 1 ppm. Strong effects on the droplet composition from product evaporation seem to be absent.
This journal is © The Royal Society of Chemistry 2023 |