{"id":"6e37aac0-efac-4097-94f3-596726e2fc21","arxiv_id":"2607.24927","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Compressibility sets density contrast under patterned activity via one dimensionless parameter and thereby steers active turbulence, including a soft-confined one-dimensional vortex chain.","lead":"A continuum theory shows that compressibility alone can relocate active turbulence and pin a controllable chain of vortices at activity interfaces in nematic suspensions. That gives experimentalists a single knob—how easily density can vary—to pattern flows without hard walls.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The vortex-chain \"analytically tractable steady state\" is a one-way kinematic construction: a texture is assumed, flow is solved from it, but stationarity of Q under the computed flow and the q ~ 1/W selection rule are never verified. The Δϕ = R/√(1+R) and turbulence-relocation results are unaffecte","rationale":"The reader identified the two-fluid one-fluid reduction (SI Sec. S1) as the weakest assumption. That is a real modeling limitation but, as the reader themselves noted, it bounds experimental mapping rather than internal correctness — the claim is explicitly developed within this model class, and per good-faith reading I do not treat \"the model may not match future experiments\" as load-bearing for a theory paper. I therefore only partially agree: I locate the load-bearing soft spot elsewhere, in the internal support for the most novel sub-claim (the vortex chain), where the analytics are kinematic rather than self-consistent and the wavenumber selection is heuristic. I keep the verdict at ACCEPT because (i) the central quantitative result, Δϕ = R/√(1+R), is an exact consequence of the scalar reduction, I reproduced its algebra, and it matches parameter-free numerics across κ ∈ [1,100] with the wavelength window of validity honestly delineated in Fig. S5; (ii) the turbulence-relocation mechanism is backed by both an effective-activity argument and simulations; (iii) the vortex chain itself is independently established by the full numerical solution (Fig. 3c), with public code — my concern is about the strength of the analytical characterization, not the existence of the state. If the proposed residual/persistence test failed badly (residual O(1), m coarsening), the appropriate move would be CONDITIONAL with the chain re-framed as a long-lived transient, but there is no current evidence for that, and the paper's own language (\"approximate analytical solution\") is not overreaching enough to force a downgrade now.","tokens_in":32145,"tokens_out":7368,"duration_ms":91870,"concrete_test":"Two-part check. (1) Self-consistency residual: insert the ansatz (Q from Eqs. S138/S139/S144, ϕ from Eq. S133, v from Eqs. S154–S156 with K₁, K₂ from S170/S173) into the RHS of Eq. (3) and compute ‖λ_N A − [ω,Q] + H/γ − ∇·(vQ)‖ over the interfacial region; if this residual is comparable to the individual terms rather than small at O(θ₀), the texture is not a steady solution and the \"steady state\" claim needs re-derivation or re-labeling. (2) Dynamical persistence: initialize the full solver with the analytical ansatz at W = L/64 and W = L/32, run ≥10× the reported t_tot, and track vortex number m(t); if m drifts/coarsens or fails to scale as L/W, the selection rule q ~ 1/W and the steady-state interpretation fail even though a transient chain exists.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's strongest claim has three parts. Parts (a) and (b) — the density-contrast law and the relocation of turbulence — check out. I verified the scalar-model algebra (Eqs. S110–S127): the zero-flux solution ϕ = C/(1+Rf), the sinusoidal normalization C = √(1+R), and Δϕ = R/√(1+R) are internally correct, the dropped viscous term η∂²_xv introduces corrections of order κDηq² (~6% at κ=20 for λ=L/2, rising at κ=100), and the agreement with 2D numerics is parameter-free. Part (c), however, is the soft spot. The interfacial vortex-chain solution (Sec. S7) proceeds by (i) assuming a fixed order-parameter profile S(x) = S₀/2 in the interfacial region, (ii) assuming a bend texture θ(y) = θ₀ sin(qy) with free amplitude θ₀, (iii) taking ϕ(x) from the scalar model, and (iv) solving only the Stokes equation (S141) for the resulting vorticity. Nowhere is it checked that this texture is a fixed point of the full dynamics: inserting the ansatz into Eq. (3), the advection ∇·(vQ), flow-alignment λ_N A, and co-rotation [ω,Q] terms generated by the computed v generically act on the assumed director field at leading order in θ₀, so ∂_tQ ≠ 0 unless a cancellation holds that is never demonstrated. The wavenumber selection q ~ 2π/W (hence m ~ L/W vortices) rests on balancing Frank energies of longitudinal vs. transverse modulations with the extra assumption A_∥ ~ A_⊥ (Eq. S140) — a heuristic imported from hard-channel confinement [69], not a stability or dispersion calculation. Additional approximations sit exactly where they hurt most: the friction is linearized as ζϕ ≈ ζ across the interface even though ϕ swings from ~0 to ϕ* ≈ 2(1+R)/(2+R) there at high compressibility, and α_B(x) = α(x) is imposed although the main text uses α_B⁰ = (2/3)α⁰. None of this shows the chain does not exist — Fig. 3(c) numerics do show it — but the claim of an \"analytically tractable dynamical steady state\" currently means \"analytically computed flow around an assumed texture,\" and if the texture is not","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript develops an effective one-fluid hydrodynamic theory of a compressible active nematic, derived in the SI from a two-fluid quasi-2D suspension model in which only the active species compresses while the solvent absorbs out-of-plane flux. Three main results are presented. (1) Under spatially patterned activity, a scalar 1D reduction (Q=0) of the continuity plus Brinkman-Stokes equations yields a zero-flux steady state ϕ(x) ∝ 1/(1+Rf(x)) controlled by a single dimensionless ratio R ≃ α_B⁰κ; for a sinusoidal profile the density contrast is Δϕ = R/√(1+R), and this matches full 2D numerics with no fitting in the large-wavelength regime. (2) Tuning κ at fixed activity relocates active turbulence from high- to low-activity regions, rationalized by an effective activity |α(x)|ϕ(x) crossing the isotropic and nematic thresholds α_I^c, α_N^c from standard linear stability. (3) At sharp activity interfaces the system stabilizes a one-dimensional chain of counter-rotating vortices, for which Sec. S7 gives an approximate Green's-function solution of the Stokes equation with an assumed bend texture, and a Frank-energy argument estimating q ~ 2π/W. Linear stability, propagating density–nematic waves, and number-fluctuation crossovers (giant vs. normal) are worked out carefully in the SI.","tokens_in":32763,"tokens_out":3129,"duration_ms":57582,"significance":"If the results hold, the paper makes a substantial and timely contribution: it provides a parameter-free, closed-form prediction Δϕ = R/√(1+R) for the density contrast under patterned activity that agrees with full 2D numerics without fitting; it demonstrates a concrete control knob (compressibility) for relocating active turbulence, directly relevant to light-patterned microtubule–kinesin experiments; and it identifies a soft-confinement mechanism for one-dimensional vortex chains in a free-standing fluid. The work ships reproducible CUDA-based simulation code (Ref. [76]), the linear stability analyses are complete and consistent (transverse instability shown to preempt longitudinal ones), and the number-fluctuation crossover between dry and wet limits is cleanly derived. The falsifiable scaling m ~ L/W and the R-controlled density contrast give experimentalists direct tests.","major_comments":[{"comment":"Sec. S7, Eqs. (S138)–(S157): the vortex-chain solution is a one-way kinematic construction. A texture is assumed (fixed S(x)=S0/2, bend ansatz θ(y)=θ₀ sin(qy), ϕ from the scalar model), and only the Stokes equation (S141) is solved. Stationarity of Q under the computed flow is never checked: inserting the ansatz into Eq. (3), the advection ∇·(vQ), co-rotation [ω,Q], and flow-alignment λ_N A terms generated by the computed v act on the assumed director field at O(θ₀), so ∂_tQ ≠ 0 unless a cancellation holds that is not demonstrated. Since the abstract advertises an 'analytically tractable dynamical steady state,' a consistency check is needed — e.g., evaluating the residual of Eq. (3) on the solution and showing it is higher order, or demonstrating numerically that the texture maintains the sinusoidal form and wavenumber over long times.","section":"Sec. S7, Eqs. (S138)–(S157)"},{"comment":"Eq. (S140) and main text: the wavenumber selection q ~ 2π/W (hence m ~ L/W vortices) rests on balancing Frank energies of longitudinal vs. transverse modulations with the additional assumption A_∥ ~ A_⊥, imported from hard-wall channel confinement [69]. This is not a stability or dispersion calculation, and the paper tests only a single width (W = L/64 = 8). Since the scaling m ~ L/W is a falsifiable prediction, the authors should test it by varying W in the numerics and measuring the vortex count, or explicitly qualify the selection rule as a heuristic estimate with an uncontrolled amplitude assumption.","section":"Eq. (S140); Fig. 3(c)"},{"comment":"Sec. S7: the quantitative reach of the analytical vortex-chain solution is limited by two uncontrolled inputs that should be stated and, if possible, bounded. (i) The bend amplitude θ₀ is free, so the vorticity magnitude in Fig. 3(d) cannot be compared to Fig. 3(c) without a prescription for θ₀; the caption should state how θ₀ and q were chosen for the plotted analytic solution. (ii) The friction is linearized as ζϕ ≈ ζ across the interface (above Eq. (S141)), yet ϕ varies by O(1) there (Fig. S7); and the calculation assumes α_B(x) = α(x) (below Eq. (S146)) while the main-text numerics use α_B = (2/3)α₀. These approximations are reasonable for a first estimate, but their effect on μ_q and the vortex extent should be discussed.","section":"Sec. S7; Fig. 3(c,d)"}],"minor_comments":[{"comment":"The dropped viscous term η∂²_xv in the scalar reduction (Eq. (S108)) introduces corrections of order κDηq² relative to the retained terms; at κ=20 with λ=L/2 this is a few percent, but it grows with κ toward the κ=100 end of Fig. 2(b). A brief error estimate or regime-of-validity statement in the main text would strengthen the parameter-free claim.","section":"Eq. (S108), main-text Fig. 2 discussion"},{"comment":"Formatting: the sentence on the interfacial vortex chain reads '[66, 67][68] Although...' — the citations run into the text and should be reformatted.","section":"Main text, p. 4"},{"comment":"Fig. 3(c) vs 3(d): the figure would benefit from a shared color scale or a line-cut comparison of ω(x,y) so the reader can judge the agreement beyond the qualitative vortex pattern.","section":"Fig. 3(c,d)"},{"comment":"Main text states α_B⁰ = (2/3)α₀ as a definition; a one-line motivation (e.g., from the microstructure or from matching to experiments) would help readers who wish to generalize the activity patterning.","section":"Main text, activity patterning paragraph"},{"comment":"The SI derives giant number fluctuations (∆N ~ N) only in the friction-dominated limit, while the pure-viscous simulations in Fig. S4 find approximately normal fluctuations after rescaling. The main-text summary of Sec. S4 correctly distinguishes the two regimes, but a pointer to Fig. S4's rescaled collapse would prevent confusion.","section":"Sec. S4; main text p. 3"}],"recommendation":"major_revision","confidential_remarks":"The experimental counterpart motivating the work (Ref. [36], Cheng et al., \"in preparation\") is co-authored by several of the present authors and is cited as the system the theory \"quantitatively reproduces,\" but no data from it are shown here; the editor may wish to confirm the two submissions are appropriately coordinated. The citation pattern leans heavily on the authors' own prior work, though this largely reflects their role in developing the field."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful core is simple: they write a wet active-nematic hydrodynamics that keeps density free, show that patterned isotropic activity sets the density contrast through one number R ≃ α_B κ, and then use κ to move active turbulence from high- to low-activity regions. That fills a real gap between dry compressible models and the usual incompressible wet ones, and it is aimed straight at the light-patterned MT–kinesin experiments with large density swings.\n\nWhat they do well is the scalar reduction and the steering story. The 1D advection–diffusion solution is clean, parameter-free against the 2D runs at long wavelength, and the I–II–III crossover with R is easy to read off the effective activity αϕ versus the isotropic and nematic thresholds. Linear stability, structure factors, and the SI are careful. Code is linked. Circularity is low: R and Δϕ come from continuity plus overdamped Stokes, not from fitting the turbulence maps.\n\nThe soft spot is the interfacial vortex chain. Numerics in Fig. 3(c) show a localized chain; that part is real. The “analytically tractable dynamical steady state” is oversold. They fix S and a bend texture, pull ϕ from the scalar model, and solve Stokes. They never close the Q equation under the resulting flow, and q ~ 1/W is a hard-channel energy balance with A_∥ ~ A_⊥, not a dispersion relation. Friction is linearized across a large density jump and α_B = α is imposed in the SI while the main text uses 2/3. So you get an analytic flow around an assumed texture, not a demonstrated fixed point. That does not touch Δϕ(R) or the relocation of turbulence.\n\nOther limits are ordinary modeling ones: uniform solvent density in the two-fluid reduction, ideal-gas pressure, companion experiment still in prep. None of those are load-bearing errors inside the model class.\n\nThis is for people who do active nematics, patterned activity, or density-coupled active flows. Worth a reading-group slot and a serious referee. I would cite the compressible framework and the R-controlled contrast; I would treat the vortex-chain analytics as a useful sketch backed by numerics, not as a closed solution. Send it out.","headline":"Solid compressible wet theory with a clean density-contrast law and real control of turbulence; the vortex-chain analytics are weaker than advertised but the numerics still carry that piece.","tokens_in":33087,"tokens_out":580,"would_cite":true,"duration_ms":16212,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Compressibility steers density and turbulence in active nematics, and can pin a one-dimensional vortex chain at activity interfaces.","keywords":["compressible active nematics","active turbulence","patterned activity","density contrast","vortex chain","soft confinement","active Péclet number","lyotropic liquid crystals"],"falsifier":"In a light-patterned quasi-2D active nematic, measure the steady density contrast versus compressibility (or activity strength) under a long-wavelength sinusoidal activity profile; the claim fails if the contrast does not track R/√(1+R) or if raising compressibility does not relocate turbulence into the low-activity region.","tokens_in":32641,"feed_emoji":"🌊","tokens_out":850,"duration_ms":15081,"temperature":0.7,"pith_summary":"Light-patterned active suspensions show large density variations that standard incompressible theories cannot capture. This paper builds a continuum theory of compressible active nematics in which density and flow evolve together. Under patterned activity, extensile isotropic pressure pushes material out of high-activity regions and into low-activity ones; a solvable one-dimensional reduction shows that the density contrast is set by a single dimensionless ratio linear in compressibility. Raising that ratio relocates active turbulence from high- to low-activity regions. At sharp activity interfaces the same mechanism confines flow to a one-dimensional chain of counter-rotating vortices held by soft, activity-induced confinement rather than solid walls. The result frames compressibility as an experimentally tunable knob for organizing both density and turbulent flow.","feed_headline":"Compressibility steers active turbulence and pins vortex chains","feed_subtitle":"One dimensionless ratio sets density contrast and can relocate flow from bright to dark regions","key_machinery":"The active compressibility ratio R ≃ α_B^0 κ (activity Péclet number comparing active transport to passive relaxation). In the exactly solvable scalar 1D reduction it alone sets the steady density profile ϕ(x) = C/(1+R f(x)) and the contrast Δϕ = R/√(1+R) for a sinusoid, thereby controlling where effective activity exceeds the turbulence threshold.","core_discovery":"Compressibility is a control parameter for density variations and turbulent flow in active nematic suspensions. Under spatially patterned activity the density contrast is governed by the single dimensionless ratio R proportional to activity times compressibility, and raising R steers active turbulence from high- to low-activity regions and, at sharp interfaces, stabilizes an analytically tractable one-dimensional vortex chain held by activity-induced soft confinement.","pith_inferences":["If solvent density is not uniform, a full two-fluid treatment may shift the R threshold at which turbulence relocates, so experiments that independently vary solvent compressibility would test the reduction.","The same soft-confinement mechanism may generate controllable defect or vortex tracks in other light-addressable active systems beyond nematics.","Mapping R across activity wavelengths shorter than the vorticity correlation length would chart where the scalar density theory breaks and nematic mixing dominates."],"forward_implications":["Compressibility can be used experimentally to move active turbulence between illuminated and dark regions without changing the activity pattern itself.","Sharp activity interfaces can replace hard walls as reconfigurable boundaries that localize a one-vortex-thick chain.","The interfacial vortex chain is in principle mobile and controllable, offering a route to transport suspended objects.","Density organization and flow localization in light-responsive microtubule–kinesin suspensions become quantitatively predictable from the single parameter R."],"fun_headline_variants":["Compressibility steers turbulence and pins vortex chains","Single ratio R sets density contrast and relocates active flow","Patterned activity plus compressibility locks 1D vortex chains","Raise compressibility to shift turbulence into low-activity zones","Activity-induced soft confinement stabilizes steady vortex chains"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The model treats the suspension as an effective one-fluid system in which the solvent density stays uniform and only the active component can compress, while both species share the same in-plane velocity.","fun_headline_variants_meta":{"raw":{"variants":["Compressibility steers turbulence and pins vortex chains","Single ratio R sets density contrast and relocates active flow","Patterned activity plus compressibility locks 1D vortex chains","Raise compressibility to shift turbulence into low-activity zones","Activity-induced soft confinement stabilizes steady vortex chains"]},"model":"grok-4.5","effort":"low","cost_usd":0.004387,"raw_usage":{"total_tokens":1267,"prompt_tokens":695,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":43868000,"prompt_tokens_details":{"text_tokens":695,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":491,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":695,"tokens_out":81,"duration_ms":8395,"temperature":1.0,"reasoning_tokens":491,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T05:37:38.659909+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"In a light-patterned quasi-2D active nematic, measure the steady density contrast versus compressibility (or activity strength) under a long-wavelength sinusoidal activity profile; the claim fails if the contrast does not track R/√(1+R) or if raising compressibility does not relocate turbulence into the low-activity region.","supporting_citations":[],"review_version":1}