Predictions for the neutrino parameters in the minimal gauged U(1)$_{L_\mu-L_\tau}$ model

TL;DR

This paper analyzes the neutrino mass matrix structure in a minimal gauged U(1)$_{L___ au}$ model, predicting specific CP phases, neutrino masses, and neutrinoless double beta decay parameters consistent with current experiments.

Contribution

It provides a unique, symmetry-based prediction of neutrino parameters within the minimal gauged U(1)$_{L___ au}$ model, linking mass matrix structure to observable quantities.

Findings

Predicts Dirac CP phase around 1.59π-1.70π.

Estimates sum of neutrino masses between 0.12-0.40 eV.

Forecasts neutrinoless double beta decay effective mass 0.017-0.12 eV.

Abstract

We study the structure of the neutrino mass matrix in the minimal gauged U(1)$_{L_\mu-L_\tau}$ model, where three right-handed neutrinos are added to the Standard Model in order to obtain non-zero masses for active neutrinos. Because of the U(1)$_{L_\mu-L_\tau}$ gauge symmetry, the structure of both Dirac and Majorana mass terms of neutrinos is tightly restricted. In particular, the inverse of the neutrino mass matrix has zeros in the $(\mu,\mu)$ and $(\tau,\tau)$ components, namely, this model offers a symmetric realization of the so-called two-zero-minor structure in the neutrino mass matrix. Due to these constraints, all the CP phases-the Dirac CP phase $\delta$ and the Majorana CP phases $\alpha_2$ and $\alpha_3$-as well as the mass eigenvalues of the light neutrinos $m_i$ are uniquely determined as functions of the neutrino mixing angles $\theta_{12}$, $\theta_{23}$, and $\theta_{13}$, and the squared mass differences $\Delta m_{21}^2$ and $\Delta m_{31}^2$. We find that this model predicts the Dirac CP phase $\delta$ to be $\delta\simeq 1.59\pi$-$1.70\pi$ ($1.54\pi$-$1.78\pi$), the sum of the neutrino masses to be $\sum_{i}m_i \simeq 0.14$-0.22 eV ($0.12$-0.40 eV), and the effective mass for the neutrinoless double beta decay to be $\langle m_{\beta \beta }\rangle \simeq 0.024$-0.055 eV ($0.017$-0.12 eV) at $1\sigma$ ($2\sigma$) level, which are totally consistent with the current experimental limits. These predictions can soon be tested in future neutrino experiments. Implications for leptogenesis are also discussed.

Figures

Tables

Table 1. Table 1 : Input values for the neutrino oscillation parameters we use in this paper. We take them from Ref.  [ 51 ] .

Parameter	Best fit	1$\sigma$ range	2$\sigma$ range
$\delta m^2/{10}^{-5}{\text{eV}}^2$	7.37	7.21–7.54	7.07–7.73
$\Delta m^2/{10}^{-3}{\text{eV}}^2$	2.525	2.495–2.567	2.454–2.606
${\sin }^2{\theta }_{12}/{10}^{-1}$	2.97	2.81–3.14	2.65–3.34
${\sin }^2{\theta }_{23}/{10}^{-1}$	4.25	4.10–4.46	3.95–4.70
${\sin }^2{\theta }_{13}/{10}^{-2}$	2.15	2.08–2.22	1.99–2.31
$\delta /\pi$	1.38	1.18–1.61	1.00–1.90

Equations

$$Q_{L_{\mu }-L_{\tau }}(\text{Dirac}):\left(\begin{array}{ccc}0& 1& -1\\ -1& 0& -2\\ 1& 2& 0\end{array}\right),$$

$$Q_{L_{\mu }-L_{\tau }}(\text{Majorana}):\left(\begin{array}{ccc}0& 1& -1\\ 1& 2& 0\\ -1& 0& -2\end{array}\right),$$

$$\Delta \mathcal{L}=$$

$$-\frac{1}{2}{M}_{ee}{N}_e^c{N}_e^c-M_{\mu \tau }{N}_{\mu }^c{N}_{\tau }^c-{\lambda }_{e\mu }\sigma N_e^c{N}_{\mu }^c-{\lambda }_{e\tau }{\sigma }^{\ast }{N}_e^c{N}_{\tau }^c+\text{h.c.},$$

$${\mathcal{L}}_{\mathrm{mass}}$$

$${\mathcal{M}}_D=\frac{v}{\sqrt{2}}\left(\begin{array}{ccc}{\lambda }_e& 0& 0\\ 0& {\lambda }_{\mu }& 0\\ 0& 0& {\lambda }_{\tau }\end{array}\right),{\mathcal{M}}_R=\left(\begin{array}{ccc}{M}_{ee}& {\lambda }_{e\mu }⟨\sigma ⟩& {\lambda }_{e\tau }⟨\sigma ⟩\\ {\lambda }_{e\mu }⟨\sigma ⟩& 0& M_{\mu \tau }\\ {\lambda }_{e\tau }⟨\sigma ⟩& M_{\mu \tau }& 0\end{array}\right).$$

$${\mathcal{M}}_{{\nu }_L}\simeq -{\mathcal{M}}_D{\mathcal{M}}_R^{-1}{\mathcal{M}}_D^T.$$

$$U^T{\mathcal{M}}_{{\nu }_L}U=\text{diag}(m_1,m_2,m_3),$$

$$U=\left(\begin{array}{ccc}{c}_{12}{c}_{13}& s_{12}{c}_{13}& s_{13}{e}^{-i\delta }\\ -s_{12}{c}_{23}-c_{12}{s}_{23}{s}_{13}{e}^{i\delta }& c_{12}{c}_{23}-s_{12}{s}_{23}{s}_{13}{e}^{i\delta }& s_{23}{c}_{13}\\ s_{12}{s}_{23}-c_{12}{c}_{23}{s}_{13}{e}^{i\delta }& -c_{12}{s}_{23}-s_{12}{c}_{23}{s}_{13}{e}^{i\delta }& c_{23}{c}_{13}\end{array}\right)\left(\begin{array}{ccc}1& & \\ & e^{i\frac{{\alpha }_2}{2}}& \\ & & e^{i\frac{{\alpha }_3}{2}}\end{array}\right),$$

$${\mathcal{M}}_{{\nu }_L}^{-1}=U\text{diag}(m_1^{-1},m_2^{-1},m_3^{-1})U^T\simeq -({\mathcal{M}}_D^{-1}{)}^T{\mathcal{M}}_R{\mathcal{M}}_D^{-1}.$$

$$\frac{1}{m_1}{V}_{\mu 1}^2+\frac{1}{m_2}{V}_{\mu 2}^2{e}^{i{\alpha }_2}+\frac{1}{m_3}{V}_{\mu 3}^2{e}^{i{\alpha }_3}$$

