skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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 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.  more » « less
Award ID(s):
1654104
PAR ID:
10252229
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Atmospheric Chemistry and Physics
Volume:
21
Issue:
8
ISSN:
1680-7324
Page Range / eLocation ID:
6175 to 6198
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models. 
    more » « less
  2. 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 insights 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. 
    more » « less
  3. Because traditional paleofire archives (e.g., burn scars on trees, charcoal in lake sediments) are not available in all settings, new ways of reconstructing past fire activity are needed. We focus here on polycyclic aromatic hydrocarbons (PAHs) in stalagmites. PAHs are organic molecules composed of two or more fused aromatic rings formed through incomplete combustion of organic matter, and vary in molecular weight depending on combustion conditions. Because the use of PAHs in stalagmites as a paleofire indicator is still in its infancy and because the production, deposition, and transport of PAHs into a cave is a complex and multi-faceted system, we tested the reproducibility of PAHs in two coeval and precisely-dated aragonite stalagmites – KNI-51-F and KNI-51-G - from KNI-51 (15.3°S, 128.6°E), a shallow cave located in the Kimberley region of tropical Western Australia. KNI-51-F and KNI-51-G span 1110-1620 CE and 1310-1630 CE, respectively. Each was hand-milled for analysis in continuous sections spanning approx. 2 mm-tall intervals at Ca’ Foscari University. Owing to differences in growth rate, temporal resolutions for KNI-51-F and KNI-51-G were 3±2 and 1±0.4 yr/sample, respectively. Chemical preparations and analysis methods follow those of Argiriadis et al. (2019) Analytical Chemistry, volume 91. In order to assess replication between the two stalagmites, we compared total abundances of low molecular weight (LMW: Napthalene, Acenaphthylene, Acenaphthene, Fluorene), medium molecular weight (MMW: Phenanthrene, Anthracene, Fluoranthene, Pyrene, Benzo(a)Anthracene, Chrysene, Retene), and high molecular weight (HHM: Benzo(b)Fluoranthene, Benzo(k)Fluoranthene, Benzo(e)Pyrene, Benzo(a)Pyrene, Perylene, Benzo(ghi)Perylene, Indeno(1,2,3-c,d)Pyrene, Dibenzo(A,H)Anthracene) PAHs. Total abundances of LMW, MMW, and HMW PAHs are similar (<10 ng/g) except for HMW PAHs in KNI-51-G, which are generally <1 ng/g. Total LMW and MMW abundance time series replicate well, with multiple synchronous multidecadal periods characterized by consistently low PAH abundances, suggestive of reduced bushfire activity, punctuated by intervals of high PAH abundances, likely reflecting frequent bushfire. Less coherence exists between HMW PAHs. 
    more » « 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 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. 
    more » « less
  5. Polycyclic aromatic hydrocarbons (PAHs) are organic molecules containing adjacent aromatic rings. Infrared emission bands show that PAHs are abundant in space, but only a few specific PAHs have been detected in the interstellar medium. We detected 1-cyanopyrene, a cyano-substituted derivative of the related four-ring PAH pyrene, in radio observations of the dense cloud TMC-1, using the Green Bank Telescope. The measured column density of 1-cyanopyrene is 1 .52×10 12 cm−2, from which we estimate that pyrene contains up to 0.1% of the carbon in TMC-1. This abundance indicates that interstellar PAH chemistry favors the production of pyrene. We suggest that some of the carbon supplied to young planetary systems is carried by PAHs that originate in cold molecular clouds. 
    more » « less