🔬 condensed matter

Hexatic Order Coupled with Thermal Noise Produces Bubbles in Two-Dimensional Active Matter

This study uses particle-based simulations to demonstrate that the formation of bubbles within the dense phase of two-dimensional active matter requires both hexatic order and a non-zero magnitude of thermal translational noise, which together facilitate the cooperative motion necessary for bubble generation.

Original authors: Luke Langford, Ahmad K. Omar

Published 2026-03-19

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Luke Langford, 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

The Big Picture: A Crowd of Self-Driving Cars

Imagine a giant, flat parking lot filled with thousands of tiny, self-driving cars (these are the "active particles"). These cars have a very simple rule: they never stop moving forward, but they can bump into each other and bounce off. They don't have brakes, and they don't talk to each other.

In a world where these cars are just driving around, you might expect them to spread out evenly. But, if there are enough of them and they are driving fast enough, something strange happens: they clump together. They form a massive, dense traffic jam (the "dense phase") surrounded by empty space (the "dilute phase"). This phenomenon is called Motility-Induced Phase Separation (MIPS).

The Mystery: The Bubbles in the Traffic Jam

Scientists have noticed something weird about these traffic jams in 2D simulations. Inside the massive clump of cars, bubbles of empty space start to form. It's like looking at a solid block of concrete and seeing little pockets of air floating inside it.

For a long time, scientists didn't know why these bubbles formed.

Theory A: Maybe the traffic jam is unstable, and the bubbles are just a glitch in the physics (like a reversed version of how soap bubbles usually merge).
Theory B: Maybe the cars are organizing into little neighborhoods with specific patterns (called "hexatic order"), and the bubbles form at the borders where these neighborhoods meet.

The Experiment: Two New Ingredients

The authors of this paper decided to test two specific things to solve the mystery:

Uniformity: Are the cars all exactly the same size, or are they different sizes? (Monodisperse vs. Polydisperse).
Jitter: Are the cars driving on a perfectly smooth, frictionless road, or is the road slightly bumpy, causing them to shake and wiggle randomly? (Thermal Noise).

The Discovery: The "Perfect Storm" for Bubbles

The researchers ran thousands of simulations and found a surprising recipe for creating bubbles. You need both of the following ingredients, but in a very specific way:

1. The Cars Must Be Identical (Uniformity)

If the cars are all different sizes (some big trucks, some tiny compacts), they get stuck in a messy, disorganized pile. In this messy pile, no bubbles form. The traffic jam is just a solid, chaotic block.

The Metaphor: Think of a pile of mixed-up rocks. They fit together tightly with no gaps. But if you have a pile of identical marbles, they can arrange themselves into neat, repeating patterns. The bubbles only appear when the cars are identical marbles that can organize themselves.

2. The Road Must Be Slightly Bumpy (Thermal Noise)

This is the most shocking part. The researchers thought that the "bumpiness" of the road (thermal noise) didn't matter much because the cars were driving so fast on their own.

The Result: If the road is perfectly smooth (zero noise), the cars get stuck in a rigid grid. They can't move around enough to create bubbles.
The Twist: If the road has even the tiniest amount of bumpiness (a tiny bit of random shaking), the bubbles appear instantly!
The Analogy: Imagine a crowd of people holding hands in a tight circle. If they stand perfectly still, they are a solid wall. But if they all start doing a tiny, random "jitter" dance (even a very small one), they suddenly gain the ability to wiggle, shift, and create gaps (bubbles) between them.

The "Secret Sauce": Cooperative Motion

Why does a tiny bit of shaking make such a huge difference?
The paper found that when the cars are identical and the road is slightly bumpy, the cars start dancing together. They don't just wiggle randomly; they move in groups. A whole neighborhood of cars will shift and rotate together like a single unit.

This cooperative motion is the key. It allows the dense traffic jam to fluidly rearrange itself, creating and maintaining those bubbles. Without the tiny bit of "jitter," the cars are too stiff and rigid to cooperate.

The Conclusion

The paper solves the mystery by showing that bubbles in 2D active matter aren't just a random glitch. They are a specific phenomenon that requires:

Order: The particles must be uniform enough to form structured neighborhoods (hexatic order).
Jitter: The particles need a tiny amount of random thermal noise to unlock their ability to move cooperatively.

In short: You can't get bubbles in a 2D active system if your particles are too messy (different sizes) or too stiff (no random noise). You need a crowd of identical particles that are just "jittery" enough to dance together and make room for the bubbles.

1. Problem Statement

The paper addresses a persistent open question in the field of active matter: the origin of bubbles (gas-phase pockets) observed within the dense liquid phase of two-dimensional (2D) motility-induced phase separation (MIPS).

Context: Purely repulsive, self-propelled particles (Active Brownian Particles or ABPs) undergo phase separation into dense and dilute regions. In 2D simulations, the dense phase is often decorated with bubbles.
Existing Theories: Previous explanations have been conflicting:
- Reversed Ostwald Ripening: Theoretical models (e.g., Active Model B+) suggest negative surface tension causes bubbles to nucleate and grow, reversing the standard coarsening process.
- Topological Defects: Empirical observations suggest bubbles form at the grain boundaries of hexatic domains (quasi-crystalline order), implying a structural origin rather than a purely thermodynamic one.
The Gap: It remains unclear whether hexatic order is sufficient for bubble formation, what role thermal noise plays (often assumed negligible in active matter), and why these bubbles are absent in 3D simulations.

