Thermodynamic Dissociation Constant of the Bisulfate Ion from Raman and Ion Interaction Modeling Studies of Aqueous Sulfuric Acid at Low Temperatures

Cations and Anions Affect the Speed of Sound in Water Oppositely

Identifying the composition of a solution using acoustics remains a challenge. It is known that for low salt concentrations the speed of sound in water increases linearly with the concentration of the electrolyte, but the contribution of individual cations and anions is unknown. We introduce the concept of intrinsic sound speed Ai to quantify the contribution of ions to the speed of sound. We found that cations increase the speed of sound in water whereas anions decrease the speed of sound. Hydration layers around the ions play a major role. Because cations have a hydration layer thicker than that of anions, their contribution to the speed of sound is larger than that of anions. Experimental data on salts not used to determine the contribution of individual ions are in quantitative agreement with the predicted values. Our method can be applied to various systems containing small quantities of ions, molecules, or particles. With the knowledge that cations increase the speed of sound, we were able to explain previously unexplained data in the literature.

Effect of Acidity on Ice Nucleation by Inorganic-Organic Mixed Droplets

Aerosol acidity significantly influences heterogeneous chemical reactions and human health. Additionally, acidity may play a role in cloud formation by modifying the ice nucleation properties of inorganic and organic aerosols. In this work, we combined our well-established ice nucleation technique with Raman microspectroscopy to study ice nucleation in representative inorganic and organic aerosols across a range of pH conditions (pH -0.1 to 5.5). Homogeneous nucleation was observed in systems containing ammonium sulfate, sulfuric acid, and sucrose. In contrast, droplets containing ammonium sulfate mixed with diethyl sebacate, poly(ethylene glycol) 400, and 1,2,6-hexanetriol were found to undergo liquid-liquid phase separation, exhibiting core-shell morphologies with observed initiation of heterogeneous freezing in the cores. Our experimental findings demonstrate that an increased acidity reduces the ice nucleation ability of droplets. Changes in the ratio of bisulfate to sulfate coincided with shifts in ice nucleation temperatures, suggesting that the presence of bisulfate may decrease the ice nucleation efficiency. We also report on how the morphology and viscosity impact ice nucleation properties. This study aims to enhance our fundamental understanding of acidity's effect on ice nucleation ability, providing context for the role of acidity in atmospheric ice cloud formation.

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (∼1-6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0-18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6-7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes, and aqueous sulfate solutions also play a critical role as electrolytes for many prominent battery chemistries. While the thermodynamic literature on it is quite well-developed, comprehensive studies of the solvation structure, particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear nuclear magnetic resonance (NMR) approach to the elucidation of the solvation structure and dynamics over wide temperature (-10 to 50 °C) and concentration (0-18 M) ranges, combining the ¹⁷O shift, line width, and T₁ relaxation measurements, ³³S shift and line width measurements, and ¹H pulsed-field gradient NMR measurements of proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (∼4.5 M), with the high-concentration regime dominated by a strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H₂O/H₃O⁺ and H₂SO₄/HSO₄^-/SO₄^{2- 17}O resonances and the ¹H self-diffusivity. However, the ¹⁷O chemical shift difference between the H₂O/H₃O⁺ and H₂SO₄/HSO₄^-/SO₄^2- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling molecular dynamics simulations and density functional theory NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR, and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Inline Raman Spectroscopy and Indirect Hard Modeling for Concentration Monitoring of Dissociated Acid Species

We propose an approach for monitoring the concentration of dissociated carboxylic acid species in dilute aqueous solution. The dissociated acid species are quantified employing inline Raman spectroscopy in combination with indirect hard modeling (IHM) and multivariate curve resolution (MCR). We introduce two different titration-based hard model (HM) calibration procedures for a single mono- or polyprotic acid in water with well-known (method A) or unknown (method B) acid dissociation constants pK_a. In both methods, spectra of only one acid species in water are prepared for each acid species. These spectra are used for the construction of HMs. For method A, the HMs are calibrated with calculated ideal dissociation equilibria. For method B, we estimate pK_a values by fitting ideal acid dissociation equilibria to acid peak areas that are obtained from a spectral HM. The HM in turn is constructed on the basis of MCR data. Thus, method B on the basis of IHM is independent of a priori known pK_a values, but instead provides them as part of the calibration procedure. As a detailed example, we analyze itaconic acid in aqueous solution. For all acid species and water, we obtain low HM errors of < 2.87 × 10^-4mol mol^-1 in the cases of both methods A and B. With only four calibration samples, IHM yields more accurate results than partial least squares regression. Furthermore, we apply our approach to formic, acetic, and citric acid in water, thereby verifying its generalizability as a process analytical technology for quantitative monitoring of processes containing carboxylic acids.

Predictive models for the phase behaviour and solution properties of weak electrolytes: nitric, sulphuric, and carbonic acids

The distribution of ionic species in electrolyte systems is important in many fields of science and engineering, ranging from the study of degradation mechanisms to the design of systems for electrochemical energy storage. Often, other phenomena closely related to ionic speciation, such as ion pairing, clustering and hydrogen bonding, which are difficult to investigate experimentally, are also of interest. Here, we develop an accurate molecular approach, accounting for reactions as well as association and ion pairing, to deliver a predictive framework that helps validate experiment and guides future modelling of speciation phenomena of weak electrolytes. We extend the SAFT-VRE Mie equation of state [D. K. Eriksen et al., Mol. Phys., 2016, 114, 2724-2749] to study aqueous solutions of nitric, sulphuric, and carbonic acids, considering complete and partially dissociated models. In order to incorporate the dissociation equilibria, correlations to experimental data for the relevant thermodynamic equilibrium constants of the dissociation reactions are taken from the literature and are imposed as a boundary condition in the calculations. The models for water, the hydronium ion, and carbon dioxide are treated as transferable and are taken from our previous work. We present new molecular models for nitric acid, and the nitrate, bisulfate, sulfate, and bicarbonate anions. The resulting framework is used to predict a range of phase behaviour and solution properties of the aqueous acids over wide ranges of concentration and temperature, including the degree of dissociation, as well as the activity coefficients of the ionic species, and the activity of water and osmotic coefficient, density, and vapour pressure of the solutions. The SAFT-VRE Mie models obtained in this manner provide a means of elucidating the mechanisms of association and ion pairing in the systems studied, complementing the experimental observations reported in the literature.

Thermodynamic Modeling of Calcium Sulfate Hydrates in a CaSO4-H2SO4-H2O System from 273.15 to 473.15 K up to 5 m Sulfuric Acid

To prevent scaling and to recycle aqueous solutions in industrial processes, the thermodynamic properties of the CaSO₄-H₂SO₄-H₂O system are studied by thermodynamic modeling with the Pitzer model. The published solubility data of calcium sulfate hydrates in sulfuric acid solutions were collected and reviewed critically. Then, the CaSO₄-H₂SO₄-H₂O system was modeled using the Pitzer activity coefficient approach from critically selected experimental data to obtain optimized parameters. The model reproduces the solubility data with good accuracy up to 5 m sulfuric acid at temperatures of 283.15-368.15, 283.15-473.15, and 298.15-398.15 K for gypsum (CaSO₄·2H₂O), anhydrite (CaSO₄), and hemihydrate (CaSO₄·0.5H₂O), respectively. However, at temperatures above 398.15 K and sulfuric acid concentration above 0.5 mol/kg, the solubility of anhydrite predicted by our model deviates significantly from the literature data. Our model predicts that the solubility of anhydrite would first increase but then decrease in more concentrated sulfuric acid solutions, which is in disagreement with the experimental data showing constantly increasing solubilities as a function of increasing sulfuric acid concentration. This discrepancy has been discussed. The transformations of gypsum to anhydrite and hemihydrate were predicted in sulfuric acid solutions. With increasing H₂SO₄ concentration, the transformation temperatures of gypsum to anhydrite and hemihydrate will decrease. Thus, gypsum is stable at low temperatures in solutions of low H₂SO₄ concentrations and transforms to anhydrite at high temperatures and in concentrated H₂SO₄ solutions, while hemihydrate is always a metastable phase. Furthermore, the predicted results were compared with previous experimental studies to verify the accuracy of the model.

Thermochemical study for remediation of highly concentrated acid spill: Computational modeling and experimental validation

The release of concentrated acid solutions by chemical accidents is disastrous to our environmental integrity. Alkaline agents applied to remedy the acid spill catastrophe may lead to secondary damages such as vaporization or spread out of the fumes unless substantial amount of neutralization heat is properly controlled. Using a rigorous thermodynamic formalism proposed by Pitzer to account short-range ion interactions and various subsidiary reactions, we develop a systematic computational model enabling quantitative prediction of reaction heat and the temperature change over neutralization of strongly concentrated acid solutions. We apply this model to four acid solutions (HCl, HNO₃, H₂SO₄, and HF) of each 3 M-equivalent concentration with two neutralizing agents of calcium hydroxide (Ca(OH)₂) and sodium bicarbonate (NaHCO₃). Predicted reaction heat and temperature are remarkably consistent with the outcomes measured by our own experiments, showing a linear correlation factor R² greater than 0.98. We apply the model to extremely concentrated acid solutions as high as 50 wt% where an experimental approach is practically restricted. In contrast to the extremely exothermic Ca(OH)₂ agent, NaHCO₃ even lowers solution temperatures after neutralization reactions. Our model enables us to identify a promising neutralizer NaHCO₃ for effectively controlling concentrated acid spills and may be useful for establishment of proper strategy for other chemical accidents.

The Acidity of Atmospheric Particles and Clouds

Acidity, defined as pH, is a central component of aqueous chemistry. In the atmosphere, the acidity of condensed phases (aerosol particles, cloud water, and fog droplets) governs the phase partitioning of semi-volatile gases such as HNO3, NH3, HCl, and organic acids and bases as well as chemical reaction rates. It has implications for the atmospheric lifetime of pollutants, deposition, and human health. Despite its fundamental role in atmospheric processes, only recently has this field seen a growth in the number of studies on particle acidity. Even with this growth, many fine particle pH estimates must be based on thermodynamic model calculations since no operational techniques exist for direct measurements. Current information indicates acidic fine particles are ubiquitous, but observationally-constrained pH estimates are limited in spatial and temporal coverage. Clouds and fogs are also generally acidic, but to a lesser degree than particles, and have a range of pH that is quite sensitive to anthropogenic emissions of sulfur and nitrogen oxides, as well as ambient ammonia. Historical measurements indicate that cloud and fog droplet pH has changed in recent decades in response to controls on anthropogenic emissions, while the limited trend data for aerosol particles indicates acidity may be relatively constant due to the semi-volatile nature of the key acids and bases and buffering in particles. This paper reviews and synthesizes the current state of knowledge on the acidity of atmospheric condensed phases, specifically particles and cloud droplets. It includes recommendations for estimating acidity and pH, standard nomenclature, a synthesis of current pH estimates based on observations, and new model calculations on the local and global scale.

Measurements of Immersion Freezing and Heterogeneous Chemistry of Atmospherically Relevant Single Particles with Micro-Raman Spectroscopy

In the atmosphere, there are several different trajectories by which particles can nucleate ice; two of the major pathways are deposition and immersion freezing. Single particle depositional freezing has been widely studied with spectroscopic methods while immersion freezing has been predominantly studied either for particles within bulk aqueous solutions or using optical imaging of single particles. Of the few existing spectroscopic methods that monitor immersion freezing, there are limited opportunities for investigating the impact of heterogeneous chemistry on freezing. Herein, we describe a method that couples a confocal Raman spectrometer with an environmental cell to investigate single particle immersion freezing along with the capability to investigate in situ the impact of heterogeneous reactions with ozone and other trace gases on ice nucleation. This system, which has been rigorously calibrated (temperature and relative humidity) across a large dynamic range, is used to investigate low temperature water uptake and heterogeneous ice nucleation of atmospherically relevant single particles deposited on a substrate. The use of Raman spectroscopy provides important insights into the phase state and chemical composition of ice nuclei and, thus, insights into cloud formation.

Kinetics of Levoglucosenone Isomerization

We studied the acid-catalyzed isomerization of levoglucosenone (LGO) to 5-hydroxymethylfurfural (HMF) and developed a reaction kinetics model that describes the experimental data across a range of conditions (100-150 °C, 50-100 mm H₂ SO₄ , 50-150 mm LGO). LGO and its hydrated derivative exist in equilibrium under these reaction conditions. Thermal and catalytic degradation of HMF are the major sources of carbon loss. Within the range of conditions studied, higher temperatures and shorter reaction times favor the production of HMF. The yields of HMF and levulinic acid decrease monotonically as tetrahydrofuran is added to the aqueous solvent system, indicating that water plays a role in the LGO isomerization reaction. Initial-rate analyses show that HMF is produced solely from LGO rather than from the hydrated derivative of LGO. The results of this study are consistent with a mechanism for LGO isomerization that proceeds through hydration of the anhydro bridge, followed by ring rearrangement analogous to the isomerization of glucose to fructose.

Sulphur Kβ emission spectra reveal protonation states of aqueous sulfuric acid

In this paper we report an X-ray emission study of bulk aqueous sulfuric acid. Throughout the range of molarities from 1 M to 18 M the sulfur Kβ emission spectra from H2SO4 (aq) depend on the molar fractions and related deprotonation of H2SO4. We compare the experimental results with results from emission spectrum calculations based on atomic structures of single molecules and structures from ab initio molecular dynamics simulations. We show that the S Kβ emission spectrum is a sensitive probe of the protonation state of the acid molecules. Using non-negative matrix factorization we are able to extract the fractions of different protonation states in the spectra, and the results are in good agreement with the simulation for the higher part of the concentration range.

Direct Measurement of pH in Individual Particles via Raman Microspectroscopy and Variation in Acidity with Relative Humidity

Atmospheric aerosol acidity is an important characteristic of aqueous particles, which has been linked to the formation of secondary organic aerosol by catalyzing reactions of oxidized organic compounds that have partitioned to the particle phase. However, aerosol acidity is difficult to measure and traditionally estimated using indirect methods or assumptions based on composition. Ongoing disagreements between experiments and thermodynamic models of particle acidity necessitate improved fundamental understanding of pH and ion behavior in high ionic strength atmospheric particles. Herein, Raman microspectroscopy was used to determine the pH of individual particles (H2SO4+MgSO4) based on sulfate and bisulfate concentrations determined from νs(SO4(2-)) and νs(HSO4(-)), the acid dissociation constant, and activity coefficients from extended Debye-Hückel calculations. Shifts in pH and peak positions of νs(SO4(2-)) and νs(HSO4(-)) were observed as a function of relative humidity. These results indicate the potential for direct spectroscopic determination of pH in individual particles and the need to improve fundamental understanding of ion behavior in atmospheric particles.

Structure and ionization of sulfuric acid in water

High-sensitivity Raman spectra of aqueous H₂SO₄ indicate that the acid in water is H₄SO₅ dissociating to 3H₃O⁺ and 1HSO₅³⁻.

A water activity based model of heterogeneous ice nucleation kinetics for freezing of water and aqueous solution droplets

Immersion freezing of water and aqueous solutions by particles acting as ice nuclei (IN) is a common process of heterogeneous ice nucleation which occurs in many environments, especially in the atmosphere where it results in the glaciation of clouds. Here we experimentally show, using a variety of IN types suspended in various aqueous solutions, that immersion freezing temperatures and kinetics can be described solely by temperature, T, and solution water activity, a(w), which is the ratio of the vapour pressure of the solution and the saturation water vapour pressure under the same conditions and, in equilibrium, equivalent to relative humidity (RH). This allows the freezing point and corresponding heterogeneous ice nucleation rate coefficient, J(het), to be uniquely expressed by T and a(w), a result we term the a(w) based immersion freezing model (ABIFM). This method is independent of the nature of the solute and accounts for several varying parameters, including cooling rate and IN surface area, while providing a holistic description of immersion freezing and allowing prediction of freezing temperatures, J(het), frozen fractions, ice particle production rates and numbers. Our findings are based on experimental freezing data collected for various IN surface areas, A, and cooling rates, r, of droplets variously containing marine biogenic material, two soil humic acids, four mineral dusts, and one organic monolayer acting as IN. For all investigated IN types we demonstrate that droplet freezing temperatures increase as A increases. Similarly, droplet freezing temperatures increase as the cooling rate decreases. The log10(J(het)) values for the various IN types derived exclusively by Tand a(w), provide a complete description of the heterogeneous ice nucleation kinetics. Thus, the ABIFM can be applied over the entire range of T, RH, total particulate surface area, and cloud activation timescales typical of atmospheric conditions. Lastly, we demonstrate that ABIFM can be used to derive frozen fractions of droplets and ice particle production for atmospheric models of cirrus and mixed phase cloud conditions.

First and second deprotonation of H2SO4on wet hydroxylated (0001) α-quartz

We present anab initiomolecular dynamics study of deprotonation of sulfuric acid on wet quartz, a topic of atmospheric interest.

Optical trapping and Raman spectroscopy of solid particles

Stable levitation and spectroscopic interrogation of solid particles is achieved, over extended time periods, using a new optical trap design.

Disproportionation of bromous acid HOBrO by direct O-transfer and via anhydrides O(BrO)2 and BrO-BrO2. An ab initio study of the mechanism of a key step of the Belousov-Zhabotinsky oscillating reaction

The results are reported of an ab initio study of the thermochemistry and of the kinetics of the HOBrO disproportionation reaction 2HOBrO (2) ⇄ HOBr (1) + HBrO(3) (3), reaction ( R4' ), in gas phase (MP2(full)/6-311G*) and aqueous solution (SMD(MP2(full)/6-311G*)). The reaction energy of bromous acid disproportionation is discussed in the context of the coupled reaction system R2-R4 of the FKN mechanism of the Belousov-Zhabotinsky reaction and considering the acidities of HBr and HOBrO(2). The structures were determined of ten dimeric aggregates 4 of bromous acid, (HOBrO)(2), of eight mixed aggregates 5 formed between the products of disproportionation, (HOBr)(HOBrO(2)), and of four transition states structures 6 for disproportionation by direct O-transfer. It was found that the condensation of two HOBrO molecules provides facile access to bromous acid anhydride 7, O(BrO)(2). A discussion of the potential energy surface of Br(2)O(3) shows that O(BrO)(2) is prone to isomerization to the mixed anhydride 8, BrO-BrO(2), and to dissociation to 9, BrO, and 10, BrO(2), and their radical pair 11. Hence, three possible paths from O(BrO)(2) to the products of disproportionation, HOBr and HOBrO(2), are discussed: (1) hydrolysis of O(BrO)(2) along a path that differs from its formation, (2) isomerization of O(BrO)(2) to BrO-BrO(2) followed by hydrolysis, and (3) O(BrO)(2) dissociation to BrO and BrO(2) and their reactions with water. The results of the potential energy surface analysis show that the rate-limiting step in the disproportionation of HOBrO consists of the formation of the hydrate 12a of bromous acid anhydride 7 via transition state structure 14a. The computed activation free enthalpy ΔG(act)(SMD) = 13.6 kcal/mol for the process 2·2a → [14a](‡) → 12a corresponds to the reaction rate constant k(4) = 667.5 M(-1) s(-1) and is in very good agreement with experimental measurements. The potential energy surface analysis further shows that anhydride 7 is kinetically and thermodynamically unstable with regard to hydrolysis to HOBr and HOBrO(2) via transition state structure 14b. The transition state structure 14b is much more stable than 14a, and, hence, the formation of the "symmetrical anhydride" from bromous acid becomes an irreversible reaction for all practical purposes because 7 will instead be hydrolyzed as a "mixed anhydride" to afford HOBr and HOBrO(2). The mixed anhydride 8, BrO-BrO(2), does not play a significant role in bromous acid disproportionation.

Galena weathering under simulated calcareous soil conditions

Exploitation of polymetallic deposits from calcareous mining sites exposes galena and others sulfides to weathering factors. Galena weathering leads to the formation of lead phases (e.g., PbSO(4), PbCO(3)) with a higher bioaccessibility than galena, thus increasing the mobility and toxicity of lead. Despite the environmental impacts of these lead phases, the mechanisms of galena oxidation and the transformation of lead secondary phases, under neutral-alkaline carbonated conditions, have rarely been studied. In this work, an experimental approach, combining electrochemical and spectroscopic techniques, was developed to examine the interfacial processes involved in the galena weathering under simulated calcareous conditions. The results showed an initial oxidation stage with the formation of an anglesite-like phase leading to the partial mineral passivation. Under neutral-alkaline carbonated conditions, the stability of this phase was limited as it transformed into a cerussite-like one. Based on the surface characterization and the formation of secondary species, the weathering mechanisms of galena in calcareous soil and its environmental implications were suggested.

Freezing of water and aqueous NaCl droplets coated by organic monolayers as a function of surfactant properties and water activity

This study presents heterogeneous ice nucleation from water and aqueous NaCl droplets coated by 1-nonadecanol and 1-nonadecanoic acid monolayers as a function of water activity (a(w)) from 0.8 to 1 accompanied by measurements of the corresponding pressure-area isotherms and equilibrium spreading pressures. For water and aqueous NaCl solutions of ~0-20 wt % in concentration, 1-nonadecanol exhibits a condensed phase, whereas the phase of 1-nonadecanoic acid changes from an expanded to a condensed state with increasing NaCl content of the aqueous subphase. 1-Nonadecanol-coated aqueous droplets exhibit the highest median freezing temperatures that can be described by a shift in a(w) of the ice melting curve by 0.098 according to the a(w)-based ice nucleation approach. This freezing curve represents a heterogeneous ice nucleation rate coefficient (J(het)) of 0.85 ± 0.30 cm(-2) s(-1). The median freezing temperatures of 1-nonadecanoic acid-coated aqueous droplets decrease less with increasing NaCl content compared to the homogeneous freezing temperatures. This trend in freezing temperature is best described by a linear function in a(w) and not by the a(w)-based ice nucleation approach most likely due to an increased ice nucleation efficiency of 1-nonadecanoic acid governed by the monolayer state. This freezing curve represents J(het) = 0.46 ± 0.16 cm(-2) s(-1). Contact angles (α) for 1-nonadecanol- and 1-nonadecanoic acid-coated aqueous droplets increase as temperature decreases for each droplet composition, but absolute values depend on employed water diffusivity and the interfacial energies of the ice embryo. A parametrization of log[J(het)(Δa(w))] is presented which allows prediction of freezing temperatures and heterogeneous ice nucleation rate coefficients for water and aqueous NaCl droplets coated by 1-nonadecanol without knowledge of the droplet's composition and α.

Homogeneous Ice Nucleation in Water and Aqueous Solutions

This review provides an introduction to ice nucleation processes in supercooled water and aqueous solutions. Concepts for experimental techniques suitable to study homogeneous ice nucleation are addressed, in particular differential scanning calorimetry of inverse emulsions. Ice nucleation data from aqueous solutions have been analyzed using two approaches, and the interrelations between those are examined. It is argued that the ice nucleation process is driven entirely by thermodynamic quantities and how this can be understood in the context of three proposed theories for supercooled liquid water. Ice nucleation data for pure water droplets surrounded by a gas have been compiled and evaluated; within experimental uncertainty neither a volume dependent nucleation process nor a surface dependent nucleation process is convincingly supported by the analysis. Finally, open questions in the area of supercooled aqueous solutions and ice nucleation are discussed.