Twenty years after the discovery of $\mu-\tau$ symmetry

TL;DR

This paper reviews 20 years of $- au$ symmetry in neutrino physics, emphasizing systematic model extensions and their integration into grand unified theories, amidst experimental progress.

Contribution

It provides a systematic approach to extending $- au$ symmetry models and explores their application within grand unified theories.

Findings

Systematic extensions of $- au$ symmetry models are proposed.

Application of these models to GUT frameworks is discussed.

The paper highlights the need for more fundamental theoretical outlooks.

Abstract

It has passed 20 years after we proposed $\mu-\tau$ symmetry in light neutrino mass matrix. This model is simple but reproduced the characterestic properties of lepton sector. After that, during the experimental developments, there have appeared so many extensions but most of those phenomenological models are lacking systematic outlooks towards more fundamental theories. In this paper, we try to consider rather systematic model extensions and application to GUT model.

Equations

$$M_{\nu}=\left(\begin{array}[]{ccc}0&A&\pm A\\ A&B&C\\ \pm A&C&B\end{array}\right).$$

$${\sin }^2{\theta }_{12}=0.304\pm 0.014$$

$$\Delta m_{21}^2=(7.53\pm 0.18)\times 10^{-5}\text{mbox}{eV}^2$$

$${\sin }^2{\theta }_{23}=0.51\pm 0.05\text{mbox}(normalmasshierarchy)$$

$${\sin }^2{\theta }_{23}=0.50\pm 0.05\text{mbox}(invertedmasshierarchy)$$

$$\Delta m_{32}^2=(2.44\pm 0.06)\times 10^{-3}\text{mbox}{eV}^2\text{mbox}(normalmasshierarchy)$$

$$\Delta m_{32}^2=(2.51\pm 0.06)\times 10^{-3}\text{mbox}{eV}^2\text{mbox}(invertedmasshierarchy)$$

$${\sin }^2{\theta }_{13}=(2.19\pm 0.12)\times 10^{-2}$$

$$U=\left(\begin{array}[]{ccc}c_{13}c_{12},&c_{13}s_{12},&s_{13}e^{-i\delta}\\ -c_{23}s_{12}-s_{23}c_{12}s_{13}e^{i\delta},&c_{23}c_{12}-s_{23}s_{12}s_{13}e^{i\delta},&s_{23}c_{13}\\ s_{23}s_{12}-c_{23}c_{12}s_{13}e^{i\delta},&-s_{23}c_{12}-c_{23}s_{12}s_{13}e^{i\delta},&c_{23}c_{13}\end{array}\right)\equiv\left(\begin{array}[]{ccc}c_{13}c_{12},&c_{13}s_{12},&s_{13}e^{-i\delta}\\ U_{\mu 1},&U_{\mu 2},&s_{23}c_{13}\\ U_{\tau 1},&U_{\tau 2},&c_{23}c_{13}\end{array}\right)$$

$$M_{\nu}=U\left(\begin{array}[]{ccc}-m_{1},&0,&0\\ 0,&m_{2},&\\ 0,&0,&m_{3}\end{array}\right)U^{T}$$

$${\left(M_{\nu }\right)}_{11}=-m_1{c}_{12}^2{c}_{13}^2+m_2{s}_{12}^2{c}_{13}^2+m_3{s}_{13}^2{e}^{-2i\delta },$$

$${\left(M_{\nu }\right)}_{12}$$

$${\left(M_{\nu }\right)}_{13}$$

$${\left(M_{\nu }\right)}_{22}$$

$${\left(M_{\nu }\right)}_{33}$$

$$M_{\nu }$$

$$=$$

$$-c_{12}^2{m}_1+s_{12}^2{m}_2=0,$$

$$-m_1$$

$$m_2$$

$$m_3$$

$$⟨m{⟩}_{ee}\equiv \sum _i{U}_{ei}{m}_i{U}_{ei}=0.$$

$$(a)=\left(\begin{array}[]{ccc}0&A&A\\ A&B&C\\ A&C&C\end{array}\right)~{}~{}(b)=\left(\begin{array}[]{ccc}0&A&A\\ A&B&B\\ A&B&C\end{array}\right)~{}~{}(c)=\left(\begin{array}[]{ccc}0&A&A\\ A&B&C\\ A&C&B\end{array}\right)~{}~{}(d)=\left(\begin{array}[]{ccc}0&A&-A\\ A&B&C\\ -A&C&B\end{array}\right)$$

