Amorphous phase state diagrams and viscosity of ternary aqueous organic/organic and inorganic/organic mixtures.

A Dimer Coagulation, Isolation and Coalescence (DCIC) technique is used to probe the phase behaviour and glass transition temperatures of ternary aerosol mixtures. The DCIC technique is used to perform temperature and relative humidity dependent viscosity measurements at viscosities near 5 (cid:2) 10 6 Pa s. Measurements include organic–organic and organic–inorganic mixtures composed of sucrose–citric acid and sucrose–sodium nitrate. The data reported here add additional insight into the wide discrepancies in glass transition temperatures reported for pure sodium nitrate. The phase diagram model used in the work of Rothfuss and Petters ( Phys. Chem. Chem. Phys. , 2017, 19 , 6532–6545) is expanded to include multiple solute components. Data and model predictions of the mixtures are in good agreement with the modified model. These measurements are compared with values from Holographic Optical Tweezer (HOT) measurements taken at room temperature. Overall, the viscosities determined from the DCIC and HOT techniques are in good agreement.


I. Introduction
A significant fraction of ambient submicron non-refractory aerosol mass can be attributed to organic components. 1 The organic fraction evolves through functionalisation, fragmentation and oligomerisation reactions, 2 can become internally mixed with inorganic aerosol components and responds to changes in relative humidity (RH) by absorbing/desorbing water. The oxidation of volatile organic compounds in the gas phase leads to sufficiently high molecular weight/low volatility compounds to form secondary organic aerosol (SOA). Recent studies suggests that SOA can exist in a highly viscous or even an amorphous solid phase, with significant implications for aerosol light interactions in the atmosphere, cloud droplet and ice crystal formation, and further heterogeneous chemical processing. [3][4][5][6][7][8] A number of studies report that SOA can have a viscosity between 10 4 -10 9 Pa s at ambient temperatures, dependent on SOA formation conditions, i.e. the identity of the precursor volatile organic compound, and both the concentration of the VOC and RH. [7][8][9][10] This range in viscosity suggests that semi-solid viscosities may be relatively common, particularly at the lower temperatures found in the middle and upper troposphere. 11,12 Aerosol viscosity influences diffusion constants of molecules within a particle, 8,13,14 rates of heterogeneous reaction, 15 the photodegradation of molecules within an aerosol particle 16 and the apparent volatility of the aerosol components. 15 In addition, highly viscous particles may shield reactive species in the particle bulk, allowing for enhanced long-range transport of environmental contaminants. 17,18 Under dry conditions, aerosol composition can achieve supersaturated states with respect to solute concentration, existing in a metastable non-crystalline amorphous or subcooled liquid form. 19 Consequently, aerosol-based techniques are a necessity for measuring the viscosity of the amorphous, metastable states formed in supersaturated aerosol. This is challenging for a number of reasons including the small sizes of aerosol particles (1 nm-10 mm), the wide range in accessible viscosity to a metastable phase (10 À3 to 410 12 Pa s) 13,20 and the often extremely small sample volumes derived from sampling techniques (mg to mg). As a result, the majority of the approaches available for measuring aerosol viscosity are indirect, such as the poke-and-flow technique 21 and fluorescence lifetime imaging in which viscosity is inferred from the fluorescence lifetime of a molecular rotor. [22][23][24][25] Despite these experimental challenges, it is now apparent that the addition of functional groups and increasing molecular weight are generally consistent with increasing viscosity. 4,26,27 Numerous studies have now linked increasing viscosity with progressive oxidation using OH radicals or ozonolysis. 22,[28][29][30][31] In order to rationalise the behaviour of more chemically complex atmospheric organic particles, laboratory studies have used highly oxygenated species, such as sucrose and citric acid, as proxies for viscous SOA. 32 For these systems, a decrease in relative humidity (RH) has been shown to suppress the rate of semi-volatile organic loss, reaction rates or uptake of gaseous species, indicative of a concomitant increase in viscosity. 15,33,34 Aerosol coalescence techniques have been developed to infer particle viscosity, where the time-dependence of the shape of the coagulated dimer can be related to the viscosity. 35 Direct measurements of particle (3-6 mm) viscosity can be made using Holographic Optical Tweezers (HOT) via the coalescence of two optically confined particles. 27,35 Two different analysis approaches are employed dependent on the viscosity regime. At low viscosity, the elastically back scattered light may be used to determine particle viscosity with coalescence occurring on the 10's ms timescale. For coalescence events with particles of viscosity higher than B1 Pa s, brightfield imaging reveals the relaxation in shape directly, with the timescale for relaxation used to determine viscosity. 35 Pajunoja et al. 36 use a variant of the dimer coalescence technique, where naturally coagulated SOA aerosol are impacted onto a thin layer of carbon, and the degree of coalescence is analysed using scanning electron microscopy. 36 Further, Zhang et al. 37 determine the viscosity of secondary organic material that coagulated during particle formation and growth, by measuring the shape factor change via measuring particle electrical mobility and particle mass. 37 Rothfuss and Petters 38 prepare and isolated dimers B100 nm in diameter using a dimer coagulation, isolation, and coalescence (DCIC) technique where dimers are produced by coagulation of oppositely charged particles followed by removal of uncoagulated particles via electrostatic filtration. 38,39 The DCIC technique can be used to measure the temperature dependence of viscosity and the temperature dependence of the RH range where viscosity transitions from 1 Â 10 6 to 3 Â 10 7 Pa s. The primary advantage of this technique is the ability to infer the temperature dependence of the viscosity and potentially extrapolate the measurements to estimate glass transition temperatures. Data for sucrose aerosol determined with the submicron size DCIC technique were consistent with super-micron size HOT coalescence technique by Power et al. 20,35,40 Using the DCIC approach, a phase diagram model for binary aqueous sucrose has been proposed by Rothfuss and Petters 38 using a modified Vogel-Fulcher-Tammann equation, the Gordon-Taylor mixing rule, and a mass based hygroscopicity model. Here, we extend this model to multicomponent systems and test the model with two ternary model systems, aqueous mixtures of sucrose-citric acid and sucrose-sodium nitrate. The viscosity of these mixed component aerosol systems can be expected to transition across the sensitive range of the DCIC measurement approach making them appropriate model systems. Indeed, the measurements presented here are the first study in which the viscosity of ternary mixtures are measured and compared using both the DCIC and HOT coalescence techniques at different temperatures and relative humidities. The DCIC technique is used to probe the phase behaviour and glass transition temperatures of the ternary aerosol mixtures. The coagulated dimer measurements are compared with a limited number of mixture data points from HOT measurements, which provide an absolute measure of aerosol viscosity and provide additional viscosity data at 22 1C. Such measurements on mixtures of increasingly complex composition are crucial to determine if the properties of mixtures can be treated purely as additive, based on data from binary component mixtures, or if a more complex interplay of component interactions is required to predict the properties of mixtures. 20,38,39 Further, this study provides more data on the influence of an added inorganic component on mixture viscosity, for comparison with previous work of Power et al. who examined the aqueous sucrose-NaCl ternary system as well as the appropriate binary aerosol. 20

