Oliver
Reich
,
Michael J.
Gleichweit
,
Grégory
David
,
Nicole
Leemann
and
Ruth
Signorell
*
Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, CH-8093, Zurich, Switzerland. E-mail: rsignorell@ethz.ch
First published on 15th February 2023
Sea salt aerosol is among the most abundant aerosol species in Earth's atmosphere, and its hygroscopicity is an important parameter to quantify its interaction with solar radiation. Conflicting values for the hygroscopic growth have been reported in the literature, which decreases the accuracy with which their impact on Earth's climate can be modelled. Here we report new values of the hygroscopic growth for a selection of salt compositions representative of atmospheric sea salt. These values are obtained from single optically trapped aqueous droplets with dry radii between 0.3 and 2 μm, using a recently developed method for single particle mass measurement in an optical trap. We compare our results to earlier studies and propose a way to reconcile the apparent discrepancies found in the literature. Within our studies, we also observe the crystallization of CaSO4·2H2O (Gypsum) during the drying of optically trapped sea salt droplets at significantly larger relative humidity of 65–68% than the main efflorescence relative humidity at 50%. This preceding transition occurred in the absence of any contact of the particle with a surface.
Environmental significanceSea salt aerosol is among the most abundant aerosol species in Earth's atmosphere. Its hygroscopic properties influence cloud formation and its interaction with solar radiation. Despite the continuing effort, its hygroscopicity has not yet been finally quantified, resulting in large uncertainties in the predicted aerosol optical depth over oceans. New values for the hygroscopic growth of sea salt aerosols obtained from mass measurements in a single particle optical trap now allow us to reconcile apparent discrepancies of previously reported values. This has the potential to improve the accuracy with which the impact of sea salt aerosols on the climate of Earth can be modelled. |
Sea spray aerosol is produced in maritime areas over oceans and seas, and is one of the most abundant species in the atmosphere.7,8 As such, sea spray has been the subject of many studies over the last couple of decades. Its properties such as chemical composition, formation mechanisms and size distributions are explained in detail in the literature.8–10 Sea spray is composed of inorganic salts, which are usually referred to as sea salt, organic and biological matter. The organic and biological matter varies significantly for different particle sizes and is relevant in particular in the lower submicrometer size range.11,12 In contrast, the composition of the sea salt component is very similar to the composition of the sea water from which the droplets are generated,8,13 and the oceans have a very similar chemical composition across the Earth.8 In this article, we focus on the sea salt component of sea spray aerosols.
Despite the continuing effort, the hygroscopicity of sea salt is not yet quantified with high confidence. In the submicron size range, the hygroscopic growth of an aerosol is typically measured in size selected particle ensembles using hygroscopic tandem mobility analyzers (HTDMAs).14 The hygroscopic size growth factor at RH 90%, that is, the ratio of the particle size at RH 90% and at low RH (RH < 10%), has been reported between 2.11 and 2.46,15–17 thus spanning a range of approximately 15%. For particles with sizes of a few micrometers, trapping and measuring the hygroscopic growth of single levitated particles using electro-dynamic balances (EDBs)18–20 have been widely used. In this size range, the hygroscopic growth of sea salt was for a long time believed to be similar to that of pure NaCl,8 as measured by Tang et al.21 In their work, Tang et al. determined the hygroscopic mass growth factor, that is, the ratio of the particle masses at given RH and RH = 0%, of several sea salt compositions and reported values at RH 90% of 7.0 for NaCl and 6.8 for the sea salts. Recently, these values were critically evaluated by Zieger et al.,17 who found a significantly reduced hygroscopic growth for submicron and micron sized sea salt particles. For the same size range as in the work of Tang et al., a hygroscopic mass growth factor of 5.24 at RH 90% was found using EDBs. In the same study, it was shown that such a reduction of the hygroscopic growth parameter in climate models would lead to a reduction in the predicted aerosol optical depth over oceanic areas of up to 15%. The present uncertainty in the available literature data thus significantly limits the accuracy with which the impact of atmospheric sea salt on Earth's climate can be modelled.
In this work, we report new measurements of the hygroscopic mass growth of single NaCl and sea salt particles as well as values for their densities as a function of RH. In contrast to previous studies, we use optical trapping to investigate single particles. The measurement of the hygroscopic mass growth and density is facilitated by an approach combining the measurement of the particle's size and molecular composition via Raman spectroscopy with a recently developed method for mass measurement of single optically trapped aerosol particles.22 While mass measurements with EDBs are typically limited to particles with sizes of a few micrometers and above, our approach enables us to report data also for the large submicron size range. In comparison with the values reported by Tang et al.21 and Zieger et al.17 for micrometer sized particles, we find that our data are very similar to the results of Tang et al. We also take a critical look at the definition of the reference dry mass used for the derivation of the hygroscopic mass growth factors, which has to be taken into account for their correct implementation into existing climate models. This is discussed in particular with regard to the presence of hydrates in sea salt under dry conditions, which might explain the apparent discrepancies observed in the literature.
Fig. 1 Counter-propagating optical tweezers setup. The purpose of the electro-optic modulator (EOM) and the polarizing beam splitter (PBS) is explained in the main text. The beams are expanded and later focused into the trapping cell using aspherical lenses (AL). An optical isolator (OI) protects the laser cavity from the trapping beams propagating backwards towards the source. Figure adapted from ref. 30 under Copyright 2021 Society of Photo-Optical Instrumentation Engineers (SPIE). |
The trapped particles investigated in this study are aqueous droplets generated from bulk solution using an atomizer (TSI, 3076). To study the hygroscopicity of atmospheric sea salt aerosol droplets, aqueous solutions with three different salt compositions were prepared and used as representatives: pure NaCl (>99%, Merck, 1.06404.1000), a mixture of NaCl + MgCl2·6H2O (Sigma Aldrich, 63068) + Na2SO4 (Sigma Aldrich, 238597) which we denote as sea water 1 (SW 1) and a commercial sea water solution (Sigma Aldrich, SSWS) which we denote as sea water 2 (SW 2). The specific ionic compositions of SW 1 and SW 2 are listed in Table 1. Note that the absolute ion concentrations in the bulk solutions do not correspond to the absolute ion concentrations in the trapped droplet, since the latter also depend on the RH inside the trapping cell, which governs the amount of water in the droplet. The relative concentrations in the bulk solution and the droplets, however, are the same.
SW 1 | SW 2 | Zieger et al.17 | Tang et al.21 | |
---|---|---|---|---|
Na+ | 10.45 (1.00) | 11.03 (1.00) | 10.87 (1.00) | (1.00) |
Cl− | 18.87 (1.83) | 19.84 (1.80) | 19.53 (1.80) | (1.71) |
Mg2+ | 1.37 (0.13) | 1.33 (0.12) | 1.32 (0.12) | (0.12) |
SO42- | 2.68 (0.28) | 2.77 (0.25) | 2.66 (0.25) | (0.25) |
Ca2+ | — | 0.42 (0.038) | 0.40 (0.037) | — |
K+ | — | 0.39 (0.036) | 0.35 (0.032) | — |
HCO3− | — | 0.15 (0.013) | 0.17 (0.016) | — |
Br− | — | 0.067 (0.006) | — | — |
BO33- | — | 0.029 (0.003) | 0.031 (0.003) | — |
Sr2+ | — | 0.013 (0.001) | 0.008 (0.00081) | — |
F− | — | 0.0014 (0.0001) | — | — |
The droplets are trapped inside a custom made cell, shown in Fig. 2, which enables us to control the droplets's environment. The cell is filled with nitrogen gas, and the RH inside the cell is controlled by mixing wet and dry nitrogen flows with adjustable flow rates. All experiments were conducted at room temperature (T = 21.7 ± 0.2 °C) and ambient pressure. Temperature and RH are monitored by using a sensor (Sensirion, SHT35) during the experiments.
Fig. 2 Trapping cell and measurement setups. The cell is custom made with an opening for the temperature and relative humidity sensor (T/RH). Particles are introduced into the cell via the inlet at the top (Pin). The nitrogen flow enters at the bottom of the cell (Nin) and exits together with the excess particles via an exhaust (Out). The scattered light of the trapping beams (TB1 and TB2) is collected through the back window. The scattered light beam is split, with half the power being focused on the position sensing photodiode (PSP) for the tracking of the particle position. The polarizing beam splitter (PBS) effectively removes one of the two scattered light beams and improves the accuracy of the particle tracking. The other half of the scattered light is fiber coupled into a spectrometer for the measurement of the Raman spectrum. To determine the size of the particle, light from a broad-band light source is focused on the trapped particle from the bottom window (BB source) and the scattered light is collected through the front window for the measurement of the broad-band light scattering (BLS) spectrum. Part of the figure adapted from ref. 30 under Copyright 2021 Society of Photo-Optical Instrumentation Engineers (SPIE). |
The BLS spectrum is recorded in the range 300–500 nm and compared with simulations using Mie theory. The fitting of the experimental BLS spectrum with the simulated spectra proceeds as described in detail in our previous work.22 From the fit, the radius and wavelength dependent refractive index of the particle are determined with a typical size accuracy of better than 1%.
Fopt = −k·(z(t) − z0(t)) = −kz(t) + kZ0sin(2πft) | (1) |
Fdamp = FStokes + Fcorr | (2) |
FStokes = −γ·ż | (3) |
(4) |
The optical force and the damping forces combined drive the particle dynamics according to
(5) |
(6) |
(7) |
(8) |
(9) |
(10) |
The phase φ of the particle oscillation is measured by tracking the particle position along the trapping beam axis during the trap modulation as shown in Fig. 2. The scattered light from the trapping beams is collected under a scattering angle of 90 ± 24° using an objective (Mitutoyo, Plan Apo 20×) with a relatively high numerical aperture (NA 0.42). This scattered light is split equally in power to allow simultaneous Raman measurements explained later, and one branch is focused on a position sensing photodiode (PSP, Thorlabs, PDP90A). The oscillation of the particle translates into an amplified oscillation of the focused image on the photosensitive area of the PSP, and the PSP signal is measured together with the modulation signal used to drive the EOM. From these two measurements, the demodulated amplitude and phase of the particle oscillation is obtained using a lock-in detection scheme. For accurate mass measurements, the artificial increase in the phase measurement due to the impedance of the PSP is determined and corrected as described in our previous work.22 Previously, an additional systematic measurement bias of up to 3.5% in the position was reported for micrometer sized particles. This bias arises due to the superposition of the two particle images from the respective trapping beam on the PSP. Using cross-polarized trapping beams and filtering the scattered light using a PBS effectively suppresses one of the particle images and thus removes this bias.
From Γ0 determined by the measurement of φ and the particle radius R determined by BLS, the mass of the particle is obtained using eqn (8)
(11) |
As an added benefit, the density of the spherical particles is obtained as
(12) |
Many examples of quantitative Raman spectroscopy for the determination of the sulfate mass fraction of a given aqueous solution can be found in the literature.34–37 In particular, the ratio of the peak heights of ν1-SO42− and νs-H2O modes and the ratio of their peak areas have both been shown to be proportional to the sulfate concentration of the solution to a good approximation.35,37 With respect to these earlier studies on bulk samples, the analysis of single particle Raman spectra is made more complex due to the presence of morphology dependent resonances, or whispering gallery modes (WMG),38 that arise due to the particle's spherical shape. More precisely, only Raman spectra where the ν1-SO42− and νs-H2O bands do not overlap with a WGM can be used for quantitative analysis. After these spectra have been identified by visual inspection, the peak heights of the ν1-SO42− and νs-H2O modes are determined and compared with the calibration curve, which is obtained from a dilution series of SW 1. The data used to obtain the calibration curve are shown in Fig. 3. From the standard deviation of the calibration curve with respect to the known mass fractions of the bulk solutions, we find an average uncertainty of 5.1% for the determination of the SO42− mass fraction. The calibration curve is also suited for the commercial solution with composition SW 2, since the two compositions differ only in terms of trace ion species in very low abundance.
The mass of SO42− is determined from the mass measurements and the calibrated Raman spectra of each particle. The total amount of salt is determined from the SO42− mass using the relative mass fractions listed in Table 1 for the respective composition. We refer to this salt mass as mdry, the dry mass of the particle. Different definitions of the dry mass are used in the literature; in particular, the dry mass is often used to describe the mass of the particle in an environment at 0% RH, or at a RH below a certain threshold. This is not necessarily the same mass as the one used in our case, due to the presence of hydrates that persist under dry conditions. We will comment on the implications of the different definitions of the dry mass at a later stage in this work. For now, we return to the definition of the hygroscopic mass growth factor, m*, which is obtained from the dry mass as
(13) |
The resulting densities are shown in Fig. 4a. For reference, the data from previous EDB work by Tang et al.21 are also shown and compared with bulk solution densities. The RH range above the critical value of 75.3% for NaCl41 corresponds to undersaturated droplets/solutions. In this range, the bulk values coincide with the EDB measurements. In contrast to bulk solutions, for RH below the critical RH the droplets remain homogeneous and liquid in a supersaturated state until efflorescence occurs at a RH of approximately 50%. In the supersaturated range, only the comparison with EBD results is possible. Overall, we observe good agreement between our measurement data and the results of Tang et al. To further analyze the size dependency of our accuracy, the data are grouped into three size categories: “submicron” indicates droplets with dry radii Rdry < 0.5 μm, “intermediate” indicates Rdry = 0.5–1.0 μm and “supermicron” Rdry > 1.0 μm. The reference data by Tang et al. are accordingly classified as supermicron. The dry size of each of our droplets is obtained as using the dry mass mdry and the density of pure NaCl, ρNaCl = 2170 kg m−3. The dry mass is determined from the mass measurements at undersaturation and the NaCl salt mass fractions obtained from the comparison with bulk solutions with the corresponding water activities. The dry mass of a given droplet is obtained as the average over all retrieved values for that droplet. For a small number of droplets, no suitable mass measurements were performed in the undersaturated RH range. For these droplets, the dry mass was determined from the hygroscopic mass growth factors at supersaturation measured by Tang et al.21 and treated as estimates by way of precaution. In the study presented here, droplets with dry sizes in the range Rdry = 0.3–2 μm have been investigated. Intermediate and supermicron droplets are observed to exhibit average measurement uncertainties of 8% and 4%, respectively, for the droplet density. This uncertainty rises to approximately 10% for submicron droplets.
Fig. 4 Hygroscopicity of aqueous NaCl droplets. (a) Density and (b) hygroscopic mass growth factor as a function of relative humidity (RH). Our data are divided into three groups according to the dry size of the corresponding droplets. For comparison, the bulk values and data of supermicron particles from previous EDB work by Tang et al.21 are also shown as solid and dashed lines, respectively (the data by Tang et al. are represented by a fourth order polynomial using the coefficients given in their work). The dotted line represents the interpolation of our data with a weighted polynomial fit. The critical RH of NaCl at 75.3% is indicated as cRH. Vertical error bars indicate the uncertainty arising from the droplet size and mass measurement. |
Compared with the results of Tang et al. we observe a small bias towards larger densities in our measurements for all droplet sizes. For a quantitative analysis we interpolate our data using a third order polynomial fit, weighted by the inverse square of the individual measurement uncertainties. The third order represents the best compromise between the determinacy and flexibility of the interpolation. On average our values lie 2–3% above the bulk densities in the undersaturated RH range, a value that is within the indicated experimental uncertainty of the setup. This small bias may be a residual artefact of the phenomenon reported in our previous work,22 where for certain particle size-to-laser wavelength ratios, the particle oscillation has been observed to stop at specific instances during the modulation cycle. This additional delay leads to an artificial increase in the phase φ in the harmonic oscillator model and thus to an increase of the mass and density. This phenomenon is prevalent in CPTs with parallel polarized trapping beams,30 and is therefore not significant in this work where we use perpendicular polarized beams, although a residual artefact of 2–3% cannot be ruled out.
The hygroscopic mass growth factors derived from the mass measurements and the determined dry masses are shown in Fig. 4b. For the droplets, where the dry mass could not be determined from the bulk solutions due to the lack of measurements at undersaturation, no hygroscopic mass growth factors are reported. Very good agreement with the reference bulk data and the data by Tang et al. is observed. It is common to report the hygroscopic growth factors at a RH of 90%. To determine the average hygroscopic mass growth factor, we interpolate our data with an empirical two parameter hygroscopic growth function
(14) |
The densities of SW 1 and SW 2 are shown in Fig. 5. For quantitative comparison, a weighted third order polynomial fit is performed on our density data as described above for the NaCl measurements. No significant difference of the densities between SW 1 and SW 2 can be observed. We observe slightly larger densities for RH below 90% for both sea salt solutions than the corresponding values by Tang et al. The average deviations between our fit and the data by Tang et al. are 3% and 2% for RH above 90%, 5% and 4% for RH in the range 80–90%, and 5% and 6% for RH < 80% for compositions SW 1 and SW 2, respectively. For the range RH < 90%, the discrepancy is larger than the bias observed in the NaCl data. At the same time, the experimental uncertainties for measurements with SW 1 and SW 2 are slightly smaller (3–7% at average) than those with NaCl due to overall better trapping alignment during the measurement series. Even taking into account the indicated measurement uncertainty, a small but systematic discrepancy in the droplet density of 1–2% remains for RH below 90% between our measurements and the results of Tang et al.
Fig. 5 Density of sea salt as a function of RH. (a) Data of composition SW 1. (b) Data of composition SW 2. The same size categories as for the NaCl measurements are applied. Data from previous EDB work by Tang et al.21 are shown for comparison (dashed line, polynomial interpolation). The dotted line represents the interpolation of our data with a weighted polynomial fit for quantitative analysis. The error bars are defined as for the NaCl measurements. |
In addition to the size and mass measurements, a Raman spectrum is recorded every two minutes for quantitative determination of the salt mass fraction. From this data and the droplet mass, the hygroscopic mass growth factors are determined as explained in the Methods section. These data are shown in Fig. 6. For both sea salt compositions SW 1 and SW 2, overall good agreement with the reference data by Tang et al. is observed. We observe somewhat smaller growth factors for RH > 85%. As for the NaCl data, we quantify the average hygroscopic mass growth factor by performing a weighted fit with eqn (14). From this interpolation, we obtain an average value of m*(90%) = 6.3 ± 0.3 for both compositions SW 1 and SW 2. The uncertainty range is given by an uncertainty of 5.1% of the calibration of the Raman spectra as specified in the Methods section. Our value is very similar to the one reported by Tang (6.8) and approximately 20% larger than the mass growth factors published by Zieger et al. (5.24 for supermicron particles). The similarities between SW 1 and SW 2 in both the densities and the hygroscopic mass growth factors show that the presence of the ions of low abundance (Ca2+, K+, HCO3−, Br−, BO33−, Sr2+, and F−) is not significant for the hygroscopic growth of sea salt particles.
Fig. 6 Hygroscopic mass growth factors of sea salt. (a) Data of composition SW 1. (b) Data of composition SW 2. The same size categories as for the NaCl measurements are applied. Data from previous EDB work by Tang et al.21 (solid line, polynomial interpolation) and Zieger et al.17 (diamonds) are shown for comparison. The dotted line represents the interpolation of our data (eqn (14)) and the cross indicates the value at RH 90% with the systematic uncertainty of the Raman calibration in the determination of the dry mass. The error bars of our measurement data are defined as for the NaCl measurements. |
The comparison with our data implies that in the case of Tang et al. there is only an insignificant amount of hydrates in the dried salt mass, as we do not observe larger hygroscopic mass growth values as would be expected otherwise. By contrast, our data even tend towards 3–13% lower values at RH 90%, which may be due to the slightly different sea salts used in the studies, in particular with regard to the Na+/Cl− ratio of the composition. Conversely, the hygroscopic mass growth factors of Zieger et al. are approximately 20% smaller. This is a plausible deviation if one considers the presence of hydrates in the dried sample of Zieger et al. If, as an example, in a salt crystal with ionic composition SW 2 the naturally occurring hydrates MgCl2·6H2O and Na2SO4·10H2O were formed instead of the anhydrous MgCl2 and Na2SO4, the salt mass would increase by 30% and the mass growth factors would decrease by the same relative amount. At this point, the question arises as to why Tang et al. and Zieger et al. obtain significantly different hydration states in their respective dry sample, even though they employ the same measurement principles. We are not familiar with all the details of either of the used setups, and hence our explanation here must remain speculative. It is, however, worth highlighting a few features of the experimental setups in question. Zieger et al. generated single droplets from salt solution using an ink-jet cartridge,17 while a custom made “particle gun”42,43 with similar functionality was used in the study by Tang et al. The particle generation methods are comparable; hence they are likely not related to the observed differences in the measurement data. Tang et al. initially pump the trapping chamber down to pressures of approximately 10−7 mbar42 before the dry mass is determined, whereas the pressure in the setup by Zieger et al. does not go below 150 mbar.44 Although previous studies suggest that vacuum pumping without heating is not efficient in removing hydrates from macroscopic samples,45 it is plausible that it has an enhanced effect on a single micrometer sized levitated sample. The absence of hydrates in the experiments by Tang et al. could therefore be related to their removal by vacuum pumping prior to the measurement of the dry mass.
The above considerations underline the importance of the definition of the salt mass in quantifying the hygroscopic growth of atmospheric particles. In our study, the anhydrous salt mass is used as reference. In studies, where it is not possible or not desired to use the anhydrous mass as reference, quantitative determination of hydrates in a dried sample is suggested as a method to improve comparability with other literature data. Calibrated Raman spectroscopy may be used to this end, in analogy to the determination of the dry mass of a droplet presented earlier; however, this approach was not investigated further in this study. Such complementary data increase the transferability and overall reliability of determined values and might reconcile apparent discrepancies in hygroscopic growth factors reported in the literature.
(a) 8.8 g of Na+/kg H2O, 15.0 g of Cl−/kg H2O, 1.06 g of Mg2+/kg H2O, and 2.22 g of SO42−/kg H2O
(b) 88 g of Na+/kg H2O, 150 g of Cl−/kg H2O, 10.6 g of Mg2+/kg H2O, and 22.2 g of SO42−/kg H2O
Sea salt droplets made from sea water samples (a) and (b) will have exactly the same hygroscopicity, because the different ions have the same relative abundance. At any RH, the thermodynamic equilibrium of the droplet, which determines the equilibrium concentrations of the ions and water in the droplet, is only dictated by this relative abundance of the ions (which remains the same between the sea water and droplet at any RH). Hence droplets made from two bodies of water with different salinities will still have the same mass growth factor and size growth factor, as long as they have the same relative abundance of ions. Droplets made from sea water with slightly different relative abundance of ions may have a slightly different thermodynamic equilibrium and hence slightly different hygroscopicity. However our results show that small changes of the relative abundance of ions in the sea water, like between SW 1 and SW 2, have no significant influence on the particle hygroscopicity. Even NaCl and sea salt droplets have similar hygroscopic growth. Hence the hygroscopicity measurements presented in this paper appear to be representative for most inorganic sea salt droplets when the influence of hydrates is properly taken into account. Inorganic sea salt droplets may experience chloride depletion due to aging processes as explained by Su et al.,46 which are not considered here. The effects of such aging on the hygroscopicity of individual inorganic sea salt droplets have only recently begun to be studied.47 Rosati et al.‘s results47 suggest that the chloride depletion mostly affects nanometer size particles by increasing the amount of hydrates in dry particles. Therefore, chloride depletion may not affect the overall hygroscopicity of sea salt as much, especially if the influence of hydrates is carefully considered, as suggested in this article.
For particles with composition SW 1 the molecular peaks are observed to remain qualitatively the same down to a RH of approximately 50%, where efflorescence occurs. By contrast, particles with composition SW 2 undergo a preceding transition at a RH of approximately 68%, where ν1-SO4 at 976 cm−1 splits into two peaks at 976 cm−1 and 1007 cm−1. Every investigated particle with composition SW 2 showed the same behavior. The preceding transition was observed between 65 and 68% RH. When adding an amount of 0.12 g per 100 ml of CaCl2 to composition SW 1 and repeating the measurement series, the preceding transition occurred, indicating that the origin of this transition is linked to the presence of Ca2+ ions.
By comparison of the Raman signal of different CaSO4 hydrates,48 the emerging peak centered at 1007 cm−1 can be assigned to crystalline CaSO4·2H2O (Gypsum). The residual peak at 976 cm−1 is attributed to the remaining SO42− in the aqueous phase that does not partake in the crystallization reaction (there is a large excess of SO42− with respect to Ca2+ ions in the sea salt composition). The crystallization of CaSO4·2H2O has a negligible effect on the density and mass growth of the particle. The transition is not noticeable in the respective data sets (Fig. 5 and 6). This agrees with the observation that sea salt compositions SW 1 and SW 2 show very similar densities and hygroscopic mass growth factors, implying that ions in lower abundance play a negligible role for these quantities. As such, the removal of the Ca2+ ions and the corresponding amount of SO42− from the aqueous phase is not expected to have a noticeable impact on either size, mass or density of the particle. An upper estimate of the effect of the phase transition on the density can be obtained by assuming that all Ca2+ ions take part in the crystallization process. The mass fraction of crystallized CaSO4·2H2O with respect to the whole anhydrous salt mass is then approximately 4.5%. Taking into account the experimentally determined mass growth factor of 3.0 at the transition RH ≈ 65% and the density of CaSO4·2H2O of approximately 2300 kg m−3, roughly twice the density of the aqueous solution, the change in particle density upon phase transition can be estimated to be approximately 0.75%. The change in density during phase transition is therefore lower than our experimental uncertainty.
The efflorescence RH (ERH) of sea salt is reported to be approximately 45%.17,21 Experimentally, we find our sea salt droplets to go through the main efflorescence phase transition at somewhat larger values at approximately 50% (see Fig. 7). Regardless of the exact ERH, CaSO4·2H2O crystallization is observed at considerably larger RH. Experimentally, we therefore observe efflorescence in two steps at distinct RH, with the main step being the efflorescence transition at approximately 50% preceded by the crystallization of CaSO4·2H2O at around 65%. Similar phase transitions of sea salt solutions in multiple steps have previously been studied for droplets on a substrate.49,50 For samples containing Ca2+ ions, the crystallization of CaSO4 hydrates was observed at RH as large as 90%.49 The presence of a substrate implies the existence of a contact surface for heterogeneous nucleation, which is known to facilitate phase transitions at larger RH,51 but which is absent for the levitated droplets in this work. Heterogeneous nucleation processes are thought to be suppressed in the absence of a contact surface, and the associated multiple step phase transitions have not been previously reported for levitated particles to the best of our knowledge. Our findings however show that such transitions may occur even for levitated droplets with no apparent contact surface.
We note that the multistep transitions of sea salt droplets discussed here are to be distinguished from similarly denoted processes, which describe a single phase transition for which one or several precursor states have been identified,52–54 or from phase transitions which proceed over one or several intermediate states to the final product as governed by the Ostwald step rule.29,55,56 Such multistep transitions have been reported for levitated particles and other systems and, while of great interest for fundamental studies of nucleation processes, are beyond the scope of the present work.
Finally, in sea water compositions containing Ca2+ and SO42−, the crystallization of CaSO4·2H2O (Gypsum) was observed at a RH of 65–68%, which is significantly larger than the RH at which the efflorescence of the particle occurred (RH ≈ 50%). This crystallization is typically attributed to heterogeneous nucleation processes that occur when the particle is in contact with a surface, such as a substrate or the surface of a secondary particle. By contrast, we observe the crystallization of CaSO4·2H2O in the absence of any apparent contact surface.
This journal is © The Royal Society of Chemistry 2023 |