Study of atomic effects on electron spectrum in bound-muon decay process

M. Y. Kaygorodov School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2G7    Y. S. Kozhedub Department of Physics, St. Petersburg State University, Universitetskaya 7-9, 199034 St. Petersburg, Russia    A. V. Malyshev Department of Physics, St. Petersburg State University, Universitetskaya 7-9, 199034 St. Petersburg, Russia Petersburg Nuclear Physics Institute named by B.P. Konstantinov of National Research Center “Kurchatov Institute”, Orlova roscha 1, 188300 Gatchina, Leningrad region, Russia    A. O. Davydov Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2G7    Y. Wu Institute of Applied Physics and Computational Mathematics, Beijing 100088, China    S. B. Zhang School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China
Abstract

For the bound-muon decay process, the study of atomic effects on the electron spectrum near its endpoint is performed within the framework of the Fermi effective theory. The analysis takes into account for corrections due to finite-nuclear-size, nuclear-deformation, electron-screening, and vacuum-polarization effects, all of which are incorporated self-consistently into the Dirac equation. Furthermore, the nuclear-recoil correction to the muon binding energy is included. Calculations are carried out for the isotopes of C, Al, and Si, which are of a particular importance for forthcoming experiments aimed at search for the charged-lepton flavor-violating process of muon-to-electron conversion in a nuclear field.

I Introduction

In contrast to the well-established phenomena of quark-generation mixing and neutrino-flavor oscillations, there is still no experimental evidence for the charged-lepton flavor violation (CLFV), see, e.g., the reviews in Refs. [1, 2, 3, 4] and references therein. One of promising channels for investigating the CLFV is the coherent muon-to-electron conversion in the field of a nucleus. As this process may involve hypothetical lepton-quark interactions, it can be used to test various extensions of the Standard Model. The current experimental upper limit on the conversion branching ratio is reported in Ref. [5]. Several next-generation experiments [6, 7, 8], based on the grounds of Ref. [9], aim to significantly improve the sensitivity to this rare process. The prospects of these experiments to explore physics beyond the Standard Model are discussed in Ref. [10].

The bound-muon decay, also referred to in the literature as the μ−superscript𝜇\mu^{-}italic_μ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT decay in orbit , is a source of physical background which accompanies the muon-to-electron conversion. The experiments on the conversion are focused on the electron energy spectrum near its endpoint, where the background contributions are expected to be minimal. However, in the case of bound-muon decay, the electron spectrum close to the endpoint is strongly affected by atomic effects. Consequently, to reliably identify the muon-to-electron conversion signal, it is essential to develop a precise theoretical model of the Standard Model-allowed background arising from bound-muon decay, with careful inclusion of all relevant atomic effects.

