Pith. sign in

REVIEW 3 major objections 4 minor 43 references

Active Brownian particles in a 3D tube whose walls oscillate in a traveling wave move fastest at an intermediate driving frequency, a resonance-like effect that also boosts their effective diffusion.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-04 17:00 UTC pith:EG7XQS2A

load-bearing objection Solid extension of FJ to pulsating 3D tubes with a resonance peak corroborated by simulations; main weakness is the unproven moving-boundary reduction and the choice of parameters that violate the paper's own smoothness condition. the 3 major comments →

arxiv 2509.11147 v1 pith:EG7XQS2A submitted 2025-09-14 cond-mat.soft cond-mat.stat-mechphysics.bio-ph

Entropic active particle transport in pulsating 3D geometries

Rahul Sinha , Ankit Gupta , P. S. Burada This is my paper
classification cond-mat.soft cond-mat.stat-mechphysics.bio-ph MSC 82C3182C70
keywords active Brownian particlesFick-Jacobs approximationentropic barriersstochastic resonanceoscillating tubeeffective diffusionmicrofluidicsphase velocity
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper argues that active Brownian particles in a 3D tube with a pulsating, traveling-wave boundary exhibit a peak in both average velocity and effective diffusion at a specific oscillation frequency. This resonance-like enhancement grows with the particles' self-propulsion speed and the wall oscillation amplitude, and the behavior is captured by a generalized Fick-Jacobs approximation. A sympathetic reader would see this as a tunable transport mechanism: by adjusting the wall frequency, one can control how quickly and how broadly active particles spread through a confined channel.

Core claim

On the paper's own terms, the central claim is that for an active Brownian particle in a spatially periodic 3D tube with radius W(z,t)=a sin(2πz/L−ωt)+b, the ensemble-averaged axial velocity and effective diffusion coefficient attain maximum values at an optimal frequency of wall oscillation. This phenomenon is qualitatively and semi-quantitatively reproduced by a Fick-Jacobs equation that reduces the 3D dynamics to a 1D transport problem with a time-dependent entropic potential. The resonance is enhanced by increasing self-propulsion speed v0 and amplitude a. In the limit where the bottleneck radius (b−a) goes to zero, the mean velocity saturates at the phase velocity of the wall (ωL/2π) wh

What carries the argument

