Glass and liquid phase diagram of a polyamorphic monatomic system.

We perform out-of-equilibrium molecular dynamics (MD) simulations of a monatomic system with Fermi-Jagla (FJ) pair potential interactions. This model system exhibits polyamorphism both in the liquid and glass state. The two liquids, low-density (LDL) and high-density liquid (HDL), are accessible in equilibrium MD simulations and can form two glasses, low-density (LDA) and high-density amorphous (HDA) solid, upon isobaric cooling. The FJ model exhibits many of the anomalous properties observed in water and other polyamorphic liquids and thus, it is an excellent model system to explore qualitatively the thermodynamic properties of such substances. The liquid phase behavior of the FJ model system has been previously characterized. In this work, we focus on the glass behavior of the FJ system. Specifically, we perform systematic isothermal compression and decompression simulations of LDA and HDA at different temperatures and determine "phase diagrams" for the glass state; these phase diagrams varying with the compression/decompression rate used. We obtain the LDA-to-HDA and HDA-to-LDA transition pressure loci, P(LDA-HDA)(T) and P(HDA-LDA)(T), respectively. In addition, the compression-induced amorphization line, at which the low-pressure crystal (LPC) transforms to HDA, P(LPC-HDA)(T), is determined. As originally proposed by Poole et al. [Phys. Rev. E 48, 4605 (1993)] simulations suggest that the P(LDA-HDA)(T) and P(HDA-LDA)(T) loci are extensions of the LDL-to-HDL and HDL-to-LDL spinodal lines into the glass domain. Interestingly, our simulations indicate that the P(LPC-HDA)(T) locus is an extension, into the glass domain, of the LPC metastability limit relative to the liquid. We discuss the effects of compression/decompression rates on the behavior of the P(LDA-HDA)(T), P(HDA-LDA)(T), P(LPC-HDA)(T) loci. The competition between glass polyamorphism and crystallization is also addressed. At our "fast rate," crystallization can be partially suppressed and the glass phase diagram can be related directly with the liquid phase diagram. However, at our "slow rate," crystallization cannot be prevented at intermediate temperatures, within the glass region. In these cases, multiple crystal-crystal transformations are found upon compression/decompression (polymorphism).