Estimation of secondary organic aerosol viscosity from explicit modeling of gas-phase oxidation of isoprene and <i>α</i>-pinene

Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we apply the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model to conduct explicit chemical modeling of isoprene photooxidation and α-pinene ozonolysis and their subsequent SOA formation. The detailed gas-phase chemical schemes from GECKO-A are implemented into a box model and coupled to our recently developed glass transition temperature parameterizations, allowing us to predict SOA viscosity. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model, including autoxidation and particle-phase reactions, leading to the formation of high-molar-mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

Estimation of Secondary Organic Aerosol Viscosity from Explicit Modeling of Gas-Phase Oxidation of Isoprene and α-pinene

Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we conduct explicit modeling of isoprene photooxidation and α-pinene ozonolysis and subsequent SOA formation using the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model. Our recently-developed parameterizations to predict glass transition temperature of organic compounds are implemented into a box model with explicit gas-phase chemical mechanisms to simulate viscosity of SOA. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model including gas-phase dimerization and particle-phase reactions leading to the formation of high molar mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may also play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous 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.

Mass Accommodation and Gas-Particle Partitioning in Secondary Organic Aerosols: Dependence on Diffusivity, Volatility, Particle-phase Reactions, and Penetration Depth

Mass accommodation is an essential process for gas-particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface-bulk exchange, it can be important for viscous semisolid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas-particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs-Sutugin approximation of mass-transport kinetics at the gas-particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (KM-GAP, Shiraiwa et al., 2012) and two-film model solutions (MOSAIC, Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semisolid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in case of low-volatile compounds, whereas it can decrease by several orders of magnitude in case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs-Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models.

The influence of particle composition upon the evolution of urban ultrafine diesel particles on the neighbourhood scale

A recent study demonstrated that diesel particles in urban air undergo evaporative shrinkage when advected to a cleaner atmosphere (Harrison et al., 2016). We explore, in a structured and systematic way, the sensitivity of nucleation-mode diesel particles (diameter < 30 nm) to changes in particle composition, saturation vapour pressure, and the mass accommodation coefficient. We use a multicomponent aerosol microphysics model based on surrogate molecule (C16−C32 n-alkane) volatilities. For standard atmospheric conditions (298 K, 1013.25 hPa), and over timescales (ca. 100 s) relevant for dispersion on the neighbourhood scale (up to 1 km), the choice of a particular vapour pressure dataset changes the range of compounds that are appreciably volatile by two to six carbon numbers. The nucleation-mode peak diameter, after 100 s of model runtime, is sensitive to the vapour pressure parameterisations for particles with compositions centred on surrogate molecules between C22H46 and C24H50. The vapour pressure range, derived from published data, is between 9.23 × 10−3 and 8.94 × 10−6 Pa for C22H46 and between 2.26 × 10−3 and 2.46 × 10−7 Pa for C24H50. Therefore, the vapour pressures of components in this range are critical for the modelling of nucleation-mode aerosol dynamics on the neighbourhood scale and need to be better constrained. Laboratory studies have shown this carbon number fraction to derive predominantly from engine lubricating oil. The accuracy of vapour pressure data for other (more and less volatile) components from laboratory experiments is less critical. The influence of a core of non-volatile material is also considered; non-volatile core fractions of more than 5 % are inconsistent with the field measurements that we test the model against. We consider mass accommodation coefficient values less than unity and find that model runs with more volatile vapour pressure parameterisations and lower accommodation coefficients are similar to runs with less volatile vapour pressure parameterisations and higher accommodation coefficients. The new findings of this study may also be used to identify semi-volatile organic compound (SVOC) compositions that play dominating roles in the evaporative shrinkage of the nucleation mode observed in field measurements (Dall'Osto et al., 2011).

Studying volatility from composition, dilution, and heating measurements of secondary organic aerosols formed during <i>α</i>-pinene ozonolysis

Traditional yield curve analysis shows that semi-volatile organic compounds are a major component of secondary organic aerosols (SOAs). We investigated the volatility distribution of SOAs from α-pinene ozonolysis using positive electrospray ionization mass analysis and dilution- and heat-induced evaporation measurements. Laboratory chamber experiments were conducted on α-pinene ozonolysis, in the presence and absence of OH scavengers. Among these, we identified not only semi-volatile products, but also less volatile highly oxygenated molecules (HOMs) and dimers. Ozonolysis products were further exposed to OH radicals to check the effects of photochemical aging. HOMs were also formed during OH-initiated photochemical aging. Most HOMs that formed from ozonolysis and photochemical aging had 10 or fewer carbons. SOA particle evaporation after instantaneous dilution was measured at < 1 and ∼ 40 % relative humidity. The volume fraction remaining of SOAs decreased with time and the equilibration timescale was determined to be 24–46 min for SOA evaporation. The experimental results of the equilibration timescale can be explained when the mass accommodation coefficient is assumed to be 0.1, suggesting that the existence of low-volatility materials in SOAs, kinetic inhibition, or some combined effect may affect the equilibration timescale measured in this study.

Estimates of the organic aerosol volatility in a boreal forest using two independent methods

The volatility distribution of secondary organic aerosols that formed and had undergone aging – i.e., the particle mass fractions of semi-volatile, low-volatility and extremely low volatility organic compounds in the particle phase – was characterized in a boreal forest environment of Hyytiälä, southern Finland. This was done by interpreting field measurements using a volatility tandem differential mobility analyzer (VTDMA) with a kinetic evaporation model. The field measurements were performed during April and May 2014. On average, 40 % of the organics in particles were semi-volatile, 34 % were low-volatility organics and 26 % were extremely low volatility organics. The model was, however, very sensitive to the vaporization enthalpies assumed for the organics (ΔHVAP). The best agreement between the observed and modeled temperature dependence of the evaporation was obtained when effective vaporization enthalpy values of 80 kJ mol−1 were assumed. There are several potential reasons for the low effective enthalpy value, including molecular decomposition or dissociation that might occur in the particle phase upon heating, mixture effects and compound-dependent uncertainties in the mass accommodation coefficient. In addition to the VTDMA-based analysis, semi-volatile and low-volatility organic mass fractions were independently determined by applying positive matrix factorization (PMF) to high-resolution aerosol mass spectrometer (HR-AMS) data. The factor separation was based on the oxygenation levels of organics, specifically the relative abundance of mass ions at m∕z 43 (f43) and m∕z 44 (f44). The mass fractions of these two organic groups were compared against the VTDMA-based results. In general, the best agreement between the VTDMA results and the PMF-derived mass fractions of organics was obtained when ΔHVAP = 80 kJ mol−1 was set for all organic groups in the model, with a linear correlation coefficient of around 0.4. However, this still indicates that only about 16 % (R2) of the variation can be explained by the linear regression between the results from these two methods. The prospect of determining of extremely low volatility organic aerosols (ELVOAs) from AMS data using the PMF analysis should be assessed in future studies.

The effect of viscosity and diffusion on the HO₂ uptake by sucrose and secondary organic aerosol particles

We report the first measurements of HO2 uptake coefficients, γ, for secondary organic aerosol (SOA) particles and for the well-studied model compound sucrose which we doped with copper(II). Above 65 % relative humidity (RH), γ for copper(II)-doped sucrose aerosol particles equalled the surface mass accommodation coefficient α = 0.22 ± 0.06, but it decreased to γ = 0.012 ± 0.007 upon decreasing the RH to 17 %. The trend of γ with RH can be explained by an increase in aerosol viscosity and the contribution of a surface reaction, as demonstrated using the kinetic multilayer model of aerosol surface and bulk chemistry (KM-SUB). At high RH the total uptake was driven by reaction in the near-surface bulk and limited by mass accommodation, whilst at low RH it was limited by surface reaction. SOA from two different precursors, α-pinene and 1,3,5-trimethylbenzene (TMB), was investigated, yielding low uptake coefficients of γ < 0.001 and γ = 0.004 ± 0.002, respectively. It is postulated that the larger values measured for TMB-derived SOA compared to α-pinene-derived SOA are either due to differing viscosity, a different liquid water content of the aerosol particles, or an HO2 + RO2 reaction occurring within the aerosol particles.