{"id":"c29d2147-829d-43aa-9920-e2a392261225","arxiv_id":"2505.17539","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Tangential polar activity amplifies shear thinning of flexible polymers, changing the viscosity exponent from 1/2 to 4/3 and the tumbling-frequency exponent from 2/3 to 1/3.","lead":"Using computer simulations and an analytical model, the authors show that tangentially self-propelled (polar) flexible polymers under linear shear flow stretch more along the flow, shrink more in the transverse direction, and shear-thin much more strongly than passive polymers. The shear viscosity drops with a power-law exponent of about 4/3 in an intermediate regime, compared with 1/2 for passive polymers, before crossing over to passive behavior at very high shear rates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The virial stress in Sec. VI/Eq. (E15) includes the non-zero net active force with lab-frame positions; the resulting eta_p may be origin-dependent, so the central -4/3 viscosity exponent could be an artifact. A center-of-mass or force-free check is needed.","rationale":"I read the paper as claiming a model-specific but generalizable scaling law: in a dry, overdamped, excluded-volume active polar polymer under shear, the intermediate-activity regime has eta_p/eta_p^0 ~ Wi_Pe^{-4/3} and nu_t tau_r ~ Wi_Pe^{1/3}, with activity-independent collapse. The simulation data are internally consistent, the tumbling scaling follows from the measured <G_yy> via the stated argument, and the analytical model provides partial independent support, although its Lagrange multiplier mu is not self-consistently determined. The acknowledged lack of hydrodynamics is a transferability caveat, not an internal flaw, so I do not base the attack on it. The most load-bearing internal soft spot is the virial stress definition: Eq. (E15) sums lab-frame r_i against forces whose total does not vanish. For passive polymers the cancellation makes the formula safe; with tangential active forces the net active force is nonzero and shear-aligned, so the active term can inject an origin-dependent monopole contribution. If that term is present, the central intrinsic-viscosity exponent is not a material property. The proposed check is decisive and requires no new physics: recompute the stress with center-of-mass-relative coordinates and with a shifted box. If the exponent survives, the concern is retired; if not, the paper needs at least a major revision or a rederivation of the stress. The reader's weakest_assumption (dry model / hydrodynamics) is related but not identical, so I mark partial agreement. Since the reader already made the verdict conditional on model-transfer questions, and my concern is about internal correctness, I keep the conditional stance rather than escalate.","tokens_in":16103,"tokens_out":24724,"duration_ms":311481,"concrete_test":"Recompute eta_p from existing trajectories using the center-of-mass-relative virial sigma_xy = -(1/V) sum_i <(F_i^int + F_a_i)(r_i,y - R_cm,y)>, and also shift the periodic box in y by half a box length, for Pe=100, N=200, and Wi_Pe in [10, 10^3]. If the log-log slope moves away from -4/3 by more than about 0.1, or if the shifted-box value changes by a comparable amount, the central viscosity exponent is not a well-defined intrinsic property; if the slopes and values are unchanged, the virial implementation is acceptable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the computed intrinsic viscosity. Section VI defines sigma_xy = -sum_i <(F_i^x + F_a_i^x) r_i^y>/V, and App. E Eq. (E15) repeats this with lab-frame positions. For passive chains, sum_i F_i^int = 0, so this virial expression is origin-independent. Here, however, sum_i F_a_i = f_a sum_i t_i = f_a (R_N - R_1)/l_0, which is not zero and, under shear, has a nonzero mean along x. The active-force term in the virial therefore contains an origin-dependent monopole contribution proportional to the center-of-mass y-coordinate. The text does not state that r_i is measured relative to the center of mass (contrast Eq. (E13)-(E14) for the gyration tensor). If the code uses absolute coordinates, eta_p depends on the arbitrary placement or size of the periodic box, and the reported eta_p/eta_p^0 ~ Wi_Pe^{-4/3}, as well as its similarity to <G_yy>, could be a coordinate-convention artifact rather than a rheological exponent. This is a concrete internal issue distinct from the acknowledged absence of hydrodynamics, and it can be checked from the same trajectories by recomputing the stress in two coordinate conventions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses coarse-grained Brownian dynamics simulations of a dry, overdamped, self-avoiding flexible polymer whose monomers are driven by tangential active forces, and compares the results with a Gaussian-chain analytical model developed in Ref. [49] and rederived in Appendix E. For linear shear flow, it reports that in an intermediate, activity-dominated Weissenberg-number regime the polymer gradient size shrinks as ⟨G_yy⟩/⟨G_yy^0⟩ ∼ Wi_Pe^{-4/3}, the intrinsic viscosity thins as η_p/η_p^0 ∼ Wi_Pe^{-4/3}, the alignment angle behaves as tan(2χ) ∼ Wi_Pe^{-1}, and the tumbling frequency grows only as ν_t τ_r ∼ Wi_Pe^{1/3}, before a crossover to the passive exponents (1/2 for viscosity, 2/3 for tumbling) at large Wi_Pe. The zero-shear viscosity is found to decrease with activity and to collapse as a function of N_m Pe.","tokens_in":16440,"tokens_out":9397,"duration_ms":69335,"significance":"If the reported exponents are correct, the paper establishes a qualitatively new scaling regime for dilute polar active polymers: tangential activity changes the rheological and dynamical exponents by large factors relative to passive chains, and the onset shear rate is controlled by Pe. The claim is concrete and falsifiable, and the manuscript includes several internal consistency checks: two independent estimators of the tumbling time agree; the main power laws are shown for N_m = 200, 400, and 1000; and the near-proportionality of η_p and ⟨G_yy⟩ over four decades in Wi is used as a cross-check. The explicit statement that hydrodynamic interactions are neglected is an appropriate caveat. The main weaknesses are that the central stress calculation has a possible origin-dependence issue and that the exponent values are quoted without quantitative uncertainty; both are fixable without changing the scope of the paper.","major_comments":[{"comment":"The shear stress is defined as σ_xy = -∑_i ⟨(F_i^x + F_{a,i}^x) r_i^y⟩/V, and the same defining relation is used in Eq. (E15). For the active forces in this model, ∑_i F_{a,i}^x = f_a (R_N - R_1)_x/l_0, which is not zero and has a nonzero mean under shear. Consequently, the active contribution to this virial expression is not invariant under a shift of the origin of the y-coordinate unless r_i is understood to be measured relative to the polymer center of mass. The manuscript explicitly uses CM-relative coordinates for the gyration tensor in Eq. (4) and Eq. (E13), but no such statement is made for the stress. Because the central η_p/η_p^0 ∼ Wi_Pe^{-4/3} result is extracted from this σ_xy, please state the coordinate convention actually used in the simulations and recompute η_p using r_i - r_cm (or otherwise prove the origin independence of the reported values). This is a load-bearing check and should be reported explicitly.","section":"Section VI, Eq. (E15)"},{"comment":"The paper quotes the asymptotic exponents 4/3 and 1/3 without regression intervals or any uncertainty estimate. Since the central claim is specifically that the exponents differ from the passive values 1/2 and 2/3, please include quantitative power-law fits for each Pe, with confidence intervals and the Wi_Pe ranges over which the fit is performed, and state the statistical uncertainty of the extracted exponents.","section":"Figs. 2(a), 5(a), 6"},{"comment":"The analytical comparison is performed at fixed r_d = 25 for Pe = 150 rather than by solving the Lagrange multiplier μ(γ̇) self-consistently through Eqs. (E4) and (E10). The theory therefore provides a consistency check at chosen parameter values, not an independent derivation of the −4/3 exponent or the 1/3 tumbling exponent. Please state this limitation wherever the analytical model is said to 'support' or 'confirm' the simulation results.","section":"Appendix E"}],"minor_comments":[{"comment":"The phrase 'aside form end effects' should be corrected to 'aside from end effects'.","section":"Section V"},{"comment":"The phrase 'suggest a root to the development' should be corrected to 'suggest a route to the development'.","section":"Section VIII"},{"comment":"In Eq. (E15), 'viral stress' should be 'virial stress'.","section":"Appendix E"},{"comment":"The duplicated phrase 'for for both properties' should be reduced to 'for both properties'.","section":"Section VI"},{"comment":"The inset uses N_m Pe as the horizontal axis; a brief sentence in the caption explaining why this is the natural scaling variable would improve readability.","section":"Figure 5(b)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe headline result here—η_p/η_p^0 ∼ Wi^{-4/3} and ν_t τ_r ∼ Wi^{1/3} for tangentially driven polar polymers under shear—is genuinely new as far as I can tell, and it deserves a serious referee. I read this as a simulation paper with a supporting analytical model, not the reverse.\n\nWhat is strong: three chain lengths, five activities, two independent tumbling definitions that agree, a fluctuation-in-activity check, and a four-decade comparison between ⟨G_yy⟩ and η_p. The data collapse in the intermediate regime and the crossover back to passive behavior are convincing. The scaling argument linking ν_t to ⟨G_yy⟩ is internally consistent: if ⟨G_yy⟩ ∼ Wi^{-4/3}, the tumbling exponent 1/3 follows from τ_t ∼ 1/(γ√⟨G_yy⟩). That part looks solid. The analytic model, though, is not an independent prediction: it carries over Ref [49] and fixes the Lagrange multiplier by hand (rd = 0, 25) instead of solving the constraint, so it is not a first-principles derivation of 4/3.\n\nThe soft spot that bothers me most is the virial stress. Main text and Eq. (E15) write σ_xy = -(1/V) Σ_i ⟨(F_i^x + F_ai^x) r_i^y⟩ with lab-frame r_i. For passive chains Σ_i F_i = 0 and the expression is origin-independent. Here Σ_i F_ai = fa(r_N - r_1)/l0, which is nonzero and has a nonzero mean along x under shear, so the active term in the virial contains an origin-dependent contribution proportional to the center-of-mass y-coordinate. The paper nowhere says r_i is measured relative to the COM, and it contrasts with the gyration tensor definition in Eqs. (E13)-(E14), which explicitly use Δr_i. If the code used absolute coordinates, η_p depends on where the polymer sits in the box, and the −4/3 viscosity law could be a coordinate artifact. This is checkable from existing trajectories by recomputing the stress in two conventions, and it must be checked before the central claim is taken as established. The analytical expression E16 also seems to drop the active-force term, which is either an inconsistency or a hint that the simulations should not be including it either.\n\nMissing error bars on exponents and no code/data availability are minor but real. The dry-polymer approximation is acknowledged in the Summary and is a model limitation, not a disqualifier. Citation pattern is fine; Ref [49] is genuinely the basis of the theory.\n\nAudience: active-polymer rheology and single-polymer dynamics. I would send it to review, with a request for the coordinate-invariance check, error bars, and a clear statement of what goes into the stress. The tumbling and conformation results may well hold; the viscosity exponent is the one I would not yet bet on.","headline":"Genuinely new active-polar-polymer shear exponents, and the simulations look mostly careful, but the viscosity route has an origin-dependence problem that needs a coordinate-invariance check before the 4/3 exponent is trusted.","tokens_in":16952,"tokens_out":8224,"would_cite":true,"duration_ms":68806,"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":"Tangentially propelled polar polymers under shear flow shear-thin with viscosity exponent 4/3 and tumble with exponent 1/3 until shear dominates.","keywords":["active polymers","shear thinning","tumbling dynamics","polar activity","Weissenberg number","Péclet number","intrinsic viscosity","excluded volume"],"falsifier":"Measure the intrinsic viscosity and tumbling frequency of a dilute suspension of tangentially propelled flexible polar filaments (or simulate the same model with full hydrodynamic interactions) over $10 < Wi_{Pe} < 1000$ at $Pe > 5$; finding viscosity slopes other than $-4/3$ or tumbling slopes other than $1/3$ in the window where passive chains show $-1/2$ and $2/3$ would settle the claim negatively. A sharper internal check is the paper's prediction that $\\langle G_{yy}\\rangle$ and $\\eta_p$ track each other within 25–35% across four decades of $Wi_{Pe}$; a measurement where those two curves diverge would falsify the proposed mechanism.","tokens_in":15917,"feed_emoji":"🌀","tokens_out":7559,"duration_ms":61709,"temperature":0.7,"pith_summary":"This paper asks what happens when a flexible polymer is actively propelled along its own contour while a linear shear flow is applied. Combining coarse-grained simulations with an analytical model, it establishes a new activity-dominated rheological regime: the polymer's intrinsic viscosity falls as $\\eta_p/\\eta_p^0 \\sim Wi_{Pe}^{-4/3}$ with shear rate, far steeper than the $Wi_{Pe}^{-1/2}$ of a passive chain, while its tumbling frequency grows only as $\\nu_t \\tau_r \\sim Wi_{Pe}^{1/3}$ rather than $Wi_{Pe}^{2/3}$. The mechanism is that polar activity amplifies shear-induced compression in the gradient direction, so the polymer presents a smaller cross-section to the flow and overturns more slowly. At sufficiently high shear rates, ordinary flow dominates and the passive exponents return. If the claim is right, tangential polar activity is a generic route to much stronger shear thinning in dilute polymer solutions.","feed_headline":"Polar activity steepens polymer shear-thinning to a 4/3 power law","feed_subtitle":"Tangentially driven chains shrink across the flow, cutting viscosity faster and slowing tumbling until strong shear takes over.","key_machinery":"The load-bearing machinery is a coarse-grained overdamped Langevin dynamics in which each monomer experiences an active force $f_a(\\hat{t}_{i+1}+\\hat{t}_i)/2$ along the bond directions together with a linear shear flow and excluded-volume forces. The companion analytical model solves a discretized Gaussian active polar polymer by expanding in the eigenfunctions of a non-symmetric matrix; its eigenvalues $\\xi_n$ encode the activity-dependent mode relaxation, and both the gyration tensor and the virial stress (hence $\\eta_p$) are written as sums over these modes. The paper's central scaling argument is geometric: during a tumbling event a monomer is convected at velocity $v_x \\approx \\dot{\\gamma} y_t$ with $y_t \\approx \\sqrt{\\langle G_{yy}\\rangle}$, so $\\tau_t \\sim (\\dot{\\gamma}\\sqrt{\\langle G_{yy}\\rangle})^{-1}$; inserting $\\langle G_{yy}\\rangle \\sim Wi_{Pe}^{-4/3}$ gives $\\nu_t \\tau_r \\sim Wi_{Pe}^{1/3}$. Here $Wi_{Pe} = \\dot{\\gamma}\\tau_r(Pe)$ is the Weissenberg number built on the activity-dependent relaxation time and $Pe = f_a l_0/(k_B T)$ measures the active force.","core_discovery":"The central discovery is a qualitative change in the shear-rate scaling of dilute tangentially driven active polar polymers. In the activity-dominated intermediate window, the normalized intrinsic viscosity obeys $\\eta_p/\\eta_p^0 \\sim Wi_{Pe}^{-4/3}$, the gradient-direction radius of gyration shrinks as $\\langle G_{yy}\\rangle/\\langle G_{yy}^0\\rangle \\sim Wi_{Pe}^{-4/3}$, the flow alignment satisfies $\\tan(2\\chi) \\sim Wi_{Pe}^{-1}$, and the tumbling frequency obeys $\\nu_t \\tau_r \\sim Wi_{Pe}^{1/3}$; all quantities then cross over to the passive exponents ($-1/2$, $-1/3$, and $2/3$) once shear dominates over activity. The paper also finds that the zero-shear viscosity of self-avoiding active polymers is reduced by roughly a factor of two relative to passive chains, that this reduction requires excluded-volume interactions, and that zero-shear data collapse onto a universal curve as a function of $N_m Pe$.","pith_inferences":["A direct test would be to measure single-chain viscosity or tumbling of motor-driven cytoskeletal filaments in a microfluidic shear cell; slopes other than $-4/3$ and $1/3$ would signal that hydrodynamics or boundary conditions change the mechanism.","If full hydrodynamic simulations shift one exponent but not the other, the tight link the paper finds between $\\langle G_{yy}\\rangle$ and $\\eta_p$ would break, suggesting the stress and conformation are less directly coupled in real solvents.","The same tangential-propulsion coupling might produce steepened response in extensional or oscillatory flows, but the geometric tumbling argument would need re-derivation for those kinematics."],"forward_implications":["Dilute solutions of tangentially propelled polar polymers should exhibit a viscosity drop with shear much steeper than passive solutions, with the 4/3 exponent visible in the window $10 < Wi_{Pe} < 10^3$.","In the same window, tumbling slows relative to passive chains, making the stretch–recoil period a mechanical signature of polar activity.","For large activity, the intermediate-regime response becomes independent of $Pe$, so different activities collapse onto a single master curve and activity only sets the onset shear rate.","The crossover to passive exponents at large $Wi_{Pe}$ means the active signature can be erased by strong flow, which matters for processing flows."],"supporting_citations":[{"why":"Supplies the analytical eigenfunction model of tangentially driven active polar linear polymers used for the theory curves.","marker":"[49]"},{"why":"Provides the passive-polymer shear-flow baseline: conformation, alignment, and viscosity scaling with Wi.","marker":"[32]"},{"why":"Supplies the passive tumbling and viscosity exponents and the convective tumbling scaling argument extended here.","marker":"[36]"},{"why":"Provides the passive tumbling times and the gyration-tensor correlation method used to extract tumbling frequencies.","marker":"[39]"},{"why":"Gives the active Brownian polymer shear results ($Wi^{-3/4}$) that contrast with the steeper polar-polymer exponent.","marker":"[55]"},{"why":"Documents the activity-induced collapse of tangentially driven active polymers that underlies the zero-shear viscosity reduction.","marker":"[45]"},{"why":"Shows the fore-aft symmetry breaking and inhomogeneous stretching of polar active polymers invoked for the local stretching analysis.","marker":"[62]"},{"why":"Supplies the simulation model with tangential active forces, bond springs, and excluded volume used in the main runs.","marker":"[46]"}],"fun_headline_variants":["Active polar polymers shear-thin with a 4/3 power law","Activity amplifies polymer shear-thinning to a 4/3 exponent","Polar activity sets new shear-thinning exponent for polymers","Shear-thinning exponent of active polymers jumps to 4/3","Active polymers shear-thin steeper: 4/3 power law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes a dry, overdamped polymer with no hydrodynamic interactions and with driving tangential to the bonds; if solvent-mediated flow around the chain materially changes how it aligns, shrinks, or dissipates stress, the 4/3 and 1/3 exponents may not survive in real fluids.","fun_headline_variants_meta":{"raw":{"variants":["Active polar polymers shear-thin with a 4/3 power law","Activity amplifies polymer shear-thinning to a 4/3 exponent","Polar activity sets new shear-thinning exponent for polymers","Shear-thinning exponent of active polymers jumps to 4/3","Active polymers shear-thin steeper: 4/3 power law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000896,"raw_usage":{"total_tokens":3817,"prompt_tokens":856,"completion_tokens":2961,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":2868}},"tokens_in":472,"tokens_out":2961,"duration_ms":17678,"temperature":1.0,"reasoning_tokens":2868,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:45:00.928261+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the intrinsic viscosity and tumbling frequency of a dilute suspension of tangentially propelled flexible polar filaments (or simulate the same model with full hydrodynamic interactions) over $10 < Wi_{Pe} < 1000$ at $Pe > 5$; finding viscosity slopes other than $-4/3$ or tumbling slopes other than $1/3$ in the window where passive chains show $-1/2$ and $2/3$ would settle the claim negatively. A sharper internal check is the paper's prediction that $\\langle G_{yy}\\rangle$ and $\\eta_p$ track each other within 25–35% across four decades of $Wi_{Pe}$; a measurement where those two curves diverge would falsify the proposed mechanism.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the passive-polymer shear-flow baseline: conformation, alignment, and viscosity scaling with Wi."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the passive tumbling and viscosity exponents and the convective tumbling scaling argument extended here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the passive tumbling times and the gyration-tensor correlation method used to extract tumbling frequencies."},{"cited_title":"Mart ´ ın-G´ omez, G","cited_arxiv_id":null,"evidence_quote":"Gives the active Brownian polymer shear results ($Wi^{-3/4}$) that contrast with the steeper polar-polymer exponent."},{"cited_title":"Patra, K","cited_arxiv_id":null,"evidence_quote":"Documents the activity-induced collapse of tangentially driven active polymers that underlies the zero-shear viscosity reduction."},{"cited_title":"Fazelzadeh, E","cited_arxiv_id":null,"evidence_quote":"Shows the fore-aft symmetry breaking and inhomogeneous stretching of polar active polymers invoked for the local stretching analysis."},{"cited_title":"Bianco, E","cited_arxiv_id":null,"evidence_quote":"Supplies the simulation model with tangential active forces, bond springs, and excluded volume used in the main runs."}],"review_version":1}