The theoretical study of the bound-muon decay process was initiated in Refs. [11, 12]. Later, various aspects of the problem were investigated in Refs. [13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. The most recent and comprehensive analysis of the complete electron spectrum across different nuclei was carried out in Refs. [23, 24], where fully relativistic wave functions constructed for the finite-nuclear-charge distribution model were employed. The specific case of Al1327superscriptsubscriptAl1327{}^{27}_{13}\mathrm{Al}start_FLOATSUPERSCRIPT 27 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 13 end_POSTSUBSCRIPT roman_Al was scrutinized in Refs. [25, 26, 27]. The isotope dependence of the electron spectrum near the endpoint was explored in Ref. [28].

The present study has two primary goals. The first objective is to derive a general expression for the electron spectrum corresponding to an arbitrary bound state of the muon. The second one is to study the influence of various atomic effects, including finite nuclear size, nuclear deformation, muon nuclear recoil, and electron screening, on both the muon binding energy and the resulting electron spectrum. In addition, corrections arising from quantum electrodynamics (QED), specifically the vacuum-polarization effect, are also considered. The analysis of atomic effects on the electron spectrum is performed for several isotopes: C   612superscriptsubscriptC612{}^{12}_{\,\,\,6}\mathrm{C}start_FLOATSUPERSCRIPT 12 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 6 end_POSTSUBSCRIPT roman_C and Si1428superscriptsubscriptSi1428{}^{28}_{14}\mathrm{Si}start_FLOATSUPERSCRIPT 28 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 14 end_POSTSUBSCRIPT roman_Si, which will be used in the experiment described in Ref. [6], and Al1327superscriptsubscriptAl1327{}^{27}_{13}\mathrm{Al}start_FLOATSUPERSCRIPT 27 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 13 end_POSTSUBSCRIPT roman_Al, which will be employed in the experiments discussed in Refs. [7, 8].

The paper is organized as follows. In Sec. II, a brief derivation of a general expression for the electron spectrum within the central-field approximation is presented, employing two independent but equivalent approaches to verify the consistency of the results. Sec. III outlines the details of the numerical calculations, presents and discusses the results, and provides a comparison with existing data from the literature. Sec. IV is reserved for the conclusions. The manuscript also contains seven appendices, each addressing specific technical aspects related to the derivations in Sec. II.

The relativistic units are used throughout the paper unless specified otherwise.

II Theory

The bound-muon decay process is described by the reaction

μ−→e−+ν¯e+νμ,→superscript𝜇superscript𝑒subscript¯𝜈𝑒subscript𝜈𝜇\mu^{-}\to e^{-}+\bar{\nu}_{e}+\nu_{\mu},italic_μ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT → italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT , (1)

where the initial-state muon μ−superscript𝜇\mu^{-}italic_μ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT is bound in a spherically symmetric nuclear potential, final-state electron e−superscript𝑒e^{-}italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT is unbound, and the final-state neutrinos ν¯esubscript¯𝜈𝑒\bar{\nu}_{e}over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and νμsubscript𝜈𝜇\nu_{\mu}italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT can be regarded as free massless particles. Instead of treating the reaction as a conventional three-body decay, it is advantageous to reformulate it as an equivalent two-to-two scattering process:

μ−+νe⁢(pνe)→e−+νμ⁢(pνμ),→superscript𝜇subscript𝜈𝑒subscript𝑝subscript𝜈𝑒superscript𝑒subscript𝜈𝜇subscript𝑝subscript𝜈𝜇\mu^{-}+\nu_{e}\left(p_{\nu_{e}}\right)\to e^{-}+\nu_{\mu}\left(p_{\nu_{\mu}}% \right),italic_μ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) → italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) , (2)

where the outgoing electron anti-neutrino ν¯esubscript¯𝜈𝑒\bar{\nu}_{e}over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT, characterized by an asymptotic four-momentum pν¯esubscript𝑝subscript¯𝜈𝑒p_{\bar{\nu}_{e}}italic_p start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT, is reinterpreted as an incoming neutrino with opposite four-momentum pνe=−pν¯esubscript𝑝subscript𝜈𝑒subscript𝑝subscript¯𝜈𝑒p_{\nu_{e}}=-p_{\bar{\nu}_{e}}italic_p start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - italic_p start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT. In Eq. (2), pνμsubscript𝑝subscript𝜈𝜇p_{\nu_{\mu}}italic_p start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT is an asymptotic four-momentum of the muonic neutrino νμsubscript𝜈𝜇\nu_{\mu}italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT, and its definition stays unchanged compared to Eq. (1). This formulation facilitates the application of standard in quantum-field-theory techniques.

Within the framework of the Fermi effective theory, the lepton-neutrino interaction operator governing the process described in Eq. (2) takes the following form:

VF⁢(1,2)=superscript𝑉F12absent\displaystyle V^{\mathrm{F}}\left(1,2\right)=italic_V start_POSTSUPERSCRIPT roman_F end_POSTSUPERSCRIPT ( 1 , 2 ) = GF2⁢δ⁢(r→1−r→2)subscript𝐺F2𝛿subscript→𝑟1subscript→𝑟2\displaystyle\frac{G_{\mathrm{F}}}{\sqrt{2}}\delta\left(\vec{r}_{1}-\vec{r}_{2% }\right)divide start_ARG italic_G start_POSTSUBSCRIPT roman_F end_POSTSUBSCRIPT end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG italic_δ ( over→ start_ARG italic_r end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - over→ start_ARG italic_r end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) (3)
×\displaystyle\times× [γ0⁢γρ⁢(1−γ5)]⁢(1)⁢[γ0⁢γρ⁢(1−γ5)]⁢(2),delimited-[]superscript𝛾0superscript𝛾𝜌1superscript𝛾51delimited-[]superscript𝛾0subscript𝛾𝜌1superscript𝛾52\displaystyle\left[\gamma^{0}\gamma^{\rho}\left(1-\gamma^{5}\right)\right]% \left(1\right)\left[\gamma^{0}\gamma_{\rho}\left(1-\gamma^{5}\right)\right]% \left(2\right),[ italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_ρ end_POSTSUPERSCRIPT ( 1 - italic_γ start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT ) ] ( 1 ) [ italic_γ start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( 1 - italic_γ start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT ) ] ( 2 ) ,

where the indices (1)1(1)( 1 ) and (2)2(2)( 2 ) here and in what follows label the particles the operator acts on, GFsubscript𝐺FG_{\mathrm{F}}italic_G start_POSTSUBSCRIPT roman_F end_POSTSUBSCRIPT is the Fermi constant, γρsuperscript𝛾𝜌\gamma^{\rho}italic_γ start_POSTSUPERSCRIPT italic_ρ end_POSTSUPERSCRIPT are the Dirac gamma matrices, and the summation over the repeated Lorentz indices is implied. The tree-level amplitude of the process is given by

A=⟨e⁢νμ|VF|μ⁢νe⟩,𝐴quantum-operator-product𝑒subscript𝜈𝜇superscript𝑉F𝜇subscript𝜈𝑒A=\braket{e\nu_{\mu}}{V^{\mathrm{F}}}{\mu\nu_{e}},italic_A = ⟨ start_ARG italic_e italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_ARG | start_ARG italic_V start_POSTSUPERSCRIPT roman_F end_POSTSUPERSCRIPT end_ARG | start_ARG italic_μ italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG ⟩ , (4)

where the bra-ket notation is employed and the written matrix element implies integration over the spatial and spin coordinates of the particles involved.

Given the spherical symmetry of the system, it is natural to adopt the central-field approximation. To obtain the electron spectrum, one has to integrate over the quantum numbers of neutrinos, which are not observed in the experiment, accounting for all their possible configurations consistent with the energy-conservation law, Eμ+Eνe=Ee+Eνμsubscript𝐸𝜇subscript𝐸subscript𝜈𝑒subscript𝐸𝑒subscript𝐸subscript𝜈𝜇E_{\mu}+E_{\nu_{e}}=E_{e}+E_{\nu_{\mu}}italic_E start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_E start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT. The operator in Eq. (3) ensures that only neutrino states with the appropriate helicity contribute to the transition amplitude. As a result, the integration domain over the neutrino variables in Eq. (4) can be formally extended without altering the outcome. By employing different representations of the neutrino states in Eq. (4), two different but formally equivalent approaches for evaluating the electron spectrum emerge. These alternative formulations lead to different final expressions, the equivalence of which can be demonstrated explicitly. For completeness, both approaches are presented and discussed in detail in this work.

The first method is a general, brute-force approach in which all particles are represented in the spherical-wave basis. For each particle, the basis functions are described by a set of relativistic central-field quantum numbers: ζ≡(E,ϰ,μ)𝜁𝐸italic-ϰ𝜇\zeta\equiv\left(E,\varkappa,\mu\right)italic_ζ ≡ ( italic_E , italic_ϰ , italic_μ ), where E𝐸Eitalic_E denotes the energy (or an equivalent quantum number), ϰitalic-ϰ\varkappaitalic_ϰ is the relativistic angular quantum number, and μ𝜇\muitalic_μ is the projection of the total angular momentum j=|ϰ|−1/2𝑗italic-ϰ12j=|\varkappa|-1/2italic_j = | italic_ϰ | - 1 / 2. For the bound muon, the principal quantum number n𝑛nitalic_n is used instead of the energy E𝐸Eitalic_E. The differential decay rate in this approach is given by

