🔬 condensed matter

Insertion space in repulsive active matter

This paper extends the statistical geometry of insertion space from equilibrium hard spheres to repulsive active matter in one and two dimensions, deriving closed-form expressions that reveal how activity increases total insertion volume, enhances connectivity, and encodes signatures of collective phase behaviors.

Original authors: Luke K. Davis, Karel Proesmans

Published 2026-07-21

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Luke K. Davis, Karel Proesmans

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 trying to find a spot to stand without bumping into their neighbors. In the quiet, orderly world of physics, we've long known how to calculate the "empty space" available for a new dancer to join the party. This is called the "insertion space," and for decades, scientists have used it to predict how gases and liquids behave when they are calm and in balance. But what happens when the dancers aren't just standing still or moving randomly? What if they are self-propelled robots, swimming bacteria, or tiny particles that generate their own energy to zoom around? This is the world of "active matter." Unlike a passive crowd that just jiggles with heat, active matter is alive with motion, constantly burning energy to push itself forward. Understanding how much room is left in a crowd of these energetic, self-driving particles is a huge puzzle. It matters because these systems are everywhere, from the way our cells organize themselves to the future of swarms of tiny robots. If we can figure out how much "room to move" exists in a chaotic, self-driving crowd, we might finally understand how to store them, control them, or even predict when they will suddenly clump together.

In this new study, researchers Luke K. Davis and Karel Proesmans decided to tackle this puzzle by asking a simple question: If you have a box full of these energetic, repulsive particles that hate to touch each other, how much space is actually left for you to slip one more particle in? They looked at this "insertion space" in one-dimensional lines and two-dimensional flat surfaces, using both mathematical formulas and computer simulations.

The team discovered something surprising and counter-intuitive. In a calm, passive crowd, adding more energy (like shaking the box) usually just makes the particles jitter more, but it doesn't necessarily create more room. However, in the world of active matter, the act of the particles zooming around and pushing themselves actually increases the total amount of empty space available. It's as if the dancers, by constantly running and turning, accidentally create larger gaps between them than they would if they were just standing still. The researchers found that as the particles become more active, the "insertion space" doesn't just get bigger; it also tends to stay more connected. Instead of the empty space breaking up into many tiny, isolated islands, the active particles seem to keep the gaps linked together, forming a more open network.

They tested this idea using two different models of active particles: "run-and-tumble" particles (which zoom in a straight line, stop, and then randomly pick a new direction) and "active Ornstein-Uhlenbeck" particles (which move with a smooth, wobbly persistence). In both cases, whether they were looking at a single line of particles or a flat plane, the results were consistent. The more active the particles were, the more total empty volume existed. Interestingly, while the total space grew, the number of separate empty pockets actually went down. This suggests that activity merges small gaps into fewer, larger, and more connected open areas.

The researchers also found that these changes in empty space leave a clear "signature" that hints at when the particles are about to change their behavior. For instance, in two dimensions, they noticed that the way the empty space changed suddenly shifted at specific packing densities—exactly where scientists know these particles usually start to freeze or form crystal-like structures. This suggests that measuring the "room to move" could be a new way to detect when a swarm of active particles is about to undergo a phase change, like turning from a fluid into a solid.

While the math worked perfectly for the one-dimensional line, the two-dimensional simulations showed that the relationship gets a bit more complex as the particles get crowded, with the empty space sometimes behaving in non-linear ways. However, the core finding remains robust: activity makes the crowd more spacious. The authors emphasize that these results come from simulations and theoretical derivations, not physical experiments yet, but they believe the findings are solid enough to be tested in real-world labs with things like self-propelled colloids. This work provides the first quantitative map of the "statistical geometry" of active matter, offering a new lens to see how energy and motion reshape the very space these particles occupy.

Technical Summary: Insertion Space in Repulsive Active Matter

Problem Statement
The paper addresses the fundamental question of how to characterize the "room" available to insert an additional object into a system of repulsive active particles. While the statistical geometry of insertion space (IS) is well-established for equilibrium hard spheres—where exact relations link geometric quantities (like total insertion volume and cavity size) to thermodynamic state variables (compressibility and chemical potential)—this framework has not been extended to active matter. Active systems, such as self-propelled colloids or biological cells, are far from equilibrium and exhibit phenomena like motility-induced phase separation (MIPS). The authors investigate whether the geometric metrics of insertion space can provide a quantitative foundation for understanding the thermodynamics and phase behavior of these non-equilibrium systems.

