Diffusion of gravitactic chiral active Brownian particles in an asymmetric channel

1. Problem Statement

The study addresses the transport and diffusive behavior of chiral active Brownian particles (self-propelled micro- and nanoswimmers) confined within a two-dimensional asymmetric channel (an entropic ratchet). Specifically, the authors investigate how an external gravitational field, which exerts a torque on the particles (gravitaxis), influences their motion.

While the diffusive behavior of active particles in confinement and the effects of external fields have been studied separately, the interplay between chirality (intrinsic rotation), gravitational torque, and geometric confinement in an asymmetric structure remains unexplored. The central question is how these factors combine to create resonant phenomena, rectification, and enhanced diffusion, particularly near the bifurcation point where the external torque matches the particle's intrinsic angular velocity.

2. Methodology

The authors employ a numerical approach based on overdamped Langevin dynamics to model the system.

• Physical Model:
  • Particles: Chiral active particles with self-propelled velocity $V_0$, intrinsic angular velocity $\Omega_0$, and subject to a gravitational torque proportional to $\Omega \sin \theta$.
  • Environment: A 2D asymmetric channel with periodicity $L$. The channel width varies linearly from a minimum $w_{min}$ to a maximum $w_{max}$, defined by an aspect ratio $\epsilon = w_{min}/w_{max}$.
  • Boundary Conditions: Sliding-reflecting conditions are applied at the walls. Upon collision, the normal component of the velocity is reversed (elastic reflection), while the tangential component is preserved, allowing particles to slide along the walls until reoriented by fluctuations.
• Governing Equations:
  The dynamics are described by coupled Langevin equations for position $\mathbf{r}$ and orientation angle $\theta$:
  $$\dot{\mathbf{r}} = V_0 \hat{n} + \sqrt{2D_t} \boldsymbol{\xi}(t)$$
  $$\dot{\theta} = \Omega_0 - \Omega \sin \theta + \sqrt{2D_r} \zeta(t)$$
  Where $D_t$ and $D_r$ are translational and rotational diffusion constants, and $\xi, \zeta$ are Gaussian white noises.
• Dimensionless Analysis:
  The system is rescaled using characteristic time $\tau = L^2/D_t$ and length $L$. Key dimensionless parameters include:
  • $v_0$: Dimensionless self-propulsion speed.
  • $\omega_0$: Dimensionless intrinsic angular velocity.
  • $\omega$: Dimensionless external torque (gravitational).
  • $\alpha$: Dimensionless rotational diffusion rate.
  • $\epsilon$: Channel aspect ratio.
• Simulation:
  The authors solve the equations numerically using the stochastic Euler algorithm. Transport characteristics are calculated by averaging over $10^4$ stochastic trajectories in the steady state.
  • Average Velocity ($v$): $\lim_{t\to\infty} \langle x(t) \rangle / t$
  • Effective Diffusion Coefficient ($D_{eff}$): $\lim_{t\to\infty} (\langle x^2(t) \rangle - \langle x(t) \rangle^2) / 2t$

3. Key Contributions

1. Discovery of Resonant Diffusion: The paper identifies a specific regime where the external torque $\omega$ approaches the intrinsic angular velocity $\omega_0$ ($\omega \approx \omega_0$), leading to a phenomenon termed resonant diffusion.
2. Mechanism Elucidation: The authors map the angular dynamics to an overdamped driven pendulum. They demonstrate that the resonance arises from a saddle-node bifurcation in this pendulum system. At the critical point $\omega = \omega_0$, the effective potential softens, allowing for large fluctuations that enhance diffusion.
3. Role of Rotational Diffusion: The study clarifies that resonant diffusion is a low-noise phenomenon. It is prominent only for low rotational diffusion rates ($\alpha$). As $\alpha$ increases, the resonance vanishes, and the system transitions to passive Brownian behavior.
4. Geometric Dependence: The work quantifies how the channel's aspect ratio ($\epsilon$) modulates transport, showing that wider bottlenecks enhance both the average velocity and the magnitude of the diffusion peak.

4. Key Results

A. Steady-State Spatial Distribution

• No Torque ($\omega = 0$): Particles accumulate inhomogeneously, primarily near the left corners of the channel cells due to chirality and geometric asymmetry.
• Resonance ($\omega \approx \omega_0$): A pronounced accumulation of particles occurs near the upper-left corner of the channel. This is caused by the emergence of a stable fixed point in the orientation angle ($\theta^* \approx \pi/2$) combined with the channel geometry.
• High Torque ($\omega \gg \omega_0$): Particles align horizontally and accumulate near the channel center and exit, leading to a symmetric distribution.
• Effect of $\alpha$: At the resonance point, low $\alpha$ yields strong upper-left accumulation. As $\alpha$ increases, particles also accumulate in the lower-left corner, eventually symmetrizing the distribution and destroying the resonance effect.

B. Transport Characteristics

• Average Velocity ($v$):
  • Rectification: For $\omega \ll \omega_0$, particles exhibit directed motion (rectification) in the positive $x$-direction due to the channel shape.
  • Saturation: For $\omega > \omega_0$ (locked state), $v$ increases monotonically with $\omega$ and saturates at a value determined by the channel aspect ratio.
  • Non-monotonicity: Near the bifurcation ($\omega \approx \omega_0$), $v$ shows a non-monotonic dependence on the rotational diffusion rate $\alpha$, peaking at an optimal $\alpha$.
• Effective Diffusion Coefficient ($D_{eff}$):
  • Resonant Peak: $D_{eff}$ exhibits a sharp, enhanced peak exactly at the bifurcation point ($\omega/\omega_0 = 1$).
  • Dependence on Parameters:
    • The peak is most pronounced for low $\omega_0$ and low $\alpha$.
    • The peak height increases with the channel aspect ratio ($\epsilon$).
    • For high $\alpha$ or high $\omega_0$, the peak disappears.
  • Physical Origin: The enhancement is driven by the "softening" of the effective potential $U(\theta) = -\omega_0 \theta - \omega \cos \theta$ near the bifurcation, which allows for larger orientational fluctuations and particle accumulation in high-concentration regions.

5. Significance and Applications

The findings provide a fundamental understanding of how external fields and geometric confinement can be tuned to control active matter transport.

• Microfluidics and Lab-on-a-Chip: The resonant diffusion and rectification effects offer new strategies for particle separation and sorting based on chirality or rotational diffusion rates without external chemical gradients.
• Targeted Drug Delivery: Understanding how gravitational or magnetic torques affect the accumulation of microswimmers in confined biological structures (e.g., capillaries) could improve the efficiency of targeted delivery systems.
• Active Materials: The results aid in the design of advanced active materials where transport properties can be dynamically switched by tuning external fields relative to intrinsic particle properties.
• Theoretical Insight: The work bridges the gap between the physics of driven pendulums (bifurcation theory) and active matter, providing a minimal model for resonant transport in entropic ratchets.

In summary, the paper demonstrates that resonant diffusion is a robust phenomenon in chiral active systems confined in asymmetric channels, occurring specifically when the external torque matches the intrinsic rotation, provided that rotational noise is low. This resonance leads to significant particle accumulation and enhanced transport, controllable via channel geometry and field strength.