II.a DCIC experimental method
This experimental setup and dimer coalescence technique has been extensively discussed in previous publications. However, the experimental set up utilised in this work differs slightly from the previous literature and the modified instrument design is shown in Fig. 1. 26,38,39 Aerosol was generated using a dual syringe pump system feeding two constant output atomisers (TSI 3076), generating two separate aerosol streams, which are then passed through two sets of tandem Nafion dryers (PD-50T-24MSS, PermaPure). The two streams of dried aerosol are then passed through 210 Po radiation sources to yield equilibrium distributions of charged particles before injection into two separate differential mobility analysers (DMA), a TSI-DMA (TSI 3071) and a high flow DMA referred to as the HFDMA. 41,42 As depicted in Fig. 1, the TSI-DMA is connected to a power supply applying a negative bias potential and, hence, selects positively charged particles. The HFDMA is connected to a power supply applying a positive bias potential and selectively filters for particles of negative polarity. Both the TSI-DMA and HFDMA were configured to produce a stream of 80 nm aerosol and operated with sheath-to-sample flow ratios of 3 L min À1 : 1.5 L min À1 .
After the DMAs, additional drying was implemented via the use of one or more cold traps positioned within each size selected stream. In all experiments, the two dried, size selected streams of aerosol were routed through two separate 12 00 long, 1/2 00 outer diameter stainless steel tubes positioned on a cold plate (ThermoElectric Cooling America AHP-1200CV) set to a temperature of -30 1C. When further drying was necessary, a second set of identical tubes were also positioned inline; these tubes were inserted between blocks of dry ice (À79 1C). Next, the two dry aerosol streams enter a commercial chest freezer which housed most of the remaining instrument elements. The freezer temperature could be controlled between À25 1C and +30 1C using a chilled bath circulator (Thermo Scientific C40P) interfaced with a heat exchanger placed inside the freezer. The temperature inside the freezer is referred to as the system temperature, and this value is tabulated alongside the measurement results (ESI †). The aerosol streams were combined inside the freezer (for the first time) in a coagulator tube (approximate residence time of 5-6 s with 3 L min À1 flow rate). An optional third port introduces particle-free pre-humidified air into the coagulator. This flow was generated by routing particle-free zero air (Teledyne model 701) through a Nafion membrane submerged in a temperature controlled water bath. A combination of flow rate (0-1 L min À1 ) and temperature adjustment of the Nafion water bath (0-25 1C) was used to control the dew-point temperature inside the coagulator.
After exiting the coagulator, a flow of 0.3 L min À1 was diverted to monitor the total particle concentration measured using a condensation nuclei counter (TSI 3020) operated in bulk scattering mode. This measurement was used to monitor stability of the particle generation system. A second flow, also 0.3 L min À1 , was routed through an electrostatic filter to remove all charged particles. The remaining neutral aerosol, comprising dimers from coagulation of oppositely charged particles and particles that spontaneously discharge en-route between the DMA and the coagulator, then pass through a temperaturecontrolled chamber. Temperature control was achieved using a thermal conditioner (temperature control unit) which is a Peltier-type thermoelectric heat exchanger (CP-200TT, TE Technology) and the temperature reached by the temperature control unit was measured using a thermistor (TE Technology MP-3193). Residence time in the loop was 5 s. A chilled-mirror dew point hygrometer (General Eastern Optisonde) was used to determine the water vapour concentration in the gas-phase.
The charge neutral particles exiting the electrostatic filter were brought to charge equilibrium using another 210 Po radiation source. The size distribution was measured using a radial DMA 43 (RDMA, sheath-to-sample flow ratio 2 : 0.3 L min À1 ) operated in scanning mobility particle sizer mode. 44 The particle count from RDMA monodisperse output was measured using a second condensation nuclei counter (TSI 3020) operated in pulse counting mode. Data reduction from the RDMA output is as described in Rothfuss and Petters, 38 and involves peak fitting to find the mode diameter of the dimer size distribution.
In a few positions throughout the instrument shown in Fig. 1, the RH was measured using temperature and humidity sensors (Rotronic HygroClip2). However, the RH used for data analysis was the water vapour concentration measured by the Dew Point Hygrometer and the RH was determined by the ratio of the vapour pressure and saturation vapour pressure evaluated at the temperature inside the conditioning loop using the parameterizations of Murphy and Koop, (eqn (8) and (10)). 45 Experiments can be separated into two distinct types: isothermal humidification and heating cycles performed under dry conditions. These experiment types are detailed in Rothfuss and Petters. 39 Briefly, isothermal humidification experiments fix the temperature of the temperature control unit inside the freezer, while slowly raising the water vapour content through the humidification line until the dimers relax to a spherical shape on the timescale of the measurement. In contrast, heating cycle experiments are performed under dry conditions (corresponding to the lowest dew-point that can be achieved with drying and the baseline system temperature inside the freezer), raising the temperature of the temperature control unit inside the freezer until the dimers relax to a spherical shape on the measurement timescale.
Some differences between the present experiments and the referenced work are summarized below. In isothermal humidification experiments, the humidified flow and the two size-selected particle flows were merged simultaneously via a Swagelok cross piece positioned at the inlet of the coagulation chamber; by contrast, the humidified flow was first introduced into the outflow of the TSI-DMA in the cited work. Additionally, under some circumstances the humidifier itself was bypassed, directing up to 1.2 L min À1 of zero air directly into the system. This was sometimes necessary at colder temperatures because the minimum conditioner RH value facilitated by the Nafion dryers and cold traps was not low enough to fully resolve the fully uncoagulated dimer mobility diameter. In heating cycle experiments the temperature control unit was sufficient.
The measured SMPS mode diameters from an isothermal humidification scan or a temperature scan are converted to the particle geometry factor x using the expression given by Rothfuss and Petters (2016), where D p is the SMPS peak diameter, D uc the fully uncoalesced and D c is the fully coalesced dimer diameters. An example of the resulting experimental data as a function of RH is reported for an isothermal humidification in Fig. 2. The solid black line represents a fit to eqn (2). In eqn (2) k is a steepness parameter, RH is the measured RH from the dew point hygrometer and RH r is the relative humidity representative of the midpoint of coalescence relaxation. For the heating cycle measurement, the terms RH and RH r are replaced with T and T r , respectively where T r is representative of the temperature at which the midpoint of coalescence relaxation occurs, and T is the temperature measured by the temperature control unit. Mapping from dimer mobility diameter to viscosity is based on the numerical procedure described Rothfuss and Petters. 39 The relationship takes into account the monomer diameter, the residence time in the temperature control unit where relaxation takes place and the surface tension of the particle. For 80 nm monomers and conditioner residence time of 5 s, measurements of RH r and T r correspond to a viscosity of B5 Â 10 6 Pa s, which is slightly modified by potential changes in surface tension with composition. Both the DCIC and HOT coalescence technique require an estimate of surface tension to compute viscosity. A number of methods have been derived for combining surface tensions of simple binary solutions to estimate the surface tensions of complex mixtures at variable RH. 46,47 In this manuscript we consider ternary mixture surface tensions which have been determined using E-AIM Model III 47 and parametrisation for surface tension values used have been provided in the ESI † (S.2).

