Search for a liquid-liquid critical point in models of silica

Crystallization and the liquid-liquid critical point in nonbonded modified-WAC models

For decades, it has been known that Liquid-Liquid Critical Points (LLCPs) can exist in one-component liquids, yet a comprehensive understanding of the conditions under which they arise remains elusive. To better comprehend the possible interplay between the LLCP and the crystalline phase, we conduct molecular dynamics simulations using the nonbonded family of modified-WAC (mWAC) models, which are known to exhibit a LLCP for certain parameter values. By comparing different versions of the mWAC model-those featuring a LLCP and those lacking one-we identify several key differences between the models relating to crystallization. Those models that do have a LLCP are found to have multiple stable crystalline phases, one of them being a solid-state ionic conductor similar to superionic ice. Moreover, we find that for models that do not have a LLCP, the liquid becomes a glass at a larger range of temperatures, possibly preventing the occurrence of a LLCP. Further studies are required to determine if these results are general or model-specific.

Evolution of Short-Range Order of Amorphous GeTe Upon Structural Relaxation Obtained by TEM Diffractometry and RMC Methods

Glasses frequently reveal structural relaxation that leads to changes in their physical properties including enthalpy, specific volume, and resistivity. Analyzing the short-range order (SRO) obtained from electron diffraction by transmission electron microscopy (TEM) in combination with Reverse-Monte-Carlo (RMC) simulations is shown to provide information on the atomic arrangement. The technique elaborated here features several benefits including reliability, accessibility, and allows for obtaining detailed structural data quickly. This is demonstrated with a detailed view of the structural changes in the as-deposited amorphous phase change material (PCM) GeTe. The data show a significant increase in the average bond angle upon thermal treatment. At the same time the fraction of tetrahedrally coordinated Ge atoms decreases due to an increase in octahedrally distorted and pyramidal motifs. This finding provides further evidence for the atomic processes that govern structural relaxation in amorphous GeTe and other PCMs. A thorough literature review finally unveils possible origins of the large discrepancies reported on the structure of amorphous GeTe.

Thermodynamic anomalies, polyamorphism and all that

The appearance and evolution of thermodynamics anomalies, and related properties, are studied for two classes of system, modelling those dominated by covalent and ionic interactions, respectively. Such anomalies are most familiar in the density but are also present in other thermodynamics variables such as the compressibility and heat capacity. By systematically varying key model parameters the emergence and evolution of these anomalies can be tracked across the phase space. The interaction of the anomalies can often be rationalized by thermodynamics 'rules'. The emergence of these anomalies may also be correlated with the appearance of polyamorphism, the existence of multiple amorphous states which differ in density and entropy. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

Thermodynamic modeling of fluid polyamorphism in hydrogen at extreme conditions

Fluid polyamorphism, the existence of multiple amorphous fluid states in a single-component system, has been observed or predicted in a variety of substances. A remarkable example of this phenomenon is the fluid–fluid phase transition (FFPT) in high-pressure hydrogen between insulating and conducting high-density fluids. This transition is induced by the reversible dimerization/dissociation of the molecular and atomistic states of hydrogen. In this work, we present the first attempt to thermodynamically model the FFPT in hydrogen at extreme conditions. Our predictions for the phase coexistence and the reaction equilibrium of the two alternative forms of fluid hydrogen are based on experimental data and supported by the results of simulations. Remarkably, we find that the law of corresponding states can be utilized to construct a unified equation of state combining the available computational results for different models of hydrogen and the experimental data.

Modeling fluid polyamorphism through a maximum-valence approach

We suggest a simple model to describe polyamorphism in single-component fluids using a maximum-valence approach. The model contains three types of interactions: (i) Atoms attract each other by van der Waals forces that generate a liquid-gas transition at low pressures, (ii) atoms may form covalent bonds that induce association, and (iii) atoms with maximal valence attract or repel each other stronger than other atoms, thus generating liquid-liquid separation. As an example, we qualitatively compare this model with the behavior of liquid sulfur and show that condition (iii) generates a liquid-liquid phase transition in addition to the liquid-gas phase transition.

Phase transitions affected by natural and forceful molecular interconversion

If a binary liquid mixture, composed of two alternative species with equal amounts, is quenched from a high temperature to a low temperature, below the critical point of demixing, then the mixture will phase separate through a process known as spinodal decomposition. However, if the two alternative species are allowed to interconvert, either naturally (e.g., the equilibrium interconversion of enantiomers) or forcefully (e.g., via an external source of energy or matter), then the process of phase separation may drastically change. In this case, depending on the nature of interconversion, two phenomena could be observed: either phase amplification, the growth of one phase at the expense of another stable phase, or microphase separation, the formation of nongrowing (steady-state) microphase domains. In this work, we phenomenologically generalize the Cahn–Hilliard theory of spinodal decomposition to include the molecular interconversion of species and describe the physical properties of systems undergoing either phase amplification or microphase separation. We apply the developed phenomenology to accurately describe the simulation results of three atomistic models that demonstrate phase amplification and/or microphase separation. We also discuss the application of our approach to phase transitions in polyamorphic liquids. Finally, we describe the effects of fluctuations of the order parameter in the critical region on phase amplification and microphase separation.

