Origin of Persistent Boundary Motion in Confined Active Matter

Technical Summary: Origin of Persistent Boundary Motion in Confined Active Matter

Problem Statement
Active matter systems, characterized by constituents that consume energy to generate directed motion, exhibit distinct behaviors under confinement, including strong boundary affinity and persistent surface motion. While extensive research has characterized the positional dynamics of confined active particles—such as the accumulation near walls and the crossover from exponential to power-law decay in positional distributions—corresponding orientational behavior has received significantly less attention. Existing studies often integrate over particle orientation or focus on average orientational order in stationary states, leaving the role of orientational fluctuations and the stochastic switching dynamics underlying persistent boundary motion largely unexplored. This work addresses the gap in understanding how orientational fluctuations couple with positional accumulation to govern the dynamics of active particles in confined geometries.

Methodology
The authors investigate a single Active Brownian Particle (ABP) confined within a two-dimensional hard circular boundary of radius $R$. The study employs two complementary approaches:

1. Molecular Simulations: Numerical integration of coupled Langevin equations for particle position ($\mathbf{r}$) and orientation ($\theta$). The study focuses on the high-Péclet limit where translational diffusion is vanishingly small ($D_t \to 0$), isolating the effects of activity and confinement. The dynamics are governed by the ratio of persistence length ($l_P = u_0/D_r$) to the confinement radius, defined as the dimensionless confinement strength $\tilde{u} = u_0/(R D_r)$.
2. Analytical Framework: A Fokker-Planck equation is derived for the joint probability density $\Pi(r, \chi; t)$, where $\chi = \theta - \phi$ is the relative orientation between the propulsion direction and the local radial vector. The analysis focuses on the stationary state ($\partial \Pi / \partial t = 0$) and utilizes asymptotic expansions and scaling arguments to derive analytical expressions for positional and orientational distributions.

Key variables include the scaled radial coordinate $z = 1 - r/R$ and the conditional orientational distribution $F(\chi|z)$.

Key Contributions and Results

• Coupling of Position and Orientation: The study establishes that the positional distribution $\phi(z)$ is directly coupled to orientational fluctuations through the conditional orientational distribution. The authors derive an exact expression showing that $\phi(z)$ depends on the variance of $\cos \chi$ (Eq. 17).
• Power-Law Decay of Positional Distribution: In the absence of translational diffusion, the positional distribution exhibits a non-trivial power-law decay near the boundary, $\phi(z) \propto z^{-2/3}$, before crossing over to a uniform bulk profile. This contrasts with previous findings in systems with non-zero translational diffusion, which showed an initial exponential decay.
• Curvature-Induced Bistability: The conditional orientational distribution $F(\chi|z)$ reveals the emergence of two preferred tangential states ($\chi \approx \pm \pi/2$) associated with clockwise and anticlockwise sliding along the boundary. These states are bistable and connected by stochastic flipping pathways.
  • Near the boundary-contact region, the distribution is unimodal and centered at $\chi = 0$ (outward radial alignment).
  • In the near-boundary region ($z > z_0$), the distribution becomes bimodal with peaks at $\chi \approx \pm \pi/2$.
• Scaling Laws:
  • The conditional variance of $\cos \chi$, $\langle \cos^2 \chi \rangle_z$, scales as $(z/\tilde{u})^{2/3}$ near the boundary and saturates to $1/2$ in the bulk.
  • The modal angle shift $\eta = |\chi_p - \pi/2|$ scales as $z^{1/3}$.
  • These scaling behaviors are universal and governed solely by the confinement strength $\tilde{u}$.
• Stochastic Flipping Dynamics: The paper identifies two distinct pathways for switching between the two tangential sliding states:
  1. Rapid boundary-localized flips: Direct transitions between tangential states without traversing the bulk.
  2. Slower bulk-mediated excursions: Transitions involving reorientation in the bulk.
     The mean waiting time $\tau_w$ between flips exhibits a power-law dependence on confinement strength, $D_r \tau_w \propto \tilde{u}^{-3/4}$, for large $\tilde{u}$. This relationship is derived using an effective Kramers escape theory over an orientational potential landscape.

Significance and Claims
The paper claims to provide a fundamental framework for understanding transport, exploration, and escape processes in confined active systems. By demonstrating that the interplay between orientational fluctuations, bistability, positional accumulation, and stochastic switching governs the observed dynamics, the authors argue that:

• The persistent boundary motion is not merely a result of hydrodynamic interactions or simple trapping but is driven by curvature-induced bistable orientational states.
• The stochastic switching between these states is the key dynamical mechanism determining residence times and escape trajectories.
• The derived scaling laws and the coupling between position and orientation offer a predictive tool for analyzing active matter in complex geometries, potentially extending to sorting mechanisms (e.g., ratchets, funnels) and escape kinetics.

The work remains modest in its scope, focusing on a minimal model (single ABP in a hard circle with vanishing translational diffusion) to isolate the mechanism. It suggests that future work could explore the effects of hydrodynamic interactions, three-dimensional geometries, and heterogeneous boundaries, but does not propose specific experimental validations or applications beyond the theoretical insights provided.