{"id":"2140a874-b2e8-451a-9036-e9ac0550efff","arxiv_id":"2606.17595","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In active nematic turbulence, all ±1/2 nematic and principal-stress defects localize on the zero contour of the in-plane deviatoric stress invariant, independent of activity sign and strength.","lead":"Simulations of active nematics show that topological defects sit on the zero contour of a simple stress anisotropy measure, and that principal stress lines up with or against cell orientation depending on whether activity is extensile or contractile. That gives experimentalists a stress-only way to locate defects in cell monolayers without needing detailed cell-shape maps.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the reader's already-flagged active-stress dominance premise.","rationale":"The central claim splits cleanly into a definitional part (principal-stress defects live at J2=0) and an emergent part (nematic defects share those loci). The definitional part needs no further scrutiny. The emergent part rests exactly on the active-stress dominance premise identified by the reader; the manuscript supplies only visual evidence and an appeal to the constitutive law. No additional internal inconsistency, missing control, or hidden assumption of comparable weight appears. Parameter count, lack of code, and visual rather than statistical quantification are already noted by the reader and do not alter the logical structure. Consequently the CONDITIONAL verdict stands: the observation is useful and publishable once quantitative coincidence metrics and a clearer separation of definitional versus emergent content are added. No change to the reader's assessment is required.","tokens_in":11656,"tokens_out":596,"duration_ms":5362,"concrete_test":"Extract the full stress tensor on a dense set of frames in the turbulent regime (ζ = ±0.2, ±0.5, ±0.8). Compute the spatial offset distribution between every nematic ±1/2 core (S minimum) and the nearest J2 minimum / (σ_xx−σ_yy=0, σ_xy=0) point; report mean offset in units of the microscopic length L0=√(L/C). Simultaneously evaluate the local ratio |σ^a| / |σ^p+σ^v| inside a 3L0 radius of each core. If mean offset exceeds ~L0 or the active fraction falls below ~0.7, the coincidence claim weakens.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest_assumption correctly isolates the only non-definitional hinge: that zeros of nematic order and of total-stress anisotropy coincide because active stress σ^a_ij = −ζ Q_ij dominates the anisotropic part of the total hydrodynamic stress. Principal-stress defects sit at J2 minima by construction (eigenvectors undefined when the tensor is isotropic). Nematic coincidence is therefore the sole substantive claim. The paper itself states (Sec. III.D) that total stress is not strictly proportional to Q, yet asserts coincidence “because of the dominant role of active stress.” No quantitative measure of that dominance (e.g., relative magnitude of |σ^a| versus |σ^p + σ^v| near defect cores, or measured spatial offset between S=0 and J2=0 loci) is supplied; support remains visual (Figs. 4, 7). If passive or viscous contributions systematically displace the zeros, the claimed localization of nematic defects on the J2 isoline fails. That is precisely the concern already raised by the reader; I find no deeper or independent load-bearing flaw.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript reports hybrid Lattice-Boltzmann simulations of the standard continuum active-nematic equations (Beris–Edwards hydrodynamics with active stress −ζQij) in two dimensions. It shows that the orientation of maximal principal stress aligns perpendicular (parallel) to the nematic director for extensile (contractile) activity, quantified by angle PDFs and cross-correlations. In the turbulent regime the authors introduce the rotation-invariant scalar J2 (second invariant of the in-plane deviatoric stress) and claim that its zero-level isolines (equivalently the loci of vanishing anisotropy) coincide with the cores of all ±1/2 topological defects of both the nematic director and the principal-stress orientation field; the coincidence is stated to be robust to the magnitude and sign of activity. The results are motivated by the experimental accessibility of stress fields in cellular monolayers via traction-force microscopy.","tokens_in":11909,"tokens_out":1063,"duration_ms":19153,"significance":"If the claimed localization of nematic defects on the J2 isoline is quantitatively robust, the work supplies a practical, stress-only criterion for locating topological defects—an experimentally useful route when cell-shape or order-parameter fields are harder to extract than principal stresses. The alignment result follows directly from the structure of active stress and is consistent with existing observations. The simulations themselves are standard and the parameter set is fully stated, so the findings are in principle reproducible. The advance is incremental rather than transformative: it re-examines known active-nematic turbulence through the lens of principal stresses and adds a geometric observation whose practical value depends on the strength of the active-stress dominance assumption.","major_comments":[{"comment":"Sec. III.D and Figs. 4, 7: Localization of principal-stress defects at J2 minima is definitional (principal directions become undefined once the stress tensor is isotropic). The only non-trivial claim is therefore the coincidence of nematic ±1/2 defects with the same loci. The text acknowledges that total stress is not strictly proportional to Q yet asserts coincidence “because of the dominant role of active stress.” No quantitative diagnostic is supplied—neither the relative magnitude |σa| / |σp + σv| near defect cores nor a histogram (or mean/variance) of spatial offsets between S = 0 and J2 = 0 points across activity strengths and signs. Visual overlay of a few snapshots is insufficient to establish the robustness asserted in the abstract and conclusion; without such a measure the central claim remains untested against the possibility that passive or viscous contributions systematical","section":"Sec. III.D, Figs. 4 and 7"},{"comment":"Sec. III.D (and Appendix B): The argument that (q1, p1) constructed from the anisotropic stress components “closely resemble” the nematic directors is again only visual (Fig. 8). A direct comparison of the two order-parameter fields (e.g., spatial correlation of their magnitudes or of their defect-core positions as a function of ζ) is needed to quantify how completely active stress dominates the anisotropy that sets the defect locations.","section":"Sec. III.D, Fig. 8"}],"minor_comments":[{"comment":"Title inconsistency: the arXiv title and abstract use “Vanishing Deviatoric Stress,” while the manuscript body title uses “Stress Anisotropy.” Choose one and keep it consistent.","section":"Title / Abstract"},{"comment":"Abstract grammar: “whose zero-level contour coincides with the locations of all ±1/2 topological defects \to are localized” is incomplete; rewrite for clarity.","section":"Abstract"},{"comment":"Fig. 1 caption and text: the angle difference Δθ is reported in radians, yet the peaks are described as π/2 and 0; a brief note that the director is headless (so Δθ ∈ [0, π/2]) would avoid confusion.","section":"Fig. 1"},{"comment":"Notation: the maximal-principal-stress director is variously called ns, np and n_s; unify the symbol.","section":"Sec. III.A–B"},{"comment":"Appendix A, Fig. 5: the time axis is labeled t/100 without stating the unit; clarify lattice-Boltzmann time units.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The load-bearing scientific issue is narrow and fixable with additional post-processing of existing simulation data; the rest of the manuscript is solid. Novelty relative to recent experimental papers on stress defects in cell monolayers (e.g., Nejad et al. 2024) is modest, so the quantitative strengthening requested above is essential for a high-impact soft-matter journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing worth knowing is that in standard active-nematic turbulence both the nematic ±1/2 defects and the principal-stress defects sit on the zero contour of the rotation-invariant deviatoric stress J2 (and on σxx−σyy=0). That localization is robust to the sign and magnitude of ζ. The alignment result—principal stress perpendicular to the director for extensile activity, parallel for contractile—is almost immediate from σa=−ζQ, but the paper documents it cleanly with angle PDFs and cross-correlations.\n\nWhat is actually new is the concrete numerical claim that the zeros coincide for both fields across the turbulent regime. The experimental motivation is real: traction-force and monolayer-stress maps are often easier to obtain than full cell-shape fields, so a stress-based defect locator has practical value. The simulations are the usual hybrid Lattice-Boltzmann Beris–Edwards setup with stated parameters; the figures (especially 4 and 7) make the visual case, and the authors correctly note that principal-stress defects sit at J2 minima by construction because the eigenvectors become undefined.\n\nThe soft spot is exactly the one the reader flagged. Nematic coincidence is not definitional; it leans on active stress dominating the anisotropic part of the total stress so that zeros of S and of J2 line up. The paper itself says total stress is not strictly proportional to Q, yet offers only visual overlays—no relative-magnitude check near cores, no measured spatial offset distribution. That is a genuine but limited gap, not a collapse of the argument. There are also a couple of formula typos in the J2 expression in the figure captions, and no code or quantitative coincidence metric is supplied. None of these are fatal for a soft-matter simulation paper.\n\nThis is for people who already work with active nematics or with stress microscopy in cell monolayers. It will not rewrite the foundations, but it gives a usable mechanical diagnostic. I would send it to peer review; a referee can simply ask for a quantitative offset histogram and a short dominance check. Worth citing if you are writing about stress–defect relations or TFM interpretation.","headline":"Clean numerical observation that both nematic and principal-stress defects sit on the vanishing contour of deviatoric stress J2, useful for TFM experiments, with the main non-definitional claim resting on active-stress dominance that is only shown visually.","tokens_in":12551,"tokens_out":551,"would_cite":true,"duration_ms":5022,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"In active nematic turbulence, all half-integer defects sit where the in-plane deviatoric stress vanishes.","keywords":["active nematics","topological defects","deviatoric stress","principal stress","active turbulence","cell monolayers","defect localization"],"falsifier":"Measure both the nematic director field and the full stress tensor in a turbulent active nematic (simulation or cell monolayer) and check whether any half-integer defect lies off the J2 = 0 contour once passive or viscous stresses are made comparable to the active stress.","tokens_in":12491,"feed_emoji":"🔬","tokens_out":907,"duration_ms":9493,"temperature":0.7,"pith_summary":"The paper studies how stresses organize themselves in a standard continuum model of active nematics, the theoretical description often used for dense cell monolayers. It shows that the direction of maximal principal stress lines up perpendicular to the cell orientation for extensile activity and parallel for contractile activity. In the chaotic turbulent regime the authors isolate a simple scalar built from the anisotropic part of the stress tensor. Its zero contour is an organizing backbone: every topological defect of charge plus or minus one-half, whether defined from the orientation field or from the stress field itself, sits exactly on that contour. The localization holds for both signs of activity and across a range of activity strengths. Because traction-force and monolayer-stress microscopy already measure stress more easily than cell shape, the result offers a practical mechanical marker for defect locations in living tissues.","feed_headline":"Defects in active fluids sit where stress anisotropy vanishes","feed_subtitle":"A simple stress scalar marks every half-integer defect, for both extensile and contractile activity","key_machinery":"The second invariant of the deviatoric stress, J2 = sqrt(((σxx−σyy)/2)^2 + σxy^2), whose vanishing requires the simultaneous conditions σxx−σyy = 0 and σxy = 0 and thereby forces the principal directions to become undefined; the isoline σxx−σyy = 0 therefore supplies a geometric backbone on which all half-integer defects must sit.","core_discovery":"In the fully developed turbulent state of an incompressible active nematic, a rotation-invariant scalar measure of the in-plane deviatoric stress (the second invariant J2, equivalently half the difference of the principal stresses) reaches zero along a continuous isoline; every plus-or-minus one-half topological defect of both the nematic director and the maximal-principal-stress director is localized on that isoline. The coincidence is independent of the magnitude and of the sign (extensile or contractile) of the activity parameter.","pith_inferences":["If the same J2 isoline continues to host defects when substrate friction or multi-layer geometry is added, the localization rule may serve as a general diagnostic for active turbulence beyond the ideal continuum model.","Experimental groups that already extract principal-stress maps from monolayer stress microscopy could test the prediction by simply contouring J2 and overlaying independently identified nematic defects.","The offset between the isotropic-stress maximum and the core of a stress defect (versus a nematic defect) suggests that flow-induced stress singularities may be used to distinguish passive from active contributions in mixed tissues."],"forward_implications":["Principal-stress defects of charge ±1/2 exist and share the same locations as the ordinary nematic defects.","The isoline of vanishing normal-stress difference can be used as a mechanical proxy to locate all half-integer defects without reconstructing the full orientation field.","Alignment of maximal principal stress with the director flips from perpendicular (extensile) to parallel (contractile), giving a local diagnostic of activity type.","Stress measurements already available from traction-force microscopy become sufficient to map defect topology in confluent monolayers."],"fun_headline_variants":["Active-nematic half-defects sit on zero isolines of deviatoric stress","Turbulent ±1/2 defects localize where in-plane stress anisotropy vanishes","J2 = 0 contours pin every nematic and principal-stress half-defect","Vanishing stress anisotropy marks all ±1/2 defects independent of activity","Rotation-invariant stress scalar zeros locate topological defects"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The claim rests on the active stress being strong enough that the zeros of stress anisotropy still coincide with the zeros of nematic order even though the total stress also contains passive and viscous contributions.","fun_headline_variants_meta":{"raw":{"variants":["Active-nematic half-defects sit on zero isolines of deviatoric stress","Turbulent ±1/2 defects localize where in-plane stress anisotropy vanishes","J2 = 0 contours pin every nematic and principal-stress half-defect","Vanishing stress anisotropy marks all ±1/2 defects independent of activity","Rotation-invariant stress scalar zeros locate topological defects"]},"model":"grok-4.5","effort":"low","cost_usd":0.006734,"raw_usage":{"total_tokens":1689,"prompt_tokens":754,"num_sources_used":0,"completion_tokens":101,"cost_in_usd_ticks":67340000,"prompt_tokens_details":{"text_tokens":754,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":834,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":754,"tokens_out":101,"duration_ms":6682,"temperature":1.0,"reasoning_tokens":834,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T13:37:53.321977+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure both the nematic director field and the full stress tensor in a turbulent active nematic (simulation or cell monolayer) and check whether any half-integer defect lies off the J2 = 0 contour once passive or viscous stresses are made comparable to the active stress.","supporting_citations":[],"review_version":2}