High Homogeneous Freezing Onsets of Sulfuric Acid Aerosol at Cirrus Temperatures

. Homogeneous freezing of aqueous solution aerosol particles is an important process for cloud ice formation in the upper troposphere. There the air temperature is low, the ice supersaturation can be high, and the concentration of ice-nucleating particles is too low to initiate and dominate cirrus cloud formation by heterogeneous ice nucleation processes. The most common description to quantify homogeneous freezing processes is based on the water activity criterion (WAC) as proposed by Koop et al. (2000). The WAC describes the homogeneous nucleation rate coefﬁcients only as a function of the 5 water activity, which makes this approach well applicable in numerical models. In this study, we investigate the homogeneous freezing behaviour of aqueous sulfuric acid aerosol particles by means of a comprehensive collection of laboratory-based homogeneous freezing experiments conducted at the AIDA (Aerosol Interaction and Dynamics in the Atmosphere) cloud simulation chamber, which were conducted as part of 17 measurement campaigns since 2007. The most recent experiments were conducted during October 2020 with special emphasis on temperatures below 200 K . Aqueous sulfuric acid aerosol 10 particles of high purity were generated by particle nucleation in a gas ﬂow composed of clean synthetic air and sulfuric acid vapor, which was added to the AIDA The resulting chamber aerosol had number concentrations from 30 cm − 3 up to several thousand cm − 3 with particle diameters ranging from about 30 nm to 1.1 µ m . Homogeneous freezing of the aerosol particles was measured at simulated cirrus formation conditions in a wide range of temperatures between 185 K and 230 K with a steady increase of relative humidity during each experiment. At temperatures between about 205 K and about 230 K , the AIDA results agree well with the WAC-based predictions of homogeneous freezing onsets. At lower temperatures, however, the AIDA results show an increasing deviation from the WAC-based predictions towards higher freezing onsets. For temperatures between 185 K and 205 K , the WAC-based ice saturation ratios for homogeneous freezing onsets increase from about 1.6 to 1.7, whereas the AIDA measurements show an increase from about 1.7 to 2.0 in the same temperature range. Based on the experimental results of our direct measurements, we suggest a new ﬁt line to formulate the onset conditions of homogeneous freezing of sulfuric acid aerosol particles as an isoline for nucleation rate coefﬁcients between 5 · 10 8 cm − 3 s − 1 and 10 13 cm − 3 s − 1 . The potential signiﬁcant impacts of the higher homogeneous freezing thresholds, as directly observed in the AIDA


Introduction
For predicting the weather and climate on our planet, it is of particular importance to understand and describe the cloud processes in the atmosphere. The evolution of mixed-phase and cirrus clouds is strongly influenced by the presence of cloud ice (Lynch, 2002). Cirrus clouds generally form in the cold upper troposphere and occur in particular in tropical regions but also in mid and high latitudes. By being able to scatter and absorb light, cirrus clouds are an important factor in the global radiation 30 budget. Solid aerosol particles often act as promoter of cloud ice formation in the atmosphere via heterogeneous nucleation at low ice supersaturations. In regions with low concentrations of ice-nucleating particles, the homogeneous nucleation of aqueous solution aerosol particles becomes more relevant for cirrus clouds at conditions of lower temperatures and higher ice supersaturations. (Cziczo and Froyd, 2014;Cziczo et al., 2013;Lynch, 2002). Koop et al. (2000) suggested the water activity criterion (WAC) to predict isolines for the homogeneous freezing rate coeffi-35 cients of aqueous solutions. This approach (hereafter called Koop2000) was deduced from measurements of the homogeneous freezing temperature observed for different aqueous solutions in dependence on their composition. According to their observations, the ice nucleation behavior only depends on the water activity of the solute, which in thermodynamic equilibrium equals the relative humidity of the surrounding environment (Koop, 2015). Thus, the parameterization of the homogeneous nucleation rate coefficients and related freezing onsets does not depend on the type of the solute, which makes this approach very suitable 40 for atmospheric models.
As aqueous sulfuric acid solution (H 2 SO 4 /H 2 O) aerosol particles are frequently present in the upper troposphere and lower stratosphere, these particles are of particular importance for freezing processes at these altitudes. Therefore, the experiments in this study focus on this aerosol type. Homogeneous freezing onsets of sulfuric acid available from literature are compiled in Table 1 and plotted in Fig. 1 in the water activity-temperature-space. In cases where water activities or relative humidities of J V = 10 13 cm −3 s −1 (∆a w = 0.32, dashed line) and J V = 5 · 10 8 cm −3 s −1 (∆a w = 0.3, dotted line)  and corrected according to the revised version of the homogeneous freezing lines by Koop and Zobrist (2009). 55 The earliest freezing experiments with aqueous sulfuric acid solutions at variable solute concentrations and temperatures were performed with macroscopic samples (Ohtake, 1993;Beyer et al., 1994;Song, 1994). The observed onsets deviate significantly from Koop2000 with low water activities close to the melting point curve. More recent flow tube experiments, where Fourier transform infrared (FTIR) spectroscopy was used to derive the concentration of freezing sulfuric acid aerosol (Bertram et al., 1996;Clapp et al., 1997), reported lower water activities especially for the freezing onsets at low temperatures (T < 60 220 K). However, the majority of the compiled freezing onsets in Fig. 1 agree with the Koop2000 predictions considering their measurement uncertainties. Homogeneous freezing experiments performed in the AIDA (Aerosol Interaction and Dynamics in the Atmosphere) cloud simulation chamber Mangold et al., 2005) show a tendency of freezing onsets towards higher water activities compared to Koop2000 at low temperatures (T < 200 K). For T < 200 K, Koop2000 predicts homogeneous freezing thresholds ranging between ice saturation ratios of 1.6 and 1.7. However, in an aircraft campaign in 65 2004, enhanced ice saturation ratios up to 2.3 at about 187 K were observed in the upper troposphere (Jensen et al., 2005), clearly exceeding the Koop2000 freezing thresholds. Also in a later aircraft campaign in 2006, more than half of the relative humidity measurements showed values exceeding Koop2000 at upper tropospheric temperatures (Lawson et al., 2008). Other atmospheric observations at T < 200 K reported atmospheric relative humidities predominantly below Koop2000, but in a few cases the threshold was also exceeded (Krämer et al., 2009;Krämer et al., 2020). These atmospheric observations seemingly 70 contradict Koop2000 if the assumption holds that ice saturation ratios are unlikely to exceed the homogeneous freezing thresholds in the atmosphere. To explain this discrepancy, the formation of metastable ice phases with higher saturation pressure than stable forms of ice was suggested by various authors (Malkin et al., 2012;Russo et al., 2014). Also, the presence of organic films on the aerosol particles has been discussed, which would slow down the water vapor uptake and thus the aerosol equilibration to the surrounding conditions (Jensen et al., 2005). It was also suggested that parameterizations for the saturation 75 pressure over supercooled liquid water at low temperatures are at least 20% too low (Jensen et al., 2005). Other authors argued that ice nucleation in solutions is accompanied and influenced by the expulsion of ions from the ice lattice (Bogdan and Molina, 2010;Shafique et al., 2016), which form a highly concentrated residual solution. It was suggested that the formation of such mixed-phase cloud particles could produce higher ice supersaturations than expected according to Koop2000 (Bogdan and Molina, 2010).

