Pith. sign in

REVIEW 2 major objections 4 minor 43 references

Looking Back: Field theory of transiently chiral active particles

T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper derives a Doi-Peliti field theory for transiently chiral active particles and shows that the mean squared displacement, reorientation-angle distribution, and orientation correlation functions follow in closed form from a…

desk verdict Solid self-contained field-theoretic derivation for a new solvable active-matter model; the exactness claim holds, but some printed formulas have removable singularities at the parameter values used in two of the figures. read the letter →

arxiv 2507.01503 v1 pith:M3YS6I6S submitted 2025-07-02 cond-mat.stat-mech

Callum Britton , Gunnar Pruessner , Thibault Bertrand This is my paper · ORCID
classification cond-mat.stat-mech PACS 05.40.-a
keywords transientlychiralactiveparticlesDoi-Pelitifieldtheoryrun-and-tumbleBrownianmeansquareddisplacementorientationautocorrelatorreorientationmemorystochastic
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper derives a Doi-Peliti field theory for transiently chiral active particles (TCAPs): active Brownian particles whose reorientation angle at each tumble is itself a diffusing degree of freedom, so successive tumble angles remain correlated over a memory time $1/D_\psi$. The formal claim is that the field-theory action exactly represents the Fokker-Planck equation of the underlying Langevin process, so every observable computed from it is an exact property of the model rather than a perturbative approximation. From that action the paper obtains closed-form mean squared displacements for uniformly random and fixed initial reorientation angles, the reorientation-angle probability density, the orientation autocorrelator, and two position-orientation cross-correlators. The results show how slow reorientation-angle diffusion produces transiently chiral trajectories, with particles tumbling in the same handedness for roughly $\gamma/D_\psi$ tumbles before the memory fades, and how that chiral signature disappears once the initial reorientation angle is averaged out. Exact solvability of a correlated-reorientation active-particle model supplies simulation benchmarks and analytic predictions for single-particle tracking of bacteria with tumble memory.

What carries the argument

The central object is the Doi-Peliti field theory built from the annihilation field $\chi(x,\phi,\psi,t)$ and its Doi-shifted response field $\tilde\chi=1+\tilde\chi$, whose action (7) is an exact rewrite of the Fokker-Planck equation (6). The working machinery is the Fourier expansion in modes $n$ (director angle) and $m$ (reorientation angle), which diagonalises the harmonic part and converts the two difficult interactions into simple momentum-space vertices: self-propulsion shifts $n$ by $\pm1$ with weight proportional to $k$, and tumbling shifts $m$ by $n$ with weight $\gamma$. Because the observables are evaluated at $k=0$ with at most two derivatives with respect to $k$, only diagrams with zero, one, or two self-propulsion vertices survive, and arbitrary numbers of tumble vertices resum as geometric series. This selection rule is what turns the perturbative expansion into exact closed-form expressions.

What would settle it

Record a long tumble sequence from a single swimmer in a regime where the reorientation angle decorrelates slowly compared with the tumble interval, and compare the empirical probability that the next tumble is same-handed with Eq. (47b) and the waiting-time distribution with an exponential; a systematic mismatch in either comparison would show that the model's core premise does not describe the measured trajectories.

Watch

Extended reading notes

Core claim

The paper's central claim is that the stochastic process defined by the Langevin equations for a TCAP—constant-speed self-propulsion along a diffusing director $\phi$, translational diffusion $D_x$, Poissonian tumbles at rate $\gamma$, and a reorientation angle $\psi$ that itself diffuses with diffusivity $D_\psi$—is exactly represented by the Doi-Peliti action (7), which is a direct transcription of the Fokker-Planck equation (6). In Fourier modes $n$ (for $\phi$) and $m$ (for $\psi$), the harmonic part of the action is diagonal, with bare propagator $G(k,n,m,\omega)=1/(-i\omega+D_x k^2+D_\phi n^2+D_\psi m^2+\gamma+r)$, and the perturbations are just three vertices: two self-propulsion vertices that shift $n$ by $\pm1$ and one tumble vertex that shifts $m$ by $n$. The paper shows that for each observable considered the diagrammatic series terminates or resums exactly, because the $k$-derivatives used to form position moments kill all diagrams with more than two self-propulsion vertices while chains of tumble vertices collapse into geometric series. All final expressions—the mean squared displacement (33) and (39), the reorientation-angle density (46), the right-handed-tumble probability (47b), the orientation autocorrelator (54), and the position-orientation cross-correlators (61) and (65)—are therefore exact consequences of the model.

