Single-particle measurements and estimations of activity coe ﬃ cients for semi-volatile organic compounds in organic aerosol of known chemical speciation †

Organic aerosol is a highly complex mixture of ∼ 10 4 to 10 5 unique compounds all possessing their own set of physico-chemical properties such as saturation vapour pressure and hygroscopicity. Most of these properties have not been experimentally measured and so must be estimated, resulting in large uncertainties in their predicted gas-particle partitioning and subsequent e ﬀ ect on human health and the climate. Dicarboxylic acids (DCAs) have been used extensively to represent oxidised organic aerosol due to their ubiquity in the atmosphere and commercial availability. Pure component saturation vapour pressures have been obtained by di ﬀ erent techniques to enable accurate treatment of their partitioning in atmospheric models. However, an understanding of the synergistic interactions between molecules in multi-component droplets containing DCAs under atmospherically relevant conditions remains limited, relying on group contribution models to estimate physicochemical properties such as activity coe ﬃ cients in complex organic multicomponent matrices. In this work, a method for extracting activity coe ﬃ cients of organic species in binary organic component droplets has been developed for investigating the e ﬀ ect of matrix character (such as functionality and molecular weight) on component volatility. Droplets containing a compound of interest and a low-volatility organic liquid are suspended in an electro-dynamic balance for hours to days and the evaporation rate is estimated as the droplet evolves in composition. Comparison with a liquid-like evaporation model allows for the extraction of experimental activity coe ﬃ cients for a range of mole fractions. Activity coe ﬃ cients extracted from simple binary organic systems have been compared to and have shown good agreement with an activity coe ﬃ cient estimator, Aerosol Inorganic – Organic Mixtures Functional groups Activity Coe ﬃ cients (AIOMFAC). The activity coe ﬃ cients of a range of DCAs dissolved in a mixed-component organic droplet have also been measured and compared to AIOMFAC. The changes in DCA activity coe ﬃ cient with evolving composition from experiments mostly show reasonable agreement with AIOMFAC for all systems investigated. This showcases the ability of AIOMFAC to predict activities for multi-organic systems containing no water.


