General perturbative framework for kinetics of rare transitions in 1-dimensional active particle systems
Pith reviewed 2026-05-10 07:46 UTC · model grok-4.3
The pith
A projection-operator approach yields an analytic transition rate for active particles valid across all persistence times.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using a projection-operator formalism, perturbative rates are obtained by integrating out the fast variable in each asymptotic regime: the activity for small persistence times and the position for large persistence times. These two expansions are combined into a rational approximation that supplies an analytic rate expression valid for all persistence times and activity strengths in the rare-event limit, in excellent agreement with simulations. The framework applies to a broad class of driven systems, with the thermal AOUP appearing as a special case.
What carries the argument
Projection-operator formalism that isolates the slow variable in each timescale limit and produces a rational approximation from the two resulting perturbative expansions.
If this is right
- The rate formula reproduces known results for small persistence times by integrating out the activity.
- New analytic rates are obtained for large persistence times by integrating out the position instead.
- The approximation supplies a parameter-free expression that holds for any activity strength in the rare-event regime.
- The same construction applies to general driven systems beyond the specific AOUP model.
Where Pith is reading between the lines
- The rational-approximation strategy may extend to other driven systems that possess two well-separated timescales.
- If the form remains accurate, it could be used to predict how activity strength and persistence together control escape rates without running separate simulations for each parameter set.
Load-bearing premise
The two asymptotic expansions can be combined into a rational approximation that remains accurate at intermediate persistence times without additional fitting or higher-order corrections.
What would settle it
Direct numerical simulations of the particle dynamics that show large deviations from the rational approximation at moderate persistence times would falsify the central claim.
Figures
read the original abstract
We present a theoretical framework that enables investigating rare transitions in a general model of an active particle in an external potential, with the thermal Active Ornstein-Uhlenbeck Particle (AOUP) appearing as a special case. Using a projection-operator formalism, we compute transition rates perturbatively in two distinct asymptotic regimes. In the regime of small persistence times-where the activity evolves much faster than the particle's position-integrating out the activity reproduces the rates previously reported in the literature. In the opposite regime of large persistence times, we instead integrate out the position and obtain the corresponding rates analytically. Together, these asymptotic expansions uniquely specify a rational approximation that remains accurate across intermediate persistence times. As a result, we obtain an analytic expression for the rate valid across all persistence times and activity strengths in the rare-event limit, which are in excellent agreement with numerical simulations. The presented framework applies to rare transitions in a broad class of driven systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a projection-operator formalism for perturbative computation of rare transition rates in a general class of one-dimensional active particles in an external potential (with AOUP as a special case). Small-persistence rates are obtained by integrating out activity and recover prior literature; large-persistence rates are obtained by integrating out position. These two leading-order expansions are asserted to uniquely determine a rational approximant that remains accurate at all intermediate persistence times in the rare-event limit, yielding an analytic rate expression that agrees excellently with simulations and applies broadly to driven systems.
Significance. If the rational interpolation is rigorously controlled, the work would supply a useful analytic tool for rare-event kinetics in active matter across the full persistence range, extending beyond the usual asymptotic limits. The systematic projection-operator approach and recovery of known small-tau results are strengths; the framework's generality to other driven systems is also positive.
major comments (2)
- [Abstract] Abstract: the central claim that the small- and large-persistence expansions 'uniquely specify a rational approximation that remains accurate across intermediate persistence times' is load-bearing yet unsupported by a controlled error bound, remainder estimate, or demonstration that the projection-operator kernels produce a meromorphic structure free of non-analytic contributions at crossover. Without this, systematic deviation remains possible where neither asymptotic applies.
- [Numerical validation] The numerical validation section: 'excellent agreement with numerical simulations' is stated without quantitative error metrics, the precise range of persistence times and activity strengths tested, or confirmation that the rational form was not adjusted post-hoc; this leaves the accuracy claim unverifiable from the given information.
minor comments (2)
- The explicit algebraic form of the rational approximant should be written out in the main text (rather than described) so that readers can reproduce the expression without reconstructing it from the asymptotics.
- Notation for the projection operators and the resulting kernels could be introduced more explicitly at first use to improve readability for readers outside the immediate subfield.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable feedback on our manuscript. We address the major comments below and will revise the manuscript to incorporate improvements where appropriate.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the small- and large-persistence expansions 'uniquely specify a rational approximation that remains accurate across intermediate persistence times' is load-bearing yet unsupported by a controlled error bound, remainder estimate, or demonstration that the projection-operator kernels produce a meromorphic structure free of non-analytic contributions at crossover. Without this, systematic deviation remains possible where neither asymptotic applies.
Authors: We acknowledge that the manuscript does not provide a rigorous controlled error bound or a proof that the kernels lead to a meromorphic structure without non-analytic terms at the crossover. The rational approximation is constructed to match the leading-order asymptotics from the projection-operator formalism in both limits, and its form is motivated by the structure of the effective rate equations. While we cannot rule out systematic deviations in principle, extensive numerical tests (detailed in the manuscript) show close agreement. In the revision, we will add a section discussing the approximation's limitations and the absence of a full error estimate, and we will moderate the claim in the abstract to reflect that it is an interpolation supported by asymptotics and numerics rather than a rigorously controlled result. revision: yes
-
Referee: [Numerical validation] The numerical validation section: 'excellent agreement with numerical simulations' is stated without quantitative error metrics, the precise range of persistence times and activity strengths tested, or confirmation that the rational form was not adjusted post-hoc; this leaves the accuracy claim unverifiable from the given information.
Authors: The referee is correct that the current presentation lacks quantitative metrics and precise parameter ranges. The rational approximant is uniquely determined by matching the two asymptotic expansions and was not fitted to intermediate data. In the revised version, we will include a table summarizing the relative errors for various persistence times and activity strengths, specify the tested ranges (e.g., persistence times from 10^{-2} to 10^2 and activity parameters up to the rare-event regime), and add a statement confirming the a priori determination of the approximant. This will make the validation section more transparent and verifiable. revision: yes
Circularity Check
No significant circularity; asymptotics independently derived and combined via interpolation validated externally
full rationale
The paper derives small-persistence rates via projection-operator integration of activity, reproducing external prior literature. Large-persistence rates follow from integrating out position in the same formalism. The rational approximation is explicitly constructed to match the two leading-order asymptotic expansions at their respective limits. Because both expansions are obtained from the projection-operator framework (one matching known results, one newly analytic) and the final expression is cross-checked against independent numerical simulations rather than being tautological with the inputs, the derivation chain does not reduce to self-definition or fitted inputs by construction. No load-bearing self-citations or uniqueness theorems from the authors' prior work are invoked.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Projection-operator formalism can be applied to separate fast and slow variables in the active-particle dynamics.
- domain assumption The rare-event limit allows a well-defined perturbative expansion for transition rates.
Reference graph
Works this paper leans on
- [1]
-
[2]
S. Ramaswamy, The mechanics and statistics of active matter, Annual Review of Condensed Matter Physics1, 323 (2010)
work page 2010
-
[3]
M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Hydrody- namics of soft active matter, Reviews of Modern Physics 85, 1143 (2013)
work page 2013
-
[4]
C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe, Active particles in complex and crowded environments, Reviews of Modern Physics88, 045006 (2016)
work page 2016
-
[5]
D. Needleman and Z. Dogic, Active matter at the inter- face between biology and physics, Nature Reviews Mate- rials2, 17048 (2017)
work page 2017
- [6]
-
[7]
M. E. Cates and J. Tailleur, Motility-induced phase sep- aration, Annual Review of Condensed Matter Physics6, 219 (2015)
work page 2015
-
[8]
P. H¨ anggi, P. Talkner, and M. Borkovec, Reaction-rate theory: fifty years after kramers, Rev. Mod. Phys.62, 251 (1990)
work page 1990
-
[9]
E. Woillez, Y. Kafri, and V. Lecomte, Nonlocal stationary probability distributions and escape rates for an active ornstein–uhlenbeck particle, Journal of Statistical Me- chanics: Theory and Experiment2020, 063204 (2020)
work page 2020
-
[10]
D. Martin and T. Arnoulx de Pirey, Aoup in the pres- ence of brownian noise: a perturbative approach, Journal of Statistical Mechanics: Theory and Experiment2021, 043205 (2021)
work page 2021
-
[11]
A. Militaru, M. Innerbichler, M. Frimmer, F. Tebbenjo- hanns, L. Novotny, and C. Dellago, Escape dynamics of active particles in multistable potentials, Nature Com- munications12, 10.1038/s41467-021-22647-6 (2021)
-
[12]
L. Caprini, F. Cecconi, and U. Marini Bettolo Marconi, Correlated escape of active particles across a poten- tial barrier, The Journal of Chemical Physics155, 234902 (2021), https://pubs.aip.org/aip/jcp/article- pdf/doi/10.1063/5.0074072/13625527/234902 1 online.pdf
-
[13]
J. O’Byrne, A. Solon, J. Tailleur, and Y. Zhao, An intro- duction to motility-induced phase separation, inOut-of- equilibrium Soft Matter, edited by C. Kurzthaler, L. Gen- tile, and H. A. Stone (The Royal Society of Chemistry, 2023)
work page 2023
-
[14]
L. Angelani, R. Leonardo, and M. Paoluzzi, First-passage time of run-and-tumble particles, The European physical journal. E, Soft matter37, 15 (2014)
work page 2014
-
[15]
E. Woillez, Y. Zhao, Y. Kafri, V. Lecomte, and J. Tailleur, Activated escape of a self-propelled particle from a metastable state, Phys. Rev. Lett.122, 258001 (2019)
work page 2019
-
[16]
V. Bach, J. Fr¨ ohlich, and I. M. Sigal, Quantum electro- dynamics of confined nonrelativistic particles, Advances in Mathematics137, 299 (1998)
work page 1998
-
[17]
C. J. Pethick and H. Smith,Bose–Einstein Condensation in Dilute Gases, 2nd ed. (Cambridge University Press, 2008)
work page 2008
-
[18]
P. Ashwin and A. S. von der Heydt, Extreme sensitivity and climate tipping points, Journal of Statistical Physics 179, 1531 (2020)
work page 2020
-
[19]
G. V. Vinci and M. Mattia, Escape time in bistable neu- ronal populations driven by colored synaptic noise, Phys. Rev. Res.7, 023172 (2025)
work page 2025
-
[20]
C. W. Gardiner,Handbook of stochastic methods for physics, chemistry and the natural sciences, 3rd ed., Springer Series in Synergetics, Vol. 13 (Springer-Verlag, Berlin, 2004) pp. xviii+415
work page 2004
- [21]
-
[22]
J. F. Mejias, H. J. Kappen, and J. J. Torres, Irregular dynamics in up and down cortical states, PLOS ONE5, 1 (2010)
work page 2010
- [23]
-
[24]
R. F. Fox, Mean first-passage times and colored noise, Phys. Rev. A37, 911 (1988). 7 NON-DIMENSIONALIZATION We start with the dimensional version of the generalized active particle Langevin equations: dx∗ dt∗ =− k∗ γ F(x ∗)− λ∗ γ Fa(x∗, a) + s 2kBT γ ηx(t∗),(14) da dt∗ =− 1 τ ∗ v′(a) + r 2 τ ∗ ηa(t∗),(15) where⟨η x(t∗ 1)ηx(t∗ 2)⟩=δ(t ∗ 1 −t ∗
work page 1988
-
[25]
We have chosenato be dimensionless, such that the non-dimensionalization will not affect it
and⟨η a(t∗ 1)ηa(t∗ 2)⟩=δ(t ∗ 1 −t ∗ 2). We have chosenato be dimensionless, such that the non-dimensionalization will not affect it. The corresponding Fokker-Planck equation takes the form ∂t∗ p∗ =∂ x∗ k∗ γ F(x ∗) + λ∗ γ Fa(x∗, a) p∗ + kBT γ ∂2 x∗ p∗ + 1 τ ∗ ∂a[v′(a)p∗] + 1 τ ∗ ∂2 ap∗.(16) Here,p ∗ denotes the dimensionful density. To render the equations...
-
[26]
terms. We can write theA ′ term in terms ofAthrough integration by parts Z ∞ −∞ 2A′e−v(a)da= 2Ae −v(a) ∞ −∞| {z } =0 − Z ∞ −∞ 2A(a)[−v′(a)e−v(a)]da,(66) such that we finally obtain r≈ Z ∞ −∞ 1 Za e−v(a) R0(a) + 1 τ v′(a)A(a)−B(a) da,(67) as mentioned in the main text in Eqs. (10) and (11). We also could have obtained Eq. (67) by simply holding that 1 τ La...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.