When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behaviour with no equilibrium counterpart. We study a two-dimensional particle model initially proposed to describe driven non-Brownian suspensions undergoing nonequilibrium absorbing phase transitions. We show that when the transition occurs at large density, the dynamics produces long-range crystalline order. In the ordered phase, long-range translational order is observed because equipartition of energy is lacking, phonons are suppressed, and density fluctuations are hyperuniform. Our study offers an explicit microscopic model where nonequilibrium violations of the Mermin-Wagner theorem stabilize crystalline order in two dimensions.

We use numerical simulations to study the dynamics of dense assemblies of self-propelled particles in the limit of extremely large, but finite, persistence times. In this limit, the system evolves intermittently between mechanical equilibria where active forces balance interparticle interactions. We develop an efficient numerical strategy allowing us to resolve the statistical properties of elastic and plastic relaxation events caused by activity-driven fluctuations. The system relaxes via a succession of scale-free elastic events and broadly distributed plastic events that both depend on the system size. Correlations between plastic events lead to emergent dynamic facilitation and heterogeneous relaxation dynamics. Our results show that dynamical behaviour in extremely persistent active systems is qualitatively similar to that of sheared amorphous solids, yet with some important differences.

The criticality of the jamming transition responsible for amorphous solidification has been theoretically linked to the marginal stability of a thermodynamic Gardner phase. While the critical exponents of jamming appear independent of the preparation history, the pertinence of Gardner physics far from equilibrium is an open question. To fill this gap, we numerically study the nonequilibrium dynamics of hard disks compressed toward the jamming transition using a broad variety of protocols. We show that dynamic signatures of Gardner physics can be disentangled from the aging relaxation dynamics. We thus define a generic dynamic Gardner cross-over regardless of the history. Our results show that the jamming transition is always accessed by exploring increasingly complex landscape, resulting in anomalous microscopic relaxation dynamics that remains to be understood theoretically.

The discovery of ultrastable glasses raises novel challenges about glassy systems. Recent experiments studied the macroscopic devitrification of ultrastable glasses into liquids upon heating but lacked microscopic resolution. We use molecular dynamics simulations to analyze the kinetics of this transformation. In the most stable systems, devitrification occurs after a very large time, but the liquid emerges in two steps. At short times, we observe the rare nucleation and slow growth of isolated droplets containing a liquid maintained under pressure by the rigidity of the surrounding glass. At large times, pressure is released after the droplets coalesce into large domains, which accelerates devitrification. This two-step process produces pronounced deviations from the classical Avrami kinetics and explains the emergence of a giant lengthscale characterizing the devitrification of bulk ultrastable glasses. Our study elucidates the nonequilibrium kinetics of glasses following a large temperature jump, which differs from both equilibrium relaxation and aging dynamics, and will guide future experimental studies.

The physics of the glass transition and amorphous materials continues to attract the attention of a wide research community after decades of effort. Supercooled liquids and glasses have been studied numerically since the advent of molecular dynamics and Monte Carlo simulations in the last century. Computer studies have greatly enhanced both experimental discoveries and theoretical developments and constitute an active and continually expanding research field. Our goal in this review is to provide a modern perspective on this area. We describe the need to go beyond canonical methods to attack a problem that is notoriously difficult in terms of time scales, length scales, and physical observables. We first summarise recent algorithmic developments to achieve enhanced sampling and faster equilibration using replica exchange methods, cluster and swap Monte Carlo algorithms, and other techniques. We then review some major recent advances afforded by these novel tools regarding the statistical mechanical description of the liquid-to-glass transition as well as the mechanical, vibrational and thermal properties of the glassy solid. We finally describe some important challenges for future research.