Polyamorphism of ice at low temperatures from constant-pressure simulations

We report results of MD simulations of amorphous ice in the pressure range 0 - 22.5 kbar. The high-density amorphous ice (HDA) prepared by compression of Ih ice at T = 80 K is annealed to T = 170 K at intermediate pressures in order to generate relaxed states. We confirm the existence of recently observed phenomena, the very high-density amorphous ice and a continuum of HDA forms. We suggest that both phenomena have their origin in the evolution of the network topology of the annealed HDA phase with decreasing volume, resulting at low temperatures in the metastability of a range of densities.Comment: 11 pages, 5 postscript figures. To be published in Physical Review Letter

A Potential Energy Landscape Study of the Amorphous-Amorphous Transformation in H2_2O

We study the potential energy landscape explored during a compression-decompression cycle for the SPC/E (extended simple point charge) model of water. During the cycle, the system changes from low density amorphous ice (LDA) to high density amorphous ice (HDA). After the cycle, the system does not return to the same region of the landscape, supporting the interesting possibility that more than one significantly different configuration corresponds to LDA. We find that the regions of the landscape explored during this transition have properties remarkably different from those explored in thermal equilibrium in the liquid phase

Comment on "First order amorphous-amorphous transformation in silica"

In a recent letter (Phys. Rev. Lett. 84, 4629 (2000)), Lacks presents evidence of a first order amorphous-amorphous transition in silica at T=0. He calculates the free energy along a path of compression and successive decompression of a sample of 108 SiO2 units. The free energy of the two branches cross each other, and this is interpreted as evidence of a first order transition. We show that this conclusion does not follow from the shown data, since qualitatively the same phenomenology is obtained in a model where a first order transition does not exist.Comment: 1 page, 1 figure, to appear in Phys. Rev. Let

Ice XII in its second regime of metastability

We present neutron powder diffraction results which give unambiguous evidence for the formation of the recently identified new crystalline ice phase[Lobban et al.,Nature, 391, 268, (1998)], labeled ice XII, at completely different conditions. Ice XII is produced here by compressing hexagonal ice I_h at T = 77, 100, 140 and 160 K up to 1.8 GPa. It can be maintained at ambient pressure in the temperature range 1.5 < T < 135 K. High resolution diffraction is carried out at T = 1.5 K and ambient pressure on ice XII and accurate structural properties are obtained from Rietveld refinement. At T = 140 and 160 K additionally ice III/IX is formed. The increasing amount of ice III/IX with increasing temperature gives an upper limit of T ~ 150 K for the successful formation of ice XII with the presented procedure.Comment: 3 Pages of RevTeX, 3 tables, 3 figures (submitted to Physical Review Letters

Amorphous-amorphous transition and the two-step replica symmetry breaking phase

The nature of polyamorphism and amorphous-to-amorphous transition is investigated by means of an exactly solvable model with quenched disorder, the spherical s+p multi-spin interaction model. The analysis is carried out in the framework of Replica Symmetry Breaking theory and leads to the identification of low temperature glass phases of different kinds. Besides the usual `one-step' solution, known to reproduce all basic properties of structural glasses, also a physically consistent `two-step' solution arises. More complicated phases are found as well, as temperature is further decreased, expressing a complex variety of metastable states structures for amorphous systems.Comment: 8 pages, 7 figures, longer version, new references adde

Relation between the High Density Phase and the Very-High Density Phase of Amorphous Solid Water

It has been suggested that high-density amorphous (HDA) ice is a structurally arrested form of high-density liquid (HDL) water, while low-density amorphous (LDA) ice is a structurally arrested form of low-density liquid (LDL) water. Recent experiments and simulations have been interpreted to support the possibility of a second "distinct" high-density structural state, named very high-density amorphous (VHDA) ice, questioning the LDL-HDL hypothesis. We test this interpretation using extensive computer simulations, and find that VHDA is a more stable form of HDA and that in fact VHDA should be considered as the amorphous ice of the quenched HDL.Comment: 5 pages, 4 fig

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 a 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.Comment: 5 pages, 4 figure

Interplay Between Time-Temperature-Transformation and the Liquid-Liquid Phase Transition in Water

We study the TIP5P water model proposed by Mahoney and Jorgensen, which is closer to real water than previously-proposed classical pairwise additive potentials. We simulate the model in a wide range of deeply supercooled states and find (i) the existence of a non-monotonic ``nose-shaped'' temperature of maximum density line and a non-reentrant spinodal, (ii) the presence of a low temperature phase transition, (iii) the free evolution of bulk water to ice, and (iv) the time-temperature-transformation curves at different densities.Comment: RevTeX4, 4 pages, 4 eps figure

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

Liquid-liquid equilibrium for monodisperse spherical particles

A system of identical particles interacting through an isotropic potential that allows for two preferred interparticle distances is numerically studied. When the parameters of the interaction potential are adequately chosen, the system exhibits coexistence between two different liquid phases (in addition to the usual liquid-gas coexistence). It is shown that this coexistence can occur at equilibrium, namely, in the region where the liquid is thermodynamically stable.Comment: 6 pages, 8 figures. Published versio