Modeling the gas-particle partitioning of secondary organic aerosol: the importance of liquid-liquid phase separation

The partitioning of semivolatile organic compounds between the gas phase and aerosol particles is an important source of secondary organic aerosol (SOA). Gas-particle partitioning of organic and inorganic species is influenced by the physical state and water content of aerosols, and therefore ambient relative humidity (RH), as well as temperature and organic loading levels. We introduce a novel combination of the thermodynamic models AIOMFAC (for liquid mixture non-ideality) and EVAPORATION (for pure compound vapor pressures) with oxidation product information from the Master Chemical Mechanism (MCM) for the computation of gas-particle partitioning of organic compounds and water. The presence and impact of a liquid-liquid phase separation in the condensed phase is calculated as a function of variations in relative humidity, organic loading levels, and associated changes in aerosol composition. We show that a complex system of water, ammonium sulfate, and SOA from the ozonolysis of α-pinene exhibits liquid-liquid phase separation over a wide range of relative humidities (simulated from 30% to 99% RH). Since fully coupled phase separation and gas-particle partitioning calculations are computationally expensive, several simplified model approaches are tested with regard to computational costs and accuracy of predictions compared to the benchmark calculation. It is shown that forcing a liquid one-phase aerosol with or without consideration of non-ideal mixing bears the potential for vastly incorrect partitioning predictions. Assuming an ideal mixture leads to substantial overestimation of the particulate organic mass, by more than 100% at RH values of 80% and by more than 200% at RH values of 95%. Moreover, the simplified one-phase cases stress two key points for accurate gas-particle partitioning calculations: (1) non-ideality in the condensed phase needs to be considered and (2) liquid-liquid phase separation is a consequence of considerable deviations from ideal mixing in solutions containing inorganic ions and organics that cannot be ignored. Computationally much more efficient calculations relying on the assumption of a complete organic/electrolyte phase separation below a certain RH successfully reproduce gas-particle partitioning in systems in which the average oxygen-to-carbon (O:C) ratio is lower than ~0.6, as in the case of α-pinene SOA, and bear the potential for implementation in atmospheric chemical transport models and chemistry-climate models. A full equilibrium calculation is the method of choice for accurate offline (box model) computations, where high computational costs are acceptable. Such a calculation enables the most detailed predictions of phase compositions and provides necessary information on whether assuming a complete organic/electrolyte phase separation is a good approximation for a given aerosol system. Based on the group-contribution concept of AIOMFAC and O:C ratios as a proxy for polarity and hygroscopicity of organic mixtures, the results from the α-pinene system are also discussed from a more general point of view

Computation of liquid-liquid equilibria and phase stabilities: implications for RH-dependent gas/particle partitioning of organic-inorganic aerosols

Semivolatile organic and inorganic aerosol species partition between the gas and aerosol particle phases to maintain thermodynamic equilibrium. Liquid-liquid phase separation into an organic-rich and an aqueous electrolyte phase can occur in the aerosol as a result of the salting-out effect. Such liquid-liquid equilibria (LLE) affect the gas/particle partitioning of the different semivolatile compounds and might significantly alter both particle mass and composition as compared to a one-phase particle. We present a new liquid-liquid equilibrium and gas/particle partitioning model, using as a basis the group-contribution model AIOMFAC (Zuend et al., 2008). This model allows the reliable computation of the liquid-liquid coexistence curve (binodal), corresponding tie-lines, the limit of stability/metastability (spinodal), and further thermodynamic properties of multicomponent systems. Calculations for ternary and multicomponent alcohol/polyol-water-salt mixtures suggest that LLE are a prevalent feature of organic-inorganic aerosol systems. A six-component polyol-water-ammonium sulphate system is used to simulate effects of relative humidity (RH) and the presence of liquid-liquid phase separation on the gas/particle partitioning. RH, salt concentration, and hydrophilicity (water-solubility) are identified as key features in defining the region of a miscibility gap and govern the extent to which compound partitioning is affected by changes in RH. The model predicts that liquid-liquid phase separation can lead to either an increase or decrease in total particulate mass, depending on the overall composition of a system and the particle water content, which is related to the hydrophilicity of the different organic and inorganic compounds. Neglecting non-ideality and liquid-liquid phase separations by assuming an ideal mixture leads to an overestimation of the total particulate mass by up to 30% for the composition and RH range considered in the six-component system simulation. For simplified partitioning parametrizations, we suggest a modified definition of the effective saturation concentration, C_j^*, by including water and other inorganics in the absorbing phase. Such a C_j^* definition reduces the RH-dependency of the gas/particle partitioning of semivolatile organics in organic-inorganic aerosols by an order of magnitude as compared to the currently accepted definition, which considers the organic species only

