Ivan Saika-Voivod

Fragile-to-strong transition and polyamorphism in the energy landscape of liquid silica.

Liquid silica is the archetypal glass former, and compounds based on silica are ubiquitous as natural and man-made amorphous materials. Liquid silica is also the extreme case of a ‘strong’ liquid, in that the variation of viscosity with temperature closely follows the Arrhenius law as the liquid is cooled toward its glass transition temperature. In contrast, most liquids are to some degree ‘fragile’, showing significantly faster increases in their viscosity as the glass transition temperature is approached. Recent studies have demonstrated the controlling influence of the potential energy hypersurface (or ‘energy landscape’) of the liquid on the transport properties near the glass transition. But the origin of strong liquid behaviour in terms of the energy landscape has not yet been resolved. Here we study the static and dynamic properties of liquid silica over a wide range of temperature and density using computer simulations. The results reveal a change in the energy landscape with decreasing temperature, which underlies a transition from a fragile liquid at high temperature to a strong liquid at low temperature. We also show that a specific heat anomaly is associated with this fragile-to-strong transition, and suggest that this anomaly is related to the polyamorphic behaviour of amorphous solid silica.

Computer simulations of liquid silica: equation of state and liquid-liquid phase transition.

We conduct extensive molecular dynamics computer simulations of two models for liquid silica [the model of Woodcock, Angell and Cheeseman, J. Phys. Chem. 65, 1565 (1976); and that of van Beest, Kramer, and van Santen, Phys. Rev. Lett. 64, 1955 (1990)] to determine their thermodynamic properties at low temperature T across a wide density range. We find for both models a wide range of states in which isochores of the potential energy U are a linear function of T(3/5), as recently proposed for simple liquids [Rosenfeld and P. Tarazona, Mol. Phys. 95, 141 (1998)]. We exploit this behavior to fit an accurate equation of state to our thermodynamic data. Extrapolation of this equation of state to low T predicts the occurrence of a liquid-liquid phase transition for both models. We conduct simulations in the region of the predicted phase transition, and confirm its existence by direct observation of phase separating droplets of atoms with distinct local density and coordination environments.

Density minimum and liquid?liquid phase transition

We present a high-resolution computer simulation study of the equation of state of ST2 water, evaluating the liquid-state properties at 2718 state points, and precisely locating the liquid–liquid critical point (LLCP) occurring in this model. We are thereby able to reveal the interconnected set of density anomalies, spinodal instabilities and response function extrema that occur in the vicinity of an LLCP for the case of a realistic, off-lattice model of a liquid with local tetrahedral order. In particular, we unambiguously identify a density minimum in the liquid state, define its relationship to other anomalies, and show that it arises due to the approach of the liquid structure to a defect-free random tetrahedral network of hydrogen bonds.

Model for reversible colloidal gelation

We report a numerical study, covering a wide range of packing fraction Phi and temperature T, for a system of particles interacting via a square well potential supplemented by an additional constraint on the maximum number n(max) of bonded interactions. We show that, when n(max)<6, the liquid-gas coexistence region shrinks, giving access to regions of low Phi where dynamics can be followed down to low T without an intervening phase separation. We characterize these arrested states at low densities (gel states) in terms of structure and dynamical slowing down, pointing out features which are very different from the standard glassy states observed at high Phi values.

Free energy and configurational entropy of liquid silica: Fragile-to-strong crossover and polyamorphism

Recent molecular dynamics (MD) simulations of liquid silica, using the BKS model [Phys. Rev. Lett. 64, 1955 (1990)]], have demonstrated that the liquid undergoes a dynamical crossover from super-Arrhenius, or fragile behavior, to Arrhenius, or strong behavior, as temperature T is decreased. From extensive MD simulations, we show that this fragile-to-strong crossover (FSC) can be connected to changes in the properties of the potential energy landscape, or surface (PES), of the liquid. To achieve this, we use thermodynamic integration to evaluate the absolute free energy of the liquid over a wide range of density and T. We use this free energy data, along with the concept of inherent structures of the PES, to evaluate the absolute configurational entropy S(c) of the liquid. We find that the temperature dependence of the diffusion coefficient and of S(c) are consistent with the prediction of Adam and Gibbs, including in the region where we observe the FSC to occur. We find that the FSC is related to a change in the properties of the PES explored by the liquid, specifically an inflection in the T dependence of the average inherent structure energy. In addition, we find that the high T behavior of S(c) suggests that the liquid entropy might approach zero at finite T, behavior associated with the so-called Kauzmann paradox. However, we find that the change in the PES that underlies the FSC is associated with a change in the T dependence of S(c) that elucidates how the Kauzmann paradox is avoided in this system. Finally, we also explore the relation of the observed PES changes to the recently discussed possibility that BKS silica exhibits a liquid-liquid phase transition, a behavior that has been proposed to underlie the observed polyamorphism of amorphous solid silica.