I Introduction
Organic aerosol (OA) is ubiquitous in the atmosphere and is estimated to contribute a major fraction in mass to the total sub-micrometre sized aerosol in the Northern Hemisphere. 1 It is comprised of a vast number of species in the order 10 4 to 10 5 , all of which possess unique physicochemical properties, including pure component vapour pressures, refractive indices, and solubilities in water. 2 These properties combine, oen synergistically, to determine the hygroscopicity, optical scattering and gas-particle partitioning of semi-volatile components in aerosols of complex chemical composition, impacting on the climate, weather systems, visibility and health. 3For most atmospheric organic components, very few of the physicochemical properties have been experimentally determined for individual components.For a comprehensive exploration of the role of organic components in the atmosphere, they must instead be estimated, which can impart large uncertainties into their atmospheric impacts. [5] Equilibrium partitioning in ambient OA can be modelled in two ways, a top-down or bottom-up approach.Top-down approaches separate compositionally complex OA into volatility bins to generate volatility basis sets which can then be readily manipulated in models. 6Conversely, bottom-up models input concentrations of individual organic species and then estimate their partitioning based on their physicochemical properties, such as saturation vapour pressure, p sat .Bottom-up models assume that OA is comprised of homogeneous liquid droplets where the timescale for internal mixing is instantaneous relative to mass ux to and from the droplet.This means that, from eqn (1), the partitioning of organics is dependent on p i,sat and g i .It is therefore imperative to use accurate values of these properties, especially for compounds whose partitioning has the highest sensitivity to organic mass loading, specically compounds with similar C * i to C OA , i.e. semi-volatile organic compounds (SVOCs).2][13][14] These models are tted from literature data which are most abundant for small molecules with simple functionality and little atmospheric relevance.For compounds with increasing molecular weight and complexity of functionality, data are oen highly extrapolated or missing, resulting in poorer estimation of their p sat .
Increasing the complexity from single-to multi-component systems introduces the need for an additional term, the activity coefficient (g), which describes the heteromolecular interactivity.This term quantitatively describes the deviation from Raoult's law, as expressed in expression 2, where, for a compound "i", p i is the vapour pressure at the surface of a solution and x i is the mole fraction.
Many studies have measured the activity coefficient of water (g w ) for a range of different organics of varying functionalities and structures due to the importance of quantifying water uptake on to organic aerosol. 15,16In binary aqueous systems, organic activity coefficients can be calculated via the Gibbs-Duhem relationship, provided the water activity is correctly measured.Activity coefficients of solution species can be predicted using group contribution method-based models.One such model is the Aerosol Inorganic-Organic Mixtures Functional groups Activity Coefficients (AIOMFAC), which has parameterised the interactivity of different functional groups based on a large body of experimental data. 17It is also consistent with the Gibbs-Duhem relationship, meaning accurate predictions of water activity give accurate predictions of organic g. 18 Most of the organic species in the tting set of AIOMFAC are small with simple functionality, meaning some interactivity may not be well understood for species with more complex functionality, leading to higher uncertainty in the predicted activity coefficient.Additionally, organics in multi-organic systems in lower RH environments may exhibit strong interactions with one another.The effects of these organic-organic interactions for atmospherically relevant organic species must be explored experimentally for future inclusion in group contribution models.
Only a handful of studies, to the authors' knowledge, have directly measured g of the organic components in aerosol droplets, despite the importance of knowing departures from ideality when determining the equilibrium partitioning.In particular, there are even fewer studies that have systematically explored the dependence of the volatility of a single component on the character of the solubilizing organic matrix (e.g. chemical functionality, O : C ratio, molecular mass).Instead, most atmospheric models use estimated activity coefficients such as those generated by online group contribution estimators, e.g.AIOMFAC, for inclusion in box models, or simply assume g org values are unity. 17,19One such experimental study of g org was conducted by Saleh and Khlystov whereby an integrated volume method was used to determine the effect of polarity on the volatility of adipic acid in binary droplets. 20The particle phase was not measured and the internal mixing was inferred from the partitioning.Cappa et al. measured g org values of dicarboxylic acids (DCAs) in a multi-component solid-state sample using thermal-desorption mass spectrometry (TD-MS). 21The temperature was ramped resulting in DCAs evaporating in order of decreasing volatility.This technique quantied bulk volatility; however, the changing composition of the sample with successive component desorption means all DCA activity coef-cients are retrieved from the experimental data at different compositions.As component activity coefficients can be strongly dependent on composition, they cannot be directly compared with one another in these experiments and instead must be quoted for their own individual matrix compositions.This same reasoning challenges assumptions made when estimating the volatility of organic components in FIGAERO-CIMS measurements, especially for compounds likely to have activity coefficients signicantly different to unity. 22Most recently, Liu et al. used isothermal evaporation experiments to identify the change in partitioning of multi-component organic aerosol. 23econdary organic aerosol (SOA) was generated in a chamber from common atmospheric precursors and introduced to a range of atmospherically relevant seed particles.Bulk SOA "activity coefficients" were then inferred through comparison of experiment with a box model.The benet of this method is that it examines the partitioning of realistic samples on realistic seeds.However, it is not possible to draw conclusions on a molecular level about the factors affecting g org .Only bulk values for activities are calculated which are not applicable if the composition were to be adapted.
An improved understanding of the interactivity of a wide range of different species in well-described mixed-component organic systems at a molecular level is essential, such that results can be integrated directly into bottom-up box models more accurately.These experimental values can then be directly compared with estimations from tools such as AIOMFAC.In this work, a method has been developed for estimating g org for species dissolved in a low-volatility organic matrix in a spherical, homogeneous droplet.The evaporation rate of a binary organic droplet is measured in an electrodynamic balance (EDB) and compared to an ideal continuum ux model to yield organic activity coefficient (g org ) estimates.Using this method, the effect of functionality and isomerism on the interactivity of the matrix can be probed.Section II details the methodology to obtain g org from experimental data and model simulations.Section III then describes the sensitivities of the results to a range of experimental and model uncertainties.For the purposes of demonstrating the methodology and the uncertainties in Sections II and III, a binary organic droplet containing Tween 20 (T20) and diethylene glycol (DEG) has been used as a model system, where DEG is the volatile evaporating from the droplet.Section III also contains data from 4 different simple binary organic droplets used to validate the new method by comparison with AIOMFAC.The approach is then applied to mixtures of DCAs with T20 in Section IV, to explore synergistic interactions between typical atmospheric SVOCs and a complex organic matrix.
II Method for estimating g org in binary organic droplets

