Equilibration times in viscous and viscoelastic aerosol particles †

Secondary organic aerosol (SOA) particles in Earth's atmosphere can exist in phase states where mass transport and chemical transformations are greatly impeded. Gas-particle partitioning in both semi-solid and solid ( e.g. glassy phase state) SOA particles can be kinetically limited, and these equilibration timescales are often evaluated using a Fickian model with a concentration-dependent di ﬀ usivity. Within that framework, particles are considered to be inviscid ﬂ uids and thus their rheology is ignored. In this work, rheological properties are not neglected and particle equilibration is investigated using both a viscous model (an incompressible Newtonian ﬂ uid) and a viscoelastic model (a Maxwell material). We derive analytic expressions for gas-particle equilibration times in viscous and viscoelastic aerosol particles and compare these times to the well-known Fickian equilibration time. These expressions are then applied to the important example of the equilibration of atmospheric aerosol particles with the surrounding relative humidity ( i.e. water uptake and loss). The two systems studied here are aqueous sucrose, which is often used as a surrogate for secondary organic material found in the atmosphere, and aqueous oxidized a -pinene, which is an atmospherically relevant secondary organic material. We show that viscous e ﬀ ects (i) are likely unimportant in single particle experiments, except under very dry conditions, and (ii) can result in water vapor equilibration times of minutes to hours at any altitude in the troposphere for accumulation mode particles.


Introduction
Binary diffusion in the liquid state is generally well-described by Fick's laws of diffusion.When assessing the timescale for establishing gas-particle equilibrium, a characteristic equilibration time, s D , based on the Fickian model is oen used in atmospheric aerosol science.For a particle with a binary (Fickian) diffusivity, D AB , in the condensed phase and radius, s, the equilibration timescale is This equation has been applied to a variety of mass transport calculations where gas-particle partitioning is limited by condensed phase diffusion, such as equilibration in aqueous droplets [1][2][3] and mixing in high viscosity secondary organic aerosol (SOA) particles. 4,51][12] These non-liquid states can have very high viscosities. 13,14For instance, glassy states that form in the SOA condensed phase can have viscosities on the order of 10 12 Pa s. 10 Therefore, one would anticipate that a Fickian framework, which treats the condensed phase as being an inviscid uid insofar as mass transport is concerned, would eventually become inaccurate.6][17][18][19][20][21] In these high viscosity systems, viscoelastic or purely viscous effects can readily dominate Fickian diffusion and very oen cannot be ignored. 22However, it would be difficult to directly apply this previous research to SOA particles as it has almost exclusively focused on polymer-penetrant systems where physical quantities are determined using models appropriate for polymer solutions. 19ig. 1 presents laboratory measurements of diffusivity and viscosity for two well-studied systems across a large range of water activities.4][35][36][37][38] Here, the Stokes-Einstein (SE) relation is used to convert diffusivity to viscosity or vice versa.Two trends are immediately apparent from Fig. 1.First, both viscosity and diffusivity vary by orders of magnitude across atmospherically relevant water activities.Second, the discrepancy between diffusivity (viscosity) predicted by the SE relation using viscosity (diffusivity) measurements can be quite large.At lower water activities, the predictions from the SE relation can be off by several orders magnitude.For instance, for aqueous sucrose at a water activity of 0.25, viscosity and diffusivity predicted using the SE relation differ by a factor of 8.2 Â 10 5 from their measured values.In addition to the breakdown of the SE relation, it is also clear from Fig. 1 that there will be a large range of water activities where the equilibrated condensed phase will not be a low-viscosity liquid.Typically, liquids have viscosities less than 10 2 Pa s and the phase state between 10 2 to 10 12 Pa s is characterized as a semi-solid. 12In both panels of Fig. 1, this semi-solid region covers a signicant portion of atmospherically relevant water activities.
In the current work, we study the inuence of uid rheology on the characteristic equilibration time of aerosol particles with their surroundings.We consider both viscous and viscoelastic uids as models for phase states relevant to atmospheric aerosol particles.These two models differ from the Fickian framework as viscosity is no longer zero during mass transport.The main question that we are concerned with is: How do the rheological properties of these phase states impact gas-particle partitioning?This is answered by deriving expressions for the characteristic equilibration time, analogous to eqn (1), for both viscous and viscoelastic spherical particles.4][35][36][37][38] Conversions from diffusivity to viscosity or vice versa were calculated using the SE relation with a temperature of 293 K and a molecular radius of 0.2 nm. 60Diffusivity measurements were performed on either binary (H 2 O + solute) or ternary (H 2 O + D 2 O + solute) systems.See referenced papers for details.
equilibration time for a range of atmospherically relevant particle sizes and conditions.We discuss the implications for single particle experiments and also examine timescales for representative conditions that exist in the planetary boundary layer, middle troposphere and upper troposphere.

