Diagnosing broken ergodicity using an energy fluctuation metric

Abstract

The Mountain and Thirumalai energy fluctuation metric, Omega(t), has been used to study the effective ergodicity of 60- and 256-atom binary Lennard-Jones mixtures in order to determine the reliability of the calculated diffusion constants at different energies. A plot of Omega(t) against 1time allows the identification of two distinct regimes: ergodic supercooled liquids, where Omega(t) approaches zero, and nonergodic glassy states, where Omega(t) asymptotically approaches a nonzero value on the molecular dynamics time scale. This approach seems to be more appropriate than attempting to define a threshold value for Omega(t)/Omega(0). The behavior of systems between these two limits, which are nonergodic on the time scale considered but may be approaching ergodicity, was examined for a range of simulation times. The calculated diffusion constants change as effective ergodicity is approached, moving closer to the Vogel-Tammann-Fulcher fit defined by higher-energy systems that are already considered to be effectively ergodic. Using the form of the decay of the metric as a measure of ergodicity, we have been able to reproduce the trend in fragility obtained by Sastry for a 256-atom system [Nature (London) 409, 164 (2001)], correcting some of our earlier results [J. Chem. Phys. 120, 8314 (2004)].