80
There are only few in situ measurements of cirrus and their microphysical properties in the field, e.g. in aircraft campaigns (Krämer et al., 2020). These measurements are often limited, as the collection of cloud samples in the atmosphere usually takes place after the onset of ice formation, so that the specific conditions at freezing onset are rarely matched in aircraft studies. However, an improved understanding of atmospheric homogeneous freezing processes can be obtained from laboratory experiments under well controlled conditions. Prompted by the observed tendency of homogeneous freezing onsets towards 85 higher ice saturation ratios than predicted by Koop2000 in previous AIDA studies at low temperatures Mangold et al., 2005), an extended set of homogeneous freezing AIDA experiments focused on the freezing of H 2 SO 4 /H 2 O aerosol particles conducted since 2007 was re-analysed in this study. To quantify the observed trend and the potential deviation from Koop2000 especially at low temperatures, the re-analysed experiments have been complemented by a more recent series of AIDA measurements focusing on temperatures below about 200 K. With this extensive data set of H 2 SO 4 /H 2 O freezing 90 onsets we can provide valuable addition to existing data sets and contribute to a better understanding of homogeneous freezing processes in the atmosphere.

Experimental setup
All experiments have been performed according to experimental procedures as described by Möhler et al. (2003), Benz et al. 95 (2005) and Wagner et al. (2008). Therefore, we only provide a brief summary of the experimental setup and procedures here.
The AIDA facility basically consists of a cloud chamber with 84 m 3 volume, which is situated in a thermal housing and can be cooled to temperatures as low as 183 K. Before starting an experiment, the AIDA cloud chamber is carefully cleaned by evacuating it to a pressure below about 0.1 hPa followed by flushing with dry and particle-free synthetic air several times between 10 hPa and 1 hPa. Then, pure water vapor from a heated liquid water reservoir is added to renew the ice layer on the 100 inner chamber walls and the AIDA is filled with dry synthetic air to ambient pressure. In the cleaned and re-filled chamber, the background aerosol particle concentration is usually lower than 0.1 cm −3 .
Fresh sulfuric acid aerosol particles were generated by passing synthetic air over a small amount of concentrated sulfuric acid located in a heated glass tube. Passive cooling of the mixture in the inlet tube was sufficient to induce the homogeneous nucleation of pure H 2 SO 4 /H 2 O particles. The number concentration and size of the particles was controlled by the setting of 105 the tube temperature and the air flow (Wagner et al., 2008). For experiments since 2017, an additional dilution flow of synthetic air was added downstream of the heated glass tube in order to better control the number concentration of generated particles.
For example, a reservoir temperature of about 383 K, an air flow of about 0.75 standard l min −1 and a dilution flow of about 14 standard l min −1 resulted in a median diameter of about 350 nm for the aerosol particles added to the cloud chamber.
A homogeneous distribution of aerosol particles throughout the cloud chamber is achieved by a mixing fan located one 110 meter above the bottom of the chamber. A few minutes after the end of injection, the aerosol particle number concentration and size distribution were measured with a Condensation Particle Counter (CPC 3010, TSI GmbH), an Aerodynamic Particle Sizer (APS 3321, TSI Inc.), and a Scanning Mobility Particle Sizer (SMPS 3934, TSI Inc.), see Sect. 2.2 for details. A cloud expansion experiment is then started by opening the valves to a vacuum pump operated at controlled pump speed. Immediately after the pump starts, the pressure and thus the temperature inside the chamber start to drop. At the beginning of the experiment, 115 the gas temperature drop is close to adiabatic. However, the increasing temperature difference between the well-mixed chamber interior and the chamber walls, which stay at almost constant temperature, causes an increasing heat flux from the walls to the chamber volume and thereby causes a deviation from the adiabatic temperature profile towards higher gas temperatures. The gas temperature was measured by five vertically and five horizontally aligned temperature sensors. These sensors are calibrated to an accuracy of 0.1 K by direct comparison to a certified platinum resistance temperature sensor (Lake Shore Cryotronics Inc., 120 sensor type PT-103-14L with calibration certificate). The sensor calibration runs are conducted in a chiller bath at temperatures between 293 K and 197 K. We account for an additional uncertainty due to the temperature variability throughout the chamber.
The analysis of measurement uncertainties is presented in more detail in Sect. 2.4.
A range of instruments was available for the measurement of gas phase and total water (Wagner et al., 2009). The water vapor pressure was measured inside AIDA by Tunable Diode Laser (TDL) absorption spectroscopy with a time resolution of 125 1 s using a laser wavelength in the near-infrared of 1370 nm (Fahey et al., 2014). The measurement is done by direct long path absorption in a plane 4.2 m above the bottom of the AIDA chamber. The accuracy of the water vapor measurement is better than ±5 % (Fahey et al., 2014). The total water concentration (which is the sum of water vapor and the amount of condensed and frozen water contained in aerosol and/or cloud particles) was measured with a dew point mirror (373LX, MBW) by sampling air from the AIDA chamber via a heated stainless steel tube (accuracy ±3 %). By subtracting the gas phase water concentration 130 from the total water concentration, the mass concentration of the condensed and frozen water can be calculated. The water vapor measurement by the AIDA TDL system has been compared with a range of other hygrometers during the AquaVIT-1 (intercomparison of atmospheric water vapor measurement techniques) campaign in 2007 and showed good agreement with established methods (Fahey et al., 2014). The ice saturation ratio, S ice , in the AIDA chamber is calculated based on the gas temperature, the TDL water vapor measurements and the saturation vapor pressure over ice given by Murphy and Koop (2005).