Liquid-liquid phase separation and morphology of internally mixed dicarboxylic acids/ammonium sulfate/water particles

Knowledge of the physical state and morphology of internally mixed organic/inorganic aerosol particles is still largely uncertain. To obtain more detailed information on liquid-liquid phase separation (LLPS) and morphology of the particles, we investigated complex mixtures of atmospherically relevant dicarboxylic acids containing 5, 6, and 7 carbon atoms (C5, C6 and C7) having oxygen-to-carbon atomic ratios (O:C) of 0.80, 0.67, and 0.57, respectively, mixed with ammonium sulfate (AS). With micrometer-sized particles of C5/AS/H_2O, C6/AS/H_2O and C7/AS/H_2O as model systems deposited on a hydrophobically coated substrate, laboratory experiments were conducted for various organic-to-inorganic dry mass ratios (OIR) using optical microscopy and Raman spectroscopy. When exposed to cycles of relative humidity (RH), each system showed significantly different phase transitions. While the C5/AS/H_2O particles showed no LLPS with OIR = 2:1, 1:1 and 1:4 down to 20% RH, the C6/AS/H_2O and C7/AS/H_2O particles exhibit LLPS upon drying at RH 50 to 85% and ~90%, respectively, via spinodal decomposition, growth of a second phase from the particle surface or nucleation-and-growth mechanisms depending on the OIR. This suggests that LLPS commonly occurs within the range of O:C < 0.7 in tropospheric organic/inorganic aerosols. To support the comparison and interpretation of the experimentally observed phase transitions, thermodynamic equilibrium calculations were performed with the AIOMFAC model. For the C7/AS/H_2O and C6/AS/H_2O systems, the calculated phase diagrams agree well with the observations while for the C5/AS/H_2O system LLPS is predicted by the model at RH below 60% and higher AS concentration, but was not observed in the experiments. Both core-shell structures and partially engulfed structures were observed for the investigated particles, suggesting that such morphologies might also exist in tropospheric aerosols

Improved AIOMFAC model parameterisation of the temperature dependence of activity coefficients for aqueous organic mixtures

