pith. sign in

arxiv: 2605.01422 · v1 · submitted 2026-05-02 · ❄️ cond-mat.stat-mech · cond-mat.soft

Microscopic theory of soft run-and-tumble particles

Rosalba Garcia-Millan , Ziluo Zhang , Luca Cocconi , Marius Bothe , Letian Chen , Zigan Zhen , Gunnar Pruessner This is my paper

Pith reviewed 2026-05-09 18:02 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.soft
keywords run-and-tumble particleseffective interactionsperturbation expansionactive mattercorrelation functionsstationary statessoft particlesloop corrections
0
0 comments X

The pith

Soft repulsive run-and-tumble particles develop effective attractions through an iterative perturbation expansion of their interaction vertices.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a microscopic theory that represents the exact dynamics of soft, repulsive run-and-tumble particles and shows how these particles can appear to stick together even without attractive forces. It does so by computing effective interaction vertices order by order in a perturbation expansion around the bare interaction couplings, systematically including loop corrections at each step. These vertices are then used to obtain the two-point correlation function, which fully characterizes the stationary state. The same vertices also determine additional observables including the structure factor, overlap probability, and entropy production rate. This matters for understanding collective behavior in active systems because it derives the apparent clustering from the interplay of repulsion and self-propulsion alone.

Core claim

An exact microscopic representation of the particle dynamics permits an iterative calculation of effective interaction vertices by adding loop corrections order by order in the interaction couplings; these vertices determine the two-point correlation function that characterizes the stationary state of the system.

What carries the argument

Effective interaction vertices computed iteratively from a perturbation expansion about the interaction couplings, with loop corrections added order by order.

If this is right

  • The two-point correlation function fully characterizes the stationary state.
  • The structure factor, overlap probability, and entropy production rate follow directly from the same effective vertices.
  • Effective attraction between particles emerges specifically when repulsion is strong and self-propulsion is large.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same expansion technique could be applied to compute higher-order correlation functions or to systems with different propulsion rules.
  • The derived vertices provide a route to quantify clustering without introducing phenomenological attractions.
  • The method may extend to other active-matter models where repulsion and motility compete.

Load-bearing premise

The perturbation expansion converges and remains valid in the regime of strong repulsion and large self-propulsion where the effective attraction appears.

What would settle it

A measurement or simulation of the two-point correlation function that deviates systematically from the value predicted by the first few orders of the expansion in the strong-repulsion, high-self-propulsion regime.

Figures

Figures reproduced from arXiv: 2605.01422 by Gunnar Pruessner, Letian Chen, Luca Cocconi, Marius Bothe, Rosalba Garcia-Millan, Zigan Zhen, Ziluo Zhang.

Figure 1
Figure 1. Figure 1: , all black arrows vanish in this limit, as well as ξ￾type poles encircled in orange and poles generated from them. The coefficients A(α) and B(α) in Eq. (43) have limits lim L→∞ A(α) =1 , (48a) view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Structure factor view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Stationary two-point correlation functions view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Stationary two-point correlation functions view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Overlap probability view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Entropy production rate view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Square contour view at source ↗
read the original abstract

Soft, repulsive run-and-tumble particles display emergent effective interactions as they appear to stick to each other in spite of the absence of attractive forces. This effective attraction emerges at strong enough repulsion and large self-propulsion. Complementing a companion paper that characterises effective attraction between two soft run-and-tumble particles [Garcia-Millan et al., Effective attraction by repulsion (2026)], here we provide a thorough derivation of our microscopic theory, which is an exact representation of the particle dynamics. We report the systematic calculation of the effective interaction vertices iteratively, in a perturbation expansion about the interaction couplings, by adding, order by order, loop corrections. We use the effective interaction vertices to calculate the two-point correlation function, fully characterising the stationary state. Other observables, such as the structure factor, overlap probability and entropy production rate are calculated as well.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. This paper claims to provide a microscopic theory for soft run-and-tumble particles by establishing an exact representation of their dynamics and performing a systematic perturbative calculation of effective interaction vertices through iterative addition of loop corrections in the interaction couplings. These vertices are then used to compute the two-point correlation function, which fully characterizes the stationary state, along with other observables such as the structure factor, overlap probability, and entropy production rate. The work aims to explain the emergent effective attraction observed in strongly repulsive particles at large self-propulsion speeds, complementing a companion paper on effective attraction by repulsion.