Load-bearing premise

The load-bearing premise is the model choice that a real particle's reorientation angle diffuses freely between tumbles and that tumbles occur as independent events at a fixed average rate; if biological or experimental reorientation statistics violate either part, the exact formulas describe the model rather than the organism.

Editorial extensions

If this is right

  • After averaging over the initial reorientation angle, the mean squared displacement is exactly $\langle |\mathbf{x}|^2(t)\rangle = 4D_x t + \frac{2v^2}{(D_\phi+\gamma)^2}\left[e^{-(D_\phi+\gamma)t}-1+(D_\phi+\gamma)t\right]$, so this observable cannot distinguish rotational diffusion from tumbling.
  • With a fixed initial reorientation angle, the MSD gains corrections proportional to $\cos(m\psi_0)\gamma^m$; because $\gamma^m$ carries the number of tumbles, the corrections quantify how many correlated tumbles occur before the reorientation memory fades, and in the slow-decorrelation regime $\tau_p \ll t \ll \tau_\psi$ the trajectory can be subdiffusive for small $\psi_0$.
  • The reorientation-angle distribution is diffusion on a ring, $P(\psi,t|\psi_0)=\frac{1}{2\pi}\vartheta_3((\psi-\psi_0)/2,e^{-D_\psi t})$, so the probability of a next right-handed tumble is $P(R|\psi_0)=\frac12+\frac{2}{\pi}\sum_k \frac{\sin[(2k+1)\psi_0]}{(2k+1)(1+(2k+1)^2D_\psi/\gamma)}$, which for $D_\psi/\gamma\ll1$ predicts about $\gamma/D_\psi$ consecutive same-handed tumbles.
  • The position-perpendicular orientation cross-correlator is antisymmetric in $\psi_0$ and decays at long times, giving a direct, time-resolved signature of the direction and duration of transient chirality.
  • Because the derivation is exact, all formulas serve as benchmarks for numerical simulations of the Langevin model without adjustable parameters.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An extension the paper leaves implicit is that fitting the closed-form MSD and orientation autocorrelator to single-particle tracking data could estimate $D_\psi/\gamma$ directly from the length of the run of consecutive same-handed tumbles, without resolving $\psi$ at every tumble.
  • The exactness mechanism—moment observables evaluated at zero wavevector truncate the self-propulsion vertex expansion—suggests that other low-order spatial moments, such as velocity autocorrelations or confinement-modified first-passage statistics, may also be computable in closed form for this model.
  • If the predicted intermediate subdiffusive regime for small fixed $\psi_0$ is found in simulations or experiments, it would distinguish transient chirality from ordinary run-and-tumble or active Brownian motion, which do not show that regime.
  • Adding a harmonic potential for $\psi$, as the paper suggests, is the simplest stress test of the formalism: the tumble vertex would be unchanged and the same diagram selection should still yield exact MSD corrections.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript introduces a model of transiently chiral active particles (TCAPs), in which the reorientation angle ψ undergoes free diffusion and tumble events occur as a Poisson process with rate γ. The authors derive a Doi-Peliti field theory from the Langevin/Fokker-Planck equations and use diagrammatic perturbation theory to compute closed-form expressions for the mean squared displacement (for uniform and fixed initial ψ0), the reorientation-angle probability density, the probability of a right-handed tumble, the orientation autocorrelator, and the position-orientation and position-perpendicular-orientation cross-correlators. The analytic results are compared with numerical simulations in several figures. The paper claims that all observables are exact, arising from an exact representation of the stochastic process, and that no parameters are fitted.