$$\frac{1}{m_1}{V}_{\tau 1}^2+\frac{1}{m_2}{V}_{\tau 2}^2{e}^{i{\alpha }_2}+\frac{1}{m_3}{V}_{\tau 3}^2{e}^{i{\alpha }_3}$$

$$e^{i{\alpha }_2}$$

$$R_2$$

$$R_3$$

$$\frac{m_2}{m_1}=\frac{1}{|R_2(\delta )|},\frac{m_3}{m_1}=\frac{1}{|R_3(\delta )|}.$$

$$\delta m^2=\Delta m_{21}^2,\Delta m^2=\Delta m_{31}^2-\frac{1}{2}\Delta m_{21}^2.$$

$$\delta m^2$$

$$\Delta m^2$$

$$\delta m^2$$

$$\Delta m^2+\frac{\delta m^2}{2}$$

$$\cos \delta \simeq \frac{\cot 2{\theta }_{12}\cot 2{\theta }_{23}}{\sin {\theta }_{13}}.$$

$${\mathcal{L}}_{\mathrm{eff}}=\frac{1}{2}{C}_{\alpha \beta }(L_{\alpha }\cdot H)(L_{\beta }\cdot H)+\text{h.c.},$$

$$\mu \frac{dC}{d\mu }=-\frac{3}{32{\pi }^2}\left[{\left(Y_e^{†}{Y}_e\right)}^{T}C+C\left(Y_e^{†}{Y}_e\right)\right]+\frac{K}{16{\pi }^2}C,$$

$$K=-3g_2^2+2\mathrm{Tr}\left(3Y_u^{†}{Y}_u+3Y_d^{†}{Y}_d+Y_e^{†}{Y}_e\right)+2\lambda ,$$

$$C(t)=I_K(t)\mathcal{I}(t)C(0)\mathcal{I}(t),$$

$$I_K(t)$$

$${\theta }_{12},{\theta }_{23},{\theta }_{13},\delta m^2,\Delta m^2,$$

$$⟨m_{\beta \beta }⟩\equiv \left|\sum _i{U}_{ei}^2{m}_i\right|=\left|c_{12}^2{c}_{13}^2{m}_1+s_{12}^2{c}_{13}^2{e}^{i{\alpha }_2}{m}_2+s_{13}^2{e}^{i({\alpha }_3-2\delta )}{m}_3\right|.$$

$$\Delta \mathcal{L}=$$

$${\mathcal{M}}_R$$

$${\hat{N}}_i^c$$

$${\hat{\lambda }}_{i\alpha }$$

$${ϵ}_1$$

Full text

UT-17-17

IPMU-17-0072

**Predictions for the neutrino parameters

in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model

**

Kento Asai1, Koichi Hamaguchi1,2, Natsumi Nagata1

1Department of Physics, University of Tokyo, Bunkyo-ku, Tokyo 133–0033, Japan

2Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), University of Tokyo, Kashiwa 277–8583, Japan

We study the structure of the neutrino mass matrix in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model, where three right-handed neutrinos are added to the Standard Model in order to obtain non-zero masses for active neutrinos. Because of the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry, the structure of both Dirac and Majorana mass terms of neutrinos is tightly restricted. In particular, the inverse of the neutrino mass matrix has zeros in the $(\mu,\mu)$ and $(\tau,\tau)$ components, namely, this model offers a symmetric realization of the so-called two-zero-minor structure in the neutrino mass matrix. Due to these constraints, all the CP phases—the Dirac CP phase $\delta$ and the Majorana CP phases $\alpha_{2}$ and $\alpha_{3}$—as well as the mass eigenvalues of the light neutrinos $m_{i}$ are uniquely determined as functions of the neutrino mixing angles $\theta_{12}$, $\theta_{23}$, and $\theta_{13}$, and the squared mass differences $\Delta m_{21}^{2}$ and $\Delta m_{31}^{2}$. We find that this model predicts the Dirac CP phase $\delta$ to be $\delta\simeq 1.59\pi$–$1.70\pi$ ($1.54\pi$–$1.78\pi$), the sum of the neutrino masses to be $\sum_{i}m_{i}\simeq 0.14$–0.22 eV ($0.12$–0.40 eV), and the effective mass for the neutrinoless double beta decay to be $\langle m_{\beta\beta}\rangle\simeq 0.024$–0.055 eV ($0.017$–0.12 eV) at $1\sigma$ ($2\sigma$) level, which are totally consistent with the current experimental limits. These predictions can soon be tested in future neutrino experiments. Implications for leptogenesis are also discussed.

1 Introduction

A gauged U(1)${}_{L_{\mu}-L_{\tau}}$ symmetric extension of the Standard Model (SM), where $L_{\mu}$ and $L_{\tau}$ stand for the $\mu$ and $\tau$ numbers, respectively, has widely been considered so far, as it is one of the several possibilities of gauging an accidental U(1) symmetry in the SM [1, 2, 3, 4]. An attractive feature of this class of models is that the muon $g-2$ anomaly [5, 6, 7, 8] may be explained by the loop contribution of the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge boson if its mass lies around the weak scale or lower [9, 10, 11, 12], though this possibility is severely restricted by the searches of the neutrino trident production process [13, 14, 15, 16]. These models may also explain anomalies in flavor physics [15, 17] and offer promising candidates for dark matter in the Universe [18, 19, 20, 21, 22]. For other recent studies on gauged U(1)${}_{L_{\mu}-L_{\tau}}$ models, see Refs. [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35].

On the other hand, neutrino oscillation data show that at least two of the active neutrinos have non-zero masses and there is sizable mixing among these neutrinos. This can be accounted for if we add right-handed neutrinos to the theory and couple them to the SM lepton and Higgs fields through the Yukawa couplings. After the Higgs field acquires a vacuum expectation value (VEV), these terms lead to the Dirac mass terms for the neutrinos. In addition, the right-handed neutrinos can have Majorana mass terms, and if these masses are much larger than the electroweak scale, then the smallness of the neutrino masses can naturally be explained by the seesaw mechanism [36, 37, 38, 39].

In gauged U(1)${}_{L_{\mu}-L_{\tau}}$ models, however, the neutrino mass structure is tightly restricted since the second and third generation leptons are charged under the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry [40, 41, 42, 43, 44]. In fact, in the U(1)${}_{L_{\mu}-L_{\tau}}$ preserving limit, the Dirac mass matrix is diagonal, while in the Majorana mass matrix only the $(e,e)$, $(\mu,\tau)$, and $(\tau,\mu)$ components can be non-zero. Such a simple structure cannot explain neutrino oscillation data, and therefore we need to break the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry. To that end, one usually introduces a SM singlet scalar that has a non-zero U(1)${}_{L_{\mu}-L_{\tau}}$ charge such that its VEV spontaneously breaks the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry and gives a mass to the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge boson. We however note that even though the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry is spontaneously broken, the structure of the neutrino mass matrix is still highly constrained if we introduce just one U(1)${}_{L_{\mu}-L_{\tau}}$-breaking scalar field and consider only renormalizable interactions. Therefore, it is interesting to study if this minimal setup can accommodate the neutrino mass structure that is consistent with the present neutrino oscillation data.

