{"id":"7e8a6d19-0ddc-41bf-b9e6-38d841dc5595","arxiv_id":"1908.08157","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An active Brownian particle in an optical trap develops a double-peaked position histogram and a low-frequency rise in its power spectrum, describable by convolving thermal and active contributions.","lead":"This paper traps a single light-driven Janus particle in an optical trap and shows its position distribution changes from a bell curve to two peaks when the particle swims. It also connects the spatial distribution to the temporal power spectrum and warns that simple 'effective temperature' measures can mislead in active systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 10 sets the low-frequency PSD corner at 2πD_r, but the ABP model implies D_r; the claimed quantitative PSD–rotational-diffusion connection is unsupported as written.","rationale":"I read the paper in good faith. The experimental observations—bimodal position histograms, low-frequency PSD rise, and the HPP/PSD energy integral—are plausible, and the simulation is genuinely parameter-free for the Langevin model. The central quantitative claim, however, is the PSD-ratio formula Eq. 10, and the paper itself labels it as a fitted ansatz ('We can then write in the following format') without derivation. Deriving the same quantity from Eq. 2 shows that the active contribution to the PSD ratio is a Lorentzian in ω with corner at D_r, not at 2πD_r. This is not a disagreement with a consensus model; it is an internal inconsistency in the paper's own mathematical setup. The reader's weakest-assumption analysis focused on drift contamination of the low-frequency PSD plateau, which is a legitimate experimental risk that the authors themselves flag. My concern is more fundamental: even if the experimental plateau is clean, Eq. 10 as written does not correctly represent the predicted spectrum. That said, the issue is correctable and does not invalidate the qualitative conclusions or the histogram analysis, so I do not move the verdict beyond conditional. The paper should either derive Eq. 10 explicitly, correct the frequency convention, or reframe the connection as qualitative. A direct numerical test of the corner frequency will resolve the matter quickly.","tokens_in":11790,"tokens_out":8271,"duration_ms":88138,"concrete_test":"Derive S_ratio(ω) from Eq. 2 as outlined and compare with Eq. 10. Concretely, run the paper's parameter-free simulation with D_r = 0.05/s, v0 = 3.5 μm/s, D_t = 0.053 μm²/s, and k = 1.5 pN/μm, then compute the PSD ratio from the simulated trajectory. If the angular frequency at which the active plateau falls to half its zero-frequency value is ω ≈ 0.05 rad/s (f ≈ 0.008 Hz), Eq. 10's corner at 2πD_r ≈ 0.314 rad/s (f ≈ 0.05 Hz) is wrong; if it is instead ω ≈ 0.314 rad/s, then the derivation from Eq. 2 must be re-examined. This single check settles whether Eq. 10 can support the claimed quantitative link between the low-frequency PSD rise and rotational diffusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest claim is the quantitative PSD-ratio formula, Eq. 10, said to connect the low-frequency rise to v0, τ_r, and D_t. The paper offers no derivation, and a direct derivation from its own Langevin model, Eq. 2, contradicts the stated corner frequency. For the passive case, S_passive(ω) = 2ξ k_BT / (k^2 + ξ^2ω^2). The active orientational driving has autocorrelation ⟨cosθ(t)cosθ(0)⟩ = (1/2)e^(−D_r|t|), so its PSD is ξ^2 v0^2 D_r / (D_r^2 + ω^2). Since the same harmonic trap filters active and passive position signals, the ratio is S_active/S_passive = 1 + [v0^2/(2D_t)] [D_r^2/(D_r^2 + ω^2)]? Actually the active ratio term is 1 + [v0^2D_r/(2D_t(ω^2 + D_r^2))], which at ω=0 gives the stated 1 + v0^2τ_r/(2D_t). The spectral shape is therefore 1/[1 + (ω/D_r)^2], with corner at ω = D_r (equivalently f = D_r/2π if ω is angular). Equation 10 instead writes (ω/ω_r)^2 with ω_r = 2πD_r, placing the corner a factor 2π higher. Since Eq. 8 defines ω through e^{iω(t−t′)} and is therefore an angular frequency, this is an internal inconsistency, not a matter of convention. If ω in Eq. 10 is instead meant to be a cyclic frequency, then Eq. 10 is inconsistent with the Fourier convention used in Eq. 8. Either way, the quantitative claim that the PSD reveals the rotational diffusion time scale τ_r through Eq. 10 is unsupported. This is independent of, and more decisive than, the drift concern: even with perfect data, the Lorentzian corner is mis-specified. The drift limitation noted in Section 3.3 further weakens any experimental determination of the plateau, but the analytic error is the more load-bearing problem for the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports experiments and numerical simulations of a single induced-charge electrophoretic Janus particle confined in an optical trap, used as a model active Brownian particle. The authors measure position histograms and power spectral densities at controlled active speeds and trap stiffnesses, observing a Gaussian-to-bimodal transition in the histogram and a low-frequency rise in the PSD. They propose a convolution decomposition of thermal and active contributions to the histogram, introduce characteristic length scales and a critical active speed, and claim that the PSD ratio between active and Brownian particles has the Lorentzian form of Eq. (10), with a low-frequency plateau related to the free-space swim speed and a corner frequency set by the rotational diffusion time. They also compare effective temperatures defined from spatial and temporal fluctuations and discuss limitations of classical thermodynamic interpretations in active systems.","tokens_in":12250,"tokens_out":5884,"duration_ms":58789,"significance":"If the claims hold, the paper provides a well-controlled experimental model system with a large dynamic range, and its convolution method for separating thermal and active contributions is a clean consequence of the linearity of the overdamped Langevin equation. The simulations are presented as parameter-free, using independently measured v0, D_t, D_r, and trap stiffness, and the reported Gaussian-to-bimodal HPP transition is consistent with the ABP physics in confinement. The proposed PSD ratio formula, if correctly specified, would be a useful quantitative tool connecting confined temporal fluctuations to the free-space active parameters. These strengths make the core experimental platform valuable. However, the central quantitative PSD claim in Eq. (10) is currently undermined by an internal frequency-convention inconsistency, and the paper itself admits a drift sensitivity that directly affects the low-frequency data used to test that formula.","major_comments":[{"comment":"The corner frequency in Eq. (10) is internally inconsistent with the manuscript's own Langevin model. Equation (8) defines omega as an angular frequency through the exponential e^{i omega (t-t')}. For the overdamped equation Eq. (2), the orientational noise has autocorrelation <cos(theta(t)) cos(theta(0))> = (1/2) e^{-D_r |t|}, whose PSD is proportional to D_r/(D_r^2 + omega^2). Since the same harmonic trap filters both active and passive position signals, the ratio of active to Brownian PSD should be S_ratio(omega) = 1 + [v0^2/(2 D_t)] [D_r^2/(D_r^2 + omega^2)], with the corner at angular frequency omega = D_r. Equation (10) instead places the corner at omega_r = 2 pi D_r, a factor 2 pi higher. If omega in Eq. (10) were intended to be a cyclic frequency, that would contradict the convention used in Eq. (8). Either way, the quantitative statement that the PSD reveals the rotational diffusion time scale through Eq. (10) is not supported as written. This is a load-bearing issue because the 'characteristic frequency' claim and the extraction of tau_r from the PSD depend on this formula.","section":"§3.3, Eq. (10)"},{"comment":"The manuscript states in Section 3.3 that 'PSD measurement is highly sensitive to drift noise and at the low frequency (long time scale), the slow drift bias may affect the experimental precision [reference needed].' The low-frequency plateau of the PSD is exactly the region used to test Eq. (10) and to define the zero-frequency effective temperature in Fig. 6. Because drift contamination can raise the low-frequency spectral power, the observed agreement between the experimental PSD ratio and Eq. (10) does not by itself establish the quantitative connection to v0, tau_r, and D_t. The authors need either to demonstrate that drift is negligible in their measured low-frequency data, to correct for drift, or to explicitly restrict the quantitative claim to drift-free simulations.","section":"§3.3, Fig. 5"},{"comment":"Equation (10) is introduced as an ansatz ('We can then write') with no derivation from Eq. (2). The manuscript cites Szamel's work for a similar expression, but that work is for the active component only and does not include the thermal baseline. A direct derivation from the linear Langevin equation is straightforward, as outlined in the first comment, and it yields a different corner frequency. The authors should provide the derivation or at least correct the frequency conventions and verify that the simulation data actually match the corrected Lorentzian form. Without this, the claimed quantitative agreement between Eq. (10) and the simulation data is not credible.","section":"§3.3, Eq. (10) and Fig. 5D"}],"minor_comments":[{"comment":"There are numerous typographical and grammatical errors that should be fixed: 'Langevine' for Langevin, 'Jans' for Janus, 'serval' for several, 'microecology' likely for micro-rheology, 'the MSD vs time graph appear s', and 'showing in Fig.1 B)' should be 'as shown in Fig. 1B'. These do not affect the science but reduce readability.","section":"Throughout"},{"comment":"Two places in Section 3.3 contain the placeholder '[reference needed]' for statements about drift sensitivity and integration time. These should be supplied or the statements should be removed.","section":"§3.3"},{"comment":"Equation (9) appears to be dimensionally inconsistent for a position PSD. The numerator 12 pi eta R k_B times T and the denominator ((k_OT/R)^2 + (eta omega)^2) do not yield units of position power spectral density (m^2 s for angular frequency). This may be a typesetting error (perhaps a missing factor or an incorrect grouping), but it should be corrected because Eq. (9) is presented as the theoretical Brownian PSD reference.","section":"§3.3, Eq. (9)"},{"comment":"Equation (2) contains 'v0 D' in the active term, which is likely a typo for 'v0 theta_hat' or 'F0 theta_hat' given the form of Eq. (1). In addition, the description of the MSD slopes in Section 3.1 is confusing: the text says the short-time slope is 2 while also saying thermal diffusion gives a slope of 1 at very short times. This should be clarified.","section":"§2.4 and §3.1"},{"comment":"The text references 'Fig 3 C and D' and discusses simulated HPPs, but Fig. 3 appears to contain only panels A and B, and panel C's caption contains a placeholder 'active speed of xxx'. The figure must be completed or the text adjusted accordingly.","section":"§3.2, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of cond-mat.soft and the experimental system is promising, but the central quantitative PSD claim needs a substantive correction. The frequency-convention error in Eq. (10) is not a cosmetic issue; it directly affects the claimed extraction of the rotational diffusion time from the PSD. I would encourage the editor to request a revision that derives or corrects Eq. (10), addresses the drift sensitivity explicitly, and cleans up the numerous placeholders and typos. The novelty claim about being the first to observe boundary accumulation in an optical trap should also be toned down unless previous literature has been searched more carefully."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is worth a look for the experiment, but the central quantitative PSD claim is wrong. The ICEP Janus particle in a quadratic optical trap is a genuinely nice system: large dynamic range, clean optical forces, and a parameter-free Langevin simulation that reproduces the bimodal histograms. The convolution decomposition in Eqs 5–7 is exact for the linear overdamped equation, and the demonstration that the full HPP is the convolution of thermal and active contributions is a solid, well-illustrated point. I give them full credit for that portion.\n\nThe problem is Eq 10. From their own Eq 2, the active orientational noise has autocorrelation ~ e^{-D_r t}, so in angular frequency the active contribution to the PSD ratio is Lorentzian with corner at ω = D_r, not at 2πD_r. The paper's ω_r = 2πD_r places the corner a factor 2π too high. This is not a convention issue: Eq 8 uses an angular-frequency Fourier transform, so Eq 10 is internally inconsistent. Even with perfect data and unaffected by drift, the claimed quantitative connection between the PSD corner and the rotational diffusion time would not hold. That is a load-bearing error for the paper's 'two characteristic frequencies' narrative.\n\nOther soft spots are secondary. The low-frequency PSD is drift-prone, as the authors themselves note (with two dangling '[reference needed]' placeholders, which should have been filled before submission). The claim to be the first to see boundary accumulation in an optical trap is overstated, given prior bacterial and acoustic-trap work. The dissipation-fraction interpretation in Fig 6B is muddled. And the figures lack error bars throughout.\n\nWhat is actually new: the experimental combination of ICEP with simultaneous optical trapping and QPD tracking, and the explicit convolution-based decomposition of the HPP into thermal and active components. That part is publishable. The PSD analysis, as written, is not.\n\nFor a peer-review decision: I would send it out despite the flaw, because the system and the HPP convolution are worth refereeing and the PSD error is correctable. But I would tell the authors that Eq 10 must be fixed (ω_r → D_r or an appropriate redefinition and re-fit), the drift issue must be addressed quantitatively, and the novelty claim toned down. If the experiment genuinely shows a corner at the corrected location, the paper becomes a solid contribution. As is, I would not cite the PSD formula, though I might cite the system.\n\nBring to reading group? Maybe, as a cautionary tale in internal consistency checks.\n\nBest,\n[Your name]","headline":"A well-constructed experimental ABP-in-trap study with a correct HPP convolution analysis, but Eq 10's factor-2π corner error invalidates the paper's central quantitative PSD claim.","tokens_in":12803,"tokens_out":6993,"would_cite":true,"duration_ms":65175,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The power spectrum of a trapped active particle is a Brownian spectrum plus a Lorentzian bump whose height encodes swim speed and rotational diffusion.","keywords":["active Brownian particle","optical trap","power spectral density","position histogram","induced-charge electrophoresis","Janus particle","effective temperature","rotational diffusion"],"falsifier":"Perform the same PSD measurement at zero applied voltage while imposing a slow stage drift; if the low-frequency plateau rises above the Brownian baseline without any active swim speed, then drift alone can mimic the Eq. 10 signature. A cleaner quantitative test is to compare the swim speed fitted from the plateau height with the independently measured free-space MSD swim speed: systematic disagreement at low voltage would show the plateau is contaminated.","tokens_in":11579,"feed_emoji":"🔬","tokens_out":10428,"duration_ms":93991,"temperature":0.7,"pith_summary":"Using a single metallic Janus particle driven by induced-charge electrophoresis and held in an optical trap, the paper claims that a trapped active Brownian particle's fluctuations factor into a thermal part and an active part, in both space and time. Spatially, the position histogram is the convolution of a Gaussian thermal distribution with an active-only distribution, which is why the histogram splits from one peak to two as the swim speed grows. Temporally, the power spectral density is the Brownian trap spectrum times a one-plus-Lorentzian ratio whose plateau height is set by $v_0^2\\tau_r/(2D_t)$ and whose corner frequency is set by the rotational diffusion rate. If the claim is right, a single trapped-particle measurement can yield the swim speed and rotational relaxation directly, and the integrated-spectrum effective temperature should not be treated as a universal measure of dissipation because it depends on trap stiffness.","feed_headline":"Noise spectrum reveals the swim speed of a trapped active particle","feed_subtitle":"Plateau height encodes swim speed and rotational diffusion; one measurement gives both.","key_machinery":"The load-bearing object is the normalized power-spectral-density ratio of Eq. 10, a one-plus-Lorentzian form that attaches the active swim speed to the low-frequency plateau and the rotational diffusion coefficient to the corner frequency. Around it, the paper builds a three-length-scale criterion: the force-balance length $l_1 = v_0\\eta/k_{\\mathrm{OT}}$, the thermal width $l_2 = \\sqrt{k_BT/k_{\\mathrm{OT}}}$, and the persistence length $l_3 = v_0/D_r$; their ratios decide whether the position histogram stays Gaussian, flattens, or becomes bimodal. The convolution identity $P(r) = P_A * P_T(r)$ is what lets the paper decouple active and thermal contributions in space, while the PSD ratio decouples them in time.","core_discovery":"The paper's central result is that the steady-state fluctuations of an active Brownian particle in a quadratic trap separate into a passive part and an active part. In the frequency domain, the normalized power spectrum is\n\n$$\\frac{S_{\\mathrm{active}}(\\omega)}{S_{\\mathrm{Brownian}}(\\omega)} = 1 + \\frac{$v_0^{2}$\\tau_r/(2D_t)}{(\\omega/\\omega_r)^2 + 1}, \\qquad \\omega_r = 2\\pi D_r,$$\n\nso the active contribution is a Lorentzian added to the Brownian baseline, with the low-frequency plateau height controlled by the free-space swim speed $v_0$ and rotational relaxation time $\\tau_r = 1/D_r$. In real space, the full position distribution is the convolution of the thermal Gaussian and the distribution of a purely active particle in the same trap; the crossover from a single peak to a bimodal distribution is governed by three length scales ($v_0\\eta/k_{\\mathrm{OT}}$, $\\sqrt{k_BT/k_{\\mathrm{OT}}}$, $v_0/D_r$) and a critical swim speed. The paper also argues that the integrated PSD and the integrated histogram give the same mean potential energy, but that this integrated quantity is trap-stiffness dependent and therefore should not be interpreted as a universal dissipation or effective temperature.","pith_inferences":["If Eq. 10 survives drift-control checks, the same trapped-particle setup could serve as a self-calibrating probe: two PSD ratio values at known frequencies would determine both $v_0$ and $D_r$ in a single run.","The one-plus-Lorentzian structure may persist for weakly anharmonic traps, with the effective trap curvature sampled by the particle replacing $k_{\\mathrm{OT}}$; checking this would separate geometry effects from genuine active dynamics.","The convolution decomposition suggests a practical assay for weak activity in biological samples: calibrate the passive trap response, then infer the active-only distribution from the difference between measured and convolved histograms."],"forward_implications":["The low-frequency plateau of the PSD ratio directly measures $v_0^2\\tau_r/(2D_t)$, so the swim speed can be read from the spectrum without tracking the particle's free-space trajectory.","Because the Brownian baseline is identical at high frequency, the active contribution can be isolated by subtracting the passive spectrum from the active spectrum at each frequency.","The shape of the position histogram is governed by the two ratios $l_1/l_2$ and $l_1/l_3$, so increasing trap stiffness at fixed swim speed pushes a bimodal histogram back toward Gaussian.","The integrated PSD and the integrated histogram yield the same mean potential energy, giving a common measure that links spatial and temporal fluctuation data.","The trap-stiffness dependence of the integrated PSD means it cannot be used as a universal dissipation or effective temperature for active matter without specifying the confinement."],"supporting_citations":[{"why":"derives the Brownian harmonic-trap power spectrum that serves as the baseline in Eq. 11","marker":"[34]"},{"why":"gives the active-only PSD form that Eq. 10 extends by adding the thermal baseline","marker":"[35]"},{"why":"establishes induced-charge electrophoresis as the propulsion mechanism that sets the swim speed","marker":"[22]"},{"why":"provides the overdamped Langevin simulation approach used for the numerical histograms and PSDs","marker":"[28]"},{"why":"provides the acoustic-trap active-particle experiment whose boundary accumulation is compared with the optical-trap histograms","marker":"[17]"},{"why":"provides the thermal-fluctuation-dominated regime of a bacteria-driven active particle in a small trap, used to place the new data in parameter space","marker":"[18]"}],"fun_headline_variants":["Trap noise splits into passive and active parts","Active particle in trap: bimodal position, Lorentzian noise","Spatial and temporal fluctuations connect in a trap","Swim speed and rotation from a single noise spectrum","Active Brownian particle's trap fluctuations decoded"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim rests on the low-frequency plateau being real active-particle signal; if slow instrumental drift seeps into that plateau, the fitted swim speed and the effective-temperature comparison lose quantitative support.","fun_headline_variants_meta":{"raw":{"variants":["Trap noise splits into passive and active parts","Active particle in trap: bimodal position, Lorentzian noise","Spatial and temporal fluctuations connect in a trap","Swim speed and rotation from a single noise spectrum","Active Brownian particle's trap fluctuations decoded"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1791,"prompt_tokens":1024,"completion_tokens":767,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":692}},"tokens_in":640,"tokens_out":767,"duration_ms":7619,"temperature":1.0,"reasoning_tokens":692,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:47:35.210471+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same PSD measurement at zero applied voltage while imposing a slow stage drift; if the low-frequency plateau rises above the Brownian baseline without any active swim speed, then drift alone can mimic the Eq. 10 signature. A cleaner quantitative test is to compare the swim speed fitted from the plateau height with the independently measured free-space MSD swim speed: systematic disagreement at low voltage would show the plateau is contaminated.","supporting_citations":[{"cited_title":"Berg-Sø rensen and H","cited_arxiv_id":null,"evidence_quote":"derives the Brownian harmonic-trap power spectrum that serves as the baseline in Eq. 11"},{"cited_title":"Szamel, Phys","cited_arxiv_id":null,"evidence_quote":"gives the active-only PSD form that Eq. 10 extends by adding the thermal baseline"},{"cited_title":"Gangwal, O","cited_arxiv_id":null,"evidence_quote":"establishes induced-charge electrophoresis as the propulsion mechanism that sets the swim speed"},{"cited_title":"Volpe, S","cited_arxiv_id":null,"evidence_quote":"provides the overdamped Langevin simulation approach used for the numerical histograms and PSDs"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the acoustic-trap active-particle experiment whose boundary accumulation is compared with the optical-trap histograms"},{"cited_title":"Argun, A","cited_arxiv_id":null,"evidence_quote":"provides the thermal-fluctuation-dominated regime of a bacteria-driven active particle in a small trap, used to place the new data in parameter space"}],"review_version":1}