Influence of particle viscosity on mass transfer and heterogeneous ozonolysis kinetics in aqueous-sucrose-maleic acid aerosol.

Mass transfer between the gas and condensed phases in aerosols can be limited by slow bulk diffusion within viscous particles. During the heterogeneous and multiphase reactions of viscous organic aerosol particles, it is necessary to consider the interplay of numerous mass transfer processes and how they are impacted by viscosity, including the partitioning kinetics of semi-volatile organic reactants, water and oxidants. To constrain kinetic models of the heterogeneous chemistry, measurements must provide information on as many observables as possible. Here, the ozonolysis of maleic acid (MA) in ternary aerosol particles containing water and sucrose is used as a model system. By varying the mass ratio of sucrose to MA and by performing reactions over a wide range of relative humidity, direct measurements show that the viscosity of the particle can be varied over 7 orders of magnitude. Measurements of the volatilisation kinetics of MA show that this range in viscosity leads to a suppression in the effective vapour pressure of MA of 3-4 orders of magnitude. The inferred values of the diffusion coefficient of MA in the particle phase closely mirror the expected change in diffusion coefficient from the Stokes-Einstein equation and the change in viscosity. The kinetics of ozonolysis show a similar dependence on particle viscosity that can be further investigated using the kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB). Two scenarios, one constraining the diffusion coefficients for MA to those expected based on the Stokes-Einstein equation and the other including the diffusion coefficients as a fit parameter, yield similarly adequate representations of the ozonolysis kinetics, as inferred from the experimental decay in the signature of the vinylic C-H stretching vibration of MA. However, these two scenarios provide very different parameterisations of the compositional dependence of the diffusion coefficients of ozone within the condensed phase, yielding qualitatively different time-dependent internal concentration profiles. We suggest that this highlights the importance of providing additional experimental observables (e.g. particle size, heterogeneity in composition) if measurements and models are to be universally reconciled.

Isotherm-Based Thermodynamic Models for Solute Activities of Organic Acids with Consideration of Partial Dissociation.

Organic acids make up a significant fraction of the organic mass in atmospheric aerosol particles. The calculation of gas-liquid-solid equilibrium partitioning of the organic acid is therefore critical for accurate determination of atmospheric aerosol physicochemical properties and processes such as new particle formation and activation to cloud condensation nuclei. Previously, an adsorption isotherm-based statistical thermodynamic model was developed for capturing solute concentration-activity relationships for multicomponent aqueous solutions over the entire concentration range (Dutcher et al. J. Phys. Chem. C/A 2011, 2012, 2013), with model parameters for energies of adsorption successfully related to dipole-dipole electrostatic forces in solute-solvent and solvent-solvent interactions for both electrolytes and organics (Ohm et al. J. Phys. Chem. A 2015). However, careful attention is needed for weakly dissociating semivolatile organic acids. Dicarboxylic acids, such as malonic acid and glutaric acid are treated here as a mixture of nondissociated organic solute (HA) and dissociated solute (H(+) + A(-)). It was found that the apparent dissociation was greater than that predicted by known dissociation constants alone, emphasizing the effect of dissociation on osmotic and activity coefficient predictions. To avoid additional parametrization from the mixture approach, an expression was used to relate the Debye-Hückel hard-core collision diameter to the adjustable solute-solvent intermolecular distance. An improved reference state treatment for electrolyte-organic aqueous mixtures, such as that observed here with partial dissociation, has also been proposed. This work results in predictive correlations for estimation of organic acid and water activities for which there is little or no activity data.

Isotherm-based thermodynamic model for electrolyte and nonelectrolyte solutions incorporating long- and short-range electrostatic interactions.

The activities of solutes and solvents in solutions govern numerous physical phenomena in a wide range of practical applications. In prior work, we used statistical mechanics and multilayer adsorption isotherms to develop a transformative model for capturing thermodynamic properties of multicomponent aqueous solutions over the entire concentration range (Dutcher et al. J. Phys. Chem. 2011, 2012, 2013). That model needed only a few adsorption energy values to represent the solution thermodynamics of each solute. In the current work, we posit that the adsorption energies are due to dipole-dipole electrostatic forces in solute-solvent and solvent-solvent interactions. This hypothesis was tested in aqueous solutions on (a) 37 1:1 electrolytes, over a range of cation sizes, from H(+) to tetrabutylammonium, for common anions including Cl(-), Br(-), I(-), NO3(-), OH(-), ClO4(-), and (b) 20 water-soluble organic molecules including alcohols and polyols. For both electrolytes and organic solutions, the energies of adsorption can be calculated with the dipole moments of the solvent, molecular size of the solvent and solute, and the solvent-solvent and solvent-solute intermolecular bond lengths. Many of these physical properties are available in the literature, with the exception of the solute-solvent intermolecular bond lengths. For those, predictive correlations developed here enable estimation of solute and solvent solution activities for which there are little or no activity data.

Diffusion and reactivity in ultraviscous aerosol and the correlation with particle viscosity.

The slow transport of water, organic species and oxidants in viscous aerosol can lead to aerosol existing in transient states that are not solely governed by thermodynamic principles but by the kinetics of gas-particle partitioning. The relationship between molecular diffusion constants and particle viscosity (for example, as reflected in the Stokes-Einstein equation) is frequently considered to provide an approximate guide to relate the kinetics of aerosol transformation with a material property of the aerosol. We report direct studies of both molecular diffusion and viscosity in the aerosol phase for the ternary system water/maleic acid/sucrose, considering the relationship between the hygroscopic response associated with the change in water partitioning, the volatilisation of maleic acid, the ozonolysis kinetics of maleic acid and the particle viscosity. Although water clearly acts as a plasticiser, the addition of minor fractions of other organic moieties can similarly lead to significant changes in the viscosity from that expected for the dominant component forming the organic matrix (sucrose). Here we highlight that the Stokes-Einstein relationship between the diffusion constant of water and the viscosity of the particle may be more than an order of magnitude in error, even at viscosities as low as 1 Pa s. We show that the thermodynamic relationships of hygroscopic response that underpin such kinetic determinations must be accurately known to retrieve accurate values for diffusion constants; such data are often not available. Further, we show that scaling of the diffusion constants of organic molecules of similar size to those forming the matrix with particle viscosity may be well represented by the Stokes-Einstein equation, suppressing the apparent volatility of semi-volatile components. Finally, the variation in uptake coefficients and diffusion constants for oxidants and small weakly interacting molecules may be much less dependent on viscosity than the diffusion constants of more strongly interacting molecules such as water.