No Arabic abstract
A novel liquid-liquid phase transition has been proposed and investigated in a wide variety of pure substances recently, including water, silica and silicon. From computer simulations using the Stillinger-Weber classical empirical potential, Sastry and Angell [1] demonstrated a first order liquid-liquid transition in supercooled silicon, subsequently supported by experimental and simulation studies. Here, we report evidence for a liquid-liquid critical end point at negative pressures, from computer simulations using the SW potential. Compressibilities exhibit a growing maximum upon lowering temperature below 1500 K and isotherms exhibit density discontinuities below 1120 K, at negative pressure. Below 1120 K, isotherms obtained from constant volume-temperature simulations exhibit non-monotonic, van der Waals-like behavior signaling a first order transition. We identify Tc ~ 1120 +/- 12 K, Pc -0.60 +/- 0.15 GPa as the critical temperature and pressure for the liquid-liquid critical point. The structure of the liquid changes dramatically upon decreasing the temperature and pressure. Diffusivities vary over 4 orders of magnitude, and exhibit anomalous pressure dependence near the critical point. A strong relationship between local geometry quantified by the coordination number, and diffusivity, is seen, suggesting that atomic mobility in both low and high density liquids can usefully be analyzed in terms of defects in the tetrahedral network structure. We have constructed the phase diagram of supercooled silicon. We identify the lines of compressibility, density extrema (maxima and minima) and the spinodal which reveal the interconnection between thermodynamic anomalies and the phase behaviour of the system as suggested in previous works [2-9]
We present a review on the study of metastable silicon, primarily focusing mainly on the aspects of liquid-liquid transition, critical point and phase behaviour, structural and dynamic properties of liquid phase as well as crystal nucleation. We begin with an extensive survey of the investigations of liquid silicon pursued over three decades, with salient experimental, theoretical and simulation results. Following which we present various scenarios put forward to rationalize the density and related anomalies often observed in water and other network forming liquids. After which we present the more recent investigations (both simulation and experimental works) of the phase behavior of Silicon. Since a significant part of metastable silicon work is on a classical empirical potential an important question to address is the reliability of this potential in describing the behavior of silicon. To provide a critical assessment of the applicability of classical simulation results to real silicon we present a comparison of the structural, dynamical, and thermodynamic quantities obtained from the SW potential with those from ab initio simulations and with available experimental data. We also discuss the sensitivity of the thermodynamic properties to model parameters.
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.
The study of liquid-liquid phase transition has attracted considerable attention. One interesting example of such phenomenon is phosphorus for which the existence a first-order phase transition between a low density insulating molecular phase and a conducting polymeric phase has been experimentally established. In this paper, we model this transition by an ab-initio quality molecular dynamics simulation and explore a large portion of the liquid section of the phase diagram. We draw the liquid-liquid coexistence curve and determine that it terminates into a second-order critical point. Close to the critical point, large coupled structure and electronic structure fluctuations are observed.
Development of reliable interatomic potentials is crucial for theoretical studies of relationship between chemical composition, structure and observable properties in glass-forming metallic alloys. Due to ambiguity of potential parametrization procedure, some crucial properties of the system, such as crystallization stability or symmetry of the ground state crystal phase, may not be correctly reproduced in computer simulations. Here we address this issue for ${rm Cu_{64.5}Zr_{35.5}}$ alloy described by two modifications of embedded atom model potential as well as by textit{ab initio} molecular dynamics. We observe that, at low supercooling, both models provide very similar liquid structure, which agrees with that obtained by textit{ab initio} simulations. Hoverer, deeply supercooled liquids demonstrate essentially different local structure and so different crystallization stability. The system, which demonstrate more pronounced icosahedral sort-range order, is more stable to crystallization that is in agreement with Frank hypothesis.
Water shows intriguing thermodynamic and dynamic anomalies in the supercooled liquid state. One possible explanation of the origin of these anomalies lies in the existence of a metastable liquid-liquid phase transition (LLPT) between two (high and low density) forms of water. While the anomalies are observed in experiments on bulk and confined water and by computer simulation studies of water-like models, the existence of a LLPT in water is still debated. Unambiguous experimental proof of the existence of a LLPT in bulk supercooled water is hampered by fast ice nucleation which is a precursor of the hypothesized LLPT. Moreover, the hypothesized LLPT, being metastable, in principle cannot exist in the thermodynamic limit (infinite size, infinite time). Therefore, computer simulations of water models are crucial for exploring the possibility of the metastable LLPT and the nature of the anomalies. In this work, we present new simulation results in the NVT ensemble for one of the most accurate classical molecular models of water, TIP4P/2005. To describe the computed properties and explore the possibility of a LLPT we have applied two-structure thermodynamics, viewing water as a non-ideal mixture of two interconvertible local structures (states). The results suggest the presence of a liquid-liquid critical point and a LLPT in this model for the simulated length and time scales. We have compared the behavior of TIP4P/2005 with other popular water-like models, namely mW and ST2, and with real water, all of which are well described by two-state thermodynamics. In view of the current debate involving different studies of TIP4P/2005, we discuss consequences of metastability and finite size in observing the liquid-liquid separation. We also address the relationship between the phenomenological order parameter of two-structure thermodynamics and the microscopic nature of the low-density structure.