II.b HOT experimental method
The DCIC measurements are compared with a limited number of viscosity measurements using holographic optical tweezers (HOT). The HOT technique has been thoroughly described in previous publications. 27,35 Thus, it is only briefly discussed here as the implementation was consistent with our earlier work. These measurements provide the dependence of viscosity on RH at 22 1C (ambient temperature). A gradient force optical trap is formed by tightly focusing a laser (532 nm, Laser Quantum Opus 3 W) through a microscope objective with a high numerical aperture. This trapping beam is split using a Hamamatsu spatial light modulator (SLM, X10468) to form two optical traps. The RH of the trapping region is controlled using a ratio of dry and wet nitrogen gas and measured using a precell capacitance probe that is regularly calibrated (Honeywell, HIH-4202A). Optically confined droplets (3-7 mm) equilibrate with the gas phase RH of the chamber. Droplet equilibration (e.g., stable droplet size) is monitored via stimulated Raman scattering appearing at wavelengths commensurate with whispering gallery modes and is used to retrieve droplet size and refractive index with an accuracy of oAE2 nm and oAE0.0005 respectively. 48 The two optically confined droplets are brought close enough together to coalesce. For droplets with a viscosity 41 Pa s, the timescale of coalescence is determined using direct brightfield imaging with a high frame rate camera (Vision Research, Phantom v. 7.3) with relaxation times 410 ms. Alternatively, for droplets with viscosity up to 100 mPa s and relaxation timescales of o1 ms, elastic light scattering is used to determine the shape relaxation timescale which is then used to calculate the RH dependent viscosity. 35,40 III. Phase diagram model In this section, a model is introduced that can be used to represent isopleths of constant viscosity (eqn (3)-(7)) and provide a way of representing the T and RH dependence of the viscosity of the ternary mixtures. All symbols are defined in Table 1. The model is similar to that of Rothfuss and Petters 38 who based this on the work of Berkemeier et al. 49 Key differences to the previous version of the model include: first, the extension of the Gordon-Taylor equation and water activity model to multiple solute components; second, the treatment of water content dependent hygroscopicity parameters, k m , and the definition of the Gordon-Taylor constants as the inverse k i = 1/k (GT,i) of that used previously. The inverse definition was introduced for simplicity. w s = (1 + k m (a w /(1 À a w ))) À1 (3) T = (B/log 10 Z À A) + T 0 Eqn (3) represents a mass based hygroscopicity parameter, k m , determined using the total mass fraction of solute in the droplet, w s and the water activity, a w . 50,51 The mass based hygroscopicity parameter for a mixture containing multiple solute components is determined using the k m,i of the individual components and the mass fraction of the ith component in the dry particle, e i , as shown in eqn (4). Note that a w = RH/(100%) and the effect of particle surface curvature on the relative humidity at which the particle composition is equilibrated (the Kelvin effect) is not considered. ZSR mixing is assumed to model the water content of multicomponent mixtures. Eqn (5) introduces the Gordon-Taylor mixing rules for multiple dry solutes, which is re-arranged from eqn (4) in Koop et al. 4 Eqn (6) is used to determine parameter T 0 that introduces a temperature dependence to the viscosity with slope and curvature parameters (A, B) given by the VFT equation, log 10 Z = A + B/(T À T 0 ), and a dry glass transition temperature T g defined by eqn (4). This is the same assumption made by Rothfuss and Petters. 39 Eqn (7) is used to determine the temperature at which the particle has a viscosity, Z. A and B represent fitted slope and curvature parameters describing the temperature dependence of viscosity determined by Rothfuss and Petters. 39 Values uses to initialize the model are summarized in Table 2. These values are taken from the literature or from this work. Values derived from this work are explained later in the text. Fig. 3(a) shows a heating cycle experiment for a dry 60 : 40 sucrose-sodium nitrate mixture. Fig. 3(b) shows how the shape factor data are used to derive particle viscosity inferred as a function of temperature. The calculations are based on Rothfuss and Petters 38 and modified Frenkel theory 60 and relate x to Z via a lookup table. Calculations are performed for x values computed from the logistic fit for 1 o x o 4. The calculations assume a residence time of 5 s in the conditioner, a value for D mono = 80 nm, and a surface tension corresponding to an estimate for the aqueous mixture (Section S2 of the ESI †). Also shown in Fig. 3(b) is the impact of the uncertainty in surface tension on the determination of particle viscosity. The green error bars indicate the maximum uncertainty in particle viscosity assuming a value for the surface tension in the range of 0.03 o s o 0.15 J m À2 , circumscribing a generous range of possible values for surface tension of different compounds, ranging from low surface energy organic materials 61 to high surface energy inorganic melts. 62 The range demonstrates that the uncertainty on derived temperature for the 5 Â 10 6 Pa s viscosity transition Table 1 Notation for multicomponent phase diagram model, tabulated definition of all symbols used in eqn (3)-(7)