In this paper, we consider this “minimal” gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model where three right-handed neutrinos and one U(1)${}_{L_{\mu}-L_{\tau}}$-charged SM singlet scalar field are added to the SM. Then, it turns out that the observed neutrino mixing structure can be obtained only when the U(1)${}_{L_{\mu}-L_{\tau}}$-breaking scalar field has the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $\pm 1$. In this case, the $(\mu,\mu)$ and $(\tau,\tau)$ components of the Majorana mass matrix for the right-handed neutrinos remain zero even after the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry is spontaneously broken. Because of this structure of the Majorana mass matrix together with the diagonal Dirac mass matrix, the inverse of the neutrino mass matrix also has zeros in the $(\mu,\mu)$ and $(\tau,\tau)$ components [42, 45]. The minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model thus gives a concrete realization of a two-zero-minor model [46, 47]. Intriguingly, due to the condition that the $(\mu,\mu)$ and $(\tau,\tau)$ components in the inverse of the neutrino mass matrix vanish, all the CP phases in the neutrino mixing matrix—the Dirac CP phase $\delta$ and the Majorana CP phases $\alpha_{2}$ and $\alpha_{3}$—as well as the mass eigenvalues of the light neutrinos are uniquely determined as functions of the neutrino mixing angles $\theta_{12}$, $\theta_{23}$, and $\theta_{13}$, and the squared mass differences $\Delta m_{21}^{2}$ and $\Delta m_{31}^{2}$. As we see below, this prediction is independent of the U(1)${}_{L_{\mu}-L_{\tau}}$-symmetry breaking scale, and thus can be regarded as a generic prediction in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model. We find that the predicted values of the neutrino parameters are consistent with the present neutrino data and can be tested in future neutrino experiments. We also discuss the implications of our results for leptogenesis.

This paper is organized as follows. In Sec. 2, we introduce the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model and examine the neutrino mass structure in this model. We then show in Sec. 3 the predicted values of the Dirac CP phase $\delta$, the sum of the neutrino masses $\sum_{i}m_{i}$, and the effective mass for the neutrinoless double beta decay $\langle m_{\beta\beta}\rangle$ using the neutrino mixing angles and the squared mass differences obtained in neutrino oscillation experiments. In Sec. 4, we discuss the implications for the leptogenesis. Finally, our conclusions are summarized in Sec. 5. This paper ends with two Appendices which give further details of our analysis. In Appendix A, we present some miscellaneous formulae which are useful to study the neutrino structure in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model. In Appendix B, we perform similar analyses for the minimal gauged U(1)${}_{L_{e}-L_{\mu}}$ and U(1)${}_{L_{e}-L_{\tau}}$ models, and show that these minimal models fail to explain the observed neutrino oscillation data.

2 Neutrino mass structure in the minimal gauged

U(1)${}_{L_{\mu}-L_{\tau}}$ model

To begin with, we describe the minimal U(1)${}_{L_{\mu}-L_{\tau}}$ model which we discuss in this paper. The model possesses a new U(1) gauge symmetry U(1)${}_{L_{\mu}-L_{\tau}}$. Under this gauge symmetry, $\mu_{L,R}$ and $\nu_{\mu}$ have the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $+1$, $\tau_{L,R}$ and $\nu_{\tau}$ have the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $-1$, and the other SM fields have the zero U(1)${}_{L_{\mu}-L_{\tau}}$ charge. We also introduce three right-handed neutrinos $N_{e}$, $N_{\mu}$, and $N_{\tau}$ to obtain non-zero neutrino masses. The U(1)${}_{L_{\mu}-L_{\tau}}$ charges of these fields are [math], $+1$, and $-1$, respectively.111Generically speaking, the U(1)${}_{L_{\mu}-L_{\tau}}$ charges of the right-handed neutrinos can be $(0,a,-a)$ without spoiling anomaly cancellation conditions. We however find that only the $|a|=1$ case gives a neutrino mass structure that can explain the neutrino oscillations. Then we can define the right-handed neutrinos with the U(1)${}_{L_{\mu}-L_{\tau}}$ charge 0, $+1$, and $-1$, as $N_{e}$, $N_{\mu}$, and $N_{\tau}$, respectively, without loss of generality. We also note that the introduction of two right-handed neutrinos is insufficient since the observed neutrino mixing angles cannot be reproduced in this case.

In the $(e,\mu,\tau)$ basis, the U(1)${}_{L_{\mu}-L_{\tau}}$ charges of the Dirac Yukawa terms are

[TABLE]

where the $(\alpha,\beta)$ entry in the above matrix represents the U(1)${}_{L_{\mu}-L_{\tau}}$ charge of the fermion bilinear term $N_{\alpha}^{c}L_{\beta}$, with $\alpha,\beta$ the flavor indices. For the Majorana mass terms of the right-handed neutrinos, on the other hand, we have

[TABLE]

where the $(\alpha,\beta)$ component indicates the U(1)${}_{L_{\mu}-L_{\tau}}$ charge of the fermion bilinear term $N_{\alpha}N_{\beta}$. From Eq. (1), we find that the Dirac Yukawa matrix is always diagonal in the gauged U(1)${}_{L_{\mu}-L_{\tau}}$ models—for the same reason, the charged lepton Yukawa matrix is also diagonal. As long as renormalizable interactions are considered, this structure is not violated even if we introduce a U(1)${}_{L_{\mu}-L_{\tau}}$-breaking scalar field. On the other hand, the charges of the Majorana mass matrix in Eq. (2) show that only the $(e,e)$, $(\mu,\tau)$, and $(\tau,\mu)$ components can have non-zero values in the U(1)${}_{L_{\mu}-L_{\tau}}$-symmetric limit. As we mentioned above, with this simple structure, we cannot explain the required values of the neutrino mixing angles. We therefore introduce a scalar boson $\sigma$ which has a non-zero U(1)${}_{L_{\mu}-L_{\tau}}$ charge, and couple it to right-handed neutrinos. After this scalar field develops a VEV, these couplings lead to Majorana mass terms of the right-handed neutrinos. If the scalar field has the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $\pm 1$, then the $(e,\mu)$, $(e,\tau)$, $(\mu,e)$, and $(\tau,e)$ components in Eq. (2) can be induced after the scalar field acquires a VEV, while the $(\mu,\mu)$ and $(\tau,\tau)$ can be generated if the scalar has the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $\pm 2$. In the latter case, however, the Majorana mass matrix becomes block-diagonal, which makes it unable to explain the observed neutrino mixing angles. We are thus left with the case where the scalar field has the U(1)${}_{L_{\mu}-L_{\tau}}$ charge $\pm 1$, and we take it to be $+1$ in the following discussion. We refer to this model as the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model.