$$\left(\begin{array}[]{ccc}0&c_{2}\sqrt{m_{1}m_{2}}&-s_{2}\sqrt{m_{1}m_{2}}\\ c_{2}\sqrt{m_{1}m_{2}}&(-m_{1}+m_{2})c_{2}^{2}+m_{3}s_{2}^{2}&(m_{1}-m_{2}+m_{3})c_{2}s_{2}\\ -s_{2}\sqrt{m_{1}m_{2}}&(m_{1}-m_{2}+m_{3})c_{2}s_{2}&(-m_{1}+m_{2})s_{2}^{2}+m_{3}c_{2}^{2}\end{array}\right),$$

$$M_{\nu}^{\prime}=\left(\begin{array}[]{ccc}D&A&\pm A\\ A&B&C\\ \pm A&C&B\end{array}\right).$$

$$\tan 2{\theta }_{12}=\frac{2\sqrt{2}A}{B\pm C-D},$$

$$-m_1$$

$$m_2$$

$$m_3$$

$$U_{TB}=\left(\begin{array}[]{ccc}\frac{2}{\sqrt{6}}&\frac{1}{\sqrt{3}}&0\\ -\frac{1}{\sqrt{6}}&\frac{1}{\sqrt{3}}&-\frac{1}{\sqrt{2}}\\ -\frac{2}{\sqrt{6}}&\frac{1}{\sqrt{3}}&\frac{1}{\sqrt{2}}\end{array}\right).$$

$$M_{\nu}=\left(\begin{array}[]{ccc}D&A&\pm A^{*}\\ A&B&C\\ \pm A^{*}&C&B^{*}\end{array}\right)~{}~{}(\mbox{double sign corresponds}),$$

$$U=\left(\begin{array}[]{ccc}u_{1}&u_{2}&u_{3}\\ v_{1}&v_{2}&v_{3}\\ \pm v_{1}^{*}&\pm v_{2}^{*}&\pm v_{3}^{*}\end{array}\right)~{}~{}(\mbox{double sign correspond to that of (\ref{cform})})$$

$$M_{\nu }{M}_{\nu }^{†}=U{D}_{\nu }^2{U}^{†}.$$

$$M_{\nu}M_{\nu}^{\dagger}=UD_{\nu}^{2}U^{\dagger}=\left(\begin{array}[]{ccc}z&w&\pm w^{*}\\ w^{*}&x&y\\ \pm w&y^{*}&x\end{array}\right).$$

$${\theta }_{23}=\frac{\pi }{4},\delta =\pm \frac{\pi }{2}\text{mbox}forminussignature$$

$${\theta }_{23}=-\frac{\pi }{4},\delta =\pm \frac{\pi }{2}\text{mbox}forplussignature$$

$$M_{f}=P_{f}^{\dagger}\hat{M_{f}}P_{f}\equiv P_{f}^{\dagger}\left(\begin{array}[]{ccc}D_{f}&A_{f}&A_{f}\\ A^{\prime}_{f}&B_{f}&C_{f}\\ A^{\prime}_{f}&C_{f}&B_{f}\end{array}\right)P_{f}$$

Full text

Twenty years after the discovery of $\mu-\tau$ symmetry.

Takeshi Fukuyama 111E-mail:fukuyama@se.ritsumei.ac.jp

*Research Center for Nuclear Physics (RCNP),

Osaka University, Ibaraki, Osaka, 567-0047, Japan*

It has passed 20 years after we proposed $\mu-\tau$ symmetry in light neutrino mass matrix. This model is simple but reproduced the characterestic properties of lepton sector. After that, during the experimental developments, there have appeared so many extensions but most of those phenomenological models are lacking systematic outlooks towards more fundamental theories. In this paper, we try to consider rather systematic model extensions and application to GUT model.

1 Introduction

Twenty years ago we proposed first in the world $\mu-\tau$ symmetry in the neutrino mass matrix model [1],

[TABLE]

in the charged lepton diagonal base. Here $A,B,C$ are real and its components are invariant under $\mu-\tau$ exchange. (1,2) and (1,3) components are equal up to phase convention. This matrix, therefore, has been called $\mu-\tau$ symmetric mass matrix. This leads immediately to $\theta_{23}=\mp\pi/4$ and $\theta_{13}=0$ (double sign in the same order as (1)). This matrix represents the characterestic pattern of the mixing angles which is quite different from that of quark sector. The vanishing (1,1) component leads to the small mixing angle (SMA) solution on $\theta_{12}$ (See Eq.(19)), which had survived with large mixing angle (LMA) solution at that time. However, KamLAND [2] selected the larger part of solar neutrino angles, and we may set nonzero parameter $D$ in place of the vanishing (1,1) component without breaking $\mu-\tau$ symmetry. In 2013, Daya-Bay[3] made surprise the unexpectedly large $\theta_{13}$. It is impressing that our minimal SO(10) model [4] discussed in section 4 has suffered from large $\theta_{13}$ before Daya-Bay.

