{"id":"fe53fcb5-b3cc-4278-b0be-17b5954c5d1e","arxiv_id":"2602.08436","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The stationary position distribution of a run-and-tumble particle in a harmonic trap is solved exactly in 1D, 2D, and 3D; the 3D radial law is new and is not a beta distribution.","lead":"This paper derives exact closed-form formulas for where a run-and-tumble particle spends its time when trapped in a harmonic (bowl-shaped) potential, in one, two, and three dimensions, including with thermal noise. The main new result is the three-dimensional radial and joint position distribution, which is more complex than the beta distributions found in lower dimensions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: Kesten ergodicity is guaranteed and the Dirichlet-process identification is standard.","rationale":"The paper's central claim is well supported. The reduction by isotropy is exact; the Kesten recursion is provably ergodic because E[log U] = -1/α < 0 and V is bounded; the Dirichlet-process representation is standard Sethuraman stick-breaking; and the closed forms for d=1,2 match known results, while the d=3 expressions are internally consistent and agree with the reported simulations. The reader's concern about ergodicity overstates the risk because the conditions of the Kesten theorem hold automatically from the exponential run-time distribution. The remaining limitations—no code/simulation details and reliance on the cited Cifarelli-Regazzini theorem—are presentation gaps, not correctness risks. Therefore the ACCEPT verdict is unchanged.","tokens_in":37878,"tokens_out":22179,"duration_ms":208576,"concrete_test":"For W(v) = 1/(2v0) (the d=3 projected velocity law), compute the moment generating function of the density in Eq. (45) by numerical quadrature for several q and compare it to the Kesten-derived MGF in Eq. (B12) evaluated by contour integration or direct simulation. Agreement to machine precision would confirm the Cifarelli-Regazzini inversion; any mismatch would pinpoint the step where the theorem is misapplied.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I do not find a load-bearing flaw. The reader's weakest assumption is the uniqueness/invariance of the Kesten recursion and the legitimacy of the infinite-limit representation. In this model, U_n = e^{-μτ_n} ~ Beta(α,1), so E[log U] = -1/α < 0 and V_n is bounded; the classical Kesten theorem (Ref. [74]) then guarantees a unique stationary solution, and the infinite sum Σ ∏_{j<m} U_j (v_m/μ)(1-U_m) converges absolutely because the products decay exponentially. The stick-breaking weights are exactly Sethuraman's construction for a Dirichlet process with concentration α and base W; the Cifarelli-Regazzini formula (C4) applies to bounded f(v)=v/μ and finite base measure. The main formulas (45)-(47) are consistent with the known d=1,2 results and with the numerical checks reported in the paper. The only formal gap is that the paper states the Cifarelli-Regazzini theorem without proof, but this is a standard, well-established result, not a risk to the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the exact stationary state of a run-and-tumble particle (RTP) in an isotropic harmonic potential in d dimensions. It first solves a generalized one-dimensional RTP with arbitrary post-tumble velocity distribution W(v) via a Kesten recursion, representing the stationary position as an infinite stick-breaking sum equivalent to a mean functional of a Dirichlet process. Using the Cifarelli-Regazzini inversion, it obtains closed-form expressions for the single-coordinate density (Eq. 45) and moments (Eq. 47). Specializing W(v) to the projected velocity law of an isotropic RTP, the paper reconstructs the radial distribution p_R(r) and the full joint stationary density. In d=1 and d=2 the radial law is a beta distribution, while in d=3 the authors derive a new closed-form p_R(r) (Eq. 83) and joint density P(x,y,z) (Eq. 85) that are not beta. The paper also treats thermal noise D>0, showing that the stationary law is a Gaussian convolution of the D=0 law, and analyzes the resulting shape transitions, including a two-step/discontinuous maximum transition in d=3. All predictions are compared with numerical simulations.","tokens_in":38093,"tokens_out":8829,"duration_ms":96739,"significance":"This is a significant contribution. The d=3 closed-form stationary distribution closes a gap that had remained open despite several recent efforts (Refs. 39-41). The Dirichlet-process/stick-breaking route is elegant and likely transferable to other linear switching systems. The d=2 result correctly reproduces the independent beta law of Frydel, and the d=3 formulas are new and nontrivial. The method yields parameter-free exact expressions and explicit moments via Bell polynomials. The finite-D analysis, including the universal low-D scaling form near the turning surface and the first-order-like jump of the global maximum in d=3, provides concrete testable predictions. The only formal gap—the convergence of the Kesten recursion—is covered by the classical Kesten theorem because E[log U]=-1/alpha<0 and V_n is bounded; I do not regard this as a load-bearing flaw.","major_comments":[],"minor_comments":[{"comment":"The notation 'U_n = 1 - U_n, U_n ~ Beta(1, alpha)' is confusing because the same symbol is used for the original and transformed variables. Suggest using, e.g., V_n or U'_n for the transformed variable to avoid ambiguity.","section":"Sec. III B, Eq. (42)"},{"comment":"The statement that p_X(x) in Eq. (73) 'has been checked numerically' to satisfy the integro-differential equation (74) is a consistency check, not a proof. Since the derivation of Eq. (73) from the general formula is independent, this is acceptable, but it should be labeled as a numerical consistency check rather than a verification of the main result.","section":"Sec. IV C, Eq. (73)"},{"comment":"The Cifarelli-Regazzini identity is quoted without proof. This is acceptable because it is a published theorem (Ref. [75]), but the paper would be more self-contained if it stated the precise theorem and its applicability conditions (e.g., bounded f(v)=v/mu and finite base measure) in one sentence.","section":"Appendix C, Eq. (C5)"},{"comment":"References [79] and [80] appear to have identical titles. Please check whether they are distinct works or whether one citation is erroneous.","section":"References [79] and [80]"},{"comment":"References [82] and [95] are Wikipedia links. For a journal submission, consider replacing them with standard textbook or DLMF citations for Bell polynomials and Appell hypergeometric functions.","section":"References [82] and [95]"},{"comment":"The title and abstract emphasize the 'd-dimensional' solution, but explicit closed-form radial and joint distributions are worked out for d=1,2,3. For d>3, the result is an exact integral representation of the single-coordinate marginal plus explicit moments. This is not a flaw, but the abstract could be more precise by saying 'closed forms in d=1,2,3 and an exact integral representation in general d'.","section":"Abstract and Title"},{"comment":"The figure labels/caption appear to mix the symbols theta_c and mu_c. Please check that the axis labels and critical lines are consistently denoted (theta_c(alpha) for d=1,2).","section":"Fig. 8"}],"recommendation":"minor_revision","confidential_remarks":"I see no substantive correctness issues. The central derivation is sound and the d=3 result is new and significant. The minor points are presentation-level; once addressed, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read the paper carefully. The headline claim is right: 1D and 2D radial beta laws were already known, but the exact d=3 stationary radial and joint densities (Eqs. 83–85) are new, and they do not reduce to a beta law. The general stationary formula for arbitrary post-tumble velocity distribution W(v) via the stick-breaking construction is also a real conceptual step; it unifies known special cases and gives closed-form moments.\n\nThe derivation is mostly self-contained and honest. They take the Kesten recursion, solve it explicitly in the stationary limit, and identify the solution as a mean functional of a Dirichlet process. The Cifarelli–Regazzini formula then gives the density. I checked the chain: the d=2 result reproduces Frydel's beta law, and the d=3 expression is derived from the stated general formula without any obvious circular step. The numerics agree with the analytic curves in the figures.\n\nThe soft spots are minor. The paper assumes the n→∞ limit of the Kesten recursion yields a unique stationary law; it cites Kesten's theorem rather than proving the conditions. In this model that is fine: U_n ∼ Beta(α,1) so E[log U] = −1/α < 0 and V_n is bounded, so the classical theorem guarantees convergence. I'd have liked the paper to say that in one sentence, but it's not a gap in the result. The Cifarelli–Regazzini identity is also used as a black box; standard in the Bayesian nonparametrics literature, so acceptable. One specific validation — Eq. (73) satisfying the integro-differential equation (74) — is numerical rather than analytical; they state this openly. No code or simulation protocols are shipped, which is a minor reproducibility nuisance but not a flaw.\n\nI find no load-bearing objection. The central claim holds up. Who is this for? Statistical physicists working on active matter, exact nonequilibrium steady states, and anyone using Kesten recursions or Dirichlet processes. It's a good paper for a reading group: the method is instructive and the d=3 formulas are new.\n\nRecommendation: send it to peer review. A competent referee can verify the algebra; I expect the paper to be accepted with at most small revisions.","headline":"The d=3 closed-form stationary densities are genuinely new, the Kesten/Dirichlet-process method is clean, and the paper deserves serious peer review.","tokens_in":38630,"tokens_out":2430,"would_cite":true,"duration_ms":25164,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The stationary state of a run-and-tumble particle in a harmonic trap is exactly solvable in any dimension, and in three dimensions takes a closed but non-beta form.","keywords":["run-and-tumble particle","harmonic trap","stationary distribution","Dirichlet process","stick-breaking","Kesten recursion","shape transition","thermal noise"],"falsifier":"Simulate a three-dimensional run-and-tumble particle in a harmonic trap with μ=1, v0=1, and α=0.5, and measure the stationary radial histogram near the turning surface. If p_R(r) does not scale as ε^{α-1} (with log corrections) as ε = 1 - r → 0⁺, the closed form (83) is wrong.","tokens_in":37749,"feed_emoji":"🎯","tokens_out":7405,"duration_ms":64056,"temperature":0.7,"pith_summary":"Run-and-tumble particles confined by a harmonic trap reach a nonequilibrium steady state that, despite decades of study, had no exact description beyond one and special two-dimensional cases. This paper closes that gap by showing that the full stationary distribution in any dimension is encoded in the marginal distribution of a single Cartesian coordinate, which itself follows from a generalized one-dimensional model with arbitrary post-tumble velocities. Solving that model through a Kesten recursion and its stick-breaking (Dirichlet-process) representation gives closed-form densities and moments for every dimension. Specializing to the projected velocities of an isotropic tumbling particle yields explicit radial laws: beta distributions in d=1 and d=2, and a genuinely different closed form in d=3. Adding thermal noise simply convolves these results with a Gaussian, regularizing the turning-point singularities and producing a finite-temperature shape transition.","feed_headline":"Trapped run-and-tumble particles: exact steady state in any dimension","feed_subtitle":"A single coordinate's distribution determines the full radial and joint densities, and thermal noise just smooths the result.","key_machinery":"The central object is the Kesten recursion x_n = U_n x_{n-1} + V_n for the position just after each tumble, with U~Beta(α,1). Unrolling the recursion maps the stationary position to a weighted sum of the post-tumble velocities, whose random weights are exactly the stick-breaking weights of a Dirichlet process of concentration α and base measure W(v). This 'mean functional of a Dirichlet process' representation yields, via a known identity, closed-form densities and moments for arbitrary W(v). Combined with rotational invariance — the radial and joint densities are integral transforms of the single-coordinate marginal — it turns the d-dimensional stationary problem into a one-dimensional calc","core_discovery":"The paper's central claim is that the stationary position of a d-dimensional run-and-tumble particle in an isotropic harmonic trap is fully determined by the one-coordinate marginal p_X(x), and that p_X(x) for any post-tumble velocity law W(v) is exactly the density of a mean functional of a Dirichlet process. This identification turns an intractable nonlocal Fokker-Planck equation into closed-form expressions: p_X(x) is given by an integral formula involving a simple function φ_α, and all moments are Bell polynomials in the moments of W(v). For the isotropic d-dimensional RTP the projected velocity W(v) is the Beta-type law of Eq. (7); the resulting radial distribution is a beta law in d=1","pith_inferences":["The same machinery should extend to heterogeneous run speeds: if v0 is itself drawn from a distribution at each tumble, the single-coordinate formula (16) remains unchanged, so the whole d-dimensional solution carries over with the modified projected velocity law.","Because the Dirichlet-process representation is linear in the velocities, the approach likely also applies to linear switching dynamics, such as Brownian motion with a stochastically switching trap stiffness; this may yield exact two-time or multi-particle correlations.","The d=3 non-beta character suggests that for all d≥3 the radial law is not a beta distribution, with a dimension-dependent family of closed forms; the d=2 beta law may be a special coincidence linked to the planar projection being arcsine.","The finite-D two-step transition in d=3, with a discontinuous jump of the global maximum, resembles a first-order transition and could be probed experimentally in optical-tweezer setups with artificial swimmers."],"forward_implications":["In d=1 and d=2, the radial stationary law is exactly a beta distribution, with a persistence-controlled shape transition at the turning radius.","In d=3, the radial distribution is no longer beta but is given in closed form, and it still exhibits a shape transition at α=1.","Thermal noise D>0 leaves the steady state as a Gaussian convolution of the noise-free law, so all turning-point singularities are rounded and the distribution acquires Gaussian tails beyond r=v0/μ.","The Dirichlet-process representation yields exact stationary states for N-state run-and-tumble models (discrete velocity sets) with piecewise-continuous densities.","The moments of the stationary position are given in closed form for any dimension and any W(v) via Bell polynomials."],"fun_headline_variants":["Trapped run-and-tumble: exact steady state, any dimension","Dirichlet process yields exact radial law for trapped RTP","Beta distribution emerges as exact trap positional law in 1D and 2D","One-coordinate marginal solves trapped run-and-tumble in all d","Thermal noise smooths the sharp turning point: exact crossover"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The derivation assumes that the infinite stick-breaking sum representing the stationary position actually converges to the unique invariant law of the Kesten recursion for every post-tumble velocity distribution W(v); the paper invokes the representation but does not prove this ergodicity step for arbitrary W(v).","fun_headline_variants_meta":{"raw":{"variants":["Trapped run-and-tumble: exact steady state, any dimension","Dirichlet process yields exact radial law for trapped RTP","Beta distribution emerges as exact trap positional law in 1D and 2D","One-coordinate marginal solves trapped run-and-tumble in all d","Thermal noise smooths the sharp turning point: exact crossover"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00028,"raw_usage":{"total_tokens":1572,"prompt_tokens":893,"completion_tokens":679,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":598}},"tokens_in":637,"tokens_out":679,"duration_ms":7709,"temperature":1.0,"reasoning_tokens":598,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:16:20.954343+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a three-dimensional run-and-tumble particle in a harmonic trap with μ=1, v0=1, and α=0.5, and measure the stationary radial histogram near the turning surface. If p_R(r) does not scale as ε^{α-1} (with log corrections) as ε = 1 - r → 0⁺, the closed form (83) is wrong.","supporting_citations":[],"review_version":1}