🔬 condensed matter

Theory of Nonequilibrium Crystallization and the Phase Diagram of Active Brownian Spheres

This paper presents a statistical mechanical framework and equations of state, guided by computer simulations, to construct the full phase diagram of active Brownian spheres, successfully predicting how activity shifts the crystallization transition and reproducing key coexistence features like the triple point.

Original authors: Daniel Evans, Ahmad K. Omar

Published 2026-07-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Daniel Evans, Ahmad K. Omar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a crowded dance floor where everyone is just standing still. If you push them together, they naturally jumble into a messy pile. But if you push them even tighter, they suddenly snap into a perfect, orderly grid, like soldiers standing at attention. In the world of physics, this is called "crystallization," and for a long time, scientists thought they understood the rules: it's all about how much space the particles have and how much they bump into each other. This is the story of "passive" matter—things that sit still unless you push them.

But what happens if the dancers start moving on their own? Imagine if every person on that dance floor had a tiny rocket strapped to their back, zooming around in random directions. This is "active matter," a field studying living things like bacteria or synthetic robots that consume energy to move. When these self-propelling particles try to form a crystal, the old rules break down. They don't just jam together; they behave in wild, unexpected ways. Understanding this is crucial because active matter is everywhere in nature, from the cells in your body to the swarms of bacteria in a pond. If we want to build new materials or understand life itself, we need to know how these energetic particles decide when to be a chaotic fluid and when to become an ordered solid.

This paper tackles the mystery of how "active" particles—specifically, little spheres that swim around on their own—decide to turn into a crystal. The authors, Daniel Evans and Ahmad K. Omar, realized that the old math used to predict when passive stuff freezes doesn't work for these energetic swimmers. They developed a brand new set of rules, a "statistical mechanical description," to track two things: how crowded the particles are (density) and how much they look like a neat crystal (crystallinity).

Think of the particles as a crowd of hyperactive kids in a gym. In a normal crowd, if you squeeze them enough, they naturally form a line. But these kids have energy; they run around, bump into each other, and keep moving. The authors found that as you crank up their energy (or "activity"), the point at which they decide to form a crystal shifts dramatically. Instead of forming a crystal at a moderate crowd level, they need to be packed much tighter—almost to the point of being completely jammed—before they finally snap into order. It's as if the kids are so busy running around that they ignore the urge to line up until the room is absolutely packed.

The team used computer simulations to test their new theory. They built a virtual world of these swimming spheres and watched what happened as they turned up the speed. They discovered that their new equations could predict exactly when the "solid" (the crystal) and the "fluid" (the messy crowd) would coexist. They even mapped out a complete "phase diagram," which is like a weather map for matter, showing exactly when you get a solid, a liquid, or a gas.

One of the most exciting findings is that their theory correctly predicts a "triple point." In normal physics, a triple point is a specific condition where a substance can be a solid, liquid, and gas all at the same time. The authors found that for these active swimmers, there is a specific level of energy where the solid, liquid, and gas phases can all hang out together. Their math matched the computer simulations almost perfectly, capturing the complex dance between the solid and the fluid, as well as the liquid and the gas.

Crucially, the paper argues that trying to use the old, "passive" rules to predict this behavior is not just slightly wrong; it's fundamentally broken. If you tried to use the standard formulas for these energetic particles, you would get the wrong answer, predicting that the fluid would be much denser than it actually is. The authors show that the "equilibrium" rules, which assume particles are just bumping into each other without extra energy, simply don't apply here. The energy the particles generate themselves changes the game entirely.

The authors are careful to note that while their theory works incredibly well for the range of activities they tested, it relies on some approximations, especially when the particles are moving very fast. They didn't just guess; they derived these rules from the physics of how the particles move and interact, and then checked them against thousands of computer runs. The result is a much clearer picture of how active matter organizes itself. It suggests that to understand these systems, we can't just look at how crowded they are; we have to account for the energy they are constantly pumping into the system. This work doesn't just solve a puzzle for one type of particle; it offers a new roadmap for understanding how any self-moving stuff might build structures, from the microscopic world of cells to the future of smart, self-assembling materials.

Technical Summary: Theory of Nonequilibrium Crystallization and the Phase Diagram of Active Brownian Spheres

Problem Statement
The crystallization of hard spheres at equilibrium is a canonical example of an entropically driven phase transition. However, recent observations indicate that activity—found in systems ranging from motile bacteria to self-propelling colloids—dramatically alters this order-disorder transition. Specifically, finite activity causes the solid phase density of monodisperse hard spheres to rapidly increase toward maximal face-centered cubic (fcc) packing ( $\phi_{solid} \approx 0.74$ ), a stark contrast to the equilibrium value ( $\phi_{solid} \approx 0.545$ ). Despite these observations, a theoretical description of active crystallization has remained elusive because traditional thermodynamic arguments, which underpin the understanding of passive freezing, are inapplicable to active, nonequilibrium systems. There is a critical need for a theory that elucidates the driving forces behind active crystallization and predicts the solid-fluid phase diagram without relying on equilibrium thermodynamics.

