Structural Order in Glassy Water

We investigate structural order in glassy water by performing classical molecular dynamics simulations using the extended simple point charge (SPC/E) model of water. We perform isochoric cooling simulations across the glass transition temperature at different cooling rates and densities. We quantify structural order by orientational and translational order metrics. Upon cooling the liquid into the glassy state, both the orientational order parameter $Q$ and translational order parameter $\tau$ increase. At T=0 K, the glasses fall on a line in the $Q$-$\tau$ plane or {\it order map}. The position of this line depends only on density and coincides with the location in the order map of the inherent structures (IS) sampled upon cooling. We evaluate the energy of the IS, $e_{IS}(T)$, and find that both order parameters for the IS are proportional to $e_{IS}$. We also study the structural order during the transformation of low-density amorphous ice (LDA) to high-density amorphous ice (HDA) upon isothermal compression and are able to identify distinct regions in the order map corresponding to these glasses. Comparison of the order parameters for LDA and HDA with those obtained upon isochoric cooling indicates major structural differences between glasses obtained by cooling and glasses obtained by compression. These structural differences are only weakly reflected in the pair correlation function. We also characterize the evolution of structural order upon isobaric annealing, leading at high pressure to very-high density amorphous ice (VHDA).Comment: submitte

Spinodal of supercooled polarizable water

We develop a series of molecular dynamics computer simulations of liquid water, performed with a polarizable potential model, to calculate the spinodal line and the curve of maximum density inside the metastable supercooled region. After analysing the structural properties,the liquid spinodal line is followed down to T=210 K. A monotonic decrease is found in the explored region. The curve of maximum density bends on approaching the spinodal line. These results, in agreement with similar studies on non polarizable models of water, are consistent with the existence of a second critical point for water.Comment: 8 pages, 5 figures, 2 tables. To be published in Phys. Re

Saddles in the energy landscape: extensivity and thermodynamic formalism

We formally extend the energy landscape approach for the thermodynamics of liquids to account for saddle points. By considering the extensive nature of macroscopic potential energies, we derive the scaling behavior of saddles with system size, as well as several approximations for the properties of low-order saddles (i.e., those with only a few unstable directions). We then cast the canonical partition function in a saddle-explicit form and develop, for the first time, a rigorous energy landscape approach capable of reproducing trends observed in simulations, in particular the temperature dependence of the energy and fractional order of sampled saddles.Comment: 4 pages, 1 figur

Energy landscape and rigidity

The effects of floppy modes in the thermodynamical properties of a system are studied. From thermodynamical arguments, we deduce that floppy modes are not at zero frequency and thus a modified Debye model is used to take into account this effect. The model predicts a deviation from the Debye law at low temperatures. Then, the connection between the topography of the energy landscape, the topology of the phase space and the rigidity of a glass is explored. As a result, we relate the number of constraints and floppy modes with the statistics of the landscape. We apply these ideas to a simple model for which we provide an approximate expression for the number of energy basins as a function of the rigidity. This allows to understand certains features of the glass transition, like the jump in the specific heat or the reversible window observed in chalcogenide glasses.Comment: 1 text+3 eps figure

Liquid-Liquid Phase Transitions for Soft-Core Attractive Potentials

Using event driven molecular dynamics simulations, we study a three dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point and a liquid-liquid critical point. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and give us insight on the mechanisms ruling the dependence of the two critical points on the potential's parameters. The soft core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.Comment: 29 pages, 15 figure

Crystallization Mechanism of Hard Sphere Glasses

In supercooled liquids, vitrification generally suppresses crystallization. Yet some glasses can still crystallize despite the arrest of diffusive motion. This ill-understood process may limit the stability of glasses, but its microscopic mechanism is not yet known. Here we present extensive computer simulations addressing the crystallization of monodisperse hard-sphere glasses at constant volume (as in a colloid experiment). Multiple crystalline patches appear without particles having to diffuse more than one diameter. As these patches grow, the mobility in neighbouring areas is enhanced, creating dynamic heterogeneity with positive feedback. The future crystallization pattern cannot be predicted from the coordinates alone: crystallization proceeds by a sequence of stochastic micro-nucleation events, correlated in space by emergent dynamic heterogeneity.Comment: 4 pages 4 figures Accepted for publication in Phys. Rev. Lett., April 201