Significance. If the exactness claim holds, the paper provides a useful exactly solvable field-theoretic model for active particles with correlated tumbling, with all predictions expressed in the model parameters (v, Dx, Dϕ, Dψ, γ). The derivation is self-contained and carefully documented, including diagrammatic selection rules, t→0 sanity checks, and comparisons with simulations. These are genuine strengths. However, the central exactness claim is undercut by an unaddressed regularity issue: several of the final closed-form expressions are singular at the parameter values used in the main figures, and the manuscript does not state the limiting convention by which these expressions are evaluated. Because the issue is local and fixable, it does not invalidate the field-theoretic approach, but it must be resolved before the exactness claim can be accepted as stated.

major comments (2)
  1. [Eqs. (39), (41), (60), (61), (65)] The denominators Dϕ + Dψ(ℓ² − m²) + γ vanish at m = 1, ℓ = 0 when Dψ = Dϕ + γ. The parameters used in Figs. 6–7, Dϕ = 0, Dψ = 1, γ = 1, lie exactly on this resonance. The manuscript never states that the expressions are defined by continuity, so as printed the plotted formulas are undefined at the shown parameter values. I checked the m = 1 case: the pole from the term involving 1/(Dϕ − Dψ + γ) is cancelled by the ℓ = 0 term in the sum, giving a finite limit, so the singularity is removable. Nevertheless, the text should state this limiting convention explicitly and, ideally, display the limiting form or provide code/derivation showing how the limit was evaluated. Without this, the label “exact” is not self-contained for the central observables of Sections 6 and 7.
  2. [Figs. 6–7 and Section 7] The numerical simulations used to verify Eqs. (60)/(61) and (65) are compared at exactly the resonant parameter value Dψ = Dϕ + γ, yet the paper gives no simulation details: no integration scheme, time step, ensemble size, or error bars are reported. Because the analytic curves are singular at those parameters without a stated limiting convention, it is impossible for a reader to reproduce the comparison or to judge whether the simulation data were evaluated with the same limit. Please specify the numerical evaluation of the sums (including truncation) and the handling of the removable singularities, and provide sufficient simulation details for reproducibility.
minor comments (4)
  1. [Section 4.2.2] The text says “The data shown in Equation (3)” but should refer to Figure 3.
  2. [Section 8] “reorientaion” is a typo for “reorientation”.
  3. [Eq. (24) and surrounding notation] The primed variables k′ and ω′ appear without being defined in Eq. (24) and elsewhere; the delta functions later make their meaning clear, but an early sentence defining the convention would improve readability.
  4. [Throughout] The manuscript would benefit from a brief statement of the parameter regime in which the closed forms are valid, e.g., Dψ > 0 and Dϕ, γ ≥ 0, since several expressions contain denominators such as Dψ m² and Dψ(i² − j²) that are singular at Dψ = 0 or for coincident lattice indices.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Doi-Peliti action is built explicitly from the Fokker-Planck equation, and every observable is evaluated from the stated propagator and vertices without fitted parameters.

full rationale

