{"id":"3d534e57-fb9f-47c3-aa77-c154ce83a17c","arxiv_id":"2508.18789","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Simulations show that the location of the counterforce on the fluid, not just the active monomer, determines whether an active polymer creates pusher or puller flow fields, and that stiff chains transmit forces in a way that can flip the expected flow type.","lead":"Using computer simulations, the authors show that putting the balancing fluid force in different places around an active polymer chain changes the flow pattern it creates, making it act like a puller or a pusher swimmer. The result gives a simple design rule for controlling the motion of tiny synthetic or biological swimmers.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pusher/puller inversion for d=±σ is based on near-field flows; far-field dipole is only fitted for d=±4.5σ, so the inversion may be an artifact of the finite counterforce volume, compressibility, and torque noise.","rationale":"The reader's weakest_assumption identifies the counterforce construction as potentially unfaithful to a physical swimmer, which is a real concern. My analysis sharpens this by pointing to a specific evidential gap: the central claim is supported by near-field flow visualizations for d=±σ, but the only quantitative dipole extraction is performed for d=±4.5σ, where the counterforce is far from the polymer. The finite size of the counterforce volume, the compressibility, and the lack of strict torque-freeness could all distort the near-field flow in a way that mimics pusher/puller behavior without a genuine asymptotic dipole. This concern does not refute the paper's results but does support the CONDITIONAL verdict: the method is promising and internally consistent, yet the key classification lacks the far-field or benchmark validation needed to rule out numerical artifacts. I therefore recommend keeping the verdict unchanged. My agreement with the reader is partial because I emphasize a different facet: not just the faithfulness of the counterforce model, but the missing asymptotic check for the main configurations.","tokens_in":18225,"tokens_out":5669,"duration_ms":59184,"concrete_test":"For the stiff (βκ=100) hd- and td+ chains with d=±σ, simulate in a box of side L=40σ and extract the far-field velocity along the perpendicular direction at the x3 plane of maximum flow, fitting v_sol·e_ϱ = -p/(8πη)ϱ^{-2} (Eq. 27) over ϱ ∈ [5σ, 15σ]. If the sign of p disagrees with the near-field classification in Fig. 6 (i.e., p>0 for hd- or p<0 for td+), the reported inversion is an artifact of the local counterforce implementation rather than a consequence of internal force propagation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that counterforce placement tunes the flow field rests on the near-field velocity patterns in Fig. 6 for the d=±σ models. However, the only quantitative validation of a force-dipole far field is performed in Sec. VIII for d=±4.5σ (Figs. 10-11), where the counterforce center is far from the active monomer. For d=±σ, the counterforce is distributed over a sphere of diameter 2σ centered at r_cf (Eq. 9), overlapping the active monomer and its neighbor. The reported 'puller' flow for stiff hd- and 'pusher' flow for stiff td+ are read off at finite distances where the local counterforce volume, the ~10% density inhomogeneity (Sec. V), and the admitted lack of strict torque-freeness can each contribute flow components of comparable or dominating magnitude. In a 3D Stokes flow, a spurious rotlet decays as r^-2, the same as a force dipole, so torque fluctuations could contaminate both the sign and magnitude of the measured dipole. Because no benchmark against a known swimmer is provided, and the far-field asymptotics for d=±σ are not shown, the claimed 'tuning' of pusher/puller flows by counterforce placement is not yet established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports multi-particle collision dynamics (MPCD) simulations of active polymer chains (N=10) in which a single end monomer – head or tail – experiences a constant force along the local bond direction. To keep the system force-free, the authors introduce a local counterforce distributed over solvent particles inside a sphere of diameter 2σ centered at r_cf = R_iact + d e_f (Eq. 9). Four models are defined by active site (head/tail) and counterforce displacement (d=±σ): hd+, hd−, td+, td−. The paper examines swimming velocity, backbone orientational autocorrelation, MSD, radius of gyration, bond-angle distributions, monomer-resolved force projections, and solvent flow fields. For stiff chains, the active site has little effect on structural/dynamic properties, while for flexible chains head activity induces effective stiffening and tail activity causes crumpling and a 'cat's tail' sub-diffusive regime. The main claim is that the counterforce position tunes the hydrodynamic flow fields: hd+ and td− give the naively expected puller and pusher flows, whereas stiff hd− and td+ show the opposite, an inversion attributed to force redistribution along the polymer backbone. Far-field fits for d=±4.5σ confirm dipolar flows with opposite dipole signs.","tokens_in":18641,"tokens_out":7749,"duration_ms":75527,"significance":"The central idea – that the flow topology of an active polymer can be selected by counterforce placement, and that internal force transmission can override the naive local force-dipole expectation – is interesting and would be of value to the active-matter and soft-matter communities if firmly established. The simulation study is systematic and internally consistent across several observables (swimming velocity, RG, Cf, MSD, flow fields, per-monomer forces). The monomer-resolved force analysis in Fig. 9 directly supports the interpretation that the strongest backbone force lies near the counterforce volume. However, the load-bearing evidence for the 'tuning' claim is Fig. 6, which is qualitative, whereas the only quantitative far-field dipole analysis (Sec. VIII, Figs. 10–11) is performed for d=±4.5σ. The potential artifacts from torque fluctuations, compressibility, and the finite counterforce volume are acknowledged but not quantified. The new counterforce method is also not validated against a known test case. Thus the paper is a promising contribution that requires additional quantitative support before the main claim can be accepted.","major_comments":[{"comment":"The central claim that the flow field can be tuned by counterforce placement rests on the visual classification in Fig. 6 for d=±σ models. The only quantitative evidence for a force-dipole far field is the fit to Eq. (27) for d=±4.5σ (Figs. 10–11). For d=±σ the counterforce volume (diameter 2σ centered at r_cf = R_iact ± σ e_f) overlaps the active monomer and its neighbor, so the flow at the distances shown is not necessarily in the dipole regime. The ~10% density inhomogeneity (Sec. V) and the fluctuating torque can contribute velocity components of comparable magnitude. Please provide a quantitative far-field analysis (e.g., r^{-2} fit vs. r^{-3}, or inclusion of a rotlet term) for the d=±σ cases, or otherwise demonstrate that the classification is robust to these local effects.","section":"Sec. V and Sec. VIII (Fig. 6, Figs. 10–11)"},{"comment":"The system is not torque-free: distributing -F_a over n_cf solvent particles inside a sphere produces a stochastic torque whose magnitude fluctuates as n_cf and the particle positions fluctuate. The authors note that the torque fluctuates around zero, but they do not quantify it. In a 3D Stokes flow a rotlet decays as r^{-2}, the same as a force dipole, so a spurious torque could contaminate both the sign and the magnitude of the measured dipole in Fig. 6 and the fits of Sec. VIII. Please report the distribution and RMS of the total torque on the solvent (or the polymer+solvent system) and estimate the resulting rotlet velocity field. A simple control would be to compare the flow of a symmetric force pair (e.g., two equal and opposite forces far from the chain) with and without the stochastic torque.","section":"Sec. II.B (Eqs. 9–10) and Sec. V"},{"comment":"The new locally-tuned counterforce scheme is not validated against a known hydrodynamic test case. For example, the flow field of a single active monomer with a distant counterforce (or a rigid dumbbell) could be compared with the analytical Oseen tensor/dipole solution, and the effect of the finite counterforce sphere (diameter 2σ) and the MPCD collision-cell size could be quantified. Without such a benchmark it is difficult to rule out artifacts specific to the method, particularly for the quantitative dipole strengths p reported in Sec. VIII (p = 15.7 and −17.3 mσ²/τ²).","section":"Sec. II.B (method validation)"}],"minor_comments":[{"comment":"The phrase 'theoretical MSD' / 'prediction' is too strong: v0, τr, and Dpass are all measured from the same simulations, so Eq. (18) is a self-consistency check rather than an independent prediction. Please rephrase.","section":"Sec. III, Eq. (18)"},{"comment":"Typo: 'beding rigidity' should be 'bending rigidity'.","section":"Sec. IV"},{"comment":"The symbol lcm is used in the Reynolds number formula but the text says choosing lc = Nσ as characteristic size; please align the notation.","section":"Sec. III, Eq. (15)"},{"comment":"The force dipole strengths p are quoted without uncertainties. Please provide standard errors or confidence intervals from the fits.","section":"Sec. VIII"},{"comment":"The caption is uninformative; please specify which panels correspond to which of the four models and the stiffness values, and define the arrows/color scale in the caption.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the simulation effort is solid. The main concern is that the central claim is currently supported by a qualitative figure; the missing quantitative analysis can, in my view, be supplied with modest additional effort. I also note that no reproducibility code is provided, and the data availability statement is unusually brief; this is not a blocker but would strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper introduces a genuinely new lever: instead of spreading the counterforce uniformly over the whole fluid, they put it inside a small sphere whose center is displaced by d along the active force direction. That local displacement is the new knob, and it does change the measured flow topology around the chain. The head/tail activity distinction and the stiff-chain inversion (hd− gives puller, td+ gives pusher) are supported by per-monomer force analysis, which is the best part of the paper: it shows that for stiff chains the active force is transmitted along the backbone and the effective force dipole is set by the propagation, not by the local geometry of the applied force.\n\nThe simulations themselves are careful: 48 runs, long averaging, multiple observables (swimming speed, MSD, radius of gyration, contact maps, flow fields) that are internally consistent. The PRW comparison is fine as a consistency check, though it is not an independent test since v0, τr and Dpass all come from the same runs. That should be stated as such, not as a prediction.\n\nThe soft spots are mostly about how much weight the central \"tuning\" claim can carry. The far-field dipole fit is only done for d=±4.5σ, where the counterforce sphere is well separated from the chain. For the d=±σ cases that are the main demonstration in Fig. 6, the classification into pusher/puller is read off from near-field flows. There the counterforce volume overlaps the polymer, the solvent density is perturbed by ~10%, and the system is not strictly torque-free. The torque averages to zero, so I would not expect a systematic rotlet in the averaged fields, but the finite-size of the counterforce volume does add higher multipoles that could affect the sign of the flow at the distances shown. The claim would be much stronger if they showed a far-field fit (or at least a decay check) for the d=±σ cases, or benchmarked the local counterforce method against a known swimmer. As it stands, the inversion for stiff chains is plausible and supported by the monomer-resolved forces, but the quantitative dipole characterization is only established for the d=±4.5σ geometry.\n\nMinor points: no error bars on the fitted dipole strengths; no code/data deposit (they offer data on request, which is the field norm but not great); the fully flexible tail-active \"cat's tail\" motion is a nice observation but not analyzed in depth. These are easily addressable.\n\nWho is this for? People doing MPCD simulations of active filaments or designing synthetic swimmers. It deserves a serious referee: the method is new, the findings are interesting, and the main gap is quantitative validation of the near-field classification, not a fatal flaw. I would send it to review with a request for far-field analysis or an appropriate caveat.","headline":"New local counterforce method for MPCD active polymers is a real addition, and the stiff-chain pusher/puller inversion is plausible, but the central tuning claim needs far-field support for the d=±σ cases before I'd fully buy it.","tokens_in":19049,"tokens_out":3553,"would_cite":true,"duration_ms":37937,"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 flow field around a stiff active polymer is set by where the balancing counterforce sits relative to the chain's center, not by the local force dipole.","keywords":["active polymers","multi-particle collision dynamics","hydrodynamic interactions","pusher/puller flow fields","force dipole","force-free swimmer","bending stiffness","polymer conformation"],"falsifier":"Simulate a stiff 10-mer with d varied continuously from +4.5σ to −4.5σ and track the far-field ϱ^-2 coefficient: the force-propagation picture predicts the sign flips exactly when the counterforce center crosses the chain's center of mass, with near-zero dipole strength at the crossing; a sign flip at some other point, or no sign flip, would falsify it.","tokens_in":18162,"feed_emoji":"🌀","tokens_out":6725,"duration_ms":66617,"temperature":0.7,"pith_summary":"The paper introduces a way to simulate an active polymer chain in a fluid so that the chain's total force on the fluid is exactly zero, by pushing back on a small, movable sphere of solvent rather than on the whole fluid at once. Using this method, the authors show that the chain's flow field can be switched between a 'puller' pattern (fluid pulled in along the swimming axis) and a 'pusher' pattern (fluid pushed out along that axis) just by choosing where the counterforce sits. They then find something unexpected: for stiff chains, this switch is decided by how the active force travels through the polymer backbone to the chain's center, not by the orientation of the local force pair, so naive expectations about which end is active fail. Along the way they show that a head-active monomer straightens the chain like a stiff rod, while a tail-active monomer crumples it, and that this conformational difference survives even when hydrodynamic interactions are switched off.","feed_headline":"Counterforce placement tunes active-polymer flow from puller to pusher","feed_subtitle":"For stiff chains the flow type follows force travel through the polymer, not the local force pair.","key_machinery":"The central mechanism is the movable counterforce volume: a sphere of diameter 2σ, centered at r_cf = R_iact + d e_f, in which the reaction force −F_a is distributed equally among the n_cf solvent particles present. Varying d (chosen here as ±σ or larger values such as ±4.5σ) shifts where the fluid is pushed back, changing the sign and center of the effective force dipole. Work done: it enforces a strictly force-free polymer–fluid system while providing a single tunable parameter that selects the hydrodynamic flow type; combined with bond force transmission in the chain, it explains why stiff chains show flow fields centered at the chain's center of mass and why flexible chains behave differ","core_discovery":"By specifying the position of the counterforce volume relative to the active monomer—ahead or behind the active force—the authors obtain four models (head-active with counterforce before/behind, tail-active with counterforce before/behind) and show that the counterforce position alone can tune the hydrodynamic flow field of the active polymer between puller and pusher types. The deeper discovery is that the resulting flow type is not simply fixed by the sign of the local force dipole: for stiff chains, the active force is transmitted almost instantaneously through the bonds to all monomers, so the effective dipole emerges at the chain's center and its sign is set by the counterforce's displa","pith_inferences":["If the flow type is set by counterforce position relative to the chain center of mass for stiff chains, then continuously sweeping d through the chain's center should flip the far-field dipole sign at the crossing point—an easily testable prediction for simulations.","The same local-counterforce idea could generalise to other flexible active objects (e.g., ring polymers or filaments with distributed activity): the effective multipole order and flow direction should be predictable from the force transmission matrix along the contour, not just the local active site.","The decoupling of conformation from flow type—since counterforce placement changes swimming speed but barely changes chain shape—suggests design rules for artificial microswimmers where propulsion efficiency and far-field disturbance can be optimised separately.","Because the paper finds hydrodynamic interactions barely affect chain conformation, the pusher/puller tuning might be observable in a lattice-Boltzmann or boundary-element implementation as well, provided the reaction force is applied over a similar local volume; the 10% density variations near the dipole are an MPCD compressibility artifact that more accurate solvers could remove."],"forward_implications":["Swimming speed of stiff active polymers grows linearly with active force, while flexible polymers show a nonlinear, activity-site-dependent relation.","Head activity straightens the chain (activity-induced stiffening), increasing orientation persistence; tail activity crumples the chain locally and speeds up backbone decorrelation, including a 'cat's tail' sub-diffusive regime in flexible chains.","Conformational and dynamic differences between head- and tail-active polymers persist even when hydrodynamic interactions are turned off.","For stiff chains, placing the counterforce near a monomer slows the polymer, whereas placing it outside (hd+ and td−) yields faster swimmers, and in all cases the far field is a force dipole whose direction is set by counterforce placement relative to the chain center.","The measured force-dipole strengths for head-active d = −4.5σ and tail-active d = +4.5σ are nearly equal in magnitude but opposite in sign, confirming the rod-like force transmission picture."],"supporting_citations":[{"why":"Supplies the multi-particle collision dynamics (MPCD) solvent algorithm that includes thermal noise and hydrodynamic interactions.","marker":"[34]"},{"why":"Earlier active-filament simulation applying a counterforce within MPCD collision cells; the method this work modifies by making the counterforce local and tunable.","marker":"[52]"},{"why":"Wet active filaments where the counterforce is applied globally to the whole system; the baseline this work contrasts with local counterforce.","marker":"[53]"},{"why":"Gives the pusher/puller classification, the persistent-random-walk model, and the far-field dipole decay used to extract dipole strengths.","marker":"[70]"},{"why":"Provides the low-Reynolds swimmer flow-field description that motivates pusher/puller identification.","marker":"[71]"},{"why":"Reviews hydrodynamic interactions in polymer solutions, supporting the claim that hydrodynamic interactions can qualitatively change polymer dynamics and justifying the HI-free comparison.","marker":"[33]"}],"fun_headline_variants":["Counterforce position sets active polymer flow type","For stiff chains, counterforce site decides puller vs pusher flow","Where counterforce pushes fluid tunes polymer flow direction","Active polymer flow flips with counterforce placement"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The whole flow-tuning picture rests on treating the small sphere of solvent that receives the counterforce as a faithful stand-in for a real swimmer's reaction force, even though the system is not strictly torque-free and the solvent's density varies by about 10% near the force dipole.","fun_headline_variants_meta":{"raw":{"variants":["Counterforce position sets active polymer flow type","For stiff chains, counterforce site decides puller vs pusher flow","Where counterforce pushes fluid tunes polymer flow direction","Active polymer flow flips with counterforce placement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000337,"raw_usage":{"total_tokens":1695,"prompt_tokens":735,"completion_tokens":960,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":897}},"tokens_in":479,"tokens_out":960,"duration_ms":10194,"temperature":1.0,"reasoning_tokens":897,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T16:12:23.481356+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a stiff 10-mer with d varied continuously from +4.5σ to −4.5σ and track the far-field ϱ^-2 coefficient: the force-propagation picture predicts the sign flips exactly when the counterforce center crosses the chain's center of mass, with near-zero dipole strength at the crossing; a sign flip at some other point, or no sign flip, would falsify it.","supporting_citations":[{"cited_title":"Malevanets \\ and\\ author R","cited_arxiv_id":null,"evidence_quote":"Supplies the multi-particle collision dynamics (MPCD) solvent algorithm that includes thermal noise and hydrodynamic interactions."},{"cited_title":"Elgeti \\ and\\ author G","cited_arxiv_id":null,"evidence_quote":"Earlier active-filament simulation applying a counterforce within MPCD collision cells; the method this work modifies by making the counterforce local and tunable."},{"cited_title":"Conformation and dynamics of wet tangentially-driven active filaments","cited_arxiv_id":"2407.17602","evidence_quote":"Wet active filaments where the counterforce is applied globally to the whole system; the baseline this work contrasts with local counterforce."}],"review_version":1}