The observed data of leptonic mixing matrix nowaday are summarized as [5]

[TABLE]

Even in these refined data, our $\mu-\tau$ symmetric model does not lose its significance since such simple and real symmetric model is basically an idealized model and remains valid as the zero’th order approximation of more sophiscated models.

Indeed there have appeared a vast variety of papers during new experimental developments. Unfortunately, most of those phenomenological models are lacking systematic analyses valid for theoretical developments from phenomenological model to more fundamental one.

In this paper we reconsider $\mu-\tau$ symmetry in these experimental backgrounds and try to fill this deficit.

This paper is organized as follows. In section 2 we review the original $\mu-\tau$ symmetric model. This model is extended in section 3. In section 4 we argue the correlaton with GUTs. Section 5 is devoted to discussions.

2 $\mu-\tau$ symmetric model

Since neutrino oscillation experiments are wholly insensitive to the Majora phases, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is in general written in the form

[TABLE]

Here $c_{ij}=\mbox{cos}\theta_{ij},~{}s_{ij}=\mbox{sin}\theta_{ij}$ as usual. Neutrino mass matrix is written as

[TABLE]

in the charged lepton diagonal base, where we can set $-1$ for $m_{1}$ using the rephasing. Its explicit components are

[TABLE]

[TABLE]

[TABLE]

As is easily checked, $\theta_{23}=\pm\pi/4$ and $\theta_{13}=0$ is the unique solutions for real $\left(M_{\nu}\right)_{12}=\mp\left(M_{\nu}\right)_{13}$, $(M_{\nu})_{22}=(M_{\nu})_{33}$ relations (double sign corresponds).

If we adopted $\theta_{23}=\pi/4$ and $\theta_{13}=0$, then neutrino mass matrix becomes

[TABLE]

where $c_{1}(s_{1})\equiv\cos\theta_{12}(\sin\theta_{12})$ for brevity. In (LABEL:mixing2) we assumed further

[TABLE]

and we obtain the mass matrix of (1). Neutrino masses are expressed in terms of $A,B,C$,

[TABLE]

The double sign corresponds to (1). Eq.(19) indicates that neutrinoless double beta decay does not happen in this limit, since

[TABLE]

(19) favored the small mixing angle solution for $\theta_{12}$ which still had survived at that time. Also the vanishing (1,1) component is interested in connection with seesaw invariant mass matrix [1].

[TABLE]

(c) is transformed to (d) by the interchange of $C$ to $-C$ and these are physically equivalent as follows. If we leave $\theta_{23}$ as a free parameter and keep the assumtion (19), then $M_{\nu}$ is reduced to

[TABLE]

where $c_{2}(s_{2})\equiv\cos\theta_{23}(\sin\theta_{23})$ for brevity. Therefore, (c) and (d) are corresponding to $s_{23}=-\pi/4$ and $\pi/4$, repecetively. $\theta_{23}$ has been determined from the mixing factor $\sin^{2}2\theta_{23}$ and they are equivalent. (a) and (b) are also substantially same and from Eq.(LABEL:mixing2) they are enforced to $m_{3}\approx 0$. This is the case of inverted hierarchy. This symmetric and seesaw invariant concepts is extended to two-zero texture [6].

Eq.(19) was imposed as it enables us to fix all three masses by the same three parameters as (2). You can easily generalize this simple model as

[TABLE]

Eqs.(19) and (2) in this case are generalized to

[TABLE]

[TABLE]

If $B\pm C-D=A,~{}U$ goes to tri-bi-maximal case [7],

[TABLE]

3 Extension of $\mu-\tau$ symmetric model

The original $\mu-\tau$ symmetric model was real and can not involve CP phases. So a naive extension is to extend it to a complex and symmetric mass matrix retaining $\mu-\tau$ symmetry. The reasoning why we adhere to a symmetric matrix will be explained in section 4. Eq.(1) is a real symmetric matrix. Its naive extension is

[TABLE]