135
For some experiments, in situ Fourier Transform Infrared spectroscopy (FTIR) was used to deduce the composition and mass concentration of the H 2 SO 4 /H 2 O aerosol particles before and during expansion cooling (Wagner et al., 2008). Laser light scattering and depolarization measurements with the SIMONE instrument (Schnaiter et al., 2012) allowed for a precise determination of the freezing onset of the H 2 SO 4 /H 2 O aerosol particles during expansion cooling, as described in Sect. 2.3. In the SIMONE set-up, a polarized laser beam with an emission wavelength of 488 nm is directed horizontally across the AIDA 140 chamber. Photomultipliers detect the light scattered by the aerosol and/or cloud particles in the forward (2°) and backward (178°) direction. At 178°, the scattering intensities are measured polarization-resolved, enabling the determination of the backscattering linear depolarization ratio, which is zero for the spherical H 2 SO 4 /H 2 O solution droplets but increases as soon as aspherical ice crystals are formed after the freezing onset.
The number size distribution of aerosol and cloud particles was continuously measured by two optical particle counters 145 (OPC; model: welas, Palas GmbH), which were sampling via a vertical line from the bottom of the AIDA vessel. The size classification and measurement ranges of the two OPCs (named welas1 and welas2) depend on the particle shape and refractive index. For example, for spherical particles with a refractive index of 1.33, the size ranges of the two sensors are 0.7 to 46 µm (welas1) and 5 to 240 µm (welas2). In Sect. 2.3, it is described how the OPC records are analyzed to provide a second measure for the freezing onset of the H 2 SO 4 /H 2 O aerosol particles in addition to the SIMONE depolarization data. The H 2 SO 4 /H 2 O aerosol particles were characterized after their injection into the AIDA chamber. A CPC continuously measured the number concentration of aerosol particles, c aerosol with diameters in the size from 10 nm up to 3 µm inside the chamber. The size distribution of the particles was measured two times after the injection by APS and SMPS. The SMPS detects aerosol particles in a size range between 0.014 and 0.82 µm (mobility diameter), whereas the APS only measures larger 155 particles with diameters between 0.523 and 19.81 µm (aerodynamic diameter). As aqueous aerosol particles are spherical, the dynamic shape factor to convert the mobility diameter to a geometric diameter was set to unity. To transform the aerodynamic diameter to a geometric diameter, we applied the same shape factor together with the particle density, which was calculated with Model I of the E-AIM model (Clegg et al., 1992;Carslaw et al., 1995;Massucci et al., 1999;Wexler and Clegg, 2002;Clegg and Brimblecombe, 2005) for the specific equilibrium H 2 SO 4 /H 2 O composition. The combination of APS und SMPS data provides particle size distributions over the whole relevant size range. For most of the experiments, the measured size distribution could be well described and fitted by an unimodal lognormal distribution with the parameters N (total number concentrations), D m (median diameter) and σ (geometric standard deviation). Figure 2 shows the lognormal fits of the size distributions for all AIDA experiments for which APS/SMPS size distribution measurements were available. The median diameter ranges from about 30 nm to 1.1 µm. The measured size distribution of experiment number 14 of the ICE01 campaign (grey 165 dots) shows the typically good agreement between the measurement and the lognormal fit. For several experiments with a start temperature below about 210 K, we observed continuously increasing aerosol particle concentrations in the CPC measurements after the H 2 SO 4 aerosol injection had been stopped. The injected H 2 SO 4 vapor further nucleated small aerosol particles in the cold environment in the AIDA chamber, creating a second mode in the number size distribution at diameters between 10 nm and 30 nm. As these nanometer-sized particles from the occasional second nucleation mode retain a higher solute concentration 170 due to the Kelvin effect and have a lower freezing probability due to their small volume compared to the H 2 SO 4 /H 2 O particles of diameters > 100 nm in the principal nucleation mode, they do not significantly contribute to the observed homogeneous ice nucleation modes during expansion cooling.

Determination of freezing onsets
The freezing onsets were determined based on ice crystal number concentration measurements by the OPCs and the depolar-175 ization ratio records by SIMONE. We illustrate our approach using the measured time series of experiment number 4 of the TROPIC04 campaign as an example, see Fig. 3. At time t = 0 the expansion started with pressure and temperature immediately decreasing (black and red lines, upper panel). The relative humidity with respect to ice computed for gas phase and total water increased accordingly (middle panel). The H 2 SO 4 /H 2 O aerosol particles took up water in order to maintain thermodynamic equilibrium with the environment, which is visible in the faster increase in total water in comparison with the gas phase. This 180 observation was supported by the SIMONE light scattering measurements (bottom panel), in which the forward scattering intensity (2°, dark red) also continuously increased with the beginning of expansion cooling, which is due to the increasing size of the aerosol particles due to the water uptake. The depolarization ratio (red), however, remained constant at the zero background level in this initial time period, which indicates that the particles still had a spherical shape, meaning they were in liquid phase. After about 280 s of pumping, a sharp increase was observed in the depolarization ratio, which is indicative for 185 the formation of aspherical ice crystals in the chamber (red dashed line). The forward scattering signal also showed a strong increase at the same time because the nucleated ice crystals rapidly grow to large sizes with much higher scattering cross sections compared to the unactivated H 2 SO 4 /H 2 O aerosol particles. A few seconds later, the activated fraction N ice /c aerosol , relating the number concentration of nucleated ice crystals, N ice , to c aerosol from the CPC data, sharply increased, indicating the detection of ice particles in the OPC measurements (blue dashed line). The temperature and ice saturation ratio at the times 190 derived in this way denoted the ice onset conditions T ice and S ice,f r .
For the derivation of N ice from the OPC measurements, it was feasible for most experiments to define a specific optical threshold size to distinguish between smaller-sized H 2 SO 4 /H 2 O aerosol particles and larger sized ice crystals. This procedure was applied to experiments at starting temperatures above about 200 K, where the ice crystals rapidly grew to much larger sizes in comparison to the aerosol population. The optical size threshold was typically set to values between 2 and 10 µm, depending 195 on the aerosol particle size distribution. All particles larger than this size threshold were counted as ice, all smaller particles as H 2 SO 4 /H 2 O aerosol particles. For experiments at temperatures below about 200 K, ice crystal growth is comparatively slow and the sizes are much smaller than those of ice crystals formed at higher temperatures. In these low temperature experiments, N ice was therefore evaluated by analyzing the relative increase of the number concentrations in each size bin of the size distribution measurements instead of defining a size threshold. The ice number concentration and the resulting activated fraction 200 derived from both methods show a sharp increase as soon as ice was formed, and we define the time of this sharp increase as the time of freezing onset. An estimate for the uncertainty of the freezing onset time is given in Sect. 2.4.
Regarding the SIMONE measurements, the onset of ice nucleation was defined as the start of the sharp increase in the backscattering linear depolarization ratio and the overall back-scattering signals. As SIMONE measures the particle phase in situ, this method is not influenced by any possible artefacts induced by sampling. The time corresponding to the sharp increase of 205 the depolarization ratio was evaluated visually from the individual experimental time series. Figure A1 shows a comparison of the ice onset conditions determined from the OPC and SIMONE measurements. Both the values for T ice and S ice,f r agree well for all analyzed experiments, underlining the consistency between the SIMONE measurement, where the nucleated ice crystals are detected in situ in the AIDA chamber, and the OPC measurement, where the ice crystals are sampled from the chamber.
In general, the freezing onsets determined in the AIDA chamber are sensitive to nucleation rate coefficients between J V = 210 5·10 8 cm 3 s −1 and J V = 10 13 cm 3 s −1 , as reported in Möhler et al. (2003). These nucleation rate coefficients were determined for sulfuric acid aerosol particles with diameters between 0.5 and 2 µm with an estimated freezing probability between 0.02 and 0.5. For the definition of the freezing probability, see Equ. 2 in Möhler et al. (2003). The time delay for the detection of the nucleation peak ranged between 1 and 10 s .