The interaction terms relevant to neutrino masses are then given by

[TABLE]

where the dots indicate the contraction of the SU(2)L indices. After the Higgs field $H$ and the singlet scalar $\sigma$ acquire VEVs $\langle H\rangle=v/\sqrt{2}$ and $\langle\sigma\rangle$,222We can always take the VEV of $\sigma$ to be real by using U(1)${}_{L_{\mu}-L_{\tau}}$ transformations. respectively, these interaction terms lead to the neutrino mass terms:

[TABLE]

where

[TABLE]

The mass matrix for the light neutrinos is then given by [36, 37, 38, 39]

[TABLE]

An explicit expression for ${\cal M}_{\nu_{L}}$ can be found in Appendix A.1. We can diagonalize this mass matrix by using a unitary matrix $U$ (PMNS matrix [48]):

[TABLE]

which can be parametrized as

[TABLE]

where $c_{ij}\equiv\cos\theta_{ij}$ and $s_{ij}\equiv\sin\theta_{ij}$ for $\theta_{ij}=[0,\pi/2]$, $\delta=[0,2\pi]$, and we have ordered $m_{1}<m_{2}$ without loss of generality. We follow the convention of the Particle Data Group [48], where $m_{2}^{2}-m_{1}^{2}\ll|m_{3}^{2}-m_{1}^{2}|$ and $m_{1}<m_{2}<m_{3}$ (Normal Ordering, NO) or $m_{3}<m_{1}<m_{2}$ (Inverted Ordering, IO).

If $m_{i}=0$ ($i=1$ or 3), then $\text{det}({\cal M}_{\nu_{L}})=0$. As shown in Appendix A.1, in this case we cannot have desired mixing angles since ${\cal M}_{\nu_{L}}$ becomes block-diagonal. Thus, we focus on the $m_{i}\neq 0$ case, where we obtain from Eqs. (6) and (7)

[TABLE]

We then notice that the $(\mu,\mu)$ and $(\tau,\tau)$ components of these terms vanish since ${\cal M}_{D}$ is diagonal and ${\cal M}_{R}$ has zeros in these components. This structure is sometimes called two-zero minor [46, 47]. For other previous studies of the two-zero minor structure, see Refs. [49, 50]. These two vanishing conditions then lead to

[TABLE]

where the matrix $V$ is defined by $U=V\cdot\text{diag}(1,e^{i\alpha_{2}/2},e^{i\alpha_{3}/2})$. Notice that neither the U(1)${}_{L_{\mu}-L_{\tau}}$-breaking singlet VEV $\langle\sigma\rangle$ nor Majorana masses $M_{ee}$ and $M_{\mu\tau}$ appear in these conditions explicitly. For this reason, the following discussions based on these equations have little dependence on the U(1)${}_{L_{\mu}-L_{\tau}}$-symmetry breaking scale; it may be around the electroweak scale, or as large as $10^{(13-15)}$ GeV, which is a prime scale for the masses of right-handed neutrinos since small neutrino masses are explained with ${\cal O}(1)$ Yukawa couplings via the seesaw mechanism [36, 37, 38, 39]. It follows from Eqs. (10) and (11) that

[TABLE]

with333These expressions are consistent with the results presented in Ref. [50]. We also find from the explicit expressions (A.3) and (A.4) that the corresponding equations in Ref. [49] disagree with ours.

[TABLE]

where we have used $\widetilde{V}^{T}=V^{-1}$ and ${\rm det}V=1$ with $\widetilde{V}$ being the cofactor matrix of $V$.444The cofactor $\tilde{A}_{ij}$ of a matrix $A$ is given by the determinant of the submatrix formed by removing the $i$-th row and $j$-th column of the matrix $A$, multiplied by a factor of $(-1)^{i+j}$. The cofactor matrix of $A$, $\widetilde{A}$, is defined by $\widetilde{A}\equiv(\tilde{A}_{ij})$. We then have ${A}^{-1}=({\rm det}A)^{-1}\widetilde{A}^{T}$. In Appendix A.2, we give explicit expressions for $R_{2}$ and $R_{3}$ in terms of neutrino oscillation parameters. By taking the absolute values of the equations in (12), we find

[TABLE]

Therefore, these mass ratios are given as functions of the Dirac CP phase $\delta$. Notice that $R_{2,3}^{*}(-\delta)=R_{2,3}$ since $R_{2,3}$ contains a single CP phase $\delta$. As a consequence, $m_{2,3}/m_{1}$ are symmetric under the reflection $\delta\to-\delta$ (and $\pi+\delta\to\pi-\delta$) as we see below.

In Fig. 1, we plot the mass ratios $m_{2}/m_{1}$ and $m_{3}/m_{1}$ against the Dirac CP phase $\delta$ using Eq. (15). The bands show uncertainty coming from the 1$\sigma$ error in the neutrino mixing parameters. For input parameters of the neutrino mixing angles, we use $\sin^{2}\theta_{12}=2.97^{+0.17}_{-0.16}\times 10^{-1}$, $\sin^{2}\theta_{13}=2.15^{+0.07}_{-0.07}\times 10^{-2}$, and $\sin^{2}\theta_{23}=4.25^{+0.21}_{-0.15}\times 10^{-1}$ [51] (see Table 1).555See also Refs.[52, 53]. It is found that the resultant uncertainty mainly comes from the error in $\theta_{23}$. From this figure, we find that when $\delta\simeq 0,\,2\pi$ the observed neutrino mixing angles are incompatible with the condition $m_{1}<m_{2}$, and thus these regions are excluded. Moreover, around $\delta\simeq\pi$ we have $m_{1}<m_{3}<m_{2}$, which disagrees with the possible neutrino mass ordering: either $m_{1}<m_{2}<m_{3}$ or $m_{3}<m_{1}<m_{2}$ [48]. As a consequence, in the region where a consistent neutrino mass ordering is obtained, the neutrino mass ordering is always Quasi-Degenerate Normal Ordering with $m_{1}\lesssim m_{2}\lesssim m_{3}$, which is realized around $\delta\sim\pi/2$ and $3\pi/2$.

According to the recent global-fit analysis in Ref. [51], the NO is somewhat favored over IO at $\sim 2\sigma$ level. There are quite a few proposed experiments that may determine the neutrino mass ordering at more than 3$\sigma$ level within a decade [54, 55], such as PINGU [56], ORCA [57], and JUNO [58, 59]. To confirm NO in these future experiments would be a first consistency check of our model.