II.a Electrodynamic balance (EDB) measurements
A solution containing 2 organic components, one of which is signicantly more volatile than the other, and an appropriate solvent, such as ethanol or water, is made.The solution is dispensed using a MicroFab piezo-actuated micro-dispenser and suspended electrodynamically in an EDB in a dry nitrogen gas ow of controlled temperature.The solvent completely evaporates from the droplet in 1-2 seconds, leaving a binary organic droplet suspended in the EDB.The droplet is illuminated by a 532 nm laser and resultant elastically scattered light centred around 45°from the forwards direction is collected with a CCD camera and recorded for analysis of the radius.The light scattering pattern is referred to as the phase function.This methodology has been described in detail in previous work and will not be described in detail here. 24The EDB setup described has been used in previous studies to measure the vapour pressure of organics of a range of volatilities. 25he method used to estimate the radius, herein referred to as the "Mie peak tting method", compares the angular positions of the interference fringes in the phase function against a library of simulated fringe positions for a given refractive index or refractive index range, a method described by Davies et al. 26 For an evaporating droplet of changing n (i.e. a mixedcomponent droplet), the fringes are compared to a library of angles spanning a range of r and n at once, rather than assuming a xed n and correcting at a later stage as is done, for example, using the geometric optics approximation method in previous studies. 24,27Oen, several refractive index and radius combinations will provide sufficiently low error to be identied giving seemingly offset trends in evolving radius with time.In this case, the tted data can be constrained offline to report a single trend which most closely matches the expected evolving n.To demonstrate the importance of using a varied n approach, a comparison has been made between the xed and varied n Mie peak tting methods for an example binary droplet of T20 and DEG (outlined in Section I), with initial x DEG = 0.93, shown in Fig. 1.Compound structures are shown in Table 1.
Fig. 1(A) shows the xed n treatment in black and the varied n treatment in red.Using a xed n treatment to estimate r of a droplet of changing n will give an incorrect evaporation gradient, but with a reduced scatter in the evolution of radius with time.If Dn is sufficiently large, however, the radius can sometimes switch to an alternative trend line which has a better t for the xed value of n chosen.Fig. 1(B) shows the radius switching between 2 parallel offset radius tracks when using the xed n approach, illustrating the problem of using this approach to estimate the evolving radius of a droplet of changing n.By contrast, the varied n approach smoothly traverses between these 2 parallel tracks, giving a different, more reliable evaporation rate to that calculated from either of the 2 parallel tracks of the xed n approach.

II.b Simulating ideal evaporation in binary droplets
Simulations of evolving radius as a function of time have been performed using a numerical model based on the Maxwell equation for evaporative ux.Maxwell's equation is valid for spherical homogeneous droplets in the continuum regime, where evaporative cooling is negligible, meaning the temperature of the droplet surface is equal to the temperature of the surrounding gas.Eqn (3) combines the Maxwell equation with Raoult's law (eqn (2)) to give the mass ux of a component "i", dm i /dt, where r is the droplet radius, M i is the component molecular weight, D i is the component diffusion coefficient, T is the droplet and gas ow temperature, and p i,N is the partial pressure of "i" at an innite distance from the surface: Modied versions of this equation have been used to generate models in earlier studies. 28,29Eqn (3) can also be expressed in terms of the change in radius squared with time, dr 2 /dt, as shown in eqn (4).
The model requires an initial radius and composition, given in terms of component mole fractions, x i , along with component physicochemical properties required in eqn (3).The evaporative mass loss is then calculated and summed for all components in a designated timestep, and the droplet mass, density, radius and component mole fractions are updated.This process is then repeated for a given number of timesteps to generate an evolving mass, which can be converted to radius with the droplet density, r d .In this rst step in the analysis methodology, the interactivity is assumed to be ideal, meaning g is set to 1.For the same T20/DEG droplet measurement reported in Fig. 1, evolving r, n and r d vs. t have been simulated in Fig. 2.
It is it assumed that the rate of mass transport is slow compared to the rate of heat transport to the droplet surface and, thus, the droplet experiences negligible evaporative cooling.For evaporating "semi-volatile" species at room temperature and pressure, this is the case.It is assumed that there are no kinetic limitations to mass transport such as can be imposed by surface lms or viscosity, and there is no reactivity.Binary gas diffusion coefficients in nitrogen gas, D i,N 2 , are estimated using the Chapman-Enskog method described in Bilde et al.
1][32] r d is calculated assuming volume additivity in using the mass fraction mixing rule for density, given in eqn (6), where MF i and r i,melt are the mass fraction and melt density of component "i" respectively.This means that any deviation in mixed-component density from ideality has been omitted.This is unlikely to be a signicant factor for liquid droplets comprised of organics in which the strongest interactions are hydrogen bonds and would more likely be signicant for solid mixtures of salts.
Lastly, the droplet refractive index, n d , has been calculated using the molar refraction mixing rule, given in eqn (8).This requires the calculation of several pure and mixed-component properties, using mixing rules given in eqn ( 5)- (8).Eqn ( 5) is used to calculate the molar refraction of each component, R i , which can be calculated from the pure component properties n, M i and r i,melt .Eqn ( 7) and ( 8) are used to calculate the average molecular weight, M d , and molar refraction of the droplet, R d .Eqn ( 9) is simply a rearranged form of eqn (5), using droplet properties rather than pure component properties.While this mixing rule has been found to be the most appropriate in work by Cai et al., it requires r d which assumes volume additivity, as mentioned above. 33