Analysis of uncertainties 215
In this Section, we will discuss the sources of uncertainties for the determination of the freezing onset conditions T ice and S ice,f r .
Regarding the freezing onset temperature T ice , there are two major sources of uncertainties influencing its determination.
The first source of uncertainty is related to the onset time determined from the OPC and SIMONE records, which is described in detail in Sect. 2.3. The comparison shown in Fig. A1 has illustrated that both methods to determine the freezing onset conditions 220 agree very well. Moreover, given the steepness of the N ice and depolarization ratio curves at the onset of ice nucleation, the experimental error in the determination of the freezing onset time is assumed to be small. We determined an uncertainty of ±3 s for the onset time. Based on estimates for the diffusion length of water molecules in the H 2 SO 4 /H 2 O aerosol particles at low temperature, which are discussed in greater detail in Sect. 3, we assume that this uncertainty covers the time the initial ice nucleus needs to grow into the detection range of the OPCs. The variability of the gas temperature T in this 6 s time 225 interval is a first factor contributing to the uncertainty of T ice . The second source of uncertainty in the determination of T ice is the temperature distribution inside the chamber. Five vertically and five horizontally distributed temperature sensors inside the AIDA chamber were considered to account for temperature inhomogeneity for each experiment. Prior to the experiment start, the deviation of the temperature measurements of all sensors is usually below 0.1 K. In the course of the experiment, the measured temperatures start to deviate. Especially the deviation in the vertical is increasing, whereas the horizontal deviation 230 remains comparably small. Depending on the spread in vertical temperatures in the 6 s time interval around the freezing onset, we account for an additional uncertainty in T ice . As the deviation of the temperature sensors at the time of freezing onset varies from experiment to experiment, the uncertainty for every single experiment were analyzed individually. A typical deviation at the time of the freezing is about ±0.5 K.
Regarding the ice saturation ratio at freezing onset S ice,f r , there are three major sources contributing to the uncertainty in 235 its determination. The first source is the uncertainty in the determination of the freezing onset time from OPC and SIMONE measurements, which is assumed to be ±3 s, as described above. For calculating S ice,f r , the measured water vapor pressure p w is divided by the water saturation vapor pressure with respect to ice p sat,ice (T ). For p sat,ice , the parameterization given by Murphy and Koop (2005) was used. As this parameterization depends on temperature, the temperature uncertainty including the temperature distribution inside the chamber in the time interval of 6 s affects the uncertainty of S ice,f r via this parameter and 240 denotes the second factor contributing to the uncertainty in S ice,f r . The third source of uncertainty is the direct measurement of p w by the TDL hygrometer, with a general measurement accuracy of ±5 %. It is assumed that the TDL has a slight tendency to underestimate water vapor content at higher temperatures and to overestimate at lower temperatures (T < about 220 K) (Fahey et al., 2014). For a reasonable correction of the TDL measurements especially at low temperatures, a specific time interval after pumping was stopped and a dense cloud persisted in the chamber was considered. During this time, the interior of the chamber 245 is assumed to be saturated with respect to ice. We can therefore use this time period as a reference to correct our measurements to reasonable values. The highest correction factors were applied on the lowest temperature experiments, which resulted in a shift of up to 0.17 S ice towards lower values. For the highest temperatures, the correction factor was often near zero or only a few percent. Finally, we account for the uncertainty in S ice,f r by considering the three sources of uncertainty individually determined for each experiment.

Overview of AIDA campaigns and experiments
All experiments evaluated for this work are listed in Tables B1, B2 and B3, together with the initial pressure, p 0 , and temperature, T 0 , before the start of the expansion, as well as the temperature and ice saturation ratio, T ice and S ice,f r , and the associated cooling rate, cr, and vertical velocity, v, at the derived homogeneous freezing onset. In Figure B1, six histograms show the range and frequency of experimental conditions for all considered AIDA experiments. Except for two experiments, 255 all expansions were started at ambient pressure around 1000 hPa. Note that in the atmosphere, cirrus clouds are formed at much lower pressures. This, however, is not expected to be crucial for investigating the freezing onset conditions in terms of T and S ice . With various start temperatures between 190 K and 236 K, our experiments cover a broad range of temperatures relevant for homogeneous freezing processes. Most of the experiments started between 210 K and 215 K, but with the series of TROPIC campaigns conducted in 2019 and 2020, we have significantly extended the data set at temperatures below 200 K. The 260 majority of cooling rates determined at the associated freezing onset was around 1 K min −1 , which corresponds to a vertical velocity of about 1.67 m s −1 assuming dry adiabatic conditions. Only a few particular experiments had higher cooling rates up to 5 K min −1 (about 8.33 m s −1 ). The injected H 2 SO 4 /H 2 O aerosol particle population also covered a broad range of characteristic parameters with mean diameters between 30 nm and 1.1 µm and number concentrations from 30 cm −3 up to several thousand cm −3 . Finally, the weight percent composition of the H 2 SO 4 /H 2 O solution particles, wt%, at the beginning of each 265 experiment ranged between 32 % and 41 %. As the aerosol particles are assumed to be in thermodynamic equilibrium with the environment before starting the expansion, Model I of the E-AIM model (Clegg et al., 1992;Carslaw et al., 1995;Massucci et al., 1999;Wexler and Clegg, 2002;Clegg and Brimblecombe, 2005) was used to calculate the equilibrium weight percent composition at the measured temperature and relative humidity conditions inside the chamber.