Now that the mass ordering has been fixed, we determine $m_{1}$ and $\delta$ by using the above results. From the neutrino oscillation experiments, we can measure666 In terms of the squared mass differences $\Delta m_{21}^{2}\equiv m_{2}^{2}-m_{1}^{2}$ and $\Delta m_{31}^{2}\equiv m_{3}^{2}-m_{1}^{2}$, $\delta m^{2}$ and $\Delta m^{2}$ are expressed as

$\delta m^{2}=\Delta m_{21}^{2}~{},\qquad\Delta m^{2}=\Delta m_{31}^{2}-\frac{1}{2}\Delta m_{21}^{2}~{}.$

(16)

[TABLE]

These quantities are related to the neutrino mass ratios as

[TABLE]

By solving these equations, we can determine $m_{1}$ and $\delta$. The observed values of $\delta m^{2}$ and $\Delta m^{2}$ are given by $\delta m^{2}\simeq 7.37\times 10^{-5}~{}\text{eV}^{2}$ and $\Delta m^{2}\simeq 2.525\times 10^{-3}~{}\text{eV}^{2}$, respectively [51]. This means that the right-hand side of Eq. (19) is much smaller than that in Eq. (20). From Fig. 1, we see that such a hierarchy can be realized only when $m_{2}^{2}/m_{1}^{2}\simeq 1$, i.e., $|R_{2}(\delta)|\simeq 1$. With the explicit formula of $R_{2}(\delta)$ in Eq. (A.3), this leads to

[TABLE]

When the best-fit values of the mixing angles $\theta_{ij}$ are used, this leads to $\cos\delta\simeq 0.46$, which corresponds to $\delta\simeq 0.35\pi$ or $1.65\pi$. In Eq. (A.6) in Appendix A.3, we give a cubic equation whose solution gives an exact value of $\cos\delta$ as a function of the mixing angles $\theta_{ij}$ and the squared mass differences $\delta m^{2}$ and $\Delta m^{2}$. As discussed there, the solution (21) approximates the real solution of the cubic equation (A.6) at ${\cal O}(\delta m^{2}/\Delta m^{2})$ level. By solving Eq. (A.6) numerically, we find $\cos\delta=0.445$ for the best-fit values of $\theta_{ij}$, $\delta m^{2}$, and $\Delta m^{2}$, which means $\delta=0.353\pi$ or $1.647\pi$, and justifies the expected accuracy of the approximated formula (21).777This result disagrees with the observation in Ref. [21], where it was observed that $\delta\simeq 0$ is predicted in the minimal gauged $L_{\mu}-L_{\tau}$ model. In the analysis of Ref. [21], the Dirac Yukawa terms of neutrinos are supposed to be real in the basis where ${\cal M}_{R}$ has only one CP phase, though it is not the generic case. In addition, we find disagreement between the neutrino mass matrix shown in Ref. [21] and ours, as pointed out in Appendix A.1.

From Eqs. (13) and (19), we obtain $m_{1}$ as a function of $\delta$. An explicit expression for $m_{1}$ is given in Appendix A.2. As noted above, $|R_{2,3}(\delta)|$ are symmetric with respect to the reflection $\delta\to-\delta$. This implies that $m_{1}$ is also symmetric under this reflection and thus depends only on $\cos\delta$ (not $\sin\delta$). As a result, even though there are two solutions for $\delta$, $m_{1}$ (and thus $m_{2,3}$ as well) is uniquely determined. Finally, by substituting the above results into Eq. (12), we determine the Majorana CP phases $\alpha_{2}$ and $\alpha_{3}$. Again, since $R^{*}_{2,3}(-\delta)=R_{2,3}(\delta)$, we have $\alpha_{2,3}(-\delta)=-\alpha_{2,3}(\delta)$, as seen in Fig. 7 in Appendix A.2.

In the next section, we discuss the predictions of the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model with the recent oscillation data. Before closing this section, let us give some general remarks.

• •

If the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry breaking scale is much higher than the electroweak scale, we expect sizable quantum corrections to the neutrino mass matrix. Such quantum corrections can be taken into account by using renormalization group equations. Remarkably, it is found that the two-zero minor structure of ${\cal M}_{\nu_{L}}$ is preserved throughout the renormalization group flow [46, 47]. To see this, we first note that below the right-handed neutrino mass scale, which is around the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry breaking scale in this model, right-handed neutrinos are integrated out to give the following dimension-five effective operator:

[TABLE]

where $C_{\alpha\beta}$ has the two-zero minor structure at the right-handed neutrino mass scale. The renormalization group equation of the Wilson coefficient $C_{\alpha\beta}$ at one-loop level is [60]

[TABLE]

with

[TABLE]

where $Y_{u}$, $Y_{d}$, and $Y_{e}$ denote the up-type, down-type, and charged-lepton Yukawa matrices, respectively, $g_{2}$ is the SU(2)L gauge coupling, and $\lambda$ is the Higgs quartic coupling: ${\cal L}_{\rm quart}=-\frac{1}{2}\lambda(H^{\dagger}H)^{2}$. Now recall that the charged lepton Yukawa matrix is diagonal in our model. In this case, the above equation can readily be solved as follows [61]:

[TABLE]

where $t\equiv\ln(\mu/\mu_{0})$ with $\mu_{0}$ being the initial scale, and

[TABLE]

Note that ${\cal I}(t)$ is a diagonal matrix. Therefore, if $C_{\mu\mu}^{-1}(0)=C^{-1}_{\tau\tau}(0)=0$, then $C_{\mu\mu}^{-1}(t)=C^{-1}_{\tau\tau}(t)=0$, which proves that the two-zero minor structure of the Wilson coefficient $C$ remains at low energies. As a result, the two-zero minor neutrino-mass structure in our model is robust against quantum corrections, even if the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry breaking scale is much higher than the electroweak scale.

• •

By performing a similar analysis, we can study the minimal gauged U(1)${}_{L_{e}-L_{\mu}}$ and U(1)${}_{L_{e}-L_{\tau}}$ models. We again obtain two conditions similar to Eqs. (10) and (11), which follow from the zero components in the inverse of the neutrino mass matrix in each model. We however find that in these cases it is unable to find parameter space consistent with the observed values of the neutrino oscillation parameters, as demonstrated in Appendix B. We thus conclude that within the minimal approach discussed in this paper, the U(1)${}_{L_{\mu}-L_{\tau}}$ case is the only possibility that offers a desirable neutrino mass structure.

• •

