{"id":"1abfd703-8652-4c56-8181-3a14d8c137d5","arxiv_id":"1908.03429","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A theoretical model computes the Fermi momentum distribution and in-medium branching ratios for the non-mesonic decay of an eta-mesic helium-3 bound state through eta to 2gamma and eta to 3pi0 channels.","lead":"This paper models what happens when an eta meson bound to a helium-3 nucleus decays directly into photons, rather than being absorbed by the nucleus. It provides momentum distributions and decay branching ratios that experimenters need for Monte Carlo simulations of the WASA-at-COSY search for eta-mesic helium.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central numbers are an unconstrained function of the assumed (V0,W0) optical-potential parameters; without data constraints or uncertainty estimates, the claimed MC inputs are not robust.","rationale":"The reader's weakest-assumption analysis correctly identifies that the optical-potential parameters (V0,W0) control every quantitative output, and that these parameters are not constrained by data in the paper. My stress-test reaches the same conclusion: Table 1 shows branching ratios spanning roughly a factor of 27 across the four parameter sets, and the abstract presents these as decisive MC inputs without recommending a set or providing an error band. This makes the central claim—that the model provides the kinematic template and branching ratios crucial for experimental data interpretation—conditional on an unvalidated choice of potential. The reader's CONDITIONAL verdict is therefore appropriate; no fatal internal contradiction exists, since solving Eq. (1) for a given potential is standard and the parameter scan is transparent. I do not see a reason to move the verdict to ACCEPT or REJECT, so the verdict remains UNCHANGED. The secondary points (the asserted FSI half-reduction and the use of vacuum partial widths in Eq. (4)) reinforce the need for qualification but are less load-bearing than the unresolved parameter dependence.","tokens_in":7065,"tokens_out":5088,"duration_ms":62008,"concrete_test":"Constrain (V0,W0) by fitting the WASA-at-COSY pd -> 3He eta excitation function data (Refs. [56,57]) in the same Klein-Gordon framework, or by using the eta-4He constraint from Ref. [7] to set the plausible parameter range. Then compute the Bayesian posterior over (V0,W0) and propagate it to BR*(eta->2gamma) and BR*(eta->3pi0). If the posterior span remains larger than a factor of about 3, the paper should present the branching ratios as a continuous band with confidence intervals rather than four discrete values, and the Monte Carlo simulations should be run over that band to quantify the resulting uncertainty in detection efficiency.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative content of the paper—binding energies, absorption widths, the relative 3He-eta momentum distribution in Fig. 2, and the in-medium branching ratios in Table 1—is obtained by solving Eq. (1) with four hand-picked optical-potential parameter sets, (V0,W0)=(-75,-20), (-90,-20), (-75,-1), (-90,-1) MeV. No experimental input constrains these parameters in the 3He-eta system; the weakly absorptive value W0=-1 MeV is imported from the 4He-eta analysis of Ref. [7], while no analogous constraint is provided for the real part. The resulting branching ratios vary by more than an order of magnitude, from 2.5e-5 to 6.8e-4 for eta->2gamma and from 2.1e-5 to 5.6e-4 for eta->3pi0, i.e. a factor of about 27 depending on the assumed potential. The abstract nevertheless presents 'in-medium branching ratios' as crucial inputs for Monte Carlo simulations and experimental data interpretation, without stating which parameter set is preferred or attaching any uncertainty to the values. Since the experimental search is designed to test the existence and properties of the eta-3He bound state, using any one of these unconstrained numbers as a simulation template cannot yield robust efficiency estimates or reliable interpretation of a null or positive signal. The paper is transparent in showing the four parameter sets, but it does not offer a procedure to discriminate among them, so the central deliverable is a scan over an unconstrained model parameter rather than a predictive result. A secondary, related weakness is that the claimed final-state-interaction reduction of the signal 'to half' is asserted without derivation or reference, adding another uncontrolled factor of two to the branching-ratio estimates.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a theoretical model for the non-mesonic decay of an η-mesic 3He bound state produced in pd collisions. The authors solve the Klein-Gordon equation with an optical potential U_opt(r) = (V0 + iW0) ρ(r)/ρ0 for four parameter sets, obtaining binding energies, absorption widths, the relative 3He-η momentum distribution, and in-medium branching ratios BR* for η→2γ and η→3π0. They then describe a Monte Carlo framework, using these distributions to estimate WASA-at-COSY detection acceptances for pd→3He2γ and pd→3He6γ. The stated goal is to provide crucial inputs for experimental data interpretation in the search for η-mesic 3He.","tokens_in":7361,"tokens_out":2712,"duration_ms":29616,"significance":"If the computed momentum distributions and in-medium branching ratios were robust, the paper would fill a genuine gap: previous simulations assumed the N* momentum distribution equals the nucleon distribution, while the present model derives the bound η momentum from a quantum mechanical wavefunction. The paper also extends the study of η-mesic decays to the non-mesonic channels, for which the WASA-at-COSY measurement exists. However, the deliverable is a scan over four unconstrained optical-potential parameter sets with no data constraint and no uncertainty quantification; the resulting branching ratios vary by a factor of about 27. The central quantitative claims are therefore only as reliable as the assumed (V0,W0), and the paper does not provide the reader with a preferred set or a procedure to discriminate among them. The manuscript is transparent about the parameter dependence, but it does not currently support the abstract's assertion that these are 'crucial' inputs without qualification.","major_comments":[{"comment":"All quantitative results (binding energy, absorption width, momentum distribution, branching ratios) are obtained by solving Eq. (1) with four hand-picked (V0,W0) sets, no experimental input constrains these parameters for the 3He-η system, and the weakly absorptive value W0=-1 MeV is imported from the 4He-η analysis of Ref. [7]. The resulting BR* values span 2.5e-5 to 6.8e-4 for η→2γ and 2.1e-5 to 5.6e-4 for η→3π0, i.e. a factor of about 27. The abstract presents these as 'crucial' for Monte Carlo simulations and data interpretation without stating which parameter set is preferred or attaching any uncertainty. Since the experimental search itself aims to test the existence and properties of the η-3He bound state, using any one of these unconstrained numbers as a simulation template does not yield robust efficiency estimates or reliable interpretation of a null or positive signal. The authors should either constrain the potential from existing data (e.g. threshold scattering, dd/ pd reactions), or clearly frame the results as an illustrative sensitivity study with explicit caveats and uncertainties.","section":"Section 2.1, FSI statement"},{"comment":"The sentence 'We found that the FSI reduces the strength of the signal to half' is a load-bearing correction to the branching ratios, but no derivation, model, or reference is given. This factor directly multiplies the predicted signal in the Monte Carlo efficiency estimate, so its origin and uncertainty must be documented. If it is a qualitative estimate, that should be stated explicitly; if it is based on a calculation, the calculation should be presented or cited.","section":"Section 2.1, Eq. (4) and Table 1"},{"comment":"Equation (4) uses the in-vacuum partial widths Γ_η→2γ and Γ_η→3π0 in the numerator and only modifies the denominator by adding Γ_abs. Labeling the result as an 'in-medium branching ratio' is misleading if the decay widths themselves are modified in the nuclear medium (e.g. through in-medium η properties or N*(1535) dynamics). The authors should either justify that the in-medium modification of the partial widths is negligible for these channels, or include the modification in the calculation. As written, BR* is a vacuum-branching fraction suppressed by absorption, not a fully in-medium branching ratio.","section":"General"},{"comment":"The manuscript contains no uncertainty estimates, either statistical or systematic, for any of the reported numbers (binding energies, widths, momentum distributions, branching ratios, acceptances). Given that the primary input parameters are unconstrained, the absence of a sensitivity analysis leaves the central deliverable, the Monte Carlo template, without a documented range of validity. At minimum the authors should provide a quantitative statement of how the acceptance and efficiency estimates depend on the parameter set.","section":"General"}],"minor_comments":[{"comment":"The notation '-(75,20) MeV' is ambiguous; it should be written as (V0,W0)=(-75,-20) MeV to avoid confusion between the sign of the real and imaginary parts.","section":"Table 1"},{"comment":"The normalization condition ∫|R(p)|²p²dp=1 is stated, but the figure caption for Fig. 2 does not explicitly state whether the plotted quantity is |R(p)|²p² or |R(p)|²; the text says |R(p)|²p², which is consistent, but the y-axis label '|R(p)|²p² [MeV⁻¹]' is unconventional and should be clarified.","section":"Section 2.1, Eq. (3)"},{"comment":"The acceptance values (60% and 40%) are quoted without specifying the excess energy Q, the optical-potential parameter set, or the analysis cuts used in the Monte Carlo; these details are needed for reproducibility.","section":"Section 2.2"},{"comment":"Ref. [17] (A. Fix et al., Phys. Lett. B 772, 663 (2017)) and Ref. [7] (M. Skurzok et al., Phys. Lett. B 772, 663 (2018)) appear to refer to the same volume and page; please verify the citations and ensure they are correctly distinguished.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is essentially a model-input paper for an experimental search by the same collaboration. The unconstrained parameter scan is transparent, but the central claim that these are 'crucial' MC inputs is not supported without either data constraints or uncertainty quantification. A major revision that reframes the results as a sensitivity study, adds a derivation or reference for the FSI factor, and discusses in-medium width modifications would make the paper publishable. If the journal's scope prioritizes experimental analysis tools, the authors might also consider submitting a shorter technical note once the robustness issues are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Know this paper as a model calculation feeding the WASA-at-COSY search: it supplies the first in-medium branching ratios for eta->2gamma and eta->3pi0 from a bound 3He-eta state, plus a relative momentum distribution to replace the old assumption that the eta momentum follows the nucleon distribution. That is genuinely useful for Monte Carlo efficiency estimates, and the authors are upfront that they scan four optical-potential parameter sets rather than pretending to one answer.\n\nThe Klein-Gordon treatment is standard and the numerical work looks consistent: binding energies, absorption widths, and momentum distributions follow from solving Eq. (1) with the stated (V0,W0). Presenting the resulting branching ratios as a table with four cases is honest. The efficiency example for WASA is a sensible illustration of how the template would be used.\n\nThe soft spots are concentrated in what the paper claims versus what it actually delivers. Eq. (4) is a simple ratio of vacuum partial width to total plus absorption width; calling the result an \"in-medium branching ratio\" overstates it, because the numerator is not modified in medium. More importantly, the quantitative outputs—binding energy, width, momentum distribution, and the branching ratios—are all consequences of the assumed optical potential, and the paper provides no data constraint for the 3He-eta system. The weakly absorptive W0=-1 value is borrowed from the 4He-eta analysis, and nothing distinguishes among the four parameter sets. Since the branching ratios span more than an order of magnitude (factor ~27), any single number used as an MC input carries an uncontrolled uncertainty. The paper would be stronger if it presented results as a parameter scan with explicit caveats, or if it could point to any observable that would discriminate among the sets.\n\nThe \"FSI reduces the signal to half\" statement is also asserted without derivation or reference. That is a small but real gap, because it directly changes the expected yields. No error bars anywhere is a minor issue for a calculation like this, but worth noting.\n\nThe citation pattern is appropriate: prior optical-model work and the N* momentum model are cited, and the self-citations are to the relevant WASA analyses. I don't see a fatal contradiction, and the paper is transparent about its inputs. The central argument—that these channels need a dedicated kinematic model—holds up.\n\nWho benefits: experimentalists analyzing the WASA pd data and theorists working on eta-mesic nuclei. It deserves a serious referee; the referee should push for a clearer statement that these are model-dependent template values, not predictions, and for either a derivation of the FSI factor or removal of that factor pending calculation.\n\nRecommendation: send to review, conditional on revisions that reframe the branching ratios as parameter-dependent estimates and address the FSI factor.","headline":"Useful kinematic template for the WASA search, but its central branching ratios are unconstrained products of assumed optical-potential parameters and the FSI factor is asserted, not derived.","tokens_in":8018,"tokens_out":1955,"would_cite":false,"duration_ms":20019,"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":"The paper models the non-mesonic decay of the $\\eta$-mesic $^{3}\\mathrm{He}$ bound state, computing the relative momentum distribution of the $^{3}\\mathrm{He}$-$\\eta$ pair and the in-medium branching ratios for $\\eta\\to 2\\gamma$ and…","keywords":["mesic nuclei","non-mesonic decay","optical potential","Monte Carlo simulations","eta-mesic helium-3","eta to two gamma","eta to three neutral pions","Klein-Gordon equation"],"falsifier":"Measure the $^{3}\\mathrm{He}$ recoil momentum spectrum in $pd\\to{}^{3}\\mathrm{He}2\\gamma$ events near the $\\eta$ threshold. If the non-mesonic decay mechanism operates, the spectrum should track the narrow $|R(p)|^{2}p^{2}$ distribution computed from the Klein-Gordon wavefunction, not the broader nucleon Fermi distribution; observing a broad nucleon-like spectrum, or no bound-state signal at the predicted binding energies of roughly 4--13~MeV, would falsify the model. Alternatively, extracting the in-medium branching ratio from the $2\\gamma$ to $6\\gamma$ yield ratio and finding it outside the predicted $2\\times10^{-5}$ to $7\\times10^{-4}$ range would rule out the assumed potential strengths.","tokens_in":6846,"feed_emoji":"⚛️","tokens_out":15580,"duration_ms":139594,"temperature":0.7,"pith_summary":"This paper presents a theoretical model for how a hypothetical $\\eta$-mesic $^{3}\\mathrm{He}$ nucleus could decay without the $\\eta$ being absorbed: the bound $\\eta$ decays while still moving inside the nucleus, into two photons or into three neutral pions that then decay to six photons. The authors solve the Klein-Gordon equation with an assumed optical potential between the $\\eta$ and the $^{3}\\mathrm{He}$ nucleus, obtaining the bound-state wavefunction, the relative $^{3}\\mathrm{He}$-$\\eta$ momentum distribution, and in-medium branching ratios for the two decay channels. These outputs replace the previous practice of assuming the decaying particle's momentum distribution is the same as the nucleon distribution in the nucleus. The model supplies the kinematic and efficiency template needed to interpret the search for $\\eta$-mesic $^{3}\\mathrm{He}$ in $pd\\to{}^{3}\\mathrm{He}2\\gamma$ and $pd\\to{}^{3}\\mathrm{He}6\\gamma$ reactions.","feed_headline":"Eta bound to helium-3: decay rates and momenta now predicted","feed_subtitle":"The model yields the momentum spread and in-medium branching ratios needed to interpret the helium-3 search data.","key_machinery":"The central object is the solution of the Klein-Gordon equation for the $\\eta$-$^{3}\\mathrm{He}$ bound state, with the optical potential $U_{\\mathrm{opt}}(r)=(V_{0}+iW_{0})\\rho(r)/\\rho_{0}$ encoding the strong interaction between the meson and the nucleus. The complex energy $E_{\\mathrm{KG}}$ gives the binding energy and the nuclear absorption width $\\Gamma_{\\mathrm{abs}}=-2\\,\\mathrm{Im}(E_{\\mathrm{KG}})$, while the Fourier transform of the coordinate-space wavefunction yields the relative momentum-space wavefunction and hence the momentum distribution $|R(p)|^{2}p^{2}$. In-medium branching ratios follow from $\\mathrm{BR}^{*}=\\Gamma_{\\eta\\to X}/(\\Gamma_{\\mathrm{tot}}^{\\eta}+\\Gamma_{\\mathrm{abs}})$, using vacuum partial widths. This machinery converts the assumed potential parameters into the experimentally usable kinematic distributions and decay probabilities.","core_discovery":"The central claim is that the non-mesonic decay of the $\\eta$-mesic $^{3}\\mathrm{He}$ bound state can be modeled by treating the $\\eta$ as a particle bound in an optical potential $U_{\\mathrm{opt}}(r)=(V_{0}+iW_{0})\\rho(r)/\\rho_{0}$, with its Fermi momentum obtained from the Fourier transform of the Klein-Gordon wavefunction. For four assumed parameter sets $(V_{0},W_{0})=(-75,-20),(-90,-20),(-75,-1),(-90,-1)$~MeV, the calculation yields binding energies between about 4 and 13~MeV, absorption widths between about 0.8 and 21~MeV, and in-medium branching ratios from about $2\\times10^{-5}$ to $7\\times10^{-4}$ for both $\\eta\\to2\\gamma$ and $\\eta\\to3\\pi^{0}$. The resulting momentum distribution is generally narrower than the nucleon Fermi distribution in $^{3}\\mathrm{He}$, so it changes how simulated events populate the detector. The paper also finds that pion final-state interactions reduce the $3\\pi^{0}$ signal by half, and that the geometrical acceptance of the detector is about 60% for the $2\\gamma$ channel and 40% for the $6\\gamma$ channel. These numbers are offered as the inputs that Monte Carlo simulations of the search reactions need to interpret the experimental data.","pith_inferences":["An immediate consequence the authors do not spell out: a null result in the $2\\gamma$/$6\\gamma$ search cannot by itself exclude an $\\eta$-mesic $^{3}\\mathrm{He}$ state, because the predicted branching ratios are so small; the same data would need to be reinterpreted with the production cross-section and acceptance folded in.","The same Klein-Gordon machinery could be applied to the $\\eta$-mesic $^{4}\\mathrm{He}$ non-mesonic decay channels, using the density and potential parameters already constrained in the deuteron-deuteron experiments, giving comparable in-medium branching ratios and momentum distributions testable with existing data.","Because the predicted $^{3}\\mathrm{He}$-$\\eta$ momentum distribution is much narrower than the nucleon Fermi distribution for the deeper potentials, the measured $^{3}\\mathrm{He}$ recoil momentum spectrum in the $2\\gamma$ channel could distinguish $V_{0}=-75$~MeV from $V_{0}=-90$~MeV even with modest statistics."],"forward_implications":["Monte Carlo simulations of the $pd\\to{}^{3}\\mathrm{He}2\\gamma$ and $pd\\to{}^{3}\\mathrm{He}6\\gamma$ searches can now use a bound-$\\eta$ momentum distribution derived from the wavefunction instead of the nucleon Fermi distribution, changing the expected $^{3}\\mathrm{He}$ recoil spectrum and the detector efficiency.","The in-medium branching ratios are suppressed by factors of hundreds to tens of thousands relative to vacuum because the nuclear absorption width dominates the total width, so the non-mesonic channels are rare and observed signals would require sizeable production cross sections.","Pion final-state interactions halve the $6\\gamma$ signal strength, so analyses of that channel must include the correction when estimating yields.","Measuring the binding energy and width from the excitation function, together with the observed $2\\gamma$/$6\\gamma$ yield ratio, would discriminate among the four assumed optical potential parameter sets."],"supporting_citations":[{"why":"Supplies the weakly absorptive $W_0=-1$~MeV potential parameter imported from the $\\eta$-$^{4}\\mathrm{He}$ analysis.","marker":"[7]"},{"why":"Provides the model with which measured excitation functions were compared to constrain the $\\eta$-nucleus optical potential parameters.","marker":"[14]"},{"why":"Documents the previous Monte Carlo assumption that the decaying particle's momentum distribution equals the nucleon distribution, which this model replaces.","marker":"[58]"},{"why":"Supplies the vacuum $\\eta$ total width and the $\\eta\\to2\\gamma$ and $\\eta\\to3\\pi^{0}$ branching ratios used in the in-medium branching ratio formula.","marker":"[61]"},{"why":"Gives the theoretical $^{3}\\mathrm{He}$ density distribution that enters the optical potential.","marker":"[62]"},{"why":"Provides the $^{3}\\mathrm{He}$ density distribution used in the calculation.","marker":"[63]"},{"why":"Describes the detector geometry used for the acceptance simulation.","marker":"[65]"}],"fun_headline_variants":["Eta in helium-3: decay widths and momentum now calculated","Predicting eta-mesic 3He non-mesonic decays","Model gives decay rates for eta bound in helium-3","Helium-3 eta decay: new branching ratios and momenta"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation depends on assuming the force between the eta and the helium-3 nucleus has one of four particular strengths, with the weakly absorbing value taken from a different nucleus; no data or uncertainty constrains that force, and every predicted number changes if the true force differs.","fun_headline_variants_meta":{"raw":{"variants":["Eta in helium-3: decay widths and momentum now calculated","Predicting eta-mesic 3He non-mesonic decays","Model gives decay rates for eta bound in helium-3","Helium-3 eta decay: new branching ratios and momenta"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000597,"raw_usage":{"total_tokens":2829,"prompt_tokens":1019,"completion_tokens":1810,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":635,"completion_tokens_details":{"reasoning_tokens":1736}},"tokens_in":635,"tokens_out":1810,"duration_ms":13480,"temperature":1.0,"reasoning_tokens":1736,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:13:05.240190+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $^{3}\\mathrm{He}$ recoil momentum spectrum in $pd\\to{}^{3}\\mathrm{He}2\\gamma$ events near the $\\eta$ threshold. If the non-mesonic decay mechanism operates, the spectrum should track the narrow $|R(p)|^{2}p^{2}$ distribution computed from the Klein-Gordon wavefunction, not the broader nucleon Fermi distribution; observing a broad nucleon-like spectrum, or no bound-state signal at the predicted binding energies of roughly 4--13~MeV, would falsify the model. Alternatively, extracting the in-medium branching ratio from the $2\\gamma$ to $6\\gamma$ yield ratio and finding it outside the predicted $2\\times10^{-5}$ to $7\\times10^{-4}$ range would rule out the assumed potential strengths.","supporting_citations":[{"cited_title":"Skurzok et al","cited_arxiv_id":null,"evidence_quote":"Supplies the weakly absorptive $W_0=-1$~MeV potential parameter imported from the $\\eta$-$^{4}\\mathrm{He}$ analysis."},{"cited_title":"Ikeno et al","cited_arxiv_id":null,"evidence_quote":"Provides the model with which measured excitation functions were compared to constrain the $\\eta$-nucleus optical potential parameters."},{"cited_title":"Nogga, Ph","cited_arxiv_id":null,"evidence_quote":"Documents the previous Monte Carlo assumption that the decaying particle's momentum distribution equals the nucleon distribution, which this model replaces."},{"cited_title":"Tanabashi, et al., Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the vacuum $\\eta$ total width and the $\\eta\\to2\\gamma$ and $\\eta\\to3\\pi^{0}$ branching ratios used in the in-medium branching ratio formula."},{"cited_title":"Hiyama, B","cited_arxiv_id":null,"evidence_quote":"Provides the $^{3}\\mathrm{He}$ density distribution used in the calculation."}],"review_version":1}