Mass transport in a viscous and viscoelastic binary system
Mass transport of species a in a non-reacting solution is written using the convection-diffusion equation with no source term where r a is the mass concentration of species a, v is the mass average velocity, and j a is the mass ux of species a.
For a solution containing N components, we can express the mass density of the solution, r, with the partial specic volumes, v a , and mass fractions, w a , of species a using We have previously shown that, when v a are constant, the mass average velocity can be expressed in terms of mass uxes according to 39 This equation allows for the mass average velocity to be eliminated from eqn (2).For a binary system (N ¼ 2) with species A and B, eqn (4) along with the relation j A ¼ Àj B can be applied to eqn (2) to yield Therefore, to calculate the rate of change of the mass concentration in a binary system, what is now required is an expression for j A .For binary diffusion in the absence of temperature gradients and forced diffusion, the generalized Maxwell-Stefan (MS) equation 40 simplies to where, for species A, x A is the mole fraction, a A ¼ g A x A is the activity, g A is the activity coefficient, c A ¼ x A c is the molar concentration, f A ¼ r A v A is the volume fraction, c is the total molar concentration, R is the universal gas constant, T is the temperature, p is the local pressure, and D AB is the binary MS diffusivity.As indicated in eqn (6), the gradient of the natural logarithm of a A is evaluated at a constant T and p.In a binary system, the MS diffusivity and the Fickian diffusivity, D AB , are related through 41 where G AB is the thermodynamic factor for the binary system and is dened as Therefore, when g A is constant (e.g. for an ideal solution where it would be equal to one), we will have D 12 ¼ D 12 and the two diffusivities are identical.
The pressure gradient in eqn ( 6) is determined by using the equation of motion where s 4 is the stress tensor, and g is the external force per unit mass.For simplicity, we ignore the effects of inertia by setting the le-hand side of eqn (9) to zero.As we will only consider cases where g ¼ 0, eqn (9) simplies to A viscoelastic stress-strain relationship will be used to describe the stress in eqn (10).The standard Maxwell model 42 is chosen here.In this linear model of viscoelasticity, the uid consists of a spring and a dashpot in series 19 v s where E is the Young's modulus, b ¼ E/h is the ratio of the elasticity to viscosity, h is the concentration-dependent viscosity of the dashpot in the Maxwell model, 3 is the innitesimal strain tensor, and is the rate-of-strain tensor.A relaxation time can be dened as We will restrict our analysis to spherical coordinates with radial symmetry.Eqn (10) will be and the non-zero components of the rate-of-strain tensor are The solution to eqn (11) when h and E depend on the concentration of A for a system that is initially unstressed ( s s qq ðr; tÞ ¼ s ff ðr; tÞ ¼ In these equations, the concentration history of the system, extending back to the initial unstressed state at t ¼ 0, inuences the stresses at the current time, t.The origin of this memory effect in a viscoelastic system is the presence of a nite relaxation time (s l > 0).
In addition to the viscoelastic model, we will also consider a viscous model.When E / N, the Maxwell model (eqn ( 11)) reduces to which is identical to the description of stress that is found using an incompressible Newtonian uid (when the term with dilatational viscosity is set to zero, see p. 241 from ref. 40).The nonzero components of the stress tensor will simply be In the viscous model, only the current state of the system determines the stress.The memory effect seen in the solution to the viscoelastic model is absent as relaxation is now instantaneous (s l ¼ 0).
We now proceed in our analysis by considering the case where only small changes in species A and B occur during the mass transport process (this is sometimes referred to as differential sorption). 17We will linearize equations by assuming that the coefficients that appear in them are constant.When eqn ( 13) is inserted into eqn (6), the mass ux of species A can be written as where Eqn ( 5) and ( 19) along with the appropriate stress expressions can be used to determine governing partial differential equations for both viscous and viscoelastic mass transport during differential sorption.For the viscous model, the governing equation can be found by inserting eqn ( 18) into (19) and then using eqn (4) to eliminate v r from the resulting equation.All of the uxes, j r,A , and their partial derivatives can then be replaced through the application of eqn (5) and its partial derivatives.This yields The governing equation for the viscoelastic model is determined as follows.First, we recognize that for constant h and E, eqn (15) These two equations are then inserted into eqn (19) and then eqn ( 4) is used to eliminate v r .The result will be an equation that contains integrals over the interval of 0 to t.If the partial derivative with respect to time is applied to both sides of this equation, the resulting equation will contain no integrals.Then, just like the viscous case, application of eqn ( 5) and its partial derivatives can be used to replace all uxes with partial derivatives of concentrations.When this is done, the viscoelastic governing equation is found to be It is clear from inspection that when the relaxation time goes to zero (s l / 0), the viscoelastic governing equation (eqn ( 23)) reduces to the viscous governing equation (eqn ( 20)).