If we go beyond the minimal model, we may obtain a different neutrino-mass structure in the presence of the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry. For example, Refs. [25, 33] discuss gauged U(1)${}_{L_{\mu}-L_{\tau}}$ models where neutrino masses are radiatively induced at one-loop level. It turns out that the neutrino mass matrices in these models have the form of the two-zero texture [62, 63, 64, 65, 66], which predicts IO for neutrino mass spectrum. Thus, the determination of neutrino mass hierarchy in future neutrino experiments allows us to distinguish this class of models from the minimal model, where NO is predicted as we have seen above.

3 Predictions for the neutrino parameters

Using the results obtained above, we now compute quantities relevant to neutrino experiments with the errors in the neutrino oscillation parameters taken into account. For input values, we use the values given in Ref. [51], which are summarized in Table 1. In particular, we take the three mixing angles and the two mass squared differences ,

[TABLE]

as input parameters, and evaluate the predicted values of the other parameters, including Dirac CP phase $\delta$, the absolute masses $m_{i}$, their sum $\sum_{i}m_{i}$, and the effective Majorana neutrino mass $\langle m_{\beta\beta}\rangle$. The prediction for the Majorana phases $\alpha_{2}$ and $\alpha_{3}$ is also presented in Appendix A.2.

In Fig. 2, we plot the Dirac CP phase $\delta$ as functions of $\theta_{23}$ in the red lines. We vary $\theta_{23}$ in the 2$\sigma$ renege, where the 1$\sigma$ range is in between the vertical thin dotted lines. The dark (light) red bands show the uncertainty coming from the 1$\sigma$ ($2\sigma$) errors in the other parameters $\theta_{12}$, $\theta_{13}$, $\delta m^{2}$, and $\Delta m^{2}$. We find that this uncertainty is dominated by the error in $\theta_{12}$. We also show the 1$\sigma$ (2$\sigma$) favored region of $\delta$ in the dark (light) horizontal green bands. As we discussed in the previous section, there are two solutions for $\delta$ for each value of $\theta_{23}$. Intriguingly, the upper line is right in the middle of the favored range of $\delta$; in particular, $\theta_{23}\simeq 41.5^{\circ}$ gives $\delta\simeq 1.6\pi$, both of which are within the $1\sigma$ allowed region. Consequently, this model predicts $\delta\simeq 1.59\pi$–$1.70\pi$ ($1.54\pi$–$1.78\pi$) within $1\sigma$ ($2\sigma$). Future neutrino experiments can test this prediction through precision measurements of $\theta_{23}$ and $\delta$ [67].

Next, we evaluate the neutrino masses $m_{i}$, which are shown in Fig. 3(a) as functions of $\theta_{23}$. Here, the other parameters are fixed to be their best-fit values. We see that all of these masses are predicted to be $\gtrsim\sqrt{\Delta m^{2}}\simeq 5\times 10^{-2}$ eV. We also plot the sum of these neutrino masses as a function of $\theta_{23}$ in Fig. 3(b), where the dark (light) red band shows the uncertainty coming from the 1$\sigma$ ($2\sigma$) errors in the parameters other than $\theta_{23}$. In this case, it turns out that the dominant contribution to the uncertainty (except for that from the error in $\theta_{23}$) comes from the error in $\theta_{13}$, though the error in $\Delta m^{2}$ also gives a sizable contribution. We also show in the black dashed line the present limit imposed by the Planck experiment: $\sum_{i}m_{i}<0.23$ eV (Planck TT+lowP+lensing+ext) [68].888If we exclude the Planck lensing data, we obtain a slightly stringent bound: $\sum_{i}m_{i}<0.17$ eV [68]. From this figure, we find that a wide range of the parameter region predicts a value of $\sum_{i}m_{i}$ which is below the present limit, though a part of the parameter region has already been disfavored by the Planck limit. We however note that this bound relies on the standard cosmological history, and thus if some new physics effects modify the cosmological evolution, then this bound may significantly be relaxed. In any case, our model predicts a rather large value of the sum of the neutrino masses, $\sum_{i}m_{i}\simeq 0.14$–0.22 eV (0.12–0.40 eV) at $1\sigma$ ($2\sigma$), which may be probed in future cosmological observations.

These relatively large values of $m_{i}$ open up a possibility of testing this model in neutrinoless double-beta decay experiments. The rate of neutrinoless double-beta decay is proportional to the square of the effective Majorana neutrino mass $\langle m_{\beta\beta}\rangle$, which is given by

[TABLE]

It should be emphasized that, in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model, not only the neutrino masses $m_{i}$ but also the Majorana phases $\alpha_{2,3}$ are uniquely determined as functions of the other neutrino oscillation parameters. Thus, the value of the effective mass $\langle m_{\beta\beta}\rangle$ is also predicted unambiguously. Note also that this quantity has reflection symmetry with respect to $\delta\to-\delta$ and thus depends only on $\cos\delta$. In Fig. 4, we show $\langle m_{\beta\beta}\rangle$ as a function of $\theta_{23}$, where the dark (light) red band shows the uncertainty coming from the 1$\sigma$ ($2\sigma$) errors in the parameters other than $\theta_{23}$. Currently, the KamLAND-Zen experiment gives the strongest bound on $\langle m_{\beta\beta}\rangle$: $\langle m_{\beta\beta}\rangle<0.061$–0.165 eV [69] where the uncertainty stems from the estimation of the nuclear matrix element for 136Xe. We also show in Fig. 4 the most severe bound from KamLAND-Zen, $\langle m_{\beta\beta}\rangle<0.061$ eV, in the black dashed line [69]. As can be seen, most parameter region predicts a value of $\langle m_{\beta\beta}\rangle$ lower than the strongest bound. At $1\sigma$ ($2\sigma$) level, this model predicts $\langle m_{\beta\beta}\rangle\gtrsim 0.024$ eV (0.017 eV)—this can be within the reach of future neutrinoless double-beta decay experiments [70].

4 Implications for Leptogenesis

In this section, we discuss the implications of our results for the leptogenesis scenario [71], which is one of the most attractive mechanisms to explain the origin of the baryon asymmetry of the Universe. The minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model has three right-handed neutrinos coupled to the Standard Model leptons, and therefore it contains enough ingredients for the leptogenesis: the CP-violating decay of the right-handed neutrino in the early universe can generate the lepton asymmetry, which is converted to the baryon asymmetry via the sphaleron process [72].

