REVIEW 4 major objections 5 minor 35 references
A day of muon data could probe lepton-flavor-violating scalar couplings to 10^-5.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 14:42 UTC pith:5OPE6FH2
load-bearing objection 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. the 4 major comments →
Production of lepton-flavor-violating scalars through resonant positive-muon annihilation on atomic electrons
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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
What carries the argument
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
Load-bearing premise
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.
What would settle it
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).
If this is right
- 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.
Where Pith is reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [§2.4, Eqs. (2.14) and (2.17)] 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.
- [§2.4, Eqs. (2.13)–(2.16)] 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.
- [§2.1 and §2.4, f,g functions] 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.
- [§2, invisible branching ratio] 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.
minor comments (5)
- [Eq. (2.9)–(2.10)] 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).
- [Eq. (2.6) and Eq. (2.14)] 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).
- [Fig. 4 and text] 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.
- [Table 1] 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.
- [Ref. [17]] 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.
Circularity Check
No circular reduction found; minor self-citation/reproducibility debt only.
full rationale
The derivation chain is self-contained in the relevant sense: Eq. (2.8) defines the amplitude from the effective LFV interaction, Eqs. (2.9)-(2.13) square it using Dirac spinor completeness with the (2j+1) degeneracy, Eq. (2.14) performs the phase-space integration, and Eq. (2.17) sums over occupied subshells with RHF occupation numbers. The projected limits in Sec. 4 come from a profile likelihood with Asimov data, not by fitting g_phi to the reach curve. No output equation is equivalent to an input by construction. The only self-referential element is the numerical f,g input: Sec. 2 states 'We further use relativistic atomic Dirac spinors leveraging the techniques from [22] as will be discussed in [17]', and the acknowledgments state Prof. Darmé supplied the f,g values. Since [17] is a companion paper with overlapping authorship containing the concrete numerical wave functions, this is an external-reproducibility debt, but it is an omitted/forthcoming input, not a circular reduction of the cross-section claim. The Br(phi -> chi chi-bar) ~ 1 assumption is a stated physical input, not a fitted output. The alleged (2j+1) double-counting between Eqs. (2.11), (2.14), and (2.17) is a possible normalization/correctness issue affecting the absolute reach, but it does not make the prediction equal to its input; it would change the limit by a constant factor. The comparison against the Compton-profile and free/rest treatments provides an independent cross-check. Thus no circular step is exhibited; score 2 reflects the minor non-load-bearing self-citation and reproducibility debt.
Axiom & Free-Parameter Ledger
free parameters (6)
- g_D (dark-sector coupling) =
1e-4 (benchmark)
- Br(phi->chi chi-bar) =
~1 (assumed)
- RPC single-plane efficiency epsilon_RPC,0 =
95%
- RPC efficiency uncertainty delta_epsilon_RPC =
0.01% and 1%
- Scintillator efficiency epsilon_S,0 =
99%
- HIAF beam parameters =
7.5e5 mu+/s, sigma_x=sigma_y=30 mm, delta E/E=2%
axioms (6)
- ad hoc to paper The only SM interaction of phi is the effective LFV coupling L_int = -g_phi phi ebar mu (plus the dark coupling to chi).
- domain assumption The bound electron is represented by momentum-space Dirac spinors U_E kappa m built from RHF radial functions f,g; the muon is a free external scattering state.
- domain assumption All occupied atomic electrons contribute incoherently and the total cross section is a sum over subshells with occupation numbers N_nk summing to Z=82.
- ad hoc to paper phi decays invisibly with Br(phi->chi chi-bar) approximately 1.
- domain assumption Geant4 with the FTFP_BERT physics list plus 89 supplementary background processes accurately predicts the relevant backgrounds.
- standard math The profile-likelihood ratio with background-only Asimov data and unit Gaussian priors yields valid 90% CL limits.
invented entities (2)
-
phi (e-mu LFV scalar)
no independent evidence
-
chi (dark-sector particle)
no independent evidence
read the original abstract
We investigate an invisible lepton-flavor-violating scalar $\phi$ with exclusive $e-\mu$ couplings and study its resonant production via $\mu^+e^-\to\phi$ in fixed-target experiments. Since the effective center-of-mass energy is determined by the momentum of the initial-state bound electrons, atomic effects can significantly affect the resonance behavior. We therefore employ relativistic bound-state electron wave functions to calculate the production cross section and reveal a material-dependent broadening of the resonance lineshape. For the proposed HIAF experiment, fewer than one day of data taking ($6\times10^{10}$ MOT) can probe couplings at the $10^{-5}$ level at 90\% confidence level near resonance, demonstrating that high-intensity muon fixed-target experiments provide a powerful complementary probe of lepton-flavor violation.
Reference graph
Works this paper leans on
-
[1]
S. N. Gninenkoet al.[NA64], Leptonic scalar portal: Origin of muong−2 anomaly and dark matter?, Phys. Rev. D106, no.1, 015003 (2022) doi:10.1103/PhysRevD.106.015003 [arXiv:2202.04410 [hep-ph]]
Pith/arXiv arXiv 2022
-
[2]
Chenet al., DREAMuS: Dark matter REsearch with Advanced Muon Source, [arXiv:2604.10257 [hep-ph]]
X. Chenet al., DREAMuS: Dark matter REsearch with Advanced Muon Source, [arXiv:2604.10257 [hep-ph]]
-
[3]
A.M. Baldini, Y. Bao, E. Baracchini,et al., Search for the lepton flavour violating decay µ+ →e +γwith the full dataset of the MEG experiment, Eur. Phys. J. C76, 434 (2016), doi:10.1140/epjc/s10052-016-4271-x
-
[4]
MEG II Collaboration, New limit on theµ + →e +γdecay with the MEG II experiment, Eur. Phys. J. C85, 1177 (2025), doi:10.1140/epjc/s10052-025-14906-3
-
[5]
W. Bertl, S. Egli, R. Eichler,et al., Search for the decayµ + →e +e+e−, Nucl. Phys. B260, 1–31 (1985), doi:10.1016/0550-3213(85)90308-6
-
[6]
W. Bertl, R. Engfer, E.A. Hermes,et al., A search forµ-e conversion in muonic gold, Eur. Phys. J. C47, 337–346 (2006), doi:10.1140/epjc/s2006-02582-x
-
[7]
L. Willmann, P.V. Schmidt, H.P. Wirtzet al., New Bounds from a Search for Muonium to Antimuonium Conversion, Phys. Rev. Lett.82, 49–52 (1999), doi:10.1103/PhysRevLett.82.49. – 14 –
-
[8]
A. Y. Bai, H. J. Cai, C. L. Chenet al., Conceptual design of the muonium-to-antimuonium conversion experiment (MACE), Nucl. Sci. Tech.37, 57 (2026). doi:10.1007/s41365-025-01876-0
-
[9]
L. Gao, C. e. Liu, Q. Li, C. Zhou, Q. Li, L. Chen, X. Zhang, Y. Xu and Z. Sun, Probing and knocking with muons, Mod. Phys. Lett. A40(2025) no.24, 2530008 doi:10.1142/S0217732325300083 [arXiv:2503.22956 [hep-ph]]
Pith/arXiv arXiv 2025
-
[10]
Probing and knocking with muons and new physics exploration,
C. e. Liu, L. Gao, Z. Wang, J. Li, M. Fan, Z. Qin, R. Zhang, Y. Xu, X. Zhang and Q. Li,et al.“Probing and knocking with muons and new physics exploration,” Chin. Sci. Bull.71 (2025) no.4, 894-903 doi:10.1360/CSB-2025-5452
-
[11]
X. Yu, Z. Wang, C. e. Liu, Y. Feng, J. Li, X. Geng, Y. Zhang, L. Gao, R. Jiang and Y. Wu, et al.Proposed Peking University muon experiment for muon tomography and dark matter search, Phys. Rev. D110(2024) no.1, 016017 doi:10.1103/PhysRevD.110.016017 [arXiv:2402.13483 [hep-ex]]
Pith/arXiv arXiv 2024
-
[12]
C. e. Liu, R. Zhang, Z. Wang, A. M. Levin, L. Gao, J. Li, M. Fan, Y. Wu, Z. Qin, Y. Banet al., Probing Cosmic Ray Composition and Muonphilic Dark Matter via Muon Tomography, Phys. Rev. Lett.136, no.15, 151001 (2026) doi:10.1103/5jh7-fxf4
-
[13]
F. An, D. Bai, H. Cai, S. Chen, X. Chen, H. Duyang, L. Gao, S. F. Ge, S. Ge and J. He,et al.High-Precision Physics Experiments at Huizhou Large-Scale Scientific Facilities, Chin. Phys. Lett.42(2025) no.11, 110102 doi:10.1088/0256-307X/42/11/110102 [arXiv:2504.21050 [hep-ph]]
arXiv 2025
-
[14]
F. Arias-Arag´ on, L. Darm´ e, G. G. di Cortona and E. Nardi, Production of Dark Sector Particles via Resonant Positron Annihilation on Atomic Electrons, Phys. Rev. Lett.132, no.26, 261801 (2024) [erratum: Phys. Rev. Lett.133, no.21, 219901 (2024)] doi:10.1103/PhysRevLett.132.261801 [arXiv:2403.15387 [hep-ph]]
Pith/arXiv arXiv 2024
-
[15]
F. Arias-Arag´ on, G. G. di Cortona and E. Nardi, Unveiling atomic electron motion effects in e+e− high energy collisions at NA64, JHEP10, 174 (2025) doi:10.1007/JHEP10(2025)174 [arXiv:2507.08941 [hep-ph]]
arXiv 2025
-
[16]
F. Arias-Arag´ on, L. Darm´ e, G. G. di Cortona and E. Nardi, Atoms as Electron Accelerators for Measuring the Cross Section ofe +e− →Hadrons, Phys. Rev. Lett.134, 061802 (2025) doi:10.1103/PhysRevLett.134.061802
-
[17]
Darm´ e, Q
L. Darm´ e, Q. Li, H.-S. Shao, J. Shen, Y. Wu, To appear (2026)
2026
-
[18]
R. Plestid and M. ˜B. Wise, Atomic binding corrections for high-energy fixed target experiments, Phys. Rev. D110, 056032 (2024) doi:10.1103/PhysRevD.110.056032
-
[19]
F. Ignatov, R. N. Pilato, T. Teubner and G. Venanzoni, An alternative evaluation of the leading-order hadronic contribution to the muong−2 with MUonE, Phys. Lett. B848 (2024), 138344 doi:10.1016/j.physletb.2023.138344 [arXiv:2309.14205 [hep-ph]]
arXiv 2024
-
[20]
Abbiendi, Status of the MUonE experiment, Phys
G. Abbiendi, Status of the MUonE experiment, Phys. Scr.97, 054007 (2022), doi:10.1088/1402-4896/ac6297
-
[21]
M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, 1st ed., Addison-Wesley (1995), isbn:978-0-201-50397-5
1995
-
[22]
S. Gardner and J. Piekarewicz, Relativistically generated asymmetry in the missing-momentum distribution from the (e,e’p) reaction, Phys. Rev. C50, 2822 (1994) doi:10.1103/PhysRevC.50.2822. – 15 –
-
[23]
F. Biggs, L. B. Mendelsohn, and J. B. Mann, Hartree-Fock Compton profiles for the elements, At. Data Nucl. Data Tables16, 201–309 (1975) doi:10.1016/0092-640X(75)90030-3
-
[24]
Huzinaga, Gaussian-Type functions for polyatomic systems
S. Huzinaga, Gaussian-Type functions for polyatomic systems. I, J. Chem. Phys.42, 1293–1302 (1965) doi:10.1063/1.1696113
-
[25]
E.O. Steinborn and K. Ruedenberg, Rotation and Translation of Regular and Irregular Solid Spherical Harmonics, Adv. Quantum Chem.7, 1–81 (1973) doi:10.1016/S0065-3276(08)60558-4
-
[26]
E. Clementi and C. Roetti, Roothaan-Hartree-Fock atomic wavefunctions: Basis functions and their coefficients for ground and certain excited states of neutral and ionized atoms, Z≤ 54, At. Data Nucl. Data Tables14, 177–478 (1974) doi:10.1016/S0092-640X(74)80016-1
-
[27]
P. Kaijser and V. H. Smith, Jr., Evaluation of Momentum Distributions and Compton Profiles for Atomic and Molecular Systems, Adv. Quantum Chem.10, 37–76 (1977) doi:10.1016/S0065-3276(08)60578-X
-
[28]
E. Filter and E. O. Steinborn, Extremely compact formulas for molecular two-center one-electron integrals and Coulomb integrals over Slater-type atomic orbitals, Phys. Rev. A 18, 1–11 (1978) doi:10.1103/PhysRevA.18.1
-
[29]
D. K. Maretis, Talmi transformation and the multicenter integrals of harmonic oscillator functions, J. Chem. Phys.71, 917–921 (1979) doi:10.1063/1.438381
-
[30]
E. J. Weniger and E. O. Steinborn, The Fourier transforms of some exponential-type basis functions and their relevance to multicenter problems, J. Chem. Phys.78, 6121–6132 (1983) doi:10.1063/1.444574
-
[31]
E. J. Weniger, Weakly convergent expansions of a plane wave and their use in Fourier integrals, J. Math. Phys.26, 276–291 (1985) doi:10.1063/1.526970
-
[32]
X. Zhou and J. Yang, Status of the High Intensity Heavy-Ion Accelerator Facility in China, AAPPS Bull.32, 35 (2022) doi:10.1007/s43673-022-00064-1
-
[33]
Y. Xuet al., Feasibility study of a GeV-energy muon source based on the High Intensity Heavy-Ion Accelerator Facility, Phys. Rev. Accel. Beams28, 053401 (2025) doi:10.1103/PhysRevAccelBeams.28.053401
-
[34]
C. Cesarottiet al., Sensitivity potential to a light flavor-changing scalar boson with DUNE and NA64µ, Eur. Phys. J. C83, 775 (2023) doi:10.1140/epjc/s10052-023-11891-3 [arXiv:2306.07405 [hep-ex]]
Pith/arXiv arXiv 2023
-
[35]
Abiet al.[Muong-2], Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm, Phys
B. Abiet al.[Muong-2], Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm, Phys. Rev. Lett.126, 141801 (2021) doi:10.1103/PhysRevLett.126.141801 [arXiv:2104.03281 [hep-ex]]. – 16 –
arXiv 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.