Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNucleation in aqueous NaCl solutions shifts from 1-step to 2-step mechanism on crossing the spinodal
114
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work, we use large-scale molecular dynamics simulations coupled to free energy calculations to identify for the first time a limit of stability (spinodal) and a change in the nucleation mechanism in aqueous NaCl solutions. This is a system of considerable atmospheric, geological and technical significance. We find that the supersaturated metastable NaCl solution reaches its limit of stability at sufficiently high salt concentrations, as indicated by the composition dependence of the salt chemical potential, indicating the transition to a phase separation by spinodal decomposition. However, the metastability limit of the NaCl solution does not correspond to spinodal decomposition with respect to crystallization. We find that beyond this spinodal, a liquid/amorphous separation occurs in the aqueous solution, whereby the ions first form disordered clusters. We term these clusters as amorphous salt. We also identify a transition from one- to two-step crystallization mechanism driven by a spinodal. In particular, crystallization from aqueous NaCl solution beyond the spinodal is a two-step process, in which the ions first phase-separate into disordered amorphous salt clusters, followed by the crystallization of ions in the amorphous salt phase. In contrast, in the aqueous NaCl solution at concentrations lower than the spinodal, crystallization occurs via a one-step process, as the ions aggregate directly into crystalline nuclei. The change of mechanism with increasing supersaturation underscores the importance of an accurate determination of the driving force for phase separation. The study has broader implications on the mechanism for nucleation of crystals from solutions at high supersaturations.
rate research
Read More
We used molecular dynamics simulations and the path sampling technique known as forward flux sampling to study homogeneous nucleation of NaCl crystals from supersaturated aqueous solutions at 298 K and 1 bar. Nucleation rates were obtained for a range of salt concentrations for the Joung-Cheatham NaCl force field combined with the SPC/E water model. The calculated nucleation rates are significantly lower than available experimental measurements. The estimates for the nucleation rates in this work do not rely on classical nucleation theory, but the pathways observed in the simulations suggest that the nucleation process is better described by classical nucleation theory than an alternative interpretation based on Ostwalds step rule, in contrast to some prior simulations of related models. In addition to the size of NaCl nucleus, we find that the crystallinity of a nascent cluster plays an important role in the nucleation process. Nuclei with high crystallinity were found to have higher growth probability and longer lifetimes, possibly because they are less exposed to hydration water.
We report a numerical simulation of the rate of crystal nucleation of sodium chloride from its melt at moderate supercooling. In this regime nucleation is too slow to be studied with brute-force Molecular Dynamics simulations. The melting temperature of (Tosi-Fumi) NaCl is $sim 1060$K. We studied crystal nucleation at $T$=800K and 825K. We observe that the critical nucleus formed during the nucleation process has the crystal structure of bulk NaCl. Interestingly, the critical nucleus is clearly faceted: the nuclei have a cubical shape. We have computed the crystal-nucleation rate using two completely different approaches, one based on an estimate of the rate of diffusive crossing of the nucleation barrier, the other based on the Forward Flux Sampling and Transition Interface Sampling (FFS-TIS) methods. We find that the two methods yield the same result to within an order of magnitude. However, when we compare the extrapolated simulation data with the only available experimental results for NaCl nucleation, we observe a discrepancy of nearly 5 orders of magnitude. We discuss the possible causes for this discrepancy.
Thermal gradients induce concentration gradients in alkali halide solutions, and the salt migrates towards hot or cold regions depending on the average temperature of the solution. This effect has been interpreted using the heat of transport, which provides a route to rationalize thermophoretic phenomena. Early theories provide estimates of the heat of transport at infinite dilution. These values are used to interpret thermodiffusion (Soret) and thermoelectric (Seebeck) effects. However, accessing heats of transport of individual ions at finite concentration remains an outstanding question both theoretically and experimentally. Here we discuss a computational approach to calculate heats of transport of aqueous solutions at finite concentrations, and apply our method to study lithium chloride solutions at concentrations $>0.5$~M. The heats of transport are significantly different for Li$^+$ and Cl$^-$ ions, unlike what is expected at infinite dilution. We find theoretical evidence for the existence of minima in the Soret coefficient of LiCl, where the magnitude of the heat of transport is maximized. The Seebeck coefficient obtained from the ionic heats of transport varies significantly with temperature and concentration. We identify thermodynamic conditions leading to a maximization of the thermoelectric response of aqueous solutions.
We consider the thermodynamic behavior of local fluctuations occurring in a stable or metastable bulk phase. For a system with three or more phases, a simple analysis based on classical nucleation theory predicts that small fluctuations always resemble the phase having the lowest surface tension with the surrounding bulk phase, regardless of the relative chemical potentials of the phases. We also identify the conditions at which a fluctuation may convert to a different phase as its size increases, referred to here as a fluctuation phase transition (FPT). We demonstrate these phenonena in simulations of a two dimensional lattice model by evaluating the free energy surface that describes the thermodynamic properties of a fluctuation as a function of its size and phase composition. We show that a FPT can occur in the fluctuations of either a stable or metastable bulk phase and that the transition is first-order. We also find that the FPT is bracketed by well-defined spinodals, which place limits on the size of fluctuations of distinct phases. Furthermore, when the FPT occurs in a metastable bulk phase, we show that the superposition of the FPT on the nucleation process results in two-step nucleation (TSN). We identify distinct regimes of TSN based on the nucleation pathway in the free energy surface, and correlate these regimes to the phase diagram of the bulk system. Our results clarify the origin of TSN, and elucidate a wide variety of phenomena associated with TSN, including the Ostwald step rule.
It has long been assumed the Earths solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earths core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
