Ergodicity Breaking in Active Run-and-Tumble Particles in a Double-Well Potential

Technical Summary: Ergodicity Breaking in Active Run-and-Tumble Particles in a Double-Well Potential

Problem Statement
This study investigates the dynamics of a Run-and-Tumble Particle (RTP) moving in a one-dimensional symmetric double-well potential, $U(x) = \frac{1}{4}(x^2 - a^2)^2$. While the barrier crossing and stationary states of passive Brownian particles in such potentials are well-understood (governed by Kramers' law and a unique Boltzmann distribution), the behavior of active particles in the absence of thermal noise remains largely unexplored. Specifically, the authors address how the intrinsic active nature of RTPs—characterized by persistent motion with speed $v_0$ and orientation switching at rate $\gamma$—affects position fluctuations, ergodicity, and barrier crossing times in landscapes with finite barriers.

Methodology
The authors model the particle's position $x(t)$ using the Langevin equation:
$$\dot{x}(t) = f(x) + v_0\sigma(t)$$
where $f(x) = -U'(x)$ is the deterministic force and $\sigma(t) = \pm 1$ is a dichotomous noise representing the particle's internal orientation, switching at rate $\gamma$.

The analysis proceeds by examining the relationship between the potential force $f(x)$ and the active speed $v_0$. The authors identify a critical barrier height (controlled by the parameter $a$) relative to the active speed. They distinguish two regimes:

1. Ergodic Regime ($a < a_c$): The force magnitude $|f(x)|$ never exceeds $v_0$ within the relevant domain, allowing the particle to access the entire region between turning points.
2. Non-Ergodic Regime ($a > a_c$): The force magnitude exceeds $v_0$ in certain regions, creating "effective walls" where the particle cannot overcome the force regardless of its orientation.

To characterize the stationary state, the authors employ:

• Analytical Derivation: Using known results for RTPs in confining potentials to derive the stationary probability density function (PDF) $P_{st}(x)$.
• Hitting Probabilities: For the non-ergodic regime, they solve the backward Fokker-Planck equation to calculate the "splitting probability" (or hitting probability) $\pi(x_0)$ that a particle starting in the intermediate region $|x_0| < \ell_1$ will eventually get trapped in the left or right well.
• Numerical Simulations: The analytical predictions are validated against numerical integrations of the Langevin equation.

Key Contributions and Results

1. Strong Ergodicity Breaking:
   The paper demonstrates that for barrier heights exceeding a critical threshold ($a > a_c$), the system exhibits strong ergodicity breaking. Unlike passive systems where the long-time distribution is unique and independent of initial conditions, the stationary distribution of the RTP depends crucially on the initial position $x_0$:

   • If $x_0$ lies within the basin of attraction of one well (specifically $|x_0| > \ell_1$), the particle becomes permanently trapped in that well.
   • If $x_0$ lies in the intermediate region ($|x_0| < \ell_1$), the particle can reach either well with a finite probability $\pi(x_0)$. The long-time distribution is a weighted sum of the distributions in the two wells, determined by $\pi(x_0)$.
   • The authors provide an exact analytical expression for $\pi(x_0)$ (Eq. 12) and the resulting non-ergodic stationary distribution (Eq. 13).

2. Critical Barrier Height and Turning Points:
   The critical value is derived as $a_c = \left(\frac{3\sqrt{3}v_0}{2}\right)^{1/3}$.

   • For $a < a_c$, there are two turning points ($\pm \ell$), and the system is ergodic with a unique stationary distribution supported on $[-\ell, \ell]$.
   • For $a > a_c$, the equation $f(x) = \pm v_0$ yields six turning points ($\pm \ell_1, \pm \ell_2, \pm \ell_3$). These points divide the phase space into disconnected regions, preventing the particle from crossing the central barrier if it starts outside the intermediate zone.

3. Barrier Crossing Time and Violation of Kramers' Law:
   In the ergodic regime ($a < a_c$), the authors estimate the characteristic barrier crossing time $\tau$ defined via the stationary density at the barrier top: $\tau = 1/[v_0 P_{st}(0)]$.

   • The study shows that $\tau$ violates the standard Kramers-Arrhenius law ($\tau \sim e^{\beta \Delta V}$).
   • As the barrier height approaches the critical threshold from below ($a \to a_c^-$), the crossing time diverges.
   • The divergence follows a Vogel-Fulcher-Tammann-like form with an anomalous exponent: $\tau \sim (a_c - a)^{-1/2}$.

Significance and Claims
The authors claim that this work provides the first analytical characterization of ergodicity breaking in active particles within a double-well potential without thermal noise. The primary significance lies in demonstrating that active dynamics can fundamentally alter the statistical mechanics of barrier crossing:

• Non-Ergodicity: The existence of a non-ergodic phase where the system fails to explore the full phase space is a direct consequence of the active nature of the dynamics, a feature absent in passive Brownian systems.
• Initial Condition Dependence: The stationary state is not unique but depends on the initial condition, challenging standard assumptions of equilibrium statistical mechanics applied to active matter.
• Anomalous Scaling: The divergence of the barrier crossing time near the critical point with an exponent of $1/2$ highlights a distinct universality class for active barrier crossing compared to thermal activation.

The paper concludes that while finite temperature allows barrier crossing via thermal fluctuations, the "pure active nature" of the dynamics introduces a mechanism for trapping and ergodicity breaking that is qualitatively different from equilibrium phenomena.