Gel to glass transition in simulation of a valence-limited colloidal system

We numerically study a simple model for thermoreversible colloidal gelation in which particles can form reversible bonds with a predefined maximum number of neighbors. We focus on three and four maximally coordinated particles, since in these two cases the low valency makes it possible to probe, in equilibrium, slow dynamics down to very low temperatures T. By studying a large region of T and packing fraction phi we are able to estimate both the location of the liquid-gas phase separation spinodal and the locus of dynamic arrest, where the system is trapped in a disordered nonergodic state. We find that there are two distinct arrest lines for the system: a glass line at high packing fraction, and a gel line at low phi and T. The former is rather vertical (phi controlled), while the latter is rather horizontal (T controlled) in the phi-T plane. Dynamics on approaching the glass line along isotherms exhibit a power-law dependence on phi, while dynamics along isochores follow an activated (Arrhenius) dependence. The gel has clearly distinct properties from those of both a repulsive and an attractive glass. A gel to glass crossover occurs in a fairly narrow range in phi along low-T isotherms, seen most strikingly in the behavior of the nonergodicity factor. Interestingly, we detect the presence of anomalous dynamics, such as subdiffusive behavior for the mean squared displacement and logarithmic decay for the density correlation functions in the region where the gel dynamics interferes with the glass dynamics.

Phase diagram of silica from computer simulation

We evaluate the phase diagram of the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)]], a model of silica widely used in molecular dynamics (MD) simulations. We conduct MD simulations of the liquid, and three crystals ( beta -quartz, coesite, and stishovite) over wide ranges of temperature and density, and evaluate the total Gibbs free energy of each phase. The phase boundaries are determined by the intersection of these free energy surfaces. Not unexpectedly for a classical pair potential, our results reveal quantitative discrepancies between the locations of the BKS and real silica phase boundaries. At the same time, we find that the topology of the real phase diagram is reproduced, confirming that the BKS model provides a satisfactory qualitative description of a silicalike material. We also compare the phase boundaries with the locations of liquid-state thermodynamic anomalies identified in previous studies of the BKS model.

Effect of bond lifetime on the dynamics of a short-range attractive colloidal system.

We perform molecular dynamics simulations of short-range attractive colloid particles modeled by a narrow (3% of the hard sphere diameter) square well potential of unit depth. We compare the dynamics of systems with the same thermodynamics but different bond lifetimes, by adding to the square well potential a thin barrier at the edge of the attractive well. For permanent bonds, the relaxation time tau diverges as the packing fraction phi approaches a threshold related to percolation, while for short-lived bonds, the phi dependence of tau is more typical of a glassy system. At intermediate bond lifetimes, the phi dependence of tau is driven by percolation at low phi , but then crosses over to glassy behavior at higher phi . We also study the wave vector dependence of the percolation dynamics.

Routes to colloidal gel formation

We discuss features of simple inter-particle potentials which are able to generate low-packing fraction arrested states, i.e. gels, in the absence of a macroscopic phase separation phenomenon. We suggest that the ratio between surface and bulk free energy is crucial in favoring ideal gel states. Two possible models for gels, one based on the competition of short range attraction and long range repulsions and the other on local constraints disfavoring packed local structures are discussed.

Energy landscape of a simple model for strong liquids

We calculate the statistical properties of the energy landscape of a minimal model for strong network-forming liquids. Dynamic and thermodynamic properties of this model can be computed with arbitrary precision even at low temperatures. A degenerate disordered ground state and logarithmic statistics for the local minima energy distribution are the landscape signatures of strong liquid behavior. Differences from fragile liquid properties are attributed to the presence of a discrete energy scale, provided by the particle bonds, and to the intrinsic degeneracy of topologically disordered networks.