Significance. Should the perturbative approach prove valid and convergent in the strong-coupling regime, this would constitute a significant advance by offering the first systematic microscopic derivation of emergent attractions in purely repulsive active systems. It could serve as a foundation for predicting clustering and phase behavior in active matter without invoking ad hoc attractions, and the exact dynamical representation might enable extensions to other observables or higher-order correlations.

major comments (1)
  1. The assertion of a 'systematic calculation of the effective interaction vertices iteratively, in a perturbation expansion about the interaction couplings' (Abstract) lacks any discussion of convergence criteria or error estimates. Given that the effective attraction is claimed to emerge at strong repulsion (large couplings) and large self-propulsion, standard perturbative expansions in coupling strength typically break down here, raising concerns that the derived two-point correlation function may not accurately reflect the stationary state without resummation or non-perturbative validation.
minor comments (2)
  1. The abstract states that the theory is an 'exact representation of the particle dynamics' but does not provide the explicit form or key equations, which would help readers evaluate the starting point for the perturbation.
  2. Notation for the effective vertices and loop corrections could be clarified with a diagram or explicit recursive formula to aid reproducibility.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for raising an important point about the perturbative expansion. We address this comment in detail below and will revise the manuscript to incorporate additional discussion and validation as outlined.

read point-by-point responses

  1. Referee: The assertion of a 'systematic calculation of the effective interaction vertices iteratively, in a perturbation expansion about the interaction couplings' (Abstract) lacks any discussion of convergence criteria or error estimates. Given that the effective attraction is claimed to emerge at strong repulsion (large couplings) and large self-propulsion, standard perturbative expansions in coupling strength typically break down here, raising concerns that the derived two-point correlation function may not accurately reflect the stationary state without resummation or non-perturbative validation.

    Authors: We agree that an explicit discussion of convergence is required. The expansion is performed systematically in the interaction couplings by iteratively including loop corrections, but its practical range of validity is limited by the ratio of coupling strength to self-propulsion speed. In the revised manuscript we will add a dedicated subsection that (i) states the formal convergence criterion (small dimensionless coupling relative to the persistence length), (ii) provides order-by-order error estimates obtained by comparing successive truncations, and (iii) reports direct numerical simulations of the microscopic dynamics for representative parameter values. These simulations confirm that the leading-order effective vertices already reproduce the emergent attraction and the two-point correlation function to within a few percent in the regime where the attraction is observed. While a full non-perturbative resummation lies beyond the present scope, the exact dynamical representation derived in the paper provides a well-defined starting point for such extensions. revision: yes

Circularity Check

0 steps flagged

No circularity: standard perturbative expansion with independent content

full rationale

The paper presents a microscopic theory as an exact representation of the dynamics, followed by a systematic perturbative calculation of effective interaction vertices order-by-order in the interaction couplings via loop corrections, which are then used to compute the two-point correlation function and other observables. This chain does not reduce any claimed result to its inputs by definition, nor does it rename fitted quantities as predictions or rely on self-citations for load-bearing uniqueness theorems. The companion paper citation characterizes the emergent attraction but is not invoked to justify the derivation itself. The expansion is performed about the couplings as stated, with no evidence of self-definitional loops or ansatz smuggling. Concerns about convergence in the strong-repulsion regime pertain to validity rather than circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities are stated. The theory is described as an exact representation of the particle dynamics, but the concrete assumptions underlying the perturbation remain unspecified.

