🔬 condensed matter

Splitting probabilities of confined active particles

This paper employs the backward Fokker-Planck equation to analytically and numerically investigate how particle activity, chirality, and channel geometry influence the splitting probabilities and escape dynamics of active particles in one-dimensional intervals and two-dimensional corrugated channels.

Original authors: Sarafa A. Iyaniwura, Zhiwei Peng

Published 2026-03-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sarafa A. Iyaniwura, Zhiwei Peng

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 hallway filled with people trying to get to one of two exits: the left door or the right door.

In the world of physics, we usually study passive particles (like dust motes or a drunk person stumbling). They move randomly, bumping into walls and each other, with no real direction. If you drop a passive particle in the middle of the hall, its chance of exiting left or right is purely based on how close it is to each door. It's a simple, linear game of chance.

But this paper is about active particles. Think of these as people who are running with a purpose. They have their own internal energy (like a battery or a motor) that makes them swim or walk in a specific direction. They don't just stumble; they propel themselves.

The authors of this paper wanted to answer a simple question: If these "self-propelled" runners are trapped in a hallway, which door are they more likely to escape through, and how does the shape of the hallway change their odds?

Here is a breakdown of their findings using everyday analogies:

1. The Straight Hallway (The 1D Interval)

Imagine a long, straight corridor with a left exit and a right exit.

The Passive Case: If you drop a passive particle in the middle, it has a 50/50 chance. If you drop it closer to the right, it's more likely to hit the right door. It's a straight line.
The Active Case (Running): Now, imagine the particles are runners.
- If they are slow runners: They still mostly follow the "closer door" rule, but their running adds a little twist. If they start on the left but are facing right, they might run past the middle and escape the right door, even though they started closer to the left.
- If they are fast runners: This is where it gets wild. If the runners are very fast and energetic, it doesn't matter where they start in the middle of the hall. Because they are zooming so fast, they are equally likely to hit the left or right door eventually. The only time their starting position matters is if they are right next to a door; then, they just run straight out.
- The "Spin" Factor (Chirality): Some of these runners are also spinning in circles (like a dog chasing its tail). If they spin very fast, they stop moving in a straight line and start wobbling in place. The paper found that if they spin fast enough, they act just like the passive, stumbling dust motes again. Their "running" power is canceled out by their spinning.

2. The Wavy Hallway (The Corrugated Channel)

Now, imagine the hallway isn't straight. The walls are wavy, like a snake sliding through a tunnel. Sometimes the tunnel is wide, and sometimes it gets very narrow.

The Bottleneck Effect: In a wavy hallway, the narrow parts act like bottlenecks.
- If the right exit is located in a narrow, tight squeeze, it becomes much harder to get out that way. Even if a runner is heading right, the narrow walls might bounce them back or slow them down.
- The paper found that the shape of the walls creates an "entropic barrier." It's like trying to push a large suitcase through a narrow door; the geometry itself fights against the particle.
The "Flat" Approximation: When the hallway is very thin (like a long, narrow tube), the authors found a clever shortcut (called the Fick-Jacobs reduction). They realized they could ignore the side-to-side wiggles and just look at the forward motion. It's like saying, "If the hallway is thin enough, the particle is just moving forward, and the walls don't matter much."
When the Shortcut Fails: However, if the hallway is wide and wavy, that shortcut stops working. The particles get stuck in the wide parts or bounce off the walls in complex ways. The "running" particles don't just move forward; they get trapped in the curves, and their path becomes much more complicated than the simple math predicted.

Why Does This Matter?

You might wonder, "Who cares about math runners in a hallway?"

This research is actually a map for understanding real life:

Biology: Inside your body, cells and bacteria are constantly moving. Some are passive, but many (like sperm or immune cells) are active runners. Understanding how they navigate narrow blood vessels or crowded tissues helps us understand how they find targets (like a virus or a tumor).
Technology: Scientists are building tiny "micro-machines" that swim through fluids to deliver medicine. This paper helps them design better channels and devices so these tiny robots don't get stuck or go the wrong way.

The Big Takeaway

The paper teaches us that activity changes the rules of the game.

Passive particles are at the mercy of random chance and distance.
Active particles are at the mercy of their own energy, their direction, and how much they spin.
The environment (the shape of the hallway) acts as a filter, amplifying or canceling out their ability to escape.

In short: If you want to know where a particle will go, you can't just look at where it started. You have to ask: Is it running? Is it spinning? And is the hallway narrow or wide?

1. Problem Statement

The paper addresses the fundamental problem of splitting probabilities (also known as exit or escape probabilities) for active particles confined within bounded domains. Unlike passive Brownian particles, active particles (such as bacteria, synthetic microswimmers, or cells) possess internal energy sources that drive self-propulsion and directional persistence. This non-equilibrium behavior fundamentally alters transport statistics, particularly in confined or structured environments like biological microenvironments and microfluidic devices.

The specific challenges addressed include:

How self-propulsion (activity) and chirality (intrinsic rotation) modify the likelihood of a particle escaping through a specific boundary before another.
The interplay between geometric confinement (e.g., corrugated channels) and active motility mechanisms.
The lack of comprehensive analytical frameworks for active particles in structured geometries compared to the well-established theory for passive Brownian motion.

2. Methodology

The authors employ the backward Fokker–Planck equation (BFPE) to model the splitting probability. This approach allows for the calculation of the probability of reaching a specific boundary starting from a given position and orientation, without needing to simulate full time trajectories.

Key Methodological Components:

