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

Efflorescence of Ammonium Sulfate and Coated Ammonium Sulfate Particles: Evidence for Surface Nucleation

Using optical microscopy, we investigated the efflorescence of ammonium sulfate (AS) in aqueous AS and in aqueous 1:1 and 8:1 (by dry weight) poly(ethylene glycol)-400 (PEG-400)/AS particles deposited on a hydrophobically coated slide. Aqueous PEG-400/AS particles exposed to decreasing relative humidity (RH) exhibit a liquid−liquid phase separation below 90% RH with the PEG-400-rich phase surrounding the aqueous AS inner phase. Pure aqueous AS particles effloresced in the RH range from 36.3% to 43.7%, in agreement with literature data (31−48% RH). In contrast, aqueous 1:1 (by dry weight) PEG-400/AS particles with diameters of the AS phase from 7.2 to 19.2 μm effloresced between 26.8% and 33.9% RH and aqueous 8:1 (by dry weight) PEG-400/AS particles with diameters of the AS phase from 1.8 to 7.3 μm between 24.3% and 29.3% RH. Such low efflorescence relative humidity (ERH) values have never been reached before for AS particles of this size range. We show that these unprecedented low ERHs of AS in PEG-400/AS particles could not possibly be explained by the presence of low amounts of PEG-400 in the aqueous AS phase, by a potential inhibition of water evaporation via anomalously slow diffusion through the PEG coating, or by different time scales between various experimental techniques. High-speed photography of the efflorescence process allowed the development of the AS crystallization fronts within the particles to be monitored with millisecond time resolution. The nucleation sites were inferred from the initial crystal growth sites. Analysis of the probability distribution of initial sites of 31 and 19 efflorescence events for pure AS and 1:1 (by dry weight) PEG-400/AS particles, respectively, showed that the particle volume can be excluded as the preferred nucleation site in the case of pure AS particles. For aqueous 1:1 (by dry weight) PEG-400/AS particles preferential AS nucleation in the PEG phase and at the PEG/AS/substrate contact line can be excluded. On the basis of this probability analysis of efflorescence events together with the AS ERH values of pure aqueous AS and aqueous PEG-400/AS particles aforementioned, we suggest that in pure aqueous AS particles nucleation starts at the surface of the particles and attribute the lower ERH values observed for aqueous PEG-400/AS particles to the suppression of the surface-induced nucleation process. Our results suggest that surface-induced nucleation is likely to also occur during the efflorescence of atmospheric AS aerosol particles, possibly constituting the dominating nucleation pathway

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

Experimental evidence for excess entropy discontinuities in glass-forming solutions

Glass transition temperatures T_g are investigated in aqueous binary and multi-component solutions consisting of citric acid, calcium nitrate (Ca(NO_3)_2), malonic acid, raffinose, and ammonium bisulfate (NH_4HSO_4) using a differential scanning calorimeter. Based on measured glass transition temperatures of binary aqueous mixtures and fitted binary coefficients, the T_g of multi-component systems can be predicted using mixing rules. However, the experimentally observed T_g in multi-component solutions show considerable deviations from two theoretical approaches considered. The deviations from these predictions are explained in terms of the molar excess mixing entropy difference between the supercooled liquid and glassy state at T_g. The multi-component mixtures involve contributions to these excess mixing entropies that the mixing rules do not take into account

Palladium-Catalyzed Oxidation of β-C(sp3)-H Bonds of Primary Alkylamines through a Rare Four-Membered Palladacycle Intermediate.

Site-selective functionalizations of C-H bonds are often achieved with a directing group that leads to five- or six-membered metallacyclic intermediates. Analogous reactions that occur through four-membered metallacycles are rare. We report a challenging palladium-catalyzed oxidation of primary C-H bonds β to nitrogen in an imine of an aliphatic amine, a process that occurs through a four-membered palladacyclc intermediate. The success of the reaction relies on the identification, by H/D exchange, of a simple directing group (salicylaldehyde) capable of inducing the formation of this small ring. To gain insight into the steps of the catalytic cycle of this unusual oxidation reaction, a series of mechanistic experiments and density functional theory (DFT) calculations were conducted. The experimental studies showed that cleavage of the C-H bond is rate-limiting and formation of the strained four-membered palladacycle is thermodynamically uphill. DFT calculations corroborated these conclusions and suggested that the presence of an intramolecular hydrogen bond between the oxygen of the directing group and hydroxyl group of the ligating acetic acid is crucial for stabilization of the palladacyclic intermediate

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

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