Methodology
The study employs a combination of analytical derivations and numerical simulations across one (d=1) and two (d=2) spatial dimensions.

Models: The authors analyze two primary models of repulsive active matter:
1. Run-and-Tumble Particles (RTPs): Studied both on a lattice (discrete sites) and off-lattice (continuous space with WCA potential).
2. Active Ornstein-Uhlenbeck Particles (AOUPs): Studied off-lattice in continuous space.
Analytical Approach (d=1): For the 1D lattice RTP model, the authors derive closed-form expressions for IS metrics. They utilize a mean-field approximation for weak activity ( $K \lesssim 1$ ) and a coalescence-fragmentation picture (based on the DCR model) for high activity ( $K \gg 1$ ). These analytical results are expressed in terms of the contact probability ( $P_C$ ) and packing fraction ( $\eta$ ).
Numerical Approach:
- d=1: Both on-lattice and off-lattice simulations are performed to validate the analytical predictions.
- d=2: Due to the complexity of defining cavities in higher dimensions (where multiple particles can bound a single cavity), the authors develop a voxel-based algorithm. This method partitions space into small voxels, removes those within exclusion zones, and calculates IS metrics (cavity size, number, total volume, and surface area) efficiently. Simulations are conducted for off-lattice RTPs.

Key Contributions and Results

Extension of Equilibrium Relations: The paper establishes the first quantitative link between the statistical geometry of insertion space and active matter. It defines geometric observables: the mean insertion cavity size ( $\langle \tilde{\Delta} \rangle$ ), the number of cavities ( $\langle N_\Delta \rangle$ ), the total insertion volume ( $\langle V_I \rangle$ ), and the geometric compressibility ( $Z_G$ ) and chemical potential ( $\mu_G$ ).
Analytical Expressions in 1D:
- The authors derive exact relations for 1D equilibrium hard spheres (Tonks gas) and extend them to active RTPs.
- Weak Activity: Closed-form expressions are derived up to $O(K^3)$ , showing corrections to the equilibrium values based on the contact probability.
- High Activity: A distinct regime is identified where the system forms long, correlated domains. Here, the contact probability and cavity sizes scale with the persistence length, leading to simplified relations where $\langle V_I \rangle/V \approx 1 - \eta$ in the limit of extreme activity.
Universal Effects of Activity:
- Increased Total Insertion Volume: Across all packing fractions and activity levels, the total insertion volume ( $\langle V_I \rangle$ ) is strictly higher in active systems compared to their passive counterparts.
- Reduced Cavity Count: Activity tends to reduce the total number of distinct insertion cavities ( $\langle N_\Delta \rangle$ ), implying that the available space becomes more connected (fewer, larger cavities).
- Geometric Metrics: The geometric compressibility ( $Z_G$ ) and geometric chemical potential ( $\mu_G$ ) are generally lower for active systems than for passive ones, suggesting a "softer" effective equation of state in terms of insertion geometry.
2D Findings and Phase Signatures:
- In 2D, the total insertion volume remains higher with increased activity, and the space remains more connected.
- Non-monotonicity: Unlike in 1D, the mean cavity size in 2D is not strictly monotonic with activity at high packing fractions, showing a "shelf" behavior.
- Phase Transition Signatures: The metrics exhibit distinct features at packing fractions corresponding to known phase transitions in hard discs (freezing, liquid-hexatic coexistence, and jamming). Specifically, the slope of the geometric compressibility ratio ( $Z_G/Z_G^{eq}$ ) changes sign near these transition points, suggesting that insertion space geometry encodes information about collective phase behaviors.
Model Independence: The findings hold for both RTPs and AOUPs, and the high-activity analytical expressions derived for the lattice RTP model quantitatively describe the off-lattice AOUP data when mapped via an effective persistence parameter.

Significance
The paper claims to provide the first quantitative foundation for the statistical geometry of active matter. By demonstrating that activity systematically alters the available space for insertion—increasing total volume while enhancing connectivity—the work offers a new geometric perspective on active thermodynamics. The authors posit that these insertion space metrics serve as robust indicators of collective phase behaviors, potentially offering a geometric route to understanding phenomena like MIPS without relying solely on mechanical definitions of pressure or temperature, which are often ambiguous in non-equilibrium systems. The results are presented as experimentally testable using monodisperse active colloids.

More like this