Characterising the Evaporation Kinetics of Water and SemiVolatile Organic Compounds from Viscous Multicomponent Organic Aerosol Particles. Physical Chemistry Physics

The physicochemical changes experienced by organic aerosol particles undergoing dehydration into the surrounding gas phase can be drastic, forcing rapid vitrification of the particle and suppressing internal diffusion. Until recently, experimental studies have concentrated on quantifying diffusional mixing of either water or non-volatile components, while relatively little attention has been paid to the role of semivolatile organic component (SVOC) diffusion and volatilisation in maintaining the equilibrium between the gas and particle phases. Here we present methods to simultaneously investigate diffusivities and volatilities in studies of evolving single ternary aerosol particle size and composition. Analysing particles of ternary composition must account for the multiple chemical species that volatilise in response to a step change in gas phase water activity. In addition, treatments of diffusion in multicomponent mixtures are necessary to represent evolving heterogeneities in particle composition. We find that the contributions to observed size behaviour from volatilisation of water and a SVOC can be decoupled and treated separately. Employing Fickian diffusion modelling, we extract the compositional dependence of the diffusion constant of water and compare the results to recently published parametrisations in binary aerosol particles. The treatment of ideality and activity

in each case is discussed, with reference to use in multicomponent core shell models. Meanwhile, the evaporation of a SVOC into an unsaturated gas flow may be treated by Maxwell's equation, with slow diffusional transport manifesting as a suppression in the extracted vapour pressure.

I. Introduction
In recent years, a diverse range of laboratory and ambient aerosol measurements have supported the existence of a moisture driven glass transition in secondary organic aerosol (SOA) particles. [1][2][3][4] Vitrification proceeds via the removal of water from the particle following a change in gas relative humidity (RH), which in turn leads to large increases in organic solute concentration close to the surface. 5,6 The increasing dehydration arrests molecular diffusion and confines the system to a metastable, amorphous 'solid-melt' type structure with no available kinetic pathway to undergo a phase transition and form an equilibrium crystal phase. 7, 8 Current work is focused on quantifying the influence of functional group chemistry on both the mesoscopic properties of single droplets (for example, the influence of aerosol composition on particle viscosity and phase), [9][10][11] and the diffusional dynamics and reactivity of molecules residing within viscous particles. 12,13 Such a fundamental understanding of the underlying physical chemistry may be essential to help qualify (or quantify) the environmental impacts of ensembles of glassy and semisolid and particles. atmospheric conditions over which these phenomena may lead to a slowing of particle phase mixing have begun to be investigated, [19][20][21] with some consensus emerging that diffusive limitations to water transport rarely occur in SOA in the planetary boundary layer. 22 In contrast, kinetically limited mass transport of SVOC in SOA could have ramifications for mass concentrations and size distributions of aerosol 23 and the lifetime of species in the particle phase. 24 Indeed, there is considerable interest from a more fundamental perspective on the relationship between diffusion on the molecular scale and viscosity on the mesoscale of particles. 8,15 The purpose of this paper is to explore methods for quantifying the competing contributions arising from the simultaneous evaporative mass transport of water and SVOCs in single-particle studies of viscous and glassy aerosol.
There are currently two single particle strategies by which the dependence of the diffusion constant of water on composition in viscous aerosol can be determined. In the first strategy, particles are exposed to step changes in water activity driving the formation of concentration gradients. The rate of water loss from the surface is extremely rapid; often, the core of the particle cannot respond rapidly enough to equalise the water activity before diffusion is effectively arrested near the surface. Such radial heterogeneity is well documented, 6,25 and modelling the subsequent internal water transport phenomena is central to extracting diffusion coefficients from experiments conducted on optically trapped and electrodynamically levitated particles. 16,26,27 Indeed, we have shown that care must be taken in interpreting kinetic measurements from radially inhomogeneous particles: the radial inhomogeneity at any moment in time is a complex function of the processing history of the particle and so very rarely even starts in a uniform homogeneous form. 28 To the best of our knowledge, most studies that assess the impact of aerosol microphysics on the composition of atmospheric aerosol do not account for this process history within an evolving gas phase environment, although some do account for diffusion using a core-shell model.
In a second strategy, the water activity is maintained at a constant value but the gas phase water is instantaneously switched from normal to deuterated water, leading to progressive isotopic exchange within the particle. 27,29 The change in isotopic ratio can then be used to extract the water activity dependence of the diffusion constant of water. This method has been used to study particles either deposited on a substrate or optically trapped. 30 Indeed, a similar approach has been extended to study the slow diffusional mixing of non-deuterated and deuterated sucrose between two contacting phases at low water activity, highlighting the conformity of sucrose diffusion constants to the SE equation. 18 More broadly, the diffusion constants of large organic molecules in viscous aerosol or through viscous shells on aerosol particles have been inferred from kinetic measurements of the rates of volatilisation of tracer compounds including, for example, pyrene. 24 In this publication, we will describe a framework for assessing the coupled evaporation kinetics of water and an SVOC from a viscous particle, allowing both a determination of the diffusion constant of water and the degree of suppression in the volatility of the SVOC with varying particle viscosity. We first review the framework for representing the compositional dependence of diffusion constants and the expected timescales for particle-gas partitioning of water and SVOCs in Section II. In Section III, we review the experimental methods used here. In Section IV we present a method for characterising the variation in particle size with time using the Kohlrausch-Williams-Watts equation, 31 separating out the timescales for water transport and for volatile loss of the SVOC. Further, we provide the first direct determinations of diffusion constants of water in ternary component aerosol. This benefits from the application of a new algorithm for routinely fitting multiple kinetic events to yield a compositionally dependent diffusion constant averaged over many events, a considerable step forward from previous work where computational demands restricted the analysis of only individual events. Finally, we examine measurements of the inferred loss rates of a prototypical SVOC (malonic acid) from viscous aerosol, exploring the variation in evaporation rate with moisture content and providing a qualitative examination of the dependence on particle viscosity.

