Ben Wodena,
Yizhou Sub,
Maximilian W. A. Skoda
c,
Adam Milsom
b and
Christian Pfrang
*bd
aDepartment of Chemistry, University of Reading, Reading, UK
bSchool of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, UK. E-mail: c.pfrang@bham.ac.uk
cISIS Neutron and Muon Source, Didcot, UK
dDepartment of Meteorology, University of Reading, Reading, UK
First published on 3rd December 2024
Aqueous aerosols are often covered in thin films of surface-active species, such as fatty acids which are prominent components of both sea spray and cooking emissions. The focus of our study is one-molecule thin layers of linoleic acid (LOA) and their behaviours when exposed to ozone in multi-component films at the air–water interface. LOA’s two double bonds allow for ozone-initiated autoxidation, a radical self-oxidation process, as well as traditional ozonolysis. Neutron reflectometry was employed as a highly sensitive technique to follow the kinetics of these films in real time in a temperature-controlled environment. We oxidised deuterated LOA (d-LOA) as a monolayer, and in mixed two-component films with either oleic acid (h-OA) or its methyl ester, methyl oleate (h-MO), at room temperature and atmospherically more realistic temperatures of 3 ± 1 °C. We found that the temperature change did not notably affect the reaction rate (ranging from 1.9 to 2.5 × 10−10 cm2 s−1) which was similar to that of pure OA. We also measured the rate coefficient for d-OA/h-LOA to be 2.0 ± 0.4 × 10−10 cm2 s−1. Kinetic multi-layer modelling using our Multilayer-Py package was subsequently carried out for further insight. Neither the change in temperature nor the introduction of co-deposited film components alongside d-LOA consistently affected the oxidation rates, but the deviation from a single process decay behaviour (indicative of autoxidation) at 98 ppb is clearest for pure d-LOA, weaker for h-MO mixtures and weakest for h-OA mixtures. As atmospheric surfactants will be present in complex, multi-component mixtures, it is important to understand the reasons for these different behaviours even in two-component mixtures of closely related species. The rates we found were fast compared to those reported earlier. Our work demonstrates clearly that it is essential to employ atmospherically realistic ozone levels as well as multi-component mixtures especially to understand LOA behaviour at low O3 in the atmosphere. While the temperature change did not play a crucial role for the kinetics, residue formation may be affected, potentially impacting on the persistence of the organic character at the surface of aqueous droplets with a wide range of atmospheric implications.
The oxidation of these organic species to form secondary organic aerosols (SOAs) has been studied extensively in the past,9–12 and the SOA thereby produced will have climatic cloud lifetime effects. However, the partitioning of these species into monolayers at the surfaces of cloud droplets has also been observed to affect their reactivity.4–6,13 Perhaps more significant, though, may be the effects that the presence of a surface monolayer has on the cloud droplets themselves. Cloud formation, growth, evaporation, and rainout processes have a critical dependence upon surface tension as a result of the Kelvin effect.14 The presence of insoluble organic films at the air–water interface has a lowering effect on the surface tension of the resultant water aerosols, with impacts on interactions with clouds, and with cloud lifetime effects.3 The lowering of the surface tension associated with the presence of a monolayer would seem, via a weakening of the Kelvin effect, to suppress the destabilisation of small droplets and resultant suppression of nucleation that the Kelvin effect produces.
On this basis, the presence of these surfactants on aerosols could promote nucleation via reducing the surface tension of the droplets thereby nucleated. Additionally, the presence of such monolayers on the surface of existing droplets enhances their stability at lower diameters. This enhancement of nucleation and comparative stabilisation of lower diameter droplets could lead to increased cloud lifetimes.
There are also several other routes via which the presence of a monolayer film on the surface of a cloud droplet has important consequences for atmospheric chemistry. On the one hand, such films can act as a barrier to mass transport between the water droplet and the atmosphere. Many studies have shown that insoluble films at the interface are able to significantly retard both the evaporation of the water droplet itself and the transport of oxygen (O2), ozone (O3), ammonia (NH3), and many other important atmospheric species between the water droplet bulk phase and the atmospheric bulk phase.15–22 This inhibition of cross-interfacial transport could hinder the progress of atmospheric reactions inside aqueous aerosols that rely heavily on mass-transport steps across the air–water boundary.23
On the other hand, the surface film itself can act as a ‘2D solvent’ for species from the atmosphere, thus allowing the dissolution of species that would not normally dissolve in a water droplet, or altering the solvation behaviour of species that normally would dissolve in the aqueous droplet.24–28 This effect may play an important role in the transport of various species by water droplets, as it allows species that would not normally be transported by water droplets to be adsorbed at the coated interface.3 These changes in solvation behaviour at the interface could impact the rates and even nature of atmospheric reactions that occur heterogeneously at the surface of such particles. Studies of some reactions have shown a reduction in rates, such as the rate of production of nitric acid (HNO3) from dinitrogen pentoxide (N2O5) and water29,30 as a result of the presence of these coatings3 which is a key atmospheric process for the redistribution of nitrogen amongst various species in the atmosphere.
In summary, the reactivity and properties of these surfactant films deserve further study, as they likely have a vast range of effects that are of climatic and meteorological significance, be that by affecting the reactivity of the film components and associated SOA production, by weakening the Kelvin effect and thus extending the lifetime of reflective clouds, or by modifying key atmospheric chemical processes that involve adsorption to the surface of or diffusion into the bulk of water droplets.
The by far most commonly studied unsaturated fatty acid in the context of atmospheric aerosol coatings is oleic acid (OA).5,6,31–34 OA will react with highly reactive gas-phase species in the atmosphere such as O3,5,32,33,35 OH34 and NO3.6,31 OA’s single double bond limits the complexity of potential products being formed making it a fairly well-understood model system for ozone interaction with atmospheric surfactants.36 The poly-unsaturated analogue of OA, linoleic acid (LOA), represents a logical next step for investigations.
The focus of the study presented here is one-molecule thin layers of linoleic acid (LOA) and their behaviours when exposed to ozone in multi-component films at the air–water interface. LOA contains two double bonds allowing for ozone-initiated autoxidation, a radical self-oxidation process, as well as traditional ozonolysis.37,38 The core of the experiment involved studying the oxidation of custom-deuterated LOA (d-LOA) as a pure monolayer, and in mixed monolayers with oleic acid (h-OA) and methyl oleate (h-MO; the methyl ester analogue of OA) at both room temperature and a more atmospherically relevant temperature of 3 ± 1 °C. The OA/LOA system was also studied in a reverse deuteration configuration (d-OA/h-LOA) in both these temperature conditions. Exploratory data from d-LOA/stearic acid (h-SA; the saturated analogue of OA) mixed monolayers was also gathered. The custom-deuterated surfactants employed in the present study are displayed in Fig. 1.
![]() | ||
Fig. 1 Custom-deuterated surfactant molecules employed in the present study: (a) deuterated linoleic acid (d-LOA); and (b) deuterated oleic acid (d-OA). |
Compared to OA reactivity, there has been fairly little previous work on the oxidation of LOA by ozone: He et al.39 studied a different LOA morphology and used much higher ozone concentrations than found in the atmosphere; Chu et al.37 suggested that competing LOA reaction mechanisms dominate under different ozone conditions.
A custom-built aluminium gas flow cell (volume of ca. 1.5 L) was fitted with a PTFE liquid trough (inner dimensions of 238 mm × 70 mm). This reaction chamber was mounted on the sample stage and interfaced with the gas delivery system. The trough was filled with 90 mL of null-reflecting water (NRW). The height of the air–liquid interface was aligned with respect to the neutron beam using a Keyence laser displacement sensor (model no. LK-G402), which was coupled into the sample chamber via a quartz window to allow automated height adjustment during the measurements. Height adjustments over a ca. 2 hours experiment were always less than 0.15 mm for a water height of ca. 5 mm in the trough. The chamber was designed to provide a controlled and confined environment in which monolayers can be oxidised by a gas-phase oxidant while under analysis by both neutron reflectometry and infra-red reflection absorption spectroscopy (IRRAS) as described in detail in a method paper.40 For the present study, a further iteration of the sample environment development was deployed to allow for cooling of the subphase, in order to access relevant atmospheric temperature conditions (compare Woden et al., 2021).33 Because of the low relative humidity (RH) used in the work presented here, we did not observe any condensation within the reaction chamber or on the windows throughout the experiments even at near-freezing temperatures (3 ± 1 °C).
Monolayers were spread using 20–40 μL of the spreading solutions in chloroform, leaving a monolayer of the dissolved species after solvent evaporation. Dry oxygen was continuously flowing (flow rate: 1.2 L min−1) into the chamber to provide a low (<10%) RH environment and to avoid build-up of any gas-phase products. Data were recorded for several minutes before ozone (O3) was admitted into the chamber. O3 was generated by exposing the O2 flow to UV light using a commercial PenRay ozoniser (UVP Ltd, Cambridge) to ozonise the stream of O2 (99.999%; BOC Ltd) regulated by an electronic mass flow controller to achieve O3 mixing ratios in the range of 98–983 ppb; the ozoniser was calibrated offline using UV-vis absorption at 254 nm and an absorption cross-section value of 1.13 × 10−17 cm−2 (see Daumont et al., 1992).42 We were working in large excess of O3 compared to the organic monolayer, and [O3] remained approximately constant during the reaction.
Neutron reflectivity (the fraction of incident neutrons reflected, R) varies as a function of the energy and reflection angle of the incident neutrons (expressed as momentum transfer, Q) and the scattering length density, SLD (ρ) and thickness (τ) of the monolayer as shown in eqn (1) (based on Lu et al., 2000):43
![]() | (1) |
The SLD and layer thickness can be inferred from the relationship between reflectivity and momentum transfer measured by the instrument. The two parameters are fitted over the whole Q range as a combined ρτ value, which corresponds to a surface concentration of scattering length, from which the surface concentration of the surfactant of interest can be determined if the scattering length is known (e.g. b = 315 fm for deuterated oleic acid). Reflectivity curves of R vs. Q were fitted using MOTOFIT44 to give ρτ values for each 20 s time slice. These fitted ρτ values can be converted into surface concentration (Γ) values for the surfactant following eqn (2):43
![]() | (2) |
The error bars for the SLD × thickness product were propagated using the calculated fitting errors for each parameter obtained from MOTOFIT.44 In order to achieve a sufficient contrast of the monolayer compared to the other phases in the experimental system, a deuterated form of oleic acid (d-OA), was used (Sigma-Aldrich at 98% atom D; 99.9%). Spreading onto an aqueous subphase will cause the acidic deuterium to be exchanged with the subphase, so we use the scattering length of oleic acid with 33 deuterium atoms (315 fm) to calculate the surface concentration from ρτ values. For custom-deuterated linoleic acid (d-LOA; Sigma-Aldrich at 98% atom D; 98%), the equivalent consideration leads to the presence of 31 deuterium atoms in the fully deuterated LOA.
Multilayer-Py45 provides a framework for constructing kinetic multi-layer models, so that the model code is produced automatically for the user, removing potential human error in typing out the model code. This code is written in a readable format, enabling the code to be shared easily and facilitating more reproducible modelling results. This is further supported through the Jupyter notebook,46 which is a document that incorporates both Python code and markdown text and is becoming an increasingly popular way of sharing and describing scientific code. The usefulness of Multilayer-Py has been demonstrated by application to the ozonolysis of OA45 employing both KM-SUB47 and KM-GAP48 modelling approaches.
A detailed description of the KM-SUB model concept used in the present study has been presented previously47 and is not repeated here. Essentially, the model splits our experimental system into a number of layers. The diffusion of reactants between each layer and the reaction of each component within each layer are resolved. Surface chemistry and the adsorption and desorption of gaseous species are also resolved. For this specific study, we introduced four different treatments of the oxidation of LOA by ozone: (i) single step oxidation considering the gradual build-up of ozone in our reaction chamber; (ii) the same assumptions as in (i), but now also assuming an inert residue remaining; (iii) two-step oxidation with a reaction product reacting with the reactant fatty acid in the second step; and (iv) the same assumptions as in (iii), but now also assuming an inert residue remaining at the interface.
Pseudo-first-order rate coefficients were fitted to these neutron reflectometry time series displayed in Fig. 2 using the stretched exponential analytical model outlined in the ESI Section S3 in Woden et al. (2021).33 Fig. 3 displays a second order plot for this reaction (omitting data from the highest ozone concentration, as these reaction conditions were judged to be too fast to be fitted, given the limitation of a 20 seconds time resolution in our experimental approach). The second order rate coefficient obtained from these data is 2.0 ± 0.4 × 10−10 cm2 s−1. This is slightly lower than the values we measured for pure OA ((2.2 ± 0.4) × 10−10 cm2 s−1 at 21 ± 1 °C and (2.2 ± 0.2) × 10−10 cm2 s−1 at 2 ± 1 °C),33 but the error bars overlap significantly.
![]() | ||
Fig. 3 Second order rate plot for the oxidation of d-OA in a mixed monolayer with h-LOA based on the data displayed in Fig. 2 (omitting the data for [O3] = 983 ppb). |
A further experiment was performed collecting neutrons over the full Q range accessible to INTER, in order to quantify any residue left behind after ozonolysis (as we have reported for low temperature ozonolysis of a pure OA monolayer; see Woden et al., 2021).33 Fig. 4 shows reflectivity curves before and after oxidation at 3 ± 1 °C where a residue can clearly be identified. Fitting monolayer parameters to these reflectivity curves to quantify the absolute amount of deuterated material adsorbed at the interface determined that 13% of the material present before ozonolysis remained afterward. This is similar to the range seen for pure OA monolayers at these temperatures (up to 11.1% deuterated residue was reported in Woden et al., 2021).33 Fig. 4 shows that the fit to the data collected before reaction deviates at high Q values. This is because the reflectivity background is not always Q-independent, as assumed in the MOTOFIT model and this becomes important at high Q where there is little signal compared to the background. This can be avoided by using a fixed background parameter derived from a clean air-NRW measurement, and this is the procedure used for fitting most of the data presented. This does, however, assume that the background will not change across the conditions studied (including temperature variations). Therefore, when carrying out experiments specifically to quantify post-oxidation residue under various conditions, in which distinguishing genuine residue from background signal is particularly important, this fixed-background method was avoided.
Due to time constraints during neutron beamtime, no analogous experiment could be performed at room temperature. However, the kinetic data shown in Fig. 2(a) and (c) – collected across a restricted Q range and therefore not suitable for precise absolute quantification of small amounts of material – suggests that a similar residual monolayer remained.
Fig. 8(a) to (c) overlay the temperature conditions, displaying only the highest and lowest ozone concentrations for visual clarity.
![]() | ||
Fig. 8 Time series of neutron reflectometry data contrasting the two temperature conditions for the highest and lowest [O3] for the oxidation (see Fig. 5–7 for individual plots for all [O3] studied) of (a) a pure d-LOA monolayer; (b) a mixed monolayer of d-LOA with h-OA; and (c) a mixed monolayer of d-LOA with h-MO. |
Fig. 9 displays a second-order plot for the LOA oxidation, and Table 1 summarises the fitted second-order rate coefficients for the three monolayer types in the two temperature conditions, as well as the rate coefficients for each monolayer type (treating temperature as irrelevant), for each temperature (treating monolayer type as irrelevant), and for all data combined. Table 1 demonstrates that neither the change in temperature nor the introduction of a co-deposited film component alongside d-LOA consistently affected the rate of reaction of LOA with ozone.
![]() | ||
Fig. 9 Second-order rate plot for the data displayed in Fig. 5–8. |
Monolayer | Rate coefficient/10−10 cm2 s−1 | ||
---|---|---|---|
21 ± 1 °C | 3 ± 1 °C | Combined | |
d-LOA | 2.1 ± 0.7 | 2.0 ± 0.7 | 2.0 ± 0.4 |
d-LOA/h-OA | 2.0 ± 0.7 | 2.5 ± 0.7 | 2.3 ± 0.4 |
d-LOA/h-MO | 1.9 ± 0.7 | 2.5 ± 0.7 | 2.0 ± 0.4 |
Combined data | 2.0 ± 0.2 | 2.3 ± 0.3 | 2.1 ± 0.2 |
The rates displayed in Table 1 are higher compared to those reported by He et al.39 However, our study uses a very different morphology (monolayer at the air–water interface) compared to that used by He et al. Furthermore, He et al. use a much higher ozone concentration than our study (ca. 10 ppm compared to ca. 100 to 1000 ppb in the present study) and Chu et al. have specifically warned that, due to the competing mechanisms for this reaction that dominate under different conditions, extrapolating from high ozone concentrations downwards is likely to be problematic.37 Other studies on this heterogenous reaction have mostly reported uptake coefficients, rather than rate coefficients, and these have varied by around an order of magnitude as reviewed by He et al.39 The use of an uptake coefficient, which is more dependent on reaction conditions and geometry than is a rate coefficient, makes these studies less useful as a guide for what to expect in our study, which uses a monolayer at the air–water interface rather than particulate phase or film-coated flow tube setups. Differences in mechanism may also be driving some of these discrepancies.
An experiment analogous to that described earlier for the d-OA/h-LOA system (compare Fig. 4) was performed in order to quantify any post-oxidation residue. Again, only the low temperature was studied due to time constraints at neutron beamtime experiments. Fig. 10(a)–(c) display reflectivity curves before and after oxidation at 3 ± 1 °C for d-LOA, d-LOA/h-OA and d-LOA/h-MO, respectively.
For d-LOA as a pure monolayer (Fig. 10(a)), a residue is clearly observed, and fitting the monolayer parameters quantifies this as ca. 7% of the initial adsorbed deuterated material. For d-LOA in a mixed monolayer with h-OA, a residue is also observed and quantified as ca. 11% of the initial adsorbed deuterated material. The data for d-LOA in a mixed monolayer with h-MO are slightly more difficult to interpret. The fitting process used throughout this work does successfully converge and fit a curve (displayed in Fig. 10(c) as the red line) that would represent ca. 13% of the initial adsorbed deuterated material. However, a visual appraisal of Fig. 10(c) suggests that the post-oxidation reflectivity curve differs mainly from the null-reflecting water background in that it is unusually noisy at very low Q. From these noisy data alone it is very difficult to quantify the amount of residual monolayer, so further experimental confirmation would be useful. As well as repeating the low temperature residue quantification experiment for d-LOA/h-MO, future work on this system could perform analogous room temperature experiments. As for the OA ozonolysis in the presence of LOA discussed earlier, the kinetic data collected over a limited Q range (see Fig. 5–8) suggests that similar residual monolayers are present, but a full Q range characterisation would be very useful to confirm this. As further experiments were not feasible within the limitations of time-constrained beamline facility access, this was one of the key motivations for the more detailed modelling analysis presented in Section 3.4.
The mechanism for LOA oxidation is far more complex than that for OA oxidation; the former involves multiple competing pathways of ozonolysis and ozone-initiated autoxidation that result in different products. One mechanism or the other has previously been observed to dominate depending on ozone concentration and relative humidity.37 The linear increase in pseudo-first-order reaction rate with increasing ozone concentration suggests that one dominant mechanism is being observed across the range of ozone concentrations used in this work (ca. 100 to 1000 ppb). Chu et al. noted that higher ozone concentrations (above 250 ppb) inhibited the build-up of autoxidation products, as did higher relative humidity.37 All our work is carried out with ozone dissolved in a stream of dry oxygen, so the relative humidity will be very low. Based on the range of ozone concentrations we employed and the low relative humidity in our system, based on this literature we would have expected to observe mainly autoxidation except for the highest [O3] conditions. The autoxidation mechanism is complex and is not presented here as we were not conducting any product analysis to directly distinguish reaction mechanism, but rather focussing on reaction kinetics and residue formation in two-component monolayers at different temperatures (a detailed explanation of LOA autoxidation can be found in Chu et al.).37
Given the abundance of saturated surfactants in atmospheric aerosols, we also carried out exploratory experiments to test the effect of co-deposited stearic acid (SA) – the saturated analogue to OA – on the LOA monolayer oxidation. Fig. 11(a) displays the oxidation of a d-LOA/h-SA monolayers at two ozone concentrations at the two temperatures of interest. While this data set was too limited to reliably calculate a second-order rate coefficient, a visual appraisal of the data suggests that temperature once again does not exert a significant effect on the reaction rate. As far as the effect of SA on the reaction rate goes, however, this preliminary data suggest that there may be an impact. Fig. 11(b) illustrates the room temperature data overlaid with the analogous data for d-LOA/h-OA oxidation. It does appear that the reaction is faster in the presence of SA. Further work would be useful to establish whether this effect is reproducible, as there may be a degree of inter-run variability.
A full Q range residue quantification experiment was performed for the d-LOA/h-SA system at 3 ± 1 °C, and the results are shown in Fig. 11(c). The post-oxidation reflectivity appears no different from the air-NRW background barring an intriguing oscillation at low Q that would merit further investigation as would the residue behaviour at 21 ± 1 °C.
![]() | ||
Fig. 12 Multilayer-Py45 modelling fits to selected [O3] for pure d-LOA oxidation (fits to the full set of [O3] conditions is presented in the ESI;† the model runs presented here consider ozone build-up in the reaction chamber, a single-stage oxidation process and residue formation). The figures display the normalised decay of d-LOA as a function of time with experimental data (red and blue symbols), the results of the global optimisation using Markov chain Monte Carlo (MCMC) sampling (pink lines close to the best fit) and the optimised model fit (solid black line). (a) 983 ppb & 21 ± 1 °C; (b) 983 ppb & 3 ± 1 °C; (c) 492 ppb & 21 ± 1 °C; and (d) 492 ppb & 3 ± 1 °C. |
![]() | ||
Fig. 13 Multilayer-Py45 modelling fits to selected [O3] for d-LOA/h-OA oxidation (fits to the full set of [O3] conditions is presented in the ESI;† the model runs presented here consider ozone build-up in the reaction chamber, a single-stage oxidation process and residue formation). The figures display the normalised decay of d-LOA as a function of time with experimental data (red and blue symbols), the results of the global optimisation using Markov chain Monte Carlo (MCMC) sampling (pink lines close to the best fit) and the optimised model fit (solid black line). (a) 983 ppb & 21 ± 1 °C; (b) 983 ppb & 3 ± 1 °C; (c) 492 ppb & 21 ± 1 °C; and (d) 492 ppb & 3 ± 1 °C. |
![]() | ||
Fig. 14 Multilayer-Py45 modelling fits to selected [O3] for d-LOA/h-MO oxidation (fits to the full set of [O3] conditions is presented in the ESI;† the model runs presented here consider ozone build-up in the reaction chamber, a single-stage oxidation process and residue formation). The figures display the normalised decay of d-LOA as a function of time with experimental data (red and blue symbols), the results of the global optimisation using Markov chain Monte Carlo (MCMC) sampling (pink lines close to the best fit) and the optimised model fit (solid black line). (a) 983 ppb & 21 ± 1 °C; (b) 983 ppb & 3 ± 1 °C; (c) 492 ppb & 21 ± 1 °C; and (d) 492 ppb & 3 ± 1 °C. |
The kinetic multi-layer modelling demonstrated that model optimisation using Markov chain Monte Carlo (MCMC) sampling yields good fits for the rate coefficients obtained with standard kinetic fitting (see fits in Fig. 12–14 which are using the rate coefficients stated in Table 1).
Furthermore, Multilayer-Py allowed optimised residue fitting. The results (see Table 2) suggest that the residues for pure d-LOA (Fig. 12) are consistently lower (11–14%) than for its mixtures with h-OA (13–27%; Fig. 13) followed by the mixtures with h-MO (19–27%; Fig. 14). The modelling showed no significant difference between residues at the two different temperatures: only for the d-LOA/h-MO mixtures (see Table 2) did we find a small difference with the residue averaging at ca. 23% for 3 ± 1 °C compared to ca. 20% for 21 ± 1 °C, however no clear difference was observed for the other monolayer compositions.
Monolayer | Relative residue/% | |||||||
---|---|---|---|---|---|---|---|---|
21 °C & 983 ppb | 21 °C & 492 ppb | 21 °C & 246 ppb | 21 °C & 98 ppb | 3 °C & 983 ppb | 3 °C & 492 ppb | 3 °C & 246 ppb | 3 °C & 98 ppb | |
d-LOA | 14.48 | 12.86 | 14.31 | 12.07 | 10.87 | 13.59 | 12.39 | 13.17 |
d-LOA/h-OA | 19.38 | 13.90 | 16.22 | 20.21 | 16.40 | 14.31 | 12.80 | 26.79 |
d-LOA/h-MO | 21.32 | 19.14 | 20.68 | 19.94 | 22.49 | 19.76 | 22.94 | 27.39 |
Overall, the modelling at room temperature results in slightly better fits across the different monolayer compositions at all [O3] except for the lowest concentration of 98 ppb (see Fig. S25 in the ESI†). For pure d-LOA at low [O3] (98 ppb) the decay shape consistently deviates from that compatible with a single decay process (see Fig. S25(a) & (b)†); this deviation is weaker, but still visible, for the mixtures with h-MO and less apparent for the h-OA mixtures (contrast Fig. S25(a)–(d)†); the data were too limited to quantify this effect, but this observation is consistent with previous work suggesting that the additional process of autoxidation becomes important at low O3.37
Kinetic multi-layer modelling using our Multilayer-Py package allowed us to optimise the residue fitting as residue formation showed a temperature dependence in our previous work for oleic acid ozonolysis (Woden et al., 2021).33 While the modelling results presented here suggest that the residues for pure d-LOA are consistently lower (11–14%) than for its mixtures with h-OA (13–27%) and h-MO (19–27%), it showed no clear difference between residues at the two different temperatures: only for the d-LOA/h-MO mixtures did we find a small difference with the residue averaging at ca. 23% for 3 ± 1 °C compared to ca. 20% for 21 ± 1 °C, but no significant difference was observed for the other monolayer compositions across the different O3 levels.
Since the low-temperature conditions used here are atmospherically realistic, it is key to understand if a product film persists and thus needs to be considered when assessing the impact of unsaturated fatty acid partitioned to the air–water interface. The presence of stable (non-oxidisable) reaction products could also lead to a build-up of inert monolayers during the aerosol life cycle with potential implications for cloud formation. Our previous work on oleic acid ozonolysis33 showed that a residual surface film (likely formed of ozonolysis products nonanoic acid and a mixture of azelaic and 9-oxononanoic acids) was retained at the interface after ozonolysis at low temperatures, but not at room temperature. For the binary mixtures studied here, we did not find such a clear temperature dependence of residue formation. As surfactants will be present in the atmosphere in complex, multi-component mixtures, it is important to understand the reasons for these different behaviours even for closely related systems. While the temperature change did not impact on the kinetics, residue formation may be affected, thus altering the persistence of the organic character at the surface of aqueous droplets.
Furthermore, it would be very useful to explore further the kinetic behaviour at lower ozone concentrations (as long as reactions can still be observed within the limited timeframe of neutron beamtime experiments) given the unusual decay shapes we found especially for pure d-LOA at 98 ppb. This would also extend the range of data available for second-order fitting increasing the robustness of the fitted parameters.
Finally, it would be worth establishing the residue remaining after reaction for a wider range of conditions also including full Q range runs at room temperature.
Kinetic multi-layer modelling was used to: (i) confirm the rate coefficients obtained with standard kinetic fitting; (ii) illustrate that for pure d-LOA at low [O3] (98 ppb) the decay shape consistently deviates from that compatible with a single decay process; this deviation was found to be weaker, but still visible, for the mixtures with h-MO and less apparent for the h-OA mixtures (this observation is consistent with previous work suggesting that the additional process of autoxidation becomes important at low O3, but the reason for the differences between pure d-LOA and the two-component mixtures merits further investigation); and (iii) optimise the residue fitting suggesting that the residues for pure d-LOA are consistently lower (11–14%) than for its mixtures with h-OA (13–27%) followed by the mixtures with h-MO (19–27%). The latter trend in the model data is consistent with full Q range NR experiments carried out at the highest [O3] conditions, although these experiments suggest a lower residue proportion ranging from ca. 7% to 13%.
In summary, neither the change in temperature nor the introduction of co-deposited film components alongside d-LOA consistently affected the LOA oxidation rates, but the deviation from a single process decay behaviour (indicative of autoxidation) at 98 ppb is clearest for pure d-LOA, weaker for the h-MO mixtures and the weakest for h-OA mixtures.
As surfactants will be present in the atmosphere in complex, multi-component mixtures, it is important to understand the reasons for these different behaviours even in two-component mixtures of closely related species. The rates we found were fast compared to those reported earlier for a different LOA morphology under much higher ozone concentrations. Our work demonstrates clearly that it is essential to employ atmospherically realistic ozone levels as well as multi-component mixtures especially to understand LOA behaviour at low O3 in the atmosphere. While the temperature change did not play a crucial role for the kinetics, residue formation may be affected, potentially impacting on the persistence of the organic character at the surface of aqueous droplets with a wide range of atmospheric implications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4fd00167b |
This journal is © The Royal Society of Chemistry 2025 |