The progression of thermodynamic anomalies in MX networks with local tetrahedral geometries.

Key thermodynamic anomalies in density and compressibility, as well as the related stability limits, are determined using an ionic model for BeF which includes many-body polarization terms. BeF is chosen as an example of an archetypal network-forming system whose structure can be rationalised in terms of connected local tetrahedral coordination polyhedra. The anion dipole polarizability (which effectively controls the bond angles linking neighbouring tetrahedra) is used as a single free parameter in order to help rationalise the changes in the anomaly locations in phase space, whilst all other potential parameters remain fixed. The anomalies and stability limits systematically shift to lower temperature and higher pressure as the anion polarizability is increased. At high dipole polarizabilities the temperature of maximum density anomaly locus becomes suppressed into the supercooled regime of the phase space. The movements of the anomaly loci are analysed in terms of the network structure and the correlation with the inter-tetrahedral bond angles is considered. The high sensitivity of the anomalies to the details of the potential models applied is discussed with reference to previous works on related systems. The relationship to analogous studies on Stillinger-Weber liquids is discussed.