This study presents a new, improved parameterisation of the temperature dependence of activity coefficients in the AIOMFAC (Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients) model applicable for aqueous as well as water-free organic solutions. For electrolyte-free organic and organic–water mixtures the AIOMFAC model uses a group-contribution approach based on UNIFAC (UNIversal quasi-chemical Functional-group Activity Coefficients). This group-contribution approach explicitly accounts for interactions among organic functional groups and between organic functional groups and water. The previous AIOMFAC version uses a simple parameterisation of the temperature dependence of activity coefficients, aimed to be applicable in the temperature range from ~ 275 to ~ 400 K. With the goal to improve the description of a wide variety of organic compounds found in atmospheric aerosols, we extend the AIOMFAC parameterisation for the functional groups carboxyl, hydroxyl, ketone, aldehyde, ether, ester, alkyl, aromatic carbon-alcohol, and aromatic hydrocarbon to atmospherically relevant low temperatures. To this end we introduce a new parameterisation for the temperature dependence. The improved temperature dependence parameterisation is derived from classical thermodynamic theory by describing effects from changes in molar enthalpy and heat capacity of a multi-component system. Thermodynamic equilibrium data of aqueous organic and water-free organic mixtures from the literature are carefully assessed and complemented with new measurements to establish a comprehensive database, covering a wide temperature range (~ 190 to ~ 440 K) for many of the functional group combinations considered. Different experimental data types and their processing for the estimation of AIOMFAC model parameters are discussed. The new AIOMFAC parameterisation for the temperature dependence of activity coefficients from low to high temperatures shows an overall improvement of 28% in comparison to the previous model version, when both versions are compared to our database of experimentally determined activity coefficients and related thermodynamic data. When comparing the previous and new AIOMFAC model parameterisations to the subsets of experimental data with all temperatures below 274 K or all temperatures above 322 K (i.e. outside a 25 K margin of the reference temperature of 298 K), applying the new parameterisation leads to 37% improvement in each of the two temperature ranges considered. The new parameterisation of AIOMFAC agrees well with a large number of experimental data sets. Larger model–measurement discrepancies were found particularly for some of the systems containing multi-functional organic compounds. The affected systems were typically also poorly represented at room temperature and further improvements will be necessary to achieve better performance of AIOMFAC in these cases (assuming the experimental data are reliable). The performance of the AIOMFAC parameterisation is typically better for systems containing relatively small organic compounds and larger deviations may occur in mixtures where molecules of high structural complexity such as highly oxygenated compounds or molecules of high molecular mass (e.g. oligomers) prevail. Nevertheless, the new parameterisation enables the calculation of activity coefficients for a wide variety of different aqueous/water-free organic solutions down to the low temperatures present in the upper troposphere

Corrigendum to "A thermodynamic model of mixed organic–inorganic aerosols to predict activity coefficients" published in Atmos. Chem. Phys., 8, 4559–4593, 2008

No abstract available

A thermodynamic model of mixed organic-inorganic aerosols to predict activity coefficients

Tropospheric aerosols contain mixtures of inorganic salts, acids, water, and a large variety of organic compounds. Interactions between these substances in liquid mixtures lead to discrepancies from ideal thermodynamic behaviour. By means of activity coefficients, non-ideal behaviour can be taken into account. We present here a thermodynamic model named AIOMFAC (Aerosol Inorganic-Organic Mixtures Functional groups Activity Coefficients) that is able to calculate activity coefficients covering inorganic, organic, and organic-inorganic interactions in aqueous solutions over a wide concentration range. This model is based on the activity coefficient model LIFAC by Yan et al. (1999) that we modified and reparametrised to better describe atmospherically relevant conditions and mixture compositions. Focusing on atmospheric applications we considered H^+, Li^+, Na^+, K^+, NH^+_4, Mg^(2+), Ca^(2+), Cl^−, Br^−, NO^−_3, HSO^−_4, and SO^(2−)_4 as cations and anions and a wide range of alcohols/polyols composed of the functional groups CH_n and OH as organic compounds. With AIOMFAC, the activities of the components within an aqueous electrolyte solution are well represented up to high ionic strength. Most notably, a semi-empirical middle-range parametrisation of direct organic-inorganic interactions in alcohol+water+salt solutions strongly improves the agreement between experimental and modelled activity coefficients. At room temperature, this novel thermodynamic model offers the possibility to compute equilibrium relative humidities, gas/particle partitioning and liquid-liquid phase separations with high accuracy. In further studies, other organic functional groups will be introduced. The model framework is not restricted to specific ions or organic compounds and is therefore also applicable for other research topics

Vapor pressures of substituted polycarboxylic acids are much lower than previously reported

