---
abstract: 'We propose a new experiment to search for a time-reversal $(T)$ symmetry breaking process in muon decay and the Majoranality of the neutrinos. In the presence of $V+A$ interactions, the Majoranality appears as a $T$-violating term in the muon decay width as shown by Doi et al [@Doi1982], while in the Standard Model such a $T$-violating term is negligibly small. The presences of $V+A$ interactions and the corresponding heavy right-handed Majorana neutrinos give us an important clue to solve two major issues in particle physics, the deficit of baryon asymmetry in the universe and the Majoranality of neutrinos. In the experiment, the polarization of positrons from muon decays is measured using a polarimeter consisting of a magnetized foil and a segmented calorimeter. According to our result of numerical calculation, a factor of ten improvement in sensitivity to the $T$-violating process is expected by a year of measurement at J-PARC Materials and Life Science Experimental Facility, compared to the most recent precursor experiment.'
author:
- Takeshi Fukuyama
- Sohtaro Kanda
- Daisuke Nomura
- Koichiro Shimomura
bibliography:
- 'bib.bib'
title: 'Hunting for $T$-violation and Majoranality of Neutrinos in Muon Decays'
---

Introduction
============

Just after the discovery of parity violation [@Lee1952; @Wu1957], possible tests of time reversal invariance ($T$-invariance) have been proposed [@Jackson1957]. Violation of $T$-invariance immediately implies $CP$-violation if invariance under $CPT$ transformation is assumed. $CP$-violation was afterward first found in kaon decays [@Christenson1964] and later in B-meson systems [@Belle2001; @Babar2001]. These phenomena are successfully explained by the Kobayashi-Maskawa phase in the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix [@Kobayashi1973; @Cabibbo1963]. Nevertheless, the observed baryon asymmetry in the universe (BAU) strongly indicates the existence of additional $CP$-violating factors beyond that of the Standard Model (SM).

Also, in the lepton sector, we have recently obtained great achievements in understanding the generation structure, including the determination of the Maki-Nakagawa-Sakata (MNS) mixing matrix elements [@Maki1962; @Pontecorvo1957]. In these situations, it is one of the most important unsolved issues in particle physics whether the neutrinos are Majorana or Dirac particles.

The most attractive idea to explain the tiny masses of the neutrinos is the seesaw mechanism via superheavy right-handed Majorana neutrinos [@Minkowski1977; @Yanagida1979; @GellMann1980; @Glashow1980; @Mohapatra1980], indicating that the light neutrinos are Majorana particles. These superheavy neutrinos may have induced the observed baryon asymmetry via leptogenesis [@Fukugita1986]. In this scenario, decays of the superheavy neutrinos induce asymmetry in the lepton number, which is later converted to the baryon number asymmetry when the temperature of the universe is ${\cal O}(10) \text{TeV}$ due to the thermal fluctuation between vacua with different $B+L$ numbers [@Kuzmin1985]. Thus, the superheavy right-handed neutrinos may solve the two major unsolved issues (the BAU and the neutrino Majoranality) mentioned above simultaneously. In this scenario, the existence of the right-handed weak current $V+A$ at very high energy is naturally expected since the right-handed neutrinos are most naturally predicted in many SO(10) grand unified theories (GUTs), where symmetry between the $V+A$ and $V-A$ currents are automatically predicted at high enough energies. This is also a robust prediction from many other GUTs.

In the presence of $V+A$ interactions, muon decays have a special implication since two emitted neutrinos interfere differently depending on whether they are Majorana or Dirac particles. Fortunately, the interference effect appears most clearly in the $CP$-violating term essentially without contamination of the SM background [@Doi1982]. Thus the search for $T$-violation in muon decay from both theoretical and experimental sides is a very attractive target to solve the above-mentioned problems, and is the motivation of this paper.

This paper is organized as follows. In Section \[section:section2\], we estimate the size of the $T$-violating term in the presence of $V+A$ current. It gives the relationship between $T$-violation and the Majorana property of the neutrinos. In Section \[section:section3\] we discuss the $\nu$SM (SM+massive neutrinos) background on the $T$-violating term which appears from radiative corrections. This is indispensable in searching for effects from new physics beyond the SM (BSM physics). In the subsequent sections, we give an experimental proposal to detect the $T$-symmetry breaking term in muon decay and the Majoranality of the neutrinos. In Section \[section:section4\], a measurement principle to detect the transverse polarization of the positrons emitted from muon decays is described. In Section \[section:section5\], we study practical methods and feasibility of the proposed experiment. Section \[section:section6\] is devoted to discussions, where the relation with another typical Majoranality experiment, neutrinoless double beta decay, is also argued. Section \[section:section7\] is the summary.

