{"id":"19f23959-b16f-4eb7-a886-afb430b7a80d","arxiv_id":"2507.18761","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Lambda hyperon spin polarization on an iso-energy density freezeout surface is derived in linear response theory, with corrections to isothermal freezeout controlled by a small thermodynamic coefficient C.","lead":"Using a local equilibrium density operator, the authors derive spin polarization formulas that account for the freezeout surface being at constant energy density instead of constant temperature. They show the resulting corrections to the standard isothermal treatment are small, at most about 10% even at low collision energies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sec. III.A's spacetime-gradient heuristic is uncontrolled, and the ≤10% bound is inferred from |C|, not from computed polarization integrals; both need a numerical freezeout test before the central claim is accepted.","rationale":"The paper's central claim has two components: the derivation of the iso-energy freezeout polarization formulas (30) and the quantitative statement that corrections to isothermal freezeout are at most about 10% down to 11.5 GeV. The mathematical structure of the derivation is internally coherent and reduces properly to the isothermal limit when C → 0, so I do not see an algebraic defect that would justify rejection. The weakest point is the heuristic in Sec. III.A: approximating the larger of δT and δμ by a first-order spacetime Taylor expansion and fixing the smaller through Eq. (23). The error in this replacement is not bounded by C; it is a second-order gradient correction that can be comparable to the retained C-suppressed terms. In addition, the 10% statement is not a computed polarization but an extrapolation from the magnitude of C. The actual correction terms in Eqs. (30) involve ratios of chemical-potential gradients to kinematic vorticity and shear, and the ζC factor in the spin-Hall term can exceed |C| when μ/T is large. These concerns are addressable with a concrete hydrodynamic test, so the appropriate verdict remains conditional rather than accept or reject. The reader's weakest-assumption identification matches this concern, with the additional quantitative-bound issue reinforcing rather than replacing it.","tokens_in":13412,"tokens_out":11581,"duration_ms":130070,"concrete_test":"Run BEShydro with EOS3 at sqrt(sNN) = 11.5 GeV and extract both freezeout surfaces from the same hydrodynamic output: the isothermal surface T = T_D and the iso-energy surface e = e0 with e0 matched to e(T_D, μ̄) at the average freezeout chemical potential. Compute the Λ mean spin vector with the isothermal formula (17) on the first surface and with the iso-energy formula (30) on the second surface. Compare the pT-integrated and first-harmonic components. If any component changes by more than 10% of the isothermal value, the abstract's quantitative claim fails; if all changes are below 10%, the claim is supported. This check also indirectly probes the Sec. III.A heuristic, because large residual differences between the exact T and μ variations and their Taylor/constrained approximations on the freezeout surface would manifest as discrepancies between the two calculations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of Eqs. (30) rests on the Sec. III.A prescription: for |C| < 1, replace the exact hypersurface variation δμ by (y−x)^α ∂_α μ and set δT = −C (y−x)^α ∂_α μ via Eq. (23). This step is not error-controlled. Eq. (23) holds only to first order in the differences between the two surface points; when C is small, the retained δT term is O(C ∂μ), while the neglected quadratic terms are O(∂²) and can be comparable to or larger than the retained correction. The paper's claim that the approximation error is 'suppressed by C' compares with the alternative of using (y−x)·∂T, but it does not bound the residual error in Eq. (26). The C-dependent terms in Eqs. (30) may therefore be of the same order as the truncation error, so the geometric correction is not rigorously established. Separately, the abstract's ≤10% bound is not computed from Eq. (30). The correction terms are C β_[ρ ∂_σ] μ and ζC ∂_σ μ, whose size relative to the isothermal terms is set by C times ratios such as (β ∂μ)/ω, plus ζC in the spin-Hall channel; it is not simply |C|. Even with |C| ≤ 0.1, the ratio (β ∂μ)/ω can be large in regions that dominate local or harmonic polarization, and the paper itself concedes that local observables may be affected. Thus the quantitative central claim is an inference from the thermodynamic coefficient, not a computed result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a linear-response formalism for Lambda spin polarization at a constant-energy-density freezeout hypersurface. Starting from the local equilibrium density operator, the authors relate variations of temperature and chemical potential along the iso-energy surface through the thermodynamic coefficient C = (d e/d mu)|_T / (d e/d T)|_mu. For |C| < 1 they propose replacing the chemical-potential variation by its spacetime Taylor expansion and determining the temperature variation from delta T = -C delta mu, obtaining modified polarization formulas in Eqs. (30) that reduce to the isothermal expressions when C tends to zero. They evaluate C for two equations of state, EOS3 and HRG, find |C| <~ 0.1 for freezeout conditions down to sqrt(s_NN) = 11.5 GeV, and conclude that corrections to isothermal freezeout are at most 10%, although they may matter for local polarization observables.","tokens_in":13761,"tokens_out":6728,"duration_ms":73439,"significance":"If correct, the formulas provide a simple, equation-of-state-dependent improvement over the isothermal freezeout prescription for the Beam Energy Scan, and the paper is careful in separating the exact hypersurface relation from the subsequent Taylor expansion. The explicit reduction to prior results and the evaluation of C for two equations of state are useful. However, the quantitative 10% claim is not derived from the polarization integrals themselves, and the approximation underlying Eqs. (30) lacks a controlled error estimate, so the significance is conditional on further numerical validation.","major_comments":[{"comment":"The approximation step is uncontrolled. Equation (23) is derived from delta e = 0 and therefore holds only for displacements tangent to the iso-energy surface, while Eq. (26) promotes it to a relation between full spacetime Taylor coefficients of mu and T, including normal components. Moreover, replacing the hypersurface variation delta mu by (y-x)^alpha d_alpha mu leaves an error of order (y-x)^2 d^2 mu that is not proportional to C and can be comparable to the retained C-dependent terms in Eqs. (30). The statement that the error is suppressed by C compares delta T with (y-x) dot dT, but it does not bound the residual relative to the correction terms. A numerical test on a hydrodynamical freezeout surface, or at least an estimate of the neglected second derivatives, is required to validate Eqs. (30).","section":"Sec. III.A, Eqs. (23)-(26)"},{"comment":"The claim that corrections to isothermal freezeout are at most 10% is inferred from the smallness of C, but it is not computed from Eq. (30). The correction terms in Eqs. (30) are C beta_[rho d_sigma] mu and zeta C d_sigma mu, so their relative size is set by C times ratios such as (beta d mu)/omega, plus zeta C in the spin-Hall channel; those ratios can be large in the regions that dominate local or harmonic polarization. The paper's own caveat that local observables may become relevant acknowledges this. To support the abstract's quantitative statement, the authors should evaluate the polarization with and without the geometric corrections on a realistic freezeout surface, or at least estimate the relevant gradient ratios.","section":"Abstract and Sec. IV, Fig. 1"}],"minor_comments":[{"comment":"The caption contains a duplicated sentence: \"Lines of constant chemical potential are also shown:\" appears twice.","section":"Fig. 1 caption"},{"comment":"In the paper organization sentence, \"the local equilibrium density operation\" should presumably be \"the local equilibrium density operator.\"","section":"Sec. I, second-to-last paragraph"},{"comment":"The text refers to the \"C > 1 regime\" where the earlier discussion used |C| > 1; the absolute value should be kept for consistency, since C can be negative.","section":"Sec. III.A, after Eq. (28)"},{"comment":"There is an index mismatch in Eq. (28a): the left-hand side is written as delta beta^mu while the right-hand side is expressed in terms of nu indices; this should be corrected for readability.","section":"Eq. (28)"},{"comment":"The abstract states that corrections are \"at most a 10% effect,\" while Sec. V says they are \"typically below 10%\"; these statements are not identical and should be reconciled.","section":"Sec. V vs Abstract"},{"comment":"The bracket notation in Eqs. (32) is difficult to parse because the antisymmetrization brackets are opened on u and closed on different gradient terms; using explicit parentheses or a displayed symmetrization convention would improve clarity.","section":"Eqs. (32a)-(32c)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the derivation is presented honestly. My main concern is that the headline quantitative claim of at most 10% is presented as a result even though it is an inference from the smallness of C, and the authors themselves defer a full hydrodynamic implementation. I would recommend that the published version either contain the numerical test or state the bound as a plausibility estimate rather than a quantitative prediction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: Palermo and Shokri extend the linear-response spin polarization formalism to freezeout on an iso-energy density hypersurface, and they do it cleanly. The new piece is a thermodynamic coefficient C that relates tangent variations of T and μ on the surface, and formulas (30) that reduce to the known isothermal ones when C→0. They also evaluate C for two equations of state and find |C| ≲ 0.1 even at √sNN=11.5 GeV. That part is solid and worth having: it gives a compact criterion for when isothermal freezeout is a good approximation.\n\nThe soft spot is the step from |C| small to \"corrections at most 10%\". The paper does not compute the polarization from Eq. (30); it infers the size of the correction from the size of C alone. That inference ignores the ratios that actually enter the integrand: the extra terms are C β_[ρ ∂σ] μ and ζ C ∂σ μ, which compare to the leading vorticity/shear terms as C times (β∂μ)/ω, plus ζC in the spin-Hall channel. In regions where ω is small, that ratio can be large even if C is 0.1. The abstract's \"at most 10%\" is therefore stronger than what the derivation supports.\n\nThere is also a more technical concern about the Sec. III.A prescription itself. The paper replaces the exact hypersurface variation δμ with a spacetime Taylor term and then sets δT = -C(y-x)·∂μ. That is a first-order-in-coordinate-difference statement, and the neglected quadratic terms are not bounded relative to the retained C term. The authors are aware that hypersurface and spacetime variations differ; their argument that the error is \"suppressed by C\" compares the C term to the naive (y-x)·∂T alternative, not to the truncation error. So the geometric correction is plausible but not rigorously error-controlled. A numerical freezeout test with an actual hydrodynamic surface would settle both issues: compute (30) against a full evaluation of (12) with exact variations.\n\nNone of this kills the paper. The formalism is coherent, the limits C→0 and C→∞ behave as expected, and the multi-chemical-potential extension is a nice addition. The result is a legitimate extension of Becattini et al. and is not circular. But the headline quantitative claim is conditional on a numerical test that the paper explicitly leaves to future work. For a referee: this deserves serious review, not a desk reject, because the framework is useful and the flaw is fixable. I would probably cite it for the formulas and the C criterion, with a caveat about the 10% bound.","headline":"A clean generalization of isothermal freezeout to constant energy density with a useful coefficient C, but the headline 10% bound is an inference from C, not a computed polarization.","tokens_in":14232,"tokens_out":2296,"would_cite":true,"duration_ms":23763,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.Nq"],"model":"deepseek-v4-flash","headline":"The paper derives Lambda spin-polarization formulas for constant-energy-density freezeout and shows the corrections to isothermal freezeout stay below ten percent at RHIC energies.","keywords":["Lambda hyperon spin polarization","heavy-ion collisions","freezeout hypersurface","iso-energy density","local equilibrium density operator","linear response theory","thermodynamic coefficient","RHIC beam-energy scan"],"falsifier":"Take a computed freezeout surface from a 3+1D hydrodynamic simulation at $\\sqrt{s_{NN}}=11.5$ GeV, evaluate the exact tangent variations of $T$ and $\\mu$ on the surface, and compare them with $\\delta T=-C\\,\\delta\\mu$ using spacetime gradients; if the mismatch is comparable to the gradient size rather than suppressed by $C$, the central formulas fail.","tokens_in":13207,"feed_emoji":"🌀","tokens_out":10181,"duration_ms":100191,"temperature":0.7,"pith_summary":"At freezeout, hydrodynamic simulations typically stop on a surface of constant energy density, while spin-polarization formulas have been built assuming constant temperature. This paper derives what changes for Lambda hyperon spin polarization when the freezeout surface is genuinely iso-energy, using the local equilibrium density operator and keeping the surface geometry explicit before expanding in gradients. The entire correction is controlled by one thermodynamic coefficient, $C = (\\partial e/\\partial\\mu)_T/(\\partial e/\\partial T)_\\mu$, and the paper finds that $C$ stays below 0.1 for the equations of state and freezeout conditions relevant to RHIC, down to $\\sqrt{s_{NN}}=11.5$ GeV. The finite-density corrections to the isothermal formulas are therefore at most about ten percent, but they can still matter for local polarization observables and for equations of state with a critical point.","feed_headline":"Constant-energy freezeout shifts Lambda spin by at most 10%","feed_subtitle":"The isothermal freezeout approximation for Lambda polarization holds to about ten percent down to RHIC's lowest energy.","key_machinery":"The central object is the thermodynamic coefficient $C=(\\partial e/\\partial\\mu)_T/(\\partial e/\\partial T)_\\mu$, which measures how much temperature must change to compensate a chemical-potential change when energy density is held fixed. The argument combines this coefficient with the surface identity $\\delta T=-C\\,\\delta\\mu$ and a guiding rule: on an iso-energy-density surface, approximate the larger of the two variations by its spacetime Taylor expansion and determine the smaller one from the constraint. Inserting these geometric variations into the linear-response expansion of the local equilibrium density operator yields the polarization formulas, and taking $C\\to 0$ returns the isothermal results.","core_discovery":"The paper's central claim is that for freezeout on an iso-energy-density hypersurface, the surface constraint forces the temperature variation and chemical-potential variation to be linked by $\\delta T = -C\\,\\delta\\mu$, where $C$ is the ratio of the energy-density derivatives. In the relevant regime $|C|<1$, the chemical-potential variation should be approximated by its spacetime gradient and the temperature variation eliminated through the constraint; this produces modified spin-vector formulas for thermal vorticity, thermal shear, and the spin-Hall effect, shown in Eqs. (30). The formulas reduce exactly to the isothermal freezeout result when $C\\to 0$. Evaluating $C$ with the EOS3 and hadron-resonance-gas equations of state, the paper finds $|C|\\lesssim 0.1$ at freezeout conditions corresponding to RHIC energies from 200 GeV down to 11.5 GeV, so the isothermal approximation is correct to about ten percent.","pith_inferences":["If a future equation of state with a critical point makes $|C|$ grow, the same formulas predict where and when isothermal freezeout breaks down; evaluating $C$ along the freezeout curve of a critical-point EoS is a direct test of that sensitivity.","The geometry-informed substitution is not limited to spin: any observable computed with the local equilibrium density operator on a constant-energy-density surface should receive analogous $O(C)$ corrections, so spectra and flow harmonics are natural places to look for the same effect.","Comparing Eq. (30) with the isothermal formulas on identical hydrodynamic freezeout surfaces would isolate whether the predicted ten-percent shifts are visible in current local-polarization measurements or hidden by their statistical uncertainties."],"forward_implications":["At RHIC beam energies from 200 GeV down to 11.5 GeV, $|C|$ remains below about 0.1 for both EOS3 and HRG, so the $|C|<1$ formulas in Eq. (30) are the relevant ones.","The isothermal freezeout approximation for Lambda polarization is safe at the ten-percent level across the RHIC energy scan, so sizable deviations in global polarization data cannot be blamed on finite-density freezeout geometry.","Finite-density corrections can still shift local polarization observables in low-energy collisions, where the relative orientation of gradients changes from cell to cell.","The prescription extends to multiple conserved charges by choosing the intensive variable with the largest $|\\partial e/\\partial X|$ as the dependent one, which fixes the form of the corrections when more than one chemical potential is present.","The $C>1$ regime, a near iso-chemical-potential freezeout, is not reached by the two equations of state considered."],"supporting_citations":[{"why":"Defines the local equilibrium density operator whose linear-response expansion is the starting point of the calculation.","marker":"[17–19]"},{"why":"Provides the earlier linear-response derivation of vorticity- and shear-induced spin polarization used as the intermediate formulas.","marker":"[28]"},{"why":"Gives the isothermal freezeout polarization formulas that are the baseline and the C → 0 limit of the new expressions.","marker":"[29]"},{"why":"Establishes the spin-Hall contribution with a chemical-potential gradient that enters the isothermal result.","marker":"[36]"},{"why":"Supplies the EOS3 equation of state used to evaluate the coefficient C at finite chemical potential.","marker":"[39, 40]"},{"why":"Supplies the hadron resonance gas equation of state for the second evaluation of C.","marker":"[41]"},{"why":"Provides freezeout values of μ/T at RHIC energies that set the scale for the ten-percent estimate.","marker":"[38]"}],"fun_headline_variants":["Constant-energy freezeout keeps Lambda spin shift under 10%","Isothermal freezeout for Lambda spin valid to 10%","Spin polarization: constant-energy freezeout close to isothermal","Lambda spin correction from freezeout geometry stays below 10%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on trusting that the larger change—in temperature or in the baryon chemical potential—along the freezeout surface is well described by its local spacetime gradient, with the smaller change then fixed by the constant-energy-density relation. If the true surface variation is not captured by that gradient, the smallness of $C$ does not by itself keep the formulas accurate.","fun_headline_variants_meta":{"raw":{"variants":["Constant-energy freezeout keeps Lambda spin shift under 10%","Isothermal freezeout for Lambda spin valid to 10%","Spin polarization: constant-energy freezeout close to isothermal","Lambda spin correction from freezeout geometry stays below 10%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000277,"raw_usage":{"total_tokens":1601,"prompt_tokens":845,"completion_tokens":756,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":686}},"tokens_in":461,"tokens_out":756,"duration_ms":7312,"temperature":1.0,"reasoning_tokens":686,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:08:33.913156+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a computed freezeout surface from a 3+1D hydrodynamic simulation at $\\sqrt{s_{NN}}=11.5$ GeV, evaluate the exact tangent variations of $T$ and $\\mu$ on the surface, and compare them with $\\delta T=-C\\,\\delta\\mu$ using spacetime gradients; if the mismatch is comparable to the gradient size rather than suppressed by $C$, the central formulas fail.","supporting_citations":[{"cited_title":"The hadron resonance gas model,","cited_arxiv_id":null,"evidence_quote":"Supplies the hadron resonance gas equation of state for the second evaluation of C."}],"review_version":2}