Kinetic modeling of formation and evaporation of secondary organic aerosol from NO3 oxidation of pure and mixed monoterpenes

Organic aerosol constitutes a major fraction of the global aerosol burden and is predominantly formed as secondary organic aerosol (SOA). Environmental chambers have been used extensively to study aerosol formation and evolution under controlled conditions similar to the atmosphere, but quantitative prediction of the outcome of these experiments is generally not achieved, which signifies our lack in understanding of these results and limits their portability to large-scale models. In general, kinetic models employing state-of-the-art explicit chemical mechanisms fail to describe the mass concentration and composition of SOA obtained from chamber experiments. Specifically, chemical reactions including the nitrate radical (NO3) are a source of major uncertainty for assessing the chemical and physical properties of oxidation products. Here, we introduce a kinetic model that treats gas-phase chemistry, gas-particle partitioning, particle-phase oligomerization, and chamber vapor wall loss and use it to describe the oxidation of the monoterpenes a-pinene and limonene with NO3. The model can reproduce aerosol mass and nitration degrees in experiments using either pure precursors or their mixtures and infers volatility distributions of products, branching ratios of reactive intermediates and particle-phase reaction rates. The gas-phase chemistry in the model is based on the Master Chemical Mechanism (MCM) but trades speciation of single compounds for the overall ability of quantitatively describing SOA formation by using a lumped chemical mechanism. The complex branching into a multitude of individual products in MCM is replaced in this model with product volatility distributions and detailed peroxy (RO2) and alkoxy (RO) radical chemistry as well as amended by a particle-phase oligomerization scheme. The kinetic parameters obtained in this study are constrained by a set of SOA formation and evaporation experiments conducted in the Georgia Tech Environmental Chamber (GTEC) facility. For both precursors, we present volatility distributions of nitrated and non-nitrated reaction products that are obtained by fitting the kinetic model systematically to the experimental data using a global optimization method, the Monte Carlo genetic algorithm (MCGA). The results presented here provide new mechanistic insight into the processes leading to formation and evaporation of SOA. Most notably, the model suggests that the observed slow evaporation of SOA could be due to reversible oligomerization reactions in the particle phase. However, the observed non-linear behavior of precursor mixtures points towards a complex interplay of reversible oligomerization and kinetic limitations of mass transport in the particle phase, which is explored in a model sensitivity study. The methodologies described in this work provide a basis for quantitative analysis of multi-source data from environmental chamber experiments but also show that a large data pool is needed to fully resolve uncertainties in model parameters.