where $D,C\in\mathcal{R},~{}A,B\in\mathcal{C}$. This matrix is diagonalized by the unitary matrix [8],

[TABLE]

as

[TABLE]

Then

[TABLE]

In Eqs. (28), (29), and (31), double sign corresponds. Here $x,z\in\mathcal{R},~{}y,w\in\mathcal{C}$. In this complex form we have, in addition to $\theta_{23}=\frac{\pi}{4},~{}\theta_{13}=0$, another solution,

[TABLE]

$\delta=-\frac{\pi}{2}$ is interesting since it is the global minimum (though $1~{}\sigma$) [9].

One strategy for extending $\mu-\tau$ symmetry is the following: The extensions not only explain the leptonic CP phase but also must include quark sector. This is because we are considering GUT as its more fundamental final correspondent. One of such examples preserves $\mu-\tau$ symmetry up to phase but breaks the symmetric property of mass matrices [10] like,

[TABLE]

where

[TABLE]

and $f$ includes up-type and down-type quark. First we diagonalize $\hat{M_{f}}$ by two orthogonal matrices $O_{f1}$ and $O_{f2}$ as

[TABLE]

[TABLE]

where $U_{TB}$ is a tri-bi-maximal mixing matrix (27). We found that these matrices are consistent with the experimental data of CKM mixing matrix. This is the extension to quark sector but is left on the same phenomenological level as the original work of lepton sector. Hereafter we restrict ourselves in symmetric mass matrices again. If we involve quark sector, it must reveal some higher symmetric (more fundamental) new character. In this sence, though there are a vast variety of these extensions, most of these phenomenological extensions have no systematic idea leading to more fundamental theoretical models. For the route from (low energy) phenomenological model to (high energy) more fundamental one, some symmetry must play an important role.

Let us explain it in the well known example: QED lagrangian has $U(1)$ and Lorentz invariances,

[TABLE]

Under T transformation, ${\bf E}\to-{\bf E}$ and ${\bf B}\to{\bf B}$, and therefore $\theta\to-\theta$, if we preserve T invariance. In another word, T invariance requires $\theta=0$. Lorentz invariance breaks to spatial rotation invariance for dielctics,

[TABLE]

Thus, permeability ($\mu$) and permittivity ($\epsilon$) characterize this symmetry breaking. Eq.(39) ((38)) may be considered as a phenomenological (a fundamental) model. According to this general idea, how the phenomenological $\mu-\tau$ symmetry is incorporated to higher symmetry group or more comprehensive model ? It is natural to incorporate quark sector in this higher symmetric world. We consider here $A_{4}$ group as a candidate for it [11]. $A_{4}$ is the four degreed symmetry group with even permutation whose elements we denote as $(a_{1},a_{2},a_{3},a_{4})\cite[cite]{[\@@bibref{}{Altarelli}{}{}]}.~{}A_{4}$ is generated by the $S$ and $T$ and their products, which satisfy

[TABLE]

The three-dimensional unitary representation, in a basis where the element $S$ is diagonal, is built up from:

[TABLE]

Let us practice the transformation, for instance, $ST$ to $V~{}\equiv~{}(a_{1},a_{2},a_{3})^{T}$

[TABLE]

The rule of the game for reading the permutation group of four degree from three dimensional vector is to make plus element change to $a_{4}$ and do minus signs interchange, and $ST$ corresponds to $(a_{4},a_{1},a_{3},a_{2})$. Thus $S$ means the $2\leftrightarrow 3$ symmetry and $T$ does cyclic permutation or equivalently $Z_{3}$. Mathematically this is the elementary example of Sylow’s theorem [13]. The order of $A_{4}$ is $12=2^{2}\times 3$, and it is the product of normal subgroup $V_{4}$, composed of $(1,2)(3,4),~{}(1,3)(2,4),~{}(1,4)(2,3),~{}1$ and $Z_{3}$. Thus $A_{4}\supset(2\leftrightarrow 3)~{}\mbox{symmetry}~{}\times~{}Z_{3}$. Namely, $\mu-\tau$ symmetry may be considered the residual symmetry broken from $A_{4}$. This fact is very important for the model buiding of more fundamental theories. So far we have not considered $Z_{3}$, and let us consider how $Z_{3}$ symmetry appears. Corresponding to this extension, we generalize from the lepton sector to the quark-lepton sector, denoting their fields as $\psi_{i}$ ($i$: generation), and call 2-3 symmetry instead of $\mu-\tau$ summetry in that case.