d⁢W2⁢p⁢(ζe,ζμ)d⁢Ee𝑑superscript𝑊2psubscript𝜁𝑒subscript𝜁𝜇𝑑subscript𝐸𝑒\displaystyle\frac{dW^{2\mathrm{p}}\left(\zeta_{e},\zeta_{\mu}\right)}{dE_{e}}divide start_ARG italic_d italic_W start_POSTSUPERSCRIPT 2 roman_p end_POSTSUPERSCRIPT ( italic_ζ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_ζ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ) end_ARG start_ARG italic_d italic_E start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG =2⁢π⁢∫δ⁢(−Ee−Eνμ+Eμ+Eνe)absent2𝜋𝛿subscript𝐸𝑒subscript𝐸subscript𝜈𝜇subscript𝐸𝜇subscript𝐸subscript𝜈𝑒\displaystyle=2\pi\int\,\delta\left(-E_{e}-E_{\nu_{\mu}}+E_{\mu}+E_{\nu_{e}}\right)= 2 italic_π ∫ italic_δ ( - italic_E start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) (5)
×|A⁢(ζe,ζνμ,ζμ,ζνe)|2⁢d⁢ζνμ⁢d⁢ζνe,absentsuperscript𝐴subscript𝜁𝑒subscript𝜁subscript𝜈𝜇subscript𝜁𝜇subscript𝜁subscript𝜈𝑒2𝑑subscript𝜁subscript𝜈𝜇𝑑subscript𝜁subscript𝜈𝑒\displaystyle\times\left|A\left(\zeta_{e},\zeta_{\nu_{\mu}},\zeta_{\mu},\zeta_% {\nu_{e}}\right)\right|^{2}d\zeta_{\nu_{\mu}}\,d\zeta_{\nu_{e}},× | italic_A ( italic_ζ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_ζ start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_ζ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT , italic_ζ start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_ζ start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_d italic_ζ start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT ,

where it is assumed that all unbound sates are normalized to the delta function in energy. This formulation will be referred below to as the two-particle approach and labeled with ”2p“.

The second approach represents only the muon and electron states in the spherical-wave basis, while the neutrinos are described using the plane-wave basis, characterized by a set of quantum numbers ξ=(p→,ms)𝜉→𝑝subscript𝑚𝑠\xi=\left(\vec{p},m_{s}\right)italic_ξ = ( over→ start_ARG italic_p end_ARG , italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ), where p→→𝑝\vec{p}over→ start_ARG italic_p end_ARG is the three-momentum and mssubscript𝑚𝑠m_{s}italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT is the spin projection onto p→→𝑝\vec{p}over→ start_ARG italic_p end_ARG. The latter wave functions are normalized to the delta function in momentum space. In this approach, the partial-wave electron spectrum can be expressed as

d⁢W1⁢p⁢(ζe,ζμ)d⁢Ee𝑑superscript𝑊1psubscript𝜁𝑒subscript𝜁𝜇𝑑subscript𝐸𝑒\displaystyle\frac{dW^{1\mathrm{p}}\left(\zeta_{e},\zeta_{\mu}\right)}{dE_{e}}divide start_ARG italic_d italic_W start_POSTSUPERSCRIPT 1 roman_p end_POSTSUPERSCRIPT ( italic_ζ start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , italic_ζ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ) end_ARG start_ARG italic_d italic_E start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG =2⁢π⁢∫δ⁢(−Ee−Eνμ+Eμ+Eνe)absent2𝜋𝛿subscript𝐸𝑒subscript𝐸subscript𝜈𝜇subscript𝐸𝜇subscript𝐸subscript𝜈𝑒\displaystyle=2\pi\int\,\delta\left(-E_{e}-E_{\nu_{\mu}}+E_{\mu}+E_{\nu_{e}}\right)= 2 italic_π ∫ italic_δ ( - italic_E start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) (6)