T-violating muon decays in New Physics beyond the SM {#section:section2}
====================================================

Before discussing muon decays in the presence of the $V+A$ current together with the usual $V-A$ one, we discuss the $V+A$ current in GUTs. As a BSM physics, GUTs are most promising and give comprehensive pictures. Among them, SO(10) GUT [@Fukuyama2013] is general enough. In the process of the symmetry breaking from SO(10) GUT to the SM, the right-handed current naturally appears in a symmetry breaking sequence like [^1] $$\begin{aligned}
\rm SO(10)\to SU(4)_C\otimes SU(2)_L\otimes SU(2)_R\to \nonumber \\
\rm \text SU(3)_C\otimes SU(2)_L\otimes SU(2)_R\otimes U(1)_{B-L}\to SM.\end{aligned}$$ This sequence shows that as the energy scale goes down, right-handed gauge bosons become massive at higher energy scale than the electroweak scale ($v_{\text{EW}}$) and the weak (left-handed) gauge bosons do so at $v_{\text{EW}}$. That is, left- and right-handed chiral Lagrangians successively appear as the energy scale goes up.

We should note that it is not clear which stage in the sequence above is directly related to the string scale where the gravitation is unified. For example, if we consider $ \text{SU}(4)_{\text C}\otimes \text{SU}(2)_{\text L}\otimes \text{SU}(2)_{\text R}\otimes \text{D}$ (Pati-Salam (PS) phase [@Pati1974] with D parity which interchanges left (L)- and right (R)-handed currents) as a remnant from SO(10), their coupling constants $g_{\text L},\ g_{\text R}$ and mixing matrices $U_{lj},\ V_{lj}$ satisfy following equations at the mass scale of GUT ($M_{\text{GUT}}$) $$\begin{aligned}
g_{\text L}=g_{\text R}, ~~U_{lj}=V_{lj}.
\label{PS}\end{aligned}$$ The mixing matrices $U_{lj}$ and $V_{lj}$ connect the weak $(l)$ and mass $(j)$ eigenstates. The former and the latter mixing matrices correspond to the left- and right-handed sectors, respectively. The constraints in Eq. [(\[PS\])]{} are realized at $M_{\text{GUT}}$ but start to be violated as the energy goes down to the SM scale by renormalization group equations.

However, if the PS phase is the stage which is directly related to the string scale, we are free from the above constraints. Furthermore, if we start from $\text{SU}(3)_{\text C}\otimes \text{SU}(2)_{\text L}\otimes \text{SU}(2)_{\text R}\otimes \text{U}(1)_{\rm{B-L}}$, the quarks and leptons are separated and there are no leptoquark-type gauge bosons which connect the quarks and leptons. If we forget about GUTs but still consider the same mixing matrix as in Eq. [(\[PS\])]{}, the lower limit of $M_{\text{WR}}$ (mass of right-handed gauge boson) is relaxed to $M_{\text{WR}}>1.6$ TeV [@Beall1982]. Thus the right-handed currents appear in many BSM physics, where the $V+A$ interactions coexist with the $V-A$ interactions.

Another theoretical focus of this paper is the Majoranality of neutrinos, which is the natural result of the seesaw mechanism of heavy right-handed neutrinos and Dirac neutrinos. Thus, as BSM physics, we consider the general effective interaction which includes the right-handed current $$\begin{aligned}
H_{\text W}=\frac{G_{\text F}}{\sqrt{2}}\left[j_{{\text L}\alpha}^\dagger j_{\text L}^\alpha
   +\lambda j_{{\text R}\alpha}^\dagger j_{\text R}^\alpha
   +\kappa\left(j_{{\text L}\alpha}^\dagger j_{\text R}^\alpha
   +j_{{\text R}\alpha}^\dagger j_{\text L}^\alpha\right)\right],
\label{JLJR}\end{aligned}$$ where $j_{{\text L}\alpha}$ and $j_{{\text R}\alpha}$ are given by $$\begin{aligned}
j_{{\text L}\alpha} &\equiv 
 \sum_{l=e,\mu,\tau}\overline{l}(x)\gamma_\alpha(1-\gamma_5)
 \nu_{l{\text L}}(x),\\
j_{{\text R}\alpha} &\equiv \sum_{l=e,\mu,\tau}\overline{l}(x)\gamma_\alpha
   (1+\gamma_5)\nu_{l{\text R}}'(x).\end{aligned}$$ Here $\nu_{\text L}~(\nu_{\text R}')$ are the left-handed (right-handed) weak eigenstates of neutrinos. The constants $\lambda$ and $\kappa$ are related to the mass eigenvalues of the weak bosons in the left- and right-handed gauge sector ($W_{\text L}, W_{\text R})$ as follows: $$\begin{aligned}
 \label{LR1}
W_{\text L}&=&W_1\cos\zeta+W_2\sin\zeta,\\
W_{\text R}&=&-W_1\sin\zeta+W_2\cos\zeta,\\
\frac{G_{\text F}}{\sqrt{2}}&=&\frac{g^2}{8}\cos^2\zeta \frac{M_1^2\tan^2\zeta+M_2^2}{M_1^2M_2^2}\label{LR2},\\
\lambda &\equiv& \frac{M_1^2+M_2^2\tan^2 \zeta}{M_1^2\tan^2\zeta+M_2^2} \label{LR3},\\
\kappa &\equiv& \frac{-(M_2^2-M_1^2)\tan \zeta}{M_1^2\tan^2\zeta+M_2^2}\label{LR4},\end{aligned}$$ where $M_1$ and $M_2$ are the masses of the mass eigenstates $W_1$ and $W_2$, respectively, and $\zeta$ is the mixing angle which relates the mass eigenstates and the gauge eigenstates.

Using the effective Hamiltonian, the general form of the muon differential decay rate in the rest frame of $\mu^\mp$ is given by [@Doi1982] $$\begin{aligned}
\frac{d\Gamma}{d{\bf q}_e}
&=\frac{m_\mu G_{\text F}^2}{3(2\pi)^4} \bigg\{ N(e)
 \pm\frac{ \left( \right. {\bf q}_e\cdot \boldsymbol{\zeta}_\mu \left. \right)}{E}P(e)
 \mp\frac{ \left( \right. {\bf q}_e\cdot \boldsymbol{\zeta}_e \left. \right)}{E}Q(e) \nonumber \\
&-\frac{({\bf q}_e\times \boldsymbol{\zeta}_\mu)\cdot
  \left( \right. {\bf q}_e\times \boldsymbol{\zeta}_e \left. \right)}{|{\bf q}_e |^2}R(e)
-\frac{ \left( \right. {\bf q}_e\cdot \boldsymbol{\zeta}_\mu \left. \right)
 \left( \right. {\bf q}_e\cdot \boldsymbol{\zeta}_e \left. \right)}{|{\bf q}_e |^2}S(e) \nonumber \\
&+ \frac{\boldsymbol{\zeta}_\mu\cdot
 \left( \right. {\bf q}_e\times\boldsymbol{\zeta}_e \left. \right)}{E}T(e) \bigg\},
\label{eq:drate}\end{aligned}$$ where $E$ and ${\bf q}_e$ are the energy and momentum of the emitted electron (positron), and $\boldsymbol{\zeta}_e$ and $\boldsymbol{\zeta}_\mu$ are the polarization vectors of $e^{\mp}$ and $\mu^{\mp}$, respectively. In Eq. (\[eq:drate\]) we assume that we detect $e^{\mp}$ in the final state but not the neutrinos. The upper and lower signs in Eq. [(\[eq:drate\])]{} correspond to $\mu^-$ and $\mu^+$, respectively. This comes from the difference of the density matrix, $$\begin{aligned}
\rho=\frac{1}{2}(\gamma^\mu q_\mu \pm m)(1+\gamma_5 \gamma^\nu a_\nu)~. \end{aligned}$$ Here the double-sign “$\pm$” corresponds to Eq. [(\[eq:drate\])]{}. The four-dimensional spin vector $a_\mu$ is [@BLP] $$\begin{aligned}
{\bf a}=\boldsymbol{\zeta}
+\frac{{\bf q}(\boldsymbol{\zeta}\cdot{\bf q})}{m(\epsilon+m)}, ~~~~
a^0=\frac{{\bf q}\cdot\boldsymbol{\zeta}}{m},\end{aligned}$$ where $\boldsymbol{\zeta}$ is twice the mean spin vector in the rest frame and $\epsilon$ is the energy. $N(e)$, $S(e)$, and $R(e)$ are $P$-even terms. $P(e)$ and $Q(e)$ are $P$-odd terms. The $T(e)$ term is $P$-odd and $T$-odd ($i.e.$, $CP$-odd) term. Its explicit form is given by Eq. [(\[Te\])]{}, whose presence indicates that the neutrinos are Majorana particles in the scheme of our $(V+A)$ $\oplus$ $(V-A)$ model. The dominant $\lambda$ term in Eq. [(\[Te\])]{} comes from the interference of the diagrams of Fig. \[fig:muondecay\_intrfr\] (a) and (b) with $(A,B,C,D)=(L,L,R,R)$ and $(R,R,L,L)$. Here $L\equiv V-A$ and $R\equiv V+A$. This interference takes place only if the neutrinos are Majorana particles: Indeed, these diagrams indicate the same reactions for Majorana neutrinos. The $T(e)$ term indicates BSM physics up to very tiny corrections due to the $\nu$SM.