Dynamical crossover and its connection to the Widom line in supercooled TIP4P/Ice water

We perform molecular dynamics simulations with the TIP4P/Ice water model to characterize the relationship between dynamics and thermodynamics of liquid water in the supercooled region. We calculate the relevant properties of the phase diagram, and we find that TIP4P/Ice presents a retracing line of density maxima, similar to what was previously found for atomistic water models and models of other tetrahedral liquids. For this model, a liquid-liquid critical point between a high-density liquid and a low-density liquid was recently found. We compute the lines of the maxima of isothermal compressibility and the minima of the coefficient of thermal expansion in the one phase region, and we show that these lines point to the liquid-liquid critical point while collapsing on the Widom line. This line is the line of the maxima of correlation length that emanates from a second order critical point in the one phase region. Supercooled water was found to follow mode coupling theory and to undergo a transition from a fragile to a strong behavior right at the crossing of the Widom line. We find here that this phenomenology also happens for TIP4P/Ice. Our results appear, therefore, to be a general characteristic of supercooled water, which does not depend on the interaction potential used, and they reinforce the idea that the dynamical crossover from a region where the relaxation mechanism is dominated by cage relaxation to a region where cages are frozen and hopping dominates is correlated in water to a phase transition between a high-density liquid and a low-density liquid.

The anomalies and criticality of liquid water

The origin of water's anomalies has been a matter of long-standing debate. A two-state model, dating back to Röntgen, relies on the dynamical coexistence of two types of local structures-locally favored tetrahedral structure (LFTS) and disordered normal-liquid structure (DNLS)-in liquid water. Phenomenologically, this model not only explains water's thermodynamic anomalies but also can rationalize the existence of a liquid-liquid critical point (LLCP) if there is a cooperative formation of LFTS. We recently found direct evidence for the coexistence of LFTS and DNLS in the experimental structure factor of liquid water. However, the existence of the LLCP and its impact on water's properties has remained elusive, leaving the origin of water's anomalies unclear. Here we propose a unique strategy to locate the LLCP of liquid water. First, we make a comprehensive analysis of a large set of experimental structural, thermodynamic, and dynamic data based on our hierarchical two-state model. This model predicts that the two thermodynamic and dynamical fluctuation maxima lines should cross at the LLCP if it exists, which we confirm by hundred-microsecond simulations for model waters. Based on recent experimental results of the compressibility and diffusivity measurements in the no man's land, we reveal that the two lines cross around 184 K and 173 MPa for real water, suggesting the presence of the LLCP around there. Nevertheless, we find that the criticality is almost negligible in the experimentally accessible region of liquid water because it is too far from the LLCP. Our findings would provide a clue to settle the long-standing debate.

A mixed radial, angular, three-body distribution function as a tool for local structure characterization: Application to single-component structures

A mixed radial, angular three-body distribution function g₃(r_BC, θ_ABC) is introduced, which allows the local atomic order to be more easily characterized in a single graph than with conventional correlation functions. It can be defined to be proportional to the probability of finding an atom C at a distance r_BC from atom B while making an angle θ_ABC with atoms A and B, under the condition that atom A is the nearest neighbor of B. As such, our correlation function contains, for example, the likelihood of angles formed between the nearest and the next-nearest-neighbor bonds. To demonstrate its use and usefulness, a visual library for many one-component crystals is produced first and then employed to characterize the local order in a diverse body of elemental condensed-matter systems. Case studies include the analysis of a grain boundary, several liquids (argon, copper, and antimony), and polyamorphism in crystalline and amorphous silicon including that obtained in a tribological interface.

Effects of the bond polarity on the structural and dynamical properties of silica-like liquids

Silica is a network-forming liquid that shares many properties with water due to its tetrahedral structure. It undergoes a transition from a fragile to a strong liquid as the temperature is decreased, which is accompanied by a structural change to lower density and higher tetrahedral order. In order to disentangle the effects of Coulomb and van der Waals interactions on the structure and dynamics of liquid silica, we modify the bond polarity by changing the partial charges assigned to each atom. Using molecular dynamics simulations, we show that density, tetrahedral order, and structural relaxation times decrease when reducing bond polarity. Moreover, we find that the density maximum and the fragile-to-strong transition move to lower temperatures until they eventually vanish when the partial charges are decreased below approximately 75% of their regular value. Irrespective of whether strong or fragile behavior exists, structural relaxation is governed by hopping motion at sufficiently low temperatures. As long as there is a strong regime, the energy barrier associated with strong dynamics decreases with decreasing partial charges, but the dependence on the bond polarity differs from that of the activation energy in the Arrhenius regime at high temperatures. We show that the fragile-to-strong transition is associated with structural changes occurring between the first and second coordination shells that lead to a decrease in density and an increase in tetrahedral order. In particular, independent of the value of the partial charges, the distribution of the local structures is the same at this dynamic crossover, but we find no evidence that the effect occurs upon crossing the Widom line. In the fragile regime at intermediate temperatures, the relaxation times are well described by a previously proposed model which decomposes the apparent activation energy into a constant single-particle contribution and a temperature-dependent collective contribution. However, our results for silica-like melts do not obey several common relations of the model parameters reported for molecular glass formers.

