{"id":"70cb6a19-46de-4bca-ba71-a94eeb07e98f","arxiv_id":"2608.07059","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"An active Brownian particle in a Maxwell-Voigt fluid has the same mean-square displacement as an active particle in a diffusing harmonic trap, demonstrated experimentally with a Janus colloid.","lead":"The authors show that an active Brownian particle in a single-relaxation viscoelastic fluid can be described as a self-propelled particle in a slowly diffusing harmonic well, and they emulate this medium with a computer-steered optical trap. Experiments on a Janus colloid reproduce the predicted mean-square displacement, whose shape is set by the ordering of three timescales.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Experimental validation is underdetermined: all four model parameters are fitted to the measured MSD, so the reported 'quantitative agreement' demonstrates fit flexibility rather than independent confirmation of the model.","rationale":"The reader's verdict of CONDITIONAL is appropriate. The analytical and simulation parts of the paper rest on well-established additive MSD results and a clearly specified propagator scheme; no internal mathematical inconsistency was identified in Eq. 2 or its derivation. The actual risk is entirely in the experimental claim of 'quantitative agreement,' which is weakened by fitting all four model parameters to the same MSD curve and by the direct contradiction between the Methods text (λ is fitted) and the Conclusions (λ is known from the trap). The paper also overstates that the environment's viscoelastic properties are 'unaffected by the strength of self-propulsion' while acknowledging that laser power couples trap stiffness to propulsion speed. These issues do not undermine the theoretical contribution, but they do mean the experimental validation, as reported, does not independently confirm the model. Hence the conditional verdict should stand, with the requirement that the authors either provide independent calibration of the parameters or soften the claim of experimental validation. My concern matches the reader's weakest assumption precisely; the strongest claim itself is not threatened by any soft spot I could identify in the theory or simulations.","tokens_in":10442,"tokens_out":2024,"duration_ms":19933,"concrete_test":"Perform independent calibration experiments on the same dynamic-trap setup: (1) With the active particle in a static trap (DHW=0), measure the stationary position distribution to extract k from equipartition (or Boltzmann inversion), and measure the mean-square displacement to extract V and τR from the known static-HW HBABP expression. (2) With a passive particle (V=0), measure the trap-center trajectory to determine DHW directly, and the passive MSD to determine τk and λ. (3) Then compute the predicted active MSD using Eq. 2 with these independently measured parameters and no free fitting, and compare to the experimental active MSD. If the predicted curve falls within the experimental uncertainty, the claim of quantitative agreement is validated; if systematic deviations appear, the model or the assumed parameter mapping is wrong.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that an ABP in a single-relaxation VE fluid is described by an HBABP with long-time diffusion—is well supported analytically and by simulation, where all parameters are set independently. The weak point is the experimental 'validation.' The Methods paragraph states: 'We compute the time-averaged MSD from the tracked trajectories and fit it to Eq. 2, obtaining τk, τR, λ, and V with fitting uncertainties below 10%.' Yet the Conclusions state 'with λ known from the trap rather than fitted.' These statements are mutually contradictory, and the paper never reports an independent measurement of any of the four fitted parameters. The trap stiffness k, the HW diffusivity DHW (hence λ = γHW/k), the propulsion speed V, and the orientational relaxation time τR are each controlled or measurable by separate experiments, but none is reported. Fitting four parameters to a single MSD curve will necessarily reproduce the observed shape, including the features shown in Fig. 4b,c. Consequently, the agreement between experiment and Eq. 2 is evidence of the fit's expressiveness, not of the model's predictive validity. The central theoretical construction may well be correct, but the experimental demonstration, as presented, cannot distinguish the proposed model from any four-parameter MSD ansatz that happens to share the same functional form.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that an active Brownian particle (ABP) in a single-relaxation Maxwell-Voigt fluid is equivalent to a harmonically bound ABP whose harmonic-well center undergoes free diffusion (HBABP with long-time diffusion). The authors derive the mean-square displacement as the sum of a static-well HBABP contribution and a free-diffusion term (Eq. 2, with the degenerate case Eq. 3), simulate the stochastic dynamics for five timescale orderings, and use a Pt-silica Janus colloid in a dynamic optical trap to test two of those orderings experimentally. They report quantitative agreement among theory, simulation, and experiment.","tokens_in":10718,"tokens_out":28520,"duration_ms":258080,"significance":"The theoretical framework is attractive: it reduces the three-timescale competition (equilibration time tau_k, persistence time tau_R, and viscoelastic relaxation time lambda) to a tractable stochastic model, and the dynamic-optical-trap realization offers a tunable experimental emulation whose viscoelastic response is not perturbed by activity. Strengths include the explicit analytic limits (V = 0 recovers HBBP with long-time diffusion; lambda -> infinity recovers static-well HBABP), the clearly specified simulation propagator, and the honest statement of limitations such as linear response and the coupling between laser power and propulsion speed. The main weakness is the experimental validation: all four parameters (tau_k, tau_R, lambda, V) are fitted to the same experimental MSD curve that Eq. 2 is supposed to predict, with no independent calibration of any of them. As a result, the reported agreement with experiment currently demonstrates the flexibility of a four-parameter fit rather than predictive validity of the model.","major_comments":[{"comment":"The Methods paragraph states: \"We compute the time-averaged MSD from the tracked trajectories and fit it to Eq. 2, obtaining tau_k, tau_R, lambda, and V with fitting uncertainties below 10%.\" The Conclusions, by contrast, claim \"Experiment tests the case lambda > tau_k, tau_R, with lambda known from the trap rather than fitted.\" These statements are mutually inconsistent. All four model parameters are obtained from the same experimental MSD curve shown in Fig. 4b,c, and the paper reports no independent measurement of the trap stiffness k (hence tau_k), the imposed HW diffusivity D_HW (hence lambda), the propulsion speed V, or the orientational relaxation time tau_R. With four free parameters, a fit of Eq. 2 to an MSD of this functional form will reproduce the observed shape essentially by construction. The quantitative agreement in Fig. 4 therefore does not validate the model unless independent calibrations are supplied. I ask the authors to add, for each of the two regimes, independently measured values of k from a passive trapped bead at the same laser power, D_HW from the prescribed trap trajectory, V from a static-trap or free-propulsion assay, and tau_R from orientation tracking, and then to refit the MSD with at most one free parameter, or to re-scope the experimental portion as a parameterization study rather than a validation.","section":"Experimental section and Conclusions"},{"comment":"There is an inconsistency between the stochastic equations and the simulation update. Eq. (1a) for x_HBABP, described as \"the ABP position relative to the HW,\" contains no coupling to the motion of the HW center. The simulation update, however, includes the term \"-Delta x_HW,i\" in the argument of the exponential, meaning the relative coordinate is driven by the trap-center displacement. For a particle in a moving harmonic well, the correct relative-coordinate equation is d x_HBABP/dt = -x_HBABP/tau_k - d x_HW/dt + V cos(phi) + noise, not Eq. (1a). These two systems have different MSDs: in the moving-well system, the contribution of the center diffusion is filtered through the trap relaxation and is not simply additive at all times. The authors should state which system Eq. (2) actually solves, add the missing -d x_HW/dt term to Eq. (1a), or explicitly identify the approximation (for instance lambda >> tau_k, so D_HW/D_s = tau_k/lambda is small) under which the additive form Eq. (2) is valid for the dynamic-trap experiment. As written, the simulation and the analytic model are not solving the same equations, and the numerical agreement with Eq. (2) is not a check of the model as stated.","section":"Model and simulation propagator (Eq. 1 and the paragraph after Eq. 4)"}],"minor_comments":[{"comment":"The first term in Eqs. (2) and (3) is typeset as 4D_HBABP tau (1 - e^{-tau/tau_k}). If this is literal, it gives a vanishing short-time slope instead of the required 4D_HBABP tau and an unbounded long-time growth instead of the plateau 4D_HBABP tau_k. The correct prefactor should be 4D_HBABP tau_k; please check the typesetting and ensure the subscript is not lost.","section":"Eqs. (2) and (3)"},{"comment":"The captions say \"solid lines denote fits from Eq. 2\" even though the simulation parameters were set independently. Please clarify whether any parameter is free in these fits, and if so, report the fitted values; otherwise, replace \"fits\" with \"Eq. 2\" plotted with the input parameters.","section":"Simulation results, Figs. 2 and 3"},{"comment":"The fitted values of tau_k, tau_R, lambda, and V are never listed, despite the statement that fitting uncertainties are below 10%. Reporting the actual values, with confidence intervals, is necessary for the reader to judge whether the fitted parameters are physically plausible and consistent with the trap and particle properties.","section":"Experimental results, Fig. 4"},{"comment":"The abstract calls tau_k and lambda \"the crossover and equilibration times of the VE fluid, tau_k and lambda, respectively.\" This is confusing because lambda is the relaxation time and tau_k is the equilibration time; consider using consistent terminology throughout.","section":"Abstract and nomenclature"},{"comment":"The data availability statement says the raw data are not public because of size. For a paper whose experimental claim rests on a four-parameter fit, it would be helpful to at least publish the fitted parameter values, the MSD data in tabular form, and the simulation code, so that the fit procedure and the agreement with Eq. 2 can be independently checked.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The theoretical and simulation parts of the manuscript are well within the scope of cond-mat.soft and are likely publishable after the experimental calibration issue is resolved. The contradiction between the Methods statement that all four parameters are fitted and the Conclusions claim that lambda is known from the trap is a concrete, fixable problem rather than a reason to reject. I would encourage the editor to request independent calibration measurements or a clear downgrade of the experimental claims. The coupling inconsistency between Eq. (1) and the simulation propagator also needs to be resolved, but it may be a presentation issue if the authors intended an approximation that is not stated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the theory section is solid, the simulations are clearly specified, and the experimental platform—a phoretically active Janus colloid in a dynamic optical trap—is genuinely new. But the experimental validation as written does not establish the model: every parameter in Eq. 2 is fitted to the same measured MSD, and the paper contradicts itself about whether λ is fitted or known from the trap.\n\nWhat is new is the extension of the diffusing-harmonic-well emulation to active Brownian particles, plus the regime classification by the ordering of τk, τR, and λ. The analytical MSD is a correct additive combination of the static-well HBABP result and free diffusion of the well center; the λ→∞ and V=0 limits reduce correctly to published results. The simulation scheme is reproducible from the text. All of that is in good shape, and I would not call the central model wrong.\n\nThe soft spot is the experimental section. The Methods say the trajectory MSD is fit to Eq. 2, obtaining τk, τR, λ, and V with uncertainties below 10%. The Conclusions say λ was known from the trap rather than fitted. Both cannot be true. No independent measurement of trap stiffness, propulsion speed, or imposed HW diffusivity is reported, so the agreement in Figs. 4b and 4c is fit quality, not prediction. That is a load-bearing weakness for the claim of quantitative experimental validation, though not for the theoretical construction. The paper should either report independent calibrations (Boltzmann-inversion stiffness, measured propulsion speed from isolated trajectories, imposed DHW) or explicitly reframe the experiment as a consistency check rather than validation.\n\nMinor point: because the paper's own limits recover previously published expressions, the genuinely new theoretical content is the three-timescale ordering and the experimental realization. That is worth stating plainly, but it is not a flaw by itself.\n\nThis paper is for people working on active matter in viscoelastic media or on optical-trap emulation of rheology. A serious referee could sort out the experimental issue. I would send it to review, with the expectation of a major revision on the experimental claims.","headline":"A clean analytic/simulation story about active Brownian motion in a single-relaxation viscoelastic fluid, with an experimental showcase whose validation is undercut by fitting all four parameters to the same MSD.","tokens_in":11244,"tokens_out":1974,"would_cite":true,"duration_ms":18901,"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":"An active Brownian particle in a single-relaxation viscoelastic fluid behaves exactly like a self-propelled colloid inside a harmonic well whose center itself undergoes ordinary Brownian motion, with the full mean-square displacement…","keywords":["Active Brownian particle","viscoelastic fluid","Maxwell-Voigt model","harmonically bound active Brownian particle","long-time diffusion","dynamic optical trap","Janus colloid","timescale ordering"],"falsifier":"Independently measure the trap stiffness $k$ (e.g., from the equipartition variance of a trapped passive bead), the imposed trap diffusivity $D_{\\mathrm{HW}}$ (from the variance of the generated Brownian trajectory), and the propulsion speed $V$ (from the short-time ballistic slope of the MSD in a static trap), then compare these values with those obtained by fitting Eq. 2 to the active-particle MSD. If the fitted and independently measured values disagree beyond experimental uncertainty in either of the two reported regimes, the claimed quantitative validation would fail.","tokens_in":10216,"feed_emoji":"🧪","tokens_out":6312,"duration_ms":49132,"temperature":0.7,"pith_summary":"The paper shows that an active Brownian particle (ABP) swimming in a single-relaxation viscoelastic fluid moves precisely like a self-propelled particle trapped in a harmonic well whose center diffuses. This equivalence reduces the complicated fluid memory to a tractable stochastic model whose mean-square displacement is a sum of three terms: confined thermal motion, active propulsion, and long-time diffusion of the well. The ordering of three timescales—equilibration in the well, orientation persistence, and medium relaxation—selects the observable dynamical regimes, from superdiffusive runs through elastic plateaus to a final diffusive regime. The authors validate the model with numerical simulations across all timescale orderings and with experiments using a Pt-coated Janus colloid in a dynamically steered optical trap that mimics the viscoelastic medium. The practical payoff is a tunable, activity-independent viscoelastic environment for studying active motion in structured fluids.","feed_headline":"Active particles in viscoelastic fluids act like swimmers in a drifting trap","feed_subtitle":"A single-relaxation fluid's memory maps onto a diffusing harmonic well, and three timescales decide how the swimmer moves.","key_machinery":"The central object is the Maxwell-Voigt (MV) model with a single spring stiffness $k$ shared by a Voigt element (spring in parallel with high-frequency dashpot $\\gamma_s$) and a Maxwell element (spring in series with low-frequency dashpot $\\gamma$). Its stochastic surrogate is a harmonic well of stiffness $k$ with solvent friction $\\gamma_s$ whose center diffuses with friction $\\gamma_{\\mathrm{HW}}$, so the particle position splits into $x(t)=x_{\\mathrm{HBABP}}(t)+x_{\\mathrm{HW}}(t)$. The MSD formula (Eq. 2) adds the confined active contribution, the active propulsive contribution, and free diffusion of the well, with the three timescales $\\tau_k=\\gamma_s/k$, $\\tau_R=1/D_R$, and $\\lambda=\\gamma_{\\mathrm{HW}}/k$ controlling each regime.","core_discovery":"The paper contends that an active Brownian particle in a single-relaxation Maxwell-Voigt fluid is dynamically equivalent to a harmonically bound active Brownian particle whose confining well itself undergoes free diffusion. The particle coordinate separates into an HBABP part (position relative to the diffusing well) and the well-center displacement, so the total mean-square displacement is the sum of the static-well HBABP result and $4D_{\\mathrm{HW}}\\tau$, as given by Eq. 2. Because both the well equilibration time $\\tau_k$ and the relaxation time $\\lambda$ stem from the same spring stiffness $k$ but different friction coefficients, the relative ordering of $\\tau_k$, $\\tau_R$, and $\\lambda$ fully determines the shape of the MSD. Simulations confirm the analytical MSD in all orderings, and experiments with a Pt-coated Janus colloid in a dynamic optical trap reproduce the persistence-dominated and confinement-dominated regimes quantitatively.","pith_inferences":["If the equivalence is quantitative, the same three-timescale competition should appear in real single-relaxation VE fluids such as wormlike micelles, so the HBABP-with-long-time-diffusion formula could be used to extract an active particle's persistence time from a single MSD measurement.","The model's linear single-relaxation restriction suggests a natural multi-mode extension: a superposition of independent diffusing wells would produce a generalized Maxwell fluid with multiple plateaus and memory-dependent diffusion.","Because the optical trap couples laser power to both stiffness and thermophoretic speed, the claim that propulsion speed is independent of the medium's rheology could be tested more cleanly with chemically powered Janus colloids, where the propulsion mechanism is decoupled from the trap parameters."],"forward_implications":["If the equivalence holds, the full mean-square displacement of an ABP in a single-relaxation VE fluid is known in closed form from four parameters, making active-probe microrheology a straightforward three-timescale analysis.","The ordering of $\\tau_k$, $\\tau_R$, and $\\lambda$ dictates which dynamical regimes appear: persistence-dominated dynamics ($\\tau_R < \\tau_k$) shows superdiffusion followed by activity-enhanced diffusion and a plateau, while confinement-dominated dynamics ($\\tau_k < \\tau_R$) shows two plateaus; if the medium relaxes before reorientation ($\\lambda < \\tau_R$), activity signatures vanish entirely.","The dynamic optical trap provides a configurable VE medium whose rheological parameters ($\\tau_k$ and $\\lambda$) are set by trap stiffness and imposed trap diffusivity and are unaffected by the propulsion speed, enabling experiments in regimes that are inaccessible with real polymeric fluids.","In the limit $\\lambda \\to \\infty$ the model reduces to the static-well HBABP, and in the passive limit $V=0$ it recovers the known HBBP with long-time diffusion, so the new result nests both existing cases as special limits."],"supporting_citations":[{"why":"Supplies the passive HBBP with long-time diffusion model that the active extension builds on.","marker":"[34]"},{"why":"Provides the dynamic optical trap method used to emulate the diffusing harmonic well in experiment.","marker":"[35]"},{"why":"Gives the static-well HBABP MSD and the two-timescale interplay that this work extends to a moving well.","marker":"[36]"},{"why":"Establishes the stable optical trapping of thermophoretically active Janus colloids used in the experiments.","marker":"[38]"},{"why":"Demonstrates self-thermophoretic propulsion of Janus particles, the active mechanism for the particle.","marker":"[42]"}],"fun_headline_variants":["Diffusing harmonic well captures viscoelastic dynamics of active particles","Active swimmers in a viscoelastic fluid: equivalent to a diffusing trap","One diffusing well explains active particle motion in viscoelastic media","Viscoelastic fluids act as drifting traps for active Brownian particles","Three timescales govern active motion in viscoelastic fluids via a diffusing trap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the four fitted parameters ($\\tau_k$, $\\tau_R$, $\\lambda$, and $V$) extracted from the measured mean-square displacement correspond to the actual physical stiffness, diffusivity, and propulsion speed of the trap and particle, rather than being merely free fitting parameters.","fun_headline_variants_meta":{"raw":{"variants":["Diffusing harmonic well captures viscoelastic dynamics of active particles","Active swimmers in a viscoelastic fluid: equivalent to a diffusing trap","One diffusing well explains active particle motion in viscoelastic media","Viscoelastic fluids act as drifting traps for active Brownian particles","Three timescales govern active motion in viscoelastic fluids via a diffusing trap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000968,"raw_usage":{"total_tokens":4120,"prompt_tokens":952,"completion_tokens":3168,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":3072}},"tokens_in":568,"tokens_out":3168,"duration_ms":20036,"temperature":1.0,"reasoning_tokens":3072,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:34:34.532292+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Independently measure the trap stiffness $k$ (e.g., from the equipartition variance of a trapped passive bead), the imposed trap diffusivity $D_{\\mathrm{HW}}$ (from the variance of the generated Brownian trajectory), and the propulsion speed $V$ (from the short-time ballistic slope of the MSD in a static trap), then compare these values with those obtained by fitting Eq. 2 to the active-particle MSD. If the fitted and independently measured values disagree beyond experimental uncertainty in either of the two reported regimes, the claimed quantitative validation would fail.","supporting_citations":[{"cited_title":"Khan and T","cited_arxiv_id":null,"evidence_quote":"Supplies the passive HBBP with long-time diffusion model that the active extension builds on."},{"cited_title":"Halder and M","cited_arxiv_id":null,"evidence_quote":"Provides the dynamic optical trap method used to emulate the diffusing harmonic well in experiment."},{"cited_title":"Halder and M","cited_arxiv_id":null,"evidence_quote":"Establishes the stable optical trapping of thermophoretically active Janus colloids used in the experiments."},{"cited_title":"Jiang, N","cited_arxiv_id":null,"evidence_quote":"Demonstrates self-thermophoretic propulsion of Janus particles, the active mechanism for the particle."}],"review_version":1}