II.c Estimation of g org
Following the rst coarse step in the analysis method described in II.b, estimation of g is performed using a bespoke Python script.The gradient of radius squared with time, dr 2 /dt, is  calculated for the experimental and model data as a function of time.Using the simulated x vs. t from the model, dr 2 /dt is then converted to a function of x.A polynomial is then used to t the experimental dr 2 /dt vs. x to remove the impact of optical radius uctuations.A polynomial order of 2 or 3 is sufficient to capture the shape of the data without inadvertently tting the optical radius uctuations.Using eqn (4), for a small change in x, dx, g can be estimated from the difference between experimental and ideally modelled dr 2 /dt at equivalent x, and, thus, compositions, as given in the expression below: The dependence of g on x can then be estimated for the whole compositional range accessed in the experiment by iterating forwards in time.g is estimated at a minimum of 100 nm radius intervals, double the Dr that the Mie peak tting technique (described in Section IIa) can resolve. 26This is demonstrated in Fig. 3 for the example T20/DEG droplet used throughout this section.
As the DEG evaporates from the droplet, the Dr required per Dx decreases (i.e. less DEG needs to evaporate to change its mole fraction).This results in adjacent g values being increasingly spaced in x at smaller x DEG .

III Sensitivity of estimated g org to uncertainties and validation of method with AIOMFAC
The uncertainty in the retrieval of activity coefficients and the corresponding mole fractions can be attributed to a combination of experimental and model uncertainties.We consider these individually using the T20/DEG droplet as a test system, before discussing the overall uncertainty.

III.a Experimental uncertainties
Firstly, there is a radius uncertainty of ±100 nm in the tting of the phase function with the Mie peak tting method. 26For the test system, this imparts an uncertainty in x which increases with decreasing x to >±9% for x < 0.5.This is shown for the T20/ DEG droplet in Fig. S1 in the ESI.† The larger the molecular weight of the evaporating compound compared to the matrix compound, and the larger the starting radius of the droplet, the lower the sensitivity of x to r.The uncertainty imparted on g from the uncertainty in optical r is only minor due to the 2 nd order polynomial t, which removes the uctuations in dr 2 /dt from optical radius uctuations.Secondly, there is an uncertainty in experimental temperature of ±1 K from the thermocouple measurement which imparts an uncertainty on the p sat of the evaporating organic used in the model through the Clausius-Clayperon relationship.Using an average literature DH vap for DEG of 63.3 kJ mol −1 , a ±1 K uncertainty in T corresponds to a 9% average uncertainty in p sat . 34This impacts the simulated evaporation rates of droplets (as per eqn (3)) which, in turn, impact estimations of g.The % uncertainty imparted on g by a constant uncertainty in T scales proportionally with DH vap , as shown in Fig, S2 in the ESI.†

