{"id":"186794a3-58e1-42ea-8dbb-9b077bf80c6c","arxiv_id":"2607.24102","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Transverse-force spectra of exchangeable active multipoles obey Δu k^4 + B k^6, with the crossover set by the density of locally unscreened moments.","lead":"This paper introduces a new way to define hyperuniform order in flowing active matter: instead of asking whether particles are evenly spaced, it asks which part of the forces the fluid actually responds to. In simulation models of exchangeable active pairs, local neutral clusters make the transverse-force spectrum decay as k^6, while unpaired carriers force a k^4 leak and set a finite screening length.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Δu + B k^2 normal form is fit only over a crossover-dominated window; the k^4→k^6 asymptote is an untested extrapolation.","rationale":"The reader's CONDITIONAL verdict already identifies the core weakness: the analytic two-term normal form is fit over k≤0.9, AICc cannot always distinguish it from a power law, and the extrapolation to k→0 is assumed rather than demonstrated. My stress-test sharpens this: the only independent microscopic anchor (ρloc_u–Δu proportionality) constrains Δu but not B, so the k^6 term and the resulting R^-4 velocity decay are not independently verified. The larger-N simulation and the independent cluster-dipole measurement would directly settle whether S1(k) really flattens to a plateau in the deep infrared and whether the k^2 curvature is controlled by neutral clusters. The paper's own SM S5.6 and S8 admit the ambiguity and the finite-range/analyticity boundary, which strengthens the case for requiring these checks before full acceptance. However, the paper already has supporting evidence: two distinct models, finite-size persistence, a local–spectral proportionality, and a velocity prediction with direct/spectral ratios near unity. These are real independent supports, so the claim is not baseless. The concern is that the universal asymptotic form is not yet proven, but the conditional framework the reader used is the appropriate response. Therefore the verdict should remain CONDITIONAL; no change to the reader's assessment is warranted, but the stated conditions should explicitly require the deep-infrared test and the independent B measurement.","tokens_in":24579,"tokens_out":11863,"duration_ms":111033,"concrete_test":"Run Model I at N=16384 (or N=4096 as a first step) at koff=0.06 and 0.4 with the same parameters; this gives L≈256 and k_min≈0.0245, i.e., more than 20× below k×. Compute S1(k) with the exact particle-sum estimator. If the lowest resolved shells lie on the Δu+Bk^2 curve extrapolated from the k≤0.9 fit, the normal form is confirmed in the deep infrared. If S1 continues to fall below the fitted plateau as k→0, the fitted Δu is a finite-window artifact and the claimed k^4 infrared leakage and ξ_scr are not established. As a complementary check, independently measure B from the instantaneous cluster dipole variance, B_ind = (1/(2V))<Σ_c |Σ_a M_a(r_a−R_c)|^2>, using the same snapshots and no spectral fit, and compare to the fitted B; agreement would confirm that the k^2 term is physically controlled by neutral clusters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central law Eq. (7) rests on S1(k)=Δu+Bk^2+O(k^4) holding to k→0. The numerical support is weaker than the abstract implies: all fits use k≤0.9, and at N=1024 the smallest resolved wave number is k_min≈0.098, only ~0.18 k× at koff=0.06. The fit window is therefore dominated by crossover behavior, not by a clean k^4 (plateau) or k^6 (analytic) asymptote. SM S5.6 explicitly says a free power law can fit individual finite-window spectra as well as the two-term crossover, so AICc does not select the normal form. The independent ρloc_u–Δu proportionality (Fig. 2b) validates only the plateau amplitude as a local defect density; nothing independently measures B or the k^2 curvature. B is extracted from the same S1(k) that the law is tested against. A non-analytic correction (e.g., c k^2 ln k or k^{2+ε}) would be absorbed into B/Δu over this window while changing the true k→0 limit and the predicted R^-4 velocity term. The paper acknowledges the analyticity/finite-range assumption (SM S8) but never tests it with a deeper infrared probe. Without that, the universal k^4/k^6 separation and ξ_scr=sqrt(B/Δu) are an extrapolation rather than a measured asymptotic statement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces 'response-selected hyperuniformity' as a property of a source–response pair rather than of a bare density or continuum field. For an incompressible active fluid whose units carry signed first moments, the transverse projector selects a force sector whose spectrum is argued to have the normal form S1(k)=Δu+Bk^2+O(k^4) (Eq. 6). This yields the central law χ̃_T^f(k)=Ŵ(k)^2(Δu k^4+B k^6) (Eq. 7), with a crossover k×=sqrt(Δu/B), a screening length ξ_scr=k×^(−1), and a strict hyperuniform velocity limit only when Δu=0. The claim is tested in two stochastic many-body models: an exchangeable valence-one active-multipole fluid and a fixed-partner active-molecule model. Both are reported to exhibit the same S1 and force-spectrum forms. The paper further argues that complete partner renewal preserves the multipole inheritance, that the plateau Δu is controlled by an independently measured local unpaired-moment density ρ_loc^u (Fig. 2b), and that the resulting velocity variance has the R^(−2)/R^(−4) crossover (Eq. 12).","tokens_in":24946,"tokens_out":3303,"duration_ms":32404,"significance":"If the central asymptotic claim holds, this is a genuinely new organizing principle for hyperuniformity in active matter: it shifts attention from the spectrum of a prescribed field to the response-selected sector of a signed source, and it makes a falsifiable prediction (the universal k^4→k^6 crossover and its defect-controlled leakage) that is unusual in this literature. The paper's strengths are substantial: the cluster-expansion derivation of Eq. (6) is clean; the simulations use independent seeds with confidence intervals, exact Fourier summation, finite-size checks over N=256–1024 (Model I) and up to 4096 (Model II); and the local residual ρ_loc^u is an independent, physically motivated observable that does correlate with the fitted plateau. The paper is also unusually candid about its own limitations, explicitly flagging in SM S5.6 that AICc cannot always distinguish the crossover form from a free power law and in SM S8 that finite cluster size and finite-range correlations are assumed. The main weakness is that the numerical verification of the asymptotic k→0 normal form is less deep than the abstract implies, as detailed below.","major_comments":[{"comment":"The central law χ̃_T^f(k)=Ŵ^2(Δu k^4+B k^6) rests on S1(k)=Δu+Bk^2+O(k^4) holding all the way to k→0. The numerical support does not establish this asymptotic statement: all fits use k≤0.9, and at N=1024 the smallest resolved wave number is k_min≈0.098, only about 0.18 k× at k_off=0.06. The fit window is therefore dominated by crossover behavior, not by a clean k^4 plateau or a clean k^6 asymptote. SM S5.6 explicitly admits that a free power law can fit individual finite-window spectra as well as the two-term crossover, so AICc does not select the normal form. Moreover, B is extracted from the same S1(k) that the law is tested against, so the collapse in Fig. 2(d) is a consistency check rather than an independent prediction. A nonanalytic correction (e.g., c k^2 ln k or k^(2+ε)) would be absorbed into B/Δu over this window while changing the true k→0 limit and the predicted R^(−4) veloci","section":"Eq. (7), SM S4.1, SM S5.6"},{"comment":"The collapse of S1/(Bk^2) onto 1+(k×/k)^2 is presented as evidence for the universal crossover form. But k× and B are fitted from each individual spectrum using the very two-term form being tested, so the collapse is guaranteed asymptotically by construction for any S1 that is smooth over the fitted window. The independent anchor is the Δu–ρ_loc^u proportionality, which validates the plateau amplitude as a local defect density, but nothing independently measures B or the k^2 curvature. I recommend stating more explicitly that Fig. 2(d) is a parameterization check, not a falsifiable prediction, and moving the main evidentiary weight to the velocity prediction and the local–spectral proportionality, which are genuinely independent.","section":"Fig. 2(d), Eq. (9)"},{"comment":"The local diagnostic ρ_loc^u includes only unpaired carriers and explicitly excludes the residual of imperfect bound pairs. The claim that unpaired carriers dominate Δu rests on the observed near-linear Δu–ρ_loc^u relation, but this is not a direct measurement of the bound-pair residual. If imperfect bound-pair residuals were proportional to ρ_loc^u, the near-linear relation could hold even if those residuals contributed a significant part of the plateau. The fixed-partner model's healing branch also maintains a small but nonzero plateau (Δu≈1.0–1.3×10^(−4), Table S7), so the 'strictly hyperuniform velocity limit' described in the abstract is never actually realized in simulation; it is a theoretical limiting statement. I would not require the authors to realize Δu=0, but the text should be careful not to present the zero-residual limit as a demonstrated simulation result rather than a c","section":"Eq. (10), SM S3.3, SM S7.4"}],"minor_comments":[{"comment":"There are several typographical artifacts: 'universallaweχf' and 'response-selectedhyperuniformity' lack spacing; 'eχf' and other calligraphic symbols use inconsistent fonts. These should be cleaned before publication.","section":"Abstract/Introduction"},{"comment":"The caption says 'shaded regions are 95% confidence intervals' but the panels in Fig. 1(e) show bands; 'shaded regions' should be 'shaded bands' for consistency with the text. Also, the gray k^4 and k^6 slope guides are useful, but it would help to state explicitly over which k range each guide is drawn.","section":"Fig. 1"},{"comment":"The parameter 'reaction stride n_r' is listed but not defined in the table caption; define it as the number of integration steps between reaction attempts for completeness.","section":"SM S5.1 / Table S1"},{"comment":"The statement that the real-space fit length ξ_v is smaller than ξ_scr is clear, but the discussion would benefit from a sentence explaining why the Gaussian window finite-k effects push ξ_v below ξ_scr; currently the reader must infer this from the figure.","section":"SM S5.7"},{"comment":"Some references appear to be future-dated (e.g., [25] '2026', [27] '2026', [11] '2026'). If these are preprints or accepted articles, please add the arXiv or DOI identifiers; if they are placeholders, update them before publication.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed paper with a genuinely new conceptual contribution and unusually honest supplementary material. The central issue is not correctness of the finite-window data but the gap between the asymptotic universal claim (Eq. 7 as k→0) and the numerical evidence, which is confined to a crossover-dominated window. The authors themselves acknowledge this in SM S5.6 and SM S8, which is a credit to them, but the main text's abstract overstates the support. With a reasonable expansion of the infrared analysis (or a modest reformulation of the claim as a finite-window crossover law), the paper could be acceptable. I recommend major revision, not rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a real conceptual advance, not just another hyperuniformity example. The idea that \"hyperuniform\" should be predicated on a source–response pair—since incompressible flow only sees transverse forces—is clean and likely to stick. The two simulation models and the independent local residual measurement give it real weight. But the paper's marquee quantitative claim, the universal k^4-to-k^6 crossover, is fit over a window that straddles the crossover, so the asymptotic purity of that law is not yet established.\n\nWhat's new: response-selected hyperuniformity, exchangeable multipole inheritance, defect-controlled infrared leakage. The cluster expansion leading to S1 = Δu + B k^2 is standard but applied to the right object (transverse selected first moments), and the k^4 from the two-derivative force map is transparent. The valence-one model with no prescribed partners shows partner memory decays yet order persists—that's a nice result. The local unpaired-moment density ρ_loc^u predicting Δu, with a slope near 1 across 34 runs, is a genuinely independent anchor. The velocity variance prediction matching the spectral integral to a few percent is also good.\n\nSoft spots, in order. First, the data only resolve k_min ~ 0.1, and k× ~ 0.55 at weak turnover, so the \"plateau\" and \"k^6\" asymptotes are never actually seen cleanly; the fits are dominated by the crossover. Their own SM S5.6 says AICc can't always distinguish the two-term form from a free power law. Second, B is extracted from the same S1(k) the law is tested against, so the collapse is a consistency check, not a prediction. Third, the \"hidden\" claim would be stronger if they had directly shown S_ρ(k) for the number density is non-hyperuniform; they assert it but don't display it. That's a small fix. Fourth, no code or data repository is provided, which matters for a paper whose quantitative claims rest on fitted parameters.\n\nI don't think any of these are fatal. The analyticity argument is solid for finite-range clusters; the infrared leakage only matters if you believe the k^2 term is truly the leading correction. The authors are careful about seed-level statistics and explicitly disclaim fractional exponents as crossovers. I'd send this to a serious referee, with the request to push for deeper infrared simulations and public data.\n\nIt will be a useful paper for the active-matter and hyperuniformity communities, and I'd cite it as a conceptual framework, not as a measured asymptotic law. My recommendation: accept conditional on that request.","headline":"Genuine conceptual advance—hyperuniformity defined relative to a response operator—with a sound analytic core and two modeling checks, but the marquee k^4-to-k^6 crossover is fit over a crossover-dominated window, so the asymptotic purity of the law is not yet proven.","tokens_in":25430,"tokens_out":2465,"would_cite":true,"duration_ms":25517,"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 paper claims that hyperuniformity in active fluids comes from locally balanced active forces, and that rare unpaired moments set the finite screening length of quiet flow.","keywords":["hyperuniformity","active matter","response-selected order","multipole screening","transverse force spectrum","valence-one fluid","infrared crossover","screening length"],"falsifier":"Measure the transverse-force spectrum at wavenumbers well below k× = sqrt(Δu/B) in a steady active fluid with a measurable finite density of unpaired moments: if chi^f_T continues as k^6 rather than bending to k^4 in that deepest infrared window, the normal form S1 = Δu + B k^2 fails. Concretely, the paper's collapse S1/(B k^2) = 1 + (k×/k)^2 and the predicted real-space velocity variance V_v(R) ≃ a Δu R^{-2} + b B R^{-4} are testable; a free-power-law fit that beats the two-term form on the smallest resolved shells and does not bend toward k^4 would contradict the claim.","tokens_in":24422,"feed_emoji":"🌊","tokens_out":11262,"duration_ms":89247,"temperature":0.7,"pith_summary":"Hyperuniformity — anomalously suppressed fluctuations at long wavelengths — in an active fluid is not determined by how uniformly the particles are arranged, but by which part of their microscopic forcing the surrounding fluid can actually respond to: in an incompressible fluid, longitudinal forces are absorbed into pressure, so only the transverse (solenoidal) part of the signed active-moment field drives flow. The authors define response-selected hyperuniformity as suppressed long-wavelength fluctuations in that transverse source sector, and derive the universal law chi^f_T(k) = W(k)^2 (Δu k^4 + B k^6) for the transverse-force spectrum, with a crossover wavenumber k× = sqrt(Δu/B). They verify it in two simulated active fluids — a valence-one carrier fluid in which partners constantly exchange, and a fixed-partner molecular model — and find that complete partner renewal leaves the k^6 window intact while the unpaired-moment density tunes the k^4 residual. A sympathetic reader should care because this reframes 'hyperuniform active matter' from a structural property into a property of a source–response pair, and identifies the rare unscreened moments that set the finite screening length of quiet flow.","feed_headline":"Partner exchange doesn't kill quiet flow in active fluids","feed_subtitle":"A finite density of unpaired active moments fixes the screening length of velocity fluctuations.","key_machinery":"The central object is the response-selected first-moment spectrum S1(k): the spectrum of the signed active-moment field after contraction with the transverse projector P^T_ij = δ_ij − k̂_i k̂_j, which selects the sector of active forcing that can actually drive incompressible flow. The identity carrying the argument is the cluster expansion: a homogeneous population of finite, locally neutral clusters contributes B k^2 + O(k^4) to S1 regardless of partner identity, while unpaired or imperfectly neutral carriers contribute a plateau Δu. The source-to-force map (divergence of the active stress followed by transverse projection) adds four powers of k, yielding chi^f_T = W^2 k^4 S1 = W^2 (Δu k^4","core_discovery":"Central claim: the response-selected first-moment spectrum obeys S1(k) = Δu + B k^2 + O(k^4), giving the transverse-force spectrum chi^f_T(k) = W(k)^2 (Δu k^4 + B k^6). The k^6 term comes from finite locally neutral clusters of opposite-signed active moments whose leading multipole cancels; the k^4 term comes from a nonzero infrared plateau Δu due to unscreened or unpaired moments. Since the Stokes Green function adds an extra k^{-2} per vector component, the velocity spectrum is ≃ (Δu + B k^2)/η^2 — strictly hyperuniform only when Δu = 0, with any finite residual restoring a velocity plateau and setting a screening length ξ_scr = sqrt(B/Δu). Two simulated active fluids (exchangeable valence","pith_inferences":["Editorial inference — the framework suggests a practical design rule for quiet active fluids: minimize the local density of unpaired selected moments, since the plateau Δu (not the bound fraction) controls the screening length.","Editorial inference — in chiral, odd-elastic, or compressible fluids, the response operator can mix longitudinal and transverse sectors, creating additional channels through which hidden residuals leak into observable motion; the paper notes this possibility but does not analyze it.","Editorial inference — the fixed-partner model's inability to select the complementary branch from every random preparation, while the exchangeable model assembles it de novo, hints that partner-exchange kinetics may be important for realizing response-selected hyperuniformity from generic initial conditions."],"forward_implications":["If Δu = 0, the fluid's velocity fluctuations are strictly hyperuniform even when the particle density is not; the order lives in the hidden transverse-moment sector, invisible to ordinary structure factors.","Any finite density of locally unscreened moments is an infrared-relevant perturbation: at sufficiently long wavelengths the k^4 term always dominates, so strict velocity hyperuniformity is lost and the quiet-flow regime has a finite range ξ_scr = sqrt(B/Δu).","Constant partner exchange does not destroy the hidden order — the k^6 window persists under complete partner renewal — because the order is carried by instantaneous local neutrality rather than permanent molecular identity.","The same universal normal form appears in two different models with different reaction chemistry and defect routes, so the leading powers and crossover are fixed by local neutrality, analyticity, and finite correlation length, not by a particular binding scheme.","Non-integer spectral exponents measured over finite wavenumber windows (for example ≈ 4.8) are the crossover signature predicted by the two-term normal form, not new asymptotic exponents."],"fun_headline_variants":["Unpaired moments set quiet-flow length in active fluids","Active fluids keep quiet flow unless moments are unpaired","Zero unpaired moments gives hyperuniform active flow","Response-selected order: why some active flows stay calm","Screening length fixed by rare unscreened active moments"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the long-wavelength first-moment spectrum has the exact analytic form Δu + B k^2 + O(k^4) all the way to k → 0; the paper's Supplemental Material (Sec. S5.6) notes the two-term fit is made only over k ≤ 0.9 and a free power law fits individual spectra as well, so the universal k^4/k^6 separation and ξ_scr = sqrt(B/Δu) presuppose that this analytic expansion remains valid in the deepest infrared.","fun_headline_variants_meta":{"raw":{"variants":["Unpaired moments set quiet-flow length in active fluids","Active fluids keep quiet flow unless moments are unpaired","Zero unpaired moments gives hyperuniform active flow","Response-selected order: why some active flows stay calm","Screening length fixed by rare unscreened active moments"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000605,"raw_usage":{"total_tokens":2698,"prompt_tokens":823,"completion_tokens":1875,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":1798}},"tokens_in":567,"tokens_out":1875,"duration_ms":11344,"temperature":1.0,"reasoning_tokens":1798,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T23:03:15.154018+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the transverse-force spectrum at wavenumbers well below k× = sqrt(Δu/B) in a steady active fluid with a measurable finite density of unpaired moments: if chi^f_T continues as k^6 rather than bending to k^4 in that deepest infrared window, the normal form S1 = Δu + B k^2 fails. Concretely, the paper's collapse S1/(B k^2) = 1 + (k×/k)^2 and the predicted real-space velocity variance V_v(R) ≃ a Δu R^{-2} + b B R^{-4} are testable; a free-power-law fit that beats the two-term form on the smallest resolved shells and does not bend toward k^4 would contradict the claim.","supporting_citations":[],"review_version":1}