{"id":"ebb8400f-565d-44bf-9678-b083600a1f41","arxiv_id":"2607.18669","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Relativistic bound-electron momentum distributions broaden the mu+ e- -> phi resonance, and one day of the proposed HIAF-PKMu exposure could exclude g_phi down to ~10^-5 in the 100-200 MeV window.","lead":"High-intensity muons hitting a lead target could produce a hypothetical particle that couples electrons to muons; the authors calculate how atomic electron motion broadens the resonance and find one day of HIAF running could probe couplings down to 10^-5. The result is a projected sensitivity, not a measurement, and rests on an invisible-decay assumption and beam parameters that have not yet been demonstrated.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Normalization of Eq. (2.14) combined with occupation numbers in Eq. (2.17) likely double-counts the (2j+1) angular degeneracy, inflating the projected HIAF reach.","rationale":"The reader's verdict is CONDITIONAL, and I agree that the paper should not be accepted without revision. My most load-bearing concern is the normalization/double-counting issue in Eqs. (2.14) and (2.17), which the reader also flagged in their rationale as a possible double-counting concern. However, the reader's stated 'weakest_assumption' emphasizes the invisible-decay branching ratio and HIAF beam parameters; those are model/beam assumptions that the paper explicitly declares, whereas the normalization issue is an internal consistency problem that directly controls all numerical sensitivity projections. If the double counting is real, the absolute cross section is inflated and the projected limits in Figs. 4 and 5 shift. This does not necessarily destroy the qualitative claim—the corrected reach may still be in the 10^-5 ballpark—but it must be fixed and the projections redone before the quantitative headline can be trusted. The unpublished f,g wavefunctions from Ref. [17] are an additional reproducibility concern, but they are secondary to the normalization ambiguity. I therefore recommend keeping the CONDITIONAL verdict (UNCHANGED) while requiring the authors to clarify the normalization and rerun the projections.","tokens_in":12049,"tokens_out":13255,"duration_ms":118515,"concrete_test":"Recompute sigma_Pb using Eq. (2.14) with the replacement sum_{nkappa} N_nkappa sigma_nkappa -> sum_{nkappa} sigma_nkappa (or, equivalently, re-derive Eq. (2.9) with an average 1/(2j+1) over m and keep N_nkappa), then regenerate the one-day exclusion curve in Fig. 5. If the g_phi limit changes by more than ~20% (or the Pb peak cross section in Fig. 2 drops by the occupation-weighted factor), the double counting is confirmed and the headline reach must be revised. An independent cross-check: for a hydrogenic 1s shell (j=1/2, N=2), compare the resulting sigma with the Compton-profile convolution normalized to one electron; the current formula yields a factor-2 excess.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central reach claim depends on the absolute normalization of the Pb cross section. Eq. (2.9) defines |M|^2_kappa with a sum over the magnetic quantum number m, and Eq. (2.11) supplies the factor (2j+1). Eq. (2.14) therefore already contains the full angular degeneracy for the subshell. Eq. (2.17) then multiplies each sigma_nkappa by the occupation number N_nkappa; for a closed subshell N_nkappa = 2j+1. The two steps count the same m-states twice, so the Pb cross section is enhanced by an extra factor of (2j+1) per subshell (e.g. 8 for the 4f subshell). The paper's consistency claim—that the momentum density is normalized to Z=82—does not resolve this: the stated spinor normalization in Eq. (2.5) is (1/(2pi)^3) sum_m ∫ U†U d^3k = 2E_A, not a number-density normalization to Z. If sigma_nkappa is intended per m-state, then Eq. (2.9) should average over m rather than sum; if it is intended per subshell, Eq. (2.17) should not include N_nkappa. Either way, the product (2j+1) × N_nkappa overcounts by a factor of 2j+1. Since the projected coupling limit scales roughly as g_phi ∝ sigma_Pb^{-1/2}, the headline '10^-5 in less than one day' could shift by up to O(sqrt(8)), and the relative subshell contributions in Fig. 2 change.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a hypothetical lepton-flavor-violating scalar φ with an off-diagonal e–μ Yukawa coupling g_φ and an additional dark-sector coupling g_D. It computes the resonant production cross section μ+ e− → φ on bound atomic electrons in a lead target, using relativistic Dirac bound-state wave functions (with RHF f,g functions supplied by a co-author and to appear in Ref. [17]). The authors find a material-dependent resonance broadening that differs from both the electron-at-rest approximation and the Compton-profile folding method. They then simulate the proposed HIAF-PKMu detector with Geant4 and derive projected 90% CL limits on g_φ, claiming that less than one day of running (6×10^10 MOT) can probe couplings at the 10^−5 level near resonance, assuming the φ decays invisibly, Br(φ→χχ̄)≈1.","tokens_in":12502,"tokens_out":8760,"duration_ms":87996,"significance":"If the normalization is correct, the paper offers a genuinely more rigorous treatment of bound-electron effects for this process than the Compton-profile convolution, and it gives a concrete, falsifiable experimental projection for a facility that is currently under construction. The strengths are the use of Dirac bound-state amplitudes, the attempt to sum all occupied subshells of lead, a detailed Geant4 simulation, and a profile-likelihood treatment of detector-efficiency systematics. However, the absolute normalization of the central cross section is not transparent, and there is a likely double-counting of the (2j+1) magnetic degeneracy. Because the headline coupling reach scales as g_φ ∝ σ_Pb^{−1/2}, this issue directly affects the paper's central claim and must be resolved before the projections can be trusted.","major_comments":[{"comment":"The absolute normalization appears to double-count the magnetic degeneracy. Eq. (2.9) defines |M|^2_κ with an explicit sum over m, and Eq. (2.11) evaluates that sum as proportional to (2j+1). Thus σ_nκ in Eq. (2.14) is already the full-subshell cross section. Eq. (2.17) then multiplies by N_nκ = 2j+1 for a closed subshell, counting the same m-states twice. The normalization in Eq. (2.5), sum_m ∫ U†U d^3k/(2π)^3 = 2E_A, is a spinor normalization, not a number-density normalization to Z; it does not repair the overcount. Either average over m in Eq. (2.9) and keep N_nκ, or drop N_nκ in Eq. (2.17). Since g_φ ∝ σ^{-1/2}, the projected limits in Figs. 4 and 5 are shifted by a factor up to O(√(2j+1)), and the subshell weighting in Fig. 2 also changes.","section":"§2.4, Eqs. (2.14) and (2.17)"},{"comment":"The jump from the squared amplitude to the integrated cross section is not shown. Eq. (2.14), including the denominator 8 p_B^2 E_A k_A, the x0 substitution, and the kinematic limits in Eq. (2.16), is asserted without derivation. The absolute normalization therefore cannot be independently checked from the text. Please provide the full phase-space integration, including the treatment of the energy δ-function and the flux factor for a bound electron, and demonstrate that the free/rest and Compton-profile limits are recovered. This is load-bearing because the projected reach depends directly on the absolute value of σ_Pb.","section":"§2.4, Eqs. (2.13)–(2.16)"},{"comment":"The numerical f and g functions are not given in the manuscript; the acknowledgments state they were supplied by Prof. Luc Darmé and will appear in companion paper Ref. [17]. Since Fig. 2 and all sensitivity projections depend on these functions, the calculation is not reproducible from the present text. Please include the RHF radial functions or a numerical parametrization, or make the values publicly available, so that Eq. (2.14) can be evaluated independently.","section":"§2.1 and §2.4, f,g functions"},{"comment":"The headline sensitivity assumes Br(φ→χχ̄)≃1, which requires g_D > g_φ and the existence of a dark state χ. If the scalar instead decays visibly to e± μ∓, the downstream veto removes signal and the quoted g_φ limits do not apply. The abstract's claim of probing 'couplings at the 10^-5 level' is therefore contingent on a dark-sector parameter. Please state this assumption explicitly in the abstract and summary, and either present the reach as a function of (g_φ, g_D, Br) or specify the benchmark g_D used. As written, the world-leading probe claim is broader than what the analysis actually constrains.","section":"§2, invisible branching ratio"}],"minor_comments":[{"comment":"The projector is miswritten: it should be sum_s v(P_B,s) ar v(P_B,s) = /P_B − m_B, not sum_s ar v(P_B,s) v(P_B,s). As printed, the first factor is a scalar and cannot be multiplied into the trace in Eq. (2.12).","section":"Eq. (2.9)–(2.10)"},{"comment":"The notation 'X Z' in Eq. (2.6) is confusing; Z is elsewhere defined as atomic number. It should be written as a sum over occupied subshells or over electrons. There is also a stray bracket/typographical imbalance in the numerator of Eq. (2.14).","section":"Eq. (2.6) and Eq. (2.14)"},{"comment":"The horizontal axis is labeled 'm (MeV)' but the figures and text use m_φ; the plotted range 105–140 MeV is narrower than the quoted 100–200 MeV window. Please make the mass label and range consistent.","section":"Fig. 4 and text"},{"comment":"The header 'Signal(×10−4)' is unclear, and the table does not state the reference coupling at which the signal sample is generated or how the reweighting to arbitrary g_φ is implemented. Please clarify.","section":"Table 1"},{"comment":"Since Ref. [17] supplies the central numerical f,g input, the placeholder 'To appear (2026)' is insufficient. Provide a fuller reference or, preferably, make the data available in this paper.","section":"Ref. [17]"}],"recommendation":"major_revision","confidential_remarks":"The (2j+1) double-counting appears likely and is directly testable by a benchmark derivation. If confirmed, the projected limits weaken by up to about √8, but the framework remains publishable after correction. I would also press the authors to make the f,g data available, because the current reliance on an unpublished companion paper prevents independent verification of the central numerical result. The 'world-leading probe' wording should be tempered unless the invisible-decay assumption is elevated to the abstract. The paper is within scope for JHEP and the experimental projection is interesting, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is a serious, well-motivated proposal for probing e-mu LFV scalars via mu+ e- -> phi at HIAF, and the atomic-binding treatment is genuinely new. But I think the absolute normalization of the lead cross section is wrong: the equations as written double-count the (2j+1) degeneracy, so the projected reach is over-optimistic by up to about a factor of three in coupling.\n\nWhat's good: applying relativistic bound-state Dirac wave functions at the amplitude level to this channel, and showing the material-dependent resonance broadening. The comparison with the Compton-profile approach is useful, and the detector simulation is unusually concrete—background processes, cutflow, and a data-driven RPC calibration scheme are all there. The authors are also honest about the invisible-decay assumption.\n\nThe main problem: Eq (2.9) defines |M|^2_kappa with a sum over m, and Eq (2.11) brings in the (2j+1) factor. So Eq (2.14) already contains the full subshell degeneracy. Then Eq (2.17) multiplies by the occupation number N_nkappa, which is 2j+1 for a closed subshell. Same m-states counted twice. The claim that the momentum density is normalized to Z=82 doesn't fix this, because Eq (2.5) is not a number-density normalization. As written, sigma_Pb is too large by a factor of 2j+1 per subshell; for the 4f shell that's 8. Since g_phi scales as sigma^-1/2, the headline '10^-5 in less than a day' shifts to roughly 3e-5, and the subshell weights in Fig 2 change. This is not a peripheral detail—the projected limit is the paper's main result.\n\nTwo smaller issues: the f,g functions come from an unpublished companion paper, so the calculation can't be independently checked from the deposited text; and the reach assumes Br(phi->chi chi-bar)=1 with fixed HIAF beam parameters, which should be treated as a systematic band rather than a fixed input.\n\nBottom line: worth refereeing, but a referee should be asked to verify the normalization in Section 2.4 carefully. I'd want the companion paper (or data/code) released before trusting the numbers.","headline":"A well-motivated LFV scalar search at HIAF, but the central cross-section normalization looks double-counted, so the headline reach is probably too optimistic by a factor of a few.","tokens_in":13047,"tokens_out":5511,"would_cite":false,"duration_ms":86809,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V10","81T80"],"pacs":["13.10.+q","14.80.-j"],"model":"deepseek-v4-flash","headline":"A day of muon data could probe lepton-flavor-violating scalar couplings to 10^-5.","keywords":["lepton flavor violation","muon fixed-target experiment","invisible scalar","resonant production","bound-state Dirac wave function","HIAF","mu+e- annihilation","dark sector"],"falsifier":"A direct calculation of the signal yield for the visible decay mode phi -> e+ mu- using the same detector geometry and veto logic: if the projected limits for the visible mode are significantly weaker than the invisible-mode limits, the headline claim of 10^-5 reach would not hold for models where the scalar has a large visible branching ratio. Alternatively, a precise measurement of the mu+ e- -> phi cross section on a thin lead target at a known beam energy would test the absolute normalization of Eq. (2.14).","tokens_in":11886,"feed_emoji":"🧲","tokens_out":1269,"duration_ms":17701,"temperature":0.7,"pith_summary":"This paper proposes a new way to search for lepton-flavor violation: a positively charged muon annihilates with an atomic electron to produce a new invisible scalar particle. The authors calculate the production rate using relativistic bound-electron wave functions, showing that atomic motion broadens the resonance line shape in a material-dependent way. For a proposed fixed-target experiment at HIAF, they find that less than one day of data could constrain the muon-electron coupling g_phi to the 10^-5 level near resonance, making this the most sensitive probe in the 100–200 MeV mass range. This would provide a direct, complementary test of physics beyond the Standard Model that is not accessible through other channels.","feed_headline":"A day of muon data could probe LFV scalar couplings to 10^-5","feed_subtitle":"Resonant mu+ e- annihilation on bound lead electrons turns HIAF into the world's most sensitive probe of e-mu flavor violation.","key_machinery":"The central machinery is the relativistic bound-state Dirac spinor treatment of target electrons. The bound electron is described by Dirac wave functions with large and small radial components (f and g) and spinor spherical harmonics, and the squared Feynman amplitude is computed by summing over bound-electron magnetic quantum numbers. This yields a first-principles cross-section formula that naturally includes the electron momentum distribution, replacing the electron-at-rest approximation. The authors compare this with the Compton-profile convolution method and show that the new treatment predicts a slightly lower peak cross section because the bound-electron momentum also modifies the amp","core_discovery":"The paper demonstrates that resonant muon annihilation on bound atomic electrons, mu+ e- -> phi, can be a powerful probe of lepton-flavor-violating scalar interactions. Using relativistic Dirac bound-state wave functions for the target electrons, the authors compute the production cross section and show that atomic-motion effects broaden the resonance line shape in a target-dependent way, reducing the peak cross section relative to the electron-at-rest approximation. They then combine this cross section with detailed Geant4 simulations of the proposed HIAF-PKMu experiment, including a four-cut selection that suppresses backgrounds by more than 99.999%, and derive projected 90% confidence-lev","pith_inferences":["The paper implicitly assumes that the invisible decay channel dominates; if the scalar visibly decays to e+ mu-, the downstream-veto topology changes and the quoted g_phi limits would need to be re-derived. A testable extension would be to compute limits for the visible decay mode using the same detector simulation.","The sensitivity claim relies on the absolute normalization of the cross-section formula in Eq. (2.14), which is not fully derived in the text. A dedicated derivation or independent numerical check could verify this normalization, as it directly sets the signal yield.","The use of RHF wave functions for all 82 electrons of lead is a strong approximation; a more complete treatment with exact Dirac-Fock wave functions or including electron-electron correlation might shift the resonance shape. This could be tested by comparing with future precision measurements if a signal were found.","The paper does not discuss the potential interference between the resonant mu+ e- -> phi production and the Standard Model background (e.g., mu+ e- -> anything). If such interference is non-negligible, it could alter the signal yield even away from the resonance peak; this would be a natural follow-up calculation."],"forward_implications":["If correct, the HIAF experiment would become the world-leading probe of electron-muon lepton-flavor-violating scalars in the 100–200 MeV mass window, a region with no current competitive limits.","With only 10 minutes of data, the projected limit already surpasses the full-dataset result of NA64 mu, and one hour rivals the 2000-hour DREAMuS projection.","The material-dependent resonance broadening implies that target choice matters: experiments can tune the target material to optimize sensitivity for a given beam energy and scalar mass.","The relativistic bound-state formalism can be applied to other processes involving bound electrons, such as precision muon-electron scattering at MUonE, potentially improving theoretical predictions there.","If no signal is seen, the experiment would place strong constraints on models that explain the muon g-2 anomaly through a leptophilic scalar, closing a favored parameter space region."],"fun_headline_variants":["Muon annihilation on bound electrons exposes lepton-flavor violation","One day of muon data probes LFV scalars to 10^-5","Resonant muon-electron annihilation probes flavor-violating scalars","Atomic motion broadens muon resonance, aiding LFV sensitivity","HIAF muon beam finds flavor violation in a single day"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The projected sensitivity assumes the produced scalar decays invisibly (Br(phi->chi chi-bar) = 1), which requires a dark-sector particle chi with a coupling larger than g_phi; if the scalar instead decays visibly to e+ mu-, the quoted limits do not apply.","fun_headline_variants_meta":{"raw":{"variants":["Muon annihilation on bound electrons exposes lepton-flavor violation","One day of muon data probes LFV scalars to 10^-5","Resonant muon-electron annihilation probes flavor-violating scalars","Atomic motion broadens muon resonance, aiding LFV sensitivity","HIAF muon beam finds flavor violation in a single day"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001065,"raw_usage":{"total_tokens":4267,"prompt_tokens":675,"completion_tokens":3592,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":3514}},"tokens_in":419,"tokens_out":3592,"duration_ms":26686,"temperature":1.0,"reasoning_tokens":3514,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T14:42:50.715756+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct calculation of the signal yield for the visible decay mode phi -> e+ mu- using the same detector geometry and veto logic: if the projected limits for the visible mode are significantly weaker than the invisible-mode limits, the headline claim of 10^-5 reach would not hold for models where the scalar has a large visible branching ratio. Alternatively, a precise measurement of the mu+ e- -> phi cross section on a thin lead target at a known beam energy would test the absolute normalization of Eq. (2.14).","supporting_citations":[],"review_version":1}