{"id":"ad108be1-f988-40f1-b04a-9b59d642a353","arxiv_id":"2509.21644","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A centrality-scaled double ratio of net electric charge to net baryon number in p+208Pb collisions is predicted to be a new, measurable probe of the lead neutron skin.","lead":"The authors propose a new way to measure the neutron skin of lead by comparing electric charge and baryon number in peripheral versus central proton-lead collisions at the LHC. Their hydrodynamic simulations predict this ratio is sensitive to the neutron-rich surface layer of lead.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neglecting charge diffusion in Eq. (2) leaves the predicted neutron-skin sensitivity of the double ratio untested, and the paper's own footnote suggests diffusion is not a small correction at relevant rapidities.","rationale":"I read the paper as a simulation-based proposal: the claim is not that the model is the true QCD evolution, but that a clean, robust observable exists. The most exposed link in that argument is the longitudinal transport of conserved charges, because the observable is explicitly designed to compare charge ratios before and after a full hydrodynamic plus hadronic evolution. The paper's own footnote concedes that diffusion can qualitatively change the net-charge distribution at midrapidity; it then asserts, without a test, that forward-rapidity double ratios are immune. That assertion is directly load-bearing: if diffusion redistributes Q and B between central and peripheral events differently, the observed R could be dominated by transport rather than by the neutron skin. I agree with the reader's identification of this as the weakest assumption. The proposed computational test, adding diffusion currents and scanning their coefficients, would settle the question in a way that the current manuscript does not. Other concerns, such as the p+Pb centrality-impact-parameter correlation or the baryon-stopping model, are also present and partly acknowledged, but they are encoded in the initial-state model and are not as directly tied to a known neglected term in the evolution equations. The verdict should remain CONDITIONAL: the idea is worth pursuing, but the robustness claim needs the diffusion test before it can be considered established.","tokens_in":9449,"tokens_out":8200,"duration_ms":82712,"concrete_test":"Recompute the p+208Pb 5.02 TeV and 72 GeV predictions with Eq. (2) replaced by the full charge-current equations including baryon and electric-charge diffusion terms (e.g., the DNMR-type V^μ_{B,Q}=κ_{B,Q}(∇^μ(μ/T)-... ) used in Ref. [30]), scanning κ_B and κ_Q over a factor of 2 around the values adopted in that reference. Track the double ratio R^{Q,B}_{c1,c2} in 4.5<y_lab<5.5 for ΔRnp=0 and 0.28 fm. If the zero-skin baseline shifts by more than the skin-induced signal, or the slope dR/dΔRnp changes by more than about 20%, the robustness claim in the abstract is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For the central claim to hold, the Q/B double ratio must track the neutron skin as in Figs. 1–2, i.e., the longitudinal transport of electric charge and baryon number during the hydrodynamic phase must be centrality-independent apart from known effects. Equation (2) enforces this by construction: J^μ_{B,Q}=n_{B,Q}u^μ contains no dissipative diffusion term. The text justifies this by asserting that diffusion effects are 'similar' in the two centrality classes, but no quantitative estimate is provided. The diffusion current is driven by gradients of μ/T and acts over backgrounds with different temperatures, gradients, and lifetimes in 0–20% and 80–100% p+Pb events, so cancellation in the double ratio is not self-evident. Footnote 2 is direct evidence that the neglect is not harmless: including baryon diffusion can drive the net electric charge through zero at midrapidity, making the ratio singular there. If the zero crossing moves with centrality, the forward-rapidity window used here (4.5<y<5.5) could experience a sizable centrality-dependent shift in Q/B rather than the clean skin signal. Because the observable is a ratio of independently diffusing charges, this is a load-bearing, untested assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a new observable for constraining the neutron skin thickness of 208Pb in p+208Pb collisions: the double ratio of net electric charge to net baryon number (or net proton number) in forward rapidity windows, taken for peripheral (80-100%) events and divided by the same ratio for central (0-20%) events. The authors compute this observable with the iEBE-MUSIC (3+1)D viscous hydrodynamic framework, using 3D-Glauber initial conditions with baryon junctions, charge-dependent neos-4D equation of state, and a hadronic afterburner. Results are presented for LHC collider mode at 5.02 TeV and for the LHCb SMOG2 fixed-target configuration at 72 GeV. The paper finds that the double ratio decreases monotonically with increasing neutron skin thickness, with particularly clean sensitivity in collider mode, and recommends experimental measurement by LHCb.","tokens_in":9843,"tokens_out":3861,"duration_ms":39328,"significance":"If the central prediction is correct, this is a genuinely new and experimentally accessible avenue for neutron-skin measurements that is complementary to parity-violating electron scattering and to Bayesian analyses of flow observables. The paper has notable strengths: the observable is a forward output of a multi-stage model with no parameter fitted to the double-ratio target; the work uses a modern (3+1)D hydrodynamic description with multiple conserved charges and a lattice-QCD-based equation of state; initial-state and final-state results are compared; and the data are stated to be openly available. The proposed observable could help resolve the tension between the PREX-II result and ab initio calculations of the 208Pb neutron skin. However, several load-bearing transport assumptions are not quantitatively supported, and the predictions are presented without uncertainty quantification, so the sensitivity claim is not yet fully established.","major_comments":[{"comment":"The neglect of dissipative diffusion in the conserved charge currents is load-bearing for the central claim, because the double ratio directly compares the longitudinal transport of electric charge and baryon number between two centrality classes. The text asserts that diffusion effects are 'similar' in the two centralities, but no quantitative estimate is provided. More seriously, footnote 2 states that including baryon diffusion can make the net electric charge vanish at midrapidity, making the ratio impractical there. If the zero crossing moves with centrality, the forward window 4.5<y<5.5 used for the main prediction could receive a centrality-dependent shift in Q/B. I ask for a quantitative test: include a baryon and charge diffusion term in Eq. (2) (even with a simple diffusion coefficient), and show that the double-ratio sensitivity to the neutron skin in the chosen rapidity window is preserved, or quantify how much the inferred skin changes when diffusion is included.","section":"Eq. (2) and footnote 2"},{"comment":"The assumption that baryon stopping is energy independent is a strong simplification for predictions at two very different collision energies, 72 GeV and 5.02 TeV. The text mentions that a more rigorous treatment would account for energy dependence, and later claims the observable is 'robust to variations in stopping', but no stopping parameter scan is shown. Please provide a systematic variation of the baryon-junction stopping parameters and show the resulting spread in the double ratio. Without this, the quantitative curves in Figs. 1 and 2 may be tied to one specific stopping model, which is particularly relevant because the fixed-target forward window already deviates from unity by about 20% at zero neutron skin.","section":"Model description, baryon stopping assumption"},{"comment":"The predicted double-ratio curves are presented as single lines with no statistical error bars and no systematic envelope from the model inputs (Woods-Saxon radii, diffuseness parameters, shear and bulk viscosities, particlization energy density, and equation of state). Since the proposed experimental test must distinguish different values of ΔRnp, the reader needs to know the expected statistical precision of the quoted rapidity windows and the model-induced systematic uncertainty. Without this information, the statement that the probe is 'sensitive and robust' cannot be evaluated quantitatively. I recommend adding error bars from event statistics and a band from plausible variations of at least the diffuseness and viscosity parameters around their default values.","section":"Figs. 1-3 and uncertainty quantification"},{"comment":"The entire method relies on the correlation between event activity (multiplicity) and impact parameter in p+A collisions, but footnote 1 only states that the correlation is 'sizable enough' without a quantitative measure. Because p+A collisions have weaker multiplicity-impact-parameter correlations than A+A collisions, it is important to show the impact-parameter distribution and average number of participants in the 0-20% and 80-100% classes for the specific centrality definitions used. I ask for this diagnostic, and also for a demonstration that the centrality selection itself does not introduce a bias that mimics or obscures the neutron-skin signal.","section":"Footnote 1 and centrality selection"}],"minor_comments":[{"comment":"The caption contains a duplicated word: 'triangle and and star markers' should read 'triangle and star markers'.","section":"Fig. 1 caption"},{"comment":"The centrality classes are defined by multiplicity in rapidity intervals 1.25<y_lab<3 (collider) and 0.6<y_lab<2.6 (fixed target), but the text does not specify whether these are pseudo-rapidity or rapidity intervals. Please clarify the exact observable used for the event-activity selection.","section":"Centrality definition, Section 3"},{"comment":"In the fixed-target forward window the double ratio deviates from unity by about 20% even at zero neutron skin; the manuscript mentions this but does not quantify the impact of this baseline on the extraction precision. A few sentences with the expected signal-to-baseline ratio would help the reader judge the practical sensitivity.","section":"Fig. 2 and Fig. 3 discussion"},{"comment":"The statement 'The data that support the findings of this article are openly available [41], embargo periods may apply' is vague. Please specify the repository or DOI and the embargo status, since reproducibility is a strength of this work.","section":"Data availability statement"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper proposes something genuinely new: a centrality-scaled double ratio of net electric charge to net baryon (or proton) number in p+Pb collisions, measured at forward rapidity, as a probe of the 208Pb neutron skin. I have not seen this specific observable before, and the logic is sound—peripheral collisions preferentially sample the neutron-rich surface, so the Q/B ratio shifts with skin thickness. The simulation study is thorough: (3+1)D hydrodynamics, charge-dependent lattice EoS, full hadronic afterburner, and results at both LHC collider and SMOG2 fixed-target energies. The paper is honest about its own limitations, which matters here.\n\nWhat the paper does well: the observable is experimentally accessible in principle (net charge and net protons in a forward rapidity window), and the predicted double ratio shows a clear, monotonic, roughly linear dependence on the skin thickness. The Q/p proxy works almost as well as Q/B, which is useful. The authors also show that the sensitivity is robust in the fragmentation region, where initial-state and particlization effects are smaller.\n\nThe soft spot is the one you flagged. The hydrodynamic currents in Eq. (2) contain no diffusion term, and the authors justify this by asserting diffusion effects are \"similar\" across centralities. That is not demonstrated. Worse, footnote 2 says that including baryon diffusion can drive the net charge through zero at midrapidity—so diffusion is not a small effect everywhere. If the zero crossing moves with centrality, the forward-rapidity window could get a centrality-dependent shift that masquerades as a skin signal. This is a load-bearing assumption for the central claim, and the authors themselves provide a hint that it may fail.\n\nA second, smaller issue: the baryon stopping is assumed energy-independent, which the authors acknowledge. And the predictions come with no uncertainty quantification—no statistical or systematic error bands on the curves. That is normal for a proposal paper, but it means the claimed 4% vs 20% deviations are not yet tied to experimental precision.\n\nBottom line: this is a promising proposal, not a validated measurement. The diffusion issue is the one thing I would want addressed before betting on it. But the idea is new, the simulation framework is state-of-the-art, and the authors are upfront about the main caveats. It deserves a serious referee and, ideally, a follow-up with diffusion included.\n\nRecommendation: engage with it, but treat the quantitative sensitivity as provisional until the charge diffusion question is settled.\n\nBest,","headline":"A genuinely new observable for the lead neutron skin, with a clean simulation-based sensitivity study whose main vulnerability is the neglected charge diffusion, explicitly acknowledged in footnote 2.","tokens_in":794,"tokens_out":958,"would_cite":true,"duration_ms":14115,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","21.10.Gv","24.10.Nz"],"model":"deepseek-v4-flash","headline":"A double ratio of charge yields in p+Pb collisions can measure the neutron skin of 208Pb.","keywords":["neutron skin","lead-208","proton-lead collisions","net electric charge","net baryon number","relativistic hydrodynamics","LHCb","symmetry energy"],"falsifier":"Compute the same double ratio with nonzero baryon and electric-charge diffusion coefficients in the hydrodynamic currents; if the centrality-differential shift in the 4.5 < ylab < 5.5 window at √sNN = 5.02 TeV is as large as the change produced by varying ΔRnp by 0.11 fm, the claimed sensitivity is not robust. Experimentally, an LHCb measurement of the net-charge-to-proton double ratio in that window that shows a flat or opposite centrality dependence would falsify the prediction.","tokens_in":9248,"feed_emoji":"⚛️","tokens_out":6561,"duration_ms":55754,"temperature":0.7,"pith_summary":"The paper proposes that the neutron skin of a heavy nucleus—the excess of neutrons at the nuclear surface—can be measured at colliders through conserved charge yields in proton–lead collisions, without parity-violating electron scattering. The central observable is a double ratio: the net-electric-charge-to-net-baryon ratio in peripheral (80–100%) events divided by the same ratio in central (0–20%) events, evaluated in a forward rapidity window on the lead-going side. Using event-by-event (3+1)-dimensional relativistic hydrodynamics, the authors show this double ratio decreases approximately linearly with the neutron skin thickness of 208Pb, and that the net-charge-to-proton version in the 4.5 < y < 5.5 window at LHC collider mode is ready for near-term experimental test. Such a measurement would give a new, independent constraint on the density dependence of the nuclear symmetry energy and could help adjudicate the tension between the PREX-II result and ab initio calculations.","feed_headline":"Lead's neutron skin shows up in charge ratios at the LHC","feed_subtitle":"A forward-rapidity double ratio in p+Pb collisions could settle the size of lead's neutron skin.","key_machinery":"The load-bearing object is the double ratio $R_{c_1,c_2}^{X,Y} = [N_X(c_1)/N_Y(c_1)] \\div [N_X(c_2)/N_Y(c_2)]$, with $c_1 = 80\\text{--}100\\%$ and $c_2 = 0\\text{--}20\\%$. It isolates centrality-dependent changes in the ratio of net electric charge ($X=Q$) to net baryon number ($Y=B$) or to net protons ($Y=p$), so that the geometry-driven effect of the neutron-rich surface survives while overall normalization cancels. The supporting machinery is the initial-state sampling with isospin-dependent Woods–Saxon profiles whose neutron diffuseness is varied to dial $\\Delta R_{np}$, followed by (3+1)-dimensional viscous hydrodynamic evolution with conserved charge currents and a charge-dependent lattice-QCD equation of state.","core_discovery":"The paper's central claim is that the centrality-dependent double ratio of net electric charge to net baryon number in p+208Pb collisions is a sensitive and robust probe of the lead neutron skin. Low-multiplicity (peripheral) events preferentially sample the nuclear surface, so they involve more proton–neutron interactions than high-multiplicity (central) events; the net electric charge per baryon is therefore suppressed in peripheral relative to central events, and this suppression grows with the neutron skin thickness. The authors compute this double ratio with a (3+1)-dimensional viscous hydrodynamic model initialized with isospin-dependent Woods–Saxon densities and baryon-junction stopping, and report an approximately linear fall of the double ratio with ΔRnp. They provide predictions at √sNN = 5.02 TeV and 72 GeV, and show that replacing net baryons by net protons, the experimentally accessible proxy, gives nearly identical results.","pith_inferences":["The calibration curve could shift if dissipative charge diffusion does not affect peripheral and central events equally; rerunning the simulation with finite baryon and electric-charge diffusion coefficients would test this directly.","Combining this collider-based extraction with PREX-II and the Bayesian heavy-ion extraction cited in the paper could yield a tighter joint posterior on the symmetry-energy slope L than any single method alone.","Because the observable is built from event-class ratios rather than absolute yields, part of the detector acceptance error may cancel; a realistic LHCb acceptance study would clarify whether the 4.5 < y < 5.5 window is indeed feasible in the near term.","The linear dependence on ΔRnp suggests that even an upper bound from a first measurement would already narrow the range of allowed symmetry-energy slopes, an implication the paper does not explicitly draw."],"forward_implications":["LHCb can test the collider-mode prediction by measuring the net-charge-to-proton double ratio in the 4.5 < ylab < 5.5 window for 0–20% and 80–100% centrality classes.","A confirmed signal would provide an independent cross-check of the PREX-II extraction of the 208Pb neutron skin.","Because the double ratio reacts almost linearly to ΔRnp, a single measurement would calibrate the skin thickness once the relation is fixed by the simulation.","The same methodology extends to other neutron-rich nuclei, such as 48Ca, and to the fixed-target 208Pb-on-hydrogen setup at SMOG2.","Measurements in the deeper fragmentation region 5.5 < ylab < 7.5 are predicted to be even cleaner but require detector coverage beyond current LHCb acceptance."],"supporting_citations":[{"why":"Supplies the PREX-II measurement of the 208Pb neutron skin that serves as the experimental comparison band.","marker":"[10]"},{"why":"Provides the ab initio neutron-skin range used as the other comparison band and establishes the tension with PREX-II.","marker":"[12]"},{"why":"Demonstrates the prior Bayesian extraction of the lead neutron skin from relativistic nuclear collisions, the approach this paper contrasts with.","marker":"[13]"},{"why":"Defines the dynamical initial-state model from which the collision geometry and charge distributions are built.","marker":"[21]"},{"why":"Adds the longitudinal dynamics, including baryon-junction stopping, that transports baryon charge toward midrapidity.","marker":"[22]"},{"why":"Sets up the conserved-charge transport treatment in isobar collisions on which the present charge currents rely.","marker":"[26]"},{"why":"Derives the net baryon diffusion terms whose omission is the paper's stated approximation.","marker":"[30]"},{"why":"Provides the four-dimensional QCD equation of state with baryon, electric-charge, and strangeness chemical potentials.","marker":"[34]"}],"fun_headline_variants":["Charge ratio in p+Pb reveals lead's neutron skin","Neutron skin measured via charge-to-baryon ratio in colliders","New probe: charge ratios size up lead's neutron skin","Lead's neutron skin extracted from charge asymmetries at LHC","Forward-rapidity charge ratio pins down lead neutron skin"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction assumes that dissipative diffusion of conserved charges, which is omitted from the hydrodynamic currents, affects peripheral and central events in the same way; if baryon or electric-charge diffusion differs between event classes at forward rapidity, the calibration of the double ratio to the neutron skin thickness could shift substantially.","fun_headline_variants_meta":{"raw":{"variants":["Charge ratio in p+Pb reveals lead's neutron skin","Neutron skin measured via charge-to-baryon ratio in colliders","New probe: charge ratios size up lead's neutron skin","Lead's neutron skin extracted from charge asymmetries at LHC","Forward-rapidity charge ratio pins down lead neutron skin"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1591,"prompt_tokens":913,"completion_tokens":678,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":529,"completion_tokens_details":{"reasoning_tokens":592}},"tokens_in":529,"tokens_out":678,"duration_ms":5607,"temperature":1.0,"reasoning_tokens":592,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:45:34.893306+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same double ratio with nonzero baryon and electric-charge diffusion coefficients in the hydrodynamic currents; if the centrality-differential shift in the 4.5 < ylab < 5.5 window at √sNN = 5.02 TeV is as large as the change produced by varying ΔRnp by 0.11 fm, the claimed sensitivity is not robust. Experimentally, an LHCb measurement of the net-charge-to-proton double ratio in that window that shows a flat or opposite centrality dependence would falsify the prediction.","supporting_citations":[{"cited_title":"Baryon Transport in Color Flux Tubes","cited_arxiv_id":"2311.17906","evidence_quote":"Adds the longitudinal dynamics, including baryon-junction stopping, that transports baryon charge toward midrapidity."},{"cited_title":"Zenihiroet al., Neutron density distributions of Pb- 204, Pb-206, Pb-208 deduced via proton elastic scattering at Ep=295 MeV, Phys","cited_arxiv_id":null,"evidence_quote":"Sets up the conserved-charge transport treatment in isobar collisions on which the present charge currents rely."}],"review_version":1}