The derivation chain is self-contained and not circular. Equation (6) is a direct transcription of the Langevin dynamics (1)-(5), and Equation (7) writes the corresponding Doi-Peliti action out explicitly, so the accompanying citation [19], though written by an author of the present paper, is not load-bearing; the mapping from Fokker-Planck to action is standard and is exhibited in the manuscript itself. All reported observables, Equations (33), (39), (41), (46), (47), (52)-(54), (60), (61) and (65), are obtained from the bare propagator (16) and the explicit vertices (17)-(18) by exact summations and partial-fraction identities; no parameter is fitted to a subset of data and then relabelled a prediction, and the numerical simulations are independent solutions of the stochastic process used only for comparison. The recoveries of the ABP-with-tumbling MSD and of diffusion on a ring are cross-checks of the formalism rather than renamings of the target results. The one substantive defect found is a regularity gap, not circularity: the closed forms contain denominators D_phi + D_psi(ell^2 - m^2) + gamma that vanish at D_psi = D_phi + gamma for (m, ell) = (1, 0), the parameters of Figures 6 and 7, with no stated limiting convention; this affects the exactness claim as printed, but it does not make any result equivalent to its input by construction.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The model introduces a new stochastic process (TCAP), but no new physical entities beyond the model's degrees of freedom (position, orientation, reorientation angle). The derivation rests on standard field theory techniques and the explicit model assumptions. There are no free parameters fitted to data; all results are parameterized by the model's physical constants.

assumptions (4)
  • domain assumption The tumble events occur as a Poisson process with rate γ, independent of the orientation and reorientation angle.
    This is part of the model definition, stated in Section 2. The field theory derivation follows from this assumption. It is not externally validated but is a standard modeling choice.
  • domain assumption The reorientation angle ψ evolves by free diffusion between tumbles, with diffusion constant Dψ.
    Stated in Section 2.1, Equation (4). This is the key new modeling ingredient, and it is used throughout the derivations.
  • standard math The Doi-Peliti action exactly represents the Fokker-Planck equation.
    This is a standard result in field theory, referenced from [13,14,37-40]. The paper relies on this to translate Equation (6) to Equation (7).
  • standard math The Fourier transform conventions and the use of the normal-ordering prescription are valid.
    The paper uses a specific Fourier representation and normal-ordering (Section 3). These are standard techniques in Doi-Peliti field theory.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Looking Back: Field theory of transiently chiral active particles." pith.science (2026). https://pith.science/paper/M3YS6I6S

@misc{pith2026250701503,
  author       = {Pith},
  title        = {Pith review of: Looking Back: Field theory of transiently chiral active particles},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/M3YS6I6S}},
  note         = {Machine review of arXiv:2507.01503}
}
read the original abstract

We derive a Doi-Peliti Field Theory for transiently chiral active particles in two dimensions, that is, active Brownian particles that undergo tumbles via a diffusing reorientation angle. Using this framework, we compute the mean squared displacement for both uniformly distributed and fixed initial reorientations. We also calculate an array of orientation-based observables, to quantify the transiently chiral behaviour observed.

Figures

Figures reproduced from arXiv: 2507.01503 by the authors.