Waterlike anomalies in a two-dimensional core-softened potential

We investigate the structural, thermodynamic, and dynamic behavior of a two-dimensional (2D) core-corona system using Langevin dynamics simulations. The particles are modeled by employing a core-softened potential which exhibits waterlike anomalies in three dimensions. In previous studies in a quasi-2D system a new region in the pressure versus temperature phase diagram of structural anomalies was observed. Here we show that for the two-dimensional case two regions in the pressure versus temperature phase diagram with structural, density, and diffusion anomalies are observed. Our findings indicate that, while the anomalous region at lower densities is due the competition between the two length scales in the potential at higher densities, the anomalous region is related to the reentrance of the melting line.

Fluctuations near the liquid–liquid transition in a model of silica

Molecular dynamics simulations reveal anomalous small-angle scattering and liquid–liquid phase separation in an ionic model of silica.

Water: A Tale of Two Liquids

Water is the most abundant liquid on earth and also the substance with the largest number of anomalies in its properties. It is a prerequisite for life and as such a most important subject of current research in chemical physics and physical chemistry. In spite of its simplicity as a liquid, it has an enormously rich phase diagram where different types of ices, amorphous phases, and anomalies disclose a path that points to unique thermodynamics of its supercooled liquid state that still hides many unraveled secrets. In this review we describe the behavior of water in the regime from ambient conditions to the deeply supercooled region. The review describes simulations and experiments on this anomalous liquid. Several scenarios have been proposed to explain the anomalous properties that become strongly enhanced in the supercooled region. Among those, the second critical-point scenario has been investigated extensively, and at present most experimental evidence point to this scenario. Starting from very low temperatures, a coexistence line between a high-density amorphous phase and a low-density amorphous phase would continue in a coexistence line between a high-density and a low-density liquid phase terminating in a liquid–liquid critical point, LLCP. On approaching this LLCP from the one-phase region, a crossover in thermodynamics and dynamics can be found. This is discussed based on a picture of a temperature-dependent balance between a high-density liquid and a low-density liquid favored by, respectively, entropy and enthalpy, leading to a consistent picture of the thermodynamics of bulk water. Ice nucleation is also discussed, since this is what severely impedes experimental investigation of the vicinity of the proposed LLCP. Experimental investigation of stretched water, i.e., water at negative pressure, gives access to a different regime of the complex water diagram. Different ways to inhibit crystallization through confinement and aqueous solutions are discussed through results from experiments and simulations using the most sophisticated and advanced techniques. These findings represent tiles of a global picture that still needs to be completed. Some of the possible experimental lines of research that are essential to complete this picture are explored.

Universal scaling in the aging of the strong glass former SiO2

We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely, the dynamical heterogeneities. We perform molecular dynamics simulations of SiO2 with van Beest-Kramer-van Santen interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time tw measured from the instant of the quench. We find that both the aging behavior of the dynamic susceptibility χ4 and the aging behavior of the probability distribution P(fs,r) of the local incoherent intermediate scattering function fs,r can be described by simple scaling forms in terms of the global incoherent intermediate scattering function C. The scaling forms are the same that have been found to describe the aging of several fragile glass formers and that, in the case of P(fs,r), have been also predicted theoretically. A thorough study of the length scales involved highlights the importance of intermediate length scales. We also analyze directly the scaling dependence on particle type and on wavevector q and find that both the average and the fluctuations of the slow aging dynamics are controlled by a unique aging clock, which is not only independent of the wavevector q, but is also the same for O and Si atoms.

Tunable Liquid-Liquid Critical Point in an Ionic Model of Silica

Recently, it was shown that the Woodcock-Angell-Cheeseman model for liquid silica [L. V. Woodcock, C. A. Angell, and P. Cheeseman, J. Chem. Phys. 65, 1565 (1976)] is remarkably close to having a liquid-liquid critical point (LLCP). We demonstrate that increasing the ion charge separates the global maxima of the response functions, while reducing the charge smoothly merges them into a LLCP, a phenomenon that might be experimentally observable with charged colloids. An analysis of the Si and O coordination numbers suggests that a sufficiently low Si/O coordination number ratio is needed to attain a LLCP.