{"id":"ba5dfcb0-b1b8-437d-8b3e-99f7a70e579d","arxiv_id":"2512.14478","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Spatially modulated tangential activity switches a semiflexible polymer between collapsed and swollen conformations, with a mode-number threshold.","lead":"This paper uses a mathematical model of a self-propelled chain to show that when the propulsion strength varies in a wave along the chain, the chain can shrink into a compact ball or stretch out, depending on the wave's wavelength. The result offers a practical lever for shaping artificial active filaments and for interpreting how cells might pattern motor proteins on the cytoskeleton.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The uncontrolled mean-field closure (Eq. 19) and constant-Θ approximation (Eq. A24) control the sign of the predicted R_G,2 correction; the theory-simulation threshold mismatch (m*=4 vs m*=2) leaves the central quantitative claim unverified at Pe→0.","rationale":"I read the paper in good faith: the idea that spatially modulated tangential activity can switch a semiflexible polymer between compact and swollen conformations is plausible, and the simulations do show a clear qualitative effect. The analytical expansion is a real contribution, and the paper is explicit about several limitations. However, the central quantitative claim—that Eqs. (30) predict which modes collapse and which swell—depends on the mean-field decoupling of Eq. (19) and the subsequent constant-Θ approximation of Eq. (A24). Neither step has a controlled error estimate, and the discrepancy between the analytical threshold m*=4 and the simulated m*=2 is a red flag that these approximations are not merely cosmetic. The reader's verdict (CONDITIONAL) already reflects this concern, and my stress test does not move that verdict; it sharpens the concrete check needed to resolve it: a low-Pe simulation of the same continuum model, not the FENE bead-spring model, with enough precision to extract the coefficient of the Pe² correction. If that check is passed, the conditional concerns would be resolved; if it fails, the quantitative form of the central claim would need to be revised.","tokens_in":20958,"tokens_out":9933,"duration_ms":458469,"concrete_test":"Simulate a bead-spring discretization of the continuum Rouse model in Eq. (1) with harmonic bonds, bending rigidity, and free-end boundary conditions, at Péclet numbers Pe = 0.02, 0.05, 0.1, 0.2 for N=200, l_p=0.5, and modes m=0,...,6. Fit the deviation R_G - R_G,0 as A(m)·Pe² and compare the fitted sign pattern of A(m) with the sign of R_G,2 from Eq. (42). If the fitted A(m) is negative for m=2 or differs near m=4, then the mean-field/constant-Θ pipeline is quantitatively wrong; if it matches, the concern is resolved. Optionally, also recompute R_G,2 using the full s-dependent ⟨Θ(s)⟩_0 of Eq. (A21) instead of the constant b² to isolate the role of Eq. (A24).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result—that Eq. (30) predicts which forcing modes shrink or swell the chain—rests on two uncontrolled approximations. Eq. (19) replaces the average of a ratio, ⟨ r_n·r_j / sqrt(Σ r_m∂ψ_m · r_n∂ψ_n) ⟩, by the ratio of averages, ⟨r_n·r_j⟩ / sqrt(⟨Θ⟩), a mean-field decoupling with no error bound. In the short-persistence-length treatment, ⟨Θ(s)⟩_0 is further replaced by the s-independent constant b² (Eq. A24), discarding the boundary-layer variation that the paper itself derives in Eq. (A22). Both approximations enter the integrals ξ^(0), ξ^(1) that determine the sign and magnitude of the second-order corrections R_G,2 and R_E,2. The observed disagreement between the analytical threshold m*=4 and the simulated m*=2 is exactly the kind of quantitative failure an uncontrolled closure would produce. The paper itself concedes (Introduction) that the simulations could not reach the small-Pe regime required by the expansion, and attributes the threshold mismatch to model idealizations. The claimed 'quantitative confirmation' is therefore not established: the sign trend is reproduced, but the coefficient and threshold are not. If the sign pattern of the Pe² coefficient from Eq. (42) is not reproduced in the Pe→0 limit of the continuum model, the central prediction of a mode-dependent transition is called into question.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a semiflexible active polymer driven by spatially modulated tangential forces. Using a continuum Rouse model with bending rigidity, the authors perform a small-activity expansion of the mode correlations and derive second-order expressions for the gyration radius and end-to-end distance for a sinusoidal forcing profile. The central claim is that spatially structured activity breaks the self-similar scaling of passive chains and produces a mode-dependent transition from shrinking for low forcing modes to swelling for high modes. Langevin dynamics simulations are used to support the analytical predictions. The manuscript also reports a threshold mismatch m*=4 (theory) versus m*=2 (simulation) and concedes that the simulations could not reach the small-Péclet regime required by the expansion.","tokens_in":21426,"tokens_out":4246,"duration_ms":39358,"significance":"If the analytical framework is correct, it provides a useful, parameter-free route to predict how patterned tangential activity reshapes polymer conformations. The paper has clear strengths: the Péclet number is identified rather than fitted, the stretching-coefficient correction D2 is derived from a fixed-contour-length constraint, the expansions are worked out in detail, and simulation data are publicly deposited. However, the quantitative validation is not currently established. The key predictions rest on an uncontrolled mean-field closure and on a constant-⟨Θ⟩ approximation, and the simulation comparison is made at Péclet numbers far outside the expansion regime. The qualitative sign trend is plausible, but the paper's stronger claim of quantitative confirmation is not supported by the evidence presented.","major_comments":[{"comment":"The central approximation replaces the average of a ratio, ⟨ r_n / sqrt(Σ_m r_m² (∂_s ψ_m)²) ⟩, by the ratio of averages, ⟨r_n⟩ / sqrt(Σ_m ⟨r_m²⟩ (∂_s ψ_m)²). This mean-field-type decoupling is uncontrolled and enters the derivation of every higher-order correlation, including Eqs. (30), (42), and (46). No error estimate or consistency check is provided. The discrepancy m*=4 versus m*=2 reported in Sec. IV B is exactly the kind of quantitative failure this closure could produce. The authors should either justify the closure from a small parameter, test it against the exact stationary distribution in a simplified model, or quantify its error numerically.","section":"Sec. II B, Eq. (19)"},{"comment":"In the short-persistence-length treatment, ⟨Θ(s)⟩_0 is replaced by the constant b², discarding the spatial dependence derived in Eq. (A22), including boundary-layer variation. This approximation is used to evaluate ξ^(0) and ξ^(1), which determine the sign and magnitude of R_G,2 and R_E,2. Since the central claim depends on the sign of these corrections, the sensitivity of the threshold m* to this approximation must be examined. For example, using the full expression for ⟨Θ(s)⟩_0, or keeping the next term in the κN expansion, would show whether the predicted m*=4 is robust.","section":"Appendix A 2, Eq. (A24)"},{"comment":"The claimed 'quantitative confirmation' is not supported by the data shown. The simulations use f_a=0.01 and 10 (Pe≈0.14–0.20 and Pe≈141–200), while the expansion requires Pe_l≪1. The upper panel of Fig. 6 is in a marginal regime where the second-order correction is very small, and the lower panel is far outside the expansion regime. Moreover, the theoretical curves in Figs. 6(b,c) and 8(b,c) are computed at Pe=1 and are not overlaid with the simulation data. The paper should either present a direct quantitative comparison at a Péclet number where the expansion is expected to hold, or explicitly restrict the claim to qualitative agreement of the sign trend.","section":"Sec. IV A, Fig. 6"},{"comment":"The theory and simulations disagree on the leading-order sign of R_E,1 for m=3: the text states that the theory predicts a decrease of R_E for m=3, 'in contrast to simulation results.' This is attributed to higher-order corrections, but no evidence is provided that the O(Pe) prediction becomes correct in the Pe→0 limit of the discrete model. Since this is a leading-order parity-dependent prediction, it should be tested directly or the claim should be softened.","section":"Sec. IV B, end-to-end distance"}],"minor_comments":[{"comment":"The derivation of the gyration tensor appears to contain an algebraic slip: the term −⟨r0,i r0,j⟩ is written and then dropped, but the orthogonality of the cosines makes the center-of-mass contribution vanish for different reasons. Please clarify the notation and the cancellation.","section":"Eq. (39)"},{"comment":"The Péclet number is denoted Pe_l, Pe_m, and Pem in different places (Eqs. (32), (A36), and Fig. 6). Please unify the notation and define the subscript.","section":"Notation"},{"comment":"The sentence 'in both the upper and lower sub-panels' is misleading because the two panels correspond to very different Péclet numbers; the qualitative statement should be separated from the quantitative comparison.","section":"Sec. IV A"},{"comment":"The theoretical curves are shown for Pe=1, but the text sometimes reads as if they were directly comparable to the simulation data at Pe=0.141–200. Please state explicitly that the analytical plots are illustrative of the weak-activity expansion.","section":"Sec. IV A, Figs. 6-8"},{"comment":"The shift of domain from [−N/2,N/2] to [0,N] is performed without stating the corresponding boundary conditions in terms of the shifted eigenfunctions. A brief comment would help readers verify the trigonometric representation.","section":"Sec. II D, Eq. (27)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid formal structure and the qualitative phenomenon is interesting, but the load-bearing quantitative comparison is not yet established. The uncontrolled closure in Eq. (19) and the constant-Θ approximation in Eq. (A24) control the sign that is central to the paper. I would support publication if the authors either provide a concrete validation of the closure (e.g., a simplified exactly solvable test) or explicitly reframe the claims as qualitative and add a careful discussion of the approximation error. The current abstract's 'quantitatively confirm' overstates what the simulations show."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper you asked about is worth a look, but take the abstract's \"quantitatively confirm\" with a grain of salt. What's genuinely new: the authors push a continuum Rouse model with bending rigidity into a systematic weak-activity expansion, keeping the full eigenfunction basis, and they get closed forms for mode correlations, gyration radius, and end-to-end distance under a single-mode sinusoidal tangential force. The mode-dependent transition—low modes collapse the chain, high modes swell it—is physically plausible and, as far as I know, not in the literature in this analytic form. They also cleanly identify the dimensionless control parameter Pe ~ f b sqrt(N)/kBT and show parity effects: odd forcing modes give a first-order contribution to the end-to-end distance, even modes do not. That's a real step beyond the uniform-activity treatment.\n\nThe derivation is careful in its own terms. No free parameters are fit to simulation; the Péclet number follows from the expansion, and the stretching-coefficient correction is fixed by contour-length conservation. They post the simulation data. That all counts for something.\n\nThe soft spots are real but not damning. The mean-field closure in Eq. 19—replacing the average of a ratio by the ratio of averages in the tangent-vector denominator—is uncontrolled, and the companion constant-Θ approximation in App. A also discards boundary-layer structure. These enter the integrals that set the sign of the second-order correction, so the threshold m* is exactly where an uncontrolled closure would bite. The authors admit that their Langevin simulations could not reach the small-Pe regime the expansion requires; the comparison is at moderate Pe. The predicted threshold (m*=4) and the simulated one (m*=2) differ by a factor of two. I don't think this kills the central idea, but it does mean the paper has not quantitatively confirmed its own headline prediction.\n\nThe paper is for active soft matter theorists, and it deserves a serious referee. I'd send it out, with the request that the authors either error-control the closure—even numerically, by evaluating the exact average on the simulated trajectories—or tone down the confirmation language to \"qualitative agreement\" for the sign of the effect. If they address that, it's a solid contribution. As is, it's a promising paper with an overstated validation.","headline":"Solid analytic expansion for patterned active polymers, but the abstract overstates the confirmation: an uncontrolled closure and a theory-simulation threshold mismatch leave the central quantitative claim unverified.","tokens_in":21770,"tokens_out":2540,"would_cite":true,"duration_ms":22326,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.40.-a","82.35.Lr"],"model":"deepseek-v4-flash","headline":"Sinusoidally modulated tangential activity can collapse a semiflexible polymer into a globule at low mode numbers or swell it into alternating stretched and compressed segments at higher modes, with the switch near mode 4 analytically and m","keywords":["active polymers","tangential activity","Rouse model","semiflexible polymers","gyration radius","end-to-end distance","mode-dependent conformations","Péclet number"],"falsifier":"Run Langevin simulations of a tangential-active polymer at Pe_m ≪ 1 (for example f_a ≈ 10⁻³ k_BT/σ with N = 400, so Pe_m ≈ 0.02) and measure R_G² as a function of the forcing mode m. The analytic prediction is a sign change in the correction at m* = 4; observing the sign flip at m* = 2 at such small Péclet numbers would falsify the quantitative closure.","tokens_in":20890,"feed_emoji":"🧬","tokens_out":5986,"duration_ms":51110,"temperature":0.7,"pith_summary":"The paper sets out to show that spatially patterned tangential activity is a control knob for polymer conformation even when the active force is weak. Using a semiflexible Rouse model expanded in small activity, it derives mode-correlation formulas that predict a sign change in the gyration-radius correction at a threshold forcing mode m*: low-mode forcing (m = 0,1) shrinks the chain into globule-like shapes, while higher modes produce alternating stretched and compressed segments and a swollen chain. The same expansion explains how a polymer can be compact in gyration radius yet extended end-to-end, because the two metrics weight modes differently. The authors support the analytic predictions with Langevin simulations, while noting that the simulations run at finite rather than perturbatively small activity and that the analytic threshold (m* = 4) differs from the simulated one (m* = 2).","feed_headline":"Active polymers shrink or swell depending on the activity pattern","feed_subtitle":"Low spatial modes collapse a semiflexible chain into a globule; higher modes stretch it into alternating segments.","key_machinery":"The machinery is a continuum Rouse model with bending rigidity, expanded in the passive-chain eigenmodes ψ_n(s) = sqrt(2/N) cos(π n s/N), with eigenvalues ζ_n ∝ n²(1+Λ_p n²) and a single-mode tangential active force f(s) = f_m cos(π m s/N). To make the expansion tractable, the paper uses a mean-field decoupling of the normalized tangent vector: the average of r_n / |∂s r| is replaced by ⟨r_n⟩ / sqrt(⟨Θ(s)⟩), and ⟨Θ(s)⟩ is then approximated by its peak value b². This yields closed perturbative formulas for ⟨r_i·r_j⟩ at orders 0, 1, and 2 in the Péclet number Pe_m = f_m b √N / (k_B T), which is the combination that identifies the weak-activity regime uniformly in chain length.","core_discovery":"The paper's central claim is that the spatial mode of tangential self-propulsion determines whether a semiflexible polymer shrinks or swells. Under weak sinusoidally modulated propulsion, the second-order correction to the gyration radius flips sign at a threshold mode number: uniform or low-mode forcing produces compact, globule-like conformations, whereas higher modes generate alternating stretched and compressed segments that make the chain globally swollen. The analytical expressions (Eqs. 30) for mode correlations, gyration radius, and end-to-end distance show that activity breaks self-similar scaling and that different size metrics respond differently, allowing conformations that are c","pith_inferences":["If the threshold m* is governed by the normalized eigenvalues z̄_i = i²(1+Λ_p i²), one could tune persistence length or polymerization degree to place a desired spatial mode on either side of the transition, effectively programming polymer shape by choosing which modes are active.","The analytic-versus-simulated threshold discrepancy suggests that the mean-field closure underrepresents single-mode effects; a variational or renormalized treatment might shift m* and improve quantitative agreement at finite Péclet numbers.","The same weak-activity expansion could be extended to multi-mode or time-modulated forcing, where interference between modes might create regimes that neither mode produces alone, such as collapse driven by two modes that individually swell the chain.","The parity argument for rings is directly testable: simulations of ring polymers at small Péclet should show no leading-order gyration-radius change, with higher-order corrections and swelling emerging only as the activity grows."],"forward_implications":["Weak, spatially structured activity can induce both collapse and swelling of a single polymer, making force patterning a practical alternative to changing solvent quality.","Because the relevant Péclet number grows as √N at fixed force, longer chains respond much more strongly to the same active forcing.","A chain may be compact by gyration radius but extended end-to-end, so measuring either quantity alone can miss the conformational state.","For ring polymers, leading-order activity corrections vanish by parity (only even modes contribute), so topology sharply constrains the linear-polymer result.","Non-conservative activity yields a quadratic dependence of the gyration radius on force strength, distinguishing active collapse from conservative attraction, which is linear."],"fun_headline_variants":["Propulsion mode flips polymers from globule to swollen","Active polymer size follows spatial mode of propulsion","Spatial modes decide active polymer collapse or stretch","Patterned activity controls polymer compaction and expansion","Propulsion wavelength sets polymer conformation state"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing approximation is the mean-field replacement of the average of the normalized tangent ratio by the ratio of averaged quantities (Eq. 19), together with treating ⟨Θ(s)⟩ as a constant; the paper's own simulations could not reach the perturbatively small forces and show a shifted threshold (m* = 4 analytically vs m* = 2 numerically), so the closure is plausible but uncontrolled.","fun_headline_variants_meta":{"raw":{"variants":["Propulsion mode flips polymers from globule to swollen","Active polymer size follows spatial mode of propulsion","Spatial modes decide active polymer collapse or stretch","Patterned activity controls polymer compaction and expansion","Propulsion wavelength sets polymer conformation state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000682,"raw_usage":{"total_tokens":2903,"prompt_tokens":685,"completion_tokens":2218,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":429,"completion_tokens_details":{"reasoning_tokens":2148}},"tokens_in":429,"tokens_out":2218,"duration_ms":15160,"temperature":1.0,"reasoning_tokens":2148,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T16:01:09.916290+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run Langevin simulations of a tangential-active polymer at Pe_m ≪ 1 (for example f_a ≈ 10⁻³ k_BT/σ with N = 400, so Pe_m ≈ 0.02) and measure R_G² as a function of the forcing mode m. The analytic prediction is a sign change in the correction at m* = 4; observing the sign flip at m* = 2 at such small Péclet numbers would falsify the quantitative closure.","supporting_citations":[],"review_version":1}