The partitioning of compounds between the aerosol and gas phase is a primary focus in the study of the formation and fate of secondary organic aerosol. We present measurements of the vapor pressure of 2-methylmalonic (isosuccinic) acid, 2-hydroxymalonic (tartronic) acid, 2-methylglutaric acid, 3-hydroxy-3-carboxy-glutaric (citric) acid and DL-2,3-dihydroxysuccinic (DL-tartaric) acid, which were obtained from the evaporation rate of supersaturated liquid particles levitated in an electrodynamic balance. Our measurements indicate that the pure component liquid vapor pressures at 298.15 K for tartronic, citric and tartaric acids are much lower than the same quantity that was derived from solid state measurements in the only other room temperature measurement of these materials (made by Booth et al., 2010). This strongly suggests that empirical correction terms in a recent vapor pressure estimation model to account for the inexplicably high vapor pressures of these and similar compounds should be revisited, and that due caution should be used when the estimated vapor pressures of these and similar compounds are used as inputs for other studies

A thermodynamic framework for bulk–surface partitioning in finite-volume mixed organic–inorganic aerosol particles and cloud droplets

Atmospheric aerosol particles and their interactions with clouds are among the largest sources of uncertainty in global climate modeling. Aerosol particles in the ultrafine size range with diameters less than 100 nm have very high surface-area-to-volume ratios, with a substantial fraction of molecules occupying the air–droplet interface. The partitioning of surface-active species between the interior bulk of a droplet and the interface with the surrounding air plays a large role in the physicochemical properties of a particle and in the activation of ultrafine particles, especially those of less than 50 nm diameter, into cloud droplets. In this work, a novel and thermodynamically rigorous treatment of bulk–surface equilibrium partitioning is developed through the use of a framework based on the Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) model in combination with a finite-depth Guggenheim interface region on spherical, finite-volume droplets. We outline our numerical implementation of the resulting modified Butler equation, including accounting for challenging extreme cases when certain compounds have very limited solubility in either the surface or the bulk phase. This model, which uses a single, physically constrained interface thickness parameter, is capable of predicting the size-dependent surface tension of complex multicomponent solutions containing organic and inorganic species. We explore the impacts of coupled surface tension changes and changes in bulk–surface partitioning coefficients for aerosol particles ranging in diameters from several micrometers to as small as 10 nm and across atmospherically relevant relative humidity ranges. The treatment of bulk–surface equilibrium leads to deviations from classical cloud droplet activation behavior as modeled by simplified treatments of the Köhler equation that do not account for bulk–surface partitioning. The treatments for bulk–surface partitioning laid out in this work, when applied to the Köhler equation, are in agreement with measured critical supersaturations of a range of different systems. However, we also find that challenges remain in accurately modeling the growth behavior of certain systems containing small dicarboxylic acids, especially in a predictive manner. Furthermore, it was determined that the thickness of the interfacial phase is a sensitive parameter in this treatment; however, constraining it to a meaningful range allows for predictive modeling of aerosol particle activation into cloud droplets, including cases with consideration of co-condensation of semivolatile organics.</p

New and extended parameterization of the thermodynamic model AIOMFAC: calculation of activity coefficients for organic-inorganic mixtures containing carboxyl, hydroxyl, carbonyl, ether, ester, alkenyl, alkyl, and aromatic functional groups

