Search

We study nonequilibrium analogues of surface phase transitions in a minimal model of active particles in contact with a purely repulsive potential barrier that mimics a thin porous membrane. Under conditions of bulk motility-induced phase separation, the interaction strength ϵ_{w} of the barrier controls the affinity of the dense phase for the barrier region. We uncover clear signatures of a wetting phase transition as ϵ_{w} is varied. In common with its equilibrium counterpart, the character of this transition depends on the system dimensionality: a continuous transition with large density fluctuations and gas bubbles is uncovered in 2D while 3D systems exhibit a sharp transition absent of large correlations.

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behavior of equilibrium fluids with short-ranged pairwise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multibody effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of n-body effective interactions extracted from the configurational structure.

The atomic-scale structure, melting curve, and equation of state of liquid gallium has been measured to high pressure (p) and high temperature (T) up to 26 GPa and 900 K by in situ synchrotron x-ray diffraction. Ab initio molecular dynamics simulations up to 33.4 GPa and 1000 K are in excellent agreement with the experimental measurements, providing detailed insight at the level of pair distribution functions. The results reveal an absence of dimeric bonding in the liquid state and a continuous increase in average coordination number n[over ¯]_{Ga}^{Ga} from 10.4(2) at 0.1 GPa approaching ∼12 by 25 GPa. Topological cluster analysis of the simulation trajectories finds increasing fractions of fivefold symmetric and crystalline motifs at high p-T. Although the liquid progressively resembles a hard-sphere structure towards the melting curve, the deviation from this simple description remains large (≥40%) across all p-T space, with specific motifs of different geometries strongly correlating with low local two-body excess entropy at high p-T.

We model the thermodynamics of local structures within the hard sphere liquid at arbitrary volume fractions through the morphometric calculation of n-body correlations. We calculate absolute free energies of local geometric motifs in excellent quantitative agreement with molecular dynamics simulations across the liquid and supercooled liquid regimes. We find a bimodality in the density library of states where fivefold symmetric structures appear lower in free energy than fourfold symmetric structures and from a single reaction path predict a dynamical barrier which scales linearly in the compressibility factor. The method provides a new route to assess changes in the free energy landscape at volume fractions dynamically inaccessible to conventional techniques.

Among the key insights into the glass transition has been the identification of a nonequilibrium phase transition in trajectory space which reveals phase coexistence between the normal supercooled liquid (active phase) and a glassy state (inactive phase). Here, we present evidence that such a transition occurs in experiments. In colloidal hard spheres, we find a non-Gaussian distribution of trajectories leaning towards those rich in locally favored structures (LFSs), associated with the emergence of slow dynamics. This we interpret as evidence for a nonequilibrium transition to an inactive LFS-rich phase. Reweighting trajectories reveals a first-order phase transition in trajectory space between a normal liquid and a LFS-rich phase. We also find evidence for a purely dynamical transition in trajectory space.

Geometric frustration describes the inability of a local molecular arrangement, such as icosahedra found in metallic glasses and in model atomic glass formers, to tile space. Local icosahedral order, however, is strongly frustrated in Euclidean space, which obscures any causal relationship with the observed dynamical slowdown. Here we relieve frustration in a model glass-forming liquid by curving three-dimensional space onto the surface of a 4-dimensional hypersphere. For sufficient curvature, frustration vanishes and the liquid "freezes" in a fully icosahedral structure via a sharp "transition." Frustration increases upon reducing the curvature, and the transition to the icosahedral state smoothens while glassy dynamics emerge. Decreasing the curvature leads to decoupling between dynamical and structural length scales and the decrease of kinetic fragility. This sheds light on the observed glass-forming behavior in Euclidean space.