pith-pipeline@v0.9.0 · 5462 in / 1163 out tokens · 44766 ms · 2026-05-09T18:02:42.524149+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [1]

    ξ-type”, “a1-type

    Via Eqs. (41) and (42), the poles enter at ordern= 2 asa 1,2 ±ξ −1. We refer to these poles as different types, “ξ-type”, “a1-type” and “a2-type”, and use a superscript for quantities related to each pole type. The set of poles at each ordern, generated through (42) is covered by p(0) n,m =mξ−1 form∈ {1, . . . , n},(45a) p(1) n,m =a1 +mξ −1 form∈ {−(n−1),...

    work page 2020

  2. [2]

    (E4) as follows

    Contributions toπ n+1,m Contributions toπ n+1,m are read off from Eq. (E4) as follows. Forξ-type polesp (0) n,m, theA-term in Eq. (E4c) produces an up-shifted pole ofξ-type, theB-term a down-shiftedξ-type, but only ifm >1. Similarly, fora 1-type polesp (1) n,m, theA-term produces an up-shifteda 1-type and theB-term a downshifteda 1-type without exception....

    work page

  3. [3]

    (E6) and follow the same pattern as Eqs

    Contributions toζ n+1,m The result contributions toζ n+1,m from Eq. (E6) and follow the same pattern as Eqs. (F2), (F3) and (F4) with [φ], |φ1|,|φ 2|replaced by [γ],|γ 1|,|γ 2|andπreplaced byζ. From Eq. (E6c), − A(p(0) n,m) 2p(0) n,m [γ](0) n,m 7→ζ (0) n+1,m+1 form= 1, . . . , M (0) n (F5a) − B(p(0) n,m) 2p(0) n,m [γ](0) n,m 7→ζ (0) n+1,m−1 form= 2, . . ....

    work page

  4. [4]

    (E7) and follow a different pattern compared toπandζ

    Contributions toρ n+1,m Contributions toρ n+1,m derive from Eq. (E7) and follow a different pattern compared toπandζ. From Eq. (E7c), A(p(0) n,m) 2 (p(0) n,m +ξ −1)[η](0) n,m 7→ρ (0) n+1,m+1 form= 1, . . . , M (0) n (F8a) B(p(0) n,m) 2 (p(0) n,m −ξ −1)[η](0) n,m 7→ρ (0) n+1,m−1 form= 2, . . . , M (0) n (F8b) A(p(1) n,m) 2 (p(1) n,m +ξ −1)[η](1) n,m 7→ρ (1...

    work page

  5. [5]

    Garcia-Millan, L

    R. Garcia-Millan, L. Cocconi, Z. Zhang, M. Bothe, L. Chen, Z. Zhen, and G. Pruessner, Effective attraction by repulsion (2026)

    work page 2026

  6. [6]

    Fily and M

    Y. Fily and M. C. Marchetti, Athermal phase separation of self-propelled particles with no alignment, Phys. Rev. Lett.108, 235702 (2012)

    work page 2012

  7. [7]

    M. E. Cates and J. Tailleur, Motility-induced phase sepa- ration, Annu. Rev. Condens. Matter Phys.6, 219 (2015)

    work page 2015

  8. [8]

    Digregorio, D

    P. Digregorio, D. Levis, A. Suma, L. F. Cugliandolo, G. Gonnella, and I. Pagonabarraga, Full phase diagram of active brownian disks: from melting to motility-induced phase separation, Phys. Rev. Lett.121, 098003 (2018)

    work page 2018

  9. [9]

    M. E. Cates and C. Nardini, Active phase separation: new phenomenology from non-equilibrium physics, Rep. Progr. Phys.88, 056601 (2025)

    work page 2025

  10. [10]

    Slowman, M

    A. Slowman, M. Evans, and R. Blythe, Jamming and attraction of interacting run-and-tumble random walkers, Phys. Rev. Lett.116, 218101 (2016)

    work page 2016

  11. [11]

    Guillin, L

    A. Guillin, L. Hahn, and M. Michel, Long-time analy- sis of a pair of on-lattice and continuous run-and-tumble particles with jamming interactions, J. Stat. Phys.192, 123 (2025)

    work page 2025

  12. [12]

    A. Das, A. Dhar, and A. Kundu, Gap statistics of two interacting run and tumble particles in one dimension, J. Phys. A Math. Theor.53, 345003 (2020)

    work page 2020

  13. [13]

    L. Hahn, A. Guillin, and M. Michel, Jamming pair of gen- eral run-and-tumble particles: exact results, symmetries and steady-state universality classes, J. Phys. A: Math. Theor.58, 155001 (2025)

    work page 2025

  14. [14]

    M. J. Metson, M. R. Evans, and R. A. Blythe, From a microscopic solution to a continuum description of active particles with a recoil interaction in one dimension, Phys. Rev. E107, 044134 (2023)

    work page 2023

  15. [15]

    L. Hahn, A. Guillin, and M. Michel, Activity-driven clustering of jamming run-and-tumble particles: Ex- act three-body steady state by dynamical symmetry, arXiv:2509.08945 (2025)

    work page arXiv 2025

  16. [16]

    Doi, Second quantization representation for classical 25 many-particle system, J

    M. Doi, Second quantization representation for classical 25 many-particle system, J. Phys. A: Math. Gen.9, 1465 (1976)

    work page 1976

  17. [17]

    Doi, Stochastic theory of diffusion-controlled reaction, J

    M. Doi, Stochastic theory of diffusion-controlled reaction, J. Phys. A: Math. Gen.9, 1479 (1976)

    work page 1976

  18. [18]

    Peliti, Path integral approach to birth-death processes on a lattice, J

    L. Peliti, Path integral approach to birth-death processes on a lattice, J. Phys. (Paris)46, 1469 (1985)

    work page 1985

  19. [19]

    Cardy, Reaction-diffusion processes (2006), accessed 9 Apr 2026, preprint of [?]

    J. Cardy, Reaction-diffusion processes (2006), accessed 9 Apr 2026, preprint of [?]

    work page 2006

  20. [20]

    Pruessner and R

    G. Pruessner and R. Garcia-Millan, Field theories of ac- tive particle systems and their entropy production, Rep. Progr. Phys.88, 097601 (2025)

    work page 2025

  21. [21]

    Hansen and I

    J.-P. Hansen and I. R. McDonald,Theory of simple liq- uids(Academic Press, London, UK, 2006)

    work page 2006

  22. [22]

    Zhang and G

    Z. Zhang and G. Pruessner, Field theory of free run and tumble particles in d dimensions, J. Phys. A55, 045204 (2022)

    work page 2022

  23. [24]

    Bothe, L

    M. Bothe, L. Cocconi, Z. Zhen, and G. Pruessner, Parti- cle entity in the doi-peliti and response field formalisms, J. Phys. A Math. Theor.56, 175002 (2023)

    work page 2023

  24. [25]

    Garcia-Millan and G

    R. Garcia-Millan and G. Pruessner, Run-and-tumble mo- tion in a harmonic potential: field theory and entropy production, J. Stat. Mech.: Theory Exp.2021(6), 063203

    work page 2021

  25. [26]

    J. L. Schiff,The Laplace transform: theory and applica- tions(Springer Science & Business Media, 1999)

    work page 1999

  26. [27]

    Optimal ratchet po- tentials for run-and-tumble particles.arXiv:2204.04070, 2022

    Z. Zhen and G. Pruessner, Optimal ratchet potentials for run-and-tumble particles, arXiv:2204.04070 (2022)

    work page arXiv 2022

  27. [28]

    Zhang and R

    Z. Zhang and R. Garcia-Millan, Entropy production of nonreciprocal interactions, Phys. Rev. Res.5, L022033 (2023)

    work page 2023

  28. [29]

    Le Bellac,Quantum and Statistical Field Theory [Phe- nomenes critiques aux champs de jauge, English](Oxford University Press, New York, NY, USA, 1991) translated by G

    M. Le Bellac,Quantum and Statistical Field Theory [Phe- nomenes critiques aux champs de jauge, English](Oxford University Press, New York, NY, USA, 1991) translated by G. Barton

    work page 1991

  29. [30]

    Zhang, L

    Z. Zhang, L. Feh´ ert´ oi-Nagy, M. Polackova, and G. Pruessner, Field theory of active brownian particles in potentials, New J. Phys.26, 013040 (2024)

    work page 2024

  30. [31]

    Matsubara, A new approach to quantum-statistical mechanics, Prog

    T. Matsubara, A new approach to quantum-statistical mechanics, Prog. Theor. Phys.14, 351 (1955)

    work page 1955

  31. [32]

    Espinosa, On the evaluation of Matsubara sums, Math

    O. Espinosa, On the evaluation of Matsubara sums, Math. Comput.79, 1709 (2010)

    work page 2010

  32. [33]

    Spiegel, S

    M. Spiegel, S. Lipschutz, J. Schiller, and D. Spellman, Complex variables with an introduction to conformal mapping and its applications(McGraw-Hill, New York, NY, USA, 2009)

    work page 2009

  33. [34]

    J. E. Marsden and M. J. Hoffman,Basic complex analy- sis, Vol. 3 (WH Freeman, New York, NY, USA, 1973)

    work page 1973

  34. [35]

    Sekimoto, Stochastic energetics (Springer-Verlag, Berlin, Germany, 2012) pp

    K. Sekimoto, Stochastic energetics (Springer-Verlag, Berlin, Germany, 2012) pp. I–XVIII, 1–322

    work page 2012

  35. [36]

    Cocconi, R

    L. Cocconi, R. Garcia-Millan, Z. Zhen, B. Buturca, and G. Pruessner, Entropy production in exactly solvable sys- tems, Entropy22, 1252 (2020)

    work page 2020