Multiphase reactivity of polycyclic aromatic hydrocarbons is driven by phase separation and diffusion limitations.

Benzo[a]pyrene (BaP), a key polycyclic aromatic hydrocarbon (PAH) often associated with soot particles coated by organic compounds, is a known carcinogen and mutagen. When mixed with organics, the kinetics and mechanisms of chemical transformations of BaP by ozone in indoor and outdoor environments are still not fully elucidated. Using direct analysis in real-time mass spectrometry (DART-MS), kinetics studies of the ozonolysis of BaP in thin films exhibited fast initial loss of BaP followed by a slower decay at long exposure times. Kinetic multilayer modeling demonstrates that the slow decay of BaP over long times can be simulated if there is slow diffusion of BaP from the film interior to the surface, resolving long-standing unresolved observations of incomplete PAH decay upon prolonged ozone exposure. Phase separation drives the slow diffusion time scales in multicomponent systems. Specifically, thermodynamic modeling predicts that BaP phase separates from secondary organic aerosol material so that the BaP-rich layer at the surface shields the inner BaP from ozone. Also, BaP is miscible with organic oils such as squalane, linoleic acid, and cooking oil, but its oxidation products are virtually immiscible, resulting in the formation of a viscous surface crust that hinders diffusion of BaP from the film interior to the surface. These findings imply that phase separation and slow diffusion significantly prolong the chemical lifetime of PAHs, affecting long-range transport of PAHs in the atmosphere and their fates in indoor environments.

Polycyclic aromatic hydrocarbons (PAHs), including benzo[a]pyrene (BaP), are among the most prominent toxic air pollutants, posing a threat to human health because their metabolites and oxidation products are carcinogenic and mutagenic (1). PAHs are emitted into the atmosphere from incomplete combustion and biomass burning and by smoking and cooking in indoor environments. Due to its low vapor pressure, BaP resides mostly in the condensed phase, and heterogeneous oxidation of BaP by oxidants such as OH and O3 is a major atmospheric loss pathway (2). Laboratory measurements show rapid degradation of BaP against ozone when adsorbed to a variety of substrates such as water, ammonium sulfate, soot, and organic compounds (3–6).

Upon chemical aging in the atmosphere, PAH-containing particles are likely coated by semi- or low-volatile organic compounds, which are formed by multigenerational gas-phase oxidation of volatile organic compounds. Laboratory experiments have shown that organic coatings can retard multiphase reactions of ozone with PAHs due to kinetic limitations of bulk diffusion (3, 7, 8). The extent of coating effects depends on the phase state of organic coatings, which can be liquid, amorphous semisolid, or glassy solid, depending on chemical composition, relative humidity (RH), and temperature (9). This shielding effect of PAHs from oxidation was recently implemented into regional and global air-quality models showing how regional and global distributions and transport of BaP can be affected by the temperature and humidity dependence of diffusivity and reactivity of BaP-containing particles (10, 11). These modeling results are consistent with observations of BaP at remote terrestrial and marine sites, even in the polar regions (2, 12).

Atmospheric aerosol particles are often mixtures of organics, inorganics, and water, which are subject to complex and nonideal behavior including liquid–liquid and/or liquid–solid phase separation, as predicted by thermodynamic models (13, 14). Experimental studies have demonstrated the existence of multiple phases in model mixtures or laboratory-generated secondary organic aerosol (SOA) particles as well as in field-collected organic aerosol particles (15–17). Liquid–liquid phase separation was also observed in SOA particles free of inorganic salts (18–20). The interplay of phase state and nonideality can affect aerosol mass concentration and the characteristic time scale of gas-particle mass transfer (21). Changes in viscosity upon chemical aging of organic particles were observed (22), but it is not yet clear how the effects and interplay of nonideality and phase state evolve upon multiphase chemical interactions to affect fates and atmospheric long-range transport of organic compounds.

Although extensive laboratory measurements and modeling of PAHs have been conducted under atmospherically relevant conditions, chemical transformation of PAHs in indoor environments, where we spend ∼90% of the time, is poorly understood (23). PAHs can be transported from outside air and can also be emitted by indoor activities such as smoking, cooking, and burning of solid fuels (24). Recent findings show that BaP diol-epoxide products, which are highly carcinogenic, can be formed in indoor environments (25), and chronic exposure to these compounds can be detrimental to human health. Thus, both the chemical transformation of PAHs with gaseous ozone under indoor-relevant conditions and the PAH reactivity in mixtures of chemicals emitted from indoor sources need to be elucidated. A long-standing, unexplained issue in the field of environmental chemistry and health is the observation that although ozonolysis experiments show a rapid initial loss of condensed-phase PAHs, prolonged exposure to ozone does not necessarily lead to complete PAH decay (26–29).