![ Tree level diagrams for muon decays. In the diagrams, $A, \ldots, D$ denotes the handedness of the currents, $L$ or $R$. In the set of the diagrams with $(A,\ldots,D)=(L,L,R,R)$, the neutrinos with the subscript $j$ in (a) and (b) have the same helicity, and those with the subscript $k$ also have the same helicity, both in the massless limit. The same statement also holds for the set with $(A,\ldots,D)=(R,R,L,L)$. So the diagrams of each set interfere for the case of the Majorana neutrinos.[]{data-label="fig:muondecay_intrfr"}](muondecay_intrfr_v3.eps){width="\linewidth"}

The explicit forms of $N(e),...,T(e)$ are given in [@Doi1982] and we give the only necessary term $T(e)$ in this paper, $$\begin{aligned}
T(e)=&-\epsilon_M\lambda h^I(m_\mu-m_e^2/m_\mu) \nonumber \\
&-3\kappa\left(v^Im_e-(v^I-4\kappa z^I)m_\mu\right) \nonumber \\
\approx &- \epsilon_M\lambda h^I(m_\mu-m_e^2/m_\mu),
\label{Te}\end{aligned}$$ where $\epsilon_M=1$ and $\epsilon_M=0$ correspond to the cases that the neutrinos are Majorana and Dirac particles, respectively. The second equality comes from the fact that $v^I$ and $z^I$, whose explicit values are given in [@Doi1981], are suppressed by ${m_i}/{m_\mu}$. Here, $m_i$ is the light neutrino mass, and are safely neglected except for around the maximum positron energy (see Eq. (\[Emax\])). The parameter $h^I$ is defined as $$\begin{aligned}
h^I\equiv \epsilon_M \sum_{j,k}F_{jk}^{1/2}G_{jk}\text{Im} C_{jk}'\approx \epsilon_M \sum_{j,k}\text{Im} C_{jk}',
\label{hI}\end{aligned}$$ where $$\begin{aligned}
\label{Fjk}
F_{jk}&=[1-(m_j+m_k)^2/\Delta^2][1-(m_j-m_k)^2/\Delta^2],\\
\label{Gjk}
G_{jk}&=1+(m_j^2+m_k^2)/\Delta^2-2(m_j^2-m_k^2)^2/\Delta^4,\\
\label{Cjk}
C'_{jk}&=U_{ej}^*V_{ek}V_{\mu j}^*U_{\mu k}\end{aligned}$$ with $$\begin{aligned}
\Delta^2 \equiv 2m_\mu (W-E). \end{aligned}$$ Here $W$ is the maximum positron energy in the limit of $m_\nu=0$, $$\begin{aligned}
W=\frac{m_\mu^2+m_e^2}{2m_\mu}.
\label{Emax}\end{aligned}$$ The imaginary part of $C_{jk}'$ is the “generalized" Jarlskog invariant. The left and right-handed mixing matrices, $U$ and $V$, respectively, connect the weak eigenstates with their mass eigenstates as $$\begin{aligned}
\nu_{l{\text L}}=\sum_{j=1}^{6}U_{lj}N_{j {\text L}},~~\nu_{l{\text R}}'=\sum_{j=1}^{6}V_{lj}N_{j {\text R}}.
\label{mixing}\end{aligned}$$ That is, the Jarlskog invariant of the CKM matrix [@Jarlskog1985] is replaced by that of the MNS mixing matrices of the left and right-handed neutrinos [@Maki1962]. This $T(e)$ term violates $T$-invariance. However, more importantly, only this term among $N(e),...,T(e)$ terms eventually vanishes if the neutrinos are not Majorana and is not contaminated by the $\nu$SM contributions. The magnitude of $\lambda$ is of order $(M_{\text{WL}}/M_{\text{WR}})^2$. The present lower limit of $M_{\text {WR}}$ is obtained from $\Delta L=2$ process [@Mandal2017] and of order $5$ TeV and $\lambda\leq 2.6\times 10^{-4}$ under the assumption of $V_{Rij}=\delta_{ij}$ [@PDG2016]. There still remains much larger possibility to detect the Majoranality than in the $\nu$SM case which will be shown in Eq. [(\[SM\])]{}.

Here we add some comments on the $T$-violating term, $\boldsymbol{\zeta}_\mu\cdot({\bf q}_e\times \boldsymbol{\zeta}_e)$. This term itself may appear as a result from new interactions other than $V+A$. For instance, if we generalize the Dirac bi-spinor coupling to $S$, $P$, $V$, $A$, $T$ (scalar, pseudo-scalar, vector, axial-vector, and tensor, respectively) using the Michel parameters [@Michel1950], such a term appears even in the case of Dirac neutrinos. In our theory, a new physical entity is restricted in the heavy right-handed current and $T(e)$ term does not depend on the QCD vacuum angle. So the compatibility of the non-null result is immediately checked with other phenomena.

If there is the right-handed current, the lepton electric dipole moment (EDM) appears mainly from one loop diagram as shown in Fig. \[fig:LRedm\], whereas it appears from four loop diagrams in the SM.

![The EDM diagram of the lepton $l_i$ $(i=1,\ldots,3)$ which appears at one-loop.[]{data-label="fig:LRedm"}](EDM-1loop_v2.eps){width="\linewidth"}

The lepton EDM of Fig. \[fig:LRedm\] is estimated [@Bernreuther1991] with negligible generation mixing to be $$\begin{aligned}
|d_{i}| <  8.2\times 10^{-27}\,
   \frac{|\text{Im}(m_{\text D})_{ii}|}{\text{MeV}}~ \text{e cm}, \end{aligned}$$ where $(m_{\text D})_{ii}$ denotes the $i$-th generation of the Dirac mass of the neutrino ($m_{\text D} \overline{\nu_{\text L}} \nu_{\text R}$). If $\text{SU}(3)_{\text C} \otimes \text{SU}(2)_{\text L} \otimes \text{SU}(2)_{\text R} \otimes \text{U}(1)_{\rm{B-L}}$ appears after the breaking of SO(10), $m_{\text D}$ is of order of up-type quark masses and the muon EDM may be significantly enhanced relative to the electron EDM. The muon EDM is searched for at the muon $g-2$/EDM experiments using muon storage rings at Brookhaven National Laboratory (BNL) [@Bennett2006], Fermi National Accelerator Laboratory (FNAL) [@Grange2015], and Japan Proton Accelerator Research Complex (J-PARC) [@Abe2019].

