Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThermodynamics and kinetics of crystallisation in deeply supercooled Stillinger-Weber silicon
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the kinetics of crystallization in deeply supercooled liquid silicon employing computer simulations and the Stillinger-Weber three body potential. The free energy barriers to crystallisation are computed using umbrella sampling Monte Carlo simulations, and for selected low temperature and zero pressure state points, using unconstrained molecular dynamics simulations to reconstruct the free energy from a mean first passage time formulation. We focus on state points that have been described in earlier work [Sastry and Angell, Nature Mater., 2, 739, 2003] as straddling a first order liquid-liquid phase transition (LLPT) between two metastable liquid states. It was argued subsequently [Ricci et al., Mol. Phys., 117, 3254, 2019] that the apparent phase transition is in fact due the loss of metastability of the liquid state with respect to the globally stable crystalline state. The presence or absence of a barrier to crystallization for these state points is therefore of importance to ascertain, with due attention to ambiguities that may arise from the choice of order parameters. We discuss our choice of order parameters and also our choice of methods to calculate the free energy at deep supercooling. We find a well-defined free energy barrier to crystallisation and demonstrate that both umbrella sampling and mean first passage time methods yield results that agree quantitatively. Our results thus provide strong evidence against the possibility that the liquids at state points close to the reported LLPT exhibit slow, spontaneous crystallisation, but they do not address the existence of a LLPT (or lack thereof). We also compute the free energy barriers to crystallisation at other state points over a broad range of temperatures and pressures, and discuss the effect of changes in the microscopic structure of the metastable liquid on the free energy barrier heights.
rate research
Read More
Using a recently developed thermodynamic integration method, we compute the precise values of the excess Gibbs free energy (G^e) of the high density liquid (HDL) phase with respect to the crystalline phase at different temperatures (T) in the supercooled region of the Stillinger-Weber (SW) silicon [F. H. Stillinger and T. A. Weber, Phys. Rev. B. 32, 5262 (1985)]. Based on the slope of G^e with respect to T, we find that the absolute entropy of the HDL phase increases as its enthalpy changes from the equilibrium value at T ge 1065 K to the value corresponding to a non-equilibrium state at 1060 K. We find that the volume distribution in the equilibrium HDL phases become progressively broader as the temperature is reduced to 1060 K, exhibiting van-der-Waals (VDW) loop in the pressure-volume curves. Our results provides insight into the thermodynamic cause of the transition from the HDL phase to the low density phases in SW silicon, observed in earlier studies near 1060 K at zero pressure.
The supercooled liquid silicon, modeled by Stillinger-Weber potential, shows anomalous increase in heat capacity $C_p$, with a maximum $C_p$ value close to 1060 K at zero pressure. We study equilibration and relaxation of the supercooled SW Si, in the temperature range of 1060 K--1070 K at zero pressure. We find that as the relaxation of the metastable supercooled liquid phase initiates, a straight line region (SLR) is formed in cumulative potential energy distributions. The configurational temperature corresponding to the SLR is close to 1060 K, which was earlier identified as the freezing temperature of 4-coordinated amorphous network. The SLR is found to be tangential to the distribution of the metastable liquid phase and thus influences the broadness of the distribution. As the bath temperature is reduced from 1070 K to 1060 K, the effective temperature approaches the bath temperature which results in broadening of the metastable phase distribution. This, in turn, causes an increase in overall fluctuations of potential energy and hence an increase of heat capacity. We also find that during initial stages of relaxation, 4-coordinated atoms form 6-membered rings with a chair--like structure and other structural units that indicate crystallization. Simultaneously a strong correlation is established between the number of chair-shaped 6-membered rings and the number of 4-coordinated atoms in the system. This shows that all properties related to 4-coordinated particles are highly correlated as the SLR is formed in potential energy distributions and this can be interpreted as a consequence of `freezing of amorphous network formed by 4-coordinated particles.
Below the melting temperature $T_m$ crystals are the stable phase of typical elemental or molecular systems. However, cooling down a liquid below $T_m$, crystallization is anything but inevitable. The liquid can be supercooled, eventually forming a glass below the glass transition temperature $T_g$. Despite their long lifetimes and the presence of strong barriers that produces an apparent stability, supercooled liquids and glasses remain intrinsically metastable state and thermodynamically unstable towards the crystal. Here we investigated the isothermal crystallization kinetics of the prototypical strong glassformer GeO$_2$ in the deep supercooled liquid at 1100 K, about half-way between $T_m$ and $T_g$. The crystallization process has been observed through time-resolved neutron diffraction for about three days. Data show a continuous reorganization of the amorphous structure towards the alpha-quartz phase with the final material composed by crystalline domains plunged into a low-density, residual amorphous matrix. A quantitative analysis of the diffraction patterns allows determining the time evolution of the relative fractions of crystal and amorphous, that was interpreted through an empirical model for the crystallization kinetics. This approach provides a very good description of the experimental data and identifies a predator-prey-like mechanism between crystal and amorphous, where the density variation acts as blocking barrier.
This article summarizes some of the open questions in the field of active matter that have emerged during Active20, a nine-week program held at the Kavli Institute for Theoretical Physics (KITP) in Spring 2020. The article does not provide a review of the field, but rather a personal view of the authors, informed by contributions of all participants, on new directions in active matter research. The topics highlighted include: the ubiquitous occurrence of spontaneous flows and active turbulence and the theoretical and experimental challenges associated with controlling and harnessing such flows; the role of motile topological defects in ordered states of active matter and their possible biological relevance; the emergence of non-reciprocal effective interactions and the role of chirality in active systems and their intriguing connections to non-Hermitian quantum mechanics; the progress towards a formulation of the thermodynamics of active systems thanks to the feedback between theory and experiments; the impact of the active matter framework on our understanding of the emergent mechanics of biological tissue. These seemingly diverse phenomena all stem from the defining property of active matter - assemblies of self-driven entities that individually break time-reversal symmetry and collectively organize in a rich variety of nonequilibrium states.
A theoretical treatment of deeply supercooled liquids is difficult because their properties emerge from spatial inhomogeneities that are self-induced, transient, and nanoscopic. I use computer simulations to analyse self-induced static and dynamic heterogeneity in equilibrium systems approaching the experimental glass transition. I characterise the broad sample-to-sample fluctuations of salient dynamic and thermodynamic properties in elementary mesoscopic systems. Findings regarding local lifetimes and distributions of dynamic heterogeneity are in excellent agreement with recent single molecule studies. Surprisingly broad thermodynamic fluctuations are also found, which correlate well with dynamics fluctuations, thus providing a local test of the thermodynamic origin of slow dynamics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