We present a new and considerably extended parameterization of the thermodynamic activity coefficient model AIOMFAC (Aerosol Inorganic-Organic Mixtures Functional groups Activity Coefficients) at room temperature. AIOMFAC combines a Pitzer-like electrolyte solution model with a UNIFAC-based group-contribution approach and explicitly accounts for interactions between organic functional groups and inorganic ions. Such interactions constitute the salt-effect, may cause liquid-liquid phase separation, and affect the gas-particle partitioning of aerosols. The previous AIOMFAC version was parameterized for alkyl and hydroxyl functional groups of alcohols and polyols. With the goal to describe a wide variety of organic compounds found in atmospheric aerosols, we extend here the parameterization of AIOMFAC to include the functional groups carboxyl, hydroxyl, ketone, aldehyde, ether, ester, alkenyl, alkyl, aromatic carbon-alcohol, and aromatic hydrocarbon. Thermodynamic equilibrium data of organic-inorganic systems from the literature are critically assessed and complemented with new measurements to establish a comprehensive database. The database is used to determine simultaneously the AIOMFAC parameters describing interactions of organic functional groups with the ions H^+, Li^+, Na^+, K^+, NH_(4)^+, Mg^(2+), Ca^(2+), Cl^−, Br^−, NO_(3)^−, HSO_(4)^−, and SO_(4)^(2−). Detailed descriptions of different types of thermodynamic data, such as vapor-liquid, solid-liquid, and liquid-liquid equilibria, and their use for the model parameterization are provided. Issues regarding deficiencies of the database, types and uncertainties of experimental data, and limitations of the model, are discussed. The challenging parameter optimization problem is solved with a novel combination of powerful global minimization algorithms. A number of exemplary calculations for systems containing atmospherically relevant aerosol components are shown. Amongst others, we discuss aqueous mixtures of ammonium sulfate with dicarboxylic acids and with levoglucosan. Overall, the new parameterization of AIOMFAC agrees well with a large number of experimental datasets. However, due to various reasons, for certain mixtures important deviations can occur. The new parameterization makes AIOMFAC a versatile thermodynamic tool. It enables the calculation of activity coefficients of thousands of different organic compounds in organic-inorganic mixtures of numerous components. Models based on AIOMFAC can be used to compute deliquescence relative humidities, liquid-liquid phase separations, and gas-particle partitioning of multicomponent mixtures of relevance for atmospheric chemistry or in other scientific fields

Influence of particle-phase state on the hygroscopic behavior of mixed organic-inorganic aerosols

Recent work has demonstrated that organic and mixed organic–inorganic particles can exhibit multiple phase states depending on their chemical composition and on ambient conditions such as relative humidity (RH). To explore the extent to which water uptake varies with particle-phase behavior, hygroscopic growth factors (HGFs) of nine laboratory-generated, organic and organic–inorganic aerosol systems with physical states ranging from well-mixed liquids to phase-separated particles to viscous liquids or semi-solids were measured with the Differential Aerosol Sizing and Hygroscopicity Spectrometer Probe at RH values ranging from 40 to 90%. Water-uptake measurements were accompanied by HGF and RH-dependent thermodynamic equilibrium calculations using the Aerosol Inorganic-Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) model. In addition, AIOMFAC-predicted growth curves are compared to several simplified HGF modeling approaches: (1) representing particles as ideal, well-mixed liquids; (2) forcing a single phase but accounting for non-ideal interactions through activity coefficient calculations; and (3) a Zdanovskii–Stokes–Robinson-like calculation in which complete separation of the inorganic and organic components is assumed at all RH values, with water uptake treated separately in each of the individual phases. We observed variability in the characteristics of measured hygroscopic growth curves across aerosol systems with differing phase behaviors, with growth curves approaching smoother, more continuous water uptake with decreasing prevalence of liquid–liquid phase separation and increasing oxygen : carbon ratios of the organic aerosol components. We also observed indirect evidence for the dehydration-induced formation of highly viscous semi-solid phases and for kinetic limitations to the crystallization of ammonium sulfate at low RH for sucrose-containing particles. AIOMFAC-predicted growth curves are generally in good agreement with the HGF measurements. The performances of the simplified modeling approaches, however, differ for particles with differing phase states. This suggests that no single simplified modeling approach can be used to capture the water-uptake behavior for the diversity of particle-phase behavior expected in the atmosphere. Errors in HGFs calculated with the simplified models are of sufficient magnitude to produce substantial errors in estimates of particle optical and radiative properties, particularly for the assumption that water uptake is driven by absorptive equilibrium partitioning with ideal particle-phase mixing