skip to main content

Title: Non-equilibrium interplay between gas–particle partitioning and multiphase chemical reactions of semi-volatile compounds: mechanistic insights and practical implications for atmospheric modeling of polycyclic aromatic hydrocarbons
Abstract. Polycyclic aromatic hydrocarbons (PAHs) are carcinogenic air pollutants. The dispersion of PAHs in the atmosphere is influenced by gas–particle partitioning and chemical loss. These processes are closely interlinked and may occur at vastly differing timescales, which complicates their mathematical description in chemical transport models. Here, we use a kinetic model that explicitly resolves mass transport and chemical reactions in the gas and particle phases to describe and explore the dynamic and non-equilibrium interplay of gas–particle partitioning and chemical losses of PAHs on soot particles. We define the equilibration timescale τeq of gas–particle partitioning as the e-folding time for relaxation of the system to the partitioning equilibrium. We find this metric to span from seconds to hours depending on temperature, particle surface area, and the type of PAH. The equilibration time can be approximated using a time-independent equation, τeq≈1kdes+kads, which depends on the desorption rate coefficient kdes and adsorption rate coefficient kads, both of which can be calculated from experimentally accessible parameters. The model reveals two regimes in which different physical processes control the equilibration timescale: a desorption-controlled and an adsorption-controlled regime. In a case study with the PAH pyrene, we illustrate how chemical loss can perturb the equilibrium particulate more » fraction at typical atmospheric concentrations of O3 and OH. For the surface reaction with O3, the perturbation is significant and increases with the gas-phase concentration of O3. Conversely, perturbations are smaller for reaction with the OH radical, which reacts with pyrene on both the surface of particles and in the gas phase. Global and regional chemical transport models typically approximate gas–particle partitioning with instantaneous-equilibration approaches. We highlight scenarios in which these approximations deviate from the explicitly coupled treatment of gas–particle partitioning and chemistry presented in this study. We find that the discrepancy between solutions depends on the operator-splitting time step and the choice of time step can help to minimize the discrepancy. The findings and techniques presented in this work not only are relevant for PAHs but can also be applied to other semi-volatile substances that undergo chemical reactions and mass transport between the gas and particle phase. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Atmospheric Chemistry and Physics
Page Range or eLocation-ID:
6175 to 6198
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Secondary organic aerosols (SOA) account for a substantial fraction of airparticulate matter, and SOA formation is often modeled assuming rapidestablishment of gas–particle equilibrium. Here, we estimate thecharacteristic timescale for SOA to achieve gas–particle equilibrium undera wide range of temperatures and relative humidities using astate-of-the-art kinetic flux model. Equilibration timescales werecalculated by varying particle phase state, size, mass loadings, andvolatility of organic compounds in open and closed systems. Modelsimulations suggest that the equilibration timescale for semi-volatilecompounds is on the order of seconds or minutes for most conditions in theplanetary boundary layer, but it can be longer than 1 h if particlesadopt glassy or amorphous solid states with high glass transitiontemperatures at low relative humidity. In the free troposphere with lowertemperatures, it can be longer than hours or days, even at moderate orrelatively high relative humidities due to kinetic limitations of bulkdiffusion in highly viscous particles. The timescale of partitioning oflow-volatile compounds into highly viscous particles is shorter compared tosemi-volatile compounds in the closed system, as it is largely determined bycondensation sink due to very slow re-evaporation with relatively quickestablishment of local equilibrium between the gas phase and thenear-surface bulk. The dependence of equilibration timescales on bothvolatility and bulk diffusivity provides critical insightsmore »intothermodynamic or kinetic treatments of SOA partitioning for accuratepredictions of gas- and particle-phase concentrations of semi-volatilecompounds in regional and global chemical transport models.« less
  2. Abstract. Mass accommodation is an essential process for gas–particle partitioning oforganic compounds in secondary organic aerosols (SOA). The massaccommodation coefficient is commonly described as the probability of a gasmolecule colliding with the surface to enter the particle phase. It is oftenapplied, however, without specifying if and how deep a molecule has topenetrate beneath the surface to be regarded as being incorporated into thecondensed phase (adsorption vs. absorption). While this aspect is usuallynot critical for liquid particles with rapid surface–bulk exchange, it canbe important for viscous semi-solid or glassy solid particles to distinguishand resolve the kinetics of accommodation at the surface, transfer acrossthe gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective massaccommodation coefficient αeff that depends on penetrationdepth and is a function of surface accommodation coefficient, volatility,bulk diffusivity, and particle-phase reaction rate coefficient. Applicationof αeff in the traditional Fuchs–Sutugin approximation ofmass-transport kinetics at the gas–particle interface yields SOApartitioning results that are consistent with a detailed kinetic multilayermodel (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Modelfor Simulating Aerosol Interactions and Chemistry, MOSAIC;Zaveri et al., 2014) but deviate substantially frommore »earlier modelingapproaches not considering the influence of penetration depth and relatedparameters. For highly viscous or semi-solid particles, we show that the effective massaccommodation coefficient remains similar to the surface accommodationcoefficient in the case of low-volatility compounds, whereas it can decrease byseveral orders of magnitude in the case of semi-volatile compounds. Such effectscan explain apparent inconsistencies between earlier studies deriving massaccommodation coefficients from experimental data or from molecular dynamicssimulations. Our findings challenge the approach of traditional SOA models using theFuchs–Sutugin approximation of mass transfer kinetics with a fixed massaccommodation coefficient, regardless of particle phase state and penetrationdepth. The effective mass accommodation coefficient introduced in this studyprovides an efficient new way of accounting for the influence of volatility,diffusivity, and particle-phase reactions on SOA partitioning in processmodels as well as in regional and global air quality models. While kineticlimitations may not be critical for partitioning into liquid SOA particlesin the planetary boundary layer (PBL), the effects are likely important foramorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT)as well as in the PBL at low relative humidity and low temperature.« less
  3. The daytime oxidation of biogenic hydrocarbons is attributed to both OH radicals and O3, while nighttime chemistry is dominated by the reaction with O3 and NO3 radicals. Here, the diurnal pattern of Secondary Organic Aerosol (SOA) originating from biogenic hydrocarbons was intensively evaluated under varying environmental conditions (temperature, humidity, sunlight intensity, NOx levels, and seed conditions) by using the UNIfied Partitioning Aerosol phase Reaction (UNIPAR) model, which comprises multiphase gas-particle partitioning and in-particle chemistry. The oxidized products of three different hydrocarbons (isoprene, α-pinene, and β-caryophyllene) were predicted by using near explicit gas mechanisms for four different oxidation paths (OH, O3, NO3, and O(3P)) during day and night. The gas mechanisms implemented the Master Chemical Mechanism (MCM v3.3.1), the reactions that formed low volatility products via peroxy radical (RO2) autoxidation, and self- and cross-reactions of nitrate-origin RO2. In the model, oxygenated products were then classified into volatility-reactivity base lumping species, which were dynamically constructed under varying NOx levels and aging scales. To increase feasibility, the UNIPAR model that equipped mathematical equations for stoichiometric coefficients and physicochemical parameters of lumping species was integrated with the SAPRC gas mechanism. The predictability of the UNIPAR model was demonstrated by simulating chamber-generated SOA data undermore »varying environments day and night. Overall, the SOA simulation decoupled to each oxidation path indicated that the nighttime isoprene SOA formation was dominated by the NO3-driven oxidation, regardless of NOx levels. However, the oxidation path to produce the nighttime α-pinene SOA gradually transited from the NO3-initiated reaction to ozonolysis as NOx levels decreased. For daytime SOA formation, both isoprene and α-pinene were dominated by the OH-radical initiated oxidation. The contribution of the O(3P) path to all biogenic SOA formation was negligible in daytime. Sunlight during daytime promotes the decomposition of oxidized products via photolysis and thus, reduces SOA yields. Nighttime α-pinene SOA yields were significantly higher than daytime SOA yields, although the nighttime α-pinene SOA yields gradually decreased with decreasing NOx levels. For isoprene, nighttime chemistry yielded higher SOA mass than daytime at the higher NOx level (isoprene/NOx > 5 ppbC/ppb). The daytime isoprene oxidation at the low NOx level formed epoxy-diols that significantly contributed SOA formation via heterogeneous chemistry. For isoprene and α-pinene, daytime SOA yields gradually increased with decreasing NOx levels. The daytime SOA produced more highly oxidized multifunctional products and thus, it was generally more sensitive to the aqueous reactions than the nighttime SOA. β-Caryophyllene, which rapidly oxidized and produced SOA with high yields, showed a relatively small variation in SOA yields from changes in environmental conditions (i.e., NOx levels, seed conditions, and diurnal pattern), and its SOA formation was mainly attributed to ozonolysis day and night. To mimic the nighttime α-pinene SOA formation under the polluted urban atmosphere, α-pinene SOA formation was simulated in the presence of gasoline fuel. The simulation suggested the growth of α-pinene SOA in the presence of gasoline fuel gas by the enhancement of the ozonolysis path under the excess amount of ozone, which is typical in urban air. We concluded that the oxidation of the biogenic hydrocarbon with O3 or NO3 radicals is a source to produce a sizable amount of nocturnal SOA, despite of the low emission at night.« less
  4. 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 filmmore »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.« less
  5. Atmospheric aerosol particles with a high viscosity may become inhomogeneously mixed during chemical processing. Models have predicted gradients in condensed phase reactant concentration throughout particles as the result of diffusion and chemical reaction limitations, termed chemical gradients. However, these have never been directly observed for atmospherically relevant particle diameters. We investigated the reaction between ozone and aerosol particles composed of xanthan gum and FeCl 2 and observed the in situ chemical reaction that oxidized Fe 2+ to Fe 3+ using X-ray spectromicroscopy. Iron oxidation state of particles as small as 0.2 μm in diameter were imaged over time with a spatial resolution of tens of nanometers. We found that the loss off Fe 2+ accelerated with increasing ozone concentration and relative humidity, RH. Concentric 2-D column integrated profiles of the Fe 2+ fraction, α , out of the total iron were derived and demonstrated that particle surfaces became oxidized while particle cores remained unreacted at RH = 0–20%. At higher RH, chemical gradients evolved over time, extended deeper from the particle surface, and Fe 2+ became more homogeneously distributed. We used the kinetic multi-layer model for aerosol surface and bulk chemistry (KM-SUB) to simulate ozone reaction constrained with our observationsmore »and inferred key parameters as a function of RH including Henry's Law constant for ozone, H O3 , and diffusion coefficients for ozone and iron, D O3 and D Fe , respectively. We found that H O3 is higher in our xanthan gum/FeCl 2 particles than for water and increases when RH decreased from about 80% to dry conditions. This coincided with a decrease in both D O3 and D Fe . In order to reproduce observed chemical gradients, our model predicted that ozone could not be present further than a few nanometers from a particle surface indicating near surface reactions were driving changes in iron oxidation state. However, the observed chemical gradients in α observed over hundreds of nanometers must have been the result of iron transport from the particle interior to the surface where ozone oxidation occurred. In the context of our results, we examine the applicability of the reacto-diffusive framework and discuss diffusion limitations for other reactive gas-aerosol systems of atmospheric importance.« less