{"id":"d97c86c5-5a84-4295-ad3d-49dec060039f","arxiv_id":"2411.13622","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A double ratio of hyperon-antihyperon production in antiproton-nucleus collisions correlates strongly with neutron skin differences between isotopes, in both a schematic model and GiBUU transport simulations.","lead":"This paper proposes a new way to measure tiny differences in the neutron skin, the extra layer of neutrons at the surface of atomic nuclei, by comparing how often antiprotons create two different kinds of particle pairs when they hit different isotopes of the same element. The method could offer a precise, independent probe of neutron skins, relevant for understanding neutron stars, and could be tested at the PANDA experiment.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed 0.999 Pearson correlation between the schematic and GiBUU double ratios is driven by the shared input density distributions; a separate check of the method's sensitivity to genuine skin differences is needed.","rationale":"The paper is a well-crafted feasibility study with clear physics motivation and a plausible observable. The strongest claim is the nearly perfect correlation of Fig. 8, which is used to argue that the schematic model captures the essential dynamics. My concern is more specific than the reader's factorization worry: the correlation is computed between two quantities that both derive from the same RMF density distributions within the same code environment. What would be needed to sustain the claim of a direct measure is an external validation—either against data, or against a transport model with different absorption/annihilation physics, or against densities from a structurally different nuclear model. Without that, the 0.999 Pearson coefficient is suggestive but not yet a validation of Eq. 10. The practical recommendation (CONDITIONAL) is unchanged from the reader: accept the proposal as promising, but require validation on independent input before claiming a precision measurement. The proposed test—swapping the density generator—is concrete and would settle whether the correlation is robust.","tokens_in":26773,"tokens_out":795,"duration_ms":12042,"concrete_test":"Run the same GiBUU comparison using density distributions from a completely different nuclear structure model (e.g., Skyrme-Hartree-Fock or ab initio results for 48Ca) that yield different neutron skins, and recompute both the Sec. III schematic DR and the GiBUU DR. If the correlation coefficient between schematic and GiBUU remains near 0.999 when the input densities are varied independently of the RMF set used to generate them, the correlation is likely physical. If it degrades significantly, the correlation is an artifact of shared input. Also check the slope: the linear fit in Fig. 8 has slope ~0.9, not 1; a first-principles derivation of Eq. 10 should predict the slope, and the deviation should be explained before claiming a direct measure.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim—that Eq. 10 makes DR a direct measure of integrated neutron-skin variation—rests on the correlation shown in Fig. 8. But both arms of that verification use the same RMF-generated density distributions as input: the schematic model of Sec. III integrates the very same proton and neutron areal densities that the GiBUU initialization code samples for its phase-space distributions. A strong correlation therefore may reflect only that both calculations respond similarly to the input densities, not that the schematic absorption formula captures the relevant dynamics. Moreover, the GiBUU model itself contains the same physical assumptions about absorption and production as the schematic picture, so the comparison is not an independent test of Eq. 10. The reader's weakest_assumption—factorization κ_II = κ_I·κ_n and κ_ΛΛ ≈ κ_Σ−Λ—is real, but the more immediate gap is that the verification is partly circular with respect to the density input. The quantitative claim that DR ∝ p_abs with slope (1 + Z/N) is not tested against any data or against an alternative transport model with different absorption treatment; it is only compared with another calculation that shares the same density model and similar reaction assumptions. Thus the precision estimate of ±10% on ∆n is not robust until the method is tested against a genuinely independent calculation or measurement.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new method to probe variations of the neutron skin thickness between isotopes using the double ratio of Σ−Λ and ΛΛ pair production in antiproton–nucleus interactions. Building on a simple geometrical model, the authors derive an approximate linear relation, DR ≈ 1 + (1 + Z/N) p_abs (Eq. 10), and verify it via high-statistics GiBUU transport simulations for several isotope pairs (Ne, Ca, Ni, Xe), finding a Pearson correlation of 0.999 between the schematic and transport double ratios (Fig. 8). The paper also discusses sensitivity to density parameters, estimates an uncertainty of about ±10% on the neutron-skin variation for a ±1% double-ratio measurement, and outlines experimental prospects for PANDA.","tokens_in":27095,"tokens_out":3712,"duration_ms":42922,"significance":"If the proposed method is validated, it would provide a new, potentially high-precision observable for neutron-skin studies that is complementary to existing probes and is particularly sensitive to small variations along isotope chains. The analytical derivation is transparent and parameter-free, and the GiBUU calculations are extensive and reproducible in principle. The paper also offers a useful discussion of systematic uncertainties and experimental feasibility. The credibility of the quantitative uncertainty claim, however, depends on the degree to which the GiBUU comparison constitutes an independent verification, which is the main concern of this report.","major_comments":[{"comment":"The verification of the central relation (Eq. 10) is partly circular with respect to the input densities. The schematic model of §III integrates the same RMF-generated proton and neutron areal densities that are used to initialize the GiBUU phase-space distributions. Therefore the Pearson correlation of 0.999 may largely reflect that both calculations respond to the same density input in a similar way, rather than that the simplified absorption formula captures the dynamics of hyperon production and absorption. The paper does not compare against an independent transport code or against synthetic data generated with different density parametrizations. This weakness directly affects the precision claim in §V B (Figure 9), where the ±10% mapping from double-ratio uncertainty to neutron-skin variation is derived from a comparison that shares the same density input and similar reaction assumptions. I recommend that the authors either provide such an independent test or explicitly qualify the correlation as a consistency check rather than a validation.","section":"§IV, Fig. 8"},{"comment":"The factorization ansatz κ_II = κ_I · κ_n and the equality κ_ΛΛ ≈ κ_Σ−Λ are load-bearing assumptions in the derivation of Eq. (10). Neither is derived from first principles, and the GiBUU comparison does not isolate their validity because GiBUU includes full dynamics. The authors should provide a quantitative assessment of how strongly the double ratio depends on these assumptions, for example by running GiBUU with artificially modified absorption cross sections for Λ and Σ−, or by comparing the schematic result with a version of the schematic model where these equalities are relaxed.","section":"§II, after Eq. (7)"},{"comment":"The sensitivity study in Figure 9 varies only the neutron distribution of 48Ca while keeping the proton distribution fixed. In reality, as the paper itself notes (e.g., §III list of deficiencies and §IV C), proton distributions also differ between isotopes and may be correlated with neutron changes. The two-parameter scan in Figure 10 does include variations of Rn and an for the neutron distribution, but it does not address correlated variations between protons and neutrons. The claimed precision of ±10% on neutron-skin variation should be justified under a broader set of density variations that includes correlated proton and neutron changes.","section":"§V, Figures 9 and 10"}],"minor_comments":[{"comment":"The caption states 'for 40Ca (blue lines) and 48Ca (red lines)', but the text in §III describes protons as red and neutrons as blue. Please correct the inconsistency.","section":"Figure 3 caption"},{"comment":"The captions refer to lines marking half-density radii, while the text in §IV A refers to vertical dashed lines indicating rms radii. Please align the description of the graphical elements.","section":"Figures 5 and 6 captions"},{"comment":"Several instances of 'double ration' should read 'double ratio' (e.g., §IV C, §V A).","section":"Throughout"},{"comment":"The derivation from Eq. (7) to the double ratio is compact; a short intermediate step showing the cancellation of the common factors would improve readability.","section":"Eq. (8)"},{"comment":"The statement that 'the statistical precision shown in Table II e.g. for calcium can be reached in about half a day running for each isotope' would benefit from a brief justification based on the production yields and the assumed detection efficiencies, as the reader otherwise cannot judge the feasibility.","section":"§VI"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-motivated proposal with a transparent analytical derivation and a significant transport effort. The main weakness is that the verification against GiBUU does not constitute an independent test of the method because both sides of the comparison use the same input densities and similar physics assumptions. This needs to be addressed, or at least explicitly and rigorously qualified, before the quantitative precision claim can be considered established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Mike, this one is worth a look. The authors propose using the ratio of Σ−Λ to ΛΛ production in antiproton–nucleus collisions as a differential probe of neutron skin variations between isotopes. The trick is that at threshold, ΛΛ comes from p–p and Σ−Λ only from p–n, so the double ratio isolates the neutron content of the periphery. The analytic result, Eq. 10, DR ≈ 1 + (1+Z/N) p_abs, is clean and not fitted. The schematic model in Section III is transparent and explicitly lists its simplifications. The GiBUU transport study is high-statistics and spans Ne, Ca, Ni, and Xe isotopes, and the 0.999 Pearson correlation between the schematic and transport double ratios is a striking consistency check.\n\nThe soft spot is that the two arms of the verification share the same input densities—the RMF-generated proton and neutron distributions feed both the schematic areal density integrals and the GiBUU initialization. So Fig. 8 shows that the simple geometric treatment reproduces the transport dynamics, but it does not independently validate the density-to-observable mapping. That matters for the claim that DR is a 'direct measure' of skin increments: the calibration is only as good as the density models. The factorization assumptions in Section II (κ_II = κ_I·κ_n, κ_ΛΛ ≈ κ_Σ−Λ) are plausible but unvalidated, and the ±10% precision estimate rests on a ±1% DR measurement that has not yet been demonstrated experimentally. These are addressable, not fatal.\n\nThe paper is honest about its limitations and the literature citation is thorough. I would send it to a referee. The right referee would push for a cross-check with a different transport code or an independent density functional, and for a more detailed experimental systematic budget. If the correlation survives that, this is a genuinely useful new probe.","headline":"A clever new double-ratio observable for neutron skin studies, with a clean analytic formula and strong transport-based evidence, though the verification is partly tied to the shared input densities.","tokens_in":27550,"tokens_out":3058,"would_cite":true,"duration_ms":32549,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35","81V05"],"pacs":["25.43.+t","21.10.Gv"],"model":"deepseek-v4-flash","headline":"A double ratio of near-threshold $\\Sigma^-\\overline{\\Lambda}$ and $\\Lambda\\overline{\\Lambda}$ production in antiproton–nucleus collisions is proposed and shown to be a direct measure of the change in neutron skin thickness between two…","keywords":["neutron skin thickness","isotope chains","antiproton-nucleus collisions","hyperon-antihyperon pair production","double ratio observable","Boltzmann-Uehling-Uhlenbeck transport","relativistic mean-field densities","nuclear periphery"],"falsifier":"Measure $\\mathrm{DR}$ for $^{40}$Ca/$^{48}$Ca and for several xenon isotope pairs at a $2.4$ GeV/$c$ antiproton beam, using the measured $\\overline{p}n$ cross section to compute $p_{\\mathrm{abs}}$ from the integrated skin density; if the data fall off the common linear relation $\\mathrm{DR}=1+(1+Z/N)\\,p_{\\mathrm{abs}}$ across pairs, or if the transport model with artificially unequal $\\Lambda$ and $\\Sigma^-$ absorption inside the skin changes the double ratio substantially, the central claim is falsified.","tokens_in":26615,"feed_emoji":"⚛️","tokens_out":8999,"duration_ms":86704,"temperature":0.7,"pith_summary":"The paper proposes that the evolution of the neutron periphery between two isotopes can be read from a single dimensionless double ratio: the production of $\\Sigma^-\\overline{\\Lambda}$ pairs (which require antiproton–neutron collisions) divided by $\\Lambda\\overline{\\Lambda}$ pairs (which require antiproton–proton collisions), measured for both isotopes. Because antiprotons are strongly absorbed in the nuclear periphery, the extra neutrons of the heavier isotope act as an absorber for the $\\Lambda\\overline{\\Lambda}$ channel and as an extra source for the $\\Sigma^-\\overline{\\Lambda}$ channel, and the paper derives the approximate linear relation $\\mathrm{DR}\\approx 1+(1+Z/N)\\,p_{\\mathrm{abs}}$ linking the double ratio to the absorption probability in the added neutron layer. Full transport simulations for neon, calcium, nickel, and xenon isotope pairs confirm that this schematic double ratio tracks the transport double ratio with a Pearson correlation of $0.999$ over a wide mass range. If the relation holds, the method offers a precision tool for tiny neutron-skin variations along isotope chains that is complementary to parity-violating electron scattering and could help settle the current tension between the results for $^{208}$Pb and $^{48}$Ca.","feed_headline":"Double ratio reads neutron skin variation between isotopes","feed_subtitle":"Near threshold, the ratio of Σ−Λ to ΛΛ pairs across isotopes tracks the added neutron layer to 0.999 correlation.","key_machinery":"The load-bearing object is the double ratio $\\mathrm{DR}=(\\Sigma^-\\overline{\\Lambda}/\\Lambda\\overline{\\Lambda})_{\\mathrm{II}}/(\\Sigma^-\\overline{\\Lambda}/\\Lambda\\overline{\\Lambda})_{\\mathrm{I}}$ for a heavier isotope II relative to a reference isotope I. In the schematic picture, the extra neutron layer of II suppresses $\\Lambda\\overline{\\Lambda}$ production through absorption of the incoming antiproton with probability $p_{\\mathrm{abs}}=1-\\exp(-\\sigma_{pn}\\int_{\\Delta n}\\rho_n\\,d^3r)$, while adding a new neutron-only production region for $\\Sigma^-\\overline{\\Lambda}$; expanding in small $p_{\\mathrm{abs}}$ yields $\\mathrm{DR}\\approx 1+(1+Z/N)p_{\\mathrm{abs}}$. The finite-impact-parameter version replaces the analytic formula by the ratio of the neutron and proton areal densities within one interaction length, integrated over all impact parameters, and this schematic quantity is what correlates with the transport result.","core_discovery":"Close to the production threshold, $\\Lambda\\overline{\\Lambda}$ pairs are produced only in $\\overline{p}p$ sub-collisions and $\\Sigma^-\\overline{\\Lambda}$ pairs only in $\\overline{p}n$ sub-collisions, so the proton and neutron content of the nuclear periphery are imprinted differently on the two channels. The paper's central claim is that the double ratio $\\mathrm{DR}$ between two isotopes, defined by Eq. (1), is a direct measure of the increment of the integrated neutron skin thickness: in the small-absorption limit $\\mathrm{DR}\\approx 1+(1+Z/N)\\,p_{\\mathrm{abs}}$, where $p_{\\mathrm{abs}}$ is the probability that the antiproton is absorbed in the additional outer neutron layer of the heavier isotope. The claim is supported by a Boltzmann–Uehling–Uhlenbeck transport study of $\\overline{p}+A$ reactions at $2.4$ GeV/$c$: the double ratio computed from the neutron-to-proton content of the periphery within one interaction length, integrated over all impact parameters, is linearly correlated with the full transport double ratio with Pearson coefficient $0.999$, and the authors show that a $1\\%$ uncertainty on $\\mathrm{DR}$ for the $^{40}$Ca/$^{48}$Ca pair translates into roughly $10\\%$ uncertainty on the neutron-skin variation.","pith_inferences":["Because the isolines of $\\mathrm{DR}$ and of the rms neutron-skin thickness in the $(R_n,a_n)$ plane are not parallel, a joint analysis of $\\mathrm{DR}$ with an independent peripheral probe such as antiprotonic-atom x-rays could separate changes in the half-density radius from changes in the surface diffuseness; the paper notes the combination would be valuable but does not construct such a fit.","If precise $\\mathrm{DR}$ values were obtained across a long isotope chain, the chain itself would provide a differential map of how the neutron periphery grows with neutron number, indirectly constraining the density dependence of the symmetry energy that enters neutron-star equations of state; the paper only invokes the known general correlation.","The factorisation assumption $\\kappa_{\\mathrm{II}}=\\kappa_{\\mathrm{I}}\\kappa_n$ could be tested directly in the transport model by artificially setting the $\\Lambda$ and $\\Sigma^-$ absorption cross sections equal and then unequal inside the added skin layer; the change in the resulting $\\mathrm{DR}$ would quantify the leading systematic of the method, a test the paper does not perform."],"forward_implications":["A measurement of $\\mathrm{DR}$ for $^{40}$Ca and $^{48}$Ca at the percent level would constrain the neutron-skin variation between the two isotopes to about $10\\%$, roughly a factor of three better than the CREX uncertainty.","Absolute cross sections are not needed, and the different energy dependences of the $\\Lambda\\overline{\\Lambda}$ and $\\Sigma^-\\overline{\\Lambda}$ channels cancel in the double ratio, reducing many experimental systematics.","For the xenon chain $^{129}$Xe to $^{136}$Xe, where the neutron-skin variation is only about $0.07$ fm, the predicted double ratio rises from about $1.02$ to $1.17$, so small skin increments are experimentally visible.","The incident antiproton momentum has only a small effect on the double ratio, so the method is robust around the $2.4$ GeV/$c$ operating point.","An analogous construction applied to isotone chains could probe the evolution of proton skins rather than neutron skins."],"supporting_citations":[{"why":"Supplies the full Boltzmann–Uehling–Uhlenbeck transport model whose double ratios are the benchmark for the schematic correlation.","marker":"[56]"},{"why":"Provides the relativistic mean-field parameter sets that generate the proton and neutron densities used in both the schematic and transport calculations.","marker":"[58]"},{"why":"Gives the $\\overline{p}p$ and $\\overline{p}n$ cross sections (about 55 mb at 2.4 GeV/$c$) that fix the absorption probability and interaction length.","marker":"[63]"},{"why":"Establishes that $^{40}$Ca and $^{48}$Ca have essentially equal charge radii, making that pair the paper's showcase case.","marker":"[60]"},{"why":"Provides the $^{48}$Ca parity-violating electron scattering result whose uncertainty the proposed method would improve by about a factor of three.","marker":"[25]"},{"why":"Provides the $^{208}$Pb parity-violating electron scattering result that, together with the $^{48}$Ca result, frames the neutron-skin tension the method is intended to address.","marker":"[24]"}],"fun_headline_variants":["Spot neutron skin change with hyperon pair double ratio","Hyperon double ratio measures neutron skin variation","Antiproton + nucleus: new probe of neutron skin","Hyperon pairs map neutron skin difference"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method requires that inside the added outer neutron layer the produced $\\Lambda$ and $\\Sigma^-$ hyperons are absorbed with the same probability, and that this extra-layer absorption factorizes as a simple multiplier on the core absorption; if the two hyperon species are absorbed differently in the skin, $\\mathrm{DR}$ is no longer a clean function of the absorption probability.","fun_headline_variants_meta":{"raw":{"variants":["Spot neutron skin change with hyperon pair double ratio","Hyperon double ratio measures neutron skin variation","Antiproton + nucleus: new probe of neutron skin","Hyperon pairs map neutron skin difference"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000716,"raw_usage":{"total_tokens":3271,"prompt_tokens":1049,"completion_tokens":2222,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":2163}},"tokens_in":665,"tokens_out":2222,"duration_ms":17827,"temperature":1.0,"reasoning_tokens":2163,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:39:13.259592+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\mathrm{DR}$ for $^{40}$Ca/$^{48}$Ca and for several xenon isotope pairs at a $2.4$ GeV/$c$ antiproton beam, using the measured $\\overline{p}n$ cross section to compute $p_{\\mathrm{abs}}$ from the integrated skin density; if the data fall off the common linear relation $\\mathrm{DR}=1+(1+Z/N)\\,p_{\\mathrm{abs}}$ across pairs, or if the transport model with artificially unequal $\\Lambda$ and $\\Sigma^-$ absorption inside the skin changes the double ratio substantially, the central claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the full Boltzmann–Uehling–Uhlenbeck transport model whose double ratios are the benchmark for the schematic correlation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the relativistic mean-field parameter sets that generate the proton and neutron densities used in both the schematic and transport calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the $\\overline{p}p$ and $\\overline{p}n$ cross sections (about 55 mb at 2.4 GeV/$c$) that fix the absorption probability and interaction length."},{"cited_title":"Reinhard, X","cited_arxiv_id":null,"evidence_quote":"Establishes that $^{40}$Ca and $^{48}$Ca have essentially equal charge radii, making that pair the paper's showcase case."}],"review_version":1}