The nucleation of ice from aqueous solutions is a process essential to myriad environmental and industrial processes, but the physical factors affecting the capacity of different solutes to depress the homogeneous nucleation temperature of ice are yet poorly understood. In this work, we demonstrate that for many binary aqueous solutions of non-ionic solutes, this depression is dominated by the entropy of the liquid phase. Employing the classic Turnbull interpretation of the interfacial free energy γ∼TSliquid−Ssolid and estimating solution entropies with a Flory-style modification of the ideal entropy of mixing that accounts for solute size effects, we demonstrate that mixing entropy alone predicts experimental homogeneous nucleation temperatures across a wide variety of non-ionic solutions. We anticipate that this physical insight will not only enhance a fundamental understanding of homogeneous nucleation processes across fields but also open new avenues to the rational design of aqueous solutions for desired nucleation behaviors.
more »
« less
A size-dependent ideal solution model for liquid–solid phase equilibria prediction in aqueous organic solutions
Predictive synthesis of aqueous organic solutions with desired liquid–solid phase equilibria could drive progress in industrial chemistry, cryopreservation, and beyond, but is limited by the predictive power of current solution thermodynamics models. In particular, few analytical models enable accurate liquidus and eutectic prediction based only on bulk thermodynamic properties of the pure components, requiring instead either direct measurement or costly simulation of solution properties. In this work, we demonstrate that a simple modification to the canonical ideal solution theory accounting for the entropic effects of dissimilar molecule sizes can transform its predictive power. Incorporating a Flory-style entropy of mixing term that includes both the mole and volume fractions of each component, we derive size-dependent equations for the ideal chemical potential and liquidus temperature, and use them to predict the binary phase diagrams of water and 10 organic solutes of varying sizes. We show that size-dependent prediction outperforms the ideal model in all cases, reducing average error in the predicted liquidus temperature by 59% (to 5.6 K), eutectic temperature by 45% (to 9.7 K), and eutectic composition by 43% (to 4.7 mol%), as compared to experimental data. Furthermore, by retaining the ideal assumption that the enthalpy of mixing is zero, we demonstrate that, for aqueous organic solutions, much of the deviation from ideality that is typically attributed to molecular interactions may in fact be explained by simple entropic size effects. These results suggest an underappreciated dominance of mixing entropy in these solutions, and provide a simple approach to predicting their phase equilibria.
more »
« less
- Award ID(s):
- 1941543
- PAR ID:
- 10665219
- Publisher / Repository:
- PNAS
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 122
- Issue:
- 8
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Entropy-stabilized oxides are single-phase, multicomponent oxides that are stabilized by a large entropy of mixing, ΔS, overcoming a positive enthalpy. Due to the −TΔS term in the Gibbs' free energy, G, it can be hypothesized that entropy-stabilized oxides demonstrate a robust thermal stability. Here, we investigate the high temperature stability (1300–1700 °C) of the prototypical entropy-stabilized rocksalt oxide (MgCoNiCuZn)0.2O in air. We find that at temperatures >1300 °C, the material gradually loses Cu and Zn with increasing temperature. Cu is lost through a selective melting as a Cu-rich liquid phase is formed. Zn is sublimed from the rocksalt phase at approximately similar temperatures to those corresponding to the Cu loss, significantly below both the melting temperature of ZnO and its solubility limit in a rocksalt phase. The elemental loss progressively reduces the entropy of mixing and results in a multiphase solid upon quenching to room temperature. We posit that the high-temperature solubility of Cu and Zn is correlated providing further evidence for entropic stabilization over general solubility arguments.more » « less
-
Rhabdophane is a hydrous phosphate that commonly replaces monazite as a weathering product in critical mineral deposits during the alteration of rare earth elements (REE) bearing carbonatites and alkaline igneous complexes. It is an important host to the light (L)REE (i.e., La to Gd) but the stability and structure of binary solid solutions between the Ce and the other LREE endmembers have not yet been determined experimentally. Here we present room temperature calorimetric experiments that were used to measure the enthalpy of precipitation of rhabdophane (Ce1−xREExPO4·nH2O; REE = La, Pr, Nd, Sm, Eu, and Gd). The solids were characterized using X-ray diffraction, scanning electron microscopy, Raman spectroscopy, and the role of water in the rhabdophane structure was further determined using thermogravimetric analysis coupled with differential scanning calorimetry. The calorimetric experiments indicate a non-ideal behavior for all of the binary solid solutions investigated with an excess enthalpy of mixing (ΔHex) described by a 2- to 3-term Guggenheim parameters equation. The solid solutions were categorized into three groups: (1) binary Ce-La and Ce-Pr which display positive ΔHex values with a slight asymmetry; (2) binary Ce-Nd and Ce-Sm which display negative ΔHex values with a nearly symmetric shape; (3) Ce-Eu and Ce-Gd which display both negative and positive ΔHex values with nearly symmetric shape. The excess Gibbs energy (ΔGex) of the solid solutions was further investigated using a thermodynamic analysis approach of aqueous-solid solution equilibria and the optimization programs GEMS and GEMSFITS. The resulting ΔGex values combined with the calorimetric ΔHex values indicate that there is likely an excess entropy contribution implying important short-range structural modifications in the solid solutions dependent on the deviation of the REE ionic radii from the size of Ce3+. These observations corroborate with the trends in the Raman v1 stretching bands of the PO4-site. The excess molar volumes determined from X-ray diffraction analysis further indicate an overall asymmetric behavior in all of the studied binary solid solutions, which becomes increasingly important from La to Gd. The pronounced short-range order–disorder occurring in groups 2 and 3 solid solutions mimics some of the behavior observed from previous studies in anhydrous monazite solid solutions. This study highlights the potential to use the chemistry and the structural modifications of rhabdophane as potential indicators of formation conditions in geologic systems and permits improving our modeling capabilities of REE partitioning in critical minerals systems.more » « less
-
We introduce a notion called entropic independence that is an entropic analog of spectral notions of high-dimensional expansion. Informally, entropic independence of a background distribution $$\mu$$ on $$k$$-sized subsets of a ground set of elements says that for any (possibly randomly chosen) set $$S$$, the relative entropy of a single element of $$S$$ drawn uniformly at random carries at most $O(1/k)$ fraction of the relative entropy of $$S$$. Entropic independence is the analog of the notion of spectral independence, if one replaces variance by entropy. We use entropic independence to derive tight mixing time bounds, overcoming the lossy nature of spectral analysis of Markov chains on exponential-sized state spaces. In our main technical result, we show a general way of deriving entropy contraction, a.k.a. modified log-Sobolev inequalities, for down-up random walks from spectral notions. We show that spectral independence of a distribution under arbitrary external fields automatically implies entropic independence. We furthermore extend our theory to the case where spectral independence does not hold under arbitrary external fields. To do this, we introduce a framework for obtaining tight mixing time bounds for Markov chains based on what we call restricted modified log-Sobolev inequalities, which guarantee entropy contraction not for all distributions, but for those in a sufficiently large neighborhood of the stationary distribution. To derive our results, we relate entropic independence to properties of polynomials: $$\mu$$ is entropically independent exactly when a transformed version of the generating polynomial of $$\mu$$ is upper bounded by its linear tangent; this property is implied by concavity of the said transformation, which was shown by prior work to be locally equivalent to spectral independence. We apply our results to obtain (1) tight modified log-Sobolev inequalities and mixing times for multi-step down-up walks on fractionally log-concave distributions, (2) the tight mixing time of $$O(n\log n)$$ for Glauber dynamics on Ising models whose interaction matrix has eigenspectrum lying within an interval of length smaller than $$1$$, improving upon the prior quadratic dependence on $$n$$, and (3) nearly-linear time $$\widetilde O_{\delta}(n)$$ samplers for the hardcore and Ising models on $$n$$-node graphs that have $$\delta$$-relative gap to the tree-uniqueness threshold. In the last application, our bound on the running time does not depend on the maximum degree $$\Delta$$ of the graph, and is therefore optimal even for high-degree graphs, and in fact, is sublinear in the size of the graph for high-degree graphs.more » « less
-
null (Ed.)The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N, the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N, and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.more » « less
An official website of the United States government