Methodology
The authors develop a multi-scale statistical mechanical theory for the crystallization of 3D active Brownian particles (ABPs). The approach involves the following key steps:

Field Formulation: The system is modeled using two order parameter fields: a conserved number density field $\rho(z)$ (or volume fraction $\phi(z)$ ) and a non-conserved, dimensionless "crystallinity" field $\psi_N(z)$ (or crystallinity density $\psi_V(z) = \psi_N \rho$ ). This allows for a description where both fluid and crystal phases possess uniform density fields, with crystallinity captured solely by the order parameter.
Derivation of Coexistence Criteria: Recognizing the absence of a nonequilibrium variational principle, the authors utilize a dynamical framework based on steady-state conditions ( $\dot{\rho}=0, \dot{\psi}_V=0$ $\overset{ρ}{˙} = 0, \dot{ψ}_{V} = 0$ ). They derive coexistence criteria by analyzing the balance of fluxes and generation terms for the density and crystallinity fields.
- They define a "crystallinity pseudopotential" ( $u_\psi$ ) and a "chemical pseudopotential" ( $u_\rho$ ).
- Unlike equilibrium systems where the Gibbs-Duhem relation allows for a rigorous Maxwell construction, the authors find that the nonequilibrium Gibbs-Duhem relation cannot be exactly satisfied for ABPs.
- Consequently, they propose an approximate generalized Maxwell construction. This involves solving for a transformation tensor to establish uniformity in the pseudopotentials and invoking an interfacial nonequilibrium Gibbs-Duhem relation to approximate the final coexistence criterion.
Equations of State (EOS): Guided by computer simulations (Brownian dynamics), the authors develop physics-guided empirical EOSs for the conservative interaction pressure ( $p^C_{bulk}$ ), the active pressure contribution ( $U$ ), and the equilibrium crystallinity path ( $\psi^*_V(\rho)$ ). These EOSs capture the dependence of the system on packing fraction ( $\phi$ ) and the ratio of active run-length to particle diameter ( $\ell_0/D$ ).

Key Contributions and Results

Shift in Order-Disorder Transition: The derived equations of state reveal that increasing activity shifts the order-disorder transition to higher packing fractions. The order-disorder volume fraction ( $\phi_{ODT}$ ) increases monotonically with activity, saturating at the random-close packing limit ( $\phi_{RCP} \approx 0.645$ ) at relatively low activities ( $\ell_0/D \approx 5$ ).
Reconstruction of the Phase Diagram: By applying the approximate coexistence criteria to the developed EOSs, the authors construct the full phase diagram of active Brownian spheres. The theory quantitatively recapitulates:
- The solid-fluid coexistence curve, correctly predicting the rapid approach of solid density to close-packing as activity increases.
- The liquid-gas coexistence curve (Motility-Induced Phase Separation, or MIPS).
- The solid-liquid-gas triple point and the subsequent transition to solid-gas coexistence at high activities.
Validation Against Simulations: The theoretical predictions show strong agreement with simulation data across a range of intermediate and high activities. The theory successfully captures the binodals for both MIPS and active crystallization, a feat not previously achieved without appealing to equilibrium thermodynamics.
Failure of Equilibrium Criteria: The paper demonstrates that naively applying equilibrium coexistence criteria (equal chemical potential and pressure via the standard Gibbs-Duhem relation) leads to qualitative errors at finite activity. Specifically, equilibrium criteria vastly overestimate fluid densities and incorrectly predict solid-liquid coexistence above the triple point, whereas the system actually exhibits solid-gas coexistence.

Significance
The authors claim that this work provides the first statistical mechanical description of active crystallization that quantitatively reproduces the full phase diagram of active Brownian spheres without relying on equilibrium thermodynamic arguments. By identifying the specific nonequilibrium coexistence criteria and developing appropriate equations of state, the theory elucidates the driving forces behind the shift in the order-disorder transition. The work offers a concrete path toward a general theory for nonequilibrium crystallization, highlighting that the equilibrium Gibbs-Duhem relation is formally incorrect for active systems and that distinct nonequilibrium criteria are required to accurately predict phase behavior.

Technical Summary: Theory of Nonequilibrium Crystallization and the Phase Diagram of Active Brownian Spheres

More like this