Figure 1
Figure 1. Example trajectories for free TCAP particles — (a) schematic of successive tumbling events with reorientation angles {ψ(t1), ψ(t2), . . .}. Representative trajectories are shown for two free TCAP with Dx = 0, Dϕ = 1, γ = 1000, v = 1, ψ0 = π/12 and (b) Dψ = 0.0001, (c) Dψ = 1 and (d) Dψ = 10000. 2.2. Transient chirality In principle, there are six regimes to consider, depending on the relative ordering of the three t… view at source ↗
Figure 2
Figure 2. Mean squared displacement with uniformly chosen initial ψ0 — The MSD, Equation (33), with uniform initial ψ0 for τp = 1, Dx = 0.01 and a range of self￾propulsion velocities v, as indicated. Symbols indicate numerical simulations, solid lines Equation (33) [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Mean squared displacement with fixed initial ψ0 — The MSD for a TCAP with fixed initial ψ0, Equations (39) and (41), is plotted in the limit of slow decorrelation of ψ with (a) low tumble rate and high diffusivity of the director γ = 0.1, Dϕ = 10, Dψ = 0.1, Dx = 0.01 and (b) high tumble rate and low diffusivity of the director γ = 10, Dϕ = 0.1, Dψ = 0.1, Dx = 0.01 for various values of ψ0. Symbols indicate numerical… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Emergent transient chirality — The probability P(R|ψ0), Equation (47b), to observe a right-handed tumble following a tumble with reorientation ψ0 as (a) a function of ψ0 for a range of Dψ/γ and (b) a function of Dψ/γ for a range of ψ0. the desired result for diffusion …
Figure 5
Figure 5. Figure 5: Orientation autocorrelator — ⟨e(ϕ(t))·e(ϕ0)⟩ plotted as a function of time, t, for a variety of initial reorientation angles ψ0 with Dϕ = 0, Dψ = 1, γ = 1. Symbols indicate numerical simulations, solid lines Equation (54). Similarly, for ⟨sin (ϕ)⟩ we find ⟨sin (ϕ)⟩ = e…
Figure 6
Figure 6. Figure 6: Position-orientation cross-correlator — ⟨x(t)·e(ϕ(t))⟩ plotted as a function of time, t, for a range of initial reorientation angles ψ0 with Dϕ = 0, Dψ = 1, γ = 1. Symbols indicate numerical simulations, solid lines Equation (61). Integrating and taking r → 0 we find ⟨…
Figure 7
Figure 7. Figure 7: Position-perpendicular orientation cross-correlator — ⟨x(t) · e ⊥(ϕ(t))⟩ plotted as a function of time, t, for a variety of initial reorientation angles ψ0 with Dϕ = 0, Dψ = 1, γ = 1. Symbols indicate numerical simulations, solid lines Equation (65). an observable that…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 42 canonical work pages

  1. [1]

    Cates M E 2012Rep. Progr. Phys. 75(4) 042601

    work page

  2. [2]

    Bechinger C, Di Leonardo R, Löwen H, Reichhardt C, Volpe G and Volpe G 2016Rev. Mod. Phys. 88(4) 045006

    work page

  3. [3]

    Schnyder S K, Molina J J, Tanaka Y and Yamamoto R 2017Sci Rep 7(1) 5163

    work page

  4. [4]

    Bothe M, Lardet E, Poliakov A, Pruessner G, Bertrand T and Bordeu I 2025 Eph-ephrin-mediated differential persistence as a mechanism for cell sorting arXiv:2503.13714

    work page Pith review arXiv 2025

  5. [5]

    Kurzthaler C, Zhao Y, Zhou N, Schwarz-Linek J, Devailly C, Arlt J, Huang J D, Poon W C, Franosch T, Tailleur J and Martinez V A 2024Phys. Rev. Lett. 132(3) 038302

    work page

  6. [6]

    Howse J R, Jones R A L, Ryan A J, Gough T, Vafabakhsh R and Golestanian R 2007Phys. Rev. Lett. 99(4) 048102

    work page

  7. [7]

    Tailleur J and Cates M E 2008Phys. Rev. Lett. 100(21) 218103

    work page

  8. [8]

    Cates M E and Tailleur J 2013Europhys. Lett. 101(2) 20010

    work page

Show all 43 references
  1. [9]

    Berg H C and Brown D A 1972Nature 239 500–504

  2. [10]

    Berg H C 1993Random Walks in Biology (Princeton University Press)

  3. [11]

    coli in Motion (Springer)

    Berg H C 2004E. coli in Motion (Springer)

  4. [12]

    Patteson A E, Gopinath A, Goulian M and Arratia P E 2015Sci Rep 5(1) 15761

  5. [13]

    Doi M 1976J. Phys. A: Math. Gen. 9(9) 1465–1477

  6. [14]

    Peliti L 1985J. Phys. (Paris) 46(9) 1469–1483

  7. [15]

    Martin P C, Siggia E D and Rose H A 1973Phys. Rev. A 8(1) 423–437

  8. [16]

    Janssen H K 1976Z. Phys. B 23 377–380

  9. [17]

    de Dominicis C 1976J. Phys. (Paris) Colloque C1 37 C1247–C1253

  10. [18]

    Zhang Z, Fehértói-Nagy L, Polackova M and Pruessner G 2024New J. Phys. 26(1) 013040 8 CONCLUSION 26

  11. [19]

    Pruessner G and Garcia-Millan R 2025 Field theories of active particle systems and their entropy production arXiv:2211.11906

    2025 arXiv

  12. [20]

    Santra I, Basu U and Sabhapandit S 2020Phys. Rev. E 101(6) 062120

  13. [21]

    Banerjee T, Majumdar S N, Rosso A and Schehr G 2020Phys. Rev. E 101(5) 052101

  14. [22]

    Redig F and Wiechen H V 2024Markov Process Relat 30(2024 №2 (30)) 297–331

    2024

  15. [23]

    Friedrich B M and Jülicher F 2007PNAS 104(33) 13256–13261

  16. [24]

    Bertrand T, Zhao Y, Bénichou O, Tailleur J and Voituriez R 2018Phys. Rev. Lett. 120(19) 198103

  17. [25]

    Rashid S, Long Z, Singh S, Kohram M, Vashistha H, Navlakha S, Salman H, Oltvai Z N and Bar-Joseph Z 2019PNAS 116(24) 11770–11775

  18. [26]

    Zhao Y, Kurzthaler C, Zhou N, Schwarz-Linek J, Devailly C, Arlt J, Huang J D, Poon W C K, Franosch T, Martinez V A and Tailleur J 2024Phys. Rev. E 109(1) 014612

  19. [27]

    Berg H and Purcell E 1977Biophys. J. 20(2) 193–219

  20. [28]

    Keller E F and Segel L A 1971J. Theor. Biol. 30(2) 225–234

  21. [29]

    Endres R G and Wingreen N S 2008PNAS 105(41) 15749–15754

  22. [30]

    Knight J, Garcia-Galinda P, Pausch J and Pruessner G 2024J. R. Soc. Interface. 21(216) 20240100

  23. [31]

    Miranda S W, Asfahl K L, Dandekar A A and Greenberg E P 2022Adv Exp Med Biol 1386 95–115

  24. [32]

    d’Alessandro J, Barbier-Chebbah A, Cellerin V, Benichou O, Mège R, Voituriez R and Ladoux B 2021 Nature Communications 12 4118

    2021

  25. [33]

    Marconi U M B and Maggi C 2015Soft Matter 11(45) 8768–8781

  26. [34]

    Potdar A A, Jeon J, Weaver A M, Quaranta V and Cummings P T 2010PLoS ONE 5(3) e9636

  27. [35]

    Purcell E M 1977Am. J. Phys 45(1) 3–11

  28. [36]

    thesis Imperial College London, London, United Kingdom

    Pausch J 2019 Topics in Statistical Mechanics Ph.D. thesis Imperial College London, London, United Kingdom

    2019

  29. [37]

    Täuber U C, Howard M and Vollmayr-Lee B P 2005J. Phys. A: Math. Gen. 38(17) R79–R131

  30. [38]

    Cardy J 2008 Reaction-diffusion processesNon-equilibrium Statistical Mechanics and Turbulence ed Nazarenko S and Zaboronski O V (Cambridge, UK: Cambridge University Press) pp 108– 161 London Mathematical Society Lecture Note Series: 355, preprint available from (http: //www-th...

    2008

  31. [39]

    Pruessner G 2011 Lecture notes on non-equilibrium statistical mechanics, ch 6 accessed 26 June 2025

    2011

  32. [40]

    Täuber U C 2014 Critical Dynamics: A Field Theory Approach to Equilibrium and Non- Equilibrium Scaling Behavior 1st ed (Cambridge, UK: Cambridge University Press)

    2014

  33. [41]

    Zhang Z and Pruessner G 2022J. Phys. A 55(4) 045204

  34. [42]

    Zhang Z, Fehértói-Nagy L, Polackova M and Pruessner G 2024New J. Phys. 26(1) 013040

  35. [43]

    Gradshteyn I S and Ryzhik I M 2007 (1943) 8.18–8.19 theta functionsTable of Integrals, Series, and Products (New York, USA: Academic Press) pp 877–883 7th ed

    1943

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.