In this study, the decay kinetics of BaP upon exposure to ozone were measured using direct analysis in real-time mass spectrometry (DART-MS). Three different sets of experimental and modeling scenarios were explored: (i) the ozonolysis kinetics of pure BaP films (containing a small amount of an internal standard); (ii) the effects of the reaction substrate by mixing BaP with α-pinene SOA material, a mixture of different oxidized organic compounds; and (iii) the effects of mixing BaP with liquid substrates such as squalane, linoleic acid, and cooking oil. The experimental conditions were simulated using the kinetic multilayer model of aerosol surface and bulk chemistry (KM-SUB), which resolves mass transport and chemical reactions at the surface and in the condensed phase explicitly (30), combined with the Aerosol Inorganic-Organic Mixtures Functional Groups Activity Coefficients (AIOMFAC) thermodynamic model (31–33) at the core of a thermodynamic equilibrium framework that predicts chemical composition of different coexisting phases. Together, the 3 sets of simulations and experimental results pinpoint the critical factors that control PAH ozonolysis loss rates in a variety of organic media.

Results and Discussion

Fig. 1 shows the ozonolysis decays of pure BaP in thin films that only contain traces of bis(2-ethylhexyl) sebacate (BES) as an internal standard. There is a 1-h exposure to 15 to 1,000 ppb ozone (O3) at <5% RH and 296 K. If the BaP were uniformly deposited on the capillary, then the resulting film would be roughly 1 monolayer thick. Much more likely, the film has inhomogeneous thickness, leading to most of the BaP being present in multilayer amounts. The BaP concentration decayed slower at 15 ppb O3 than at 100 ppb O3, while BaP decay rates became saturated and very similar at 500 and 1,000 ppb O3. Interestingly, BaP was not entirely reacted away, leaving an unreacted fraction of ∼20% at high O3 mixing ratios (500 and 1000 ppb) after 1 h, which has also been observed previously, but not explained for many decades (26–29). KM-SUB was applied to simulate these data by considering surface adsorption of O3, decomposition of O3 into reactive oxygen intermediates (ROIs; e.g., sorbed O atoms) and subsequent reactions with BaP (34), bulk diffusion of O3 and BaP, and bulk reactions between O3 and BaP. The time and concentration dependence of BaP decay can be reproduced by KM-SUB with the kinetic parameter values listed in , as consistent with our previous work (34). The film thickness was set in the model to be 0.61 nm, representing roughly a double layer of molecules. Note that qualitatively similar modeling results were obtained using a film thickness 10 times larger (see for details).

Fig. 1.

Decay of BaP concentration in BaP/BES films exposed to different mixing ratios of ozone (15 to 1,000 ppb) at 296 K and RH <5%. Filled circles show the experimental data, with error bars representing the SD of 10 measurements. The solid lines are KM-SUB simulation results with a single kinetic parameter set shown in .

The modeling results indicate that the decay of BaP is initially controlled by ROI formation and surface reactions with BaP, leading to fast initial BaP decay. Once surface BaP is depleted, the key aspect of the model prediction is that the oxidation of BaP becomes limited by bulk diffusion of BaP from the film interior to the surface at longer reaction times. This explains the subsequent slow BaP decay after ∼10 min, as well as the saturation of BaP decay rates at higher O3 mixing ratios. Bulk diffusivity of BaP was predicted to be in the range of ∼10−17 to 10−19 cm2⋅s−1, which is consistent with an amorphous semisolid state of the film. Note that the modeled kinetics are insensitive to the bulk diffusivity of O3, as bulk reactions play a minor role in overall BaP degradation (3, 7). As a consequence, BaP in the deep bulk is protected from ozonolysis due to this shielding effect arising from slow BaP diffusion to the surface, which explains the unreacted BaP fraction at long reaction times in these and previous experiments (26–29).

