{"id":"b456f020-dda0-4fb0-9bae-1fdb73bb6207","arxiv_id":"2505.01814","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Second-order spin hydrodynamics is extended with nonlocal memory corrections, adding new acceleration-dependent terms to diffusion, rotational stress, and heat-current equations.","lead":"This paper adds new correction terms to the equations used to model spinning quark-gluon plasma as a fluid. The terms come from memory effects in the quantum statistical operator and depend on how fast the fluid accelerates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Curie-theorem decompositions (84)-(89) assume isotropy despite a nonzero spin density; unexamined S-dependent tensor structures could alter the claimed second-order terms.","rationale":"The reader identified the isotropy assumption in the presence of spin density as the weakest assumption, and my reading agrees: this is the point where the central claim is least secure. Every claimed new second-order term, including the headline Eq. (107), depends on the tensor decompositions (84)-(89); those decompositions require full rotational invariance, which a nonzero S^{\\mu\\nu} breaks. The paper provides no argument excluding S-dependent tensor structures, so the constitutive relations and Kubo formulas are potentially incomplete. This is a genuine load-bearing concern rather than a manufactured one: spin is the novel ingredient of the paper, and the anisotropy it induces is the natural first place where the standard Curie-theorem step can fail. The reader's CONDITIONAL verdict already reflects this uncertainty, including the need to verify imported identities from Ref. [1] and the lack of independent numerical checks. Since my concern matches the reader's weakest assumption and does not push the verdict to a different category, I leave the verdict unchanged. A concrete microscopic calculation of the correlator components in a spin-polarized system would settle whether the concern lands.","tokens_in":37114,"tokens_out":3111,"duration_ms":35907,"concrete_test":"Perform a microscopic test of Eq. (84) in a weakly interacting spin-polarized system, e.g. a free or weakly coupled Dirac gas at finite spin chemical potential \\omega, with the fluid at rest and the spin density along the z-axis. Compute the two-point functions C^{xy,xy}(x,x_1) and C^{xz,xz}(x,x_1) (with the appropriate Zubarev correlator at local equilibrium). If these two components are unequal, the single-form-factor ansatz (84) fails and the claims built on it, including Eq. (107), are not established. Alternatively, enumerate the full covariant tensor basis for the correlator under the residual symmetry group preserving S^{\\mu\\nu} and u^{\\mu}; if any independent S-containing projector appears, the paper must either show it vanishes or include its contribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new results, e.g. Eq. (107), are reached by substituting the one-form-factor decompositions (84)-(89) into the nonlocal correction formula (76). These decompositions are justified by 'medium isotropy' plus the orthogonality constraints (19). But the local-equilibrium state considered in this paper carries a nonzero spin density S^{\\mu\\nu}, which is treated as a zeroth-order thermodynamic variable, and a nonzero S^{\\mu\\nu} breaks full rotational invariance down to the subgroup that preserves S. Curie's theorem in its textbook form applies to isotropic media; in a medium with a preferred spin direction, additional tensor structures built from S^{\\mu\\nu} (and compatible with transversality and symmetry) are allowed in every two-point correlator of rank-two tensors. The paper never argues that these extra structures vanish. For example, in the local rest frame with the spin vector along z, decomposition (84) forces the correlator components C^{xy,xy} and C^{xz,xz} to be equal, whereas only O(2) rotations around S are required to leave the state invariant, so these components may differ. If such S-dependent form factors are present, the Kubo formulas (95)-(101), the first-order relations (90)-(94), and every nonlocal second-order correction, including Eq. (107), are incomplete: transport coefficients would become spin-dependent and additional spin-coupled gradient terms would appear. The issue is internal rather than a disagreement with external consensus, because the paper's own setup includes a nonvanishing zeroth-order spin density, and the relevant argument is simply absent. In addition, the paper imports several key identities, such as (104), from the appendix of Ref. [1] without re-derivation; this is a secondary concern, but the anisotropy of the local-equilibrium ensemble is the load-bearing gap for the claimed novelty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Zubarev nonequilibrium statistical operator approach to relativistic spin hydrodynamics. Starting from the local-equilibrium operator with a spin chemical potential and from a second-order expansion of the statistical operator, the authors compute nonlocal corrections to the dissipative fluxes that arise from two-point correlation functions of tensors evaluated at distinct spacetime points. They obtain new second-order terms in the shear stress tensor, bulk viscous pressure, charge-diffusion currents, rotational stress tensor, and boost heat vector, all involving the comoving derivative of the flow velocity, and they express the corresponding transport coefficients through Kubo-type formulas in terms of retarded Green's functions. The paper also writes relaxation-type equations for the dissipative currents and states that the nonlocal correction to the fully antisymmetric spin flux vanishes.","tokens_in":37356,"tokens_out":13509,"duration_ms":137720,"significance":"If the results are correct, they would show that the second-order spin hydrodynamics derived in the authors' earlier work [2] is incomplete, adding acceleration-coupled terms analogous to those found for spinless fluids in Ref. [1]. The manuscript's strengths are its systematic use of the Zubarev formalism, the explicit Kubo expressions for two- and three-point transport coefficients, and the clear separation of first-order, nonlocal second-order, and three-point-correlation contributions. The central claim is, however, conditional on the isotropy assumption used to reduce two-point correlators to single form factors, and on frequency-derivative identities taken from the spinless companion paper; both points need to be resolved before the results can be regarded as complete.","major_comments":[{"comment":"The Curie-theorem decompositions (84)-(89) assume full rotational isotropy of the local-equilibrium state. This assumption is inconsistent with the paper's own setup, in which the spin density S^{μν} is a zeroth-order thermodynamic variable (Section II.A) and therefore defines a preferred direction. In a local rest frame with spin along the z-axis the symmetry is only O(2), so, for example, Eq. (84) forces the correlator components C^{xy,xy} and C^{xz,xz} to be equal, whereas S-dependent tensor structures compatible with the orthogonality and symmetry conditions (19) could split these components. The paper never shows that such S-dependent form factors vanish. Because Eqs. (90)-(101), (103), (107), and the analogous second-order results all rely on (84)-(89), this is a load-bearing gap rather than a presentation issue.","section":"Section III.A, Eqs. (84)-(89)"},{"comment":"Equation (104) as printed is not a valid tensor identity: the left side has free indices μνρσ and τ, while the right side has only τ. Since Eq. (107) is obtained by substituting Eq. (104) into Eq. (103), the authors must provide the corrected statement (presumably a first-moment identity for the scalar correlator ∫ d^4x_1 (π̂_{λη},π̂^{λη})(x_1-x)^τ) and show how the projector in Eq. (103) is handled.","section":"Section III.B.1, Eq. (104)"},{"comment":"Several load-bearing identities that convert first moments of two-point correlators into the transport coefficients η~, ζ~, χ~, γ~, and λ~ are quoted from Appendix B of Ref. [1], which deals with spinless fluids. In the present paper the local-equilibrium operator contains the spin-chemical-potential term in Eq. (5), and the operators p̂*, Ĵ^μ_a, and related quantities are modified by spin contributions, e.g. p̂* = p̂*_spinless - K_{αβ}Ŝ^{αβ} in Eq. (46). The authors do not demonstrate that these modifications leave the first-moment identities unchanged at the required order. Without such a demonstration, the nonlocal second-order transport coefficients in Eqs. (107), (147), (176), (197), and (225) are not fully derived.","section":"Section III.B, Eqs. (104), (134)-(136), (165)-(167), (194), (215)-(217)"}],"minor_comments":[{"comment":"The title and body contain typographical errors, including 'spin hydrodynamic s' in the title and 'dissiptive' and 'hydrodyn amics' in the body.","section":"Title and body text"},{"comment":"The symmetry property of the three-point correlator is stated without proof; since it is used repeatedly, a short derivation or a precise reference would improve readability.","section":"Eq. (15)"},{"comment":"The notation for Z is inconsistent: Eq. (57) uses Z^{μν} while the definition in Eq. (58) has lower indices. Aligning the index positions would avoid confusion.","section":"Eqs. (57)-(58)"},{"comment":"The definitions of Γ~ and δ~_a contain explicit θ^{-1} factors. Although the final products in Eq. (147) are regular, the intermediate quantities are formally singular at θ=0; presenting the original regular expressions (144)-(145) directly would avoid this formal singularity.","section":"Eq. (146)"},{"comment":"The sentence 'Symmetry considerations allow the omission of the term γ_{φqq}M^{[μ}M^{ν]}' should say explicitly that the antisymmetrized product of a symmetric tensor vanishes identically, rather than invoking an unspecified symmetry consideration.","section":"After Eq. (213)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a direct extension of Refs. [1] and [2] and relies heavily on those works for the derivation of many identities. The novelty is incremental but suitable for a specialist journal. The main obstacle is the isotropy assumption in the presence of a nonzero spin density; I would encourage the editor to require the authors either to justify the vanishing of S-dependent tensor structures or to extend the tensor decompositions accordingly before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know: the paper does a real piece of work, extending Harutyunyan and Sedrakian's nonlocal second-order corrections to spin degrees of freedom; and it has a gap that could change the results—the authors assume Curie-style isotropy for two-point functions in a medium with nonzero spin density, and never justify why spin-dependent tensor structures should vanish.\n\nFirst the credit. The derivation is careful under the stated power counting. The new two-point nonlocal corrections to the rotational stress tensor and the boost heat vector, including the \\dot u^\\mu \\nabla^\\nu \\alpha terms, are genuinely new compared to both Ref [1] and their own earlier work [2]. Expressing every transport coefficient as a zero-frequency derivative of a retarded Green's function is a clean result. The algebra is dense but internally consistent, and they don't try to hide the reliance on the appendix of Ref [1] for the frequency-derivative identities—though that reliance does make the paper's technical core partly a black box.\n\nThe soft spot is real. Their local equilibrium state carries a nonzero, zeroth-order spin density S^{\\mu\\nu}. That breaks rotational invariance down to the subgroup preserving S. Curie's theorem in the form used in Eqs. (84)-(89) requires full isotropy, so additional S-dependent tensor structures are allowed in every two-point correlator. The equality e.g. C^{xy,xy}=C^{xz,xz} for the shear correlator would not hold in general. If such structures survive, the Kubo formulas (95)-(101), the first-order relations (90)-(94), and every nonlocal second-order correction—including the headline Eq. (107)—are incomplete. The paper never argues that these extra terms vanish or are subleading. That's not a disagreement with an external consensus; it's an internal gap.\n\nThere's also a minor reproducibility issue: many load-bearing identities are quoted from Ref [1] without derivation. That's acceptable in a companion paper, but it makes independent verification harder. No numerics or magnitude estimates, so we don't know how big the new terms are.\n\nWho's it for: people working on relativistic spin hydrodynamics, especially Zubarev-formalism users. The field is active and this is a serious derivation, so it deserves referee time. I'd send it to peer review, but the referee should push on the isotropy assumption. The authors need to either justify the assumption in their setup or redo the decompositions to include S-dependent form factors.","headline":"Serious formal extension of Harutyunyan–Sedrakian to spin hydrodynamics, with new nonlocal terms, but a load-bearing isotropy assumption is left unjustified.","tokens_in":37971,"tokens_out":4952,"would_cite":false,"duration_ms":51960,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Second-order spin hydrodynamics is incomplete without nonlocal two-point correlation terms.","keywords":["spin hydrodynamics","Zubarev nonequilibrium statistical operator","second-order hydrodynamics","nonlocal correlations","transport coefficients","retarded Green's functions","quark-gluon plasma","spin polarization"],"falsifier":"Compute the two-point retarded Green's functions of the dissipative currents in a background with nonzero spin density or spin chemical potential: if any correlator contains independent tensor structures built from $S^{\\mu\\nu}$ beyond the isotropic projectors, Curie's-theorem reduction fails. A cheaper check is to derive the second-order shear relaxation equation from kinetic theory with a nonlocal collision term and see whether the term $2\\tilde{\\eta}Th^{-1}\\sum_a n_a\\dot{u}^{\\langle\\mu}\\nabla^{\\nu\\rangle}\\alpha_a$ appears with the same coefficient.","tokens_in":36901,"feed_emoji":"🌀","tokens_out":7983,"duration_ms":75083,"temperature":0.7,"pith_summary":"This paper extends second-order relativistic spin hydrodynamics to include nonlocal corrections that arise when two-point correlation functions are evaluated at distinct spacetime points. It claims such correlations produce memory-effect terms, quadratic in gradients and proportional to the comoving derivative of the flow velocity, that prior derivations omitted. The result is a complete set of second-order constitutive relations for the shear stress tensor, bulk viscous pressure, charge diffusion currents, rotational stress tensor, and boost heat vector, with every new transport coefficient expressed through retarded Green's functions. If correct, the earlier second-order spin hydrodynamic equations are missing these acceleration-coupled contributions.","feed_headline":"Spin hydrodynamics gains acceleration-coupled second-order terms","feed_subtitle":"Nonlocal memory terms add acceleration couplings to spin-fluid equations; Kubo formulas fix their coefficients.","key_machinery":"The central object is Zubarev's nonequilibrium statistical operator $\\hat{\\rho} = Q^{-1}e^{-\\hat{A}+\\hat{B}}$, where $\\hat{A}$ fixes the local equilibrium thermodynamic parameters and $\\hat{B}$ carries a memory integral over the dissipative operator $\\hat{C}$. The argument works by expanding $\\hat{C}(x_1)$ about $x$ through first order in $x_1-x$, inserting that expansion into two-point correlators, and using Curie's theorem to reduce each correlator to one scalar form factor times an isotropic projector. This turns nonlocality into comoving derivatives of thermodynamic forces, producing the tilded coefficients $\\tilde{\\eta}$, $\\tilde{\\zeta}$, $\\tilde{\\chi}$, $\\tilde{\\gamma}$, $\\tilde{\\lambda}$ as first frequency derivatives of the corresponding retarded Green's functions. The three-point correlator terms from the previous formulation are retained, so the final constitutive and relaxation equations contain both the old nonlinear terms and the new nonlocal ones.","core_discovery":"The paper derives, from Zubarev's nonequilibrium statistical operator, that nonlocal two-point correlations of different-rank tensor operators generate new second-order terms in spin hydrodynamics. Its central example is Eq. (107): $\\langle \\pi^{\\mu\\nu}\\rangle^1_2 = 2\\tilde{\\eta}\\,(\\Delta^{\\mu\\nu\\rho\\sigma}D\\sigma_{\\rho\\sigma}+\\Gamma\\theta\\sigma^{\\mu\\nu}) + 2\\tilde{\\eta}T h^{-1}\\sum_a n_a \\dot{u}^{\\langle\\mu}\\nabla^{\\nu\\rangle}\\alpha_a$, where $\\tilde{\\eta}$ is the frequency derivative of the shear Kubo function. Analogous terms appear in the bulk pressure, charge diffusion currents, rotational stress tensor, and boost heat vector, and each coefficient is given as a derivative of a retarded Green's function. The same mechanism yields no nonlocal two-point correction to the fully antisymmetric spin flux $\\varpi^{\\lambda\\mu\\nu}$.","pith_inferences":["Inference: the same memory mechanism should appear in kinetic-theory or entropy-current derivations of second-order spin hydrodynamics; deriving those equations with a nonlocal collision term would test whether the acceleration couplings are scheme-independent.","Inference: in heavy-ion collisions, the new $\\dot{u}^{\\mu}\\nabla^{\\nu}\\alpha_a$ terms may be sizable in the expanding fireball where acceleration and baryon chemical potential gradients coexist, possibly affecting the rapidity dependence of $\\Lambda$ and $\\bar{\\Lambda}$ polarization.","Inference: the Curie-theorem reduction is the fragile step; if a nonzero spin density $S^{\\mu\\nu}$ forces additional tensor structures into the correlators, further terms would appear at the same order. A systematic expansion of the correlators in powers of $S^{\\mu\\nu}$ would quantify that correction.","Inference: because the coefficients are frequency derivatives of Green's functions, weakly coupled QCD or holographic computations could provide concrete values and reveal which new terms are numerically dominant."],"forward_implications":["The relaxation-type equations for $\\pi^{\\mu\\nu}$, $\\Pi$, $J_c^{\\mu}$, $\\varphi^{\\mu\\nu}$, and $q^{\\mu}$ acquire new source terms involving $\\dot{u}^{\\mu}$, the acceleration of the fluid, multiplying gradients of thermal potentials.","All new coefficients are compute-facing: they are second frequency derivatives of retarded Green's functions, so thermal field theory, lattice QCD, or holographic methods can in principle evaluate them.","Stability and causality analyses of second-order spin hydrodynamics should be redone with these terms included, since the added $\\dot{u}^{\\mu}$ couplings change the linear mode structure.","The nonlocal correction to the fully antisymmetric spin flux vanishes, so that sector of the theory is unchanged by the mechanism.","The new terms survive in the spinless limit only for shear, bulk, and diffusion; the rotational-stress and boost-heat-vector terms are specific to spin fluids."],"supporting_citations":[{"why":"Supplies the method: two-point correlations of different-rank tensors at distinct points generate second-order nonlocal corrections, which this paper generalizes to spin.","marker":"[1]"},{"why":"The authors' earlier second-order spin hydrodynamics, built on Zubarev's operator, that this paper extends and whose equations gain the new terms.","marker":"[2]"},{"why":"Provides the second-order expansion of Zubarev's statistical operator in powers of the dissipative operator, used throughout the derivation.","marker":"[47]"},{"why":"Equilibrium density-matrix prediction of vorticity-induced hyperon polarization, the physical motivation for spin hydrodynamics.","marker":"[5]"},{"why":"Supports the canonical pseudogauge choice by matching the thermodynamics of axial current interactions to Zubarev's local equilibrium operator.","marker":"[46]"}],"fun_headline_variants":["Nonlocal spin hydrodynamics gains acceleration-coupled terms","Spin fluid equations get new second-order nonlocal terms","Zubarev method yields extra spin transport coefficients","New spin hydrodynamics terms from nonlocal correlations","Acceleration couplings appear in spin hydrodynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that local equilibrium is isotropic enough for Curie's theorem to collapse every two-point correlator into a single scalar form factor, even though a nonzero spin density $S^{\\mu\\nu}$ itself breaks rotational invariance; if spin-dependent tensor structures appear, the constitutive relations and their Kubo formulas are incomplete.","fun_headline_variants_meta":{"raw":{"variants":["Nonlocal spin hydrodynamics gains acceleration-coupled terms","Spin fluid equations get new second-order nonlocal terms","Zubarev method yields extra spin transport coefficients","New spin hydrodynamics terms from nonlocal correlations","Acceleration couplings appear in spin hydrodynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000162,"raw_usage":{"total_tokens":1191,"prompt_tokens":850,"completion_tokens":341,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":269}},"tokens_in":466,"tokens_out":341,"duration_ms":3569,"temperature":1.0,"reasoning_tokens":269,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:09:50.206501+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two-point retarded Green's functions of the dissipative currents in a background with nonzero spin density or spin chemical potential: if any correlator contains independent tensor structures built from $S^{\\mu\\nu}$ beyond the isotropic projectors, Curie's-theorem reduction fails. A cheaper check is to derive the second-order shear relaxation equation from kinetic theory with a nonlocal collision term and see whether the term $2\\tilde{\\eta}Th^{-1}\\sum_a n_a\\dot{u}^{\\langle\\mu}\\nabla^{\\nu\\rangle}\\alpha_a$ appears with the same coefficient.","supporting_citations":[{"cited_title":"( 71) and ( 84) into Eq","cited_arxiv_id":null,"evidence_quote":"Supplies the method: two-point correlations of different-rank tensors at distinct points generate second-order nonlocal corrections, which this paper generalizes to spin."},{"cited_title":"This discrepancy arises from fluid expansion or compr ession and is mathematically expressed as Π = ⟨ˆp⟩ − p ( ǫ, na, Sαβ) = ⟨ˆp⟩l + ⟨ˆp⟩1 + ⟨ˆp⟩2 − p ( ǫ, na, Sαβ)","cited_arxiv_id":null,"evidence_quote":"The authors' earlier second-order spin hydrodynamics, built on Zubarev's operator, that this paper extends and whose equations gain the new terms."},{"cited_title":"( 71) into Eq","cited_arxiv_id":null,"evidence_quote":"Equilibrium density-matrix prediction of vorticity-induced hyperon polarization, the physical motivation for spin hydrodynamics."}],"review_version":1}