As we have seen in the previous sections, the light neutrino mass matrix ${\cal M}_{\nu_{L}}\simeq U^{*}\text{diag}(m_{1},m_{2},m_{3})U^{\dagger}$ in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model is uniquely determined for a given set of neutrino oscillation parameters $\theta_{12},~{}\theta_{23},~{}\theta_{13},~{}\delta m^{2},~{}\Delta m^{2}$, and $\text{sign}(\sin\delta)$. Therefore, the mass matrix of the right-handed neutrino, ${\cal M}_{R}\simeq-{\cal M}_{D}{\cal M}_{\nu_{L}}^{-1}{\cal M}_{D}$, is also tightly constrained, having only three additional free parameters, $\lambda_{e}$, $\lambda_{\mu}$ and $\lambda_{\tau}$. Note that these neutrino Yukawa couplings can be taken to be real and positive by field redefinitions. Therefore, there is no additional phase parameter in the model. By diagonalizing the masses of the right-handed neutrinos, we can rewrite the Lagrangian (3) as

[TABLE]

where $\widehat{N}_{i}^{c}$ are the right-handed neutrino fields in the basis where the masses are diagonalized, $\widehat{\lambda}_{i\alpha}$ are the neutrino Yukawa couplings in that basis, and $M_{i}$ are the masses which are taken to be real and positive. Explicitly, they are given by

[TABLE]

We emphasize again that both masses $M_{i}$ and the couplings $\widehat{\lambda}_{i\alpha}$ are completely determined by the three real parameters $\lambda_{e,\mu,\tau}$ and the oscillation parameters $\theta_{12},~{}\theta_{23},~{}\theta_{13},~{}\delta m^{2},~{}\Delta m^{2}$, $\text{sign}(\sin\delta)$.

One of the most important parameters in the leptogenesis scenario is the asymmetry parameter $\epsilon_{1}$, which represents the lepton asymmetry generated by the decay of the lightest right-handed neutrino. At the leading order, it is given by [73, 74, 75]

[TABLE]

Here, one can see that there is a correlation between the sign of the Dirac phase $\delta$ and the baryon asymmetry of the Universe. Suppose that the sign of the Dirac phase $\delta$ is flipped, $\delta\to-\delta$, while the other input oscillation parameters $\theta_{12},~{}\theta_{23},~{}\theta_{13},~{}\delta m^{2},~{}\Delta m^{2}$ as well as the neutrino Yukawa couplings $\lambda_{e,\mu,\tau}$ are fixed. As discussed in the previous sections, the Majorana phases then flip the sign, $\alpha_{2,3}\to-\alpha_{2,3}$, while the absolute masses of light neutrinos do not change, $m_{i}\to m_{i}$. Thus, the PMNS matrix transforms as $U\to U^{*}$. This then results in ${\cal M}_{\nu_{L}}\to{\cal M}_{\nu_{L}}^{*}$, leading to ${\cal M}_{R}\to{\cal M}_{R}^{*}$, $\Omega\to\Omega^{*}$, $\widehat{\lambda}\to\widehat{\lambda}^{*}$, and eventually $\epsilon_{1}\to-\epsilon_{1}$. This means that, for a given input oscillation parameters $\theta_{12},~{}\theta_{23},~{}\theta_{13},~{}\delta m^{2},~{}\Delta m^{2}$ and the neutrino Yukawa couplings $\lambda_{e,\mu,\tau}$, there is one-to-one correspondence between the sign of the Dirac phase $\delta$ and the sign of the baryon asymmetry of the Universe.

In Fig. 5, we show the asymmetry parameter $\epsilon_{1}/\lambda^{2}$ as a function of the neutrino Yukawa couplings, which are parametrized as $(\lambda_{e},\lambda_{\mu},\lambda_{\tau})=\lambda(\cos\theta,\sin\theta\cos\phi,\sin\theta\sin\phi)$.999We can take account of the one-loop renormalization-group effects by instead regarding $I_{K}^{1/2}(t)\lambda_{\alpha}{\cal I}_{\alpha\alpha}$ ($\alpha=e,\mu,\tau$) as input parameters (see Eqs. (25) and (26)). Thus, our consequence here is robust against one-loop quantum corrections. Note that the parameter $\epsilon_{1}$ scales as $\epsilon_{1}\propto\lambda^{2}$ in this parametrization, as shown in Eq. (33). The oscillation parameters $\theta_{12},~{}\theta_{23},~{}\theta_{13},~{}\delta m^{2},~{}\Delta m^{2}$ are taken to be their best-fit values in Table 1. The sign of the Dirac phase $\delta$ is taken to be negative, $\delta<0$ (or $\delta>\pi$), as it is favored at 2$\sigma$ level. Note that the observed baryon asymmetry of the Universe requires $\epsilon_{1}<0$, because the sphaleron process predicts $n_{B}/n_{L}<0$ [76]. Surprisingly, negative $\epsilon_{1}$ is realized only in the limited regions of the parameter space. This is more clearly seen in Figs. 6, where we show the regions of $\epsilon_{1}<0$ in the $(\theta,\phi)$ plane as red shaded areas. We also show the contours of the right-handed neutrino mass ratios $M_{2}/M_{1}$ and $M_{3}/M_{1}$ in the left and right panels, respectively. As we see, the asymmetry parameter $\epsilon_{1}$ can be negative when some of the right-handed neutrinos are degenerate in mass—$M_{1}\simeq M_{2}$ in the negative $\epsilon_{1}$ region around $\phi\simeq\pi/4$ while $M_{2}\simeq M_{3}$ for $\theta\simeq\pi/2$ and $\phi\simeq 0,\pi/2$. In the latter regions, $\theta\simeq\pi/2$ and $\phi\simeq 0,\pi/2$, the lightest right-handed neutrino mass is much smaller than the other ones, and $\lambda_{e,\tau}\ll\lambda_{\mu}$ ($\lambda_{e,\mu}\ll\lambda_{\tau}$) for $\phi\simeq 0$ ($\phi\simeq\pi/2$). In these cases, the absolute value of the asymmetry parameter $|\epsilon_{1}|$ is found to be quite suppressed, and thus it is rather difficult to obtain a sizable value, say $|\epsilon_{1}|\gtrsim 10^{-6}$.101010Note that the function $f(x)$ in Eq. (34) goes as $f(x)\simeq-3/(2\sqrt{x})$ for $x\gg 1$. Similarly, $|\epsilon_{1}|$ gets small for $\phi\simeq\pi/4$ and $\theta\ll\pi/4$. These observations further restrict the promising parameter region for leptogenesis to be $(\theta,\phi)\simeq(\pi/4,\pi/4)$.

The final baryon asymmetry depends on the production mechanism of the right-handed neutrino. In the case of thermal leptogenesis, the predicted range of the lightest neutrino mass in the present model, $m_{1}\gtrsim 0.03~{}\text{eV}$, corresponds to the so-called strong wash-out region (see, e.g., [77]). Moreover, as we see from Fig. 6, around $(\theta,\phi)\simeq(\pi/4,\pi/4)$, right-handed neutrino masses are close to each other. In this case, to properly estimate the net baryon asymmetry generated via the decay of right-handed neutrinos, we may need to include not only the contribution of the lightest right-handed neutrino but those of the other ones. Flavor effects [78, 79, 80] may also become important. A more detailed study will be given elsewhere [81].

