{"id":"84f2cae1-6538-49aa-a534-667291116780","arxiv_id":"1908.10231","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Including the neutron skin of zirconium in isobar collision simulations halves the predicted Ru/Zr magnetic-field strength difference from 10% to 5% in peripheral events.","lead":"This paper simulates Ru+Ru and Zr+Zr collisions at RHIC with a transport model, adding a measured neutron-skin halo to zirconium and deformation to ruthenium. It finds the expected difference in magnetic field strength between the two isobars shrinks from 10% to 5% in peripheral collisions, which would shrink the predicted signal for the chiral magnetic effect relative to background.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-of-2 reduction in the magnetic-field ratio is computed at the upper limit of the measured neutron skin and would likely weaken at the central value; the quantitative claim needs an uncertainty propagation before it can be taken at face value.","rationale":"The reader's weakest assumption correctly identifies the most load-bearing point: the quantitative headline depends on evaluating the neutron skin at the upper limit of the measured constraint. My reading of the paper confirms that Eq. 5 quotes Δr_np = 0.12 ± 0.03 fm, Table II realizes Δr_np = 0.15 fm, and the abstract reports the resulting factor-of-2 reduction without that caveat. This is not an internal inconsistency, and the direction of the effect is physically plausible: a more diffuse neutron distribution in 96Zr increases the relative weight of neutron-driven peripheral interactions and correspondingly changes the proton contribution to the magnetic field. The concern is quantitative precision, not qualitative soundness. I found no stronger objection: the transport setup is standard, the SMASH code is public, the comparison to the no-skin baseline is within the same model, and the eccentricity finding agrees with prior v2 results. The missing uncertainty propagation and the extremal choice of Δr_np do not invalidate the study, but they do mean the abstract's factor-of-2 language overstates the robustness of the central estimate. This fully supports the reader's CONDITIONAL verdict; no change in verdict is needed.","tokens_in":9851,"tokens_out":4924,"duration_ms":58114,"concrete_test":"Repeat the SMASH runs for 96Zr+96Zr with Δr_np = 0.12 fm and Δr_np = 0.09 fm, keeping the proton Woods-Saxon parameters from Table II and recomputing the neutron diffusiveness via the Appendix A procedure, then recompute the ratio ⟨B²⟩Ru/⟨B²⟩Zr for 8 < b < 12 fm. If the reduction at the central value is substantially less than a factor of 2 (for example, if the ratio moves from 1.10 only to 1.07), the abstract should be revised to present the 5% value as an upper-limit scenario; if the factor-of-2 reduction persists at the central value, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that including the neutron skin of 96Zr reduces the Ru+Ru to Zr+Zr magnetic-field strength ratio from 10% to 5% in peripheral collisions. This reduction is obtained by fixing Δr_np = 0.15 fm, the upper edge of the measured 0.12 ± 0.03 fm from Eq. 5. The paper itself states that this value is chosen 'with the purpose of studying the neutron-skin impact at its extreme,' but the abstract and summary present the factor-of-2 reduction as the main result without this caveat. Since the no-skin baseline is a 10% difference and the skin effect is what moves the ratio toward unity, the size of the effect is plausibly roughly proportional to Δr_np over this range; at the central value 0.12 fm or the lower value 0.09 fm the reduction would likely be smaller than a factor of 2. The simulation also does not propagate the experimental uncertainty on Δr_np into the reported B-field ratio, so the headline number has no error bar and is best read as an extremal sensitivity estimate rather than a central prediction. The physical mechanism itself is plausible and the in-text caveat is honest, but the quantitative claim, as stated in the abstract, is not yet supported at the level claimed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the SMASH hadronic transport code to quantify how nuclear-structure details of the isobar pair 96Zr and 96Ru affect two CME-relevant observables in 200 GeV collisions: the participant eccentricity (a flow background proxy) and the early-time magnetic field strength (a CME signal proxy). The authors introduce isospin-dependent Woods-Saxon profiles for 96Zr with a neutron skin at the upper experimental limit (Δr_np = 0.15 fm), a deformation for 96Ru (β2 = 0.158), and nucleon-nucleon short-range correlations. Their main claim is that the neutron skin reduces the Ru+Ru to Zr+Zr magnetic field ratio difference from 10% to 5% in peripheral collisions (12 > b > 8 fm), which would imply a smaller CME signal-to-background ratio than previously expected. They also report an up-to-10% deformation-driven enhancement of the eccentricity ratio in ultra-central collisions.","tokens_in":10083,"tokens_out":6717,"duration_ms":67719,"significance":"The isobar program at RHIC was designed to separate CME signal from background by comparing Ru+Ru and Zr+Zr collisions; this paper is the first transport-model study to include the measured neutron skin of 96Zr in that context. If the central claim holds, it materially changes the expected sensitivity of the isobar run and informs the centrality selection for CME analyses. The eccentricity-deformation result is also a useful cross-check of earlier hydrodynamic work. The study is based on the publicly available SMASH code, and the nuclear-structure inputs are traceable to experimental data, which is a strength. However, the headline magnetic-field result is computed at the upper limit of the measured neutron skin and no uncertainties are propagated, so the quantitative claim is an extremal estimate rather than a central prediction.","major_comments":[{"comment":"The headline claim of a factor-of-2 reduction, from 10% to 5%, is computed only for Δr_np = 0.15 fm, the upper edge of the measured 0.12 ± 0.03 fm quoted in Eq. (5). The text honestly states that this is chosen 'with the purpose of studying the neutron-skin impact at its extreme,' but the abstract and summary present the factor-of-2 reduction without this caveat. Since the effect is expected to scale with Δr_np, the reduction at the central value 0.12 fm, or at the lower value 0.09 fm, would likely be smaller than a factor of 2. The authors should either compute the ratio for at least the central and lower values of Δr_np, or propagate the experimental uncertainty and rephrase the abstract to state that the factor-of-2 is an upper-limit sensitivity estimate.","section":"Abstract; Framework (Eq. 5)"},{"comment":"The ratio of magnetic field strengths, which is the central observable, is shown without statistical uncertainties or error bands. Given that the claimed effect is a change from a ratio of about 1.10 to about 1.05, it is essential to know whether this difference is significant relative to event-by-event fluctuations. Please add error bars or confidence bands to the ratio plot and, if possible, provide the statistical significance of the reduction in the peripheral bin 12 > b > 8 fm.","section":"Fig. 3 (bottom panel)"},{"comment":"The implementation of the neutron skin is of the 'neutron-halo' type (R0,n = R0,p, d_n > d_p), but the paper does not explicitly verify that the parameters listed in Table II indeed yield the intended Δr_np = 0.15 fm when used in Eqs. (9)-(12). Adding the resulting root-mean-square radii and the corresponding Δr_np to Table II would provide a useful consistency check that the sampled distributions realize the claimed extreme value.","section":"Table II; Appendix A"}],"minor_comments":[{"comment":"There are several typos and grammatical errors, including 'assumptios' (Framework), 'pannel' (Fig. 2 caption), 'diffusiveness' (should be 'diffuseness'), and 'refer to' (should be 'referred to'). The paper would benefit from a careful proofread.","section":"Throughout"},{"comment":"Reference [14] is listed as '(2019), arXiv:1906.03373 [nucl-ex]' without an author list; this appears to be the STAR isobar paper and the citation should be completed.","section":"Introduction"},{"comment":"The sentence explaining the mechanism states that the neutron skin 'enhances the number of neutron-neutron interactions in peripheral collisions or, equivalently, the concentration of protons in the central point that contribute to Eq. 8, leading to a larger B field.' This is confusing because neutron-neutron interactions do not directly contribute to the magnetic field. Please clarify the causal chain connecting the neutron skin to the enhanced B field at the central point.","section":"Signal: Magnetic field strength (Eq. 8)"},{"comment":"The cutoff R_i < 0.3 fm used to avoid singularities in the Lienard-Wiechert sum is a free parameter. Please provide a sensitivity check with respect to this cutoff, or at least justify the chosen value quantitatively.","section":"Signal: Magnetic field strength (Eq. 8)"},{"comment":"In the top panels of Figs. 2 and 3, the solid and dashed curves for the two systems overlap considerably, making them hard to distinguish. Adding markers or different line colors/widths would improve readability.","section":"Figs. 2 and 3"},{"comment":"Equations (11) and (12) use R0 and d on the right-hand side for the charge distribution, while R0,p and d_p denote point-proton values; this notation should be stated explicitly to avoid confusion.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and important question for the RHIC isobar program, and the physical mechanism is plausible. The main concern is that the quantitative headline result is presented for the upper-limit value of the neutron skin without uncertainty propagation, which overstates the certainty of the factor-of-2 reduction. This is fixable by additional simulations at the central and lower values or a proper error analysis, so I recommend major revision rather than rejection. The eccentricity part is solid but secondary."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Genuinely new: they put the measured neutron skin of 96Zr into the geometry used by SMASH and find it shrinks the Ru/Zr magnetic-field difference in peripheral collisions from 10% to 5%. The mechanism is clear and physically sensible: with neutrons pushed to the periphery, peripheral Zr+Zr collisions become more neutron-dominated, so Ru's proton excess matters less. The calculation is straightforward, SMASH is public, and the inputs come from independent nuclear data, not from fits to CME observables. The eccentricity result is also useful: Ru deformation raises the Ru/Zr eccentricity ratio by up to 10% in ultra-central collisions, while neutron skin and short-range correlations barely matter. That separation of effects is clean.\n\nSoft spot: the factor-of-2 headline is computed at the upper limit of the measured skin, Δr_np = 0.15 fm, while the central value is 0.12 ± 0.03 fm. The text says this, but the abstract and summary present the factor of 2 without the caveat. The effect is plausibly roughly linear in the skin, so at 0.12 or 0.09 fm the reduction would be smaller. No error bars in Fig. 3, no propagation of the experimental uncertainty into the ratio. That is a real, proportionate concern, not fatal.\n\nMinor points: the deformation choice (Ru deformed, Zr not) follows one set of measurements; the paper acknowledges an alternative scenario from the literature and picks one to isolate effects, which is fine but still an assumption. The Lienard-Wiechert cutoff at 0.3 fm is standard regularization, not a red flag.\n\nOverall, the direction is plausible and the simulation supports it. The abstract overstates the certainty of the quantitative claim. A rerun at central and lower skin values, or at least an uncertainty band, would fix that. The paper deserves serious refereeing; the referee should press on that point.","headline":"A focused, physically sensible transport study showing that the measured neutron skin of 96Zr can cut the Ru/Zr magnetic-field difference in half—but the headline number rests on the upper-limit skin value and needs an uncertainty band before it should be quoted as the central result.","tokens_in":665,"tokens_out":1396,"would_cite":true,"duration_ms":25543,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q"],"model":"deepseek-v4-flash","headline":"Neutron skin halves magnetic-field gap between isobar collision systems.","keywords":["neutron skin","isobar collisions","chiral magnetic effect","magnetic field strength","nuclear deformation","eccentricity","relativistic heavy-ion collisions","Liénard-Wiechert potentials"],"falsifier":"Rerun the same calculation with the central measured value $\\Delta r_{np} = 0.12$ fm or with a newer precision measurement of the skin; if the predicted Ru/Zr magnetic-field ratio in peripheral collisions stays near the original 10% excess rather than dropping to about 5%, the claimed cancellation fails. The experimental charge-separation data from the isobar run in peripheral centralities can serve as a direct check.","tokens_in":9608,"feed_emoji":"🧲","tokens_out":7389,"duration_ms":76685,"temperature":0.7,"pith_summary":"This paper studies collisions of two mirror nuclei with the same mass, ruthenium-96 and zirconium-96, at 200 GeV per nucleon pair, and asks whether an experimentally measured nuclear-structure feature changes the magnetic field that drives the chiral magnetic effect. The central finding is that including the neutron skin of zirconium-96—the fact that its neutrons extend farther from the center than its protons—cuts the expected magnetic-field difference between the two systems by half, from about 10% to about 5% in peripheral collisions. Because the CME signal is proportional to the field, the paper concludes that the isobar run will show a smaller charge-separation signal relative to background than earlier estimates suggested. It also finds that the deformation of ruthenium-96 raises the eccentricity ratio between the systems by up to 10% in ultra-central collisions, an effect that must be matched when comparing backgrounds.","feed_headline":"Neutron skin halves magnetic-field gap in isobar collisions","feed_subtitle":"Zirconium's neutron skin offsets ruthenium's protons, shrinking the predicted charge-separation signal in peripheral events.","key_machinery":"The mechanism is an isospin-dependent nuclear density profile: protons and neutrons are no longer sampled from one common Woods-Saxon distribution but are given separate radii and diffusivities, tuned so that zirconium-96 satisfies the measured skin $\\Delta r_{np} = 0.15$ fm while keeping the overall nucleus size fixed. This shifts neutrons to large radii and, in peripheral collisions, changes which nucleons interact, thereby altering the charge distribution that generates the magnetic field via Liénard-Wiechert potentials. The same geometric setup also includes a quadrupole deformation parameter $\\beta_2 = 0.158$ for ruthenium-96 and nucleon-nucleon short-range correlations, allowing the paper to isolate the skin's effect on both the magnetic field and the participant eccentricity.","core_discovery":"The central discovery is a cancellation. Ruthenium has more protons than zirconium and would normally generate stronger magnetic fields in peripheral collisions, but zirconium has a measured neutron skin that places extra neutrons in the nuclear surface. In peripheral collisions those neutrons dominate the interaction region, and the proton charge that contributes to the magnetic field concentrates near the center of the overlap, boosting the zirconium system's field. The paper finds that with the upper limit of the measured skin, $\\Delta r_{np} = 0.15$ fm, the ratio of squared magnetic-field strengths between Ru+Ru and Zr+Zr drops from a 10% excess to about 5% in peripheral collisions ($12 > b > 8$ fm), approaching unity. The same calculation shows that the eccentricity ratio between the two systems differs by up to 10% in ultra-central collisions, driven by ruthenium's deformation rather than by the neutron skin or nucleon correlations.","pith_inferences":["The authors deliberately use the upper limit of the measured skin; if a more precise measurement settles on the central value $\\Delta r_{np} = 0.12$ fm or lower, the same calculation would likely give a ratio closer to the original 10% difference, so the headline factor-of-two reduction is an extremal scenario rather than a central estimate.","The same cancellation logic should apply to other isobar pairs with a neutron-rich member that carries a skin; a systematic comparison across nickel or tin isobar pairs could test whether the suppression scales with skin size.","A ratio observable comparing charge separation in peripheral bins, where the field ratio is near one, against mid-central bins, where it is not, could isolate the CME contribution without relying on absolute magnetic-field predictions.","Because the neutron skin changes which nucleons participate in the collision, it should also alter photon and dilepton yields that track the proton fraction, offering independent experimental checks of the geometrical change."],"forward_implications":["Peripheral isobar events are expected to show a smaller CME signal relative to background than previously thought, because the magnetic-field difference driving the signal drops from 10% to about 5%.","The centrality window for a clean CME search shifts toward more peripheral events with $b > 12$ fm, where a sizeable field difference between the two systems reappears.","The eccentricity ratio between Ru+Ru and Zr+Zr is up to 10% larger in ultra-central collisions because of ruthenium's deformation, so background comparisons in that region must correct for flow differences; the paper suggests using $b > 6$ fm for matched backgrounds.","Nuclear-structure details such as the neutron skin and deformation need to be included in the interpretation of isobar data, while short-range correlations by themselves have a negligible effect on these observables."],"supporting_citations":[{"why":"Supplies the experimentally measured neutron skin of zirconium-96, $\\Delta r_{np} = 0.12 \\pm 0.03$ fm, which is the input for the isospin-dependent density profile.","marker":"[27, 28]"},{"why":"Provides the prior prediction of a 10% larger magnetic field for Ru+Ru than Zr+Zr in peripheral collisions, the baseline that the neutron skin is found to halve.","marker":"[23, 24]"},{"why":"Provides the relativistic hadronic transport simulation used to generate the nuclear collisions and compute the observables.","marker":"[29]"},{"why":"Defines the standard Woods-Saxon density profile that the paper modifies by introducing separate proton and neutron parameters.","marker":"[30]"},{"why":"Details the unfolding procedure that converts charge distributions into point-nucleon distributions, used to implement the neutron skin while keeping the nucleus size fixed.","marker":"[37, 40, 41]"},{"why":"Provides the previous hydrodynamic prediction of up to 10% flow differences due to deformation, which the eccentricity results are compared against.","marker":"[25]"}],"fun_headline_variants":["Neutron skin halves isobar magnetic gap","Neutron skin shrinks isobar field gap to 5%","Isobar magnetism gap cut in half by neutron skin","Neutron skin dulls ruthenium's magnetic edge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline factor-of-two reduction is computed by taking the experimentally measured neutron skin of zirconium-96 at its upper limit, 0.15 fm; if the true skin is at the central value 0.12 fm or below, the reduction weakens.","fun_headline_variants_meta":{"raw":{"variants":["Neutron skin halves isobar magnetic gap","Neutron skin shrinks isobar field gap to 5%","Isobar magnetism gap cut in half by neutron skin","Neutron skin dulls ruthenium's magnetic edge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000904,"raw_usage":{"total_tokens":3943,"prompt_tokens":1053,"completion_tokens":2890,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":2821}},"tokens_in":669,"tokens_out":2890,"duration_ms":21632,"temperature":1.0,"reasoning_tokens":2821,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:49:47.299034+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the same calculation with the central measured value $\\Delta r_{np} = 0.12$ fm or with a newer precision measurement of the skin; if the predicted Ru/Zr magnetic-field ratio in peripheral collisions stays near the original 10% excess rather than dropping to about 5%, the claimed cancellation fails. The experimental charge-separation data from the isobar run in peripheral centralities can serve as a direct check.","supporting_citations":[],"review_version":1}