Unifying fluctuation-dissipation temperatures of slow-evolving nonequilibrium systems from the perspective of inherent structures

Understanding the glassy dynamics from melting temperatures in binary glass-forming liquids

It is natural to expect that small particles in binary mixtures move faster than large ones. However, in binary glass-forming liquids with soft-core particle interactions, we observe the counterintuitive dynamic reversal between large and small particles along with the increase of pressure by performing molecular dynamics simulations. The structural relaxation (dynamic heterogeneity) of small particles is faster (weaker) than large ones at low pressures, but becomes slower (stronger) above a crossover pressure. In contrast, this dynamic reversal never happens in glass-forming liquids with hard-core interactions. We find that the difference of the effective melting temperatures felt by large and small particles can be used to understand the dynamic reversal. In binary mixtures, we derive effective melting temperatures of large and small particles simply from the conversion of units and find that particles with a higher effective melting temperature usually undergo a slower and more heterogeneous relaxation. The presence (absence) of the dynamic reversal in soft-core (hard-core) systems is simply due to the non-monotonic (monotonic) behavior of the melting temperature as a function of pressure. Interestingly, by manipulating the relative softness between large and small particles, we obtain a special case of soft-core systems, in which large particles always have higher effective melting temperatures than small ones. As a result, the dynamic reversal is totally eliminated. Our work provides another piece of evidence of the underlying connections between the properties of non-equilibrium glass-formers and equilibrium crystal-formers.

Discontinuous instabilities in disordered solids

Under a sufficiently large load, a solid material will flow via rearrangements, where particles change neighbors. Such plasticity is most easily described in the athermal, quasistatic limit of zero temperature and infinitesimal loading rate, where rearrangements occur only when the system becomes mechanically unstable. For disordered solids, the instabilities marking the onset of rearrangements have long been believed to be fold instabilities, in which an energy barrier disappears and the frequency of a normal mode of vibration vanishes continuously. Here, we report that there exists another, anomalous, type of instability caused by the breaking of a "stabilizing bond," whose removal creates an unstable vibrational mode. For commonly studied systems, such as those with harmonic finite-range interparticle interactions, such "discontinuous instabilities" are not only inevitable, they often dominate the modes of failure. Stabilizing bonds are a subset of all the bonds in the system and are prevalent in disordered solids generally. Although they do not trigger discontinuous instabilities in systems with vanishing stiffness at the interaction cutoff, they are, even in those cases, local indicators of incipient mechanical failure. They therefore provide an accurate structural predictor of instabilities not only of the discontinuous type but of the fold type as well.

Equilibrium states corresponding to targeted hyperuniform nonequilibrium pair statistics

The Zhang-Torquato conjecture [G. Zhang and S. Torquato, Phys. Rev. E, 2020, 101, 032124.] states that any realizable pair correlation function g₂(r) or structure factor S(k) of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To further test and study this conjecture, we consider two singular nonequilibrium models of recent interest that also have the exotic hyperuniformity property: a 2D "perfect glass" and a 3D critical absorbing-state model. We find that each nonequilibrium target can be achieved accurately by equilibrium states with effective one- and two-body potentials, lending further support to the conjecture. To characterize the structural degeneracy of such a nonequilibrium-equilibrium correspondence, we compute higher-order statistics for both models, as well as those for a hyperuniform 3D uniformly randomized lattice (URL), whose higher-order statistics can be very precisely ascertained. Interestingly, we find that the differences in the higher-order statistics between nonequilibrium and equilibrium systems with matching pair statistics, as measured by the "hole" probability distribution, provide measures of the degree to which a system is out of equilibrium. We show that all three systems studied possess the bounded-hole property and that holes near the maximum hole size in the URL are much rarer than those in the underlying simple cubic lattice. Remarkably, upon quenching, the effective potentials for all three systems possess local energy minima (i.e., inherent structures) with stronger forms of hyperuniformity compared to their target counterparts. Our methods are expected to facilitate the self-assembly of tunable hyperuniform soft-matter systems.