The carrying tool is the generalized Fick-Jacobs equation, which reduces the full 3D Fokker-Planck dynamics to an effective 1D equation along the tube axis z. The central objects are the time- and position-dependent entropy potential A(z,t)=−k_B T ln(π W^2/L^2) and an effective diffusion coefficient D(z,t) = (D0 + v0^2/[d(d−1)D_r])[1+(h')^2]^{-α}, with α=1/2 in 3D. This maps the moving corrugated tube to a 1D advection-diffusion problem with a traveling entropic barrier, allowing semi-analytical computation of <v> and D_eff.

Load-bearing premise

The Fick-Jacobs reduction assumes that transverse equilibration remains fast when the boundary moves, so the static form of the entropy potential and diffusivity applies with no additional drift from the moving wall; the paper's own parameters give a wall slope of unity, beyond the stated |W'|<<1 condition under which that assumption is justified.

What would settle it

Compute the full 3D Fokker-Planck dynamics with the moving-wall boundary condition, including the ∂_t ln h(z,t) term that the static Fick-Jacobs form omits, and compare the predicted <v>(ω) at the paper's parameters. If the moving-wall drift term is non-negligible at the resonance frequency, the peak position or height will shift; alternatively, a microfluidic experiment tracking particle trajectories in an oscillating tube at the same dimensionless parameters would reveal the same mismatch.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Microfluidic and nanofluidic devices could use the optimal oscillation frequency to maximize throughput of active particles, with the frequency setting the characteristic speed.
  • Because the peak frequency is independent of self-propulsion speed, a fixed device geometry could control transport of particles with different activity levels at the same optimal frequency.
  • Near the closed-bottleneck limit, particle velocity locks to the wall phase velocity, offering a way to synchronize particle motion with an external mechanical drive.
  • The generalized Fick-Jacobs approach is extendable to other time-dependent confinements, such as peristaltic channels or ciliated surfaces, making it a basis for predicting transport in biological and synthetic systems.

Where Pith is reading between the lines

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

  • If valid, the resonance implies that even passive Brownian particles would show a (weaker) peak at a similar frequency; measuring the v0-dependence of the peak height would directly test whether activity is the amplifier the paper claims.
  • The saturation of <v> to the wall phase velocity suggests an experimental probe: tracking inert particles in a pulsating tube could measure wall oscillation speed without optical access to the wall itself.
  • A standing-wave boundary (sin(2πz/L)cos(ωt)) should produce zero net velocity by symmetry, providing a control experiment to confirm that the traveling wave is the source of the directed transport.
  • The paper's parameter choice (a=1/2π, b=2/2π) yields |W'|=1, violating the stated smoothness condition; testing the resonance at gentler slopes (smaller a/L) would reveal whether the peak position and height are robust or artifacts of the approximation regime.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. This manuscript studies transport of non-interacting active Brownian particles (ABPs) in a 3D tube whose radius is a traveling wave W(z,t)=a sin(2πz/L−ωt)+b. The authors reduce the 3D Fokker–Planck equation to a 1D Fick–Jacobs (FJ) equation (Eq. 4) by assuming fast transverse equilibration and replacing the static width h(z) by h(z,t). They solve this 1D equation to obtain the mean velocity <v> and the effective diffusion coefficient D_eff (Eqs. 10–11), and compare with 3D Brownian dynamics simulations. The central claims are a resonance-like maximum in <v> and D_eff at an optimal frequency ω, enhancement with activity v0 and amplitude a, and saturation of <v> to the phase velocity v_wall=ωL/(2π) with D_eff→0 when the bottleneck closes (a→b).

Significance. If the FJ reduction is valid for moving boundaries, the paper provides a simple semi-analytical framework and a falsifiable prediction: active particles in a pulsating periodic tube can be transported at a tunable velocity with a maximum at a specific driving frequency. The use of a parameter-free semi-analytical solution together with independent 3D Brownian dynamics simulations is a genuine strength; the resonance is predicted by the model and confirmed in simulation, not fitted. The saturation of <v> to the phase velocity and the vanishing of D_eff in the closed-bottleneck limit are conceptually appealing. However, the result rests on an unproved extension of the static FJ approximation to time-dependent boundaries, and the parameters used violate the paper's own smoothness condition, so the quantitative validity of the central prediction needs additional support.

major comments (3)
  1. [Semi-analytical approach, Fick–Jacobs equation, Eq. (4) and (8)] The reduction from the 3D Fokker–Planck equation (3) to Eq. (4) is taken from the static-channel FJ result (Ref. [21]) by replacing h(z) with h(z,t), but the moving-wall no-flux condition is not treated. After integrating the 3D continuity equation over the moving cross-section and imposing no-flux, the boundary velocity cancels in the integrated continuity equation, so the specific ∂_t ln h term suggested by the stress-test note is not automatic. Nevertheless, the constitutive law for the integrated longitudinal flux J(z,t) in Eq. (8) is not derived for time-dependent boundaries, and corrections of order ∂_t h or a transverse-relaxation lag may enter. Since all semi-analytical curves in Figs. 2 and 4 come from Eq. (4), this missing derivation is load-bearing. The authors should either derive the time-dependent FJ equation from Eq. (3) with the moving-wall boundary condition, or provide
  2. [Parameters used in Figs. 2–4 vs. stated FJ validity condition] The paper states that 'the FJ approximation requires smooth boundary variations, i.e., |W′(z,t)|≪1.' For the parameters used throughout, L=1 and a=1/(2π), we have max|W′| = 2πa/L = 1. This violates the stated condition by an order of magnitude. The later claim that 'FJ approximation is valid even in the case of oscillatory boundaries' is therefore not supported by the parameter set in Figs. 2–4. The authors should demonstrate the approximation for at least one parameter set with max|W′| ≪ 1, and explain why the very good agreement persists at max|W′|=1.
  3. [Figs. 2–4 and simulation details] The numerical corroboration is limited to active particles (v0>0), with no error bars, a single geometry, and no passive (v0=0) case. Because Eq. (7) is the passive limit of the moving-boundary FJ equation, a v0=0 comparison would directly test whether the 1D reduction is correct without the complication of orientational dynamics. Please include error bars or a statement of statistical uncertainty, and show the passive case and at least one additional amplitude/frequency set.
minor comments (4)
  1. [Eqs. (14)–(15)] The long-time definitions of ⟨v⟩ and D_eff should state explicitly that the averages are over the 3D Brownian dynamics ensemble and that t is the simulation time; also clarify how the long-time limit is taken in practice.
  2. [Fig. 4 caption] 'The left arrows indicate the values of saturated mean velocity ⟨v⟩ = vw = ωL/2π beyond a = b' is confusing because a=b is the closed-bottleneck limit; 'as a→b' would be clearer.
  3. [Numerical simulations] The text reads '2 × 104 trajectories' and '10 −4'; these should be typeset as 2×10^4 and 10^{-4}.
  4. [General] Minor grammatical issues: 'The dynamics of ABP is governed' should be 'ABPs'; 'Stochastic Resonance' should be 'stochastic resonance'; 'project-wide CRG/2023/006186' should be 'project no.'.

Circularity Check

0 steps flagged

No circularity: the FJ predictions are solved semi-analytically and validated against independent BD simulations; self-citations are to prior published results and no fitted parameter is renamed as a prediction.

full rationale

The paper's transport predictions are not equivalent to its inputs by construction. The generalized Fick-Jacobs equation (Eq. 4) is asserted by adapting the static FJ reduction from Ref. [21] to a time-dependent boundary, but the subsequent calculation is a genuine solution of that 1D PDE: the mean velocity and effective diffusivity in Eqs. (10)-(11) are obtained by quadrature from the solved probability density, not by fitting to the simulation data. The Brownian dynamics simulations in Figs. 2-4 use the same independently specified parameters (D0=Dr=1, tube geometry a=1/2pi, b=2/2pi, L=1, v0 as input) and are not used to determine any coefficient in the FJ equation; the agreement is an a posteriori validation. The active enhancement enters through published expressions for D(z,t) (Eq. 5, Refs. [20,22]) and the advective v0 cos(theta) term, and the resonance peak is an emergent feature of the solution rather than a fitted input. The limit a->b, where <v> approaches the wall phase velocity, is explicitly said to occur where the FJ approximation fails; it follows kinematically from the traveling no-flux boundary condition, not from the 1D reduction, so it is not a self-definitional prediction. Self-citations such as Ref. [21] (static FJ) are to published, externally falsifiable prior work, and no uniqueness theorem is invoked to forbid alternatives. One caveat is flagged: the paper states the FJ condition |W'|<<1, but the chosen parameters give max|W'|=1, and the time-dependent FJ form is not derived from the exact moving-wall boundary condition; these are correctness/validity concerns rather than circularity, since they do not make the output equal to the input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

No parameters are fitted to data; all model parameters (a, b, omega, v0, D0, Dr) are physical inputs. The central claim rests on the adiabatic Fick-Jacobs reduction, the active diffusion coefficient expression from Refs. [20,22], and the point-particle no-hydrodynamics idealization.

axioms (5)
  • domain assumption Fast transverse equilibration: the 3D Fokker-Planck equation can be integrated over x and y because the transverse diffusion time is much shorter than the axial diffusion time.
    Invoked after Eq. (3) to obtain Eq. (4); standard FJ statement from Ref. [21]. It breaks at high omega and near a = b.
  • ad hoc to paper The static Fick-Jacobs form with D(z,t) from Eq. (5) and A(z,t) from Eq. (6) remains valid when the boundary moves, i.e., no additional wall-motion terms appear.
    This is the paper's key modeling assumption; it is not derived for time-dependent boundaries and the paper notes FJ fails at high omega.
  • domain assumption Particles are point-like and walls are purely reflecting; hydrodynamic interactions are neglected.
    Stated in the numerical simulations section; standard for entropic transport models.
  • domain assumption Tube parameters satisfy b > a so bottlenecks remain open and particles can pass between cells.
    Stated after Eq. (2); needed for transport along the tube.
  • domain assumption The active Brownian particle dynamics of Eq. (1) with constant self-propulsion speed v0 and orientation diffusion Dr is the correct minimal model.
    Standard ABP model (Refs. [16-18]); assumed throughout.

pith-pipeline@v1.3.0-alltime-deepseek · 7844 in / 13907 out tokens · 154526 ms · 2026-08-04T17:00:50.633294+00:00 · methodology

0 comments
read the original abstract

We study the transport of active Brownian particles (ABPs) in three-dimensional (3D) oscillatory geometries, which are spatially periodic. We establish a generalized Fick-Jacobs approach, which reduces a 3D system to an effective 1D system based on the assumption that a fast equilibration of particles along the transversal directions of the geometry. The transport characteristics of ABPs are computed semi-analytically and corroborated by numerical simulations. At the optimal frequency of the geometry oscillation, particles exhibit higher average velocity $\langle v \rangle$ and effective diffusion coefficient $D_{\text{eff}}$, resembling the phenomena of stochastic resonance. This effect is further enhanced by the self-propelled velocity of ABPs and the amplitude of geometry oscillations. These findings have significant implications for the development of micro- and nanofluidic devices with enhanced control over particle transport and precise manipulation of small-scale biomedical devices.

Figures

Figures reproduced from arXiv: 2509.11147 by Ankit Gupta, P. S. Burada, Rahul Sinha.

Figure 1
Figure 1. Figure 1: FIG. 1. (Left) schematic illustration of the shape a 3D oscillating tube, with time-dependent [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. ( [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Snapshots of active particles (with [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) The magnitude of the mean velocity [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

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

43 extracted references

  1. [1]

    X. Long, L. Chen, X. Xiao, X. Min, Y. Wu, Z. Yang, and X. Wen, Front. Cell Dev. Biol. 12 (2024)

    2024

  2. [2]

    5 10 ⟨v⟩ 0 20 40 60 0 0. 1 0. 2

  3. [3]

    Kumral and A

    D. Kumral and A. M. Zfass, Dig. Dis. Sci. 63, 2500 (2018)

    2018

  4. [4]

    5 a Deff FIG. 4. (a) The magnitude of the mean velocity ⟨v⟩. The left arrows indicate the values of saturated mean velocity ⟨v⟩ = vw = ωL 2π beyond a = b and (b) the effective diffusion coefficient Deff as a function of the corrugation amplitude a for different values of boundary oscillation frequency ω. The right arrows denote the value of Deff calculated...

    2023

  5. [5]

    D. Wang, L. Vahala, and Z. Hao, Phys. Rev. E 98, 032402 (2018)

    2018

  6. [6]

    Marcos, H. C. Fu, T. R. Powers, and R. Stocker, Proc. Natl. Acad. Sci. U.S.A. 109, 4780 (2012)

    2012

  7. [7]

    5 10 v0 ⟨v⟩ 0 10 200 20 40 Deff ⟨v⟩ Deff 10−1 100 101 102 103 104 0 10 20 ω Deff 0 6 12 FIG. 2. ( a) The mean velocity ⟨v⟩ and ( b) the effective diffusion coefficient Deff as a function of the tube oscillation frequency ω for various values of the self-propulsion speed v0 of the particles. Symbols correspond to numerical simulation results (Eqs. 14 & 15),...

  8. [8]

    P. R. Corridon, Sci. Rep. 11 (2021)

    2021

  9. [9]

    2 a = b 0

    55(a) (b) vw = ωL 2π ω vw b a = 0 a = 0. 2 a = b 0

  10. [10]

    Reguera, A

    D. Reguera, A. Luque, P. S. Burada, G. Schmid, J. M. Rub ´ ı, and P. H¨ anggi, Phys. Rev. Lett. 108, 020604 (2012)

    2012

  11. [11]

    P. S. Burada, P. H¨ anggi, F. Marchesoni, G. Schmid, and P. Talkner, ChemPhysChem 10, 45 (2009)

    2009

  12. [12]

    P. S. Burada, G. Schmid, D. Reguera, M. H. Vainstein, J. M. Rubi, and P. H¨ anggi, Phys. Rev. Lett. 101, 130602 (2008)

    2008

  13. [13]

    Gupta and P

    A. Gupta and P. S. Burada, Phys. Rev. E 108, 034605 (2023)

    2023

  14. [14]

    Reguera and J

    D. Reguera and J. M. Rub ´ ı, Phys. Rev. E64, 061106 (2001)

    2001

  15. [15]

    Kalinay and J

    P. Kalinay and J. K. Percus, Phys. Rev. E 74, 041203 (2006)

    2006

  16. [16]

    Malgaretti, I

    P. Malgaretti, I. Pagonabarraga, and J. M. Rubi, Front. Phys. 1 (2013)

    2013

  17. [17]

    Arango-Restrepo, J

    A. Arango-Restrepo, J. M. Rubi, S. Kjelstrup, B. A. J. Angelsen, and C. de Lange Davies, Biophys. J. 120, 5255 (2021)

    2021

  18. [18]

    X. Yang, C. Liu, Y. Li, F. Marchesoni, P. H¨ anggi, and H. P. Zhang, Proc. Natl. Acad. Sci. U.S.A. 114, 9564 (2017)

    2017

  19. [19]

    J. M. Rub ´ ı, A. Lervik, D. Bedeaux, and S. Kjelstrup, J. Chem. Phys. 146 (2017)

    2017

  20. [20]

    Vicsek, A

    T. Vicsek, A. Czir´ ok, E. Ben-Jacob, I. Cohen, and O. Shochet, Phys. Rev. Lett. 75, 1226 (1995)

    1995

  21. [21]

    A. K. Omar, K. Klymko, T. GrandPre, and P. L. Geissler, Phys. Rev. Lett. 126, 188002 (2021)

    2021

  22. [22]

    R. N. Valani, B. Harding, and Y. M. Stokes, Phys. Rev. E 110, 034603 (2024)

    2024

  23. [23]

    Vachier and M

    J. Vachier and M. G. Mazza, Eur. Phys. J. E 42, 1 (2019)

    2019

  24. [24]

    Sandoval and L

    M. Sandoval and L. Dagdug, Phys. Rev. E 90, 062711 (2014)

    2014

  25. [25]

    Reguera, G

    D. Reguera, G. Schmid, P. S. Burada, J. M. Rub ´ ı, P. Reimann, and P. H¨ anggi, Phys. Rev. 12 Lett. 96, 130603 (2006)

    2006

  26. [26]

    K. J. Modica, A. K. Omar, and S. C. Takatori, Soft Matter 19, 1890 (2023)

    2023

  27. [27]

    M. F. Carusela and J. M. Rub ´ ı, J. Chem. Phys.146 (2017)

    2017

  28. [28]

    X. Ao, P. Ghosh, Y. Li, G. Schmid, P. H¨ anggi, and F. Marchesoni, Eur. Phys. J. Spec. Top. 223, 3227 (2014)

    2014

  29. [29]

    Khatri and P

    N. Khatri and P. S. Burada, Phys. Rev. E 105, 024604 (2022)

    2022

  30. [30]

    Volpe, I

    G. Volpe, I. Buttinoni, D. Vogt, H.-J. K¨ ummerer, and C. Bechinger, Soft Matter 7, 8810 (2011)

    2011

  31. [31]

    Gupta and P

    A. Gupta and P. S. Burada, J. Chem. Phys. 161 (2024)

    2024

  32. [32]

    Gammaitoni, P

    L. Gammaitoni, P. H¨ anggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998)

    1998

  33. [33]

    Z¨ ottl and H

    A. Z¨ ottl and H. Stark, J. Phys.: Condens. Matter 28, 253001 (2016)

    2016

  34. [34]

    P. S. Burada, Y. Li, W. Riefler, and G. Schmid, Chem. Phys. 375, 514 (2010)

    2010

  35. [35]

    Y. Wang, D. S. Dean, S. Marbach, and R. Zakine, J. Fluid Mech. 972, A8 (2023)

    2023

  36. [36]

    Sarfati, C

    R. Sarfati, C. P. Calderon, and D. K. Schwartz, ACS nano 15, 7392 (2021)

    2021

  37. [37]

    S. G. Nagella and S. C. Takatori, AIChE Journal 69, e18215 (2023)

    2023

  38. [38]

    Romanczuk, F

    P. Romanczuk, F. M¨ uller, and L. Schimansky-Geier, Phys. Rev. E 81, 061120 (2010)

    2010

  39. [39]

    J. D. Torrenegra-Rico, A. Arango-Restrepo, and J. Rub ´ ı, J. Chem. Phys.156 (2022)

    2022

  40. [40]

    Stubenrauch and R

    C. Stubenrauch and R. v. Klitzing, J. Phys.: Condens. Matter 15, R1197 (2003)

    2003

  41. [41]

    Marbach, D

    S. Marbach, D. S. Dean, and L. Bocquet, Nat. Phys. 14, 1108 (2018)

    2018

  42. [42]

    Bleil, P

    S. Bleil, P. Reimann, and C. Bechinger, Phys. Rev. E 75, 031117 (2007)

    2007

  43. [43]

    S. W. Kowalczyk, D. B. Wells, A. Aksimentiev, and C. Dekker, Nano Lett. 12, 1038 (2012). 13

    2012