270
The homogeneous freezing onsets in AIDA were determined as described in Sect. 2.3 using the OPC ice count and SIMONE data. Figure 4 shows the freezing onsets color-coded for the respective AIDA campaign and the corresponding year. For comparison, the WAC-based temperature dependent Koop2000 homogeneous freezing onset lines for nucleation rate coefficients of J V = 5 · 10 8 cm −3 s −1 (dotted line) and J V = 10 13 cm −3 s −1 (dashed line) are shown , as well as the liquid water saturation ratios computed using the parameterization by Murphy and Koop (2005)  Converting the solute concentration into water activity resulted in a close match of the freezing temperatures for the different solutes. It was therefore suggested to formulate the freezing nucleation rate coefficients of aqueous solution particles only as function of the water activity and the temperature. More recent studies suggested that the homogeneous nucleation rate 285 coefficients are about 2 orders of magnitude lower than those given by Koop et al. (2000) (Knopf and Rigg, 2011;Riechers et al., 2013). In terms of water activity, the experimentally derived freezing curve by Knopf and Rigg (2011) is about 0.01 lower than that predicted by Koop et al. (2000). This slightly delayed ice nucleation onset, however, does not account for the much larger deviation of the AIDA results to Koop2000. In the present study, the homogeneous freezing of water with H 2 SO 4 as the solute was measured as a function of the steady increase of the S ice and the related increase of the water content or 290 the dilution of the aerosol particles during the AIDA experiments. Thus, the observed differences of the AIDA homogeneous freezing onsets and the Koop2000 derived values may be related to the conversion of the water activity a w of the solute used by Koop et al. (2000) to S ice used in our studies as the basic water concentration unit in the aerosol system, or the assumption of thermodynamic equilibrium composition of the aerosol particles at the time of freezing onset. Therefore, the measurement uncertainties of the water vapor concentration and the gas temperature are the first key aspect to 295 discuss and interpret the AIDA data in comparison to Koop2000. In the present experiments, the water vapor concentration was directly measured in situ with a TDL setup (see Sect. 2.1) which was not yet available during the previous AIDA experiments Mangold et al., 2005). For the latter experiments, gas phase water concentration was derived from the total water measurements with the Lyman-α hygrometer FISH (Meyer et al., 2015) by subtracting the aerosol water content, which was calculated from the sulfuric acid aerosol mass concentration and the molar ratio of sulfuric acid and water in the 300 aerosol particles at thermodynamic equilibrium composition . The sulfuric acid mass concentration was determined by ion chromatography analysis of aerosol filter samples. The contribution of the aerosol water content to the total water content was below 1 % at temperatures around 235 K and increased to about 10 % at temperatures around 195 K. Therefore, this correction had only a minor contribution to the overall uncertainty of S ice,f r , which varied between ±0.08 and ±0.13 (see Table 1 in Möhler et al. (2003)). The uncertainty of the S ice,f r values derived in this manner varies from about 305 ±0.06 at the higher temperatures to about ±0.13 at the lower temperatures. Major contributions to this uncertainty come from the determination of the freezing onset time from OPC particle number and SIMONE scattering intensity and depolarization measurements, the direct measurements of the water vapor partial pressure, and the gas temperature uncertainty and chamber internal variability needed to calculate the water vapor saturation pressure with respect to the ice phase (see Sect. 2.4). Within these uncertainties, the two data sets from previous and more recent AIDA cloud simulation chamber experiments agree well 310 with each other, and show both a clear trend of increasing freezing onsets to lower temperatures in comparison to Koop2000.
The second important aspect for discussing the deviation of AIDA freezing onsets from Koop2000 is the assumption of thermodynamic equilibrium composition of the aerosol particles in the dynamic environment of changing temperature and relative humidity during the AIDA experiments. The process modeling results by Haag et al. (2003) showed that larger aerosol particles may have a slightly enhanced sulfuric acid concentration, in particular during experiments at high cooling rates and 315 low temperatures, but that the majority of smaller particles can well follow thermodynamic equilibrium in typical AIDA cloud expansion experiments. A similar result was obtained by Wagner et al. (2008) with a thorough analysis and application of the process model MAID  for AIDA experiments with aqueous sulfuric acid aerosol particles performed at temperatures between 230 K and 205 K. For lower temperatures, the viscosity of sulfuric acid increases with decreasing temperature (Williams and Long, 1995) and it can even undergo a transition to a glassy state (Koop, 2004). An enhanced 320 viscosity or a glassy state may reduce the water uptake of the aerosol particles and therefore particle growth and dilution. Such effects were discussed to e.g. inhibit homogeneous freezing of secondary organic aerosol particles at T < 200 K (Fowler et al., 2020). For aqueous sulfuric acid solutions, Koop (2004) summarized the conditions for transitions to ultraviscous and glassy states in dependence on temperature and aerosol composition. To get information about the phase state of the H 2 SO 4 /H 2 O aerosol particles in the AIDA experiments, we used Model I of the E-AIM model (Clegg et al., 1992;Carslaw et al., 1995;325 Massucci et al., 1999;Wexler and Clegg, 2002;Clegg and Brimblecombe, 2005) to calculate the sulfuric acid weight percentage for the temperature and humidity conditions before experiment start and at the conditions of the observed ice onsets (see Tables   B1, B2 and B3). For all the experiments, the calculated weight percentages are above the values for a transition to ultraviscous or glassy particles according to Koop (2004) and references therein (see Fig. D1). Especially the glass transition occurs at significantly lower temperatures compared to the AIDA starting and ice onset conditions. This is also supported by the aerosol particle forward scattering intensity measurements with the SIMONE instrument, which showed no evidence of delayed water uptake that would be expected if initially glassy particles were gradually transformed to aqueous solution droplets. To calculate the viscosity of the H 2 SO 4 /H 2 O aerosol particles in the AIDA chamber we applied the parameterization by Williams and Long (1995). This parameterization is only valid for temperatures > 200 K and for a H 2 SO 4 /H 2 O weight percentage between 30 and 80 wt% and can therefore only be applied to the starting conditions of AIDA experiments with T 0 > 200 K (see Fig.   335 D2, data points for experiments with T 0 < 200 K are displayed in a grey shaded box). The Stokes-Einstein-equation was then used to calculate the diffusion coefficient of water molecules in the aqueous sulfuric acid aerosol particles. From this, we then calculated the diffusion length of water molecules in the aerosol particles on a time scale of 6 s, which is much larger than the mean aerosol diameter. This indicates that the water diffusion is still fast enough to keep the aerosol particles in thermodynamic equilibrium.