The results from the second reaction scenario are shown in Fig. 2, which illustrates the time-dependent oxidative decay of BaP mixed with α-pinene SOA exposed to 500 ppb O3 at dry conditions, 50% RH, and 85% RH. The BaP concentration decays at very similar rates at different RHs for up to ∼10 min and then decays with faster rates and higher reacted fractions at higher RH. In theory, this behavior could be due to a phase change of α-pinene SOA, which is expected to exist as an amorphous solid at a temperature of ∼295 K and at dry conditions but becomes substantially less viscous at higher RH due to water uptake (35). According to previous viscosity measurements (35) and kinetics experiments (7), the bulk diffusivity of PAH in α-pinene SOA at 85% RH is expected to be as high as ∼10−11 cm2⋅s−1. However, when PAH is assumed to be well mixed and with such a high bulk diffusivity in the α-pinene SOA mixture, the KM-SUB simulation failed to reproduce the experimental data, as shown by the red dashed line in Fig. 2. Under such assumptions, KM-SUB predicts complete depletion of BaP within less than 1 h, indicating a lack of substantial kinetic limitations, even at 85% RH. There is clearly an aspect to this modeling scenario that is not accurately representing the experimentally observed kinetics.

Fig. 2.

(A) Decay of BaP embedded in the films of α-pinene SOA upon exposure to 500 ppb O3 at dry conditions (blue), 50% RH (green), and 85% RH (red) at 296 K. The circle markers are the experimental data, and error bars represent the SD of 10 measurements. KM-SUB modeling results are presented with consideration of liquid–liquid phase separation (solid lines) and without (dashed line for 85% RH). (B) Thermodynamic modeling results by AIOMFAC predicting phase compositions of the α-pinene SOA + BaP/BES system as a function of water activity (i.e., RH) for an initial amount of 10 ng of α-pinene SOA, 1 ng of BaP, and 0.5 ng of BES in the particle phase; i.e., before reaction with ozone. In particular, each frame corresponds to a different phase, within which the relative proportions of different components in the phase are indicated. Specifically, predicted phase compositions are shown as stacked mass fractions of the individual components for the phase α enriched by SOA compounds (Lower) and the phase β dominated by BaP and BES (Upper). Note that α-pinene SOA is treated as comprising 14 representative oxidation products (see for explanation of labels).

To probe this reaction system more deeply, the likelihood of phase separation was assessed using the AIOMFAC-based liquid–liquid equilibrium (LLE) model (31–33), where α-pinene SOA was simulated by a mixture of 14 representative oxidation products () (13, 21). Contrary to the assumption of full mixing which led to the dashed line in Fig. 2, the LLE model predicts that the mixtures separate into BaP-rich (Fig. 2 , Upper) and SOA-rich (Fig. 2 , Lower) phases of distinct compositions, dependent on the water activity. Under dry conditions, ∼65% of the total particle-bound BaP mass is predicted to partition to the phase rich in BES and BaP, with the remaining ∼35% BaP present in the SOA-rich phase. BaP’s phase preference becomes more distinct with increasing RH so that at 50% RH, about 90% of its mass resides in the BaP-rich phase, and at 85% RH, the partitioning to the BaP-rich phase approaches 100%; that is, virtually complete phase separation. The BaP-rich phase is most likely in the form of a film on top of the SOA-rich phase (i.e., full engulfing as opposed to partial engulfing), considering that BaP should have lower surface tension than polar SOA materials.

By implementing the predicted phase compositions arising from phase separation in KM-SUB, the model reproduces experimental data very well, as shown by the solid lines in Fig. 2. The kinetic model simulations reveal that the BaP decay is initially controlled by surface reactions leading to fast BaP decay, followed by diffusion of BaP molecules from the bulk to its surface. The diffusion coefficients of BaP within the surface crust (i.e., a semisolid surface layer of high viscosity) were found to have the most sensitivity in the modeling. Reactions at the surface are found to be faster than in the bulk (7), so BaP in the film interior needs to diffuse through the bulk to the surface to react. The BaP-rich layer at the surface hinders the BaP in the film interior to diffuse to the surface, resulting in significant fractions of BaP remaining unreacted after 1 h of ozone exposure.

The final set of kinetics experiments measured the decay behavior of BaP embedded in other organic liquids, including squalane, linoleic acid, and cooking oil, as shown in Fig. 3. Contrary to the SOA case, BaP is miscible with these organic compounds, as evident by visual inspection and as predicted by the thermodynamic phase equilibrium model. So, BaP should initially be homogenously mixed within the film. The KM-SUB modeling results for BaP decay in squalane upon exposure to 500 ppb O3 are shown as the dashed red line in Fig. 3. The model captures the initial BaP decay up to ∼10 min but overpredicts the long-term BaP decay in the absence of kinetic limitations.

Fig. 3.