II. The Compositional Dependence of Diffusion Constants and the Relative Timescales for SVOC and Water Partitioning
A comparison of the reported water activity dependencies of the diffusion constant of water inferred from isotope exchange and RH step experiments suggests that they are incompatible with qualitatively different shapes. 30 Recently, it has been common to adopt the Vignes parametrisation to express the diffusion constant of water in the mixture, Dw(xw), where it has a mole fraction of xw, in terms of the limiting diffusion constants of water in pure water, Dw,w, and at infinite dilution in pure solute, Dw,org: 12,16,27,30 ) . (1) Each of these studies involves fitting equation (1) to datasets of D values inferred from measurements at steady gas-phase RHs (i.e. constant water activities) in isotope exchange measurements, 27,30 or from D values inferred for a single particle during the transition in particle size and heterogeneity in water activity following a step change in gas phase RH. 16 Here, we explicitly include the Vignes dependence within a Fickian diffusion model to reproduce experimental data in which the water activity is changing within a particle following a step change in gas phase RH.
When defining the dependence of the diffusion constant on composition, as in equation (1), the dimensionless constant α is dependent on composition and is frequently referred to as an activity coefficient. 26,27 We have chosen to express this compositional dependence using a cubic dependence on mole fraction, recently introduced by Davies and Wilson 30 as a modification of the temperature dependent Margules type 32 binomial expansion of activity: C and D are fit coefficients for the different terms in the expansion. Although equation (3) is a formulation that was intended for use in treating binary solutions, 33 and its initial use in the parametrisation of aqueous sucrose activity, 34 it has shown some applicability for representing water diffusion in organic aerosol by accurately fitting water diffusion in citric acid 26 and highly complex secondary organic material. 27 Compared to the binomial expansion, the cubic parametrisation of α produces an almost identical shape for Dw(xw), in Eq (1), but with a much less pronounced plateau at the bulk water value of Dw at high water activity. Many techniques, including aerosol optical tweezers experiments, are insensitive to transitions within this region of activity, as the timescale for water transport is comparable to the timescale for the instrumentally limited change of RH. As such, it cannot be said with any confidence which of these parametrisations, (2) or (3), is a more realistic approximation to the binary diffusion coefficients in such particles. Instead, we are seeking the most accurate method to represent data recorded for highly dehydrated droplets that experience large kinetic limitations, and access compositions governed by the shape of the curve at much lower water mole fractions.
The utility, drawbacks, and recent uses of Fickian modelling for characterising diffusional transport in binary viscous aerosol have been discussed in previous publications. 31 We have chosen to use the recently developed partial differential model, Fi-PaD, described by O'Meara et al. 35 and a brief description is provided below. Fi-PaD represents the particle by a vector containing compositional information, with each element describing a radial shell within the particle bulk; here 400 shells are used, corresponding to a shell resolution of around 10 nm for typical particle radii in our single-particle measurements, 4 m. The model assumes an initial homogeneous concentration profile at equilibrium with a user specified gas phase RH, which can then undergo an immediate transition to another (or The simplified treatment of α (Equation 2) described above greatly decreases the computational expense of the data analysis as it removes one dimension from the fitting space when compared with fitting to the more complex dependence of equation (3). By iteratively varying log(Dw,org) and C it is possible to improve the overlay between the simulated radius responses and the data. This was achieved using the Matlab interior-point algorithm, within the fmincon function, starting from four randomly selected starting points in (Dw,org, C) space. The extent of agreement between each simulation and the experiment is quantified by calculation of the residual between the two response functions at each simulated timestep: where is the radial response function, as discussed in our previous publication. 31  Running on a single Intel i7-4790 core, simulation of a typical experimental RH step change (8,000-10,000 s in length) was found to take approximately 500 s, with typically around one hundred simulations required to converge the fit parameters and reproduce the measured response in size. The model was considered to have fit when all four independent fitting runs had found the nearest minimum in residual space, and the mean deviation from experiment was less than 3% of the magnitude of the experimental data. We provide an example of the fitting in Figure 1 for a ternary aqueous-sucroseglucose droplet experiencing a change in RH from 20 to 6 % RH. As in our previous work, we represent the relaxation in size following the RH step by a response function (see below) that is calculated from the time-dependence in wavelength of a whispering gallery mode (WGM) in the Raman spectrum. 31 The initial droplet size prior to the step in RH is calculated explicitly from the fingerprint of WGMs appearing in the Raman signature. The value of Dw,org is the limiting value of the diffusion coefficient of water in the fully dry limit of an isolated water molecule diffusing in a viscous matrix of sucrose and glucose present in a 1:1 mass ratio. The methods applied here offer improvements over similar diffusion models both in terms of accuracy (varying the fit parameters and parallelising over multiple CPU cores reduces the likelihood of becoming trapped in a local minimum in the parameter space) and throughput (see above).
Most studies establish experimental conditions such that the volatilisation of water or an SVOC can be studied in isolation, allowing the estimation of diffusion constants of water or the diffusion constants of the organic species. For water transport studies, this has often been achieved by examining the evaporation of water or condensation of water from a binary component particle that contains an involatile solute as the second component, for example sucrose or citric acid. 30,36 In measurements of the suppression in volatility of an organic component due to slow bulk diffusion, the focus has usually been on the rates of evaporation of a volatile under dry conditions where fluctuations in RH and water content (i.e. the presence of water as a plasticiser) can be avoided. 37-39 Alternatively, the volatilities of the organic components have been assumed to be so low that water transport can be studied in isolation, for example in the studies of water transport in SOA samples. 27 However, it should be expected that the mass transfer of water and SVOC would happen concurrently in any ambient aerosol. For example, both water and SVOC evaporation would be expected to be convoluted in the size and compositional response when the aerosol contains these two volatile components. Indeed, we have recently examined the evaporation kinetics of water (vapour pressure ~2 kPa) and maleic acid (vapour pressure ~10 -3 Pa) from ternary mixed component aerosol consisting of water, maleic acid and sucrose. 13 In this first attempt to examine the coupled loss of water, SVOCs, and the correlation with explicit measurements of particle viscosity, we showed that water loss remained fast while the SVOC loss could be suppressed by many orders of magnitude relative to its pure component vapour pressure.
An analytical treatment of the gas-particle partitioning of semi-volatile species has been described by Mai et al., 23 allowing an estimate of the timescale for evaporation of components of varying volatility.
It is appropriate at this point to us this model to assess the relative timescales of water and SVOC volatilisation from an aerosol particle. Such an assessment should allow us to identify the conditions under which the water loss rate is dominant, thereby justifying the use of the treatment described above for retrieving the diffusion constants of water. Conversely, this would also allow us to identify the timescales over which the water fraction equilibrates and the loss of the SVOC drives a measured size change. Thus, the Mai et al. 23 model was used to simulate the depletion kinetics of water or a hypothetical SVOC (with the same molecular mass as malonic acid) from a particle experiencing a small, but nonzero, change in external gas phase saturation of each species. Timescales were estimated for values of Dw and DSVOC varied logarithmically within atmospherically relevant bounds (10 -13 -10 -19 m 2 s -1 ). Timescales τSVOC and τwat are determined as the e-folding time of the concentration of each component remaining in the particle. The results indicate that for a 200nm radius particle, SVOC evaporation timescales will only become separable from those of water (τSVOC ≥ 10τwat) at around Dw = 10 -17 m 2 s -1 ( Figure 2a). Interestingly, τw continues to decrease as Dw increases above this value, whereas the SVOC volatilisation remains constant. This is due to the gas diffusional flux being substantially more limiting for an SVOC than water owing to the lower equilibrium vapour pressure of the SVOC (captured by the inverse dependence of τ on partial pressure, p, via the Henry's Law coefficient in equation (14) in the original publication 23 ). We also consider the equilibration timescales for particles of radius 5 μm in Figure 2b, a size of more relevance to particles that may be trapped during many single particle laboratory based studies, such as optical tweezers. The self-diffusion coefficient below which the timescales are found to be indistinguishable is 10 -16 m 2 s -1 , with the same divergence of regimes as observed in the 200nm case: the SVOC becomes gas transport limited while water retains a bulk kinetic limitation up to higher bulk diffusion constant.
For both particle sizes, the timescales for SVOC and water transport (Figure 2c) remain widely separated at all RHs and particle moisture contents. In producing Figure 2c we have assumed that Dw Other model parameters are as follows: temperature = 298 K, accommodation coefficient = 1, particle density = 1.4 gcm -3 , gas phase diffusion coefficient = 10 -5 m 2 s -1 . Water transport is consistently more rapid than SVOC transport as a consequence of the higher vapour pressure for water at all RHs and a bulk condensed phase diffusion constant for water that is considerably faster than expected based on the Stokes-Einstein equation. By contrast, the larger organic molecule is impeded more readily by increases in viscosity (Dw > DSVOC). 17 An additional consequence of this is that the water content adjusts rapidly to compositional change arising from volatilisation/condensation of an SVOC.

III. Experimental Description
We have provided comprehensive accounts of the time-dependent measurements of particle size and refractive index using aerosol optical tweezers (AOTs) in numerous previous publications 40, 41 and we will provide only a short overview of the instruments and procedures here. Optical tweezers are formed of a tightly focused 532nm (Laser Quantum, Opus 2W) laser beam, guided onto a vertically aligned oilimmersion objective and into a trapping cell. A mixed flow of ultrapure (dry) and humidified (wet) nitrogen is directed into the trapping cell to control the RH by altering the relative flow rates of both.
The RH is measured by a calibrated probe (HUMICAP HMT 330, Vaisala, error = ±2% RH) placed several centimetres before the trapping cell. An aerosol plume is generated using a medical nebuliser and drawn through the cell. Once a particle is captured in the optical trap, brightfield imaging may be used to visually assess the trapping stability. Inelastic backscattered light is diverted to a spectrograph (0.5 m focal length) coupled to a CCD (1 s time resolution, 0.05 nm/pixel).
Once the droplet equilibrates with the cell conditions and the radius is stable, a series of step changes of RH can be conducted. Throughout the measurement series, the droplet size and RI can be determined from the fingerprint of cavity enhanced Raman scattering at whispering gallery mode wavelengths. 42,43 The kinetic responses to changes in RH are examined either by considering the change in particle radius or the change in wavelength of a selected WGM as described in our earlier work. 31 The time intervals between RH steps are programmed in advance and were generally chosen to increase as the particles became more kinetically arrested (i.e. longer periods of RH stability at lower RH).
Here, the time-dependent wavelength shifts of individual WGMs in a Raman spectrum are used to infer the relaxation timescale in particle size following a step change in RH. Such data can be represented as a response function F( ( )), as in equation (5), where ( ) is the time dependent wavelength in nanometers. The response function is then fit to a modified stretched exponential function similar in form to the Kohlrausch-Williams-Watts equation: When compared with our previous fitting procedure, 31  The pre-exponential factor, s, then scales the data such that the fit and response functions both begin at unity. Here the e-folding timescale, , represents the kinetic limitation to the evaporation of water within the droplet. 41 For binary systems studied in previous publications, 31,44 the value of  increases as the initial and final RH are lowered.
Particle viscosity is determined using a dual trap AOTs, which has also been described in previous publications. 17,45 Two particles are independently captured in separate traps, which are formed using a holographic optical tweezer with a 532 nm wavelength laser (Laser Quantum, Opus 3W). By altering the kinoform applied to a spatial light modulator, a diffractive optical component that allows the generation of complex optical landscapes, 40 it is possible to manipulate the relative positions of the particles, and therefore to coalesce them. The relaxation of the aspect ratio of the resulting particle (from adjacent spheres to a single sphere) follows an exponential decay with a time constant, τ, that can be related to viscosity by 46 when the particles have a viscosity higher than the critical damping viscosity (~20 mPa s). σ is the surface tension and rfinal is the radius of the fully relaxed sphere.

IV. Results and Discussion
We first consider an approach for characterising and assessing the timescales of water and SVOC loss from viscous aerosol before presenting the analysis of water diffusion constants in organic aerosol of ternary composition. This is followed by a discussion of the effective suppression in volatility of the SVOC from viscous aerosol and the correlation with particle viscosity.

IV.a. Kinetics of Water Evaporation and Condensation in Ternary Aerosol with One Volatile Component
One hundred and seventy-two RH transitions have been conducted in the atmosphere surrounding twenty-seven ternary component aerosol particles consisting of water, sucrose and glucose. It should be noted that water is the only volatile component in these measurements; sucrose and glucose can be considered involatile. By capturing the non-exponential nature of relaxation in size (i.e. water partitioned to particle) following the RH step 47 using equation (6), the complex multi-exponential timedependence in droplet radius has been converted into a simple estimate of the timescale of internal water transport. 31 Time-constant values for all changes in RHs may be found in in Figure 3a. For the first time, we have also unconstrained the value of  (in Eq 6) during the fitting process (Figure 3b). The chemical significance of this parameter is somewhat less clear: In the glass physics community, where stretched exponential relaxation of macroscopic variables is commonplace, there is some consensus [48][49][50] that  is representative of the heterogeneity of the system, i.e. that there is an underlying distribution of mono-exponential events whose width becomes larger as  approaches zero. Previously, we have shown that fixed values of  between 0.4 and 0.6 represent the relaxation of binary sucrose water aerosol at low humidities; 31 similarly, a multi-exponential description of ternary systems is also appropriate ( In evaporation experiments, the time constant broadly increases as the final state of dehydration of the particle increases (i.e. as the RH change moves towards the bottom left of Figure 3(a)), unless the particle begins in a very low viscosity state at high RH. In addition, the slow release of water over long time is consistent with a value of  <1. Tentatively we can say that this is due to the larger variance in molecular evaporation timescales that together make up the observed radius response. However, we must be careful: stretched exponential behaviour is observed in every transition, including those that take place for viscous liquid aerosol above the moisture driven glass transition RH of the organic components. This suggests that we are using the stretched exponential to represent a fundamentally different physical process to that usually just ascribed to glasses.
The condensation experiments are somewhat more complex to characterise, involving the condensation of a pure water layer on the dehydrated particle surface. 44 Induction behaviour is observed in the first several hundred seconds of condensation for many of the measurements, where a thin water layer is likely to have formed on, but not yet been absorbed by, the dehydrated particle. This leads to values of  in excess of 1. 31 For purely gas-diffusion limited condensation, the size response should be independent of the internal composition with a value of  close to unity. 51 In order to successfully determine diffusion coefficients of water in the ternary mixture using the Fi_PaD model it is essential that an accurate radius response is inferred. The fully equilibrated droplet size and refractive index are determined by fitting the WGM fingerprint at a high initial RH (typically 70 -80 %) at which the droplet can be assumed to be in equilibrium with the gas phase and to be homogeneous in refractive index. Then, the droplet radius response following an RH change is inferred from the relative change in wavelength of a WGM; in such a kinetically limited system, it is inappropriate to rely on radius data produced by Mie theory fitting to a radially homogeneous sphere. 42 Instead, we must assume that the refractive index of the outer micron or so of the droplet is approximately constant throughout the transition and use the approach of Krieger and coworkers, who estimate radius changes from experimental wavelength shifts. 52 We consider, as an example, the diffusion constants of water in ternary particles containing glucose and sucrose at a 1:1 mass ratio, along with transitions in pure sucrose and pure glucose. Each transition (see Table 1) has been simulated with the Fi_PaD model to produce diffusion dependences, modelling the droplet as a two component system with the non-volatile component molecular mass taken as the weighted average of glucose and sucrose. As diffusion coefficients are assumed symmetrical (the same for all components at any given composition), no extra information may be gained by explicitly modelling the separate movements of the two organic components. The constants (Dw,org, C) within equation (1) were floated within the model as described earlier, until simulated radius responses agreed with experiment. The sigmoidal dependencies of Dw on mole fraction resulting from these individual optimal fits for binary aqueous/sucrose aerosol are shown in Figure 4 retrieved from 5 transitions in

RH.
A new method of error analysis has been employed to complement the Fickian modelling, allowing the uncertainty in Dw,org within the sensitivity of the tweezers measurement to be determined over many droplets and step changes in RH, as shown in Fig. 4. As total dehydration of the droplets is often not achieved, one must consider that significant deviations in the values of equation (1) returning the mean deviation, in log space, between two Vignes equations: the experimentally retrieved dependence Dw fit and an arbitrary function iteratively varying over ranges in physically meaningful values of Dw,org and C. RHinit and RHfin are the initial and final relative humidities over which the particle transitions, respectively. x w,i is the i th element of the mole fraction vector at each intermediate RH, over intervals of 1%. α(C) refers to equation (2). The division by nxw, the mole fraction range over which the particle transitions, converts m into a deviation per unit mole fraction, so that measurements with differing ranges RHfin -RHinit may be analysed in a consistent way. Figure S3 shows Here, we have used the merit function method described to estimate the diffusion constants Dw,org, with associated uncertainties in ternary sucrose/glucose/water aerosol prepared at a mass ratio of 1:1 for the two saccharides ( Figure 5 and Table 1). As might be expected, the fitted Dw,org coefficients and uncertainty envelopes for water the ternary mixture are found to be above the values for water in binary water/sucrose aerosol particles also determined in this work (Figures 5). In addition, the value for water in sucrose agrees well with recently reported values inferred from aerosol size response in an electrodynamic balance 6 (1×10 -12 cm 2 s -1 ) and more recent H/D isotope exchange measurements 30 (2.6×10 -17 cm 2 s -1 ). Additionally, increasing the relative molecular mass of the non-volatile component (i.e. comparing pure glucose to 1:1 glucose : sucrose to pure sucrose) inhibits water diffusion more greatly and the diffusion constant of water decreases.
A comprehensive comparison of our results applying this method to binary sucrose aerosol with previous optical tweezers and isotope exchange parametrisations is also shown in Figure 6. We have chosen to plot diffusion coefficients against mole fraction rather than gas phase activity, as previously reported, in recognition of the fact that this is the composition parameter implicitly used internally to evaluate D within Fi-PaD and within any Fickian based framework. Additionally, it is noted that activity does not appear anywhere in Eq. 1, and so a choice of thermodynamic model must then be used to implicitly relate particle phase mole fractions, xw, to water activity, in order to determine the range over which the model output is valid in each case. In addition, it is not clear whether current core-shell models are applicable in multicomponent mixtures with more than 3 components that are non-ideal. To clarify this, Fickian based diffusion models can account for non-ideal behaviour by combining Maxwell-Stefan diffusion coefficients with 'thermodynamic correction' factors through the following: 53 where [ ] is a matrix of Maxwell-Stefan diffusivities given by: and [Γ]is a matrix of thermodynamic factors: where is the Kronecker delta, Whilst Maxwell Stefan diffusion coefficients are symmetric, Fickian diffusion coefficients are not, unless the mixture is ideal or binary. By representing, and fitting, diffusion coefficients in terms of water activity, their use to then account for the compositional dependence of diffusion in non-ideal multicomponent mixtures is questionable. Expressing diffusion coefficients as a function of mole fractions allow the above expressions to be used, provided the matrix of activity coefficient derivatives can be calculated. In the above expressions, Fickian diffusion coefficients can be negative if the chemical potential gradient is such that solubility limitations drive diffusion. The efficacy of available numerical methods in modelling such multicomponent mixtures needs to be reviewed, and will form the basis of future work. For now, presenting diffusion coefficients as a function of mole fraction rather than water activity facilitates future investigations and is consistent with the assumption of ideality.
Following on from this, we stress that representing the diffusion coefficient as a function of mole fraction is distinct from presenting the Vignes dependence as a function of relative humidity, as previously reported. 30 We suggest that, based on the previous discussion, this latter form is less rigorous and necessarily requires an assumption that the thermodynamic relationships used for activity coefficient are entirely accurate. Therefore, the choice of activity coefficients can have a substantial influence on the appearance of Dw data if (either particle or gas phase) activity is chosen as the x axis.
For example, the frequently observed difference in curvature of Dw parametrisations at low activity between isotope exchange and optical tweezers data can be accounted for if it is recognised that α and xw(RH) models are different; in mole fraction space, the most recent sucrose isotope exchange data 30 appears similar in shape to that presented in Figure 6 although the irreconcilable difference in the Dw,org intercept remains. By contrast, the Fi-PaD model treats activity coefficients as unity, and so particle phase activity is equal to xw.

IV.b. Kinetics of Water Loss in ternary aerosol with 2 volatile components
As we previously reported for sucrose : maleic acid mixtures, 13 a complicating factor in assessing the time dependent size response when SVOCs are present is the exact deconvolution of the evaporation rates of competing volatiles from the droplet, water and the SVOC. In the limiting case in which the vapour pressure of the SVOC differs from that of water by a factor of 10 6 and when there is no kinetic impairment of water loss from the aerosol particle, rapid water evaporation will dominate the first several tends of seconds, after which an underlying linear time-dependent loss in r 2 becomes apparent. 54 Indeed, Figure 2 confirms the relative ranges in diffusion coefficients in which this separation in timescale occurs.
By simply fitting the stretched exponential function to the radius data, equation (6), normalised to the initial and final radii, the long tail can be assumed to have a constant gradient in r 2 with time once the water evaporation has ceased. This causes the r 2 loss to appear exactly in the residual between the data and the best fit, from which the gradient can easily be determined. In order to determine this gradient, it has been found that removing the data up to 30% of the way through an experiment provides a good compromise between fitting the largest number of points and ensuring the water loss is excluded (see Supplementary Information). If we assume that the SVOC loss begins immediately at the RH transition, and is masked by the water evaporation, a new function can be defined that takes into account the linear combination of the two contributions. A comparison of the two functions is shown in Figure 7, with both fit to typical optical tweezers data for ternary aerosol particles, with one or two volatile components. The rate at which the radius recedes, 2 , is kept constant at the value fit to the residual in the previous step, and A, r(∞), τ and β are allowed to vary. r(∞) is the final particle size containing the involatile solute and solvating water at equilibrium with the gas phase.
The deviation between this function and the experimental data is found to decrease by several orders of magnitude when compared to fits carried out using a standard stretched exponential on the same dataset (Equation 6).
The accuracy with which equation (13) reproduces radius responses for particles containing SVOCs (Fig 7b and section IVc) suggests it is valid to assume that transitions between differing semi-volatile loss rates occur rapidly upon relative humidity changes within the trapping cell. This is a consequence of the typical disparities in the time constants for SVOC and water loss illustrated in Figure 2 (typically more than a factor of 20), effectively allowing the two processes to be decoupled in the analysis.
Therefore, it is possible to probe the water loss response occurring simultaneously with the SVOC loss in the early stages of the particle response by simply subtracting the known r 2 profile from the experimental data. This produces a hypothetical 'water loss only' radius response. These data may then be modelled using Fi-PaD, making this the first instance of precise water diffusion coefficients being extracted from a ternary system with this distribution of volatilities. From a dataset of seventeen sucrose : malonic acid transitions, prepared at a 5:2 molar ratio, four were deemed to not be limited by the timescale of RH changes within the trapping cell (τ > 200s), and therefore ideal candidates for this analysis. The results are tabulated in Table 2. The assumed involatile molecular mass was once again the weighted average of the two organic components, with the composition assumed to be constant at its initial value throughout the experiment.
The compositional dependence of the diffusion constant ( Figure 8a) indicate a substantial plasticising effect of malonic acidjust under a 30% molar component within the organic matter enhances Dw,org by the same amount (relative to pure sucrose) as a one to one mass ratio of glucose (cf. figure 5). This is also borne out by presenting the diffusion dependences against viscosity data, collected via the particle coalescence methods described in Section II. Nine coalescence events were conducted between

IV.c. Kinetics of SVOC loss in Ternary
Finally, we consider the volatilisation kinetics of a malonic acid from an aqueous sucrose/malonic acid droplet, rather than water. Malonic acid is similar in size and chemical functionality to many organic molecules of atmospheric importance, as well as exhibiting a vapour pressure typical of an SVOC that may partition into the particle phase (10 -4 Pa). 55 The gradient determined for the SVOC loss can be converted into an effective vapour pressure using the equation 56 (we refer the reader to the SI for further details): This equation is derived from the Maxwell equation for a single component droplet, corrected by Raoult's law. Here R is the gas constant, T is the temperature, and MMA and Dg are the molecular mass and gas phase diffusion coefficient of malonic acid, taken as 104.1 Da and 7.2•10 -6 m 2 s -1 respectively.
ρdroplet is the density of the droplet, taken in this case to be the binary sucrose value, parametrised as a cubic function of mass fraction by Cai et. al. 57 , the malonic acid vapour pressure values determined from equation (14), are related to the pure component vapour pressure, p ○ MA, by xMA, the mole fraction in the droplet, and fMA, the activity coefficient at that mole fraction. These quantities, in addition to the mass fraction of malonic acid, FMA, were calculated using an activity/mole fraction relationship generated from the UManSysProp prediction tools 58 . Once pMA is corrected by these terms, we designate the calculated vapour pressure value as an 'effective pure component' value.
A gas flow removes evaporated SVOC from the trapping cell and maintains a background gas phase partial pressure that can be considered to be zero in the vicinity of the droplet. The gas velocity is sufficiently slow that no enhancement in evaporation rate from the droplet need be considered (i.e. the Sherwood number is ~2).Thus, in Figure 9 we report the RH dependence of the effective pure pressure of MA as the particle is dried. This is a broadly similar level of variance at low RHs to that reported in our previous publication for the volatility suppression of maleic acid in aqueous maleic acid/sucrose particles and reflects the increasing kinetic limitation imposed on volatilisation as the particle becomes more viscous. 13 One potential inaccuracy that may arise when examining inferred vapour pressures arises from the fact that particle composition is continuously changing as volatilisation proceeds. Indeed, the magnitude of the diffusive limitation to volatilisation will change as the relative proportion of remaining SVOC changes. This necessarily produces a large degree of variation of both the pure component vapour pressure estimate and retrieved refractive index at any given RH. Figure 9 shows an example of this: Each point is a unique single particle transition, with some droplets trapped and consecutively dried over multiple RH steps, while others were caught and rapidly dried just once. Data are coloured by the final retrieved refractive index to give an indication of droplet composition (see Supplementary   Information). Particles in which less malonic acid has previously volatilised exhibit refractive indices further from the pure sucrose value (nmelt = 1.5462 57 ) and higher effective vapour pressures, due to a decreased diffusive limitation. This trend leads to the spread in observed volatility for the five experiments with final RHs between 20 and 30%.

V. Conclusion
In the immediate aftermath of a rapid change in surrounding RH, a single optically tweezed organic aerosol particle will begin to vitrify, impeding the dynamics of its constituent molecules, both in the particle phase and at the droplet interface. Using Fickian diffusion modelling, we have determined the compositional dependence of the diffusion constants of water; from analysis of the particle radius Further studies suggest that the slow evaporation of malonic acid from ternary droplets is correlated with the conditioning history that the particle experiences, in addition to the composition of the starting solution and the ambient RH. Repeat experiments reveal substantial variation in this volatilisation.
Future publications will apply these methodologies to probing diffusion and volatility in increasingly complex systems, both in a 'bottom-up' sense, using multicomponent aerosol with well-defined starting compositions, and in a 'top-down' sense, through analysis of secondary organic matter produced via chamber oxidation of common VOCs. We will also begin to quantify the regional impacts of kinetic limitations to gas particle partitioning in dry air parcels, by embedding Fickian diffusion methods explicitly within a box model, to probe the evolving dynamics of an ensemble of particles.       Table 1. saccharide aerosol (as in Figure 4a). Each curve is produced from the mean Dw,org and C values of the transitions tabulated in Table 1.    Figure 4a, composition dependent diffusion coefficients fit to four radius responses using Fi-PaD (further details can be found in Table 2). (b) Single particle viscosity measurements, with a loglinear parametrisation. (c) Combination plot of data from (a) and (b) showing correlation between D and η. A comparison of the Stokes-Einstein predicted diffusivity for malonic acid (green) is also provided in panels (b) and (c).  (13) and (14) (closed circles). Points are coloured by the mean refractive index recorded after the initial water loss had ceased.