A, B
Fitted slope and curvature parameters describing the temperature dependence of viscosity T Temperature where the particle has viscosity Z a w Water activity of solution T g Glass transition temperature of the mixed particle (water + dry solutes) e i Mass fraction of the ith component in the dry particle T g,i Glass transition temperature of ith solute; i = w denotes water k m,i (a w ) Mass weighted hygroscopicity parameter of the ith solute, with parameterized dependence on water content T 0 Asymptotic temperature that produces a viscosity temperature dependence that has a slope and curvature parameters (A, B) given by the VFT equation and a dry glass transition temperature T g k i Binary Gordon-Taylor constant for the mixture of the ith component with water due to the surface tension assumption is small. The range in modeled actual surface tensions for the different mixtures is slightly smaller and ranges between 0.06 J m À2 (pure citric acid) and 0.130 J m À2 (pure sodium nitrate melt). Modelled surface tension values are used for the analysis. The resulting uncertainty in the derived temperature transition due to the uncertainty of the model surface tension estimate of the mixture (ESI, † S.2) is much smaller than the indication in figure in 3(b). Further, the glass transition temperature is estimated by extrapolating the VFT fit to a viscosity corresponding to 10 12 Pa s 38 and shown in Fig. 3(b) by the purple point. The uncertainty of the extrapolation is estimated to be AE10 1C based on results discussed later (Fig. 5).
To date extrapolation from DCIC has only been applied to two cases where glass transition temperature is known, pure citric acid and pure sucrose. For both substances the error is within the AE10 1C estimate. Fig. 4 shows a typical dataset obtained from the HOT instrument; using coalescence measurements across a range of relative humidities, a parametrisation of the viscosity of the mixture can be generated. Measurements for mixed ternary aqueous sucrose-citric acid aerosol particles are compared with the previously published parameterisations of the simple binary systems. 27 The parametrisations are generated using a polynomial fitting through the data points of viscosity collected using the HOT; the uncertainty associated with the parametrisation is a result of the uncertainty in fitting to the experimental data. In Fig. 4 each measurement (data point) for the 60 : 40 sucrose : citric acid mixture represents an experimentally determined viscosity from a single coalescence event. Data points are not shown for the binary aqueous-sucrose and aqueous-citric mixtures for overall figure clarity and because they are presented in the work of Song et al. 27 For comparison with DCIC measurements, HOT measurements are all taken at  22 AE 2 1C and the corresponding RH r value at which the DCIC measurement is made is estimated assuming a viscosity of 5 Â 10 6 Pa s. Advantageously, this technique provides measurements over a much wider range in viscosity than the DCIC technique. As a disadvantage, there is no temperature control, and this means that the HOT can provide an additional measurement only at room temperature. Coincidentally, the DCIC isothermal humidification experiments at room temperature are challenging due to the large amount of moisture that needs to be added. Thus the HOT data extends the temperature range for representing the viscous properties of the aerosol and complements the DCIC isothermal humidification approach. Fig. 5(a) and (b) summarise the data taken in this study for heating cycles under constant dry conditions for the aqueous sucrose-citric acid mixtures. In Fig. 5(a) the values of temperature correspond to the temperature at which the viscosity of the mixture is 5 Â 10 6 Pa s. The data show a gradual increase in the relaxation temperature with increasing sucrose mass fraction. The mixing rule is not linear in mass fraction space, though the convexity is relatively small.
The values in Fig. 5(b) correspond to the glass transition temperatures estimated from the modified VFT fit described in Fig. 3 from each cycle. Available literature values from Rothfuss and Petters 26 for the binary aqueous-solute system (i.e. aqueoussucrose) are also shown here, with good agreement with the DCIC measurements and VFT model extrapolation. The data in Fig. 5(b) were used to construct a binary mixing model for the glass transition temperature of sucrose-citric acid mixtures, which is shown by the solid line. The model was obtained by a nonlinear least square fit to the data of the full model (eqn (5)-(6)). The fit was performed as follows. First glass transition temperatures for pure sucrose, pure citric acid and k i for sucrose/water ( Table 2) are fixed. Next the RH in the model is set to zero. The remaining k i for citric acid/water is varied such that the error between the model and the data in Fig. 5(b) is minimized. Finally, the VFT parameters for pure sucrose are used to estimate the mixing rule for the relaxation temperature in Fig. 5(a). Therefore, the model lines in Fig. 5(a) and (b) are effectively a one parameter fit to a model whose remaining parameters were constrained from pure component data identified in the literature. The resulting model provides an acceptable representation of dry mixtures of sucrose and citric acid. Note that the scatter in the data relative to the model Fig. 5(b) is approximately twice as large as the scatter in Fig. 5(a), which is due to the added uncertainty introduced by the extrapolation of the data to T g using the VFT fit. The excellent agreement of the model and data in Fig. 5(a) provide an encouraging first  validation of the model assumptions, since the data are fit to the extrapolated data in Fig. 5(b) and since pure component VFT data were applied to the mixture. Fig. 6(a) and (b) are the corresponding plots to Fig. 5 but for the sucrose-NaNO 3 mixtures. The fitting procedure to construct model lines is identical to that used for the sucrose/citric acid data shown in Fig. 5. Note that there is significant disagreement in the literature for the glass transition temperature for NaNO 3 . Measurements by Dette and Koop 53 (open triangles in Fig. 6b) use differential scanning calorimetry (DSC) to determine the glass transition temperature of several organic/inorganic mixtures including sucrose and NaNO 3 . In their technique, aerosols are impacted on a surface and crystallisation is observed for several of the organic/inorganic mixtures considered. They report glass transition at 293 K and 287 K, for pure NaNO 3 , equivalent to placing a limit of 10 12 Pa s on the viscosity of NaNO 3 under fully dry conditions at room temperature. 53 However, this value of the glass transition temperature is not directly measured but is extrapolated from the values for the dry sucrose/NaNO 3 mixtures from a NaNO 3 dry mass fraction of 0.8.
In contrast, previous reports (open square in Fig. 6b) suggest that the glass transition temperature of sodium nitrate is in the range 216-225 K, [54][55][56][57][58] considerably lower than the value suggested by Dette and Koop. 53 Therefore, two models were constructed, one with T g = 292 K and one with T g = 220 K. Either model can explain our observations. DCIC measurements for sucrose mass fractions o0.5 were attempted. However, these particles were fully coalesced at room temperature. Suspecting a possible measurement artefact, these experiments were repeated several months later with a slightly reconfigured system. Repeated points at sucrose mass fraction 0.6 and 0.7 are in excellent agreement with the initial dataset. Again, no coalescence was observable for sucrose mass fraction o0.5 consistent with full coalescence on the measurement timescale. The DCIC technique requires the dimers to be sufficiently viscous at the cold/dry limit of the measurement scan to retain their hourglass shape. The fact that we could not produce DCIC data for sucrose mass fractions o0.5 is consistent with the cooler T g model because otherwise we would expect viscosity to be higher and measurements to be feasible at higher NaNO 3 fractions.
The DCIC measurements for sucrose mass fractions o0.5 were attempted for RH as low as B2%. Baldelli et al. 63 report a viscosity of B1 Pa s at a mass fraction water of B0.1 for aqueous NaNO 3 , a composition approaching a pure liquid NaNO 3 melt and corresponding to an RH of 24%. Sucrose has a much higher viscosity in solution (410 6 Pa s) at the same water mass fraction, 27,35 which is what allowed us to measure viscosity for sucrose mass fractions 40.5. Note the 12-order of magnitude difference reported for NaNO 3 at room temperature, between 1 Pa s at an RH of 24% and 10 12 Pa s (glass transition) at an RH of 0%. At RH B2%, DCIC measurements were not possible. Although this may be attributable to retained humidity in the coalescence chamber, these observations are also consistent with a viscosity of NaNO 3 that remains lower than the 5 Â 10 6 Pa s threshold at all T and RH studied. [54][55][56][57][58]63 The discrepancy in estimated glass transition temperature of NaNO 3 between near room temperature and B220 K has also been extensively discussed by Dette and Koop. 53 They note a contradiction where in some cases atomized NaNO 3 particles do not effloresce [64][65][66][67][68] and form amorphous particles with near zero water content at 0.05% RH, and in other cases efflorescence is observed. [66][67][68] Dette and Koop 53 argue that the inconsistency in NaNO 3 efflorescence data could be consistent with a T g near room temperature, with the ambient temperature of 290 K well below the melting point of crystalline NaNO 3 . However, a DSC measurement does not directly imply a viscosity. Given that the DSC and viscosity observations operate on different physical principles, the apparent discrepancy could indeed be explained with all cited observations being correct.   7 shows the relationship between T and RH for a value of viscosity 5 Â 10 6 Pa s for mixtures of citric acid and sucrose and for six different mixture compositions. Lines represent isopleths of constant viscosity (5 Â 10 6 Pa s) and the relationship between temperature and RH, for the sucrose-citric acid, as determined from the modified VFT model. Data from both the DCIC temperature and RH scans are presented alongside additional measurements using the HOT technique at 22 AE 2 1C. For the DCIC measurements, the exact temperature of the 5 Â 10 6 Pa s transition is recorded during the experiment. Error bars show the confidence interval of the fit for the DCIC measurements using eqn (1). For the HOT measurements, the RH-dependent viscosity is directly measured and used to generate a RH-dependent viscosity parametrisation (examples in Fig. 4). The error in the HOT measurement is determined from the range in RH values at a viscosity of 5 Â 10 6 Pa s, corresponding to that of the DCIC measurement. The large error for the HOT measurement in Fig. 7(c) is attributable to wide RH range at a comparable viscosity of 5 Â 10 6 as shown in Fig. 4. Data for pure sucrose are taken from Rothfuss and Petters 39 and are shown for reference. Points at RH = 0% are identical to those in Fig. 5(a).
The data in Fig. 7 show that the increase of the citric acid mixing fraction shifts the entire viscosity isopleth to colder temperatures and lower relative humidities. At dry conditions, the effect is due to the binary mixing of two substances with different T g 's, as shown in Fig. 5 and 6. At elevated RH, the water content of the internally mixed particle is also altered by the changed hygroscopic response of the mixture. The net effect is predicted by the modified VFT model. We note that the modified VFT model was only constrained with data from Fig. 5 and material properties listed in Table 2 that have been previously reported. These values include water activity dependent surface tension and hygroscopicity parametrisations that are available in the ESI. † Thus, the agreement between the model and data is an independent test of the adequacy of the model framework for predicting the behaviour of ternary mixtures.
In Fig. 7, DCIC and HOT measurements for the sucrosecitric acid mixture are in reasonably good agreement with the model prediction. We believe that discrepancies between the data and the model are mostly caused by measurement errors that may be larger than those indicated by the error bars. For example, the RH measurement in the present version of the DCIC system is more uncertain when dewpoint temperatures that approach the limit of the dewpoint hygrometer (À15 to À20 1C). Similarly, the inference of the precise 5 Â 10 6 Pa s from HOT measurements requires interpolation from a sometimes relatively noisy datasets as shown in Fig. 4 and the supplement of Rothfuss and Petters. 39 However, model uncertainties also exist. Specifically, the location of the viscosity isopleth is sensitive to the hygroscopicity of the aerosol. 39 The modelled hygroscopicity relies on accurate representation of the water activity vs. water content relationship for binary mixtures of water and solute below room temperature and in highly concentrated solutions. Such data are not widely available. In addition, the hygroscopicity of the mixed particle is represented by the ZSR mixing rule, 69 which assumes that no solute-solute interactions exist. This assumption is likely questionable for these highly concentrated ternary mixtures. Finally, the conjecture that the VFT parameters A and B (Table 2), are only weakly dependent of solute composition over a wide range of compositions has not been rigorously evaluated. Failure of this assumption could lead to deviations between the model and the data. Considering these possible uncertainties, we conclude that the agreement between data and model is satisfactory. More ternary systems should be evaluated before the modified VFT model is more widely applied in atmospheric models. Fig. 8 shows the relationship between T and RH for a mixture with viscosity of 5 Â 10 6 Pa s for mixtures of sucrose and NaNO 3 at four different mixing fractions. The plot is constructed in the same manner as Fig. 7, with one difference. The modified VFT model is applied for T g,NaNO 3 = 220 K (solid line) and T g,NaNO 3 = 292 K, (dashed lines), subsuming the two dry mixing models shown in Fig. 5. As would be expected, the difference between the two models increases with increasing NaNO 3 mass fraction in the particle. The T g,NaNO 3 = 220 K model is in better agreement with the data for the 70 : 30 and 60 : 40 mixture. We note that the model viscosity isopleth is also nonlinear for these mixtures. This non-linearity is caused by a strong dependence of on water activity for NaNO 3 at low RH (ESI, † S.1). We note that no water activity vs. composition data is available for RH o 10%, and k m is assumed to be equal to the value at RH 10% for this region of the phase diagram. The other caveats noted in the discussion of Fig. 7 also apply to the comparison between model observations for the sucrose-NaNO 3 mixtures. Given these uncertainties, we consider the agreement of model and data encouraging, adding to the confidence that a properly constrained modified VFT model has potential to predict the viscosity isopleth for ternary mixtures over a wide range of RH and temperature.

