{"id":"a33bd3f9-501d-4d87-b79c-be18bacd395b","arxiv_id":"2606.06496","paper_version":1,"verdict":"CONDITIONAL","confidence":"UNKNOWN","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":5,"one_line_summary":"A source-extended Lie-group algorithm derives Ward identities for elastic Burgulence, showing viscoelastic turbulence has weaker symmetry constraints than Navier-Stokes turbulence.","lead":"The paper formulates elastic turbulence in a path-integral framework and derives symmetry constraints (Ward identities) for a simplified 1D viscoelastic Burgers model. It provides a structural foundation for future nonperturbative renormalization group calculations of polymer-driven turbulence.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The X5 Ward identity (Eq. 149–152) is ill-defined for Schwartz-class test functions, and its regulator compatibility is acknowledged insufficient — the closure program in §7.3 rests on an unverified identity.","rationale":"The reader correctly identifies the X5 Ward identity convergence and regulator compatibility as the most load-bearing concern. The paper's central structural claim — that elastic turbulence has fewer symmetry constraints than Navier-Stokes — does not depend on X5 and stands on its own merits, supported by the systematic derivation of X1–X4 and the clear absence of Galilean-type constraints on the stress sector. However, the constructive part of the paper (closure schemes in §7.2.1 and §7.3) does depend on X5, and the identity's mathematical well-definedness, its bare-level validity, and its compatibility with the Wetterich regulator are all unverified. The CONDITIONAL verdict is appropriate: the framework is sound as a structural analysis, but the closure program cannot proceed without resolving the X5 issue. The paper itself is appropriately cautious, acknowledging that closure evaluation is 'outside the scope of the present work.' No adjustment to the verdict is needed.","tokens_in":27941,"tokens_out":3391,"duration_ms":136721,"concrete_test":"Compute both sides of the X5 Ward identity (Eq. 203, the δ′-coefficient relation: 0 = iΓ^{(u,2)}(0,iχ; q,Ω) − (q/χ)[Γ^{(2)}(q,Ω+iχ) − Γ^{(2)}(q,Ω)]) at the bare scale k=Λ using the leading-order derivative expansion ansatz (Eq. 176) with the power-law regulator from Eq. 157. If the two sides do not match at k=Λ, the regulator breaks X5 from the outset. Then evolve the Wetterich flow numerically in this truncation and check whether the identity is restored as k→0. If it is not restored, the assumption in §7.2.1 is unfounded and the §7.3 closure loses its basis.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the weakest link. The X5 generator (Eq. 111/115) produces f(t) = exp(χt) with χ potentially positive (e.g., A=0, B=1, α=1 gives χ = 1/(3W) > 0). The resulting Ward identity (Eq. 150) involves integrals ∫ φ(t) e^{χt} dt that diverge for φ ∈ S (Eq. 151). The authors' regularization — factoring out the divergence as δ(ω−i) — is analogous to the δ(0) factor in effective potential calculations, but that analogy is imperfect: δ(0) arises from a volume divergence in a well-defined integral, whereas here the integrand itself grows exponentially, making the pairing undefined before any volume factor appears. The physical argument (Eq. 152) that correlations decay as e^{−t/W} conflates properties of the full interacting theory with the bare-level identity: Ward identities are operator relations that must hold at the level of the regularized action, where bare correlations need not exhibit the physical decay rate. Furthermore, even if convergence is granted, the regulator preserving X5 (Eq. 157, a power-law R_k ∝ (|p|/k)^C) is acknowledged to be insufficient to suppress IR singularities of the propagator for all truncations (end of §7). In §7.2.1, the authors sidestep this by assuming the identity is restored at k→0 without enforcement along the flow — an unverified assumption. If X5 fails, the BMW-type closure in §7.3 (Eqs. 203–210) loses its key non-trivial constraint, reducing to the velocity-sector-only closure already available from X1/X2. The structural result about symmetry deficit (the comparison with Navier-Stokes) survives, as it depends on the absence of Galilean-type constraints on the stress sector, not on X5's validity.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript formulates elastic and elasto-inertial turbulence (Oldroyd-B model) within the Martin-Siggia-Rose (MSR) path-integral framework and develops a systematic, source-extended Lie-group algorithm to derive Ward identities directly from the Euler-Lagrange equations. As a dimensionally reduced model, the authors introduce an extended Burgers equation that retains the characteristic coupling between the extra stress and velocity gradient. The paper derives the Ward identities for this model, analyzes the resulting constraints on closure schemes, and outlines two complementary truncation strategies for the Wetterich flow equation: a leading-order derivative expansion and a momentum-resolving BMW-type closure. The central structural finding is that, unlike in Navier-Stokes turbulence, the stress and response-stress sectors in viscoelastic turbulence carry far fewer protecting symmetries, leading to a larger running theory space and potential mean-field instabilities.","tokens_in":28301,"tokens_out":1899,"duration_ms":141491,"significance":"The application of the functional renormalization group (fRG) to elastic turbulence is a genuinely novel undertaking with clear motivation. The source-extended Lie-group algorithm for systematically deriving Ward identities is a methodological contribution in its own right and is presented in a self-contained, algorithmic manner. The identification of the specific structural obstruction—namely, the comparative lack of constraints on the stress sector compared to the Galilean-protected velocity sector in Navier-Stokes—is a valuable insight for the community. The scaling analysis yielding z=0 from scale-invariant energy injection is a concrete, falsifiable prediction. The two proposed closure schemes (Eqs. 176–210) provide a concrete roadmap for future quantitative calculations.","major_comments":[{"comment":"§7, Eqs. (149)–(152): The Ward identity generated by X5 involves the function f(t) = exp(χt) (Eq. 115), where χ can be positive (e.g., A=0, B=1, α=1 gives χ = 1/(3W) > 0). The resulting identity (Eq. 150) involves integrals of the form ∫ φ(t) e^{χt} dt, which diverge for Schwartz-class test functions φ ∈ S (Eq. 151). The authors propose to regularize this by factoring out the divergence as an ill-defined δ(ω−i), drawing an analogy to the δ(0) factor in effective potential calculations. This analogy is imperfect: δ(0) in the effective potential arises from a volume divergence in a well-defined integral, whereas here the integrand itself grows exponentially, making the distributional pairing undefined before any volume factor appears. Furthermore, the physical argument (Eq. 152) that correlations decay as e^{−t/W} conflates properties of the full interacting theory with the bare-level Ward","section":null},{"comment":"identity, which must hold at the level of the regularized action where bare correlations need not exhibit the physical decay rate. This issue is load-bearing: the BMW-type closure in §7.3 (Eqs. 203–210) relies on the X5 Ward identity as its key non-trivial constraint on the stress sector. If the identity is ill-defined, the closure reduces to the velocity-sector-only constraints already available from X1/X2. The authors should either provide a rigorous justification for the regularization (e.g., by demonstrating that the Ward identity can be defined as a distributional limit in a suitable test-function space) or explicitly delineate the status of the X5 identity as a formal/asymptotic relation whose validity must be verified a posteriori.","section":null},{"comment":"§7.2.1, Eqs. (176)–(177) and surrounding text: The leading-order derivative expansion ansatz implements the X5 Ward identity (149) 'in its unmodified state,' explicitly assuming that 'the Ward identity be restored for k→0 and that the final fixed point that is approached is not influenced by the perturbation of the trajectory through enforcement along the entire flow.' This is an unverified assumption. The regulator preserving X5 (Eq. 157, a power-law R_k ∝ (|p|/k)^C) is acknowledged at the end of §7 to be insufficient to suppress all IR singularities of the propagator for all truncations. The authors should clarify whether the fixed point reached under this ansatz is robust to the regulator breaking, or at minimum discuss the potential sensitivity of the stability analysis (Eqs. 179–183) to this assumption.","section":null},{"comment":"§7.3, Eqs. (203)–(210): The BMW-type closure relies on evaluating vertex functions at a complex frequency shift ω = iχ (Eqs. 203, 204). The closure in Eq. (206) then assumes |Ω| ≫ |χ| to replace σ-sector vertices by their zero-momentum values. The regime of validity of this approximation is not specified. Since χ = 1/(3W) for the Oldroyd-B case (A=0, B=1, α=1), the condition |Ω| ≫ |χ| may restrict the closure to frequencies well above the polymer relaxation rate, potentially excluding the physically relevant elastic subrange. The authors should discuss whether this closure is applicable in the scaling regime of interest or whether it is limited to a UV asymptotic regime.","section":null}],"minor_comments":[{"comment":"§1: 'teh econd' should be 'the second' (in the sentence beginning 'Within the Lundgren [23], Monin [26], Novikov [27] hierarchy...').","section":null},{"comment":"§1: 'clssical Navier Stokes turublence' should be 'classical Navier-Stokes turbulence'.","section":null},{"comment":"§1: 'nomber of degrees of freemdom' should be 'number of degrees of freedom'.","section":null},{"comment":"§7: 'study pf' should be 'study of'.","section":null},{"comment":"§7: 'fuction' should be 'function' (in the text following Eq. 149).","section":null},{"comment":"§6.0.2, Eq. (111): The generator X5 is written with terms involving σ and σ̄, but the notation for the response stress field is inconsistent with the rest of the paper (sometimes σ̄, sometimes σ̃). Please standardize.","section":null},{"comment":"§5.1, Eq. (60): The δ(0) prefactor is stated to be absorbable into a normalization defining a spatial density. It would help the reader to state explicitly that this is the standard volume factor V/(2π)^d.","section":null},{"comment":"§7.1, Eq. (165): The dimension [ε] = 2z is derived, and z=0 follows from [ε]=0. It would strengthen the presentation to explicitly note that this result (z=0 for elastic turbulence) holds in arbitrary spatial dimension d, as mentioned in the parenthetical but not emphasized.","section":null},{"comment":"§6.0.3: The symmetry X7 (Eq. 133) and its variation δX7(L) = d·f(t)·L (Eq. 134) are introduced but their implications are not discussed. If these are not used subsequently, a brief comment on their (ir)relevance would be helpful.","section":null},{"comment":"References: Several entries have minor formatting issues (e.g., [16] uses a non-standard DOI format '10.1103/pbtf-rn7d').","section":null}],"recommendation":"major_revision","confidential_remarks":"The reader's report and stress-test note correctly identify the X5 Ward identity convergence issue as the central weakness. On reading the manuscript, I confirm that this concern lands: the regularization procedure in Eqs. (151)–(152) is not rigorously justified, and the conflation of bare-level operator identities with physical decay properties of the full theory is a genuine gap. However, I note that the paper's central structural claim—the comparative lack of symmetry constraints on the stress sector—does not depend on X5 alone; it holds even if X5 is discarded, since the remaining identities (X1–X4) only protect the velocity sector. The X5 identity is needed specifically for the closure program in §7.3, not for the structural conclusions of §6. This suggests the paper could be restructured to present the X5-based closure as a conditional proposal pending rigorous justification, rather than as an established result. The paper is a methods/foundation paper with no numerical results or data comparison, which is appropriate for the journal's scope, but the authors should be asked to clarify the status of their central tool (the X5 identity) rather than presenting it as established."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and substantive reading of the manuscript. The referee's comments identify genuine technical gaps in the treatment of the X5 Ward identity and its consequences for both closure schemes. We address each point below and indicate revisions where they are warranted.","responses":[{"response":"The referee is correct that the regularization argument as currently stated is not rigorous. We acknowledge two distinct issues: (1) The pairing ⟨φ, e^{χt}⟩ is not defined for φ ∈ S when χ > 0, and the analogy to δ(0) in effective potential calculations is indeed imperfect—the δ(0) there arises from a volume divergence in an otherwise well-defined integral, whereas here the integrand itself grows exponentially. (2) The physical argument in Eq. (152) that correlations decay as e^{−t/W} is an argument about the full interacting theory, not about the bare-level Ward identity, which must hold at the level of the regularized action where bare correlations need not exhibit the physical decay rate. We cannot resolve this by simply invoking physical decay at the bare level. We will revise the manuscript to explicitly delineate the status of the X5 identity as a formal relation whose validity must be verified a posteriori—specifically, by checking at the level of the truncated effective action whether the full (dressed) correlators decay sufficiently fast to render the distributional pairing well-defined. We will also add a remark that a rigorous justification would require demonstrating that the Ward identity can be defined as a distributional limit in a test-function space restricted to functions with sufficient exponential decay (e.g., the Gelfand–Shilov space S_α rather than the full Schwartz space S), and we will note this as an open problem. We agree that the consequences for the BMW closure are significant: if the X5 identity is only formal, then the stress-sector constraints it provides are not guaranteed, and the closure does reduce to the velocity-sector-only constraints from X1/X2 unless and until the identity is validated at the dressed level.","revision_made":"partial","referee_comment":"§7, Eqs. (149)–(152): The Ward identity generated by X5 involves f(t) = exp(χt) with χ potentially positive, leading to divergent integrals. The δ(0) analogy is imperfect, and the physical argument conflates bare-level and interacting correlations. If the identity is ill-defined, the BMW closure reduces to velocity-sector-only constraints."},{"response":"The referee correctly identifies that the assumption stated in the text surrounding Eqs. (176)–(177) is unverified. We agree that this should be made explicit rather than presented as a working hypothesis without qualification. In the revised manuscript, we will add a discussion of the potential sensitivity of the stability analysis (Eqs. 179–183) to this assumption. Specifically, we will note that: (a) the power-law regulator R_k ∝ (|p|/k)^C that preserves X5 is, as we already acknowledge at the end of §7, insufficient to suppress all IR singularities of the propagator for general truncations; (b) the fixed point reached under the ansatz (177) may therefore depend on the manner in which the X5 identity is enforced (or approximately enforced) along the flow; and (c) a robustness check would require varying the regulator shape and verifying that the fixed-point structure is stable, which we leave as a numerical task for future work. We will reframe the presentation to make clear that Eqs. (176)–(177) constitute a working ansatz whose self-consistency must be verified, not a rigorously justified truncation.","revision_made":"partial","referee_comment":"§7.2.1, Eqs. (176)–(177): The assumption that the Ward identity is restored for k→0 and that the fixed point is not influenced by enforcement along the entire flow is unverified. The regulator preserving X5 is acknowledged to be insufficient to suppress all IR singularities."},{"response":"The referee raises a valid and important point. For the Oldroyd-B case (A=0, B=1, α=1), we have χ = 1/(3W), so the condition |Ω| ≫ |χ| restricts the closure to frequencies well above the polymer relaxation rate 1/W. This is indeed the regime of the elastic subrange's UV tail, not the elastic subrange itself where Ω ~ 1/W. We will add an explicit discussion of the regime of validity in the revised manuscript. We note that the closure as written is an asymptotic approximation valid in the UV regime |Ω| ≫ |χ|, and that extending it to the physically relevant elastic subrange would require retaining the full frequency dependence of the σ-sector vertices rather than replacing them by their zero-momentum values. This is a nontrivial extension that increases the computational cost substantially, as it prevents the algebraic reduction of the stress-sector vertices. We will state this limitation clearly and note that the closure in its present form should be understood as a UV asymptotic scheme, with the elastic subrange requiring either a modified closure strategy or a full numerical treatment of the frequency-dependent vertices.","revision_made":"partial","referee_comment":"§7.3, Eqs. (203)–(210): The BMW-type closure assumes |Ω| ≫ |χ|, which may restrict the closure to frequencies well above the polymer relaxation rate, potentially excluding the physically relevant elastic subrange."}],"tokens_in":28115,"tokens_out":1449,"duration_ms":50724,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper develops a source-extended Lie-group algorithm for systematically deriving Ward identities from Euler-Lagrange equations with source terms, and applies it to viscoelastic Burgulence. The genuinely new result is twofold — the method itself (Section 6, Eqs. 86–91), and the structural finding that elastic turbulence has substantially weaker symmetry constraints than Navier-Stokes turbulence, because the stress sector lacks anything like Galilean invariance. That diagnostic result is real and important for the subfield. The response-velocity sector is protected from renormalization (Eq. 147), and the constitutive couplings are restricted in a derivative expansion (Eqs. 178–190). The MSR construction (Sections 3–5) is standard but correctly executed, and the scaling analysis giving z=0 from scale-invariant energy injection is clean and not circular — it's a standard physical postulate, not an input dressed up as an output. Credit is due for laying out the closure program concretely rather than just gesturing at it. The derivative expansion ansatz (Eq. 176) and the BMW-type closure (Eqs. 203–210) are worked out in real detail, even if not numerically evaluated. Now the soft spot, which the reader and stress-test correctly identify: the X5 Ward identity (Eqs. 149–152) involves integrals with exponential growth in time. The authors acknowledge this and propose factoring out the divergence as an ill-defined delta function, analogous to delta(0) in effective potential calculations. The stress-test note is right that this analogy is imperfect — delta(0) arises from a volume divergence in a well-defined integral, whereas here the integrand itself grows exponentially, making the pairing undefined before any volume factor appears. The physical argument (Eq. 152) that correlations decay as exp(-t/W) conflates properties of the full interacting theory with the bare-level identity, which must hold at the level of the regularized action. This is a real gap. If X5 fails, the BMW closure in Section 7.3 loses its key non-trivial constraint and reduces to the velocity-sector-only closure already available from X1/X2. The structural result about symmetry deficit survives regardless, since it depends on the absence of Galilean-type constraints on the stress sector, not on X5's validity. The regulator compatibility issue (end of Section 7) is a secondary concern — the authors are upfront that the power-law regulator preserving X5 is insufficient for all truncations, and they assume restoration at k→0 without enforcement along the flow. This is a known limitation of the approach, acknowledged honestly. This paper is for researchers working on fRG approaches to turbulence who need the symmetry infrastructure before quantitative calculations are possible. It deserves a serious referee — the methodological contribution is genuine and the structural result about symmetry deficit is correct and useful. The X5 issue needs to be addressed before the closure program can proceed, but that's a problem for the next paper, not a reason to reject this one.","headline":"New systematic Ward-identity method for viscoelastic turbulence; one load-bearing identity needs rigorous justification","tokens_in":29000,"tokens_out":695,"would_cite":true,"duration_ms":85597,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.57.-s","47.27.ef","05.10.Cc"],"model":"glm-5.2","headline":"Ward identities expose why elastic turbulence resists standard renormalization","keywords":[],"falsifier":"If the exponentially growing Ward identity (X_5) cannot be consistently regularized within the Wetterich flow, the momentum-resolving closure of Section 7.3 loses its primary constraint, and the derivative-expansion scheme of Section 7.2 must rely solely on the weaker identities X_1 and X_2, which do not constrain the stress sector.","tokens_in":28068,"feed_emoji":"🌀","tokens_out":918,"duration_ms":64632,"temperature":0.7,"pith_summary":"The paper develops a systematic algorithm to derive Ward identities—exact nonperturbative constraints on correlation functions—for viscoelastic turbulence by extending Lie-group symmetry analysis to the source-extended Euler-Lagrange equations of the Martin-Siggia-Rose path integral. Applied to a dimensionally reduced model called elastic Burgulence (a Burgers equation coupled to an elastic stress tensor via an Oldroyd-type constitutive law), the authors find that the symmetry structure is substantially weaker than in Navier-Stokes turbulence. In Navier-Stokes, Galilean invariance produces strong Ward identities that pin down the zero-momentum sector and enable controlled closure schemes. In the viscoelastic case, the stress and response-stress sectors carry no comparable protecting symmetry: the mean stress need not vanish, the effective action has a larger space of running couplings, and dynamical instabilities of the mean field can arise. The paper identifies which sectors are protected (the response-velocity sector does not renormalize) and which remain unconstrained, then uses these constraints to restrict the functional form of constitutive couplings in a derivative expansion and to outline two complementary closure strategies for the Wetterich renormalization-group flow equation.","feed_headline":"Symmetry deficit found in elastic turbulence's renormalization structure","feed_subtitle":"Ward identities show viscoelastic stress lacks the Galilean protection that makes Navier-Stokes closures work, blocking standard approx","key_machinery":"Ward identities derived via a source-extended Lie-group symmetry algorithm applied to the MSR action for viscoelastic Burgers equations; the Wetterich flow equation for the scale-dependent effective action; a derivative expansion ansatz for the effective action constrained by the derived identities.","core_discovery":"The central discovery is that the source-extended Lie-group algorithm systematically produces Ward identities for elastic Burgulence, and these identities reveal a structural deficit: unlike Navier-Stokes turbulence, where Galilean invariance constrains the zero-momentum sector and enables closure, the viscoelastic stress sector is comparatively unprotected. The response-velocity sector is shielded from renormalization (Eq. 147), and the constitutive couplings are restricted in a derivative expansion (Eqs. 178-190), but the stress and response-stress sectors remain largely unconstrained, making nonperturbative closures significantly harder.","pith_inferences":[],"forward_implications":["The symmetry deficit identified here means that closure schemes for elastic turbulence cannot directly borrow the zero-momentum-sector strategy that proved successful for Navier-Stokes and KPZ; new approximation architectures are needed for the stress sector.","The derivative-expansion truncation (Eqs. 178-190) provides a concrete starting ansatz for numerical computation of the renormalization-group flow, including the stability boundary of mean-field stress configurations (Eq. 172).","The scaling exponent z=0 derived from constant energy-injection rate (Eq. 165) applies not only to 1D Burgulence but to elastic turbulence in arbitrary spatial dimensions, offering a universal constraint on the fixed-point structure.","The momentum-resolving closure outlined in Section 7.3, if numerically tractable, would yield the first nonperturbative prediction of kinetic and polymeric energy spectra for elastic turbulence, for which no analytical scaling exponents currently exist."],"fun_headline_variants":["Elastic turbulence lacks symmetry protection needed for standard closures","Ward identities expose unprotected stress sector in elastic Burgulence","Viscoelastic stress lacks Galilean protection, complicating closures","Structural deficit in elastic turbulence blocks nonperturbative closures","Elastic Burgulence stress sector remains unconstrained by symmetries"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The Ward identity generated by the symmetry X_5 (Eqs. 149-150) involves a time integral with exponential growth that is formally divergent. The authors argue it exists because physical stress correlations decay exponentially in time, and propose factoring out the divergence as an ill-defined delta function. This regularization is not rigorously justified, and if it fails, the key constraint on the constitutive-sector closure in Section 7.3 is lost.","fun_headline_variants_meta":{"raw":{"variants":["Elastic turbulence lacks symmetry protection needed for standard closures","Ward identities expose unprotected stress sector in elastic Burgulence","Viscoelastic stress lacks Galilean protection, complicating closures","Structural deficit in elastic turbulence blocks nonperturbative closures","Elastic Burgulence stress sector remains unconstrained by symmetries"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":885,"prompt_tokens":412,"completion_tokens":473,"prompt_tokens_details":null},"tokens_in":412,"tokens_out":473,"duration_ms":21126,"temperature":1.0,"reasoning_tokens":474,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-05T12:26:21.551508+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the exponentially growing Ward identity (X_5) cannot be consistently regularized within the Wetterich flow, the momentum-resolving closure of Section 7.3 loses its primary constraint, and the derivative-expansion scheme of Section 7.2 must rely solely on the weaker identities X_1 and X_2, which do not constrain the stress sector.","supporting_citations":[],"review_version":1}