It was shown roughly thirty years ago that the density correlations of eigenvalues of large random matrices display a universal form, independent of most of the details of the distribution of the random matrix itself. We show that when the matrix elements evolve according to a Dyson Brownian motion, dynamical correlations retain a large degree of the universality found at equal times when expressed in terms of the characteristics of some partial differential equation in the complex plane.

The historical Mpemba effect involves a first-order phase transition. This has prompted the experimental realization of microscopic proxies in the form of a colloidal particle trapped in an asymmetric double well, for which the Mpemba effect has indeed been observed. We establish that the existence of the one-dimensional Mpemba effect for a polynomial potential is driven solely by the presence of a hard enough boundary, irrespective of the potential's double-well shape. We then show that the physics of the underlying Mpemba effect is governed not only by single-well physics but also by the high-temperature initial regime.

The Mpemba effect, a counterintuitive phenomenon where a hotter system relaxes faster than a colder one, has been widely observed in various nonequilibrium systems. Despite this progress, the fundamental structural features of the energy landscape required for its emergence remain a subject of debate. In this study, we investigate the conditions for the Mpemba effect within one-dimensional overdamped Langevin dynamics. We classify the potential landscapes based on the presence of single or double wells, their symmetry properties, and the existence of walls. We establish that the existence of the effect is primarily driven by the presence of boundaries, either hard or soft, rather than the specific internal structure of the potential landscape, such as metastability or the number of minima. By employing a spectral decomposition of the Fokker-Planck operator, we analyze the behavior of the first nontrivial eigenmode and demonstrate that its derivative acts as a Dirac delta peak in the low-temperature regime. This helps us elucidate the mechanism underlying the Mpemba effect: it appears as the interplay between this behavior and the initial population dynamics in a non-trivial way induced by the presence of the wall. Our analysis provides a unified classification across single- and double-well potentials, highlighting the crucial role of boundary conditions and asymmetry. Furthermore, we demonstrate that this framework allows for the engineering of potential landscapes capable of producing multistage Mpemba transitions.

We explore the edge flows that emerge at boundaries in nonequilibrium passive and active chiral colloidal fluids. We show that these complex interface currents obey an equation of state that relates their fluxes to bulk observables. For confined fluids, the edge flux is given by the average odd stress in the fluid. In phase-separated systems, the flux along the interface is given by the jump of the odd stress across the interface. We then use the equation of state to reveal, and contrast, the microscopic origins of the edge currents in passive and active systems.

The ultraslow dynamics of glass-formers has been explained by two views often considered as mutually exclusive: One invokes locally hindered mobility, and the other rests on the complexity of the configuration space. Here, we show that time evolution responds strongly to details of the dynamics by changing the speed of time flow: It has time-reparameterization softness. This finding reconciles both views: While local constraints reparameterize the flow of time, the global landscape determines relationships between different correlations at the same times. We show that modern algorithms developed to accelerate the relaxation to equilibrium act by changing the time reparameterization. Their success thus relies on their ability to exploit reparameterization softness. We conjecture that these results extend beyond the realm of glasses to the optimization of more general constraint satisfaction problems and to broader classes of algorithms.

We study the dynamics of weakly deformed interfaces separating two stable phases, starting from the fluctuating hydrodynamics of the phase-separating fields. Using a well-chosen definition for the interface and the dynamical-action formalism to represent path probabilities, we derive the linear relaxation of the interface and the fluctuations around it for a large class of models. Our method applies to equilibrium dynamics, where it recovers and complements existing results, but also extends to their non-equilibrium counterparts. We explain how non-linear terms can be systematically computed and illustrate their derivations in the case of (active) model A. We highlight the danger of a popular ansatz used to derive interface dynamics, which was rigorously established in equilibrium but is uncontrolled for active field theories.

Active matter systems comprise self-propelled particles that move on a substrate while leaving chemical trails that influence other particles through chemotaxis (e.g., slime-depositing bacteria). Orientational chemotaxis manifests as a torque that steers the particle toward the chemical gradient. As each particle is coupled to its own trail, the dynamics exhibits an instability: when the particle gently diffuses, it abruptly transitions to trajectories with a radius of curvature comparable to its own size, becoming apparently trapped. We argue that this trajectory instability occurs for any chemotactic coupling strength, either directly above some threshold or via an activated process below it.