(A) Decay of BaP embedded in the films of squalane (red), linoleic acid (yellow), or cooking oil (green) upon exposure to 500 ppb O3 under dry conditions. The circle markers are the experimental data, and error bars represent the SD of 10 measurements. The solid lines represent KM-SUB simulations with consideration of surface crust formation by BaP oxidation products and composition-dependent bulk diffusivity using a Vignes-type equation, while the dashed line for the squalane case does not consider this phase separation effect, failing to reproduce experiments. (B) Predicted phase compositions in mass fractions for mixtures of BaP, BaP oxidation products (BaP-6,2-dione, BaP-derived carboxylic acid, and 6,12-dihydroxy BaP), and squalane after half of the BaP is degraded (see Fig. 2 caption for explanation of the plots). One phase consists largely of squalene, BES, and BaP (β, Lower), whereas the other consists largely of the BaP oxidation products (α, Upper). LLE calculations with AIOMFAC predict for all RH (or water activities) that BaP oxidation products are virtually immiscible in squalane.

Although BaP may be miscible initially, its oxidation products may phase separate from the oily liquid reaction substrates. We considered 3 major BaP ozonolysis products, including 6,12-dihydroxy-BaP, BaP-6,12-dione, and a BaP-derived dicarboxylic acid (36), and cooking oil was treated as triolein (). The AIOMFAC thermodynamic equilibrium modeling revealed that these oxidation products are immiscible with all 3 organic liquids and the system should therefore undergo phase separation upon formation of such oxidation products (Fig. 3 and ). This suggests that the oxidation process, which leads to the formation of oxidation products that are more oxygenated than the BaP reactant, promotes phase separation. A viscous surface crust forms, which can hinder diffusion of BaP from the film interior to the surface. We note that a crusting effect has been observed in the uptake of organic nitrates (37) or ammonia to SOA, where carboxylates may also form a viscous crust at the surface of the reaction substrate (38).

To consider this effect in KM-SUB, bulk diffusion coefficients of BaP were treated as composition dependent using the Vignes-type equation (39) (). Linoleic acid and cooking oil compounds are unsaturated (i.e., with C=C double bonds), and ozone reacts with these compounds (). Previous studies have shown that high-molar-mass compounds such as peroxidic oligomers can be formed, especially under low RH (40, 41), causing an increase of viscosity (22). Thus, these products may also contribute to form the surface crust, which was explicitly treated in KM-SUB. With consideration of these effects, KM-SUB was able to capture a slow BaP decay at longer reaction times to fully reproduce the experimental data, as shown by the solid lines in Fig. 3. Retardation of oleic acid ozonolysis from surface crust formation has been observed in mixed-component particles in previous studies (42–44). These results demonstrate how phase separation leading to crust formation at the surface can affect the chemical transformation of PAH in multicomponent mixtures.

Implications for Atmospheric and Indoor Chemistry.

PAHs are often associated with soot particles that are subject to coating by semivolatile primary organic emissions and by secondary organic materials formed by oxidation of anthropogenic and biogenic precursors. Our experimental and modeling study demonstrates that the multiphase reactivity of PAHs depends strongly on the interplay of phase state and nonideal mixing with such organic materials.

In all 3 sets of reaction systems, it was found that slow diffusion through viscous phases is rate limiting at long reaction times because the BaP must diffuse to the surface to react. In the case of pure BaP films, the semisolid substrate acts as the diffusion barrier, whereas in the mixtures with oils such as squalane, linoleic, or cooking oils, the reaction products phase separate into a phase with low diffusivity. When mixed into SOA, the PAHs do not mix homogeneously with α-pinene SOA but are instead separated into different phases, affecting the multiphase reactivity of PAHs against ozone. RH controls the degree of nonideality and phase separation and also determines the phase state of PAH–SOA mixtures. In phase-separated particles with the PAH-rich phase on the particle surface, slow diffusion of PAH from the particle bulk to the surface limits PAH degradation by ozone. Moreover, the SOA-rich phase, which embeds a fraction of the PAHs, adopts an amorphous (semi)solid state at low RH, posing additional kinetic limitations of bulk diffusion and retarding ozonolysis kinetics. Overall, these effects will prolong chemical lifetimes of PAHs, facilitating long-range transport of PAHs and affecting regional and global distributions in the atmosphere and their effects on air quality and public health (10, 11).

PAHs are also found in indoor environments, and a BaP concentration of 38.8 µg/g of dust has been reported (45). PAHs are transported from outdoors to indoors and are directly emitted by indoor activities such as cooking and smoking. With cooking, BaP may be deposited and embedded into an organic film on indoor surfaces. The fate of low-volatility PAHs will be controlled by surface oxidative processes, considering the high surface-to-volume ratios in indoor environments and that organic films can persist for an indefinite time unless they are wiped off (46).