III.b Model uncertainties
There are uncertainties in the properties used in the evaporation model which also impact the estimated value of g, including the component melt densities, calculated droplet density from volume additivity and estimated diffusion coefficients.To quantify the effects of these uncertainties on the uncertainty in g, sensitivity analyses have been conducted using reasonable expected uncertainties in the properties.Simulated evaporation curves for T20/DEG droplets, such as shown in Fig. 2, were generated using the evaporation model and the input properties were modied to assess their impact on g.An uncertainty in liquid/melt density will impart an uncertainty on the estimated g, through altering the modelled dr 2 /dt, and the simulated mole fraction for which the estimated g is assigned.
To quantify this, a ±5% uncertainty in DEG liquid density, r DEG , was imposed, and the resulting uncertainty in both g and x were identied as a function of x.The results are shown in Fig. S3 in the ESI.† The resulting uncertainty in x increases as x decreases, in a similar way to the impact of optical radius uncertainty.For x DEG < 0.5, where very small changes in radius correspond to large changes in x DEG , large decreases in calculated g DEG could instead be explained by an incorrect melt r DEG resulting in a lower modelled radius corresponding to complete DEG loss.As no composition information is retrieved, such as would be from Raman spectroscopy, the amount of DEG remaining in the droplet is not known.Thus, it is best to avoid estimating g DEG at radii that could feasibly correspond to near or complete loss of DEG (and low mole fractions).A ±5% uncertainty in liquid/melt r imparts a 5-7% uncertainty in g with a minor dependence on x.
The uncertainty arising from deviation from ideality in the estimate of the mixed-component density has been investigated.For compounds that have similar molecular weights (M i ), the maximum deviation in density from ideality would be expected when there are equal amounts of both compounds. 33owever, some systems will have components with a signicantly different M i to one another, meaning that there is a signicant difference in mass fraction (MFS)-dependent or xdependent non-ideality.Additionally, binary systems comprised of signicantly different sized molecules are likely to exhibit higher density than predicted by assuming volume additivity.To investigate this, x and MFS-dependent non-ideality were investigated separately.In both systems, a maximum of 5% higher density than volume additivity was considered, which is close to the maximum departure from volume additivity of aqueous sodium chloride (∼6%).The maximum resulting uncertainties in g and x from both experiments were ±5% and 0.4% respectively, which are signicantly lower than the experimental uncertainties, and so can be ignored.
An uncertainty of ±6% in gas diffusion coefficient, D i , typical of the error from the Chapman-Enskog method, imparts an almost equivalent uncertainty on g (±6%) but imparts no uncertainty in x.In the same way, g depends upon the p sat value used in the model, which is not always possible to measure directly with this methodology.In the case where it is not possible to measure the pure component p sat values with this technique, values must be taken from the literature, which have quoted associated uncertainties.These uncertainties can be quite large, especially for semi-volatile species that are solid at room temperature which, in the example of DCAs, can range from ±20% for glutaric acid to ±75% for some functionalised DCAs.In this case, the uncertainty in literature p sat dominates the resulting uncertainty in g.

III.c Total uncertainty
The uncertainty in x is dominated by the uncertainties in optical radius and liquid/melt densities of the components, with increase as x decreases in both cases.Liquid/melt densities of species that have been parameterised from aqueous solutions typically have an uncertainty of <1%, which would result in a much lower uncertainty in x than calculated in this study.The uncertainty in g is derived from uncertainties in T (through its impact on p sat ), D i , r d and optical r.The optical r and r d have the smallest impact on the uncertainty of g.If the uncertainty in pure component vapour pressure is also small, the uncertainties in D i and T determine the overall uncertainty in g at ∼6-15% each.However, p sat values of atmospherically relevant semi-volatile at 298 K are rarely quoted to within ±15% uncertainty in the literature, as mentioned in the previous subsection.This means that in most cases, the uncertainty in g will be dominated by the uncertainty in pure component p sat .