Characteristic equilibration times
The viscous and viscoelastic equilibration times for the binary system where species A is the gas-phase species/penetrant and species B is the condensed/particle phase species can be found by solving the linear boundary problem with the appropriate boundary conditions.At the centre of the particle we have the usual symmetry condition For nanometer-sized particles with low condensed phase diffusivity, it is usually an excellent approximation to assume that the concentration of A at the surface can be xed at its nal value. 5,43,44Therefore, in this particle-diffusion limited case, the boundary condition at the particle surface will be where r A,N is the bulk gas-phase mass concentration of species A in the condensed-phase.We remark that the approximation that yields eqn (25) has only been validated within a Fickian framework. 5,43,44he particle radius is assumed to be constant here so that we do not have a moving boundary problem.The initial condition for the problem will be where r A,0 is the concentration of species A in the particle at t ¼ 0. For the viscous case, we can solve the boundary-value problem using separation of variables to get where the eigenvalues are The characteristic equilibration time, s N , for the viscous incompressible Newtonian uid is then For the case when v A ¼ v B there is no viscous ow and s N will reduce to the Fickian equilibration time (eqn (1)), provided that the solution is ideal (g A ¼ 1).However, if v A s v B then s N will not reduce to s D , even if the viscosity is zero.This occurs because the typical Fickian treatment neglects convective mass transport due to density changes.Also, we see that for s 2 [ p 2 w, the viscous pressure gradients will not signicantly affect the equilibration time.Fig. 2 shows concentration proles calculated using eqn (27) during both uptake and loss in the case when Fickian diffusion dominates (a and c) and viscous effects dominate (b and d).In panels b and d, the concentration of A is always uniform (except near r ¼ s) and equilibration is determined by material response during swelling or shrinkage (rheological limitations).
An approximate expression for s N that is more straightforward to apply in many situations of interest, and also useful in terms of physical insight, can be obtained through the following assumptions: the solution is ideal, w A ( w B , and the s 2 /p 2 4 term can be replaced by s D .Then, eqn (29) becomes In Fig. 3, eqn ( 29) and ( 30) are compared across a range of water activities and radii for the system of aqueous sucrose.Even when the condition w A ( w B is no longer true, the accuracy of eqn (30) does not begin to deteriorate.This is because s D will dominate in that case, so the second term on right-hand side of eqn ( 30) is negligible.Aside from being slightly simpler than eqn (29), one other advantage of eqn ( 30) is that M B is not needed if w A is already known.So, for instance, in a calculation involving water transport in SOA material, the molar mass of the SOA, which is typically dened as an average, will not be required.
For the viscoelastic case, we require a second initial condition to formulate the boundary-value problem due to the second order time derivative in eqn (23).We will use With this additional condition, the solution to the viscoelastic boundary-value problem is where The parameters w, 4, and l n have the same denitions as in the viscous case.The difference between the viscous and viscoelastic solutions is that the time-dependent part of eqn (32)  in the viscoelastic solution is that of a damped harmonic oscillator whereas the viscous solution was that of a homogeneous rst-order differential equation.Therefore, d n is the damping parameter and u 0,n is the natural frequency.As should be expected, in the limit where s l / 0, eqn (32) reduces to the time-dependent viscous concentration prole (eqn ( 27)).
For large t in the overdamped case (d n 2 > u 0,n 2 ), eqn (32) can be approximated as where L n ¼ d n À a n .Fig. 4 shows the parameters d 1 2 and u 0,1 2 across an activity range of zero to one for the model system of aqueous sucrose at several different radii.In all cases, the overdamped condition is satised.We will therefore assume that the overdamped solution is the only one that is relevant for atmospheric systems and that the characteristic time, s V , for the viscoelastic uid can subsequently be written as The viscoelastic equilibration time given here does not exactly reduce to the viscous equilibration time (eqn (29)) when the relaxation time goes to zero due to the approximations used to arrive at eqn (35).

Calculated characteristic equilibration times
Fig. 5 shows the characteristic equilibration times for aqueous sucrose particles calculated across a range of water activities and radii using the expressions for viscous and viscoelastic equilibration times found in Section 2.2.Although the equilibration times were evaluated from 10 mm down to 1 nm, it should be anticipated that as this lower size limit is approached, the continuum assumption used throughout this work will fail.The Fickian equilibration time (eqn (1)) is also Fig. 3 Comparison between the viscous characteristic equilibration time calculated using the exact (eqn ( 29)) and approximate (eqn ( 30)) equations as a function of radius at several different activities, a A (or mass fractions w A ). Species A is water and species B is sucrose.Parameters for aqueous sucrose at T ¼ 293 K were used for the calculations.presented for comparison.The parameterizations for diffusivity and activity for aqueous sucrose were taken from ref. 23 and viscosity was calculated using AIOMFAC-VISC. 45,46For the Young's modulus, a value of 33 GPa was used in all calculations, which is representative of reported measurements for sucrose crystals at room temperature. 47,48Note that regardless of particle size and water activity, there will always be a small difference between the Fickian equilibration time and the other two equilibration times because v A s v B and g A s 1 in the viscous and viscoelastic calculations.The origin of this discrepancy was discussed in Section 2.2.
In Fig. 5 we see that at lower RHs (drier conditions), the radius at which noticeable separation between the three curves rst appears becomes larger.The most signicant difference across all panels is always between the Fickian curve and the viscous and viscoelastic curves.In the driest example shown (panel a), the Fickian times are several orders of magnitude faster than the viscous or viscoelastic times across sizes that would be relevant to atmospheric SOA (e.g.size distributions that are peaked below 100 nm (ref.49-51)).In contrast, both the viscous and viscoelastic curves are qualitatively similar.Even quantitatively, the difference between the two curves is small.As the Young's modulus is always positive, the viscoelastic equilibration time will always be greater than the viscous equilibration time.Overall, as the RH increases, the region where this difference occurs moves to lower radii.At the highest RH (panel d), the three curves are nearly indistinguishable from each other over the plotted range of radii.
The oen large difference between the Fickian curve and the other two curves in Fig. 5 can be understood by comparing the expression for the Fickian equilibration time (eqn (1)) to the expression for the viscous equilibration time (eqn ( 29)).When it is no longer true that s 2 [ p 2 w, then viscous pressure gradients play a signicant role during equilibration.For a non-zero viscosity, this can always be achieved by either (i) decreasing the particle radius or (ii) increasing w.In the approximate equation for the viscous equilibration time (eqn ( 30)), the viscous equilibration time is separated into a sum of the Fickian equilibration time and a product involving the viscosity.It is apparent that simply increasing the viscosity does not guarantee that viscous effects will dominate, as the SE relation states that the diffusivity will simultaneously decrease (resulting in the Fickian equilibration time increasing).If the SE relation is perfectly satised, increasing the viscosity will not affect the difference between the viscous and Fickian times for radii greater than 1 nm in the case of aqueous sucrose.However, as was discussed in Section 1, the SE relation is oen inaccurate for lower water activities.In Fig. 1 it is clear that as water activity decreases (and thus viscosity increases), deviations from the SE relation become more signicant.The measured viscosity becomes orders of magnitude larger than that predicted from the diffusivity using the SE relation.Therefore, the breakdown of the SE relation provides an explanation as to why the radius at which deviations begin to occur is larger for lower water activities.
In Fig. 6 and 7 the calculations from Fig. 5 are repeated but using parameters for aqueous oxidized a-pinene particles.Viscosity was calculated with AIOMFAC-VISC using a system of 15 SOA surrogate components.Tabulated AIOMFAC-VISC calculations are provided in the ESI.† As there are no measurements of the Young's modulus for oxidized a-pinene, only viscous and Fickian calculations are shown here.To our knowledge, there are only two activity-dependent diffusivity data sets for oxidized a-pinene (with the data set from ref. 7 also being temperature-dependent).Both sets of measurements along with their ts are shown in Fig. 1b, where it can be Fig. 5 The Fickian (eqn (1)), viscous (eqn ( 29)), and viscoelastic (eqn ( 35)) characteristic equilibration time as a function of radius for an aqueous sucrose particle at a RH of (a) 10%, (b) 20%, (c) 30%, and (d) 50%.Parameters for aqueous sucrose at T ¼ 293 K were used for the calculations.
seen that, at lower water activities, the discrepancy between the two data sets is around one to two orders of magnitude.This type of difference is not seen with any of the ve aqueous sucrose data sets shown in Fig. 1a.The most straightforward explanation as to why there are such large differences between the oxidized a-pinene measurements is that the oxidation products in the two cases are signicantly different.The chemical composition of the oxidized a-pinene products is sensitive to the conditions under which they are produced, 38 and the two samples were generated under very different conditions; e.g. the samples in ref. 7 were produced using oxidation with OH radicals, while the samples in ref. 33 were Fig. 6 The Fickian (eqn (1)) and viscous (eqn ( 29)) characteristic equilibration time as a function of radius for an aqueous oxidized a-pinene particle at a RH of (a) 10%, (b) 20%, (c) 30%, and (d) 50%.Parameters for aqueous oxidized a-pinene at T ¼ 293 K were used for the calculations, with the diffusivity being calculated using the parameterization from ref. 7. Fig. 7 The Fickian (eqn (1)) and viscous (eqn (29)) characteristic equilibration time as a function of radius for an aqueous oxidized a-pinene particle at a RH of (a) 10%, (b) 20%, (c) 30%, and (d) 50%.Parameters for aqueous oxidized a-pinene at T ¼ 280 K were used for the calculations, with the diffusivity being calculated using a polynomial fit to the measurements from ref. 33.produced using O 3 oxidation.This cannot be neglected, as it can be seen that the choice of the diffusivity for oxidized apinene has a signicant impact on the calculated equilibration times: at a radius of 100 nm the Fickian and viscous times are nearly identical in every panel in Fig. 6, whereas in Fig. 7, large differences are seen in both panels (a) and (b).The implications of this will be discussed in Section 3.3.

Implications for single particle measurements
For more than a decade, single particle measurements have been used to determine the condensed phase diffusivity in aerosol particles.Six of the eight diffusivity data sets presented in Fig. 1 were measured by holding a particle in an electrodynamic balance or with optical tweezers (the other two data sets used a liquid disk on a hydrophobic slide).In all cases, the dimensions of the particles were on the order of microns.Analysis in these previous experiments was always performed within the Fickian framework.From Section 3.1, we see that such an approach should be satisfactory as, for particles in this size regime, viscous and viscoelastic effects will not affect water transport and can be safely neglected for all but the driest conditions.

Atmospheric implications
The calculations in Fig. 5-7 demonstrate that treating SOA particles, especially those with sizes in the sub-100 nm size regime (e.g.glass formers), as inviscid uids can potentially result in a many orders of magnitude error in the characteristic equilibration time.The implications of this result to Earth's atmosphere are explored by calculating characteristic equilibration times for particles with the surrounding RH using simulated RHs and temperatures from the global chemistry climate model EMAC 52 at the Earth's surface, 850 hPa ($1.35 km), and 500 hPa ($5.5 km).The RH and temperature data that was used is plotted in the ESI.† Fig. 8 shows the Fickian and viscous equilibration times for aqueous sucrose and Fig. 9 shows the Fickian and viscous equilibration times for aqueous oxidized a-pinene.Viscoelastic equilibration times were not calculated as there is no temperature-dependent measurements of the Young's modulus available for either substance.However, based on calculations in Fig. 5, viscoelastic times should not differ signicantly from viscous times.A particle radius of 100 nm was chosen as being representative of accumulation mode particles.Prior to our discussion, we note that sharp decreases in viscosity have been predicted for radii below 50 nm. 53While this effect will not impact the calculations shown in this section, a viscosity decrease would certainly need to be considered for smaller radii.
For the surrogate system of aqueous sucrose particles, shown in Fig. 8, the Fickian model never predicts equilibration times of more than a few minutes at the surface or 850 hPa.This is consistent with earlier calculations using the same Fickian Fig. 8 Characteristic equilibration time for an aqueous sucrose particle with a radius of 100 nm with the surrounding RH at (a and b) the Earth's surface, (c and d) 850 hPa, and (e and f) 500 hPa.Fickian equilibration times (a, c and e) were calculated using eqn (1) and viscous equilibration times (b, d and f) were calculated using eqn (29).Parameters for aqueous sucrose were used for the calculations.RHs and temperatures that were used as input for these calculations were taken from the global chemistry climate model EMAC, 52 and are identical to those used in ref. 8. model and similar parameters, 8 where longer times (tens of minutes or hours) were only found to occur at 500 hPa.In contrast to the Fickian calculations, the viscous model predicts equilibration times of several hours for certain locations at both the surface and 850 hPa.The implications of these longer equilibration timescales for water uptake and loss are the following: Predictions concerning heterogeneous chemistry will be affected as the particle's phase state may not be that which is found by assuming equilibrium with the surrounding RH. 4,54,55 This also applies to calculations of equilibrium particle size and the partitioning of semivolatile organics into the condensed phase. 56,57The phase state of the particle is also key to ice nucleation pathways.In typical adiabatic updras of 0.1 m s À1 , glassy or solid-containing aerosol particles with equilibration times that exceed several minutes (longer than typical CCN activation time periods), may initiate ice nucleation heterogeneously under sufficiently cold conditions. 6,9nlike aqueous sucrose particles, the equilibration times from both the Fickian and viscous calculations for aqueous oxidized a-pinene particles, shown in Fig. 9, are indistinguishable from each other at the surface and 850 hPa and are always less than a few minutes.Here, timescales longer than several minutes only occurs at 500 hPa.The similarity between Fickian and viscous times for aqueous oxidized a-pinene was anticipated based on Fig. 6, where the two times were nearly identical at a radius of 100 nm for all RHs that were shown.Of course, as the products of oxidized a-pinene can vary signicantly depending on the conditions under which the oxidation takes place, it is not possible to state that the results in Fig. 9 will apply to all oxidized a-pinene particles.In Section 3.1 we discussed how oxidized a-pinene generated under different conditions can have viscous times that are much longer than Fickian times for radii of 100 nm; i.e. the results shown in Fig. 7.However, the diffusivity function used in that example does not have a temperature dependence (it is only valid for 280 K) and cannot be used in the atmospheric calculations performed here.
It appears that the answer to the question of whether or not viscous effects can signicantly effect the equlibration of SOA particles with surrounding water vapor is dependent on the organic material that forms particles.In Fig. 8, the effect of viscosity was signicant and equilibration times much longer than Fickian ones were observed at both the surface and 850 hPa.In contrast, in Fig. 9, the difference between the viscous and Fickian times was small everywhere except for a few locations at 500 hPa.Unfortunately, there is currently a dearth of temperature-dependent diffusivity measurements of oxidized aqueous organic material.More laboratory-based studies will be required to determine whether the situation in Fig. 8 or 9 is dominant in the atmosphere.
Finally, although the calculations here were restricted to binary systems consisting of water and an organic solute, inorganic species can also be a signicant component of Fig. 9 Characteristic equilibration time for an aqueous oxidized a-pinene particle with a radius of 100 nm with the surrounding RH at (a and b) the Earth's surface, (c and d) 850 hPa, and (e and f) 500 hPa.Fickian equilibration times (a, c and e) were calculated using eqn (1) and viscous equilibration times (b, d and f) were calculated using eqn (29).Parameters for aqueous oxidized a-pinene were used for the calculations, with the diffusivity being calculated using the parameterization from ref. 7. RHs and temperatures that were used as input for these calculations were taken from the global chemistry climate model EMAC, 52 and are identical to those used in ref. 8.  organic aerosol. 58One model ternary system of atmospheric relevance where equilibration times have been measured is water + citric acid + ammonium sulfate. 59In that case, the Fickian equilibration time initially increased as ammonium sulfate was added to the mixture of water and citric acid.For RH between 10 and 40%, it was only aer a 1 : 1 molar ratio of ammonium sulfate to citric acid was exceeded that the Fickian equilibration time became faster than that of the binary case of water and citric acid.For RH between 10 and 40%, AIOMFAC-VISC calculations for this ternary system predict only a decrease in viscosity with the addition of ammonium sulfate.Therefore, for this system, eqn (30) predicts that the viscous effects become less important as an inorganic solute is added to the organic and water mixture.

Summary
The main results from this work were the formulation and solution to the problem of mass transport in both a viscous and viscoelastic spherical particle.The derived analytic expressions for gas-particle equilibration times in viscous (eqn ( 29)) and viscoelastic (eqn (35)) particles can potentially have a broad range of applications beyond those discussed here (e.g.differential sorption in polymer spheres).When applied to systems of atmospheric relevance, it was shown that viscous and viscoelastic effects will be unimportant in laboratory-based studies of water transport involving super-micron particles.In contrast, in the sub-micron regime, and in particular below 100 nm, the effects become signicant and both the viscous and viscoelastic equilibration times can be several orders of magnitude slower than the Fickian time for the same environmental conditions.For the model system of aqueous sucrose, this result was demonstrated to have signicant implications for the equilibration with water vapor in the atmosphere.However, for aqueous oxidized a-pinene, the viscous equilibration time showed little difference from the Fickian equilibration time.

Fig. 2
Fig. 2 Concentration profiles for species A during uptake (a and b) and loss (c and d) calculated using eqn (27).In panels (a) and (c) s 2 [ p 2 w and transport is dominated by Fickian diffusion.In panels (b) and (d) s 2 ( p 2 w and transport is dominated by non-Fickian effects (which are viscous in this case).

Fig. 4
Fig. 4 Viscoelastic parameters d 1 2 and u 0,1 2 as a function of water activity.Eqn (33) and parameters for aqueous sucrose at T ¼ 293 K were used for the calculations.