skip to content
 

When designing practical applications of materials, it is rather useful to know to what extent they are thermodynamically stable. For example, the use of tin to seal fuel cans has been alleged to have contributed to the failure of Robert Scott's second expedition to the South Pole: metallic β-tin is thermodynamically unstable with respect to the brittle α-tin below 13 °C. The transition does not generally occur at 13 °C however, since β-tin can remain metastable to significantly lower temperatures. However, the probability of a rare nucleation event occurring increases rapidly with supercooling, leading to the problem of ‘tin pest’. A good understanding of thermodynamics is sometimes quite literally a matter of life and death.

Indeed, knowing the conditions under which a substance is stable, but also where it is metastable, can be important in many other situations, even if they are perhaps a little less immediately deadly. Separation processes rely on knowledge of the underlying phase diagram, and pharmaceuticals must usually be crystallised into a particular polymorph, since this determines their solubility and hence activity. Predicting the phase stabilities of substances is therefore central to computational materials science and chemistry. However, it is a less straightforward matter than it seems, since thermodynamic stability is controlled not only by energetic considerations, but also by entropy, which is much more difficult to deal with, especially for systems whose configurational entropy is inaccessible to typical simulation timescales (such as proton disorder in ice or orientational ordering in quasicrystals).

Predicting the phase behaviour of substances is an especially challenging test of our computational modelling because very small differences in the underlying potential energy can drastically affect the phase behaviour. This makes it challenging to get right on the one hand, but at the same time this means that the prediction of phase diagrams can serve as a useful means of parameterising potentials: if an empirical potential is good enough to give rise to phase behaviour consistent with experiment, it is likely that it can also describe many other properties of the material of interest well, and so this can give us considerable confidence in the predictions of such a potential.

Related Publications 

Phase diagrams-Why they matter and how to predict them
PY Chew, A Reinhardt – Journal of Chemical Physics (2022) 158, 030902
A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals
A Reinhardt, PY Chew, B Cheng – J Chem Phys (2023) 159, 184110
Predicting the phase diagram of titanium dioxide with random search and pattern recognition.
A Reinhardt, CJ Pickard, B Cheng – Physical Chemistry Chemical Physics (2020) 22, 12697
Numerical calculation of free-energy barriers for entangled polymer nucleation
X Tang, F Tian, T Xu, L Li, A Reinhardt – J Chem Phys (2020) 152, 224904
Phase behavior of empirical potentials of titanium dioxide
A Reinhardt – Journal of Chemical Physics (2019) 151, 064505
Self-assembly of two-dimensional binary quasicrystals: a possible route to a DNA quasicrystal
A Reinhardt, JS Schreck, F Romano, JPK Doye – J Phys Condens Matter (2017) 29, 014006
Computing phase diagrams for a quasicrystal-forming patchy-particle system
A Reinhardt, F Romano, JPK Doye – Physical review letters (2013) 110, 255503
Re-entrant phase behavior for systems with competition between phase separation and self-assembly.
A Reinhardt, AJ Williamson, JPK Doye, J Carrete, LM Varela, AA Louis – Journal of Chemical Physics (2011) 134, 104905