The Young-Dupr & eacute; equation is a cornerstone of the theory of capillary and wetting phenomena, which governs the shape of droplets adsorbed on surfaces in equilibrium. However, many living and synthetic materials, such as swarming bacteria and active colloids, are composed of self-propelled particles that are inherently out of thermal equilibrium. The description of the wetting of surfaces by such active fluids thus requires a new framework. Here we develop an analogue to the Young-Dupr & eacute; equation for systems made of self-propelled particles. A key step is to define the liquid-gas surface tension of active fluids as the force exerted along the interface, which we show from first principles to be negative, even when active materials separate into stable liquid and gas phases. Our active Young-Dupr & eacute; equation explains why partial wetting appears in simulations where the surface tensions do not balance and reveals the underlying feedback mechanism: the interface is stable only because of steady flows, which are themselves generated by the parity symmetry-breaking interface. Unlike in passive fluids, where the droplets are scale-free, this feedback loop selects the sizes and shapes of adsorbed droplets in active materials. Our results outline a framework for understanding how active matter wets surfaces.

We rely on a hydrodynamic description of living tissues to describe the interface separating two of them with distinct constitutive properties. Using the difference in their homeostatic pressures as a control parameter, we show that the interface generates an emergent capillary surface tension that depends on the hydrodynamic scale, and that makes it very stable to large wavelengths perturbations. Using the difference of active forces the two tissues experience as a control parameter, we not only find that the front propagation mechanism shifts from the pushed wave to the Burgers wave, but we also find that the emergent surface tension is not sufficient to stabilize the interface at large enough drive and low enough viscosity.

We extend to soft repulsive interaction potentials a recently proposed irreversible swap algorithm originally designed for polydisperse hard spheres. The original algorithm performs rejection-free, irreversible, collective swap moves. We show that event-driven cluster updates of particle diameters can also be performed in continuous potentials by introducing a factorized Metropolis probability. However, the Metropolis factorization needed to deal with continuous potentials decreases the efficiency of the algorithm and mitigates the benefits of breaking detailed balance. This leads us to propose another irreversible swap algorithm using the standard Metropolis probability that accelerates the relaxation of soft sphere glasses at low temperatures compared to the original swap algorithm. We apply these efficient swap algorithms to produce very stable inherent structures with vibrational density of states lacking the quasilocalized excitations observed in conventional glasses.

When out-of-equilibrium particles interact by means of pairwise forces, their stationary distribution in general exhibits many-body interactions. In the particular case of active particles, it has been shown numerically that the Motility Induced Phase Separation cannot be explained by the effective attraction emerging from two isolated particles, thereby highlighting the role of multibody interactions. In this work, we study the thermodynamics of the Fox-UCNA approximation for active particles interacting by means of pairwise repulsive forces. Working at large space dimension we establish that multibody interactions up to infinite order are instrumental in giving rise to such collective phenomena as phase transitions. We recover a MIPS-like first order transition, but also find a liquid-liquid transition at somewhat lower persistence times. This new transition is connected to a spin glass phase of orientational-like degrees of freedom with disordered interactions set by the particle positions themselves.

We endow the elements of a random matrix drawn from the Gaussian Unitary Ensemble with a Dyson Brownian motion dynamics. We initialize the dynamics of the eigenvalues with all of them lumped at the origin, but one outlier. We solve the dynamics exactly which gives us a window on the dynamical scaling behavior at and around the Baik-Ben Arous-Péché transition. Amusingly, while the statics is well-known and accessible via the Hikami-Brézin integrals, our approach for the dynamics is explicitly based on the use of orthogonal polynomials.

We present a detailed derivation of the Langevin dynamics obeyed by a massive rigid body immersed in a chiral active bath. We show how the antisymmetric nature of the noise leads to an unusual relationship between the Langevin equation describing stochastic trajectories and the Fokker-Planck equation governing the time-evolution of the probability density. The chirality of the bath endows the object dynamics with odd diffusivity, odd mobility, and rotational ratchet effects that depend on the object symmetries. For rotationally-symmetric objects, we show that a hidden time-reversal symmetry leads to separate effective equilibrium descriptions for the translational and rotational degrees of freedom. Finally, starting from the bath dynamics, we construct a multipole expansion to quadrupolar order that allows predicting the far-field current and density modulation induced by the object on the bath.

When submerged in a chiral active bath, a passive object becomes a spinning ratchet imbued with odd transport properties. We present the most general Langevin dynamics for a rigid body in a chiral active bath, in the adiabatic limit of large object mass. For rotationally symmetric objects, odd diffusion and odd mobility are connected by an Einstein relation, that we show numerically to break down outside the adiabatic limit. As the object symmetry decreases, its dynamics becomes increasingly irreversible: a massive disk exhibits an effective equilibrium dynamics, while a rod admits distinct translational and rotational temperatures, and a wedge is fully irreversible.Conversely, this departure from equilibrium can be read in universal far-field currents and density modulations of the bath, which we measure numerically and derive analytically.

