Search for a liquid-liquid critical point in models of silica

Erik Lascaris, Mahin Hemmati, Sergey V. Buldyrev, H. Eugene Stanley, Charles Angell

Research output: Contribution to journalArticlepeer-review

61 Scopus citations


Previous research has indicated the possible existence of a liquid-liquid critical point (LLCP) in models of silica at high pressure. To clarify this interesting question we run extended molecular dynamics simulations of two different silica models (WAC and BKS) and perform a detailed analysis of the liquid at temperatures much lower than those previously simulated. We find no LLCP in either model within the accessible temperature range, although it is closely approached in the case of the WAC potential near 4000 K and 5 GPa. Comparing our results with those obtained for other tetrahedral liquids, and relating the average Si-O-Si bond angle and liquid density at the model glass temperature to those of the ice-like β-cristobalite structure, we conclude that the absence of a critical point can be attributed to insufficient "stiffness" in the bond angle. We hypothesize that a modification of the potential to mildly favor larger average bond angles will generate a LLCP in a temperature range that is accessible to simulation. The tendency to crystallize in these models is extremely weak in the pressure range studied, although this tendency will undoubtedly increase with increasing stiffness.

Original languageEnglish (US)
Article number224502
JournalJournal of Chemical Physics
Issue number22
StatePublished - Jun 14 2014

ASJC Scopus subject areas

  • General Physics and Astronomy
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Search for a liquid-liquid critical point in models of silica'. Together they form a unique fingerprint.

Cite this