Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe barrier to ice nucleation in monatomic water
297
0
0.0
(
0
)
Authors
Santi Prestipino
Ask ChatGPT about the research
No Arabic abstract
Crystallization from a supercooled liquid initially proceeds via the formation of a small solid embryo (nucleus), which requires surmounting an activation barrier. This phenomenon is most easily studied by numerical simulation, using specialized biased-sampling techniques to overcome the limitations imposed by the rarity of nucleation events. Here, I focus on the barrier to homogeneous ice nucleation in supercooled water, as represented by the monatomic-water model, which in the bulk exhibits a complex interplay between different ice structures. I consider various protocols to identify solidlike particles on a computer, which perform well enough for the Lennard-Jones model, and compare their respective impact on the shape and height of the nucleation barrier. It turns out that the effect is stronger on the nucleus size than on the barrier height. As a by-product of the analysis, I determine the structure of the nucleation cluster, finding that the relative amount of ice phases in the cluster heavily depends on the method used for classifying solidlike particles. Moreover, the phase which is most favored during the earlier stages of crystallization may happen, depending on the nucleation coordinate adopted, to be different from the stable polymorph. Therefore, the quality of a reaction coordinate cannot be assessed simply on the basis of the barrier height obtained. I explain how this outcome is possible and why it just points out the shortcoming of collective variables appropriate to simple fluids in providing a robust method of particle classification for monatomic water.
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.
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.
Ice-water, water-vapor interfaces and ice surface are studied by molecular dynamics simulations with the SPC/E model of water molecules having the purpose to estimate the profiles of electrostatic potential across the interfaces. We have proposed a methodology for calculating the profiles of electrostatic potential based on a trial particle, which showed good agreement for the case of electrostatic potential profile of the water-vapor interface of TIP4P model calculated in another way. The measured profile of electrostatic potential for the pure ice-water interface decreases towards the liquid bulk region, which is in agreement with simulations of preferential direction of motion of Li$^{+}$ and F$^{-}$ solute ions at the liquid side of the ice-water interface. These results are discussed in connection with the Workman-Reynolds effect.
A new theoretical model for self dynamic response is developed using Vibration-Transit (V-T) theory, and is applied to liquid sodium at all wavevectors q from the hydrodynamic regime to the free particle limit. In this theory the zeroth-order Hamiltonian describes the vibrational motion in a single random valley harmonically extended to infinity. This Hamiltonian is tractable, is evaluated a priori for monatomic liquids, and the same Hamiltonian (the same set of eigenvalues and eigenvectors) is used for equilibrium and nonequlibrium theory. Here, for the self intermediate scattering function Fself(q,t) we find the vibrational contribution is in near perfect agreement with molecular dynamics (MD) through short and intermediate times, at all q. This is direct confirmation that normal mode vibrational correlations are present in the motion of the liquid state. The primary transit effect is diffusive motion of the vibrational equilibrium positions, as the liquid transits rapidly among random valleys. This motion is modeled as a standard random walk, and the resulting theoretical Fself(q,t) is in excellent agreement with MD results at all q and t. In the limit for q to infinity, the theory automatically exhibits the correct approach to the free-particle limit. Also in the limit for q to zero, the hydrodynamic limit emerges as well. In contrast to the benchmark theories of generalized hydrodynamics and mode coupling, the present theory is near a priori, while achieving modestly better accuracy. Therefore, in our view, it constitutes an improvement over the traditional theories.
Crystals and other condensed matter systems described by density waves often exhibit dislocations. Here we show, by considering the topology of the ground state manifolds (GSMs) of such systems, that dislocations in the density phase field always split into disclinations, and that the disclinations themselves are constrained to sit at particular points in the GSM. Consequently, the topology of the GSM forbids zero-energy dislocation glide, giving rise to a Peirels-Nabarro barrier.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