Governing Equation: The BFPE for a chiral active particle (CAP) is derived, incorporating translational diffusion ( $D_x$ ), rotational diffusion ( $D_\theta$ ), swim speed ( $U_s$ ), and chirality ( $\Omega$ ).
Dimensionless Parameters: The system is characterized by:
- Péclet number ( $Pe = U_s L / D_x$ ): Ratio of active to diffusive transport.
- Dimensionless chirality ( $\chi = \Omega L^2 / D_x$ ): Ratio of rotational drift to diffusion.
- Rotational diffusivity parameter ( $\gamma = L^2 D_\theta / D_x$ ).
Geometries Studied:
1. One-Dimensional (1D) Interval: $[-L, L]$ with absorbing boundaries.
2. Two-Dimensional (2D) Corrugated Channel: A channel with periodic width modulation $H(x) = H_0(1 + \alpha \cos(\pi x/L))$ .
Analytical Techniques:
- Asymptotic Analysis: Perturbation expansions for weak activity ( $Pe \ll 1$ ) and singular perturbation theory for strong activity ( $Pe \gg 1$ ) to derive boundary layer structures.
- Fick–Jacobs Reduction: A dimensional reduction technique applied to the 2D channel in the small aspect ratio limit ( $\epsilon \ll 1$ ) to derive effective 1D transport equations.
Numerical Methods: The full backward Fokker–Planck equation is solved numerically using the Finite Element Method (FEM) implemented in FreeFem++ for finite aspect ratios and arbitrary activity levels.

3. Key Contributions

Analytical Solutions for 1D Active Particles: The authors derived exact analytical expressions for splitting probabilities in the weak-activity regime (up to $O(Pe^2)$ ) and leading-order asymptotic solutions for the strong-activity regime ( $Pe \to \infty$ ).
Characterization of Chirality Effects: The study quantifies how intrinsic rotation (chirality) competes with translational persistence, effectively acting as a "diffusive" mechanism that suppresses activity-induced biases.
Extension to Corrugated Channels: The work extends the splitting probability framework to 2D structured domains, validating the Fick–Jacobs approximation for narrow channels and identifying its limitations in wider channels.
Boundary Layer Analysis: A rigorous characterization of boundary layers in high-activity regimes, showing how the splitting probability transitions from a linear profile (passive) to a uniform bulk profile with sharp boundary adjustments.

4. Key Results

A. One-Dimensional Interval

Weak Activity ( $Pe \ll 1$ ): The splitting probability $\langle p_R \rangle$ $⟨ p_{R} ⟩$ deviates from the linear passive profile ( $p_0 = (x+1)/2$ $p_{0} = (x + 1) /2$ ) at order $O(Pe^2)$ $O (P e^{2})$ .
- For particles starting in the left half ( $x<0$ ), activity increases the probability of escaping the right boundary.
- For particles starting in the right half ( $x>0$ ), activity decreases the probability of escaping the right boundary (due to the tendency to swim back toward the center or left).
- Chirality Effect: As chirality ( $\chi$ ) increases, the magnitude of this correction decreases. In the limit $\chi \to \infty$ , the motion becomes effectively diffusive, and the splitting probability reverts to the passive linear profile.
Strong Activity ( $Pe \gg 1$ ):
- The splitting probability becomes nearly uniform ( $\approx 0.5$ ) across the bulk of the domain, regardless of the starting position.
- Boundary Layers: Sharp boundary layers form near the exits ( $x = \pm 1$ ) where the probability rapidly adjusts to satisfy boundary conditions ($0$ or $1$). The thickness of these layers scales with $1/Pe$.
- The analytical asymptotic solution matches numerical FEM results with high precision.

B. Corrugated Channels (2D)

Fick–Jacobs Approximation: For small aspect ratios ( $\epsilon \ll 1$ ), the 2D problem reduces to an effective 1D equation. This approximation accurately predicts splitting probabilities and captures the effect of geometric confinement.
Geometric Bias: In corrugated channels, geometric confinement creates entropic barriers.
- Particles starting on the side of the "narrow" exit face higher resistance, reducing their escape probability compared to a flat channel.
- The splitting probability is biased toward the wider exit, even for active particles.
Finite Aspect Ratios: As the aspect ratio ( $\epsilon$ $ϵ$ ) increases, the Fick–Jacobs approximation begins to deviate from the full solution.
- Transverse Non-uniformity: In wider channels, the splitting probability is no longer uniform across the transverse direction ( $y$ ). It exhibits maxima near the channel center ( $y=0$ ) because active particles with persistent motion do not fully equilibrate transversely before escaping.
Chirality in 2D: Similar to the 1D case, chirality suppresses the influence of activity, making the splitting profile more linear and less sensitive to geometric modulation.

5. Significance and Implications

Theoretical Framework: The paper provides a robust predictive framework for first-passage phenomena in active matter, bridging the gap between classical diffusion theory and non-equilibrium active transport.
Biological Relevance: The findings are crucial for understanding intracellular transport, immune cell surveillance, and bacterial chemotaxis in complex, confined environments (e.g., tissue pores, microfluidic traps).
Engineering Applications: The results offer design principles for microfluidic devices intended to sort or separate active particles based on their motility (speed, chirality) or to control the release of active cargo.
Limitations and Future Work: The current model assumes idealized point particles without hydrodynamic interactions or crowding. The authors suggest future extensions to include time-dependent boundaries, particle-particle interactions, and partially reactive boundaries to better model dense biological systems.

In summary, this work demonstrates that activity and geometry are coupled in determining escape statistics. While passive diffusion is governed solely by geometry and distance, active transport is a complex interplay of self-propulsion, rotational dynamics, and confinement, where chirality serves as a critical regulator that can tune the system between ballistic and diffusive regimes.