III.d Validation of method by comparison with AIOMFAC
To show that this method can be used to accurately retrieve activity coefficients of organics, four binary organic systems of different functionality and O : C ratio were tested and the results were compared with AIOMFAC.Systems were chosen that are likely to be well characterised by AIOMFAC, such that AIOMFAC could be assumed to give correct values.The volatile compound was the same for all binary droplets and the matrix compound was changed (see Table 1).Although DEG has been used as the volatile compound for the previous analyses, its p sat is within an order of magnitude to some of the proposed matrix compounds, meaning the assumption of negligible matrix compound evaporation would be invalid.A volatile organic such as ethanol is volatile enough to lead to signicant evaporative cooling at 298 K, and uncertainties in evaporative cooling would lead to greater uncertainties in the simulated evaporation rate and subsequent estimation of g, thus ethanol was not chosen either.Instead, carbitol was selected to be the volatile for these experiments.With a p sat over an order of magnitude larger than all matrix compounds, carbitol can be assumed to be responsible for all mass loss in a carbitol/matrix droplet, without causing noticeable evaporative cooling, meaning it was a perfect candidate for the volatile compound.All matrix compounds are either in the drop-down menu of AIOMFAC or are chemically similar to compounds that are.All matrix compounds are liquids at room temperature with known pure refractive indices and densities, removing the uncertainties associated with these properties.To avoid an uncertainty in carbitol vapour pressure from temperature uctuations, a droplet of pure carbitol was injected prior to every droplet of carbitol/matrix solution.The vapour pressure of the pure carbitol was then calculated using the Maxwell equation in terms of radius squared, given in eqn (4), where dr 2 /dt is constant for a pure component droplet.A similar comparative kinetics approach is commonly used in EDB-based hygroscopicity experiments, where the evaporation kinetics of a water droplet is used to infer the relative humidity of the gas ow prior to the injection of a solute of unknown hygroscopicity. 35In the case of all experimental droplets in this study, p i,N was assumed to be 0 due to the dry N 2 gas ow removing organics from the surface of the droplet.Additionally, the uncertainty in carbitol diffusion coefficient is also removed by this comparative kinetics approach, as it is kept constant for both pure and mixed-component calculations in eqn ( 3) and ( 4).During the series of experiments, the extracted p sat for the pure carbitol varied from 16-18 Pa, which agrees very well with the literature value of 17 Pa. 36Table 1 contains the structure of the volatile, carbitol, and the different matrix compounds used in this set of experiments, glycerol, triethyl citrate, pentaethylene glycol and Tween 20.
Fig. 4 shows the experimental estimations of g vs. x for carbitol in carbitol/matrix droplets, compared to AIOMFAC predictions.For clarity, data points have been selected every 500 nm rather than every 100 nm.The experimental estimations of g vs. x show good agreement with AIOMFAC for all systems tested, within the combined experimental and model error.Perhaps surprisingly, the experimental g vs. x of carbitol in T20 (Fig. 4(D)) also shows good agreement with AIOMFAC, despite the high molecular weight and complex functionality of T20.This supports the use of AIOMFAC to predict the interactivity of complex molecules in the atmosphere.Overall, these data validate the method for estimation of organic activity coefficients in binary organic droplets.

