REVIEW 2 major objections 4 minor 20 references
Momentum distribution of the electron pair from the charged lepton flavor violating process $\mu^-e^-\to e^-e^-$ in muonic atoms with a polarized muon
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read A polarized muon in a muonic atom yields a parity-violating electron asymmetry whose sign encodes the chirality of the lepton-flavor-violating interaction, and whose energy-angle pattern distinguishes photonic from contact operators.
desk verdict Solid, useful extension of their unpolarized rate work to spin-dependent asymmetries; the chirality-sign claim needs the single-operator caveat made explicit. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The organizing object is the spin-correlation decomposition of the differential decay rate: $$\frac{d\Gamma}{d\epsilon_1 d\Omega_1 d\Omega_2} = \frac{1}{8\$pi^{2}$}\, \frac{d\Gamma_{\rm unpol.}}{d\epsilon_1 dc_{12}}\left[1 + F_S\,P\cdot \hat{p}_{12} + F_A\,P\cdot \hat{q}_{12} + F_D\,P\cdot (\hat{p}_1\times \hat{p}_2)\right].$$ The coefficients are built from multipole-expanded Dirac amplitudes for bound and scattered leptons in the Coulomb potential of a finite nuclear charge distribution; in the plane-wave limit the asymmetry functions reduce to radial integrals ($I_{gg}$, $I_{fg}$, $I_{gf}$, $I_{ff}$) and are proportional to the handedness parameter $h_a$ of the operator. That proportionality is what carries the argument: the sign of the asymmetry directly reflects the chirality of the muon coupling.
What would settle it
A polarized-muon muonic-atom experiment (for example on a heavy nucleus such as $^{208}$Pb) that measures $F_S$, $F_A$, and $F_D$ as functions of $\epsilon_1$ and $c_{12}$ would settle the claim: if the sign of the parity-violating asymmetry for the assumed dominant operator is opposite to the template, or if $F_D$ falls far outside the $10^{-2}$ to $10^{-1}$ range for real couplings, the chiral-sign relation and the final-state-interaction estimate would be refuted. A calculation that includes 2s or 2p bound electrons and checks whether the $g_1$-versus-$g_2$ sign relation survives would also directly test the paper's core prediction.
Extended reading notes
Core claim
The paper establishes, using Dirac wave functions in the Coulomb field of a finite nuclear charge distribution, that the differential decay rate for $\mu^- e^- \to e^- e^-$ with a polarized muon factorizes into the unpolarized rate times three polarization terms. The coefficients $F_S$ and $F_A$ are parity-violating asymmetries; their sign flips when the chirality of the muon operator flips, and their pattern in $\epsilon_1$ and $c_{12}$ differs between contact and photonic operators, with the photonic case giving the largest asymmetries. The third coefficient $F_D$ is parity-even but motion-reversal-odd; although the couplings are taken real, final-state Coulomb distortion generates $F_D$ at the level $10^{-2}$ to $10^{-1}$, largest for the photonic interaction. The calculation uses a 1s muon and a 1s electron in the initial state, and treats the emitted electrons as distorted waves including the finite nuclear size.
Load-bearing premise
The calculation assumes that only the 1s bound muon and 1s bound electron contribute to the initial state and that electron-electron rescattering is negligible, so if higher-n bound electrons or electron-electron rescattering contribute significantly, the computed asymmetry curves would shift.
Editorial extensions
If this is right
- The measured sign of $F_S$ and $F_A$ would assign a left- or right-handed muon coupling to the dominant CLFV operator, information that the decay rate alone cannot provide.
- The $\epsilon_1$ and $c_{12}$ dependence separates the photonic dipole case from contact scalar and vector cases, and distinguishes same-chirality contact operators ($g_1$/$g_3$) from opposite-chirality ones ($g_5$).
- A finite nuclear charge distribution is not a minor correction: for $g_1$–$g_4$ type interactions the asymmetry would essentially vanish for a point-like nucleus, so realistic nuclear wave functions are required for any asymmetry measurement.
- The motion-reversal-odd coefficient $F_D$ is predicted to be nonzero at $10^{-2}$ to $10^{-1}$ even for real couplings, with the largest value for photonic interactions, so any CP-violation search using $F_D$ must first subtract this final-state-interaction background.
- Combining the asymmetry observables with the atomic-number dependence of the decay rate from the authors' previous work gives a two-pronged strategy for identifying the new-physics operator behind charged lepton flavor violation.
- When the next-generation muon experiments search for $\mu^- e^- \to e^- e^-$, the polarization asymmetry provides a concrete, background-labeled observable to include in the experimental design.
Reading between the lines
- Inference: The same formalism could be used to optimize the choice of target nucleus, since lighter nuclei would suppress the final-state-interaction-induced $F_D$ background while heavier nuclei amplify the finite-size effect that makes the parity-violating asymmetry visible.
- Inference: The chiral-sign relation may extend to other bound-state CLFV transitions initiated by polarized muons, such as muon-to-electron conversion in atoms, but the authors do not perform that calculation here; testing it would require an analogous distorted-wave analysis.
- Inference: The paper's templates assume a single dominant operator; a practical experimental analysis could fit mixtures of operators and treat 2s or 2p bound-electron contributions as a systematic uncertainty, which the paper mentions as a possible extension of its multipole formulas.
- Inference: If the sign-chirality relation survives inclusion of higher-n bound electrons and electron-electron rescattering, the asymmetry could serve as a model-independent handedness diagnostic even when the overall rate is dominated by several competing operators.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives the energy-angular distribution of the two final electrons in the charged lepton flavor violating process μ−e− → e−e− in a muonic atom, for a polarized bound muon. The effective Lagrangian contains photonic dipole and six contact operators of different chiralities. The authors formulate the decay rate in terms of parity-violating asymmetry functions FS and FA and a motion-reversal-odd function FD, first in a plane-wave approximation for the contact interaction and then with a full multipole expansion using Dirac wave functions in a finite nuclear charge distribution. Numerical results for 208Pb are presented for four representative cases: g1, g3, g5, and AR, each under a single-operator dominance hypothesis. The main claims are that the sign of the asymmetry directly reflects the chirality of the muon, that the energy-angle dependence distinguishes operator types, and that FD is generated by final-state interactions even for real couplings at the level of 10^-2 to 10^-1.
Significance. If the central claims hold, this paper provides new observables beyond the total decay rate for identifying the chiral structure of CLFV interactions in muonic atoms. The analytic plane-wave formulas and the multipole expansion are internally cross-checked, and the calculation is parameter-free in the sense that no fitting is involved. The computation of the motion-reversal-odd correlation FD from final-state interactions is a useful, quantitative warning for future CP-violation searches. The analysis is limited by the single-operator dominance assumption and by the lack of truncation/convergence details, but the overall approach is credible and potentially important for COMET and similar programs.
major comments (2)
- [Section IV, finding (ii); Section III, Eqs. (35)-(38)] The claim that "the sign of the asymmetry directly reflects the chirality of the muon involved" is established only under the single-operator dominance hypothesis used in Eqs. (35)-(38). For the general effective Lagrangian of Eqs. (1)-(3), F and FD are ratios of bilinear products of amplitudes, so interference between opposite-chirality operators contributes to both numerator and denominator. A parity-conserving combination, e.g., g1=g2=1 or AR=AL=1, must have a vanishing P·p̂ term by parity, so the sign is not a direct readout of muon chirality once both chiralities are present. The paper does not quantify how large an opposite-chirality admixture can be before the sign relation or the FS/FA discrimination patterns break down. Please either state the restriction to single-operator dominance in the conclusion and abstract, or provide a robustness analysis with mixed chiralities.
- [Section II.B, Eqs. (28)-(31), Figs. 2-5] The numerical coefficients cF_l and cFD_l in Eqs. (30) and (31) are defined as infinite sums over κ1, κ2, κ1', κ2', J, J', and the expansions in Eqs. (28)-(29) are infinite Legendre series. The manuscript does not state the truncation order used for Figs. 2-5 or provide convergence tests. Because the sharp operator-dependent differences in FS and FA are the central quantitative results, the absence of this information makes the numerical predictions non-reproducible and leaves open the possibility of truncation artifacts. Please state the maximum |κ| and l included and show that the asymmetries are stable under increasing the cutoffs.
minor comments (4)
- [Section III, case 2, Eq. (36)] The condition "gj≠1 = 0" in Eq. (36) should read "gj≠3 = 0", since the case under consideration is g3 = 1.
- [After Eq. (5)] The restriction to 1s bound electrons is stated without a quantitative estimate; adding a sentence on the expected size of higher-n contributions would strengthen the independent-particle-model justification.
- [Section III, Fig. 5 discussion] The text says the effect of electron-electron rescattering is much smaller than electron-nucleus scattering because of the large nuclear charge, but no numerical estimate is given; a brief quantitative statement would be helpful for assessing the FD background.
- [Figures 2-5] The figures are only color maps; since the paper makes quantitative claims about the size of asymmetries (e.g., "large" for photonic, FD ~ 10^-2 to 10^-1), including numerical tables or contour values for representative kinematics would improve reproducibility.
Circularity Check
No circularity: the asymmetry coefficients are derived from first-principles matrix elements with couplings set by convention, not fitted to the predicted observables.
full rationale
The paper derives the polarized-muon angular distribution from the transition amplitude in Eq. (5) and the general CLFV effective Lagrangian in Eqs. (1)-(3). The asymmetry functions F, FD, FS, and FA are computed from bilinear products of amplitude matrix elements via Eqs. (9), (28), (29), and (32)-(34), with no parameter fitted to the target asymmetries. The single-operator dominance cases in Eqs. (35)-(38) set one coupling to unity and all others to zero; this is a stated calculational hypothesis for illustration, not an input that is later recovered as a prediction. The dependence of the asymmetry sign on the chirality coefficient ha emerges analytically in Eqs. (20) and (25), so the chiral-sign claim is a derived consequence of the operator structure rather than an assumption. Reliance on Refs. [11,12] for the unpolarized decay rate and radial matrix elements is normal citation of independent prior derivations, and the unpolarized formulas are reproduced in Appendix A, making the present calculation self-contained for the quantities it presents. The restriction to 1s bound electrons and the neglect of electron-electron final-state rescattering are explicit approximations, not circular inputs. No step in the derivation chain is equivalent by construction to the paper's own conclusions.
Assumptions & free parameters
assumptions (6)
- domain assumption The Dirac equation with the uniform nuclear charge distribution in Eq. (39) describes the bound and scattered lepton wave functions.
- domain assumption Only initial-state 1s bound muon and 1s bound electron are included; higher orbits are neglected.
- domain assumption Electron-electron final state interaction is neglected compared with electron-nucleus distortion.
- ad hoc to paper Single operator dominance: each numerical case sets exactly one CLFV coupling to 1 and all others to 0.
- domain assumption The effective Lagrangian in Eqs. (1)-(3) spans all relevant leading CLFV operators for this process.
- standard math Fierz identities correctly rewrite the g5 and g6 vector-vector operators as scalar-scalar operators.
Cite this review
Pith. "Pith review of Momentum distribution of the electron pair from the charged lepton flavor violating process $\mu^-e^-\to e^-e^-$ in muonic atoms with a polarized muon." pith.science (2026). https://pith.science/paper/OBRIDYIZ
@misc{pith2026190811653,
author = {Pith},
title = {Pith review of: Momentum distribution of the electron pair from the charged lepton flavor violating process $\mu^-e^-\to e^-e^-$ in muonic atoms with a polarized muon},
year = {2026},
howpublished = {\url{https://pith.science/paper/OBRIDYIZ}},
note = {Machine review of arXiv:1908.11653}
}
abstract
The $\mu^-e^-\to e^-e^-$ process in a muonic atom is one of the promising probes to study the charged lepton flavor violation (CLFV). We have investigated the angular distribution of electrons from the polarized muon of the atomic bound state. The parity violating asymmetric distribution of electrons is analyzed by using lepton wave functions under the Coulomb interaction of a finite nuclear charge distribution. It is found that the asymmetry parameters of electrons are very sensitive to the chiral structure of the CLFV interaction and the contact/photonic interaction. Therefore, together with the atomic number dependence of the decay rate studied in our previous work, the angular distribution of electrons from a polarized muon should be a very useful tool to constrain the model beyond the standard model.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Contact interaction with scalar coupling, where the electrons are emitted with the same chirality: g1 = 1, AL/R = 0, and gj⁄=1 = 0. (35)
-
[2]
Contact interaction with vector coupling, where the electrons are emitted with the same chirality: g3 = 1, AL/R = 0, and gj⁄=1 = 0. (36)
-
[3]
Contact interaction, where the electrons are emitted with the opposite chirality: g5 = 1, AL/R = 0, and gj⁄=5 = 0. (37)
-
[4]
(38) In the following, we show results of asymmetry coefficients for a polarized muon in 208Pb
Photonic interaction: AR = 1, AL = 0, and gj = 0. (38) In the following, we show results of asymmetry coefficients for a polarized muon in 208Pb. The muon and electron wave functions are obtained by solving a Dirac equation with a Coulomb potential of the uniform distribution of nuclear charge ρC(r), ρC(r) = 3Ze 4πR3θ(R−r), (39) 7 with R = 1.2A1/3fm, where ...
- [5]
- [6]
- [7]
- [8]
Show all 20 references
-
[9]
A. M. Baldini et al ., arXiv:1301.7225
- [10]
- [11]
- [12]
-
[13]
T. M. Nguyen, Pros. Sci. FPCP2015 (2015) 060
2015
-
[14]
Koike, Y
M. Koike, Y. Kuno, J. Sato, and M. Yamanaka, Phys. Rev. Lett. 105, 121601 (2010)
2010
-
[15]
Uesaka, Y
Y. Uesaka, Y. Kuno, J. Sato, T. Sato, and M. Yamanaka, Phys. Rev. D 93, 076006 (2016)
2016
-
[16]
Uesaka, Y
Y. Uesaka, Y. Kuno, J. Sato, T. Sato, and M. Yamanaka, Phys. Rev. D 97, 015017 (2018)
2018
-
[17]
Kuno and Y
Y. Kuno and Y. Okada, Phys. Rev. Lett. 77, 434 (1996)
1996
-
[18]
Okada, K
Y. Okada, K. I. Okumura, and Y. Shimizu, Phys. Rev. D 61, 094001 (2000)
2000
-
[19]
R. G. Sachs, The Physics of Time Reversal (University of Chicago Press, Chicago, 1985)
1985
-
[20]
S. Ando, J. A. McGovern, and T. Sato, Phys. Lett. B 677, 109 (2009)
2009
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.