We assign $Z_{3}$ charge of each generation of fermions so as to be compatible with 2-3 symmetry [14],

[TABLE]

where $\omega^{3}=+1$. Then, the bilinear terms $\overline{q}_{Li}u_{Rj},~{}\overline{q}_{Li}d_{Rj},~{}\overline{l}_{Li}\nu_{Rj},~{}\overline{l}_{Li}e_{Rj}$ and $\overline{\nu_{Ri^{c}}}\nu_{Rj}$ are transformed as follows:

[TABLE]

where

[TABLE]

Therefore, if we assume two SU(2) doublet Higgs scalars $H_{1}$ and $H_{2}$, which are transformed as

[TABLE]

the Yukawa interactions are given as follows

[TABLE]

where

[TABLE]

Therefore,

[TABLE]

In (49), the symbol $*$ denotes non-zero quantities. Here, in order to give heavy Majorana masses of the right-handed neutrinos $\nu_{R}$, we have assumed an SU(2) singlet Higgs scalar $\Phi^{0}$, which is transformed as $H_{1}$. Mass matrices are sums of $Y_{(1)}$ and $Y_{(2)}$ and their (1,1) element must be vanished:

[TABLE]

Thus zero-texture model becomes another useful character as well as symmetric property. Then how far can we go along this line of thought ? In this case, neutrino mass matrix may have the special property of seesaw mechanism [15], and the concept of the seesaw invariance plays an important role [1, 6]. Two-zero texture is interesting from the parameter counting. Mass matrix $M_{\nu}$ is determined by 9 out-put parameters, 3 masses, 3 angles and 3 CP phases (one Dirac $\delta$ and 2 Majorana phases) in the charged lepton flavour diagonal base. Two-zero texture gives four constraints and 9-4=5 in-put free parameters [16, 17, 18]. Among others, the following textures are very important,

[TABLE]

They are related by $\mu-\tau$ excahange, $M_{\nu}^{(2)}=P_{23}M_{\nu}^{(1)}P_{23}$. Here

[TABLE]

As mentioned above, five parameters remain free in two-zero texture model. Therefore, if $\theta_{ij},~{}\Delta m_{sol}^{2},~{}\Delta m_{atm}^{2}$ are determined, $\delta,~{}\rho,~{}\sigma$ are predicted.

[TABLE]

for $M_{\nu}^{(2)}$ case [18], where $t_{ij}\equiv\tan\theta_{ij}$. For $M_{\nu}^{(1)}$ case is obtained by replacing $t_{23}$ by $-1/t_{23}$. Thus $M_{\nu}^{(1)}$ and $M_{\nu}^{(2)}$ give similar results. This fact is used in GUT formulation as as will be discussed in the next section.

4 Mass matrix model and GUT

GUT models basically search the vertical structure of quark-lepton of one generation. Inter-family (horizontal) relations like Yukawa structure are not predicted. $\mu-\tau$ symmetry may be clue to this extension.

In the previous section, we considered that $2\leftrightarrow 3$ symmetry and $Z_{3}$ suggest zero texture solution. So we consider GUT implimented witbh texture. For that purpose we set two requirements

• •

GUT model itself must have few ambiguities and be predictive enough.

• •

The reliable mass texture model should be adopted.

The $SO(10)$ grand unified theory (GUT) can provide the most promising framework to unify quarks and leptons, because the entire SM matter contents of each generation (including a right-handed neutrino) can be unified in a single irreducible representation, $\bf 16$. A particular attention has been paid to the renormalizable minimal $SO(10)$ model, where two Higgs multiplets $\{{\bf 10}\oplus{\bf\overline{126}}\}$ are utilized for the Yukawa couplings with the matter representation [19]. The couplings to the 10 and $\overline{\bf 126}$ Higgs fields can reproduce realistic charged fermion and neutrino mass matrices using their phases thoroughly [20, 21]. $\overline{\bf 126}$ Higgs is selected since it includes $({\bf 10,~{}1,~{}3})$ and $({\bf\overline{10},~{}1,~{}3})$ under the Pati-Salam subgroup which induce type I and type II seesaw mechanism, respectively. Yukawa coupling is given by

[TABLE]

where ${\bf 16}_{i}$ is the matter multiplet of the $i$-th generation, $H_{10}$ and $H_{126}$ are the Higgs multiplet of 10 and $\overline{\bf 126}$ representations.