To investigate PAH degradation kinetics with real indoor air, a BaP film was exposed to room air containing 10 ppb O3 for up to 5 h. The BaP decay rate by room air is remarkably similar to the one when BaP is exposed to 15 ppb O3 in a controlled flow tube setting, and it is modeled well (the solid red line in ) with the same kinetic parameters used to model Fig. 1. These results indicate that ozone will be the most important indoor oxidant, among others (e.g., OH, NO3), for determining the fate of these compounds in indoor environments.

To illustrate the chemical fate of BaP embedded in an indoor surface film consisting of cooking oil, we simulated the molar fraction of unreacted BaP after exposure to 0.1 to 100 ppb O3 for 1 h, 12 h, 1 d, and 1 wk. The parameters used in the model are those that successfully modeled the experimental data in Fig. 3 (). The mass concentration of BaP in cooking oil was assumed to be 1,200 ng⋅cm−3 and the film thickness was set to be 8 nm based on previous observations in indoor environments (47–49). Fig. 4 shows the results of such KM-SUB simulations. At an O3 concentration of 10 ppb, ∼80% and ∼10% of the BaP remain unreacted after 1 h and 12 h, respectively, whereas the BaP is expected to be fully reacted away after 1 wk. Typical indoor O3 concentrations range from 5 ppb to over 50 ppb (depending on air exchange rates and outdoor air pollution) (50), and the unreacted fraction of BaP molecules considering surface crust formation is larger than the case with an ideal mixing assumption, especially at higher O3 concentrations. Thus, the fate of PAH in indoor environments is strongly influenced by multiphase reactivity as impacted by surface crust formation and diffusion limitations. The true BaP lifetime may be even longer and is dependent upon the mass transfer rates that prevail under genuine indoor conditions.

Fig. 4.

Fraction of BaP remaining in an 8-nm-thick cooking oil film exposed to different gas-phase ozone mixing ratios after 1 h (black), 12 h (purple), 1 d (blue), and 1 wk (red) at 295 K and dry conditions, as simulated by KM-SUB considering phase separation and crust formation by BaP oxidation products (solid lines) or assuming ideal mixing (dashed lines).

Materials and Methods

Kinetics Experiments.

A solution containing BaP and BES was deposited on the outside sealed end of glass melting point capillary tubes, and a thin film of BaP/BES was formed after the solvent evaporated. BaP was also mixed with other organics, namely squalane, cooking oil, linoleic acid, or SOA, as formed in a 1-m3 Teflon chamber from the dark reaction of α-pinene with ozone. Heterogeneous oxidation of BaP with ozone was carried out at 295 ± 3 K in a flow tube. The indoor air oxidation of BaP was performed by placing the capillaries inside a fume hood in the Chemistry Department, University of Toronto, where the temperature, RH, and ozone mixing ratio was monitored. After oxidation, the capillaries were stored in a desiccator where a flow of nitrogen was added to avoid contamination/oxidation from the room air before analysis using DART-MS, with the approach described in full previously (51). Further details are described in .

Kinetic and Thermodynamic Modeling.

KM-SUB (30) was used for simulating kinetics experiments of BaP ozonolysis. The AIOMFAC thermodynamic model (31–33) was used to estimate thermodynamic mixing and the potential for phase separation in the PAH/organics/water system. Predicted phase compositions and associated-component activity coefficients were considered in KM-SUB. KM-SUB consists of a number of compartments and layers in which ozone, BaP, and oxidized products can undergo mass transport and chemical reactions in the particle phase. Decomposition of surface-adsorbed ozone into ROIs (e.g., O atoms) followed by surface reactions with BaP was simulated explicitly. The required kinetic parameters include the surface accommodation coefficient and the desorption lifetime of ozone, bulk diffusion coefficients of ozone and BaP, and the second-order surface and bulk reaction rate coefficients between ozone and BaP as well as those between ozone and unsaturated organic oils (). Nanoscale effects are not considered in this modeling. While beyond the scope of this study, they would need to be included to more accurately describe the composition and behavior of very thin surface layers. From theoretical considerations and experimental work, it is known that the influence of liquid–liquid and air–liquid interfaces becomes increasingly important at the nanometer scale (52, 53). The proximity of molecules/phases to such interfaces can lead to enhancement or reduction of the thermodynamic drivers for phase separation and/or bulk–surface partitioning, depending on the system (54, 55). Hence, the AIOMFAC-based calculations, while performed for bulk-system conditions, also apply approximately to nanoscale systems for the case where the interfacial energy effects nearly balance. The degree to which this case applies to the systems studied here is unknown yet is of interest for future experimental and modeling studies on aerosol size and interface effects. Further method details are described in .