Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTheoretical prediction of the homogeneous ice nucleation rate: disentangling thermodynamics and kinetics
150
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Estimating the homogeneous ice nucleation rate from undercooled liquid water is at the same time crucial for understanding many important physical phenomena and technological applications, and challenging for both experiments and theory. From a theoretical point of view, difficulties arise due to the long time scales required, as well as the numerous nucleation pathways involved to form ice nuclei with different stacking disorders. We computed the homogeneous ice nucleation rate at a physically relevant undercooling for a single-site water model, taking into account the diffuse nature of ice-water interfaces, stacking disorders in ice nuclei, and the addition rate of particles to the critical nucleus.We disentangled and investigated the relative importance of all the terms, including interfacial free energy, entropic contributions and the kinetic prefactor, that contribute to the overall nucleation rate.There has been a long-standing discrepancy for the predicted homogeneous ice nucleation rates, and our estimate is faster by 9 orders of magnitude compared with previous literature values. Breaking down the problem into segments and considering each term carefully can help us understand where the discrepancy may come from and how to systematically improve the existing computational methods.
rate research
Read More
Ice nucleation is a process of great relevance in physics, chemistry, technology and environmental sciences, much theoretical and experimental efforts have been devoted to its understanding, but still it remains a topic of intense research. We shed light on this phenomenon by performing atomistic based simulations. Using metadynamics and a carefully designed set of collective variables, reversible transitions between water and ice are able to be simulated. We find that water freezes into a stacking disordered structure with the all-atom TIP4P/Ice model, and the features of the critical nucleus of nucleation at the microscopic level are revealed. Our results are in agreement with recent experimental and other theoretical works and confirm that nucleation is preceded by a large increase in tetrahedrally coordinated water molecules.
We develop a theory in order to describe the effect of relaxation in a condensed medium upon the quantum decay of a metastable liquid near the spinodal at low temperatures. We find that both the regime and the rate of quantum nucleation strongly depend on the relaxation time and its temperature behavior. The quantum nucleation rate slows down with the decrease of the relaxation time. We also discuss the low temperature experiments on cavitation in normal $^3$He and superfluid $^4$He at negative pressures. It is the sharp distinctions in the high frequency sound mode and in the temperature behavior of the relaxation time that make the quantum cavitation kinetics in $^3$He and $^4$He completely different in kind.
The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.
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.
Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical $(N, V, E)$ ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy $E$ and number of atoms $N$, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time $tau_{rm w}$ before melting occurs. We find that the probability distribution of $tau_{rm w}$ is roughly exponential, and that the mean value $<tau_{rm w} >$ depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean $<tau_{rm w}>$ also depends strongly on system size in a way that we have quantified. For very small systems of $sim 100$ atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials, and to suggest ways of correcting for the errors of the method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
