{"id":"9ac5ddf7-aaf7-4cd3-be30-ae642b6a55de","arxiv_id":"2602.21011","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"New measurements of the 144Sm(alpha,n)147Gd reaction cross section at five sub-Coulomb energies between 14 and 21 MeV, with a detailed covariance analysis.","lead":"A nuclear physics group measured how often alpha particles trigger a reaction in samarium-144 at energies below the Coulomb barrier, producing five new cross-section values with detailed uncertainty correlations. The work feeds models of how rare heavy elements are made in stars and supports production of a medical isotope.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported cross sections are assigned to the mean beam energy, but the activation yield averages the cross section over the energy distribution; at the steepest sub-Coulomb points this can bias the data upward by 10-25%.","rationale":"The paper is an internally consistent activation measurement, and the reader's verdict of CONDITIONAL is reasonable. However, the reader's weakest assumption focuses on the target-thickness scale, which is explicitly included in the quoted 15-25% uncertainties. A more serious and unaddressed systematic is the neglect of the energy-spread convolution. At sub-Coulomb energies, the cross section rises exponentially with energy (dlnσ/dE ≈ πη/E). The activation method measures the yield integrated over the full beam-energy distribution, but the paper assigns that yield to the mean energy of the distribution. Using a Gaussian approximation for the energy spread (as the paper's own fits in Figures 3-4 suggest), the resulting bias <σ>/σ(Ē) is approximately exp(½(πη/E)²σ_E²). For the lowest-energy point (14.09 MeV, σ_E=0.30 MeV) this amounts to a +25% systematic excess, several times larger than the stated cross-section uncertainty for that point and nearly equal to it. The effect diminishes at higher energies but still reaches +12% at 16.03 MeV. This bias is energy-dependent and would steepen the low-energy excitation function, which is precisely the regime most relevant for p-process nucleosynthesis and for testing Hauser-Feshbach models. The covariance analysis in Table 8 does not include this effect, and the text does not mention it. The proposed concrete test—convolving a TALYS cross-section with the GEANT4 energy distributions—would directly quantify the bias. If the bias is confirmed, the reported cross sections at 14-16 MeV are systematically too high and the paper's central data set would require correction or a redefinition of effective energies. For these reasons, the verdict remains CONDITIONAL, but the condition should include a demonstration of this convolution correction or a quantitative argument that it is negligible.","tokens_in":10703,"tokens_out":23592,"duration_ms":214951,"concrete_test":"Obtain the GEANT4 energy distributions φ_i(E) at each target position (or re-run the simulation using the reported stack geometry and physics lists). For each of the five reported energies, compute <σ>_i = ∫σ_TALYS(E)φ_i(E)dE / ∫φ_i(E)dE using a TALYS curve that reproduces the measured points (e.g., AOMP3 or AOMP6), and compare <σ>_i with σ_TALYS(Ē_i). If <σ>_i/σ_TALYS(Ē_i) exceeds 1.05 for the 16.03 and 14.09 MeV points, the reported cross sections must be corrected for the bias or the energies re-defined as the yield-averaged effective energy.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The activation yield at each target is the convolution of the (unknown) cross section with the alpha energy distribution at that foil. The paper's GEANT4 simulation produces these distributions, but the analysis uses only the fitted mean energy and its 1σ width (Table 5); the measured activity is assigned to the mean energy Ē, so the reported σ is the yield-averaged cross section <σ> = ∫σ(E)φ(E)dE / ∫φ(E)dE, not σ(Ē). For sub-Coulomb (α,n) reactions, lnσ is approximately πη/E, with η≈0.157 Z1 Z2 sqrt(μ/E); at E=14.09 MeV this gives b=dlnσ/dE≈2.3 MeV⁻¹. For a Gaussian energy distribution with σ_E=0.30 MeV, <σ>/σ(Ē) ≈ exp(½b²σ_E²) ≈ 1.25, a +25% bias. Similar estimates yield +13% at 16.03 MeV, +6% at 17.68 MeV, +4% at 19.34 MeV, and +2.5% at 20.90 MeV. This energy-dependent bias steepens the low-energy tail of the excitation function and is not included in the covariance matrix (Table 8) nor discussed in the text. It directly affects the comparison with Hauser-Feshbach predictions and the astrophysical reaction rate, and it is not covered by the quoted target-thickness uncertainty, making it the most load-bearing unquantified systematic.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports new activation measurements of the 144Sm(α,n)147Gd cross section at five sub-Coulomb energies between 14.09 and 20.90 MeV, using stacked Sm2O3 targets on aluminium backings, GEANT4-simulated beam-energy distributions, offline γ-ray spectroscopy, and a detailed covariance/correlation analysis. The results are compared with earlier measurements and with a 432-combination Hauser-Feshbach spread from TALYS-2.0. The authors state that this is the first covariance analysis for this reaction and that the data are relevant for p-process nucleosynthesis and 147Gd production.","tokens_in":11104,"tokens_out":11136,"duration_ms":116542,"significance":"If the systematic issues below are addressed, the data set fills a gap at sub-Coulomb energies for a p-process nucleus and provides a useful covariance framework for future reaction-rate evaluations. The authors are to be credited for running GEANT4 energy-straggling simulations, applying EFFTRAN summing and extended-geometry corrections, and scanning a wide model space in TALYS. The covariance matrices appear internally consistent with the quoted uncertainty budget. However, the absolute normalization and the treatment of energy spreading over the target are currently not sufficiently quantified for the reported cross sections to be used as final nuclear-data values.","major_comments":[{"comment":"Eq. (3) evaluates σ at the fitted mean energy Ē, but the measured activity is the yield average ⟨σ⟩ = ∫σ(E)φ(E)dE / ∫φ(E)dE over the GEANT4 energy distribution. For a sub-Coulomb (α,n) reaction lnσ ≈ πη/E, so at 14.09 MeV b ≈ 2.3 MeV⁻¹; with the Table 5 width 0.30 MeV this gives ⟨σ⟩/σ(Ē) ≈ 1.25, and ≈ 1.13, 1.06, 1.04, 1.025 at 16.03, 17.68, 19.34, 20.90 MeV. The bias is not included in Table 7 or the covariance matrix of Table 8, and it steepens the low-energy tail in Fig. 9. The authors should correct for this by folding a trial σ(E) through the simulated distributions, or at least quote the model-dependent correction as a systematic uncertainty.","section":"Sec. 3.1 / Table 5"},{"comment":"N_target is never explicitly defined as the areal density of 144Sm nuclei. The targets are 67% enriched Sm2O3; if the weighed 280–350 μg/cm² layer thickness were used directly in Eq. (3), all five cross sections would be incorrect by a factor related to the enrichment and stoichiometry. Because σ ∝ 1/N_target and the target thickness is the dominant 15–25% uncertainty, the paper must state the conversion from mass per area to N_target(144Sm), and include the uncertainties in enrichment, stoichiometry, and the 229Th thickness measurement. This is needed to assess the absolute scale of the data.","section":"Sec. 2.1 / Eq. (3)"},{"comment":"The target-thickness contribution is entered as a purely uncorrelated uncertainty. However, if the 67% enrichment fraction, the Sm2O3 stoichiometry, or the calibration of the 229Th thickness method is common to all five targets, these components are common-mode and should appear in the covariance matrix as a correlated systematic. As written, the reported 7–8% correlations may underestimate the normalization correlations. Separate the thickness uncertainty into random (weighing/deposition) and common (enrichment/stoichiometry/calibration) components.","section":"Sec. 3.3 / Table 7"}],"minor_comments":[{"comment":"The text says 'proton flux' but the beam is 4He2+. Clarify that φ_b is the alpha-particle flux in particles/s, not the electrical beam current, and state how the 5% beam-current uncertainty accounts for charge-state and electron-suppression effects.","section":"Eq. (3)"},{"comment":"Units are inconsistent: Table 1 lists efficiencies as fractions (0.1293), while Table 2 lists values such as 2.694 and 2.786, apparently in percent. Use the same convention and state it explicitly in the captions.","section":"Tables 1 and 2"},{"comment":"The initial beam-energy spread is generated as a uniform random distribution with 0.2 MeV FWHM, but Table 5 reports 1σ uncertainties. State how the 0.2 MeV FWHM was converted to the Gaussian 1σ widths used in the simulation and whether the energy distributions were actually Gaussian.","section":"Sec. 2.3"},{"comment":"The grey band is the max–min envelope of 432 TALYS combinations; this is not a statistical confidence interval. Define its interpretation and, if possible, add a residual or ratio panel to quantify the agreement of the present data with the selected AOMP curves and with the literature data.","section":"Sec. 4 / Fig. 9"},{"comment":"For a nuclear-data measurement paper, raw peak areas, individual target thicknesses, irradiation/counting times, and simulation outputs should be provided in a supplement. The statement that data are 'available upon request' is not sufficient for independent verification.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The measurement is potentially valuable, but the absolute normalization and the energy-yield-averaging bias must be addressed before the data can be used in astrophysical or medical-isotope evaluations. The issues are fixable within the scope of the paper, so I recommend major revision rather than rejection. Please ask the authors to provide the raw data and a target-normalization derivation in the revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The contribution here is concrete: five new sub-Coulomb cross-section points for 144Sm(α,n) and a covariance matrix for this reaction. The previous datasets were sparse, and the covariance treatment is genuinely new for this channel. The activation work is honest and reasonably careful: GEANT4 energy-loss simulation, 152Eu efficiency calibration with summing and extended-geometry corrections, and a broad TALYS comparison with 432 parameter combinations. No circular fitting to the data was done; the TALYS band is independent. Credit is earned for the uncertainty bookkeeping—the correlation matrix looks internally consistent, and the dominant target-thickness uncertainty is stated rather than hidden.\n\nThe main soft spot is the one the stress-test note identifies. The activity is assigned to the mean beam energy, but activation integrates the cross section over the energy distribution. At sub-Coulomb energies the cross section rises steeply; with 0.30 MeV straggling at the lowest point, the bias is plausibly +25% at 14.09 MeV, and still +10% or more at the next point. That bias is not in the covariance matrix and is not discussed. It steepens the low-energy tail and directly affects the comparison with Hauser-Feshbach predictions. It is not fatal, because the GEANT4 simulation yields the full energy distributions and the authors could correct for the effect or propagate it; but as published it is a real unquantified systematic.\n\nThe target-thickness uncertainty (15–25%) is the other major concern. It is honestly reported, but the determination by weighing and 229Th alpha source is described in only a few lines, and the paper does not include raw data or simulation outputs. The standard \"available upon reasonable request\" statement is acceptable but makes independent verification harder. Minor: Eq. (3) calls φb the proton flux, which should be alpha flux.\n\nThis paper deserves a serious referee. The new data and covariance analysis are useful for p-process network calculations and for 147Gd production estimates, and the experimental work is otherwise sound. I would send it to review with a request that the authors either apply the energy-distribution correction or justify why it is negligible, and that they make the energy distributions available. If they do that, the paper is a solid incremental contribution.","headline":"New sub-Coulomb 144Sm(α,n) cross sections with a careful covariance treatment; the energy-binning bias at the lowest points is real and should be addressed before publication.","tokens_in":11564,"tokens_out":1981,"would_cite":false,"duration_ms":21524,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"New cross sections for 144Sm(alpha,n)147Gd at five sub-Coulomb energies are reported, together with the first complete covariance analysis for this reaction, giving data relevant to p-process nucleosynthesis and the medical isotope 147Gd.","keywords":["144Sm(alpha,n)147Gd","sub-Coulomb cross section","activation method","covariance analysis","p-process nucleosynthesis","Hauser-Feshbach model","stacked foil technique","147Gd SPET isotope"],"falsifier":"Re-measure the target thicknesses with an independent method (e.g., ion-beam backscattering) or re-measure the cross sections at the same five energies using a different target preparation; if the new thickness values fall outside the reported 15–25% range, or if an independent cross-section point at 14.1 MeV differs from 0.79 mb by more than 0.20 mb, the absolute scale of the data would be in question.","tokens_in":10621,"feed_emoji":"⚛️","tokens_out":17350,"duration_ms":125869,"temperature":0.7,"pith_summary":"The paper seeks to establish trustworthy cross-section values for the reaction 144Sm(alpha,n)147Gd at five alpha energies between 14 and 21 MeV, a sub-Coulomb regime where measurements are hard and data are sparse. Using the stacked-foil activation technique, it obtains cross sections from about 595 millibarns at 20.9 MeV down to 0.79 millibarns at 14.1 MeV, and it provides the first complete covariance and correlation matrix for this reaction, showing the points are mutually correlated at the 7–8% level. Such data feed directly into stellar models of heavy-element production (the p-process) and into estimates of 147Gd yields for a promising medical imaging isotope. The paper also compares the measurements with a wide spread of Hauser-Feshbach statistical-model predictions, finding that the model band brackets the data and that the choice of alpha-nucleus optical potential matters most.","feed_headline":"Five nuclear cross-section points measured below the Coulomb barrier","feed_subtitle":"For the Sm-144(alpha,n) reaction, with full error correlations for stellar models.","key_machinery":"The carrying mechanism is the stacked-foil activation measurement combined with a covariance analysis. A single 28 MeV alpha beam is sent through a series of thin, isotopically enriched 144Sm2O3 targets separated by aluminium degraders, so that one irradiation produces five different sub-Coulomb bombarding energies. A Monte Carlo simulation of energy loss and straggling through the degraders and targets fixes the mean energy and its 1σ spread at each foil. Each cross section is then derived from the offline gamma-ray activity of the 147Gd product (mainly the 229.3 keV line) using the activation formula, and the uncertainties of all inputs — beam current, detector efficiency, gamma-ray intens","core_discovery":"On its own terms, the paper reports five absolute cross sections for 144Sm(alpha,n)147Gd below the Coulomb barrier: 594.61 ± 154.38 mb at 20.90 ± 0.18 MeV, 431.44 ± 112.03 mb at 19.34 ± 0.20 MeV, 124.53 ± 32.45 mb at 17.68 ± 0.22 MeV, 20.52 ± 5.35 mb at 16.03 ± 0.27 MeV, and 0.79 ± 0.20 mb at 14.09 ± 0.30 MeV. The measurements were obtained by degrading one 28 MeV alpha beam through a stack of aluminium foils and five thin 144Sm2O3 targets, with the effective energy at each target determined by Monte Carlo simulation. The paper's methodological claim is that it is the first measurement of this reaction to include a full covariance analysis, so that the uncertainties (dominated by target thic","pith_inferences":["The authors' 15–25% target-thickness uncertainty is the largest single contribution; an independent thickness measurement (for example by backscattering spectrometry) would directly test the absolute scale of all five cross sections, since the cross section is inversely proportional to target thickness.","The reported ~7–8% correlations imply that ignoring correlations when these points feed into a reaction network would slightly understate the uncertainty on the derived stellar reaction rate; the provided matrix allows that error to be handled exactly.","If the two lowest-energy points are confirmed by an independent measurement, they could become benchmarks for the alpha optical potential in a regime where theoretical predictions are strongly divergent.","The paper's use of a single 28 MeV beam and Monte Carlo energy determination could be validated by directly measuring the degraded beam energy at each target position with a detector, which would check the simulation's energy-loss and straggling model."],"forward_implications":["The five cross-section points can be inserted directly into stellar nucleosynthesis networks, updating the production/destruction balance for 144Sm and neighboring p-nuclei.","Because the full covariance matrix is given, reaction-rate calculations for this channel can correctly propagate the common systematic components (beam flux, efficiency, gamma intensity) instead of assuming the points are independent.","The finding that the alpha optical-model potential dominates the theoretical spread suggests that these data can serve as a new constraint on alpha-nucleus potentials at sub-Coulomb energies.","For applied purposes, the cross sections give a quantitative basis for predicting 147Gd yields from alpha irradiation of enriched 144Sm, relevant to SPET imaging.","The two lowest-energy points, where the cross section drops steeply, provide a sharp test for any statistical-model calculation in the far-sub-Coulomb regime."],"fun_headline_variants":["Sub-Coulomb Sm-144(alpha,n) cross sections with full covariance","Five sub-Coulomb cross sections for Sm-144(alpha,n) with correlated errors","First full covariance analysis for Sm-144(alpha,n) below barrier","Sm-144(alpha,n) measured below Coulomb barrier with error correlations","New sub-Coulomb data for Sm-144(alpha,n) with full uncertainty matrix"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire normalization of the five cross sections rests on the measured thickness (areal density) of the thin 144Sm2O3 targets; if that measurement is systematically off by some factor, every reported cross section is off by the same factor.","fun_headline_variants_meta":{"raw":{"variants":["Sub-Coulomb Sm-144(alpha,n) cross sections with full covariance","Five sub-Coulomb cross sections for Sm-144(alpha,n) with correlated errors","First full covariance analysis for Sm-144(alpha,n) below barrier","Sm-144(alpha,n) measured below Coulomb barrier with error correlations","New sub-Coulomb data for Sm-144(alpha,n) with full uncertainty matrix"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000154,"raw_usage":{"total_tokens":1062,"prompt_tokens":770,"completion_tokens":292,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":514,"completion_tokens_details":{"reasoning_tokens":189}},"tokens_in":514,"tokens_out":292,"duration_ms":3242,"temperature":1.0,"reasoning_tokens":189,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T21:08:35.117846+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the target thicknesses with an independent method (e.g., ion-beam backscattering) or re-measure the cross sections at the same five energies using a different target preparation; if the new thickness values fall outside the reported 15–25% range, or if an independent cross-section point at 14.1 MeV differs from 0.79 mb by more than 0.20 mb, the absolute scale of the data would be in question.","supporting_citations":[],"review_version":1}