{"id":"5f3f6b46-891f-43cf-97be-cce0e3d86ebf","arxiv_id":"2607.22222","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a harmonically trapped active Brownian particle, the regime with a Boltzmann-like position distribution is activity-dominated, while the bimodal regime is equivalent to a passive trapped particle at a displaced position.","lead":"This paper studies a small self-propelled particle trapped in a harmonic (spring-like) potential and asks whether the shape of its position distribution tells you whether it is still actively moving. It finds that a bell-shaped (Boltzmann-like) distribution actually hides strong residual activity, while the bimodal 'active-looking' distribution corresponds to dynamics essentially the same as a passive trapped particle displaced away from the trap center.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Regime-II 'devoid of residual activity' overstates the model's own slow azimuthal drift and is not directly validated by accessible experimental evidence.","rationale":"The reader's weakest_assumption correctly pinpoints the τ_R→∞ idealization in Appendix B, Case-2 and the paper's own admission of a 'marginally widening' vθ_res in regime II. That is the softest load-bearing point in the argument: the advertised two-regime dichotomy (activity-dominated vs. activity-depleted HBBP-like) is stronger than the asymptotic derivation supports. I agree with the reader that the central claim is defensible in its asymptotic limits and that the clean dichotomy is softened by the azimuthal drift. The reader's classification is CONDITIONAL with medium confidence; I reviewed the full text to see whether that concern actually rises to a rejection or whether it is a matter of overstatement. I find the underlying analytical results—Gaussian broadening in regime I, displaced Gaussian/bimodal in regime II—are correct in the stated limits and well supported by the simulations shown, including the clear HBBP-like collapse of the PSDs (Fig. 7(b)) and the V-linear rc with unchanged keff (Fig. 8). The flaw is not that the asymptotic analysis is wrong; it is that the categorical language 'devoid of residual activity' and 'fully devoid of activity' (Sec. III B 5) outruns the finite-τR results, since even in regime II the azimuthal propulsion term Vθ in Eq. 11 contributes to vθ_res. That concern is already explicitly asserted by the paper itself in Fig. 6(c) and Supp. Video 2, which is strong in-scope evidence for the reader's reading. Given the reader already flagged this exact assumption, the agreement is full. My only difference from the reader is that I would keep the verdict at CONDITIONAL rather than moving it further: the concern is a qualification of the headline claim, not a refutation, because the central crossover story (timescale ratio governs which signature appears) is robust. The requested concrete test is one quantitative check that would settle whether the azimuthal widening is physically negligible or whether regime II retains a measurable residual-activity signature in the tangential channel. If the test shows the excess is small relative to the Brownian channel, the current conditional acceptance stands on strengthened ground; if not, the paper should explicitly downgrade 'devoid of residual activity' to 'radial residual activity is balanced, while azimuthal drift remains.' I do not see a second concern of comparable weight: the regime-I claim (activity-dominated with Boltzmann shape) is well supported by the variance, residual-velocity, MSD, and PSD evidence; the intermediate-regime treatment is honestly labeled as lacking closed-form position distributions; the companion Letter [35] is a real limitation for experimental validation, but the paper's own simulations and analytics carry the argument; and the derivation in Appendix A/B/C, while compact, is internally consistent at the level needed here. Therefore the verdict is UNCHANGED relative to the reader, with the agreement on the weakest assumption explicit.","tokens_in":23656,"tokens_out":2663,"duration_ms":23318,"concrete_test":"Repeat the regime-II residual-velocity analysis at fixed τR ≫ τk (e.g., τR=100 s, τk=1 s) and compute the excess azimuthal mean-square residual velocity Δθ(V) = ⟨(vθ_res)²⟩_V − ⟨(vθ_res)²⟩_{V=0} over the steady-state annulus, for V = 2, 5, 10 μm/s (and, if feasible, τR=300 s to test the asymptotic trend). Then either (a) show Δθ(V) stays bounded by the Brownian contribution (e.g., ≤ 10% of the HBBP value) over the entire V range, or (b) quantify the fractional contribution of the azimuthal propulsion term to the total residual kinetic energy. If Δθ(V) is not negligible, the categorical 'devoid of residual activity' statement requires qualification; if it is negligible, the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central dichotomy—Boltzmann-shaped position distribution equals activity-dominated dynamics, bimodal equals HBBP-like and 'devoid of residual activity'—depends on treating regime II (τR ≫ τk) with the orientation frozen during the equilibration window. Appendix B, Case-2 explicitly sets τR → ∞ and sin(ϕ(t)) = 0, giving P(x) shifted by V τk and P(y) unshifted (Eqs. 6–9). But the actual steady state in regime II includes τR = 100 s, so the ABP slowly reorients and acquires a nonzero azimuthal velocity component. The paper's own Fig. 6(c) shows P(vθ_res) 'marginally widens with increasing propulsion speed' (Sec. III B 2), and Supp. Video 2 describes a 'free and slow dynamics along the azimuthal direction.' That azimuthal drift is activity-driven (Vθ feeds vθ_res in Eq. 11); it is not present for the V = 0 HBBP. If slow azimuthal propulsion counts as residual activity, then the claim that regime-II dynamics is 'essentially reduced to' HBBP dynamics is softened, exactly as the reader's weakest_assumption notes. The bimodal distribution and HBBP-like radial fluctuations are consistent with the asymptotic picture, but the absolute statement 'devoid of residual activity' is not entailed by the equations: it rests on the τR→∞ idealization rather than on a quantitative demonstration that the azimuthal residual-velocity contribution is negligible in the plotted regime (e.g., V=10 μm/s, τR=100 s). Because the companion Letter [35] is not accessible in this manuscript, the simulation data in Fig. 6(c) are the only evidence on this point; the paper does not analyze how ⟨(vθ_res)²⟩ − ⟨(vθ_res)²⟩_{V=0} scales with Vτk or τk/τR. The central claim would still hold in a weaker form, but it is stated categorically.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes a two-dimensional overdamped active Brownian particle in an isotropic harmonic trap (HBABP), combining analytic limiting calculations with Langevin dynamics simulations in three regimes: τR ≪ τk, τR ≫ τk, and τR = τk. It derives closed-form position distributions, the MSD, effective harmonic confinement parameters, and compares simulated residual velocity components, MSDs, and power spectral densities with analytic predictions. The central claim is that the shape of the steady-state position distribution is not a reliable signature of activity: the regime with a Boltzmann-like position distribution (τR ≪ τk) is actually dominated by residual activity, whereas the regime with a bimodal position distribution (τR ≫ τk) is claimed to be effectively a harmonically bound Brownian particle at a displaced position, 'devoid of residual activity.'","tokens_in":24110,"tokens_out":8279,"duration_ms":79780,"significance":"If the central claim holds, the paper provides a useful correction to the common interpretation of the Boltzmann-to-bimodal crossover in trapped active particles and offers several experimentally accessible diagnostics: the widening of the position distribution with propulsion speed, the residual-velocity distributions, the ballistic MSD rise, the activity-dependent effective stiffness, the displaced harmonic center, and the collapse of the PSD onto the HBBP Lorentzian. The analytical limiting derivations in Appendices A–C are standard and the simulations agree with them in the three representative regimes. The paper also makes falsifiable predictions, e.g., keff decreases with Pe in regime I and rc = Vτk with unchanged keff in regime II.","major_comments":[{"comment":"The claim that regime II is 'devoid of residual activity' is stronger than what the model equations imply. At the simulated parameter values (τR = 100 s, τk = 1 s), the orientation is not frozen: Eq. (11) contains the active azimuthal term Vθ, and Appendix B Case-2 obtains the displaced HBBP solution only by taking τR → ∞ and setting sinϕ(t) = 0 for t ≲ τR. The paper itself acknowledges a 'slow azimuthal drift' (Supp. Video 2) and a marginal widening of P(vθ_res) with V (Fig. 6(c)). This azimuthal drift is activity-driven and absent for V = 0, and it is removed from the PSD analysis by computing spectra from radial strips (Fig. 7(b)). To support the absolute 'devoid' statement, the authors need a quantitative bound on the azimuthal residual-velocity contribution, e.g., ⟨vθ_res²⟩_V − ⟨vθ_res²⟩_{V=0} compared with the HBBP Brownian contribution at the largest V used. Without such a bound,","section":"§III B 2, Appendix B Case-2, Eq. (11)"},{"comment":"The observable P(r) is not unambiguously defined. Eq. (5) is the 2D joint density P(x)P(y); the true marginal radial density for the Boltzmann-like regime is (r/σ²) exp(−r²/(2σ²)), not the zero-centered Gaussian in Eq. (5). Likewise, Eq. (10) omits the polar Jacobian for a marginal radial distribution, although it is a local Gaussian if P is measured along a narrow radial strip. The text in §II B mentions trajectory segments in radial strips, but the captions of Figs. 2(d) and 5(d) say 'radial position distributions P(r)' and fit them to Eqs. (5)/(10). This ambiguity matters because the Boltzmann-inversion fits for keff and rc use this quantity. Please define P(r) precisely — marginal radial density versus conditional density along a radial strip — and display the correct normalization or curve in each figure.","section":"Eqs. (5), (10); Figs. 2(d), 5(d)"}],"minor_comments":[{"comment":"The rms-sum identity ⟨(vr_res)²⟩ + ⟨(vθ_res)²⟩ ≈ V² is introduced without an analytical estimate or error bars. A short derivation or an uncertainty statement would make the claim quantitatively robust.","section":"§III A 2, Fig. 3"},{"comment":"The super-Gaussian shape parameter n is a fittable exponent; its reported values (n = 1.44 and 9.3) are descriptive. Please state explicitly that n is not derived and report fit uncertainties.","section":"§III C 1, Fig. 10"},{"comment":"The velocity averaging window of 0.1 s is a free parameter in the residual-velocity analysis. It is at least ten times shorter than τR and τk, so the choice is reasonable, but the sensitivity of P(vr_res) and P(vθ_res) to this window should be reported in an appendix or supplement.","section":"§II B"}],"recommendation":"major_revision","confidential_remarks":"The paper relies on a companion Letter for experimental validation, but the analytical and numerical core stands alone. The main fixable issue is the absolute 'devoid of residual activity' claim in regime II; providing a quantitative azimuthal residual-velocity bound or qualifying the claim should be sufficient. The P(r) normalization ambiguity also needs to be resolved, but it is likely a presentation issue rather than a fundamental flaw."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper deserves a serious referee. Its central claim—that in a harmonic trap the shape of the steady-state position distribution is not a reliable activity signature—is correct and well-supported. The diagnostics it offers (residual velocity distributions, MSD structure, PSD collapse, effective confinement) are sensible and would be useful to anyone studying confined active particles. The paper is honest that most closed-form pieces are known: the MSD expression is cited to ref. [37] and the position distributions are standard. The new contribution is interpretive: Boltzmann-like does not mean passive, and bimodal does not mean strongly active. That correction is worth making.\n\nThe analytical limiting derivations are standard, and the simulations match them in all three regimes. The effective-confinement construction via Boltzmann inversion and the PSD collapse in regime II are particularly nice diagnostics.\n\nWhere it overreaches is the categorical 'devoid of residual activity' for regime II. Appendix B, Case-2 freezes the orientation (τR→∞), but the actual simulated regime has τR=100 s, τk=1 s. During that time the orientation slowly diffuses, producing azimuthal motion. The paper's own Fig. 6(c) shows the azimuthal residual-velocity distribution widening with V, and Supp. Video 2 describes 'free and slow dynamics along the azimuthal direction.' That is activity-driven, not present for the V=0 HBBP. So 'devoid' is not entailed; 'largely suppressed' or 'not visible on the radial diagnostic' would match the evidence. The bimodal distribution and radial HBBP-like behavior are consistent with that weaker version, and the stress-test note is right that the central idea survives.\n\nMinor concerns: the experimental validation is relegated to an inaccessible companion Letter [35], and no code or data are shipped. The simulation data in Fig. 6(c) are the only evidence on the azimuthal point, but the paper does not quantify how the azimuthal velocity variance scales with Vτk or τk/τR. That is a fixable gap.\n\nBottom line: a clear, careful study with an overstated slogan. Send it to peer review; ask for a softened regime-II claim and a quantitative treatment of the azimuthal residual-velocity contribution. I would cite it for the diagnostic package.","headline":"A useful corrective to the 'distribution shape = activity' reading of trapped ABPs, but the categorical 'devoid of residual activity' claim for regime II overstates what the paper's own equations and simulations show.","tokens_in":24606,"tokens_out":2493,"would_cite":true,"duration_ms":24749,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A trapped active particle's position distribution lies about its activity","keywords":["active Brownian particle","harmonic trap","residual activity","position distribution","power spectral density","mean square displacement","effective confinement","persistence time"],"falsifier":"Measure the radial power spectrum of an HBABP with τR≫τk at several propulsion speeds, from V=0 up to V=10 μm/s; the paper predicts all curves collapse onto the passive Lorentzian. If the PSD changes with V, or if the measured effective stiffness from Boltzmann inversion changes with V, the claim that residual activity vanishes in this regime is refuted. The complementary check in regime I is that keff must follow k/(1+Pe²/2).","tokens_in":23538,"feed_emoji":"🔬","tokens_out":8779,"duration_ms":87420,"temperature":0.7,"pith_summary":"This paper asks whether the shape of the steady-state position distribution of a self-propelled particle in a harmonic trap can tell whether the particle is still actively swimming or has been effectively tamed. The answer it reaches is no: in the regime where the distribution is Boltzmann-like (persistence time much shorter than the trap equilibration time), the trapped particle is still dominated by residual activity—the distribution broadens with propulsion speed, the residual velocity components add up to the full propulsion speed, and the power spectrum departs from the equilibrium prediction. In the opposite regime, where the distribution is bimodal or annular (persistence time much longer than trap equilibration), the particle's motion becomes statistically identical to a passive Brownian particle in a displaced harmonic well, with the propulsion exactly balanced by the restoring force at a distance Vτk. The paper therefore replaces the common 'passive vs active' reading of the distribution crossover with a timescale-controlled crossover from activity-dominated to activity-depleted dynamics.","feed_headline":"Boltzmann-like trap hides activity; bimodal trap hides none","feed_subtitle":"Position distributions can't tell if a trapped swimmer is active; the timescale ratio can.","key_machinery":"The central object is the residual velocity vres of the confined particle, i.e., the actual velocity of the trap-bound particle after the restoring force acts. Decomposed into radial and azimuthal components, vres directly reveals whether propulsion survives confinement: in regime I it retains the full propulsion speed, in regime II it reduces to the Brownian velocity of a passive particle. The argument is carried by the competition between the two timescales—the persistence time τR=1/DR and the trap equilibration time τk=γT/k—and by the pair of closed-form position distributions derived in the appendix under τR→0 or τR→∞ limits, together with the effective harmonic confinement parameters (k","core_discovery":"The central claim is that the resultant dynamics of a harmonically bound active Brownian particle (HBABP) are controlled by the ratio of the orientation persistence time τR to the trap equilibration time τk, and that the shape of the steady-state position distribution is not a reliable activity diagnostic. For τR≪τk, the distribution is Gaussian/Boltzmann-like but activity-dominated: its variance grows as (kBT/k)(1+Pe²/2), residual radial and azimuthal velocities show bimodal distributions whose mean-square values sum to V², the MSD shows a ballistic rise, and the PSD deviates from fluctuation-dissipation. For τR≫τk, the 2D distribution is an annulus that looks bimodal in projection, yet the","pith_inferences":["The paper's own Fig. 6(c) shows that the azimuthal residual velocity distribution in regime II widens slightly with propulsion speed; this suggests 'devoid of residual activity' is exact only in the τR→∞ limit, and on timescales approaching τR a weak orientational-activity channel remains—a point worth testing by measuring the azimuthal PSD.","The timescale-ratio criterion likely carries over to other self-propelled models in harmonic traps, such as run-and-tumble or Ornstein-Uhlenbeck active particles, where a similar residual-velocity decomposition could expose hidden activity beneath Boltzmann-like position distributions.","A practical extension for single-particle experiments: rather than fitting only the position distribution, collect simultaneous orientation and position trajectories and compute vres(t); the sum of mean-square radial and azimuthal residual velocities should equal V² in regime I and only the Brownian value in regime II, giving a direct per-trajectory activity meter."],"forward_implications":["In weak traps or fast reorienting swimmers (τR≪τk), a Gaussian/Boltzmann-like position distribution is consistent with strong residual activity, so experiments cannot infer passivity from distribution shape alone.","In strong traps or slow reorienting swimmers (τR≫τk), an annular/bimodal distribution does not imply far-from-equilibrium dynamics; the particle is effectively a passive Brownian particle in a shifted well at r=Vτk.","Peclet number Pe alone is not a valid measure of the activity that survives confinement; the controlling dimensionless ratio is τR/τk.","The effective trap stiffness keff and center displacement rc give two clean experimental handles: keff=k/(1+Pe²/2) with rc=0 in regime I, and keff=k with rc=Vτk in regime II.","In regime II the power spectral density collapses onto the Lorentzian of the passive trapped particle at all propulsion speeds, providing a frequency-domain test that the active contribution has been neutralized."],"fun_headline_variants":["Boltzmann shape doesn't prove passivity: active signatures persist","Trap's shape misleads; timescale ratio tells if activity remains","Residual activity hidden in Gaussian-looking trap dynamics","Active motion leaves trace in dynamics despite Boltzmann profile"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The conclusion that regime II is devoid of residual activity rests on the Appendix B, Case-2 idealization that τR→∞ (sinϕ(t)=0), so the propulsion is exactly balanced by the restoring force; the paper's own steady-state picture includes slow orientational diffusion (τR=100 s) producing a slow azimuthal drift and a marginal widening of P(vθ_res) with propulsion speed, so the clean dichotomy softens if that drift counts as residual activity.","fun_headline_variants_meta":{"raw":{"variants":["Boltzmann shape doesn't prove passivity: active signatures persist","Trap's shape misleads; timescale ratio tells if activity remains","Residual activity hidden in Gaussian-looking trap dynamics","Active motion leaves trace in dynamics despite Boltzmann profile"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000703,"raw_usage":{"total_tokens":3016,"prompt_tokens":757,"completion_tokens":2259,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":2191}},"tokens_in":501,"tokens_out":2259,"duration_ms":18116,"temperature":1.0,"reasoning_tokens":2191,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T05:24:25.864502+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radial power spectrum of an HBABP with τR≫τk at several propulsion speeds, from V=0 up to V=10 μm/s; the paper predicts all curves collapse onto the passive Lorentzian. If the PSD changes with V, or if the measured effective stiffness from Boltzmann inversion changes with V, the claim that residual activity vanishes in this regime is refuted. The complementary check in regime I is that keff must follow k/(1+Pe²/2).","supporting_citations":[],"review_version":1}