We present a detailed derivation of the Langevin dynamics obeyed by a massive rigid body immersed in a chiral active bath. We show how the antisymmetric nature of the noise leads to an unusual relationship between the Langevin equation describing stochastic trajectories and the Fokker-Planck equation governing the time evolution of the probability density. The chirality of the bath endows the object dynamics with odd diffusivity, odd mobility, and rotational ratchet effects that depend on the object symmetries. For rotationally symmetric objects, we show that a hidden time-reversal symmetry leads to separate effective equilibrium descriptions for the translational and rotational degrees of freedom. Finally, starting from the bath dynamics, we construct a multipole expansion to quadrupolar order that allows predicting the far-field current and density modulation induced by the object on the bath.

Equilibrium sampling of the configuration space in disordered systems requires algorithms that bypass the glassy slowing down of the physical dynamics. Irreversible Monte Carlo algorithms breaking detailed balance successfully accelerate sampling in some systems. We first implement an irreversible event-chain Monte Carlo algorithm in a model of continuously polydisperse hard disks. The effect of collective translational moves marginally affects the dynamics and results in a modest speedup that decreases with density. We then propose an irreversible algorithm performing collective particle swaps which outperforms all known Monte Carlo algorithms. We show that these collective swaps can also be used to prepare very dense jammed packings of disks.

We studied the nonequilibrium dynamics of a cycling three-state Potts model using simulations and theory. This model can be tuned from thermal-equilibrium to far-from-equilibrium conditions. At low cycling energy, the homogeneous dominant state cycles via nucleation and growth, while spiral waves are formed at high energy. For large systems, a discontinuous transition occurs from these cyclic homogeneous phases to spiral waves, while the opposite transition is absent. Conversely, these two modes can coexist for small systems. The waves can be reproduced by a continuum theory, and the transition can be understood from the competition between nucleation and growth.

The nature of the transition to collective motion in assemblies of aligning self-propelled particles remains a long-standing matter of debate. In this article, we focus on dry active matter and show that weak fluctuations suffice to generically turn second-order mean-field transitions into a 'discontinuous' coexistence scenario. Our theory shows how fluctuations induce a density-dependence of the polar-field mass, even when this effect is absent at mean-field level. In turn, this dependency on density triggers a feedback loop between ordering and advection that ultimately leads to an inhomogeneous transition to collective motion and the emergence of inhomogeneous travelling bands. Importantly, we show that such a fluctuation-induced first order transition is present in both metric models, in which particles align with neighbors within a finite distance, and in 'topological' ones, in which alignment is based on more complex constructions of neighbor sets. We compute analytically the noise-induced renormalization of the polar-field mass using stochastic calculus, which we further back up by a one-loop field-theoretical analysis. Finally, we confirm our analytical predictions by numerical simulations of fluctuating hydrodynamics as well as of topological particle models with either k-nearest neighbors or Voronoi alignment.

Monte Carlo simulations are widely employed to measure the physical properties of glass-forming liquids in thermal equilibrium. Combined with local Monte Carlo moves, the Metropolis algorithm can also be used to simulate the relaxation dynamics, thus offering an efficient alternative to molecular dynamics. Monte Carlo simulations are, however, more versatile because carefully designed Monte Carlo algorithms can more efficiently sample the rugged free energy landscape characteristic of glassy systems. After a brief overview of Monte Carlo studies of glass-formers, we define and implement a series of Monte Carlo algorithms in a three-dimensional model of polydisperse hard spheres. We show that the standard local Metropolis algorithm is the slowest and that implementing collective moves or breaking detailed balance enhances the efficiency of the Monte Carlo sampling. We use time correlation functions to provide a microscopic interpretation of these observations. Seventy years after its invention, the Monte Carlo method remains the most efficient and versatile tool to compute low-temperature properties in supercooled liquids.

The Young-Dupré equation is a cornerstone of the equilibrium theory of capillary and wetting phenomena. In the biological world, interfacial phenomena are ubiquitous, from the spreading of bacterial colonies to tissue growth and flocking of birds, but the description of such active systems escapes the realm of equilibrium physics. Here we show how a microscopic, mechanical definition of surface tension allows us to build an Active Young-Dupré equation able to account for the partial wetting observed in simulations of active particles interacting via pairwise forces. Remarkably, the equation shows that the corresponding steady interfaces do not result from a simple balance between the surface tensions at play but instead emerge from a complex feedback mechanism. The interfaces are indeed stabilized by a drag force due to the emergence of steady currents, which are themselves a by-product of the symmetry breaking induced by the interfaces. These currents also lead to new physics by selecting the sizes and shapes of adsorbed droplets, breaking the equilibrium scale-free nature of the problem. Finally, we demonstrate a spectacular consequence of the negative value of the liquid-gas surface tensions in systems undergoing motility-induced phase separation: partially-immersed objects are expelled from the liquid phase, in stark contrast with what is observed in passive systems. All in all, our results lay the foundations for a theory of wetting in active systems.