{"id":"702d2957-2606-4f2e-ba78-ae726ed0990e","arxiv_id":"1908.00568","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Run-and-tumble particles in a one-dimensional random potential show Sinai-type ultra-slow diffusion, and interacting particles form clusters whose mean size scales as the square root of the system size.","lead":"This paper analyzes how random, fixed environments affect the motion of self-propelled 'run-and-tumble' particles in one dimension. It finds that a bumpy energy landscape makes a single particle spread extremely slowly, and many interacting particles clump into clusters whose average size grows with the system size.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Many-body strong-disorder clustering claim rests on an unvalidated lattice-fermion analogy; direct simulations on interacting RTPs are absent.","rationale":"The reader's verdict is CONDITIONAL with medium correctness risk, which is appropriate. The single-particle analysis is careful and the Sinai mapping is supported by analytic derivation and numerical collapse to ln^4 t. The weak-disorder many-body field theory is also checked against RTP simulations. The remaining weak point is precisely the strong-disorder interacting case: the abstract's headline contrast (√L cluster growth vs finite clusters) rests entirely on the lattice-fermion analogy. The paper itself flags the missing direct simulation. My proposed test uses the authors' own Sec. 4 claim that weak translational diffusion does not change the physics; if true, it provides an efficient route to the strong-disorder steady state for actual RTPs, closing the gap. Because the concern is the same as the reader's weakest_assumption, agreement is 'agree'. Since the existing verdict already conditions acceptance on this issue, I do not adjust the verdict.","tokens_in":16778,"tokens_out":6793,"duration_ms":73208,"concrete_test":"Simulate the full interacting-RTP model of Sec. 3 (harmonic repulsion, random ratchet potential of Fig. 2) with a small translational diffusion D added, as the authors argue in Sec. 4 that noise-assisted barrier hopping does not change the physics. Choose D small enough that the strong-disorder length l* remains large but the steady state is reachable in feasible time; measure the cluster-size distribution and mean cluster size for L≫l*. If P(l) decays as l^{-3/2} and ⟨l⟩∼√L, the lattice-fermion analogy is validated for actual RTPs; if the exponent or scaling changes, the central many-body claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2 derives the headline many-body result—power-law cluster-size distribution P(l)∼l^{-3/2} and mean cluster size ∼√L—from a lattice model of non-interacting fermions in a random ±1 force landscape (Eqs. 32–33, Figs. 10–11), not from simulations of interacting RTPs. The paper states these simulations are 'prohibitively slow' and relies on the analogy between RTPs and passive particles in a random force field. That analogy is established for the single-particle steady state and MFPT (Secs. 2.1.1, 2.2.1), but the many-body interacting problem is not equivalent by any proven mapping: active particles with repulsive interactions need not minimize a free energy, and the effective quasi-potential U(x) of Eq. (6) is derived for an isolated particle. The first-return argument giving l^{-3/2} is a property of the equilibrium lattice fermion model; whether it survives continuous-space RTPs with finite-range repulsive interactions is untested. Because the abstract's central claim is the contrast between √L scaling and finite clusters in clean systems, this missing validation is the most load-bearing weakness. The weak-disorder field theory was checked against RTP numerics (Fig. 8), but no analogous check exists in the strong-disorder regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the effects of quenched disorder on run-and-tumble particles (RTPs) in one dimension. For a single particle, it derives the steady-state distribution (Eq. 3) and the mean first passage time to exit an interval for three types of disorder: a bounded external potential, a space-dependent speed, and a space-dependent tumbling rate. For disordered potential, the steady-state distribution is argued to be equivalent to that of a passive particle in a random force field, and the typical MFPT grows as exp(sqrt(L)), implying Sinai-type diffusion with mean-square displacement scaling as ln^4(t) (Eq. 13), verified numerically in Fig. 3. For speed and tumbling-rate disorder, normal diffusion is generally found, with anomalous exponents when the speed distribution vanishes as v^{-beta} or the tumbling-rate distribution has a diverging mean. For many particles in a disordered potential, the paper predicts, in the weak-disorder regime, a structure factor S(q) ~ q^{-2} (Eq. 29), and in the strong-disorder regime, a cluster-size distribution P(l) ~ l^{-3/2} with mean cluster size ~ sqrt(L) (Eqs. 32-33), the latter derived from a lattice model of non-interacting fermions in a random force landscape.","tokens_in":17072,"tokens_out":7088,"duration_ms":67180,"significance":"The single-particle analysis is a genuine contribution: the derivations are explicit and the Sinai-diffusion prediction for RTPs is non-trivial and is tested by direct simulation without fitting the exponent. The weak-disorder many-body result is also compared against RTP numerics. However, the strong-disorder many-body clustering claim, which is highlighted in the abstract, is supported only by a heuristic mapping to an equilibrium lattice fermion model and has no direct simulation of interacting RTPs; this is the main weakness. If the mapping can be substantiated, the paper would make a strong contribution; as it stands, the many-body part is more suggestive than established.","major_comments":[{"comment":"The central many-body strong-disorder result, the power-law cluster-size distribution P(l) ~ l^{-3/2} and mean cluster size ~ sqrt(L), is derived entirely from a lattice model of non-interacting fermions in a random +/-1 force landscape. The paper states (Sec. 3.2) that direct RTP simulations are 'prohibitively slow' in this regime, and no test of the mapping against interacting RTP dynamics is provided. The single-particle equivalence established in Sec. 2.1.1 concerns the steady state of an isolated particle; it does not by itself imply that soft repulsive RTPs obey the same quasi-potential statistics once interactions are present. Because this is the headline claim of the abstract, please either add a direct numerical test of the cluster-size scaling for interacting RTPs in a numerically accessible parameter range, or provide a physically motivated argument (beyond the analogy) for why the fermion lattice model is the correct effective description.","section":"Sec. 3.2, Eqs. (32)-(33), Figs. 10-11"},{"comment":"The mapping to non-interacting fermions assumes (i) that the many-particle steady state is described by filling the single-particle quasi-potential U(x) of Eq. (6) up to a chemical potential and (ii) that finite-range repulsion can be replaced by hard-core exclusion on a lattice. For RTPs with soft pairwise repulsion, neither assumption is self-evident, since the system is out of equilibrium and U(x) was derived from the no-current condition for a single particle. Fig. 10 tests only the lattice fermion model, so the l^{-3/2} law and the sqrt(L) scaling are properties of that model, not of the RTP system. The authors should either relax the claim or provide evidence that the effective interaction among trapped RTPs in the strong-disorder regime is indeed hard-core and that the quasi-potential picture survives interactions.","section":"Sec. 3.2, mapping to non-interacting fermions"}],"minor_comments":[{"comment":"The numerical verification of Eq. (13) would be more convincing with error bars or a statement of the number of disorder realizations used; the caption also appears to read '104 ratchets' where '10^4 ratchets' is intended.","section":"Fig. 3"},{"comment":"Please state explicitly what the constant A in Eq. (12) depends on (e.g., disorder strength and correlation length) and how its value is chosen in the rescaled axes of Fig. 3.","section":"Sec. 2.2.1, Eq. (12)"},{"comment":"The structure factor is evaluated for q != 0 in Eq. (29); please state the treatment of the q=0 mode in the Gaussian free energy (21), since the normalization of the total density is sensitive to this mode.","section":"Sec. 3.1, Eq. (29)"},{"comment":"The closing statement that any 1d active particle model generating ratchet currents 'should exhibit the same phenomenology' goes beyond the models analyzed; consider presenting this as a conjecture rather than a conclusion.","section":"Sec. 4"},{"comment":"There are minor typographical issues, e.g., 'dependance' in Sec. 2.2.1 and 'NSF5-BSF' in the acknowledgments, which should be corrected.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The single-particle section is strong and likely publishable on its own. The main risk is that the many-body strong-disorder claim, which is central to the paper's stated conclusions, is under-supported. I would ask the authors for a numerical test of interacting RTPs (e.g., in the ratchet potential used in Fig. 8 but with weaker disorder and larger L, or with harder interactions) or a substantial theoretical argument for the fermion mapping. If this cannot be provided, the paper should be reframed as presenting the fermion model as a conjecture. The authors should also double-check the novelty of the weak-disorder structure factor result against existing literature on Gaussian free energies with Brownian disorder, since Eq. (29) follows from standard Gaussian integrals."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The single-particle part is the real contribution. The mapping from an RTP in a bounded random potential to an equilibrium walker in a random force field is convincing: the steady-state quasi-potential argument is explicit, the MFPT calculation is spelled out, and the predicted ln^4 t Sinai spreading is checked against direct RTP simulations. The anomalous diffusion exponents for singular speed and tumbling-rate distributions also look new and are backed by numerics. The weak-disorder many-body field theory is standard but competently done, and the q^-2 structure factor is verified with interacting RTPs. Credit where due: the paper is clearly written, the algebra in the appendices is checkable, and the authors are upfront about what is simulation and what is analogy.\n\nThe soft spot is exactly the one the stress test flags. Section 3.2's strong-disorder interacting result—power-law cluster sizes P(l) ~ l^-3/2 and mean cluster size ~ sqrt(L)—is derived from non-interacting lattice fermions in a random ±1 force landscape, not from RTPs. The paper says direct RTP simulations are 'prohibitively slow' and then presents only lattice-model data. That is a transparent statement, but it leaves the abstract's central contrast with finite clusters in clean systems unsupported for the actual model. The single-particle Sinai analogy does not automatically carry over to a many-body interacting system: active particles with repulsive interactions are not minimizing a free energy, and the quasi-potential U(x) is derived for an isolated particle. The first-return argument for l^-3/2 is a property of the equilibrium fermion model. It may survive in the RTP system, but it is currently a conjecture. The missing error bars in the numerics are a minor issue; the unvalidated many-body strong-disorder claim is the load-bearing one.\n\nWho is this for? Statistical physicists working on active matter in disordered environments. The single-particle results and the weak-disorder calculation are solid and useful. The strong-disorder many-body part should be framed as a conjecture or supported by a direct simulation of interacting RTPs at accessible parameters, or by a sharper analytic argument connecting the RTP dynamics to the fermion model. I would send this to a serious referee: the paper deserves peer review, and a referee can push on Section 3.2. With that section revised or explicitly demoted, I'd be happy to cite it.","headline":"A mostly careful paper with a genuinely nice Sinai-diffusion mapping for single RTPs, but the headline many-body strong-disorder clustering claim is a conjecture resting on a lattice-fermion analog that is not validated for interacting RTPs.","tokens_in":17519,"tokens_out":952,"would_cite":true,"duration_ms":11791,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C31","82C44","60J60","82C22"],"pacs":["05.40.-a","05.40.Fb","05.60.-k"],"model":"deepseek-v4-flash","headline":"This paper aims to establish that a run-and-tumble particle in a bounded random potential behaves at long times like an equilibrium particle in a random force field, spreading as $\\langle x^2(t)\\rangle \\propto \\ln^4 t$, and that many…","keywords":["run-and-tumble particles","quenched disorder","Sinai diffusion","random force landscape","active matter","one-dimensional systems","clustering","mean first passage time"],"falsifier":"Simulate many interacting run-and-tumble particles with hard-core or short-range repulsion in a strongly disordered random ratchet potential, at system sizes far beyond the crossover length $\\ell^*$, and measure the mean cluster size as the system size doubles; if it does not grow as $\\sqrt{L}$, or if the cluster-size distribution is not a power law, the central many-body claim fails. For the single-particle claim, measure the disorder-averaged mean-square displacement in the same landscape and check whether it approaches $\\ln^4 t$ rather than any power of $t$.","tokens_in":16602,"feed_emoji":"🎲","tokens_out":11884,"duration_ms":114468,"temperature":0.7,"pith_summary":"This paper wants to show that quenched spatial disorder changes the physics of run-and-tumble particles in one dimension, and that the change depends on which parameter is disordered. Its central case is a bounded random potential: a particle that runs straight and randomly reverses direction becomes, at long times, equivalent to an equilibrium particle diffusing in a random force field. That equivalence gives strongly localized steady states and ultra-slow spreading, $\\langle x^2(t)\\rangle \\propto \\ln^4 t$, and it turns the finite clusters known in one-dimensional active systems into clusters whose mean size grows as $\\sqrt{L}$. Disorder in the speed or tumbling rate is instead shown to be generically diffusive, with anomalous exponents only when the relevant distribution has diverging moments.","feed_headline":"Active particles in random potentials spread only as log⁴ time","feed_subtitle":"A single particle acts like a walker in a random force field; many interacting ones form clusters that grow with the system size.","key_machinery":"The object carrying the argument is the effective quasi-potential $U(x)\\simeq \\frac{1}{2}\\ln\\bigl[1-(\\partial_x\\tilde V)^2\\bigr]+g\\sum_{i=1}^{\\lfloor x/\\xi\\rfloor}\\eta_i$, built from the exact steady-state density (Eqs. 4-6), where $\\tilde V=\\mu V/v$, $g=\\alpha/v$, and the $\\eta_i$ are independent random increments contributed by each correlation-length cell of the potential. Because the increments are independent and identically distributed, $U$ looks like a random walk in space and grows as $\\sqrt{L}$; this is the feature that both localizes the steady state and sets the exponential first-passage-time scale. For the many-particle strong-disorder result, the second piece of machinery is the mapping of hard-core interacting particles in such a landscape to non-overlapping particles on a lattice, where a cluster is a consecutive filled run and the cluster-size distribution follows from the statistics of the first time a random walk returns to the level that separates filled from empty sites.","core_discovery":"The central claim is a mapping: a run-and-tumble particle with constant speed and tumbling rate in a bounded random potential $V(x)$ has an exact steady state whose exponent is a sum of independent random contributions, so its effective quasi-potential is that of a random force. When the system is much larger than the correlation length of $V$, the effective potential grows as $\\sqrt{L}$ by the central limit theorem, and the first-passage-time statistics match Sinai diffusion: $\\ln\\langle\\tau\\rangle\\sim A\\sqrt{L}$, hence $\\langle x^2(t)\\rangle\\propto\\ln^4 t$. Using this single-particle equivalence, the paper derives that in weak disorder the density structure factor diverges as $q^{-2}$ and real-space correlations decay linearly with $r/L$ with an amplitude linear in $L$; in strong disorder, non-interacting particles collapse around one minimum, while hard-core interacting particles, mapped to non-overlapping particles on a random-force lattice, have cluster sizes distributed as $P(\\ell)\\sim \\ell^{-3/2}$, so the mean cluster size is $\\langle\\ell\\rangle\\sim\\sqrt{L}$. The paper also shows that disorder in the speed or tumbling rate generically leaves ordinary diffusion, except for singular distributions that produce power-law anomalous exponents.","pith_inferences":["The random-force equivalence is probably generic to any one-dimensional active dynamics that rectifies asymmetric potentials, not only run-and-tumble motion; the paper suggests this for Brownian-noise-assisted hopping and other active models, but does not prove it.","In two dimensions the same effective random forces should create circulating currents and a slower-than-diffusive spread $\\langle x^2\\rangle\\sim t/\\ln t$; the paper flags this as an expectation, not a derived result.","A direct numerical test of the strong-disorder many-body claim would simulate many interacting run-and-tumble particles in a strongly disordered ratchet potential at $L\\gg\\ell^*$ and check whether the cluster-size distribution and the $\\sqrt{L}$ scaling survive outside the lattice mapping.","If the clustering result holds, it gives a clean way to distinguish disorder-induced clustering from ordinary motility-induced phase separation: the mean cluster size has no finite thermodynamic limit and grows with system size, which could be probed by controlling speckle or ratchet-like landscapes experimentally."],"forward_implications":["A single active particle in a bounded random potential is strongly localized in its steady state, with the location of the maximum depending on the whole potential profile across the system, unlike a passive particle in a bounded potential.","The disorder-averaged mean-square displacement grows only as $\\ln^4 t$ at long times, with a crossover from ordinary diffusion set by a length scale $\\ell^*\\simeq D^2/\\sigma^2$.","At weak disorder, many-particle systems develop long-range density correlations: the structure factor diverges as $q^{-2}$ and the real-space correlation function decays linearly on the scale of the system size.","At strong disorder, non-interacting particles all gather at one minimum, while interacting particles form clusters with a power-law size distribution and mean size $\\sqrt{L}$, in contrast to the finite clusters found without disorder.","Disorder in speed or tumbling rate generically gives diffusive spreading; anomalous exponents appear only when the speed distribution near zero behaves like $v^{-\\beta}$ or the tumbling-rate tail behaves like $\\alpha^{-(1+\\mu)}$."],"supporting_citations":[{"why":"establishes that active particles experience a net force from asymmetric potentials, the ratchet mechanism behind the random-force equivalence.","marker":"[9]"},{"why":"supplies the exact steady-state solution of the run-and-tumble Fokker-Planck equation used to write Eq. (3).","marker":"[5]"},{"why":"gives the activation-rate framework for nonlinear stochastic flows from which the same steady-state and first-passage calculations follow.","marker":"[17]"},{"why":"provides the Sinai-diffusion result whose first-passage-time scaling the paper reproduces for the active particle.","marker":"[18]"},{"why":"defines the random-force landscape theory, including the logarithmic spreading and the crossover length scale used in the paper.","marker":"[19]"},{"why":"documents the finite-size clusters that one-dimensional active systems form without disorder, the baseline against which the sqrt(L) clustering is contrasted.","marker":"[22]"},{"why":"establishes the statistical mechanics of interacting run-and-tumble particles and their cluster formation without disorder.","marker":"[23]"}],"fun_headline_variants":["Random potentials make active particles move logarithmically slow","Quenched disorder induces strong clustering in active particles","Disorder turns active particle motion into Sinai slow spread","Active particles in disorder: Sinai diffusion and linear-size clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction that clusters grow as the square root of the system size rests on replacing the interacting active particles by non-overlapping particles on a random energy landscape; if this analogy misses the non-equilibrium dynamics of run-and-tumble motion, the predicted power-law clusters would not form.","fun_headline_variants_meta":{"raw":{"variants":["Random potentials make active particles move logarithmically slow","Quenched disorder induces strong clustering in active particles","Disorder turns active particle motion into Sinai slow spread","Active particles in disorder: Sinai diffusion and linear-size clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000894,"raw_usage":{"total_tokens":3867,"prompt_tokens":974,"completion_tokens":2893,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":2831}},"tokens_in":590,"tokens_out":2893,"duration_ms":22529,"temperature":1.0,"reasoning_tokens":2831,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:46:33.209164+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate many interacting run-and-tumble particles with hard-core or short-range repulsion in a strongly disordered random ratchet potential, at system sizes far beyond the crossover length $\\ell^*$, and measure the mean cluster size as the system size doubles; if it does not grow as $\\sqrt{L}$, or if the cluster-size distribution is not a power law, the central many-body claim fails. For the single-particle claim, measure the disorder-averaged mean-square displacement in the same landscape and check whether it approaches $\\ln^4 t$ rather than any power of $t$.","supporting_citations":[{"cited_title":"Active particles with soft and curved walls: Equation of state, ratchets, and instabilities.Physical Review Letters, 117:098001–5, 2016","cited_arxiv_id":null,"evidence_quote":"establishes that active particles experience a net force from asymmetric potentials, the ratchet mechanism behind the random-force equivalence."},{"cited_title":"Pressure is not a state function for generic active fluids.Nature Physics, 11:673–678, 2015","cited_arxiv_id":null,"evidence_quote":"supplies the exact steady-state solution of the run-and-tumble Fokker-Planck equation used to write Eq. (3)."},{"cited_title":"Activation rates for nonlinear stochastic flows driven by non-gaussian noise","cited_arxiv_id":null,"evidence_quote":"gives the activation-rate framework for nonlinear stochastic flows from which the same steady-state and first-passage calculations follow."},{"cited_title":"Limit behaviour of one-dimensional random walks in random envi- ronments","cited_arxiv_id":null,"evidence_quote":"provides the Sinai-diffusion result whose first-passage-time scaling the paper reproduces for the active particle."},{"cited_title":"Classical diffusion of a particle in a one-dimensional random force field.Annals of Physics, 201(2):285–341, 1990","cited_arxiv_id":null,"evidence_quote":"defines the random-force landscape theory, including the logarithmic spreading and the crossover length scale used in the paper."},{"cited_title":"Motility-induced phase separation.Annual Review of Condensed Matter Physics, 6(1):219–244, 2015","cited_arxiv_id":null,"evidence_quote":"documents the finite-size clusters that one-dimensional active systems form without disorder, the baseline against which the sqrt(L) clustering is contrasted."},{"cited_title":"Statistical mechanics of interacting run-and-tumble bacteria.Phys- ical review letters, 100(21):218103, 2008","cited_arxiv_id":null,"evidence_quote":"establishes the statistical mechanics of interacting run-and-tumble particles and their cluster formation without disorder."}],"review_version":1}