5 Conclusions

We have studied the structure of the neutrino mass matrix in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model. Because of the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry, the structure of both Dirac and Majorana mass terms of neutrinos is tightly restricted, which results in a two-zero minor structure of the neutrino mass matrix. Because of this restriction, all the CP phases and the neutrino masses are uniquely determined. We find that this model gives the following prediction at $1\sigma$ ($2\sigma$) level:

• •

Dirac CP phase: $\delta\simeq 1.59\pi$–$1.70\pi$ ($1.54\pi$–$1.78\pi$),

• •

the sum of the neutrino masses: $\sum_{i}m_{i}\simeq 0.14$–0.22 eV ($0.12$–$0.40$ eV),

• •

the effective mass for the neutrinoless double beta decay: $\langle m_{\beta\beta}\rangle\simeq 0.024$–0.055 eV (0.017–0.12 eV).

They are totally consistent with the current experimental limits, and hold independently of the U(1)${}_{L_{\mu}-L_{\tau}}$ breaking scale and the Majorana mass scale. In this sense, the above predictions are the generic features of the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model. Remarkably, these predictions can be tested in various neutrino experiments in the near future—regardless of the scale of the U(1)${}_{L_{\mu}-L_{\tau}}$ symmetry breaking—which we believe sheds light on the gauge structure of physics beyond the Standard Model.

We have also discussed the implications of the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model for the leptogenesis scenario, and found that the correct sign of the baryon asymmetry of the Universe can be obtained only in the limited regions of the parameter space, as the right-handed neutrino mass structure is also severely restricted in this model. In particular, the observed value of baryon asymmetry can be realized when right-handed neutrinos are degenerate in mass, which requires a further detailed study to assess the viability of leptogenesis in this model [81].

Acknowledgments

This work is supported by the Grant-in-Aid for Scientific Research (No.26104001 [KH], No.26104009 [KH], No.26247038 [KH], No.26800123 [KH], No.16H02189 [KH], No.17K14270 [NN]), and by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. The work of KA is supported by the Program for Leading Graduate Schools, MEXT, Japan.

Appendix

Appendix A Miscellaneous formulae

Here we give formulae that are useful for the study of the neutrino mass structure in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model.

A.1 Neutrino mass matrix ${\cal M}_{\nu_{L}}$

The light neutrino mass matrix ${\cal M}_{\nu_{L}}$ in Eq. (6) can be expressed in terms of the Lagrangian parameters in Eq. (3) as111111We obtain a different result from that given in Ref. [21].

[TABLE]

The determinant of this mass matrix is given by

[TABLE]

We find that this determinant vanishes if and only if $\lambda_{\nu}=0$ ($\nu=e$, $\mu$, or $\tau$). In this case, the mass matrix ${\cal M}_{\nu_{L}}$ is block-diagonal, and thus cannot reproduce the required neutrino mixing angles.

A.2 $R_{2}$ and $R_{3}$

The functions $R_{2}(\delta)$ and $R_{3}(\delta)$ defined in Eqs. (13) and (14), respectively, are expressed in terms of neutrino oscillation parameters as

[TABLE]

Using Eq. (19) together with Eq. (A.3), we find

[TABLE]

This expression shows that $m_{1}$ depends only on $\cos\delta$ (not $\sin\delta$), as we have argued in Sec. 2.

Using Eqs. (12), (A.3), and (A.4), we determine the Majorana CP phases $\alpha_{2,3}$ as functions of $\delta$, and show them in Fig. 7. Here, we fix the other neutrino oscillation parameters to be their best-fit values. The green band depicts the 2$\sigma$ favored region of $\delta$ predicted in the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model. As can be seen, sizable Majorana CP phases are predicted in this model. In addition, we find $\alpha_{2,3}(-\delta)=-\alpha_{2,3}(\delta)$ (or $\alpha_{2,3}(\pi-\delta)=-\alpha_{2,3}(\pi+\delta)$), as discussed in Sec. 2.

A.3 Cubic equation for $\cos\delta$

Here, we show a cubic equation whose real solution in terms of $x$ gives $\cos\delta$:

[TABLE]

where

[TABLE]

In the limit of $\epsilon\to 0$, the above equation leads to

[TABLE]

or

[TABLE]

The discriminant of the quadratic equation (A.8) is given by

[TABLE]

which is negative as $\cos 4\theta_{12}+\cos 4\theta_{23}\simeq-1.63<0$. Thus, Eq. (A.8) does not give a real solution. On the other hand, Eq. (A.9) gives

[TABLE]

which agrees to Eq. (21). From the above derivation, we see that the solution (A.11) approximates the real solution of the cubic equation (A.6) with an accuracy of ${\cal O}(\epsilon)={\cal O}(\delta m^{2}/\Delta m^{2})$.

Appendix B U(1)${}_{L_{e}-L_{\mu}}$ and U(1)${}_{L_{e}-L_{\tau}}$

In this section, we examine the neutrino mass structure in the minimal gauged U(1)${}_{L_{e}-L_{\mu}}$ and U(1)${}_{L_{e}-L_{\tau}}$ models and show that it is unable to obtain a solution that is consistent with the observed values of the neutrino oscillation parameters.121212A similar conclusion was also reached in Ref. [50].

B.1 U(1)${}_{L_{e}-L_{\mu}}$

Following an analysis similar to that in Sec. 2, we find that the $(e,e)$ and $(\mu,\mu)$ components in the inverse of the neutrino mass matrix vanish in the minimal gauged U(1)${}_{L_{e}-L_{\mu}}$ model. The corresponding two vanishing conditions are

[TABLE]

Solving these equations, we have

[TABLE]

with

[TABLE]

In Fig. 8(a), we plot the mass ratios $m_{2}/m_{1}$ and $m_{3}/m_{1}$ as functions of $\delta$. As we see from this figure, $m_{2}<m_{1}$ is predicted for any value of $\delta$, and thus there is no solution which gives an allowed pattern of neutrino mass spectrum.

B.2 U(1)${}_{L_{e}-L_{\tau}}$

In this case, the $(e,e)$ and $(\tau,\tau)$ components in ${\cal M}^{-1}_{\nu_{L}}$ are zero, which leads to

[TABLE]

These equations read

[TABLE]

with

[TABLE]

Using these equations, we plot the mass ratios $m_{2}/m_{1}$ and $m_{3}/m_{1}$ as functions of $\delta$ in Fig. 8(b). Again, $m_{2}<m_{1}$ over the whole range of $\delta$, and thus this model cannot provide a desirable neutrino mass ordering.