2. Methodology

The authors employed large-scale particle-based simulations using the HOOMD-Blue software package to disentangle the effects of structural order and thermal noise.

Model System: 2D Active Brownian Particles (ABPs) governed by overdamped Langevin equations.
- Forces: Pairwise repulsive interactions via a Weeks-Chandler-Anderson (WCA) potential.
- Activity: Self-propulsion speed $U_0$ and rotational diffusion $D_R$ .
- Thermalization: Translational noise characterized by a dimensionless coefficient $T \equiv D_T / (D_R \sigma^2)$ , where $D_T$ is translational diffusivity.
Key Variables Manipulated:
1. Polydispersity: The authors compared monodisperse systems (uniform particle size) vs. polydisperse systems (Gaussian distribution of diameters, 9% standard deviation). Polydispersity was used to systematically suppress hexatic order.
2. Thermal Noise ( $T$ ): Simulations were run across a range of $T$ , including the athermal limit ( $T=0$ ) and various finite values ( $T > 0$ ), down to extremely small magnitudes ( $T = 1/100$ ).
Analysis Metrics:
- Bubble Distribution $n(A)$ : The number of bubbles of area $A$ .
- Hexatic Order: Spatial correlation of the local bond-orientational order parameter ( $\psi_6$ ).
- Dynamics: Velocity orientation correlations ( $Q(r)$ ) and Mean Squared Displacement (MSD) to quantify collective motion.

3. Key Contributions

Disentangling Order and Noise: The study isolates the specific contributions of hexatic order and translational thermal noise, demonstrating that neither is sufficient on its own; both are required.
The Necessity of Non-Zero Thermal Noise: The authors prove that even a vanishingly small amount of translational noise is critical for bubble formation. In the strict athermal limit ( $T=0$ ), bubbles become statistically rare, regardless of the presence of hexatic order.
Cooperative Motion Mechanism: The paper identifies that thermal noise triggers cooperative motion within the dense phase (specifically within hexatic domains), which appears to be the driving mechanism for bubble generation.

4. Key Results

A. The Role of Hexatic Order (Polydispersity)

Monodisperse Systems: Exhibit long-range hexatic order.
Polydisperse Systems: The introduction of size polydispersity destroys long-range hexatic order (correlations decay within a few particle diameters).
Result: In polydisperse systems, no significant bubble formation occurs, even when thermal noise is present ( $T > 0$ ). This confirms that hexatic order is a necessary condition for bubble formation.

B. The Role of Thermal Noise (Athermal vs. Thermal)

Athermal Limit ( $T=0$ ): Even in monodisperse systems with strong hexatic order, bubbles are extraordinarily rare. The bubble size distribution $n(A)$ decays sharply, and large bubbles are essentially non-existent.
Finite Thermal Noise ( $T > 0$ ): Introducing even a tiny amount of noise (e.g., $T = 1/100$ $T = 1/100$ ) leads to a dramatic increase in bubble population.
- The distribution $n(A)$ shifts to a power-law decay ( $n(A) \sim A^{-1.45}$ ), consistent with self-organized criticality.
- The fraction of gas phase within bubbles ( $\phi_B$ ) increases by three orders of magnitude compared to the athermal case.
Conclusion: Hexatic order is necessary but insufficient; translational thermal noise is the trigger for bubble generation.

C. Dynamical Insights

Collective Motion: In thermal monodisperse systems, particles within hexatic domains exhibit collective translation and reorientation. This cooperative motion is absent in athermal systems (where dynamics are sluggish and mediated by vacancies) and polydisperse systems (where transport is rapid but uncorrelated).
Velocity Correlations: The spatial range of velocity correlations increases with thermal noise, suggesting that noise induces a cooperative dynamic state essential for nucleating bubbles.

D. Size Distribution

While the number of bubbles depends heavily on both hexatic order and noise, the size distribution (probability $P(A)$ ) for small bubbles appears to follow a similar algebraic decay across thermalized systems (monodisperse and polydisperse), suggesting that noise dictates the size scaling, while order dictates the population density.

5. Significance and Implications

Resolution of the "Bubble Paradox": The findings explain why bubbles are observed in 2D MIPS but not in 3D. In 3D, the dense phase is disordered (lacking hexatic order), and without the specific coupling of order and noise, bubbles do not form.
Refutation of Purely Deterministic Models: The results challenge models that assume active matter can be treated purely deterministically (athermal). Even infinitesimal thermal noise fundamentally alters the non-equilibrium phase behavior.
Mechanism for Bubble Formation: The paper proposes that bubbles are not merely a result of reversed Ostwald ripening driven by negative surface tension, but are dynamically generated by cooperative motion within hexatic domains, which requires thermal fluctuations to activate.
Future Directions: The authors suggest that coarse-grained continuum models (like Active Model B+) must be revised to include long-range dynamical correlations and the specific interplay between structural order and thermal noise to accurately predict bubble statistics.

In summary, the paper establishes that bubbles in 2D active matter are a cooperative phenomenon arising from the coupling of hexatic structural order and translational thermal noise, where the noise acts as a catalyst for the collective dynamics required to nucleate and sustain bubbles.