IV Activity coefficients of dicarboxylic acids (g DCA ) dissolved in a Tween 20 matrix
To demonstrate the use of this technique for atmospherically relevant species, a selection of binary organic droplets containing Tween 20 and a dicarboxylic acid, herein referred to as "T20/DCA droplets", have been investigated and g org has been estimated for a range of mole fractions.The range of DCAs investigated is shown in Table 2, along with the properties used in the model simulations.Estimations of g org were also compared to AIOMFAC.Tween 20 was selected due to its liquid phase state at room temperature and pressure, its low viscosity of ∼400 mPa s −1 and its low volatility of <10 −6 Pa (from pure component EDB measurements).Thus, its contribution to droplet evaporation is negligible and it does not kinetically inhibit the evaporation of the DCA.As it is non-volatile, its values for diffusion coefficient and vapour pressure were not included in Table 2. T20 is also soluble in solvents with a range of polarities, meaning a large range of organic compounds with varying polarity can be dissolved in this matrix and studied using this technique.T20 also contains atmospherically relevant functional groups such as esters, hydrocarbon chains and hydroxyl groups, as seen in Table 1.This justies its use as an atmospheric analogue despite not being naturally present in the atmosphere.DCAs were chosen due to their atmospheric abundance, commercial availability and known physicochemical properties such as p sat .
The gas ow was dry nitrogen at 298.15 ± 1 K.The RH was assumed to be 0% but was not measured directly due to the inaccuracy of humidity probes and the lack of a suitable probe droplet composition, oen used in our EDB studies, which can be adopted under these conditions. 37It was not possible to measure the p sat values of the pure DCAs in the same way as was described in the section IIId, as pure DCA particles are solid at 298.15 K, and thus do not produce the regular interference fringes needed for radius estimation.Instead, melt p sat values were taken from Bilde et al. (2015) and pure component melt r and n values were extrapolated from bulk solution data as described in Cai et al. (2016). 8,33otal uncertainty was calculated using the sensitivity analysis described in Section III, using compound-specic uncertainties. 33The resulting uncertainty in AIOMFAC predictions from the uncertainty in temperature of ±1 K is <1% at its maximum.
The maximum x DCA studied in these experiments for each DCA corresponds to the maximum amount of DCA while keeping the droplet homogeneous, i.e. the experimentally determined solubility limit.Measurements for succinic acid were also attempted but droplets of homogeneous morphology were not formed at any composition, preventing inference of activity coefficients by this approach.This may be due to its high lattice enthalpy which aids its efflorescence.
Fig. 5 compares experimental estimations and AIOMFAC predictions of g vs. x for 5 different DCAs in T20/DCA droplets.The directional trends of decreasing g with decreasing x match AIOMFAC predictions, however, the magnitude of this decrease is larger than predicted by AIOMFAC for the malonic acids at high mole fractions of DCA.Experimental estimations of glutaric acid show highest agreement with AIOMFAC, which is unsurprising considering its extensive coverage in the literature.
As x DCA approaches 1, i.e. the droplet tends towards pure dicarboxylic acid, g should also tend towards 1.For glutaric and methylglutaric acid, this is the case, indicating good agreement of these experiments with the literature p sat values used.However, for malonic, methylmalonic and methylsuccinic acid, g tends towards a value greater than 1, which may indicate a higher pure component p sat than agreed upon in the literature.Extrapolating g to x = 1 for the malonic acids gives a value of ∼10, which would correspond to an order of magnitude higher p sat than the literature.Although the p sat values have not been directly measured in this experiment, they have been agreed upon by several different studies at 298 K with a quoted uncertainty much lower than an order of magnitude.Disagreement between p sat values inferred from this study and the literature at the same temperature are unlikely to amount to an order of magnitude, meaning the disagreement between the results in this study and AIOMFAC at higher x cannot be fully accounted for by p sat .Another possibility is that the interactivity of T20 and the malonic acids are not well understood by AIOMFAC at high x DCA for these species.For the malonic acids, experimentally estimated DCA activities become greater than 1, indicating that phase separation would be thermodynamically favourable, however, phase separation is not predicted by AIOMFAC at any x.
Overall, DCAs show reasonable agreement with AIOMFAC at lower x DCA , but show deviations both in absolute values of g and changes in g with x at high x for malonic acids.The reasons for this need to be investigated further in future work.

V Conclusions and future work
This work provides a methodology for measuring g org for an organic species dissolved in a homogeneous mixed-component organic droplet.Binary organic droplets are suspended in an electrodynamic balance and their evolving radius is measured as the droplet evaporates.The activity coefficient can then be calculated as the composition of the droplet changes through comparison with an ideal liquid-like evaporation model.Several simple binary organic systems, deemed as well understood by AIOMFAC, were tested with this method and showed good agreement with AIOMFAC predictions of activity coefficient.To demonstrate this method for atmospherically relevant OA species, g DCA was estimated for dicarboxylic acids dissolved in Tween 20 and estimations were compared to AIOMFAC predictions.In general, AIOMFAC correctly predicted the direction and magnitude of the change in activity coefficient with evolving composition, particularly at lower x, however deviations in the experimental values from AIOMFAC were observed at higher x for malonic acids.Further investigation is needed to understand the reason for these observations.Future work will expand the coverage of different functionalities and isomers to investigate the complex relationship between g org and the chemical character of the matrix in different atmospherically relevant conditions, and investigate the potential for this methodology to be used for >2 organics per droplet.

Fig. 1
Fig. 1 Experimental T20/DEG droplet evaporation at 298.15 K and 0% RH, analysed with the Mie peak fitting method with a varied n (red scatter) and a constant n (black scatter), where (A) shows the fitted n vs. r and (B) shows the corresponding fitted r vs. t.

Fig. 2
Fig. 2 Simulated droplet r, Dn and r d vs. t for a droplet of T20/DEG at 298.15 K and 0% RH.

Fig. 3
Fig. 3 Demonstration of the estimation of g through comparison of experimental and simulated evaporation of a T20/DEG droplet at 298.15 K and 0% RH. (A) Experimental (scatter), polynomial fitted experimental (dashed line) and modelled (solid line) evaporation, dr 2 /dt, as a function of x DEG .(B) The resulting estimated relationship of g vs. x for DEG, using eqn (9), using the fitted experimental data (dashed line).

Table 1
Compounds used in binary carbitol/matrix droplet experiments

Table 2
Properties of DCAs and Tween 20 used in the model to simulate the evolving droplet radius with time