340
In order to support our assumption that the aerosol particles freezing in AIDA cloud expansion experiments are at or at least close to thermodynamic equilibrium conditions, we performed an experiment starting at a gas temperature of about 197 K, as shown in Fig. 5a. This experiment was not performed with a constant pump rate as all other experiments discussed here, but with an intermediate reduction of the pump rate in order to maintain a constant relative humidity for a longer time period. The experiment was started as other experiments with decreasing temperature and increasing relative humidity, but 345 after the relative humidity passed the Koop2000 line after about 80 s, the pump rate was reduced to keep the relative humidity at an almost constant value above the Koop2000 line for about five minutes to give the particles enough time for reaching thermodynamic equilibrium. The forward scattering intensity of the SIMONE measurements (dark red) shows a clear increase at the beginning of the experiments, indicating the growing of the aerosol particles due to water uptake with increasing relative humidity. As soon as the relative humidity is kept at a constant value, no further increase in the forward scattering intensity 350 was observed. Fig. 5b shows the relative humidity and the forward scattering intensity at the beginning of the experiment. The forward scattering intensity clearly increases during the first two minutes of the experiment, which is caused by water uptake as a result of increasing relative humidity, and results in a decreasing viscosity of the particles. As soon as the relative humidity is controlled to an almost constant value by the reduction of the pump speed, no further increase in the forward scattering intensity is observed. The signal even follows the slight variations in the relative humidity, and a delayed water uptake is 355 not indicated. This supports our assumption that the aerosol particles can well follow thermodynamic equilibrium conditions during an AIDA experiment at typical pump rates and related rates of cooling and relative humidity increase. As soon as the pump rate was increased again at about 450 s experiment time, the freezing onset was observed at a relative humidity similar to other experiments with constant pump rate. The fact that no ice formation was observed during the long time period with almost constant relative humidity well above the Koop2000 line supports our assumption, that the observed shift of the freezing onset to higher relative humidity is not caused by a delayed nucleation due to kinetic limitation of water diffusion to and uptake by the aerosol particles. The result of this experiment also indicates that the observed high freezing onsets may not be caused by a delayed growth of the pristine ice nucleus embedded in a highly concentrated solution layer as suggested in previous studies (Clapp et al., 1997;Bogdan et al., 2006;Bogdan, 2006;Bogdan and Molina, 2010;Bogdan et al., 2013).
In addition, the depolarization measurements with the SIMONE instrument and the signatures of the recorded FTIR spectra 365 show no indication for the formation of sulfuric acid hydrates (Nash et al., 2001), which could occur under the experimental conditions in the AIDA chamber (Koop et al., 1997). An impact of hydrate formation on the observed ice onsets, as shown e.g.
in the case of sulfuric acid tetrahydrate (Fortin et al., 2003), is therefore not expected.
Results from Wagner et al. (2008) and more recent measurements during the TROPIC04 campaign show that the mass fraction of sulfuric acid in the H 2 SO 4 /H 2 O particles measured by FTIR extinction spectroscopy right before the homogeneous 370 freezing onset tend to be slightly lower when compared to the composition of the freezing solutions shown in Koop et al. (2000), but still overlap in the range of uncertainties. This observation supports the assumption that the conversion from the wt%-T -space into the S ice -T -space is a potential explanation for the observed deviation to Koop2000. For this conversion, the water vapor saturation pressure with respect to the supercooled liquid water phase is needed, or more precisely, the ratio of the saturation pressures with respect to the supercooled liquid and ice phases. The descriptions for the liquid water saturation 375 pressures are rather uncertain (Koop, 2004), and existing parameterizations deviate from each other (e.g. Buck, 1981;Sonntag, 1994;Tabazadeh et al., 1997;Murphy and Koop, 2005;Nachbar et al., 2019). In this work, we are using parameterizations from Murphy and Koop (2005) to calculate both liquid water and ice saturation pressures, but in a recent study Nachbar et al. (2019) suggested a new parameterization for liquid water saturation conditions which indicates that Murphy and Koop (2005) may increasingly underestimate the liquid water saturation pressure with decreasing temperature below about 220 K. Jensen et al.

