{"id":"c95926ba-f153-4836-b34e-7da48a2f3f50","arxiv_id":"2507.01454","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"RHIC isobar data are explained by different shapes of 96Ru and 96Zr, with 96Zr showing a large octupole deformation, so nuclear structure uncertainty, not the magnetic field, dominates the observed ratios.","lead":"A combined task force of nuclear structure and heavy-ion physicists compared model calculations with RHIC data on collisions of the isobars 96Ru and 96Zr. They conclude the observed flow differences are naturally explained by different intrinsic nuclear shapes, with 96Zr likely pear-shaped. This matters because it reframes the chiral magnetic effect search and proposes isobar collisions as a precision nuclear imaging tool.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The explanation of the central v3 excess hinges on the untested conversion of B(E3;3-→0+) into a static ground-state octupole deformation of 96Zr; if the ground state is instead an octupole vibrator, the comparison with STAR data loses its low-energy anchor.","rationale":"The reader's weakest-assumption identification is accurate and is independently flagged in the manuscript, so no new verdict is needed. The report's strongest claim is broad enough to survive the loss of the static-octupole conversion in its qualitative form (different radial profiles and beta2 already generate nuclear-structure backgrounds), but the specific v3/r3 evidence and the claim of qualitative agreement with low-energy structure rest on the conversion. I therefore keep the CONDITIONAL verdict rather than raising or lowering it. The proposed check is the missing calculation the manuscript itself identifies as necessary in Sec. 3.5.2 ('This assumption has to be further studied'), and it is concrete and feasible with the MCSM and PGCM tools already deployed in Sec. 3.5–3.6.","tokens_in":55712,"tokens_out":6663,"duration_ms":80082,"concrete_test":"Perform a configuration-mixing calculation (MCSM or PGCM) for 96Zr using one Hamiltonian that reproduces the 0+ and 3- energies, compute both B(E3;3-→0+) and the intrinsic beta3 of the 0+ ground-state density separately, without assuming the two states share an intrinsic state. If the ground-state beta3 is ≲0.1 while B(E3) remains ≈42–52 W.u., the static-octupole conversion in Sec. 3.5.2 is invalid; then rerun the TRENTo/hydrodynamic v3 ratio with beta3=0 for 96Zr and check whether the central v3 excess persists or disappears. This distinguishes the static-deformation and octupole-vibration interpretations and directly settles whether the v3 evidence supports the report's central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative v3 ratio explanation rests on the assumption, flagged in Sec. 3.5.2, that the 0+ ground state and the 3- state of 96Zr share a common intrinsic state with static octupole deformation beta3 ≈ 0.2–0.27. Without that conversion, the measured B(E3;3-→0+) = 42–53 W.u. does not imply an octupole-deformed ground state; Sec. 3.2 notes that no alternating-parity band is observed and that there is no simple experimental procedure to test whether the two states have equal deformations. The mean-field and ab initio results in Sec. 3.2–3.6 give ground-state beta3 values of 0 (most Skyrme/Gogny/covariant), 0.10 (SV-mas07), 0.125 (BSkG), or 0.175 (exploratory PGCM) — all below the 0.2–0.27 used in Table 3 and in the hydrodynamic/AMPT comparisons. Some of the simulation inputs are themselves derived from STAR data (Ref. [99]), so the reproduction of the central v3 excess is partly circular. If the ground state is an octupole vibrator with negligible static beta3, the large initial-state triangularity attributed to 96Zr disappears, and the ~10% central v3 excess has no demonstrated explanation within the mechanisms tested here (diffuseness, beta2, short-range correlations, free-streaming all fail to produce it). The qualitative claim that isobar data are sensitive to nuclear structure survives, but the specific claim that the differences 'qualitatively agree' with low-energy information is contingent on this unresolved conversion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper is the report of the EMMI Rapid Reaction Task Force on isobar collisions at RHIC. It combines low-energy nuclear structure calculations (mean-field, beyond-mean-field, shell model, ab initio PGCM, lattice EFT prospects) with high-energy initial-state and full dynamical simulations (T RENTo, SMASH, AMPT, iEBE-VISHNU, Trajectum) to understand the STAR measurements of 96Ru+96Ru and 96Zr+96Zr collisions. The central claim is that the measured isobar ratios are naturally explained if the two nuclei have different intrinsic shapes, in particular a large octupole deformation in 96Zr, and that the corresponding differences qualitatively agree with low-energy nuclear information. A corollary is that nuclear-structure effects mask relative magnetic-field-driven variations, so a sound conclusion about the chiral magnetic effect requires proper uncertainty quantification of nuclear ground-state properties. The report also recommends using the 136Xe-136Ce pair for future CME searches and argues that isobar collisions can serve as a precision probe of nuclear shapes and skins.","tokens_in":56148,"tokens_out":4091,"duration_ms":51723,"significance":"If the central claim holds, this report establishes a genuinely cross-disciplinary result: high-energy isobar collision ratios become a quantitative imaging tool for nuclear ground-state shapes and surface properties, and CME conclusions become contingent on nuclear-structure uncertainties. The report is strengthened by the diversity of independent model comparisons, by its explicit engagement with low-energy experimental data, and by the unusually candid self-assessment of limitations, including the unresolved static-vs-vibrational octupole interpretation and the absence of uncertainty quantification. The Taylor-expansion approach and the correlated-sampling method of Sec. 4.6 are practically valuable methodological contributions. The significance is high for both the heavy-ion and nuclear-structure communities, provided the load-bearing octupole assumption is settled or explicitly downgraded.","major_comments":[{"comment":"The quantitative explanation of the v3 excess rests on converting B(E3; 3- -> 0+) = 42-53 W.u. into a static ground-state octupole deformation beta3 ~ 0.2-0.27, but the report itself states in Sec. 3.5.2 that this assumes the ground state and the 3-1 state share the same intrinsic state with static octupole deformation, and Sec. 3.1.2 notes that no alternating-parity band built on the ground state is observed. The independent many-body results in Tables 1-2 and in Secs. 3.4 and 3.6 give ground-state beta3 values of 0, 0.05-0.125, or at most 0.175, all below the 0.2-0.27 used in Table 3 and in the hydrodynamic/AMPT comparisons. Because the central v3 excess in Secs. 4 and 5 is generated by the large beta3(96Zr) input, the conclusion that the STAR observations are naturally explained by low-energy structure is contingent on an assumption that the report itself identifies as unresolved.","section":"Sec. 3.5.2 and Sec. 3.2"},{"comment":"There is partial circularity in the validation chain: beta3(96Zr) = 0.2 is taken from Ref. [99], a value inferred from STAR isobar data, and is then used as input to simulations that reproduce those same STAR data. This applies to the Trajectum case-5 setup, the T RENTo default/cases 3-6, and the AMPT default. As a validation of the nuclear-structure explanation, the agreement with STAR data in these sections is therefore not independent. The paper should either repeat the comparisons with beta3 values derived solely from low-energy data (including the smaller values predicted by the calculations in Sec. 3), or explicitly demonstrate that the reproduced ratios are insensitive to the provenance of beta3; without this, the 'naturally explained' claim is weakened.","section":"Sec. 5.5 (Table 9), Sec. 4.2 (Cases 3-6), Sec. 5.3"},{"comment":"No uncertainty bands are attached to the theoretical isobar ratios, despite the experimental precision of about 0.4% and the report's own recommendation in Sec. 6 that quantitative conclusions about the magnetic-field effects require uncertainty quantification. Propagated uncertainties on beta2, beta3, a, R0, dmin, and the transport parameters would be needed to assess whether the agreement with STAR data is statistically supported or merely a visual match. The final conclusions should either include such uncertainty quantification or explicitly downgrade the 'naturally explained' statement to a plausibility argument.","section":"Secs. 4-6"}],"minor_comments":[{"comment":"In the sentence 'the largest effect is observed when we move from case 4 to vase 5', 'vase' should read 'case'.","section":"Sec. 4.3.3"},{"comment":"The distinction between the Woods-Saxon deformation parameter beta_WS used in the simulation tables and the multipole moment beta_l0 of Eq. (20) should be applied consistently; for the values used here the difference is about 0.02 and should be accounted for when comparing high-energy extractions with low-energy results.","section":"Sec. 3.3.1 and Tables 3, 5, 9"},{"comment":"The row 'Case 1 and 5 difference' lists Delta_beta2^2 and Delta_beta3^2 alongside Delta R0 and Delta a, but Eq. (81) uses Delta beta_n^2 notation; the table would be clearer if it explicitly stated that the tabulated deformation differences are already squared.","section":"Sec. 5.3.2, Table 7"},{"comment":"The quantity N^off_trk,raw and the meaning of the 'off' superscript are not defined in the text; a brief definition would help readers interpret the inverse-correction procedure.","section":"Sec. 5.5.2 and Fig. 37"}],"recommendation":"major_revision","confidential_remarks":"This is a strong and unusually open collaborative report, but the central bullet is somewhat stronger than the evidence presented in the body. The main issue is not the quality of the computations but the reliance on an unresolved conversion from B(E3) to a static ground-state octupole deformation, together with the partial circularity of using a STAR-derived beta3 as simulation input. I would ask the editor to require either a direct test of the octupole-vibrator alternative or a clear downgrade of the 'naturally explained' conclusion. The paper's length is also a burden; a short executive summary of the load-bearing assumptions and their status would improve accessibility without changing the scientific content."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nYou should know this for what it is: a genuinely useful task-force report, and the strongest statement so far that the RHIC isobar ratios are dominated by nuclear structure rather than the chiral magnetic effect. The central qualitative claim holds up well. No matter which way you slice the initial-state or hydrodynamic models, the observed ratios are naturally explained by different intrinsic shapes, with a larger quadrupole deformation in 96Ru and a larger skin in 96Zr.\n\nWhat is actually new here is the breadth. The report ships a large set of original calculations: a 17-parameterization Skyrme survey, BSkG searches over the full gamma range with finite octupole deformation, new Gogny PGCM and Monte Carlo shell model results, plus systematic TRENTo, free-streaming, SMASH, AMPT, and Trajectum scans. The Taylor-expansion framework in Sec 5.3.2 is a clean way to organize which nuclear parameter drives which observable. The report also deserves credit for its honesty: the authors repeatedly flag preliminary status, missing uncertainty quantification, and their own load-bearing assumption.\n\nThat load-bearing assumption is the soft spot, and it is the one you should worry about. The entire v3 ratio explanation rests on converting the measured B(E3; 3- -> 0+) in 96Zr into a static ground-state octupole deformation beta3 ~ 0.2-0.27. As the report itself says in Sec 3.5.2, this requires the 0+ ground state and the 3- state to share the same intrinsic octupole-deformed shape. No alternating-parity band is observed in 96Zr, most mean-field and ab initio results put the ground-state beta3 lower (0 to 0.175), and there is no simple experimental test of the assumption. If 96Zr is instead an octupole vibrator, the large initial-state triangularity driving the central v3 excess has no demonstrated explanation within the mechanisms the report tested.\n\nThe partial circularity is real but less damaging than it looks. Several simulation sets take beta3 = 0.2 from Ref. [99], which was itself inferred from STAR data, and then reproduce STAR data. The independent low-energy B(E3) evidence and the diversity of models mitigate this, but the report would be stronger if it clearly separated inputs derived from STAR from inputs derived from low-energy data.\n\nWho is this for? Anyone working on the RHIC isobar program, the CME search, or nuclear structure imaging from heavy-ion collisions. It is long, but as a reference document it is valuable. It absolutely deserves a serious referee; the quantitative claims need tightening, but the central conclusion that nuclear structure dominates the isobar ratios is well supported.\n\nMy recommendation: send it to peer review, with the expectation that the octupole assumption and the circularity get addressed head-on.","headline":"A comprehensive and mostly honest task-force report that makes the nuclear-structure interpretation of the isobar data credible, but rests its sharpest quantitative claim on a flagged and unresolved assumption about octupole deformation in 96Zr.","tokens_in":56983,"tokens_out":3402,"would_cite":true,"duration_ms":39104,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Isobar collision ratios are explained by nuclear shapes, not magnetic fields.","keywords":["isobar collisions","nuclear deformation","octupole deformation","chiral magnetic effect","collective flow","neutron skin","heavy-ion collisions","nuclear structure"],"falsifier":"A decisive test would be a model-independent measurement of zirconium-96's ground-state shape, for instance a precision Coulomb-excitation or laser-spectroscopy experiment that determines whether the ground state carries a static octupole deformation $\\beta_3 \\approx 0.2$ without assuming equality of the $0^+$ and $3^-$ wave functions. Alternatively, an ab initio calculation that yields a spherical or vibration-only octupole picture for the ground state, or a hydrodynamic simulation that reproduces the measured $v_3$ ratio without any static $\\beta_3$, would falsify the report's central explanation.","tokens_in":55507,"feed_emoji":"⚛️","tokens_out":5189,"duration_ms":57735,"temperature":0.7,"pith_summary":"This report argues that the measured deviations from unity in ratios of observables from collisions of the two A=96 isobars, ruthenium-96 and zirconium-96, are naturally explained by differences in the two nuclei's ground-state shapes and skin thicknesses. If that is right, the deviations are not primarily a magnetic-field effect: the nuclear-structure background masks any relative variation genuinely driven by the chiral magnetic effect, so a sound conclusion about that effect depends on knowing the ground states. The report further establishes isobar collision ratios as a quantitative imaging tool for nuclear ground-state properties, including quadrupole deformation, octupole deformation, and neutron skin. The central technical claim is that a large static octupole deformation in zirconium-96 is the source of the measured enhancement of triangular flow in its collisions.","feed_headline":"Isobar ratios expose nuclear shapes, mask magnetic field","feed_subtitle":"Flow measurements point to distinct quadrupole and octupole deformations, making CME claims hostage to nuclear structure.","key_machinery":"The load-bearing object is the deformed Woods-Saxon density of each colliding nucleus, $\\rho(r,\\theta,\\phi) = \\rho_0/(1+\\exp[(r-R(\\theta,\\phi))/a])$ with $R(\\theta,\\phi) = R_0(1+\\beta_2 Y^0_2 + \\beta_3 Y^0_3 + \\cdots)$, whose parameters $R_0$ (radius), $a$ (diffuseness or skin), $\\beta_2$ (quadrupole), $\\beta_3$ (octupole), and $\\gamma$ (triaxiality) are sampled event-by-event to build initial conditions. The argument works through the near-linear mapping from initial eccentricities ($\\varepsilon_2$, $\\varepsilon_3$) computed from these densities to final harmonic flow coefficients ($v_2$, $v_3$), combined with a Taylor-expansion relation for isobar ratios, $\\mathcal{O}_{\\mathrm{Ru}}/\\mathcal{O}_{\\mathrm{Zr}} \\approx 1 + c_1\\Delta\\beta_2^2 + c_2\\Delta\\beta_3^2 + c_3\\Delta a + c_4\\Delta R_0$, which isolates how each nuclear parameter shows up in each observable. Octupole deformation mainly feeds triangular flow in central collisions, while the diffuseness difference drives the non-monotonic multiplicity and flow-ratio trends.","core_discovery":"The paper's central claim is that the measured isobar ratios of multiplicity, elliptic flow, and triangular flow are reproduced by initial-state and hydrodynamic calculations once the two nuclei are given different intrinsic shapes: ruthenium-96 more quadrupole-deformed, zirconium-96 with a larger neutron skin and a significant octupole deformation. This picture is consistent with low-energy nuclear structure information, in particular a strong E3 transition in zirconium-96 that can be read as a static octupole deformation $\\beta_3 \\approx 0.2$–$0.27$. Consequently, the report concludes that the isobar run does not provide evidence for the chiral magnetic effect; the nuclear-structure background dominates the ratios and must be quantified before magnetic-field effects can be isolated. It also concludes that the same measurements, once understood, turn isobar collisions into a precision probe of nuclear shapes and skins.","pith_inferences":["Editorial: if the report's method is sound, isobar ratios become a new high-energy handle on ground-state properties of isotopes that are hard to study at low energy, including nuclei relevant to neutrinoless double beta decay.","Editorial: the octupole-deformation interpretation predicts specific centrality and multiplicity dependence of higher-order cumulant ratios, such as $v_3\\{4\\}/v_3\\{2\\}$, which future isobar data could test.","Editorial: the Taylor-expansion identity suggests a general experimental strategy: measuring ratios of observables across isobar pairs cancels much of the unknown bulk dynamics, isolating geometry parameters; this could be applied to other isobar pairs beyond A=96.","Editorial: the report leaves open whether the large correlation energy that beyond-mean-field methods attribute to octupole shapes is physical; if it is not, the value of $\\beta_3$ extracted from flow will need reinterpretation."],"forward_implications":["If the central claim is correct, the ratio of triangular flow between the isobars serves as a quantitative measure of octupole deformation in zirconium-96.","The measured elliptic-flow and multiplicity ratios can be used to extract the quadrupole deformation of ruthenium-96 and the difference in skin thickness between the two nuclei.","A sound conclusion about the chiral magnetic effect requires nuclear-structure uncertainty quantification; without it, isobar ratios cannot be read as magnetic-field signals.","Future isobar runs intended to probe magnetic-field effects should use a pair in which at least one nucleus is near-spherical, since the A=96 pair is particularly geometry-sensitive.","Isobar ratio data can also constrain the short-range nucleon-nucleon repulsion parameter in initial-state models, which is poorly determined by existing global analyses."],"supporting_citations":[{"why":"Supplies the measured isobar ratios of multiplicity and flow that the report sets out to explain.","marker":"[14]"},{"why":"Proposes the Woods-Saxon parameterization with octupole deformation in zirconium-96 that reproduces the triangular-flow ratio.","marker":"[99]"},{"why":"Establishes the mechanism connecting quadrupole and octupole deformations to the harmonic flow coefficients $v_2$ and $v_3$.","marker":"[57]"},{"why":"Derives the Taylor-expansion formula for isobar ratios in terms of differences in $\\beta_2$, $\\beta_3$, $a$, and $R_0$.","marker":"[144]"},{"why":"Provides a full hydrodynamic analysis showing that the measured multiplicity and flow ratios are reproduced by deformed nuclear shapes.","marker":"[123]"},{"why":"Gives an independent hydrodynamic simulation of isobar ratios using quadrupole and octupole deformations consistent with the report's picture.","marker":"[248]"},{"why":"Supplies the low-energy $B(E3; 3^-\\to 0^+)$ strength used as evidence for strong octupole collectivity in zirconium-96.","marker":"[30]"},{"why":"Provides the shell-model calculation from which the static octupole deformation $\\beta_3 = 0.27$ is inferred.","marker":"[115]"}],"fun_headline_variants":["Isobar collisions reveal shapes, hide magnetic effect","Nuclear shapes skew isobar data, not chiral magnetism","Deformed nuclei explain isobar ratios at RHIC","Zirconium's octupole shape alters flow ratios","Isobar run maps nuclear skin, dims CME signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that zirconium-96's ground state really is octupole-deformed in the static sense, inferred by assuming that the ground state and the $3^{{-}}$ excited state are formed from the same intrinsically deformed shape; if that conversion is wrong, the triangular-flow excess needs another source.","fun_headline_variants_meta":{"raw":{"variants":["Isobar collisions reveal shapes, hide magnetic effect","Nuclear shapes skew isobar data, not chiral magnetism","Deformed nuclei explain isobar ratios at RHIC","Zirconium's octupole shape alters flow ratios","Isobar run maps nuclear skin, dims CME signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000333,"raw_usage":{"total_tokens":1836,"prompt_tokens":918,"completion_tokens":918,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":837}},"tokens_in":534,"tokens_out":918,"duration_ms":8798,"temperature":1.0,"reasoning_tokens":837,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:50:59.442953+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a model-independent measurement of zirconium-96's ground-state shape, for instance a precision Coulomb-excitation or laser-spectroscopy experiment that determines whether the ground state carries a static octupole deformation $\\beta_3 \\approx 0.2$ without assuming equality of the $0^+$ and $3^-$ wave functions. Alternatively, an ab initio calculation that yields a spherical or vibration-only octupole picture for the ground state, or a hydrodynamic simulation that reproduces the measured $v_3$ ratio without any static $\\beta_3$, would falsify the report's central explanation.","supporting_citations":[],"review_version":1}