V. Conclusions
Measurements are presented for sucrose-citric acid and sucrosesodium nitrate mixtures using the DCIC technique with complementary HOT measurements, which extend the accessible temperature range to 22 1C. The modified VFT fit, Gordon-Taylor equation, and the mass based hygroscopicity parameter have been successfully applied to ternary mixtures. There is excellent agreement between experimental DCIC measurements and the modified-VFT model for the heating cycles measurements of sucrose-citric acid mixtures ( Fig. 5 and 6) and satisfactory agreement for the isothermal humidification experiments ( Fig. 7 and 8).
The viscosity and glass transition temperature for pure component sodium nitrate has been extensively discussed with available literature values of T g = 216-225 K 54-58 and T g = 292 K extrapolated from DSC measurements reproduced from Dette and Koop. 53 In summary, the viscosity of pure sodium nitrate is unlikely to reach a viscosity of 5 Â 10 6 Pa s, the viscosity required for DCIC measurements. This conclusion is further corroborated by the better agreement between the modified-VFT fit to sucrose-sodium nitrate mixtures, using a T g = 220 K, and the DCIC measurements presented in this manuscript. Further efforts are needed to link viscosity measurements, which measure the particles' resistance to flow, and glass transition temperature measurements, which probe to response of the material to heat flow. Finally, the successful application of the multicomponent phase diagram model to ternary solute component systems is encouraging and leaves room for further expansion to more compositionally complex solutions/ samples, including secondary organic aerosols.

Conflicts of interest
There are no conflicts to declare. Fig. 8 Relationship between T, RH and Z = 5 Â 10 6 Pa s for mixtures of sodium nitrate and sucrose with solutions listed combined with respect to mass percent. In (a) 90 : 10 sucrose-NaNO 3 , (b) 80 : 20 sucrose-NaNO 3 , (c) 70 : 30 sucrose-NaNO 3 and in (d) 60 : 40 sucrose-NaNO 3 . Solid line corresponds to a sodium nitrate T g of 220 K and the dashed line corresponds to a sodium nitrate T g of 292 K. Red points are data for 100% sucrose taken from Rothfuss and Petters 39 and Power et al. 40 and are repeated in all panels for reference.