We report a computer simulation study of the glass transition for water. To mimic the difference between standard and hyperquenched glass, we generate glassy configurations with different cooling rates and calculate the T dependence of the specific heat on heating. The absence of crystallization phenomena allows us, for properly annealed samples, to detect in the specific heat the simultaneous presence of a weak pre-peak ("shadow transition"), and an intense glass transition peak at higher temperature. We discuss the implications for the currently debated value of the glass transition temperature of water. We also compare our simulation results with the Tool-NarayanaswamyMoynihan phenomenological model.
PACS numbers:Much recent research has focused on the properties of glassy water, the most common form of water in the universe, which can exist in more than one distinct amorphous form [1,2,3]. The conversion between different glass structures, the different routes producing glass structures, and the relation between the liquid and the glass phases are under active debate.A particularly relevant aspect of this debate concerns the identification of the glass transition temperature T g at ambient pressure and the magnitude of the associated jump of the specific heat, an issue which has relevance also for determining the fragility of water. Extrapolation of T g in binary aqueous solutions, in the limit of vanishing solute concentration, provides the estimate T g ≈ 136 K [4]. Early differential scanning calorimetry (DSC) studies report conflicting results. Some experiments detect the glass transition [5] but others do not [6]. An exothermal peak in the specific heat of properly-annealed hyperquenched water supports the estimate T g ≈ 136 K [7], with a specific heat jump of 1.6 − 1.9 J/mol/K. This T g value [8,9] has been recently debated [10,11,12]. It has been suggested[12] that the small peak measured in Ref.[7] is a pre-peak typical of annealed hyperquenched samples preceding the true glass transition located at T g ≈ 165 K. Assigning T g ≈ 165 K would explain many of the puzzles related to the glass transition in water [9,10,12]. Unfortunately, the T g ≈ 165 K proposal can not be experimentally tested due to the homogeneous nucleation of the crystal phase at T × ≈ 150 K.Here we report a numerical study of the temperature dependence of the specific heat across the glass-to-liquid transition for the extended simple point charge (SPC/E) model for water. We analyze the effects both of the cooling rate and of annealing ("aging") before heating the glass, since both effects are important for determining T g [13,14], and both effects have been studied extensively in many materials [15,16]. Numerical studies are particularly suited since crystallization does not take place on the time scale probed in simulations. With an appropriate choice of the heating and cooling rates to mimic the experimental conditions, we show that both the shadow and the glass transition peaks can be resolved in the same heating scan. Finally, we...