380
(2005) discussed that a shift of the liquid water saturation line to higher ice saturation ratios would also shift the homogeneous freezing onset in cirrus formation processes to higher ice saturation ratios. Higher liquid water saturation pressures would also reduce the difference between the AIDA and the Koop2000 freezing onsets for aqueous sulfuric acid aerosol particles. This is shown in Fig. 4, where Koop2000 is shown in combination with the liquid water saturation pressure parameterization of Murphy and Koop (2005) (black dashed and dotted lines, MK2005) and additional in combination with the more recent line 385 by Nachbar et al. (2019) (blue dashed and dotted lines, N2019). To fully explain the differences between the AIDA and the Koop2000 onsets, the liquid water saturation pressure would need to be even higher at low temperatures than suggested by Nachbar et al. (2019).
Finally, we do not have a solid explanation for the deviation of the H 2 SO 4 /H 2 O homogeneous freezing thresholds observed in AIDA experiments from the WAC-based homogeneous freezing lines. This deviation may be related to uncertainties in the 390 formulation of physicochemical properties at low temperatures, which are required for the conversion between the Koop2000 parameter space (a w -T -space) and the AIDA parameter space (S ice -T -space), as described above.
Based on the AIDA results and the discussion above, we provide a new fit line for homogeneous freezing of H 2 SO 4 /H 2 O aerosol particles directly measured as a function of S ice and T . When plotted the SIMONE and OPC freezing onsets in an ordinary least square fit routine the data points are fitted by with fitting coefficients a and b. The coefficients of the fit shown in Fig. 6a are a = −1.40 ± 0.05 and b = 390 ± 10 K.
This fit transferred into the S ice -T -space and compared to Koop2000 and the water saturation lines is shown in 6b (red line) with the range of fit uncertainty (red shaded area). The goodness of fit is R 2 = 0.92. Applying this new freezing threshold 400 for cloud formation processes in the atmosphere, the homogeneous freezing onset in cirrus formation would be shifted to ice saturation ratios of about 2.0 at temperatures around 185 K. A higher homogeneous freezing onset may explain the high ice saturation ratios occasionally reported for low temperatures in the upper troposphere in some field studies (Jensen et al., 2005;Lawson et al., 2008;Krämer et al., 2009;Krämer et al., 2020). However, it needs to be considered that the fit line only describes homogeneous freezing of H 2 SO 4 /H 2 O aerosol particles under laboratory conditions. Other aerosol species, which 405 could be relevant for homogeneous freezing processes in the atmosphere, are not taken into account. Application to atmospheric conditions therefore needs to be done with caution.
As homogeneously formed cirrus clouds usually have a higher optical depth due to higher ice crystal number concentrations and ice water content, the reflectivity of solar radiation is increased in comparison to heterogeneously formed cirrus (Lohmann et al., 2016). However, homogeneously formed cirrus clouds also reduce the outgoing longwave radiation (OLR) from below 410 to be emitted in space due to their higher optical depth (Lohmann et al., 2016). Consequently, a precise description of homogeneous freezing processes is crucial to understand cloud radiative effects in the present climate as well as in predictions of climate change. This may in particular be relevant for cirrus clouds in the cold tropical tropopause layer (TTL). TTL cirrus clouds have the highest occurrence frequency globally and the strongest radiative warming effect, which enhances their importance for the Earth's climate. The higher homogeneous freezing onsets would suppress the TTL cirrus cloud formation 415 resulting in a decreasing cloud fraction (Jensen et al., 2005;Schoeberl et al., 2016), which would influence the global radiation budget by impacting the OLR (Mitchell and Finnegan, 2009) and reflected shortwave solar radiation.
Besides radiative effects, the presence of cirrus clouds in the tropopause region, especially in the TTL, are assumed to influence the water vapor mass flux from the troposphere to the stratosphere (Jensen et al., 1996;Hartmann et al., 2001;Corti et al., 2006). Air passing the TTL is dehydrated by crossing the cold point in the tropopause layer (Brewer, 1949), which 420 then controls the stratospheric water vapor (Dinh and Fueglistaler, 2014). The supersaturations measured above the cold point tropopause suggest that the water exchange with the stratosphere is 10−20 % higher than expected (Rollins et al., 2016;Krämer et al., 2020), which might also be related to an enhanced freezing threshold (Schoeberl et al., 2016). To understand stratospheric water contents and related stratospheric chemistry, cirrus formation thresholds (Schoeberl et al., 2016), supersaturations and temperatures (Randel and Jensen, 2013;Randel and Park, 2019) in the tropopause region are of particular importance. Several 425 studies have investigated the relation between the amount of stratospheric water vapor and the properties of the TTL in climate change scenarios (Solomon et al., 2010;Riese et al., 2012;Randel and Jensen, 2013), which puts the importance of upper tropospheric cirrus formation parameterizations into a larger context.