The Yukawa coupling, after SO(10) symmetry is broken down to the standard model, is given as follows:

[TABLE]

Here $c_{10},~{}c_{126},~{}c_{T},~{}c_{R}$ are comlex constants. It should be remarked that mass matrices are complex and symmetric matrices because of the group property of ${\bf 10,\overline{126}}$ representations. Here we proceed to incorporate two-zero texture i n this model. We adopted the two-zero texture mass $M_{\nu}^{(2)}$ of (51).

Unfortunately, the data fittings have been performed by inputting quark sector spectrum and outputted the lepton sector. In this approach, neutrino mass texture is contaminated by the special base adopted in quark sector and shows no clear texture in $M_{\nu}$. The reason why we adopted the quark sector as input data is that the leptoin sector had been more ambiguous than the quark sector. Nowaday, however, the situation changed. The lepton parameters are more accurate rather than the quark ones, and a large threshold correction is expected in the quark sector in SUSY models. In that sence, it is better to perform the fitting by inputting the parameters in the lepton sector. The formulation is presented in not only such a practical purpose, but also to make clear the property of the solution with $v_{R}\approx 10^{16}$ GeV. Here $v_{R}$ is the typical intermediate enrygy scale, and usualy adopted $10^{13}$ GeV spoils the gauge coupling unifications [22]. Real data fitting revealed that in $v_{R}\approx 10^{16}$ GeV solution, the down quark mass is smaller than the observation, and (1,1) and (1,2) elements of $(M_{\nu})$ are smaller than the other elements in the fit result. The deficit of down quark mass can be considered as the threshold correction. Under the assumtion of $(M_{\nu})_{11}=(M_{\nu})_{12}$ lead us to [23]

[TABLE]

In Fig.1 we plot the relation between $\delta_{PMNS}$ and $\theta_{23}$ in the assumption. Of course those two mass matrix elements are not exactly zero in the fits, and provides a guide to understand the fit results for the prediction of the PMNS phase depending on the mixing angle. In Fig.2, we show the plot of proton decay, $\tau(p\to K\overline{\nu})$, in the $\delta_{PMNS}-\theta_{23}$ plane. As expected , the partial proton lifetime is larger near the curve of zero-texture in Fig.1. Near the curve, the lifetime is about 10 times bigger than the current experimental bound $\tau(p\to K\overline{\nu})>0.59\times 10^{34}$ years. Please see [23, 24] for the detail.

5 Discussion

We have not connected $2-3$ symmetry or $A_{4}$ symmetry directly with GUT symmetry SO(10). Naively it seems to be natural to consider SO(10)$\times A_{4}$ symmetry. However, the merit of renormalizable minimum SO(10) GUT is that all mass matrices, as you see in (4), are represented by only two mass matrices $M_{10}$ and $M_{126}$. So it makes this model very predictive and few room to add any assumption. If we incorporate $A_{4}$ into SO(10) GUT naively, it brings about many ambiguities on how to specify $A_{4}$ to Higgs and inflavons et.al. [26]. This spoils the high predictivity of the minimal SO(10) GUT. On the other hand, two-zero texture is very useful for full data fitting scan because we know $\left(M_{\nu}\right)_{11}\approx\left(M_{\nu}\right)_{12}\approx 0$ phenomenologically. So we can scan around this neighbourhood. And indeed we have found very good fitting around this solution. It is very interesting that long proton decay is obtained along this solution like Fig.2. Such collaboration of GUT (fundamental theory) with phenomenological model is unprecedented. $2-3$ symmetry and two-zero texture make clear the GUT solution as well as the practical usefulness of comprehensive data fittings. Thus we can not only fix all mass matrices but also predict many unobserved parameters, like proton decay and lepton flavour violation. Moreover, using the SUSY breaking boundary condition indicateded in [23], we can fit the other almost all known bounds like LFV and $\Omega_{DM}$ BR$(B_{s}^{0}\to\mu^{+}\mu^{-})$ etc. [27]

Acknowledgments

This paper is based on a talk given in a mini-workshop on gquarks, leptons and family gauge bosons hDecember 26-27, 2016 Osaka. The author expresses his sincere thanks to the organizer Y.Koide. He is grateful to H.Nishiura for his excellent collaboration on $\mu-\tau$ symmetry. He is also deeply indebted to Y.Mimura and K.Ichikawa in the recent works on the SO(10) GUT. This work is supported in part by Grant-in-Aid for Science Research from Japan Ministry of Education, Science and Culture (No. 26247036).