$T$-violating muon decays in Standard Model {#section:section3}
===========================================

To detect the signal discussed in the previous section, we have to clarify the main background processes to the signal process. Let us estimate the order of $T$-violating effect in the muon decay in the $\nu$SM framework.

![(a) Tree level diagram of muon decays. (b) An example of the diagrams which violate $CP$. $CP$-violating term in muon decay rate appears as the interference term of (a) and (b).[]{data-label="fig:CPVSM"}](CPVmuondecay_v2.eps){width="\linewidth"}

The main $T$-violating effect appears e.g. from the diagram (b) of Fig. \[fig:CPVSM\], which interferes with the tree-level diagram, Fig. \[fig:CPVSM\] (a) [^2]. This interference results in a $CP$-violating term in the muon differential decay width.

It is straightforward to estimate the size of the interference term between the two diagrams in Fig. \[fig:CPVSM\]. The $T$-violating decay width is suppressed by the factor of $\delta$ relative to the usual muon decay width, where $$\delta=\frac{G_{\text F} m_t^2}{8\pi^2\sqrt{2}}
       \frac{\alpha_s}{4\pi}\frac{1}{16\pi^2}\frac{g^2}{8}J_{CP}
\approx 1\times 10^{-13}.
\label{SM}$$ Here $J_{CP}$ is the Jarlskog parameter [@Jarlskog1985], which is a measure of $CP$-violation due to the CKM matrix [@Kobayashi1973], and $J_{CP} \approx 3.0 \times 10^{-5}$ [@PDG2016]. We find that the $T$-violating effect in the muon decay in the $\nu$SM framework is extremely small.

Experiment {#section:section4}
==========

We can search for a Majorana neutrino by measuring the transverse spin polarization of a positron from muon decay. In this section, an experimental principle and theoretical background of the experiment are described.

Polarization of a positron from muon decay
------------------------------------------

Figure \[fig:exp\_pol\] illustrates the definition of the $(x,y,z)$-coordinates in the muon rest frame. The positron polarization vector $\boldsymbol{\zeta}_e$ is described by three orthogonal components $$\begin{aligned}
\boldsymbol{\zeta}_e=(P_{\text T1},P_{\text T2}, P_{\text L}) \equiv (\bm{P}_{\text T}, P_{\text L}),\end{aligned}$$ where $P_{\text L}$ is the longitudinal component, $P_{\text T1}$ and $P_{\text T2}$ are transverse components. The $z$-axis is defined in parallel to the positron momentum. In the SM, the positron polarization vector $\boldsymbol{\zeta}_e$ is coplanar with the positron momentum ${\bf q}_e$ and the muon spin $\boldsymbol{\zeta}_\mu$. In other words, only $P_{\text L}$ and $P_{\text T1}$ have non-vanishing norms. A non-zero $P_{\text T2}$ indicates the $CP$ violation in muon decay.

![Definition of the spin polarization of a positron from muon decay: (a) The muon spin $\boldsymbol {\zeta_\mu}$ and the momentum of decay positron ${\bf q}_e$ are on the $x$-$z$ plane; (b) On the transverse polarization, the $x$ and $y$ components are defined as $P_{\text T1}$ and $P_{\text T2}$, respectively. A non-zero $P_{\text T2}$ implies the $CP$ violation in muon decay.[]{data-label="fig:exp_pol"}](exp_pol.eps){width="\linewidth"}

Polarimeter
-----------

A precision polarimeter which utilizes the transverse spin-dependent interaction is required for the experiment. For positron polarimetry, five methods have been established: Mott scattering, Bhabha scattering, Compton scattering, positron annihilation-in-flight, and Compton transmission. Each technique utilizes the spin-dependent cross section of a particular electromagnetic process. The polarimeters are qualified by an analyzing power $$\begin{aligned}
{\cal{A}} = \frac{N_+ - N_-}{N_+ + N_-},\end{aligned}$$ where $N_+$ and $N_-$ are the number of positrons polarized in parallel and anti-parallel to the axis of interest, respectively. Statistical precision of the polarimetry is evaluated by the figure of merit (FOM) $$\begin{aligned}
\label{eq:FOM}
F = N {\cal A}^2,\end{aligned}$$ where $N$ is the total number of positrons, $N=N_{+} + N_{-} $. Typical kinetic energy of incident positron ranges from 30 MeV to 50 MeV.

The analyzing power of polarimeter depends on positron energy. In the energy region of interest, the polarimeters utilizing Bhabha scattering and annihilation-in-flight are feasible. The details of these polarimeters are described in following sub-sections.

Bhabha scattering
-----------------

The Bhabha polarimeter is based on the spin-spin interaction between a positron spin and an electron spin. A thin foil of magnetized iron is commonly employed as a scattering target. It is sensitive to both transverse and longitudinal polarization components. In the high energy limit, the differential cross-section of the Bhabha scattering in the center-of-mass system is [@Swartz1988] $$\begin{aligned}
\label{Bhabha}
\left(\frac{d \sigma}{d \Omega}\right)_{\rm BHA} 
&= \frac{\alpha^2 (3+\cos^2 \theta^*)^2}{4s (1-\cos \theta^*)^2}
\left( \right. 1 - P_{\rm L}^{1} P_{\rm L}^{2} A_{\rm L}(\theta^*)  \nonumber \\
&- |\bm{P}_{\rm T}^{1}| | \bm{P}_{\rm T}^{2}| A_{\rm T}(\theta^*) 
   \cos(2 \phi - \phi_1 - \phi_2) \left. \right).\end{aligned}$$ Here $\alpha$ is the fine structure constant, $\theta^*$ is the scattering angle in the center-of-mass frame, $\sqrt{s}$ is the center-of-mass energy, and $\phi$ is the azimuth of the scattered positron. The superscripts $1$ and $2$ of the polarizations indicate the ones of beam and target, respectively. Figure \[fig:exp\_scat\] depicts the scattering plane. The azimuths $\phi_1$ and $\phi_2$ are the angles of the transverse polarizations $\bm{P}_{\rm T}^{1}$ and $\bm{P}_{\rm T}^{2}$ from the $y$-axis, respectively.

![Schematic of Bhabha scattering: (a) description of the scattering plane; (b) polarization vectors of an incident positron and a target electron. The $xyz$-frame is defined in Fig. \[fig:exp\_pol\].[]{data-label="fig:exp_scat"}](exp_scat.eps){width="\linewidth"}

The longitudinal and transverse asymmetries $A_{\text L}$ and $A_{\text T}$ are determined as follows $$\begin{aligned}
A_{\text L} = \frac{(7+ \cos^2 \theta^*)\sin^2 \theta^*}{(3+\cos^2 \theta^*)^2},
~~
A_{\text T} = \frac{\sin^4 \theta^*}{(3+\cos^2 \theta^*)^2}.\end{aligned}$$ The analyzing power is energy independent and maximal at $\theta^*=90$ degrees in the center-of-mass system for both longitudinal and transverse polarizations. The ratio of the transverse asymmetry to the longitudinal one depends on $\theta^*$ and $\phi$ for given $\phi_1$ and $\phi_2$ as $$\begin{aligned}
 \frac{A_{\text T}(\theta^*) \cos (2 \phi - \phi_1 - \phi_2)}
      {A_{\text L}(\theta^*)} 
= \frac{\sin^2 \theta^*}{7 + \cos^2\theta^*}
  \cos(2 \phi - \phi_1 - \phi_2),
\label{eq:AT_AL_ratio}\end{aligned}$$ which is the maximum at $(\theta^*, \phi)=(\pi/2, (\phi_1 + \phi_2)/2)$.

In the rest frame of the initial-state electron, which we call the lab frame for short, the scattering angle $\theta_{\text{lab}}$ of the positron is related to $\theta^*$ by $$\begin{aligned}
 \theta_{\text{lab}} = \arctan \left( \frac1{\gamma} \tan \frac{\theta^*}{2} \right),\end{aligned}$$ where $\gamma \equiv \sqrt{s}/(2m_e)$. The angle $\theta^*=\pi/2$ corresponds to $\theta_{\text{lab}}=\arctan(2m_e/\sqrt{s})$. It follows that, in the lab frame, both the analyzing power and the ratio Eq. (\[eq:AT\_AL\_ratio\]) take their maximum values at $(\theta_{\text{lab}}, \phi)= (\arctan(2m_e/\sqrt{s}), (\phi_1+\phi_2)/2)$. Later in this paper, feasibility of a new experiment to detect the positron’s transverse polarization will be discussed. There, the analyzing power and FOM of the polarimeters will be evaluated by using the GEANT4 software toolkit [@GEANT4_1; @GEANT4_2; @GEANT4_3]. Since interactions relevant to the longitudinal polarization are well established in GEANT4, firstly, we compute the longitudinal asymmetry for each small solid angle element $d\Omega = d(\cos \theta^*) d\phi$ and multiply the result by the factor $A_{\text T}(\theta^*) \cos (2 \phi - \phi_1 - \phi_2)/A_{\text L}(\theta^*)$ to calculate the transverse asymmetry. In the lab frame, this conversion factor becomes largest for the very forward region, $(\theta_{\text{lab}}, \phi) \simeq  (\arctan(2m_e/\sqrt{s}), (\phi_1+\phi_2)/2)$.

Annihilation-in-flight
----------------------

The annihilation-in-flight (AIF) polarimeter is based on the spin-dependent cross section of the positron annihilation. A scattering target is magnetized iron foil similarly to the case for the Bhabha polarimeter. The differential cross-section of AIF in the center-of-mass system in high energy limit is given by $$\begin{aligned}
\left(\frac{d\sigma}{d\Omega}\right)_{\rm AIF}
&=\frac{\alpha^2(1+\cos^2\theta^*)}
{s \sin^2\theta^*}\left( \right.
  1 + P_{\text L}^{1}P_{\text L}^{2}B_{\text L}(\theta^*) \nonumber \\
&+ |\bm{P}_{\text T}^{1}||\bm{P}_{\text T}^{2}|B_{\text T}(\theta^*)\cos (2\phi-\phi_1-\phi_2)
\left. \right),\end{aligned}$$ where the definitions of $\theta^*$ and $\phi_i$ are same as those in Eq. (\[Bhabha\]). Here $\phi$ is the azimuthal angle of one of the emitted photon. Without loss of generality, we can take the range of $\phi$ as $-90<\phi/{\text {deg.}} \leq 90 $. The functions $B_{\text L}(\theta^*)$ and $B_{\text T}(\theta^*)$ are defined by the references [@Page1957; @McMaster1961] $$\begin{aligned}
B_{\text L}(\theta^*)=1, 
~~
B_{\text T}(\theta^*)=\frac{\sin^2\theta^*}{1+\cos ^2\theta^*}.\end{aligned}$$ The AIF and Bhabha measurements can be performed simultaneously with the common apparatus. In contrast to the Bhabha polarimeter, the analyzing power of the AIF polarimeter for longitudinal and transverse polarizations are comparable. The analyzing power of the AIF polarimeter was calculated by Fetscher in detail [@Fetscher2007]. However, we emphasize that $P_{\text T2}$ appears even if his $F_{\text T2}$ vanishes because of the $V+A$ interaction and the Majoranality of the neutrinos.

Similarly to the case of Bhabha scattering, the ratio of the transverse asymmetry to the longitudinal one is given by $$\begin{aligned}
 \frac{B_{\rm T}(\theta^*) \cos (2 \phi - \phi_1 - \phi_2)}
      {B_{\rm L}(\theta^*)} 
= \frac{\sin^2 \theta^*}{1 + \cos^2\theta^*}
  \cos(2 \phi - \phi_1 - \phi_2),
\label{eq:BT_BL_ratio}\end{aligned}$$ which takes the largest value of unity at $(\theta^*, \phi)=(\pi/2, (\phi_1 + \phi_2)/2)$. The scattering angle $\theta_{\text{lab}}$ of the positron in the lab frame is related to $\theta^*$ by $$\begin{aligned}
 \theta_{\text{lab}} 
= \arctan \left( \frac{\sin \theta^*}{\gamma(\beta + \cos \theta^*)} \right),\end{aligned}$$ where $\gamma \equiv \sqrt{s}/(2m_e)$ and $\beta \equiv \sqrt{1-4m_e^2/s}$. The angle $\theta^*=\pi/2$ corresponds to $\theta_{\text{lab}}=\arctan(1/(\beta\gamma))$, where the ratio Eq. (\[eq:BT\_BL\_ratio\]) becomes largest when $\phi=(\phi_1+\phi_2)/2$. When $\sqrt{s} \gg m_e$, the ratio becomes largest in the very forward region, $(\theta_{\text{lab}}, \phi) \simeq (\arctan(2m_e/\sqrt{s}), (\phi_1+\phi_2)/2)$.

Azimuthal angle dependence
--------------------------

The transverse asymmetries in both Bhabha scattering and AIF have azimuthal dependence $\cos (2 \phi - \phi_1 - \phi_2)$ where $\phi$ is the scattering/emission azimuth, $\phi_1$ is the azimuth of positron polarization, and $\phi_2$ is the azimuth of target polarization. The azimuth $\phi_1$ is defined by the ratio between $P_{ \text T1}$ and $P_{ \text T2}$. The transverse polarizations $P_{ \text T1}$ and $P_{ \text T2}$ come from the $R(e)$ and $T(e)$ terms in Eq. (\[eq:drate\]), respectively. Only the latter violates the time-reversal symmetry.

In typical left-right symmetric models, the magnitude of the ratio $|T(e)/N(e)|$ is typically expected to be on the order of $(M_{\text{WL}}/M_{\text{WR}})^2$. When we assume the right-handed weak gauge boson mass $M_{\text{WR}}$ ranges from 1 TeV to 10 TeV, $|T(e)/N(e)| = \mathcal{O}(10^{-4}) - \mathcal{O}(10^{-2})$. On the other hand, $|R(e)/N(e)|$ is expected to be on the order of $m_e/m_\mu = \mathcal{O}(10^{-2})$. These ratios give expected ranges of $\phi_1$ from 45 degrees to 135 degrees or from -135 degrees to -45 degrees. To test the $T$-symmetry, the direction of target polarization will be in parallel or anti-parallel to the $y$-axis. In the setup of Figs. \[fig:exp\_pol\] and \[fig:exp\_scat\], it is most advantageous to make the direction of the target polarization be in parallel or antiparallel to the $y$-axis.

Possibility of a new experiment at J-PARC {#section:section5}
=========================================

In this section, a new measurement of the transverse polarization of positrons from muon decays with the high-intensity pulsed muon beam at J-PARC is proposed. At Materials and Life Science Experimental Facility (MLF) Muon Science Establishment (MUSE), the world-highest intensity of pulsed muon beam has been established. The high-intensity beam is beneficial not only for the statistical precision but also a shortening of measurement time and suppression of the systematic uncertainty arising from temporal variations of the environment.

Precursor experiments at PSI
----------------------------

The most recent experimental search for the transverse polarization of muon decay positrons was performed at Paul Scherrer Institute (PSI) [@Danneberg2005]. This precursor experiment employed the AIF polarimeter and the continuous muon beam which is delivered at random timing. As an annihilation target, a magnetized metal foil with a thickness of 1 mm was placed between the muon stopping target and the photon detector. A honeycomb array calorimeter using $\rm{Bi_{12} Ge O_{20}}$ (BGO) crystals was employed for measuring photons. The target polarization was $P_e = 7.2\ \%$ and the beam polarization was $P_\mu = 82\pm 2\ \%$. The number of analyzed annihilation events was 36.4 million. The results were $$\begin{aligned}
P_{\text T1}  &= (6.3 \pm 7.7 \pm 3.4) \times 10^{-3}, \\
P_{\text T2}  &= (-3.7 \pm 7.7 \pm 3.4) \times 10^{-3},\end{aligned}$$ where the first error is statistical uncertainty and the second is systematics. The systematic error was dominated by uncertainty in energy calibration for the calorimeter and positron energy loss before an annihilation which results in the fluctuation of the positron polarization. Note that another experiment for a measurement of the Michel parameter $\xi '' $ in polarized muon decay was performed with the apparatus aforementioned [@Prieels2014]. The result was $\xi''=0.981\pm 0.045_{\text {stat}}\pm 0.003_{\text{syst}}$. The parameter $\xi''$ is concerned with the positron’s longitudinal polarization $P_{\text L}$, whereas our concern in the proposed experiment is the positron’s transverse polarization $P_{\text T2}$.

New experiment at J-PARC
------------------------

The new experiment proposed in this paper comprises three major components: high-intensity pulsed muon beam at J-PARC, an electron-polarized scattering/annihilation target, and segmented positron and photon detectors.

### Pulsed muon beam

In contrast to the continuous beam at PSI, the pulsed muon beam at J-PARC has a periodic timing structure by nature. For a case of an experiment with pulsed beam, no muon trigger counter is required because arrival of the muon beam is synchronized to the accelerator repetition. The repetition cycle is 25 Hz and the beam has double pulse structure with the full-width at half-maximum of 100 ns. The two bunches have a timing interval of 600 ns. The expected muon beam intensity at J-PARC MLF MUSE H-Line is $1\times10^8\ \mu^+/\rm{s}$ [@Kawamura2018].

In the precursor experiments with the continuous beam, a simultaneous measurement of the Bhabha scattering and the AIF was impossible because of the difference in the data taking trigger configuration. Signatures of the Bhabha and AIF events are an electron-positron pair and photon pairs, respectively. In the new experiment, both the Bhabha and AIF measurements can be performed simultaneously.

### Decay positron polarimeter

![Experimental schematic of the decay positron polarimeter: (1) the beryllium plate for muon stopping; (2) the spin-analyzing target between the silicon pixel detectors; (3) the silicon strip detectors for electron/positron tracking; (4) the segmented BGO calorimeter for gamma-ray detection. The inset in the bottom displays a magnified figure of the spin-analyzing target. The lead collimator is placed between the muon-stopping target and the spin-analyzing target. The dashed arrows indicate the magnetization directions of the iron foils. Note that the diagram is not drawn to scale.[]{data-label="fig:exp_setup"}](exp_setup.eps){width="\linewidth"}

The polarimeter is designed for both the Bhabha and AIF measurements. Figure \[fig:exp\_setup\] depicts a schematic of the polarimeter. A thin plate of beryllium is placed for a muon-stopping target. Magnetized iron foils are employed as a spin-analyzing target. Between the muon-stopping target and the spin-analyzing target, a lead collimator is placed to limit the incident angle of positrons.

In front and behind of the spin analyzing target, silicon pixel detectors are placed in order to determine the interaction point and the positron incident angle. An electron pair from Bhabha scattering and a gamma-ray pair from AIF are detected by silicon strip detectors and segmented BGO calorimeters, respectively. The muon-stopping target is placed at the angle of 45 degrees with respect to the muon beam direction for suppression of Coulomb multiple scattering inside the target. The spin-analyzing target is placed at the angle of 60 degrees because the magnetization of the iron foil is along to an in-plane direction. In the experiment, decay positrons are polarized in anti-parallel to the $z$-axis. The $z$-projection of the target polarization is in parallel or anti-parallel to the positron polarization. We denote the former and the latter cases as “positive" and “negative" configurations, respectively.

The measurement system is axial symmetric with respect to the direction of the incident muon beam except for the tilts of targets and the plurality of polarimeters can be placed to multiply statistics. A set of eight unit polarimeters under the symmetric arrangement is practical. Symmetric placement of the polarimeters is beneficial not only for statistics but also an understanding of systematic uncertainties arising from the target polarization and detector performances.

Simulation study
----------------

Feasibility of the designed polarimeter is evaluated by Monte-Carlo simulation using GEANT4. Performance of polarimeters is quantified by the figure of merit which is defined by Eq. (\[eq:FOM\]). In order to obtain the FOM, the number of signal events and angular asymmetries are calculated stepwise.

### Angular distributions

The number of signal events per incident muon to the stopping target is calculated. For the case of Bhabha scattering, the signal events are selected by a pair of charged particle tracks, which is detected by the silicon tracker. For the case of AIF, the signals are discriminated by a pair of photons, which is detected by the calorimeters.

Simulated angular distributions are shown in Fig. \[fig:exp\_theta\]. The result is obtained with the unpolarized spin-analyzing target with the thickness of 1 mm. Here and hereafter we denote the scattering/emission angle in the lab frame $\theta_{\text{lab}}$ as $\theta$ for short. In the simulation, positive muons having the kinetic energy of 4 MeV irradiate the muon-stopping target. The azimuthal angular distribution can be assumed uniform in the region of interest.

![Simulated angular distributions of Bhabha scattering and AIF: (top) scattering/emission angle ($\theta$); (bottom) azimuthal angle $\phi$. The histograms in black and red lines correspond to Bhabha scattering and AIF, respectively. The ordinates are normalized by the number of positrons on the spin-analyzing target. The ends of the distributions result from the boundary of simulated geometry.[]{data-label="fig:exp_theta"}](exp_theta_phi.eps){width="0.9\linewidth"}

### Angular asymmetries

Spin-analyzing power of the polarimeter is numerically evaluated for polarized positrons. Figure \[fig:exp\_asym\] shows typical simulated asymmetries as functions of the scattering/emission angle $\theta$. The polarization of positrons from muon decay is assumed to be 100$\%$ in the longitudinal direction for calculation of the longitudinal asymmetry. The transverse asymmetry is calculated using the longitudinal asymmetry and the conversion factor aforementioned. Analyzing power for Bhabha scattering and AIF depend on the kinetic energies of electron-positron pair or two annihilation-quanta. Optima for the energy threshold for a measurement of the longitudinal polarization have been obtained to be 5.14 MeV and 6.66 MeV for the Bhabha polarimeter and AIF polarimeter, respectively [@Corriveau1981]. In our work, energy thresholds for a transverse polarization measurement are optimized with consideration of background processes.

![Simulated asymmetries with Bhabha scattering and AIF: (top) the result of Bhabha polarimeter case; (bottom) the result of AIF polarimeter case. The open and filled circles correspond to the longitudinal and transverse asymmetries, respectively.[]{data-label="fig:exp_asym"}](exp_asym.eps){width="0.9\linewidth"}

Maximum transverse asymmetries are 0.02 for Bhabha scattering and 0.3 for AIF. This difference mainly resulted from background contamination and the factor of seven difference in the denominators in Eq. (\[eq:AT\_AL\_ratio\]) and Eq. (\[eq:BT\_BL\_ratio\]). Detection of a pair of an electron and a positron is likely to be affected by background processes such as bremsstrahlung rather than a pair of annihilation photons. Details of event selection will be described in the next subsection. It is obvious that the AIF polarimeter is suitable for a measurement of the transverse polarization. However, the Bhabha polarimeter is still useful for measurements of the positron’s longitudinal polarization and the muon beam polarization. Simultaneous measurements of the polarizations are important to understand the systematic uncertainties arising from beam quality and detector response. Hereafter, we will focus on the AIF polarimeter.

### Sensitivity to the $T$-violating effect

In order to maximize the sensitivity for the transverse polarization, essential parameters in the experiment and analysis are optimized considering the FOM $$\mathcal{F} = \frac{(N_+ - N_-)^2}{N_+ + N_-},$$ where $N_+$ and $N_-$ are the numbers of observed events with “positive" and “negative" configurations in the target polarization, respectively. The thickness of the spin-analyzing target, the acceptances for the emission angle and the energy thresholds are optimized.

The number of signal events $N_+$ and $N_-$ are numerically calculated using the asymmetry $B_L$ shown in Fig. \[fig:exp\_asym\] and Eq. (\[eq:BT\_BL\_ratio\]). The thickness of the spin-analyzing target, the energy thresholds for calorimeters, and the acceptance of the emission angle are optimized to maximize the FOM. In the optimization procedure, the thresholds are set for both the energies of individual photons ($E_1$ and $E_2$) and the sum of the photon energies ($E_1+E_2$).

Figure \[fig:exp\_tgt\] shows the target thickness dependence of the FOM. The target thickness of 3 mm is the optimum. For this case, the optimized energy thresholds are $E_1 > 2$ MeV, $E_2 > 2$ MeV, and $(E_1+E_2)>$ 15 MeV. The acceptance for $\theta$ is set $13\pm8$ degrees. Practically, the spin-analyzing target with the thickness of 3 mm will be consisting of three layers of 1 mm thick foils.

![Calculated target thickness dependence of the FOM. At each target thickness, the energy thresholds and the emission angle acceptance are optimized to maximize the FOM. The number of simulated muons is $1\times 10^7$ for each point.[]{data-label="fig:exp_tgt"}](exp_tgt.eps){width="0.95\linewidth"}

As described in Section \[section:section4\], the transverse asymmetry depends on the angles $\phi$, $\phi_1$, and $\phi_2$. Among these angles, the target polarization azimuth $\phi_2$ can be chosen at our disposal and the emission azimuth $\phi$ is obtained by the segmented calorimeter. Hence, we can extract $\phi_1$ from the observed azimuthal dependence of the transverse asymmetry. Figure \[fig:exp\_cos\] shows an example of expected azimuthal dependence for a case of $P_{\text T2} = 7.7 \times 10^{-3}$ that corresponds to statistical uncertainty in the precursor experimental result [@Danneberg2005]. The target polarization $\phi_2$ is set zero and $\phi_1$ is assumed to be 45 degrees. The number of simulated muons is $3.3 \times 10^9$, which corresponds to ten days of data acquisition at J-PARC.

![Simulated transverse asymmetry as a function of the emission azimuth $\phi$. The red solid line indicates the fitting result by a cosine function. The oscillation amplitude obtained by the fitting is $(1.5 \pm 0.2 ) \times 10^{-4}$[]{data-label="fig:exp_cos"}](exp_cos_rev.eps){width="0.95\linewidth"}

We define the statistical significance of $P_{\text{T2}}$ by $$\sigma \equiv \frac{P_{\text{T2}}}{\delta P_{\text{T2}}} ,
\label{eq:sigma_def}$$ where $\delta P_{\text{T2}}$ is the statistical uncertainty of $P_{\text{T2}}$. Since $P_{\text{T2}}$ can be written as $$P_{\text{T2}} = c a \cos(\phi_1),$$ where $c$ is a constant and $a$ is the amplitude of the fitted cosine function in the $\phi_1$ vs observed transverse asymmetry plane like Fig. \[fig:exp\_cos\], Eq. (\[eq:sigma\_def\]) can be written as $$\sigma
 = \frac{1}{ \sqrt{ (\delta a/a)^2 + (\delta \phi_1 \tan \phi_1)^2}}.$$ In the analysis, $\sigma \ge 3$ is required for a statistically significant observation of the transverse polarization $P_{\text{T2}}$. The transverse polarization has two components ($P_{\text{T1}}$ and $P_{\text{T2}}$), and $\tan\phi_1$ gives the ratio of them. A deviation of $\phi_1$ from 90 degrees together with a non-vanishing amplitude $a$ means a non-zero $P_{\text{T2}}$, $i.e.$, a $T$-violating effect in muon decay.

Amount of measurement time necessary for $\sigma \ge 3$ is estimated for several values of the transverse positron polarization $P_{\text{T2}}$. Figure \[fig:exp\_sig\] shows the sensitivity as a function of measurement time with the AIF polarimeter for $\phi_1=45, 60$ and 75 degrees. When $|\phi_1-90^\circ|$ is large enough, the sensitivity is dominated by the requirement for the amplitude. On the other hand, the requirement for $\delta \phi_1$ is decisive for smaller $|\phi_1-90^\circ|$ cases. The muon beam intensity of $1\times10^8\ \mu^+/\rm{s}$ is assumed. The polarization of spin-analyzing target is set $7.2\%$. The target polarization direction is set “positive" or “negative" configurations and measurement time is equally divided into each configuration.

![Expected sensitivity to the transverse polarization using the AIF polarimeter. The abscissa represents the amount of time for $\sigma \ge 3$. The curves in each color correspond to the respective value of $\phi_1$. The horizontal dotted line indicates the statistical uncertainty for $P_{\text{T1}}$ and $P_{\text{T2}}$ in the precursor measurement [@Danneberg2005].[]{data-label="fig:exp_sig"}](exp_sig_rev.eps){width="0.95\linewidth"}

According to the estimation, when the magnitudes of $P_{\text{T1}}$ and $P_{\text{T2}}$ are similar, it is feasible to achieve a comparable statistical precision relative to that of the precursor experiment within four days.

Discussions {#section:section6}
===========

The proposed experiment is feasible at the existing muon facility. However, further improvement is possible with the future facility upgrade. At J-PARC, construction of a new target station is being considered. At the second target station, the muon beam intensity is expected to be several tens times higher than that at the present facility. The higher beam intensity contributes only to the statistical precision. The systematic uncertainty may be improved by using the magnetic resonance technique. In an analogy to the $\pi/2$-pulse method for nuclear magnetic resonance, the muon spin rotation will occur under a longitudinal static magnetic field and a pulse of radio-frequency field in the transverse direction [@Carne1984]. This technique enables us to rotate the muon spin polarization to an arbitrary direction. By combining the results of measurements with opposite rotation directions, the systematic uncertainty arising from the target polarization will be highly suppressed.

Lastly, we add some comments on an implication of $V+A$ interactions for neutrinoless double beta decay ($0\nu\beta\beta$) which is one of the most important experiments to detect the Majoranality of the neutrino [@Bilenky2010]. In the absense of $V+A$ interactions, it depends on the effective Majorana neutrino mass $$\begin{aligned}
{\langle}m_\nu{\rangle}\equiv |\sum U_{ej}^2m_j|\end{aligned}$$ in the $\nu$SM. In the presence of $V+A$ interactions, we have additionally $V+A$ effect and the inverse half-life time is given by [@Doi1985] $$\begin{aligned}
\left(T^{0\nu}_{1/2}\right)^{-1}
&=|M^{0\nu}|^2 \biggl[ C_1\left(\frac{{\langle}m_\nu{\rangle}}{m_e}\right)^2
+C_2{\langle}\lambda {\rangle}\frac{{\langle}m_\nu{\rangle}}{m_e}\cos \psi_1 \nonumber\\
&~~+C_3{\langle}\eta {\rangle}\frac{{\langle}m_\nu{\rangle}}{m_e}\cos\psi_2+C_4{\langle}\lambda{\rangle}^2+C_5 {\langle}\eta{\rangle}^2 \nonumber\\
&~~+C_6{\langle}\lambda{\rangle}{\langle}\eta{\rangle}\cos (\psi_1-\psi_2) \biggr].
\label{eq:tmpEq42}\end{aligned}$$ Here the explicit forms of $C_i$ and $M^{0\nu}$ are given in the reference [@Doi1985]. The parameters ${\langle}\lambda{\rangle}$ and ${\langle}\eta{\rangle}$ are defined by $$\begin{aligned}
{\langle}\lambda{\rangle}=\lambda|\sum U_{ej}V_{ej}(g_{\text V}'/g_{\text V})|,
~~
{\langle}\eta{\rangle}=\eta|\sum U_{ej}V_{ej}|,\end{aligned}$$ where the mixing matrices are defined in Eq. [(\[mixing\])]{}. The parameters $g_{\text V}$ and $g_{\text V}'$ are the coupling constants of $V-A$ and $V+A$ hadronic currents, respectively.

The present lower limit of $T^{0\nu}_{1/2}$ is [@Gando2016] $$T^{0\nu}_{1/2}>1.07\times 10^{26}\ \mbox{yr at 90\% C.L.}
\label{eq:halflife_limit}$$ This implies $${\langle}m_\nu{\rangle}<61 \sim 165 ~\mbox{meV}$$ for the case without $V+A$ interaction. In the presence of $V+A$, Eq. (\[eq:halflife\_limit\]) gives a rather complex constraint from Eq. (\[eq:tmpEq42\]) and depends on the details of nuclear structures, whereas our experiment is restricted to the leptonic world and the measurement of $T(e)$ gives a rather simple relation on Eq. (\[Te\]). Anyhow, $T$-violation is out of scope in $0\nu\beta\beta$ experiments.

Summary {#section:section7}
=======

We have discussed the measurement of the $T$-violation in muon decay. The most important theoretical character is that it almost immediately indicates the Majoranality of neutrinos. This is based on the heavy right-handed neutrino, which is the necessary and robust object in most of the BSM physics models and also is the cause of the BAU via leptogenesis. Thus the $T$-violation (discovery of a non-zero $P_{\text T2}$) leads to the definitive information for the Majoranality of the neutrinos and BAU, two major unsolved problems in particle physics. We can check the distinction from the presence of the Michel parameters by estimating the order of BAU and comparing with other Majoranality experiments, as well as theoretical consistencies in many reliable models.

The precursor experiment at PSI obtained the limits on the positron transverse polarizations with the sensitivity of $8 \times 10^{-3}$. We have proposed a new experiment with the high-intensity pulsed muon beam at J-PARC. The new experiment will employ the segmented annihilation-in-flight polarimeter with high-rate capability. The sensitivity and feasibility of the polarimeter have been studied numerically. We expect a factor of ten improvement in the sensitivity by a year of measurement.

We are greatly indebted to T. Mibe and H. Nishiura for useful discussions. This work is supported by JSPS KAKENHI Grants Nos.17H01133 (K.S, T.F, and D.N) and 16K05323 (D.N).

[^1]: The symmetry breaking to naive $\text{SU}(5)\otimes \text{U}(1)_{\text X}$ leads to too fast proton decay.

[^2]: By using the same technique used in Ref. [@Shabalin1978rs], it can be shown that no $CP$ violation occurs if we neglect the gluon in Fig. \[fig:CPVSM\] (b).