Conclusions
In this study, we present and discuss a comprehensive set of homogeneous freezing measurements for H 2 SO 4 /H 2 O aerosol 430 particles. The experiments were conducted at the AIDA cloud simulation chamber and covered a wide range of temperatures, cooling rates, aerosol sizes and number concentrations. At temperatures above 205 K, the measured ice saturation ratios at homogeneous freezing onset S ice,f r agree with the water activity criterion (WAC) based predictions by Koop et al. (2000).
Towards lower temperatures, however, the AIDA homogeneous freezing onset results show an increasing deviation from Koop2000 towards higher ice saturation ratios. For temperatures between 205 K and 185 K, the WAC-based ice saturation 435 ratios increase from about 1.6 to 1.7, whereas the AIDA measurements show an increase from about 1.7 to about 2.0 in the same temperature range.
For the comparison with Koop2000, we assume the aerosol particles to be in thermodynamic equilibrium at freezing onset.
This assumption is justified by previous measurements and process model results Wagner et al., 2008). We also demonstrate that the enhanced freezing onset is not caused by a delayed water uptake and by that a delayed dilution during 440 the AIDA experiments with steady increase of S ice . For this we conducted an AIDA experiment during which S ice was kept at an almost constant value well above the Koop2000 line for about 5 minutes. No ice formation or further water uptake of the aerosol particles was observed during this time period, but when increasing S ice in the cloud chamber to higher values, the freezing onset occurred in agreement to results from experiments with steady increase of S ice .
Further, the water activities a w used by Koop2000 have to be converted into S ice directly measured during AIDA experi-445 ments for comparing AIDA results with Koop2000. This conversion requires a formulation for the saturation vapor pressures over supercooled liquid water p sat,liq , which are rather uncertain at low temperatures (Murphy and Koop, 2005). A recent study by Nachbar et al. (2019) suggested a new parameterization for p sat,liq , which, for temperatures below about 220 K, increasingly deviates from Murphy and Koop (2005) to higher values with decreasing temperature. Higher p sat,liq would shift the Koop2000 freezing onsets to higher ice saturation ratios and would therefore reduce the difference between the AIDA and 450 the Koop2000 freezing onsets. It was already discussed in Jensen et al. (2005) that a shift of the p sat,liq line to higher ice saturation ratios would also shift the homogeneous freezing onset in cirrus cloud formation processes to higher ice saturation ratios. A higher homogeneous freezing onset as derived from our experiments may also explain field observations of high clear-sky supersaturation, which should not occur according to the freezing thresholds predicted by Koop2000 (Jensen et al., 2005;Lawson et al., 2008;Krämer et al., 2009;Krämer et al., 2020 An empirical fit to the AIDA results with the form ln(S ice ) = a + 1 T · b with fit parameters a = −1.40 ± 0.05 and b = 390 ± 10 K describes the observed homogeneous freezing onsets of H 2 SO 4 /H 2 O aerosol particles in the chamber. This fit line is based on direct measurements of the freezing onset conditions at simulated cirrus formation conditions and provides an 460 isoline for homogeneous freezing of sulfuric acid aerosol particles for nucleation rate coefficients between 5 · 10 8 cm −3 s −1 and 10 13 cm −3 s −1 . The application of this fit line to atmospheric conditions requires further work on the physical behavior of H 2 SO 4 /H 2 O aerosol particles at low temperatures and on the involvement of other atmospheric aerosol particle types. Ongoing experiments in the AIDA cloud simulation chamber aim at investigating homogeneous freezing onsets of different solutes and at constraining the descriptions for liquid water saturation pressures to experimental results.  Sci., 66, 741-754, https://doi.org/10.1175/2008JAS2542.1, 2009 A new parameterization of H2SO4/H2O aerosol composition: Atmospheric implications, Geophys. Res. Lett., 24, 1931-1934, https://doi.org/10.1029/97GL01879, 1997: Homogeneous freezing nu-  (Clegg et al., 1992;Carslaw et al., 1995;Massucci et al., 1999;Wexler and Clegg, 2002;Clegg and Brimblecombe, 2005) to transfer this weight percentage data into water activity, which is assumed to be equal to the relative humidity (assumption of thermodynamic equilibrium). If ice saturation ratios were given, the parameterizations of Murphy and Koop (2005) for the saturation pressures of supercooled liquid water and ice were used to calculate water activities. Additionally, the melting point line according to Murphy and Koop (2005) (solid line) and the homogeneous freezing thresholds for two different nucleation rate coefficients according to Koop et al. (2000) (dashed and dotted lines) are shown. The grey shaded area is indicative for the uncertainties of these homogeneous freezing thresholds (±5 % in water activity) as given in Koop (2004).  as gas (red) and wall (dark red) temperature inside the AIDA chamber as a function of time since pump start at t = 0. Accordingly, the middle panel shows the course of relative humidity with respect to ice RHice and ice saturation ratio Sice of gas phase (blue) and total water (grey). The turquoise area marks the range of freezing onset saturation ratios according to Koop et al. (2000). In the lower panel, the activated fraction of H2SO4/H2O particles determined by OPC data (blue) and the forward scattering intensity (dark red) as well as the back-scattering depolarization ratio (red) derived from SIMONE light scattering measurements are presented. The time of homogeneous freezing onsets derived from the activated fraction (blue) and from the SIMONE measurements (red) are marked by the vertical dashed lines.   Fig. 3. For the experiment shown, started at about 197 K, the pump rate was controlled in such a way that the relative humidity with respect to ice stayed relatively constant for about 5 minutes at about 170 %, hence above the homogeneous freezing threshold suggested by Koop et al. (2000). For details on this experiment, see Sect. 3. Panel (b): Enlarged view of the relative humidity and the forward scattering intensity at the beginning of the experiment.  Murphy and Koop (2005) and Nachbar et al. (2019). The blue and grey shaded areas are indicative for the uncertainties of Koop2000 (±5 % in water activity) as given in Koop (2004). Table 1. Previous measurements of homogeneous freezing of H2SO4/H2O solutions. The columns show the name of the publication, the used measurement device, the sample volume V , the cooling rate cr (if available) and information about which parameters, weight percentage composition wt% or the relative humidity RH, were measured or derived to characterize the freezing onsets ('measured' = parameter was directly measured, 'given' = parameter was calculated e.g. using the E-AIM model or others Appendix A: Comparison of OPC and SIMONE derived freezing onsets Figure A1. Comparison of the homogeneous freezing onset temperature (Tice, panel a) and ice saturation ratio (S ice,f r , panel b) derived from the SIMONE and OPC data. The respective Pearson's correlation coefficients are R = 0.999 and R = 0.996. Figure B1. Overview of experimental conditions. The six histograms show the distribution and the parameter range of start pressure p0, start temperature T0, cooling rate cr, vertical velocity v, as well as the median aerosol particle diameter Dm, the aerosol particle number concentration c aerosol and weight percentage wt% of H2SO4 in the aqueous solution particles before experiment start.
Appendix B: Overview of AIDA experiments Table B1. Summary of H2SO4/H2O AIDA experiments part I. The considered experiments are listed together with the corresponding campaign, the experiment number, the date of the experiment, the start pressure p0, the start temperature T0, the start relative humidity with respect to ice RHi,0 and liquid water RHw,0, the freezing onset temperature Tice, the ice saturation ratio at freezing onset S ice,f r , the cooling rate cr and the corresponding vertical velocity v at the freezing onset.

Campaign
Exp No. date p0 (hPa) T0 (K) RHi,0 (%) RHw,0 (%) Tice (K) S ice,f r cr (K min −1 ) v (m s −1 ) Table B2. Summary of H2SO4/H2O AIDA experiments part II. The considered experiments are listed together with the corresponding campaign, the experiment number, the date of the experiment, the start pressure p0, the start temperature T0, the start relative humidity with respect to ice RHi,0 and liquid water RHw,0, the freezing onset temperature Tice, the ice saturation ratio at freezing onset S ice,f r , the cooling rate cr and the corresponding vertical velocity v at the freezing onset.

Campaign
Exp No. date p0 (hPa) T0 (K) RHi,0 (%) RHw,0 (%) Tice (K) S ice,f r cr (K min −1 ) v (m s −1 ) Table B3. Summary of H2SO4/H2O AIDA experiments part III. The considered experiments are listed together with the corresponding campaign, the experiment number, the date of the experiment, the start pressure p0, the start temperature T0, the start relative humidity with respect to ice RHi,0 and liquid water RHw,0, the freezing onset temperature Tice, the ice saturation ratio at freezing onset S ice,f r , the cooling rate cr and the corresponding vertical velocity v at the freezing onset.

Campaign
Exp No. date p0 (hPa) T0 (K) RHi,0 (%) RHw,0 (%) Tice (K) S ice,f r cr (K min −1 ) v (m s −1 ) Figure C1. Summary of homogeneous freezing measurements of H2SO4/H2O solutions in the aw-T -space. The homogeneous freezing onsets of sulfuric acid solution samples shown in Fig. 1 are complemented by the AIDA results from this study (red points). The results in the Sice-T -space were converted into the aw-T -space by assuming equilibrium conditions aw = RHw and by using the parameterizations for the water vapor saturation pressures with respect to ice and supercooled liquid water given by Murphy and Koop (2005). The uncertainties of the calculated aw values vary between ±0.036 and ±0.076.
Appendix C: AIDA results in a w -T -space Figure D1. Phase diagram of sulfuric acid adapted from Koop (2004). The phase diagram of sulfuric acid as given by Koop (2004) and references therein is shown and complemented by the calculated weight percentage concentrations of the sulfuric acid aerosol particles before the start of the AIDA experiments of this study (red dots) and at the observed ice onset during the experiments (blue dots).
Appendix D: Phase state of supercooled sulfuric acid aerosol particles diffusion coefficients (red triangles) calculated using the parameterization of Williams and Long (1995) and the Stokes-Einstein-equation.
Panel (b): Comparison of the mean particle diameter of the sulfuric acid particle population (blue diamonds) to the diffusion length of water molecules (black squares) on a time scale of 6 s. Note that the shown data points represent the starting conditions of the AIDA experiments.
As the parameterization of Williams and Long (1995) is only valid for T > 200 K, the data points of experiments at lower temperatures are displayed in a grey shaded box.