{"id":"973fc2ce-386b-490d-9d80-98ef58d188d3","arxiv_id":"2501.00963","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For fusion of 28,30,32Si with 28Si, measured proton and alpha multiplicities and energies show GEMINI++ overpredicts proton emission and that enhanced first-step alpha emission improves agreement.","lead":"This paper measures protons and alpha particles emitted after fusing silicon nuclei of three different masses, and compares the results to a statistical decay model. The comparison shows the model gets proton emission wrong in a systematic way and suggests alpha particles are emitted earlier in the decay cascade than the model assumes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model-dependent efficiency corrections may bias the multiplicity comparison in the same direction as the claimed early-alpha deficit; the residual systematic is not quantified.","rationale":"The reader's weakest-assumption analysis correctly identifies the efficiency circularity as the load-bearing issue: the multiplicity comparison and the emission-order inference depend on Eq. (1), while the energy-spectrum shape comparison (Fig. 9) is more robust but the headline conclusion combines both. The paper is honest about the assumption and provides some validation, but no error budget is attached to it. Since the reader already conditions the verdict on quantifying this systematic, the conditional verdict remains appropriate; no adjustment is needed.","tokens_in":14982,"tokens_out":4314,"duration_ms":43706,"concrete_test":"Re-extract <M_p> using efficiencies from a modified GEMINI++ sample in which the first-step alpha fraction is increased to match the best description in Fig. 15 (e.g., by weighting first-alpha post-selected cascades to reproduce the desired fractional yield), while keeping Eq. (1) fixed otherwise. If <M_p> shifts by more than the quoted statistical uncertainty, the central conclusion is not robust to the model-dependence; if it does not, the circularity concern is settled. A useful cross-check is also to recompute epsilon_ER using the measured ER angular distribution as a reweighting function.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative step that turns the data into the central claim is Eq. (1): <M_i> = (N_{i-ER}/N_{ER}) * (epsilon_ER/epsilon_{i-ER}), with both efficiencies computed from default GEMINI++ cascades (Section II). If GEMINI++ under-represents early alpha emission—the hypothesis being advanced—then its ER recoil distribution (Fig. 5d) and its proton/alpha angular distributions are biased; epsilon_ER, epsilon_p-ER, and epsilon_alpha-ER are then all evaluated with the wrong cascade. The paper's 'reasonableness' check compares energy distributions only over the HiRA acceptance (Figs. 9-10), which does not validate the model angular distributions needed for the acceptance correction, especially since the detector measures only a few angles. Consequently, the multiplicity comparison in Fig. 8 is not fully independent of the model it is meant to test: part of the apparent proton overprediction could be an artifact of efficiencies computed from a cascade that already lacks the proposed early-alpha enhancement. The authors acknowledge the assumption but do not quantify the resulting systematic error on <M_p> or <M_alpha>.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a high-statistics measurement of proton and α-particle emission following fusion of 28,30,32Si with a 28Si target, using ER identification via energy–time-of-flight and LCP identification via ΔE–E. The authors extract average proton and α multiplicities through Eq. (1) using GEMINI++-computed geometric efficiencies, compare the multiplicities and energy spectra with GEMINI++ predictions, and find that GEMINI++ overpredicts proton multiplicities and average energies while describing α observables reasonably well for 28,30Si. They then use GEMINI++ sensitivity studies in which the identity of the first emitted particle is constrained to argue that increasing the fractional yield of early α emission improves agreement for protons, with increased early neutron emission also indicated for the neutron-rich 32Si system. The central claim is that joint measurements of particle energy distributions, multiplicities, and ER angular distributions can reveal details of the particle emission sequence.","tokens_in":15152,"tokens_out":2567,"duration_ms":26053,"significance":"If the central claim holds, the paper offers a new observable-driven handle on the de-excitation cascade of compound nuclei, with potential relevance for statistical-model descriptions of fusion–evaporation and for establishing a baseline before probing dynamical α-cluster emission. The experimental methods are standard and generally carefully described, with clean particle identification and three-isotope systematics that allow a meaningful test of neutron-excess dependence. The paper is honest about its main model-dependent assumption and presents explicit sensitivity studies rather than over-claiming direct measurement of emission order. The work is a solid contribution to the fusion–evaporation literature, although the load-bearing efficiency correction and the model-selection nature of the emission-order test require further quantitative scrutiny.","major_comments":[{"comment":"The multiplicity extraction is not independent of the model being tested: both ε_ER and ε_i−ER in Eq. (1) are computed from default GEMINI++ cascades. If the true cascade differs from GEMINI++ in the direction argued in the paper — i.e., under-represented early α emission — the angular and energy distributions entering these efficiencies are biased, and part of the apparent proton overprediction in Fig. 8 could be an artifact of the efficiency correction. The text acknowledges this assumption ('Use of GEMINI++ to determine the geometric efficiency ... assumes ...') but does not quantify the resulting systematic uncertainty on ⟨M_p⟩ or ⟨M_α⟩. The authors should provide a quantitative estimate of this systematic, for example by recomputing the efficiencies with modified first-step emission probabilities or with an alternative statistical-decay code, and show how the multiplicities in Fig. 8 shift.","section":"§II, Eq. (1)"},{"comment":"The 'reasonableness' check used to support the efficiency-correction assumption compares only energy distributions of protons and α-particles over the limited HiRA angular acceptance, with each GEMINI++ spectrum normalized to the data integral. This does not validate the angular distributions of the emitted particles, which are precisely the ingredient needed for the coincidence efficiencies ε_p−ER and ε_α−ER and which cannot be directly measured with the present few-angle HiRA setup. The paper should state this limitation explicitly and estimate the sensitivity of the extracted multiplicities to plausible changes in the model angular distributions, for example by comparing GEMINI++ angular distributions with any existing data for similar compound systems.","section":"§IV, Figs. 9–10"},{"comment":"The central conclusion — that enhanced early α emission improves agreement with experimental proton multiplicities and energies — is derived from GEMINI++ calculations in which the identity of the first emitted particle is externally constrained. This is a valuable sensitivity study, but it does not constitute a direct measurement of emission order, and the conclusion is conditional on the model's internal correlations between emission step, particle energy, and multiplicity. In particular, the authors note that increased initial α emission suppresses low-energy α yield and thus worsens the description of α energy distributions, so the suggested modification is not globally consistent. The abstract and conclusions claim that the data allow one to 'deduce interesting details of the de-excitation cascade'; the authors should either temper this claim or provide a quantitative, self-consistent estimate of how much early-α enhancement is required and whether it can simultaneously reproduce the proton and α observables within uncertainties.","section":"§V, Figs. 14–16"}],"minor_comments":[{"comment":"The detector designation appears as 'M CPU S' in several places; it should be rendered consistently as, e.g., MCP_US.","section":"§II (Fig. 1 caption and text)"},{"comment":"There is a typo: 'GEMINII++' should be 'GEMINI++'.","section":"§IV, Fig. 14 caption"},{"comment":"The label 'dahsed' should be 'dashed'.","section":"§III, Fig. 6 caption"},{"comment":"The sentence 'the GEMINI++ model provides a reasonable description the ⟨Mα⟩' is missing the word 'of' before 'the ⟨Mα⟩'.","section":"§V, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"The experimental data and analysis are of good quality and the paper is within scope for a nuclear physics journal. The main concern is the circularity between the GEMINI++-based efficiency correction and the GEMINI++ comparison used to draw the emission-order conclusion; this is acknowledged but not quantified. I would encourage the editor to request that the authors provide a quantitative systematic estimate for the efficiency-correction model dependence, as this is load-bearing for the central claim. If they can do so, the paper would be suitable for publication after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper for the data: it adds the radioactive 32Si + 28Si system to the 28,30Si fusion-evaporation dataset and makes a clean three-isotope comparison of proton and alpha multiplicities and energy spectra across neutron excess. The experiment is careful, the PID is clean, and the ER angular distributions (Fig. 5d) are a genuinely nice way to show how the balance of n/p/alpha emission is encoded in recoil kinematics. That part is good, solid experimental work.\n\nThe central claim—that GEMINI++ overpredicts proton multiplicities and average energies, and that boosting early alpha emission in the model fixes much of the discrepancy—is plausible but not as firmly anchored as it looks. The efficiency correction in Eq. (1) comes from GEMINI++ itself, and the same model is then compared to the corrected data. The paper acknowledges this and checks “reasonableness” by comparing energy spectra, but that check only covers shapes within the HiRA acceptance. It does not validate the model angular distributions that drive the acceptance correction. The stress-test note is right: if the model under-represents early alphas, the efficiencies are evaluated with the wrong cascade, and the bias could plausibly run in the same direction as the claimed deficit. I think that concern is legitimate and should be quantified before the quantitative multiplicities are taken at face value.\n\nThe emission-order conclusion is also a sensitivity study, not a measurement. The authors are honest about that—they explicitly say the forced-emission cases are not meant to reproduce the data—but they then use those cases to infer that early alpha emission is underrepresented. That is a hypothesis, not a determination. Still, the qualitative trend (GEMINI++ overpredicts protons more for the neutron-deficient systems) is large enough that it probably survives even a moderate efficiency bias.\n\nMinor points: the model spectra in Figs. 9–10 are normalized to the data integrals, which discards multiplicity information; the A/6 level-density test is only one knob and doesn't fix the proton overprediction. These are secondary.\n\nBottom line: the data are valuable, the paper is readable, and the interpretation is clearly flagged as model-based. The missing piece is a systematic uncertainty on the efficiency correction—e.g., varying the input angular/energy distributions and recomputing the efficiencies. That is addressable in revision. I'd send it to a serious referee; it deserves the time, and the referee should push on exactly that point.","headline":"Solid new three-isotope fusion-evaporation data, but the emission-order conclusion rests on a model-dependent efficiency correction whose systematic uncertainty is unquantified.","tokens_in":15789,"tokens_out":2043,"would_cite":true,"duration_ms":21871,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.70.Jj","24.60.Dr"],"model":"deepseek-v4-flash","headline":"A joint measurement of protons, alpha particles, and evaporation residues from 28,30,32Si + 28Si fusion shows that the order of particle emission is observable, and that statistical decay models need more alpha particles in the first…","keywords":["fusion-evaporation","GEMINI++","statistical decay","particle emission sequence","alpha clustering","light-charged particles","evaporation residues","neutron excess"],"falsifier":"Run GEMINI++ with a tunable fraction of cascades forced to begin with an α-particle and fit that fraction to the measured proton average energy and multiplicity for all three systems; if no single such fraction reproduces all of the proton and α observables simultaneously, the proposed emission-order modification is falsified.","tokens_in":14763,"feed_emoji":"⚛️","tokens_out":8887,"duration_ms":80108,"temperature":0.7,"pith_summary":"This paper reports a coincidence measurement of protons and α-particles emitted when silicon-28, silicon-30, and silicon-32 beams fuse with a silicon-28 target, together with the evaporation residues left behind. The authors compare the measured particle energies and multiplicities with predictions of the statistical decay code GEMINI++ and find that the model systematically overpredicts the proton multiplicity and average proton energy for all three systems, while describing α-particle energies well. Their central claim is that the de-excitation cascade is not fully captured by the default model: increasing the fractional yield of α-particle emission in the initial step brings the model's proton predictions into better agreement with the data, and for the most neutron-rich system increased initial neutron emission is also indicated. This matters because it shows that the order, not just the amount, of particle emission carries measurable information, and it provides a specific target for improving statistical decay models used to interpret fusion and cluster-emission experiments.","feed_headline":"Statistical decay code emits too few early alphas","feed_subtitle":"Proton energies and multiplicities from 28,30,32Si + 28Si fusion line up only when first-step alpha emission is boosted.","key_machinery":"The argument is carried by the coincidence measurement and by an emission-order selection applied inside GEMINI++. Experimentally, evaporation residues are identified by energy–time-of-flight in annular silicon detectors, light-charged particles by the $\\Delta E$–$E$ technique, and multiplicities are extracted from coincidences via $\\langle M_i\\rangle = (N_{i-\\mathrm{ER}}/N_{\\mathrm{ER}})(\\epsilon_{\\mathrm{ER}}/\\epsilon_{i-\\mathrm{ER}})$, with GEMINI++ supplying the geometric efficiencies. The model is a Hauser–Feshbach statistical decay code that simulates the de-excitation of a compound nucleus as a sequence of binary particle emissions. The paper then deconstructs the cascade by forcing the first emitted particle to be a neutron, proton, or α-particle and comparing those selected sub-cascades with the inclusive prediction; these selections reveal that first-step α-emission lowers the average proton energy and brings proton multiplicities down toward the measured values.","core_discovery":"The paper establishes that the order in which protons, neutrons, and α-particles are emitted from a fused compound nucleus is observable, and that the standard statistical decay code GEMINI++ does not get that order right. For fusion of 28,30,32Si with 28Si, the measured proton multiplicities and average energies lie below GEMINI++ predictions for all three systems, while the measured α-particle energies are well described and the model underpredicts the yield of low-energy protons. By examining GEMINI++ decays selected on the identity of the first emitted particle, the authors show that an α-particle emitted first lowers the excitation and barrier for subsequent proton emission, enhancing low-energy protons and reducing the proton multiplicity; increasing the fractional yield of α-particle emission in the initial step brings the model's proton $\\langle E_{\\mathrm{lab}}\\rangle$ and $\\langle M\\rangle$ into better agreement with the data. For the most neutron-rich system, 32Si, the data additionally indicate increased initial neutron emission. The paper concludes that joint measurement of particle energy distributions, multiplicities, and evaporation-residue angular distributions constrains the emission sequence, and that early α-emission (and early neutron emission in neutron-rich systems) is under-represented in the default model.","pith_inferences":["Editorial extension: if early α-emission is generally under-represented in GEMINI++, then experiments that use GEMINI++ to compute α-particle detection efficiencies—including the efficiency corrections in this paper—may carry a systematic bias whose size depends on how the α-first fraction varies with neutron excess.","Editorial extension: the inference that early neutron emission increases for 32Si could be checked directly with a neutron detector; the paper's mechanism predicts a higher average neutron energy for 32Si relative to the default GEMINI++ calculation.","Editorial extension: a systematic scan of projectile–target combinations spanning a wider range of N/Z could test whether the required α-first fraction tracks the ratio of α to neutron separation energies, which would tie the emission-order effect to the same physics that governs cluster formation in neutron-rich matter."],"forward_implications":["The proton overprediction in GEMINI++ is not fixed by changing the level-density parameter from A/7 to A/6; the A/6 choice lowers proton and α energies slightly but makes the proton-multiplicity discrepancy worse.","For 28,30Si, enhanced first-step α-emission should reduce the model's underprediction of evaporation-residue yield at $\\theta_{\\mathrm{lab}} > 12^\\circ$, since the resulting higher-energy α-particles impart a larger transverse recoil to the residue.","For 32Si, the data indicate increased initial neutron emission as well, and the paper notes that changing the emission sequence alone is not sufficient for α multiplicities, suggesting nucleon emission relative to α emission is also under-represented.","Joint measurement of light-charged-particle energy distributions, multiplicities, and evaporation-residue angular distributions constitutes a practical way to constrain emission order; measuring neutron kinetic energies and multiplicities would further constrain the sequence."],"supporting_citations":[{"why":"Provides the GEMINI++ statistical decay code, the model whose default predictions are compared with the measured multiplicities, energies, and evaporation-residue angular distributions, and whose emission-order selections are used to test the alpha-first hypothesis.","marker":"[22]"},{"why":"Records the prior observation of a 4–5 fold excess of α-particle yield in 18O+12C fusion relative to statistical-model calculations, motivating the need to understand standard statistical α emission.","marker":"[9]"},{"why":"Supplies the standard statistical-model calculations used as the comparison baseline for that earlier α-excess observation.","marker":"[10]"},{"why":"Presents TDDFT/NLF calculations showing dynamically generated α-cluster states during fusion, framing the broader question of whether α emission is dynamical or statistical that this paper's statistical-emission measurement is meant to sharpen.","marker":"[13]"},{"why":"The Bass fusion model supplies the maximum angular momentum ℓmax used as input to GEMINI++, affecting the decay cascade and hence the predicted multiplicities and energies.","marker":"[23]"},{"why":"Provides the A/6 level-density parametrization tested as an alternative to the default A/7; this alternative lowers α multiplicities but does not fix the proton overprediction.","marker":"[25]"}],"fun_headline_variants":["Early alpha emission key to matching proton data","Mismatch in proton yields traces to early alphas","Statistical code needs more first-step alphas","Fusion decay ordering exposed by proton data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that GEMINI++'s simulated angular and energy distributions of emitted particles are accurate enough to compute the detector efficiencies, since the same model is then compared with the efficiency-corrected data.","fun_headline_variants_meta":{"raw":{"variants":["Early alpha emission key to matching proton data","Mismatch in proton yields traces to early alphas","Statistical code needs more first-step alphas","Fusion decay ordering exposed by proton data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000412,"raw_usage":{"total_tokens":2129,"prompt_tokens":940,"completion_tokens":1189,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":1131}},"tokens_in":556,"tokens_out":1189,"duration_ms":9169,"temperature":1.0,"reasoning_tokens":1131,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:37:45.428736+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run GEMINI++ with a tunable fraction of cascades forced to begin with an α-particle and fit that fraction to the measured proton average energy and multiplicity for all three systems; if no single such fraction reproduces all of the proton and α observables simultaneously, the proposed emission-order modification is falsified.","supporting_citations":[{"cited_title":"Charity, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the GEMINI++ statistical decay code, the model whose default predictions are compared with the measured multiplicities, energies, and evaporation-residue angular distributions, and whose emission-order selections are used to test the alpha-first hypothesis."},{"cited_title":"Vadas, T","cited_arxiv_id":null,"evidence_quote":"Records the prior observation of a 4–5 fold excess of α-particle yield in 18O+12C fusion relative to statistical-model calculations, motivating the need to understand standard statistical α emission."},{"cited_title":"Gavron, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the standard statistical-model calculations used as the comparison baseline for that earlier α-excess observation."},{"cited_title":"Schuetrumpf and W","cited_arxiv_id":null,"evidence_quote":"Presents TDDFT/NLF calculations showing dynamically generated α-cluster states during fusion, framing the broader question of whether α emission is dynamical or statistical that this paper's statistical-emission measurement is meant to sharpen."},{"cited_title":"Bass, Nucl","cited_arxiv_id":null,"evidence_quote":"The Bass fusion model supplies the maximum angular momentum ℓmax used as input to GEMINI++, affecting the decay cascade and hence the predicted multiplicities and energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the A/6 level-density parametrization tested as an alternative to the default A/7; this alternative lowers α multiplicities but does not fix the proton overprediction."}],"review_version":1}