{"id":"413859b2-72ef-4a72-8e17-90eebf49aa21","arxiv_id":"2505.07324","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Using a proper-time expansion in the Color Glass Condensate, the authors compute the glasma's vorticity and local angular momentum along the beam axis and find they are qualitatively different, with the local angular momentum's quadrupole sign matching the measured Lambda polarization.","lead":"This paper computes the angular momentum of the glasma, the gluon field formed in the earliest instant of a heavy-ion collision. It finds that the local angular momentum along the beam axis behaves very differently from the fluid vorticity, and suggests it, not vorticity, drives the observed polarization of final hadrons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed exact cancellation of the Delta^4 term rests on Eq. (5.17), yet the paper's own nonzero vorticity implies Eq. (5.17) is not exact at the order used; without a direct check, Delta^6 dominance and the sign argument are unestablished.","rationale":"The reader's weakest assumption concerns early-time applicability; that concern is real and self-admitted in Sec. 6. The concern I raise is more direct and sits inside the derivation: the paper simultaneously relies on Eq. (5.17) being exact enough to kill the Delta^4 term and on a vorticity that is nonzero, but the same relation would force the vorticity to vanish. This is not a disagreement with the hydrodynamic consensus; it is a question of internal consistency at the order of the calculation. The check I propose is cheap: the Delta^4 coefficient is already implicit in the computed T^{0i}, and comparing it with the Delta^6 term at the box size used would settle whether L_z is genuinely dominated by third derivatives of the Poynting vector. If the check passes, the sign coincidence with Lambda polarization remains an interesting observation; if it fails, the central claim that local angular momentum, not vorticity, controls polarization loses its quantitative basis. I therefore do not think the paper should be accepted without this verification, but the concern is testable and not a reason to dismiss the framework. The reader's early-time critique and the gradient-expansion caveat remain additional conditions, but they are not the primary reason I would withhold full acceptance.","tokens_in":19141,"tokens_out":14050,"duration_ms":136687,"concrete_test":"Compute the Delta^4 coefficient C_4(r) = -tau/12 (partial_x T^{0y} - partial_y T^{0x}) directly from the eighth-order energy-momentum tensor on the grid used for Fig. 12, without imposing Eq. (5.17). If |C_4| Delta^4 is not negligible compared with the Delta^6 term at Delta = 0.2 fm, or if C_4 does not vanish to numerical precision, the claimed Delta^6 dominance fails. As a cross-check, compute omega_z from V = P/T^{00} and compare with Fig. 7; any nonzero vorticity from this definition quantifies the violation of Eq. (5.17) that also enters the Delta^4 coefficient.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 5E makes the central move: the order-Delta^4 term in L_z is 'identically zero' because the Poynting vector obeys the universal-flow relation T^{0i} = -(t/2) partial_i T^{00} (Eq. 5.17). This relation has a direct consequence that the paper does not discuss: with V = P/T^{00}, it gives V = -(t/2) grad ln T^{00}, whose curl, and therefore the vorticity omega_z of Eq. (5.1), is identically zero. Yet Sec. 5C reports a nonzero quadrupole vorticity. The only consistent reading is that Eq. (5.17) is exact only through seventh order in tau, while the numerics use the eighth-order energy-momentum tensor. If so, the Delta^4 coefficient of L_z, proportional to partial_x T^{0y} - partial_y T^{0x}, receives eighth-order corrections and is not identically zero. Since Delta = 0.2 fm is small, a nonzero Delta^4 term can dominate the claimed Delta^6 term. The numerical statement that L_z/Delta^6 is approximately independent of Delta is not accompanied by the range of Delta tested or an error estimate, so it does not currently resolve the issue. The qualitative difference between L_z and vorticity, and the sign coincidence with Lambda polarization, rest entirely on Delta^6 dominance.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies angular momentum of the glasma produced in ultrarelativistic heavy-ion collisions, using the authors' CGC-based proper-time expansion (up to eighth order) at τ=0.06 fm and mid-rapidity. It computes the global angular momentum perpendicular to the reaction plane and finds that only a small fraction of the initial nuclear angular momentum is transferred to the glasma. Its main subject is the local angular momentum along the beam axis, Lz, and the vorticity ωz of the transverse velocity field V = P/T^{00}. The authors find that Lz and ωz have qualitatively different spatial patterns: ωz is quadrupole-dominated with a sign opposite to the measured Λ polarization, while the quadrupole part of Lz at large impact parameter has the same sign as the measured Λ polarization. They explain the difference by the equation of universal flow T^{0i}=−(t/2)∂i T^{00} (Eq. (5.17)), which makes the order-Δ^4 contribution to Lz vanish, leaving a Δ^6 term built from third derivatives of the Poynting vector, whereas vorticity is a first derivative of the velocity field. They argue that local angular momentum, not thermal vorticity, may control the polarization of final-state hadrons.","tokens_in":19289,"tokens_out":12241,"duration_ms":108174,"significance":"If the Δ^6-dominance claim is correct, the paper establishes a qualitative distinction between vorticity and local angular momentum in the glasma that is not present in the naive rigid-body argument, and it offers a fresh perspective on the 'spin sign puzzle.' The analysis is transparent in several respects: the proper-time expansion is carried to eighth order with comparisons of sixth- and eighth-order results for Ly; the Fourier decomposition (Sec. 5B) is explicit; and the authors provide error bars based on the gradient-expansion cutoff δ. The central phenomenological statement is falsifiable in the sense that the sign of the quadrupole local angular momentum at large impact parameter is compared directly with data. However, the key mechanism — the vanishing of the Δ^4 term — rests on a relation that is exact only up to seventh order in the proper-time expansion, while the numerical results use the eighth-order tensor, and the numerical verification of Δ^6 dominance is not shown. This gap is the main obstacle to accepting the paper's central claim.","major_comments":[{"comment":"The claim that the order-Δ^4 contribution to Lz is 'identically zero' is not established for the quantity actually computed. Equation (5.17) is stated to hold 'exactly ... up to seventh order' in the proper-time expansion, while the numerical Lz in Sec. 5D is computed from the eighth-order energy-momentum tensor. The Δ^4 coefficient is the curl ∂xT^{0y}−∂yT^{0x}; any eighth-order correction to T^{0i} contributes to it, so the cancellation is, at best, approximate in the presented calculation. The only numerical support is the sentence 'We have verified numerically that Lz(⃗ r)/Δ^6 is approximately independent of Δ,' with no data, no range of Δ, and no error estimate. Please provide a quantitative verification: e.g., a plot of Lz(⃗ r)/Δ^6 for several values of Δ (say 0.1, 0.15, 0.2, 0.25, 0.3 fm) and the separate magnitudes of the Δ^4 and Δ^6 coefficients. This check is load-bearing because the qualitative difference between Lz and vorticity, and the sign argument, rest entirely on Δ^6 dominance.","section":"Sec. 5E, Eq. (5.17)"},{"comment":"The use of Eq. (5.17) appears inconsistent with the nonzero vorticity reported in Sec. 5C. If Eq. (5.17) were exact for the T used to compute ωz, then V = −(t/2)∇ ln T^{00} and ωz would vanish identically. The resolution is that Eq. (5.17) holds only up to seventh order in the proper-time expansion; this caveat needs to be stated wherever Eq. (5.17) is invoked. Moreover, the size of the eighth-order violation of Eq. (5.17) sets the scale of the residual Δ^4 contribution to Lz. Please quantify ∂xT^{0y}−∂yT^{0x} relative to the third-derivative combinations in Eq. (5.12) to demonstrate that the residual Δ^4 term is negligible in the region and at the impact parameters used for the sign statement.","section":"Secs. 5C and 5E"},{"comment":"The phenomenological sign claim relies on the quadrupole component Φ_L^2 at b=6 fm, but the dominant contribution to Φ_L comes from large transverse distances where the gradient expansion is least trustworthy, as the authors acknowledge in Sec. 2 and in footnote 2. The error bars obtained by varying δ_max between 0.8 and 1.0 probe only the gradient-expansion cutoff; they do not test the radial weighting of the integral or the sensitivity to the proper-time truncation order. Please provide a radial-shell decomposition of Φ_L^2 or an equivalent sensitivity test (e.g., varying the upper limit Rmax in Eq. (5.15)) and state explicitly whether the sign of the quadrupole contribution is stable.","section":"Secs. 5D, 6, Figs. 12-14"},{"comment":"The abstract's statement that 'neither vorticity nor thermal vorticity but instead the local angular momentum controls the polarization of final-state hadrons' is stronger than what the body supports. Section 6 states 'Since our method only works at early times, it is unclear if our results have any relevance for a description of the later stages of the collision,' and Section 7 states that 'a quantitative approach which would allow one to relate the local angular momentum to polarization needs to be developed.' Please align the abstract with these caveats, e.g., by saying that the results 'suggest' or 'are consistent with' a role for local angular momentum, or by adding a concrete (even schematic) argument for how the early-time Lz pattern would survive to freeze-out.","section":"Abstract and Secs. 6-7"}],"minor_comments":[{"comment":"There are several typos: 'preceeds' (p. 2) should be 'precedes', 'determing' (Sec. 7) should be 'determining', 'emphasis' (Sec. 5D) should be 'emphasize', and 'local local' (Sec. 5E) should be 'local'.","section":"Throughout"},{"comment":"The symbol h_n appears without definition and the displayed formula for ωz seems to have unbalanced parentheses; please fix.","section":"Eqs. (5.8)-(5.9)"},{"comment":"The use of δ both for the gradient-expansion parameter (Sec. 2) and for the non-radial Fourier coefficient δ_n(r) (Sec. 5B) is confusing; consider renaming one of them.","section":"Sec. 2 and Sec. 5B"},{"comment":"The discussion around Fig. 5 states that 'the curves obtained at the fourth and eighth orders show rapid growth' and that 'the second and sixth order results suggest the saturation,' which is confusing because the preceding sentence says the sixth- and eighth-order results agree for τ < 0.06 fm; please clarify which orders are being compared and why the eighth-order curve departs from the sixth-order one.","section":"Sec. 4, Fig. 5"},{"comment":"The grid geometry is not fully specified: the transverse plane dimensions and the total number of boxes are not given.","section":"Sec. 5D"},{"comment":"The estimate ω ∼ 10^{−3} fm^{−1} relies on the assumption that ∫ d^2σ·ω is conserved; the authors flag this with 'If we assume,' but it would be helpful to state explicitly that this is an order-of-magnitude assumption with no derivation for the glasma stage.","section":"Sec. 6"}],"recommendation":"major_revision","confidential_remarks":"The paper's central mechanism depends on the universal-flow relation Eq. (5.17), which is taken from the authors' own Ref. [28] and is not re-derived here; a self-contained derivation or a detailed summary would strengthen the paper. If the numerical check of Δ^6 dominance is provided and passes, the paper is likely a solid contribution; without it, the main claim is unsubstantiated. The abstract overstates the phenomenological conclusion relative to the body. The manuscript is within the journal's scope and should be sent back for the requested verification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper deserves a serious referee. The genuinely new result is the local angular momentum Lz along the beam axis: they show that the order-Delta^4 term vanishes because the Poynting vector obeys the universal flow relation T^{0i} = -(t/2) partial_i T^{00}, so Lz starts at Delta^6 (third derivatives of Poynting) while vorticity starts at first derivatives. That explains why Lz and vorticity differ qualitatively, and why the quadrupole part of Lz at large impact parameter matches the measured Lambda polarization while the vorticity has the opposite sign. That is a clear, falsifiable claim, distinct from their earlier work on global Ly.\n\nThe paper is honest and careful in several respects. Ly is computed at sixth and eighth order with a visible convergence check; the Fourier decomposition of vorticity and Lz is pedagogically useful; and Section 6 explicitly states the proper-time expansion converges only to about 0.08 fm, all results are at 0.06 fm, and it is unclear whether the early-time glasma is relevant for the later stages that produce hadrons. They do not oversell the polarization link, and they call for a quantitative transport treatment.\n\nThe load-bearing soft spot is the cancellation of the Delta^4 term. Equation (5.17) is said to hold only up to seventh order in the proper-time expansion, while the numerical energy-momentum tensor is computed to eighth order. The text calls the Delta^4 term \"exactly zero,\" but that exactness belongs to the truncated theory. With Qs*tau ~ 0.6, the eighth-order terms are not negligible. If the eighth-order contribution to curl P is nonzero, the Delta^4 term is nonzero and, because Delta = 0.2 fm is small relative to the transverse gradient scale, it would typically dominate the Delta^6 term. The paper's own nonzero vorticity shows that Eq. (5.17) is not exact at the order used, so this is an order-mixing issue, not a philosophical one. The numerical check that Lz/Delta^6 is approximately independent of Delta is the evidence that matters, but the paper gives no range of Delta tested, no fit, and no error bar. The sign argument for polarization rests entirely on Delta^6 dominance.\n\nThe other soft spots are milder: the glasma-to-freeze-out connection is an order-of-magnitude estimate using vorticity conservation and a hadron-scale Delta, and the claim that local angular momentum rather than thermal vorticity controls polarization is an argument, not a derivation. Both are acknowledged in the text.\n\nWho should read this: people working on spin polarization in heavy-ion collisions and on early-time glasma dynamics. It should go to peer review, not be desk-rejected. The referee should ask for a direct scan of Lz as a function of Delta, with the Delta^4 and Delta^6 contributions separated at the numerical order used, before accepting the central cancellation.","headline":"Carrington and Mrowczynski argue the glasma's local angular momentum, not vorticity, should set the sign of Lambda polarization; the mechanism hinges on a cancellation that may not survive at the order used in the numerics.","tokens_in":20005,"tokens_out":3797,"would_cite":false,"duration_ms":36403,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","12.38.Mh"],"model":"deepseek-v4-flash","headline":"The paper argues that the beam-axis polarization of final-state hadrons is set by the glasma's local angular momentum, not by vorticity or thermal vorticity, and it identifies the exact cancellation that makes the two quantities differ.","keywords":["glasma","angular momentum","vorticity","spin polarization","heavy-ion collisions","proper time expansion","Poynting vector","Lambda polarization"],"falsifier":"Evolve the glasma energy-momentum tensor from $\\tau\\simeq 0.06$ fm to freeze-out with a matched transport or hydrodynamic model and check whether the quadrupole sign of $L_z$ at large impact parameter survives; if the sign flips or the final polarization tracks thermal vorticity instead of local angular momentum, the central claim is wrong. A purely experimental version is a measurement of the azimuthal angle dependence of beam-axis $\\Lambda$ polarization as a function of impact parameter, looking for the predicted octupole-to-quadrupole crossover around $b\\approx 4$ fm.","tokens_in":18735,"feed_emoji":"🌀","tokens_out":12690,"duration_ms":104462,"temperature":0.7,"pith_summary":"The paper studies the glasma, the dense gluon field created in the first instants of an ultrarelativistic heavy-ion collision, and asks where its angular momentum lives. It finds that only a small fraction of the incoming nuclei's large angular momentum is transferred to the glasma, so the idea of a rapidly rotating fireball does not apply. Its central claim is about the beam-axis component: the glasma's local angular momentum and its vorticity behave qualitatively differently, because the leading derivative term in the local angular momentum vanishes exactly. The surviving term, built from third derivatives of the Poynting vector, gives a quadrupole pattern at large impact parameter whose sign matches the measured Lambda polarization, while the vorticity has the opposite sign. For that reason the paper argues that local angular momentum, not thermal vorticity, drives the polarization of final-state hadrons.","feed_headline":"Local angular momentum, not vorticity, sets Lambda polarization","feed_subtitle":"A glasma calculation matches the measured Lambda-polarization sign through the sixth-order term in local angular momentum.","key_machinery":"The load-bearing identity is the equation of universal flow, $T^{0x}=-\\frac{1}{2}t\\,\\partial_x T^{00}$ and $T^{0y}=-\\frac{1}{2}t\\,\\partial_y T^{00}$ at mid-rapidity, which relates the Poynting vector to gradients of the energy density and is exactly satisfied by the proper-time-expanded energy-momentum tensor through seventh order. Plugging this identity into the expression for the local angular momentum per unit rapidity, $dL_z/d\\eta=-\\tau\\int d^2R\\,(R_y T^{0x}-R_x T^{0y})$, kills the $\\Delta^4$ derivative term that would otherwise mimic the vorticity, leaving the $\\Delta^6$ term built from third derivatives of the Poynting vector. Vorticity, by contrast, is a first derivative of the velocity field $V^i=T^{0i}/T^{00}$. The different derivative order is what decouples the two observables.","core_discovery":"The central discovery is that for the glasma at mid-rapidity the equation of universal flow, $T^{0x}=-\\frac{1}{2}t\\,\\partial_x T^{00}$ and $T^{0y}=-\\frac{1}{2}t\\,\\partial_y T^{00}$, forces the leading $\\Delta^4$ contribution to the local beam-axis angular momentum to vanish identically. As a result $L_z$ is controlled by the $\\Delta^6$ term, which involves third spatial derivatives of the Poynting vector, whereas the vorticity $\\omega_z$ is a first derivative of the velocity field. These two quantities therefore need not look alike, and in the computed glasma they do not: at small impact parameter $L_z$ is mostly octupole, at large impact parameter it develops a quadrupole component whose sign agrees with the measured $\\Lambda$ polarization along the beam, while the glasma vorticity shows a quadrupole of the opposite sign. The paper also finds that the global angular momentum perpendicular to the reaction plane carried by the glasma is much smaller than the angular momentum of the participants, implying that the initial angular momentum stays with the valence quarks.","pith_inferences":["The paper's identity suggests a testable generalization: in any early-time model that satisfies the equation of universal flow, the local angular momentum should be dominated by third derivatives of the Poynting vector, so the same octupole-to-quadrupole crossover with impact parameter should appear.","A quantitative bridge is missing: the paper notes that no calculation connects the early-time $L_z$ to freeze-out polarization, so one could match a transport or hydrodynamic code at $\\tau\\simeq 0.06$ fm and evolve the sign pattern to see if it survives.","If the quadrupole $L_z$ is the correct driver, a clean experimental discriminator is the impact-parameter and multiplicity dependence of the azimuthal pattern of $\\Lambda$ polarization: it should cross from octupole-dominated at small $b$ to quadrupole-dominated at large $b$.","The $\\Delta^4$ cancellation may also explain why hydrodynamic models need ad hoc shear corrections to reproduce the sign: first-derivative quantities such as thermal vorticity are the wrong variable in a far-from-equilibrium system."],"forward_implications":["Spin-hydrodynamic descriptions that compute hadron polarization from thermal vorticity would need to be replaced or supplemented by a mechanism based on local angular momentum; the sign mismatch in the data is the paper's motivation.","The glasma phase is not a rapidly rotating fireball: the global angular momentum imparted to it is only a small fraction of the participants' angular momentum, consistent with vanishing global polarization at LHC energies.","The spin sign puzzle is resolved in sign at the glasma stage: the quadrupole part of $L_z$ at large impact parameter matches the measured $\\Lambda$ polarization, while the vorticity has the opposite sign.","Because the leading derivative term in $L_z$ vanishes by universal flow, any early-time system with a boost-invariant, mostly diagonal energy-momentum tensor will show the same decoupling of vorticity from local angular momentum."],"supporting_citations":[{"why":"Derives and verifies the equation of universal flow, Eq. (5.17), whose exact validity through seventh order makes the $\\Delta^4$ term vanish.","marker":"[28]"},{"why":"Introduces the equation of universal flow relating the Poynting vector to the energy-density gradient.","marker":"[41]"},{"why":"Provides the proper-time-expansion computation of the glasma energy-momentum tensor and the angular momentum formula reused here.","marker":"[23]"},{"why":"Extends the proper-time expansion and flow analysis that the paper relies on for the glasma velocity field.","marker":"[27]"},{"why":"Introduces the proper-time expansion of Yang-Mills fields used to solve the glasma dynamics.","marker":"[30]"},{"why":"Supplies the Glasma Graph approximation used to average products of gauge potentials.","marker":"[36]"},{"why":"Conjectured that beam-axis polarization is generated by elliptic-flow vorticity, the argument the paper tests against glasma.","marker":"[7]"},{"why":"Reports the observed Lambda polarization along the beam whose sign the quadrupole local angular momentum matches.","marker":"[9]"},{"why":"Gives the hydrodynamic thermal vorticity whose sign agrees with glasma vorticity but not with the measurement.","marker":"[16]"}],"fun_headline_variants":["Local angular momentum, not vorticity, drives Lambda polarization","Glasma: L_z, not omega_z, controls Lambda spin","Sixth-order term explains Lambda polarization sign in glasma","Lambda polarization traced to glasma local L_z, not vorticity","Glasma local momentum triumphs over vorticity for Lambda"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation is done at a single early instant, $\\tau=0.06$ fm, with a proper-time expansion valid only below about 0.08 fm, and the paper's own conclusion is that it is unclear whether the early-time glasma is representative of the later system that produces the measured hadrons.","fun_headline_variants_meta":{"raw":{"variants":["Local angular momentum, not vorticity, drives Lambda polarization","Glasma: L_z, not omega_z, controls Lambda spin","Sixth-order term explains Lambda polarization sign in glasma","Lambda polarization traced to glasma local L_z, not vorticity","Glasma local momentum triumphs over vorticity for Lambda"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1457,"prompt_tokens":934,"completion_tokens":523,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":437}},"tokens_in":550,"tokens_out":523,"duration_ms":4724,"temperature":1.0,"reasoning_tokens":437,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:19:13.728051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evolve the glasma energy-momentum tensor from $\\tau\\simeq 0.06$ fm to freeze-out with a matched transport or hydrodynamic model and check whether the quadrupole sign of $L_z$ at large impact parameter survives; if the sign flips or the final polarization tracks thermal vorticity instead of local angular momentum, the central claim is wrong. A purely experimental version is a measurement of the azimuthal angle dependence of beam-axis $\\Lambda$ polarization as a function of impact parameter, looking for the predicted octupole-to-quadrupole crossover around $b\\approx 4$ fm.","supporting_citations":[{"cited_title":"Gelis, E","cited_arxiv_id":null,"evidence_quote":"Derives and verifies the equation of universal flow, Eq. (5.17), whose exact validity through seventh order makes the $\\Delta^4$ term vanish."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the equation of universal flow relating the Poynting vector to the energy-density gradient."},{"cited_title":"Becattini, M","cited_arxiv_id":null,"evidence_quote":"Provides the proper-time-expansion computation of the glasma energy-momentum tensor and the angular momentum formula reused here."},{"cited_title":"Contrary to expectation, Lz(⃗ r) and ωz(⃗ r) are qualitatively different","cited_arxiv_id":null,"evidence_quote":"Conjectured that beam-axis polarization is generated by elliptic-flow vorticity, the argument the paper tests against glasma."},{"cited_title":"Adam et al","cited_arxiv_id":null,"evidence_quote":"Gives the hydrodynamic thermal vorticity whose sign agrees with glasma vorticity but not with the measurement."}],"review_version":1}