Generalized Bottom-Tau unification, neutrino oscillations and dark matter: predictions from a lepton quarticity flavor approach

TL;DR

This paper presents an $A_4$ flavor model linking dark matter stability with neutrino properties, predicting specific mass relations, CP violation, and neutrino mass ordering, with implications for neutrino oscillation experiments.

Contribution

It introduces a novel $A_4$ flavor extension of the Standard Model that correlates dark matter stability with neutrino Dirac nature and predicts specific neutrino oscillation parameters.

Findings

Predicts a generalized bottom-tau mass relation involving all families.

Induces small neutrino masses via a seesaw mechanism.

Requires significant CP violation in neutrino oscillations.

Abstract

We propose an $A_4$ extension of the Standard Model with a Lepton Quarticity symmetry correlating dark matter stability with the Dirac nature of neutrinos. The flavor symmetry predicts (i) a generalized bottom-tau mass relation involving all families, (ii) small neutrino masses are induced a la seesaw, (iii) CP must be significantly violated in neutrino oscillations, (iv) the atmospheric angle $\theta_{23}$ lies in the second octant, and (v) only the normal neutrino mass ordering is realized.

Tables

Table 1. Table 1: Charge assignments for leptons, quarks, scalars ( Φ i u superscript subscript Φ 𝑖 𝑢 \Phi_{i}^{u} , Φ i d superscript subscript Φ 𝑖 𝑑 \Phi_{i}^{d} and χ i subscript 𝜒 𝑖 \chi_{i} ) as well as “dark matter sector” ( ζ 𝜁 \zeta and η 𝜂 \eta ). Here 𝐳 𝐳 \mathbf{z} is the fourth root of unity, i.e. 𝐳 4 = 1 superscript 𝐳 4 1 \mathbf{z}^{4}=1 .

Fields	$SU{\left(2\right)}_L$	$A_4$	$Z_4$	Fields	$SU{\left(2\right)}_L$	$A_4$	$Z_4$
${\overline{L}}_i$	$\mathrm{𝟐}$	$\mathrm{𝟑}$	${𝐳}^3$	${\nu }_{e,R}$	$\mathrm{𝟏}$	$\mathrm{𝟏}$	$𝐳$
${\overline{N}}_{i,L}$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	${𝐳}^3$	${\nu }_{\mu ,R}$	$\mathrm{𝟏}$	${\mathrm{𝟏}}^{\prime }$	$𝐳$
$N_{i,R}$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	$𝐳$	${\nu }_{\tau ,R}$	$\mathrm{𝟏}$	${\mathrm{𝟏}}^{\prime \prime }$	$𝐳$
$l_{i,R}$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	$𝐳$	$d_{i,R}$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	$𝐳$
${\overline{Q}}_{i,L}$	$\mathrm{𝟐}$	$\mathrm{𝟑}$	${𝐳}^3$	$u_{i,R}$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	$𝐳$
${\Phi }_1^u$	$\mathrm{𝟐}$	$\mathrm{𝟏}$	$\mathrm{𝟏}$	${\chi }_i$	$\mathrm{𝟏}$	$\mathrm{𝟑}$	$\mathrm{𝟏}$
${\Phi }_2^u$	$\mathrm{𝟐}$	${\mathrm{𝟏}}^{\prime }$	$\mathrm{𝟏}$	$\eta$	$\mathrm{𝟏}$	$\mathrm{𝟏}$	${𝐳}^2$
${\Phi }_3^u$	$\mathrm{𝟐}$	${\mathrm{𝟏}}^{\prime \prime }$	$\mathrm{𝟏}$	$\zeta$	$\mathrm{𝟏}$	$\mathrm{𝟏}$	$𝐳$
${\Phi }_i^d$	$\mathrm{𝟐}$	$\mathrm{𝟑}$	$\mathrm{𝟏}$

Equations

$$\frac{m_b}{\sqrt{m_d{m}_s}}=\frac{m_{\tau }}{\sqrt{m_e{m}_{\mu }}}\text{leavevmode}.$$

$${\mathcal{L}}_{\mathrm{Yuk},\nu }$$

$$+$$

$$\displaystyle M_{\nu,N}\,\,=\,\left(\begin{array}[]{cccccc}0&0&0&a^{\prime}_{1}&0&0\\ 0&0&0&0&a^{\prime}_{2}&0\\ 0&0&0&0&0&a^{\prime}_{3}\\ y^{\prime}_{1}u_{1}&y^{\prime}_{2}u_{1}&y^{\prime}_{3}u_{1}&M&0&0\\ y^{\prime}_{1}u_{2}&\omega y^{\prime}_{2}u_{2}&\omega^{2}y^{\prime}_{3}u_{2}&0&M&0\\ y^{\prime}_{1}u_{3}&\omega^{2}y^{\prime}_{2}u_{3}&\omega y^{\prime}_{3}u_{3}&0&0&M\\ \end{array}\right)$$

$$a_1^{\prime }$$

$$a_2^{\prime }$$

$$a_3^{\prime }$$

$$\displaystyle M_{\nu,N}\,\,=\,\left(\begin{array}[]{cc}0&\mathrm{diag}(a^{\prime}_{1},a^{\prime}_{2},a^{\prime}_{3})\\ \mathrm{diag}(u_{1},u_{2},u_{3})\sqrt{3}U_{\rm m}\mathrm{diag}(y^{\prime}_{1},y^{\prime}_{2},y^{\prime}_{3})&M\mathrm{diag}(1,1,1)\\ \end{array}\right)$$

$$\displaystyle U_{\rm m}=\dfrac{1}{\sqrt{3}}\,\left(\begin{array}[]{ccc}1&1&1\\ 1&\omega&\omega^{2}\\ 1&\omega^{2}&\omega\\ \end{array}\right)\leavevmode\nobreak\ .$$

$$M_{\nu }=\frac{1}{M}\mathrm{diag}(a_1,a_2,a_3)\sqrt{3}{U}_{\mathrm{m}}\mathrm{diag}(y_1^{\prime },y_2^{\prime },y_3^{\prime })\text{leavevmode},$$

$$a_1$$

$$a_2$$

$$a_3$$

$$U_{\nu }^{†}{M}_{\nu }{V}_{\nu }=D_{\nu },$$

$${\theta }_{23}^{\nu }$$

$${\theta }_{12}^{\nu }$$

$$U_{LM}$$

$${\mathcal{L}}_{\mathrm{Yuk},\mathrm{u}}$$

$$+$$

$$M_u$$

$${\mathcal{L}}_{\mathrm{Yuk},\mathrm{d}}$$

$$+$$

$$M_d$$

$${\mathcal{L}}_{\mathrm{Yuk},\mathrm{l}}$$

$$+$$

$$\displaystyle M_{l}=\left(\begin{array}[]{ccc}0&a_{l}\alpha&b_{l}\\ b_{l}\alpha&0&a_{l}r\\ a_{l}&b_{l}r&0\\ \end{array}\right)\leavevmode\nobreak\ .$$

$$\text{Det}S$$

$$\text{Tr}S$$

$$(\text{Tr}S{)}^2-\text{Tr}{S}^2$$

$$v_1^d\gg v_3^d\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\mathrm{and}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\frac{v_1^d}{v_2^d}\gg \frac{y_1^l+y_2^l}{y_1^l-y_2^l}\gg 1,$$

$$r\gg \alpha \text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\mathrm{and}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}\text{leavevmode}r\gg \frac{b_l}{a_l}\gg 1.$$

$$r$$

$$a_l$$

$$b_l$$

$$r$$

Full text

I Introduction

Probably the number one mystery in particle physics is the understanding of the pattern of fermion masses and mixings from first principles. Indeed, the charged fermion mass pattern is not described in the theory: the Standard Model only allows us the freedom to fit the observed charged fermion masses, while lacking the masses of neutrinos altogether. An approach towards addressing, at least partially, the charged fermion mass problem, is the possibility of relating quarks and lepton masses as a result of a flavor symmetry Morisi:2011pt , i.e.

[TABLE]

Notice that this mass relation constitutes a consistent flavor-dependent generalization of the conventional bottom-tau SU(5) prediction, but does not require grand-unification. It provides a partial solution to the charged fermion mass problem, which can be shown to hold in some theories of flavor based on the $A_{4}$ Morisi:2011pt ; King:2013hj ; Morisi:2013eca and $T_{7}$ Bonilla:2014xla symmetries.

Turning to neutrinos, the origin of their mass, the understanding of their mixing properties and the puzzle of whether they are their own anti-particles continue to defy theorists. Underpinning the solution to such neutrino puzzles may not only write a new chapter of particle physics, but also shed light on astrophysical and cosmological puzzles. One of the latter is the puzzle of Dark Matter, believed to be associated to the existence of a new absolutely or nearly stable neutral particle.

There have been attempts at formulating joint solutions to the above shortcomings of the standard model. For example, in scotogenic models dark matter is introduced as a messenger of radiative neutrino mass generation Ma:2006km ; Hirsch:2013ola ; Merle:2016scw whose stability follows from the radiative nature of the neutrino mass. Several alternative ideas have come out, invoking non-Abelian flavor symmetries Hirsch:2012ym ; Morisi:2012fg ; King:2014nza , such as the $A_{4}$ symmetry Ma:2001dn ; Babu:2002dz ; Morisi:2013qna . For example, dark matter could be stable as a result of some remnant of the flavor symmetry associated to the pattern of neutrino mixing Hirsch:2010ru ; Boucenna:2012qb . In all these models neutrinos are Majorana type.

However, recently there has been a renewed interest in Dirac neutrinos Chulia:2016ngi ; Chulia:2016giq ; Aranda:2013gga ; Bonilla:2016diq ; Ma:2016mwh ; Ma:2015raa ; Ma:2015mjd ; Bonilla:2016zef ; Valle:2016kyz ; Abbas:2016qbl ; Abbas:2013uqh ; Wang:2016lve ; Wang:2017mcy ; Okada:2014vla ; Borah:2016zbd ; Borah:2017leo which may attain naturally small masses in many scenarios. For example, it can happen that the same flavor symmetry which presumably sheds light on the pattern of neutrino oscillation parameters, also implies that neutrinos are Dirac fermions Aranda:2013gga . Specially tantalizing is the idea that the stability of dark matter can be directly traced to the Dirac nature of neutrinos Chulia:2016ngi ; Chulia:2016giq ; Bonilla:2016diq ; Ma:2016mwh . One way to realize this idea is by means of a $Z_{4}$ Lepton Quarticity symmetry Chulia:2016ngi ; Chulia:2016giq . Within such approach the same $Z_{4}$ discrete lepton number symmetry ensures the stability of dark matter and the absence of all the Majorana mass terms. Thus owing to Lepton Quarticity, the Dirac nature of neutrinos and the stability of dark matter are intimately related: the breakdown of this symmetry will simultaneously imply loss of dark matter stability as well as the Diracness of neutrinos.

Here we focus on the Lepton Quarticity models of dark matter, along the lines pursued in Chulia:2016ngi ; Chulia:2016giq . The plan of the paper is as follows. In Sect. II we sketch in some detail the extended particle content required to realize the non-Abelian flavor symmetry of the model, and show how the Dirac nature of neutrinos and the smallness of their seesaw–induced masses both follow from our non-Abelian discrete flavor symmetry. We also briefly discuss the appearance of a viable WIMP dark matter candidate in this model. In Sect. III we present our predictions for the current and future neutrino oscillation experiments. We find that the atmospheric angle $\theta_{23}$ and the CP phase $\delta_{CP}$, whose current experimental determination is still rather poor, are tightly related to each other within our model. Finally we summarize our results in Sect. IV.

II The model setup

Here we describe the model in some detail. The particle content of our model along with the $SU(2)_{L}\otimes Z_{4}\otimes A_{4}$ charge assignments of the particles are given in Table 1.

Note that in Table 1 the $L_{i}=(\nu_{i},l_{i})^{T}$, $i=e,\mu,\tau$ denote the lepton doublets, transforming as indicated under the flavor symmetry.

Apart from the Standard Model fermions, the model also includes three right–handed neutrinos $\nu_{i,R}$ which are singlets under the $\mathrm{SU(3)_{C}\otimes SU(2)_{L}\otimes U(1)_{Y}}$ gauge group, singlets under $A_{4}$, but carry charge $\mathbf{z}$ under $Z_{4}$. We also add three gauge singlet Dirac fermions $N_{i,L},N_{i,R}$; $i=1,2,3$ transforming as triplets of $A_{4}$ and with charge $\mathbf{z}$ under $Z_{4}$, as shown in Table 1. Notice that in the scalar sector we have two different sets of fields $\Phi^{u}_{i},\Phi^{d}_{i}$; $i=1,2,3$, which are all doublets under the SU(2)L gauge group, both sets transforming trivially under $Z_{4}$. Under the $A_{4}$ flavor symmetry, $\Phi_{i}^{d}$ transforms as a triplet, while $\Phi_{i}^{u}$ transform as singlets. In addition to the above symmetries we also impose an additional $Z_{2}$ symmetry111This additional $Z_{2}$ symmetry is only required in a non-supersymmetric variant. Clearly the model can be easily supersymmetrized, in which case this additional $Z_{2}$ symmetry is no longer required.. Under this $Z_{2}$ symmetry, all the fields transform as $1$ except for $\Phi^{d}_{i}$, $l_{i,R}$ and $d_{i,R}$, which transform as $-1$. The role of this $Z_{2}$ symmetry is to prevent the Higgs doublets $\Phi_{i}^{d}$ from coupling the up-type quarks and neutrino sector, and the $\Phi^{u}_{i}$ Higgs doublets from the down-type quarks and charged leptons.

In addition we need scalar singlets, for example the $\chi_{i}$, $i=1,2,3$. These are gauge singlets transforming as a triplet under the $A_{4}$ and trivially under $Z_{4}$. We also add two other gauge singlet scalars $\zeta$ and $\eta$ both of which transform trivially under $A_{4}$ but carry $Z_{4}$ charges $\mathbf{z}$ and $\mathbf{z}^{2}$ respectively. Notice that, since under the $Z_{4}$ symmetry the field $\eta$ carries a charge $\mathbf{z}^{2}=-1$, it follows that $\eta$ can be taken to be real.

As discussed in Chulia:2016ngi ; Chulia:2016giq the lepton quarticity symmetry $Z_{4}$ serves a double purpose. It not only ensures that neutrinos are Dirac particles, but also guarantees the stability of the scalar particle $\zeta$, making it a viable dark matter WIMP. If the quarticity symmetry is broken either by an explicit soft term or spontaneously, through non-zero vacuum expectation values (vevs) to any of the scalars $\eta,\zeta$ which carry a non-trivial $Z_{4}$ charge, then both the Dirac nature of neutrinos and stability of dark matter is simultaneously lost.

We now turn our attention to the Yukawa sector of our model. In the neutrino sector the Yukawa terms relevant for generating masses for the neutrinos and the heavy neutral fermions $N_{L},N_{R}$ are given by

[TABLE]

where $3S$ and $3A$ denote the symmetric and antisymmetric $A_{4}$ triplet combinations obtained from the tensor product of two $A_{4}$ triplets. Notice also that $3S$ and $3A$ are not two different irreducible representations of $A_{4}$, which only has one triplet, but simply different contractions with the same transformation rule. Also, $y_{i},y^{\prime}_{i},c_{1},c_{2}$; $i=1,2,3$ are the Yukawa couplings which, for simplicity, are taken to be real. The parameter $M$ is the gauge and flavor-invariant mass term for the heavy leptons. Here we like to highlight the important role played by the $A_{4}$ flavor symmetry. Owing to the $A_{4}$ charges of the left and right handed neutrinos, a tree level Yukawa coupling between them of type $y_{\nu}\,\bar{L}_{L}\nu_{R}\Phi^{u}_{i}$ is forbidden. Thus neutrino masses can only appear through type-I Dirac seesaw mechanism as we now discuss.

After symmetry breaking the scalars $\chi_{i}$ and $\Phi^{u}_{i}$ acquire vevs $\left\langle\chi_{i}\right\rangle=u_{i}$; $\left\langle\Phi^{u}_{i}\right\rangle=v_{i}^{u}$; $i=1,2,3$. The invariant mass term $M$ can be naturally much larger than the symmetry breaking scales, i.e. $M\gg v^{u}_{i},u_{i}$. In this limit, for any numerical purpose the last two terms in Eq. 66 can be safely neglected. Under this approximation the $6\times 6$ mass matrix for the neutrinos and the heavy neutral fermions in the basis $(\bar{\nu}_{e,L},\bar{\nu}_{\mu,L},\bar{\nu}_{\tau,L},\bar{N}_{1,L},\bar{N}_{2,L},\bar{N}_{3,L})$ and $(\nu_{e,R},\nu_{\mu,R},\nu_{\tau,R},N_{1,R},N_{2,R},N_{3,R})^{T}$ is given by

[TABLE]

where $\omega$ is the third root of unity, with $\omega^{3}=1$ and

[TABLE]

As mentioned before, owing to the $A_{4}$ symmetry, a direct coupling between $\nu_{L}$ and $\nu_{R}$ is forbidden, leading to the vanishing of all entries in the upper left quadrant of Eq. 73. The mass matrix in Eq. 73 can be rewritten in a more compact form, as

[TABLE]

Where $U_{\rm m}$ is the usual magic matrix,

[TABLE]

Note that, in the limit $M\gg v^{u}_{i},u_{i}$ the mass matrix in Eq. (73) can be easily block diagonalized by the perturbative seesaw diagonalization method given in Ref. Schechter:1981cv . The resulting $3\times 3$ mass matrix for light neutrinos can be viewed as the Dirac version of the well known type-I seesaw mechanism. The above mass generation mechanism can also be represented diagramatically as shown in Fig. 1.

The $3\times 3$ matrix for the light neutrinos is

[TABLE]

where $a_{i}=a^{\prime}_{i}u_{i}$. We take the alignment $u_{1}=u_{2}=u_{3}=u$ for the vev of the $A_{4}$ triplet scalars $\chi_{i}$ similar to Babu:2002dz ; Ma:2001dn . In this alignment limit of $A_{4}$ triplet scalars we have

[TABLE]

which simplifies the notation, although it does not change the form of the neutrino mass matrix in Eq. (82). Notice that we have not imposed any alignment for the vevs of the $A_{4}$ singlet scalars $\Phi^{u}_{i}$ 222Doing so for the different $A_{4}$ singlet scalar vevs is not very natural. Indeed, unlike the case of $A_{4}$ triplet scalars, a priori the vevs of different $A_{4}$ singlet scalars have no reason to obey any mutual alignment. The light neutrino mass matrix of Eq. 82 with the simplified $a_{i}$ of Eq. 83 can be diagonalized by a bi-unitary transformation

[TABLE]

where $D_{\nu}$ is diagonal, real and positive. Owing to the $A_{4}$ flavor symmetry, the resulting rotation matrix acting on left handed neutrinos $U_{\nu}$ in the standard parametrization (for both hierarchies), leads to $\theta^{\nu}_{23}=\frac{\pi}{4}$ and $\delta^{\nu}=\pm\frac{\pi}{2}$ while the other two angles can be arbitrary. Thus, owing to the $A_{4}$ symmetry, $U_{\nu}$ in standard parameterization leads to following mixing angles

[TABLE]

Similar features of maximal $\theta_{23}$ and $\delta$ have been obtained previously in the context of Majorana neutrinos Babu:2002dz ; Ma:2015pma . Although the angles $\theta^{\nu}_{12}$ and $\theta^{\nu}_{13}$ can take arbitrary values they are strongly correlated with each other.

We have performed an extensive numerical scan for both type of hierarchies in the whole parameter range taking all Yukawa couplings in the perturbative range of $[-1,1]$. We find that in the whole allowed range for either type of hierarchy, one cannot simultaneously fit both $\theta^{\nu}_{12}$ and $\theta^{\nu}_{13}$ in the current global experimental range obtained from neutrino oscillation experiments Forero:2014bxa . This implies that in our model $U_{\nu}$ alone cannot explain the current neutrino oscillation data.

However, the lepton mixing matrix $U_{LM}$ which is probed by neutrino oscillation experiments is the product of the charged lepton rotation matrix $U_{l}$ with the neutrino transformation matrix $U_{\nu}$ Schechter:1980gr i.e.

[TABLE]

In our model the charged lepton mixing matrix $U_{l}$ is also non-trivial and contributes to the full leptonic mixing matrix $U_{LM}$. We now move to discuss the structure of mass matrices and mixing matrices for charged leptons as well as the up and down type quarks.

We now turn to the discussion with up type quark mass matrix. The invariant Yukawa Lagrangian relevant to generating up type quark mass matrix is given by

[TABLE]

where $y^{u}_{i}$; $i=1,2,3$ are the Yukawa couplings which for simplicity we take to be all real. After spontaneous symmetry breaking Eq. 106 leads to a diagonal mass matrix given by

[TABLE]

On the other hand, the Yukawa Lagrangian relevant to down type quarks mass generation is given by

[TABLE]

where $y^{d}_{i}$; $i=1,2$ are the Yukawa couplings which for simplicity are taken to be real. The resulting mass matrix for down type quarks after spontaneous symmetry breaking is given by

[TABLE]

where $\left\langle\Phi^{d}_{i}\right\rangle=v^{d}_{i}$; $i=1,2,3$ and $a_{d}=(y_{1}^{d}-y_{2}^{d})v_{2}^{d}$, $b_{l}=(y_{1}^{d}+y_{2}^{d})v_{2}^{d}$. Moreover, $\alpha$ and $r$ are ratios of the vevs of $\Phi^{d}_{i}$ and are given as $\alpha={}^{v_{3}^{d}}\!/_{v_{2}^{d}}$ and $r={}^{v_{1}^{d}}\!/_{v_{2}^{d}}$.

Finally, the invariant Yukawa terms for the charged leptons is given by

[TABLE]

where $y^{l}_{i}$, $i=1,2$, are the Yukawa couplings which, for simplicity, we take to be real. After symmetry breaking the charged lepton mass matrix is given by

[TABLE]

where, just as in the down quark case, here also $a_{l}=(y_{1}^{l}-y_{2}^{l})v_{2}^{d}$, $b_{l}=(y_{1}^{l}+y_{2}^{l})v_{2}^{d}$. The parameters $\alpha,r$ which are the ratios of the vevs of $\Phi^{d}_{i}$ i.e. $\alpha={}^{v_{3}^{d}}\!/_{v_{2}^{d}}$ and $r={}^{v_{1}^{d}}\!/_{v_{2}^{d}}$ are the same as those defined after Eq. 134. This matrix is completely analogous to the down-type quark mass matrix. Note that while $\alpha$ and $r$ are the same both in the quark and in the lepton sector, as they are simply ratios between the vevs of $\Phi_{i}^{d}$, while $a_{f}$ and $b_{f}$, $f\in\{l,q\}$, are different.

These mass matrices for charged leptons and down-type quarks correspond to those discussed in Morisi:2011pt ; King:2013hj ; Morisi:2013eca ; Bonilla:2014xla and lead to the generalized bottom-tau relation of 1. In section III we show that there is enough freedom to fit the charged lepton and down type quark masses within their 1-$\sigma$ range. Apart from fitting all the masses as well as leading to the generalized bottom-tau relations, the charged lepton mass matrix 158 also leads to non-trivial charged lepton rotation matrix $U_{l}$. As we show in section III this non-trivial contribution from $U_{l}$ results in a lepton mixing matrix $U_{LM}$ consistent with the current global fits to neutrino oscillation data Forero:2014bxa . The lepton mixing matrix obtained from our model also implies normal hierarchy for neutrino masses and leads to an interesting correlation between the atmospheric mixing angle and CP violating phase, the two most ill–determined parameters in leptonic mixing matrix.

III Flavor predictions: numerical results

In this section we discuss the phenomenological implications of our model. The important predictions emerging in our model are: a) the flavor-dependent bottom-tau unification mass relation of Eq. (1) b) a correlation between the two poorly determined oscillation parameters: the atmospheric angle $\theta_{23}$ and $\delta_{CP}$ and c) a normal hierarchy for the neutrinos. In this section we discuss these numerical predictions in some detail, given the experimentally measured “down-type” fermion masses, solar and reactor mixing angles as well as neutrino squared mass differences.

III.1 Charged lepton and down-type quark masses

We start our discussion by looking in more detail at the down type quark and charged lepton mass matrices discussed previously in Eqs. 134 and 158. This structure for the down type quark and charged lepton mass matrices has been previously discussed in several works Morisi:2011pt ; King:2013hj ; Morisi:2013eca ; Bonilla:2014xla . In this section for illustration purpose we first discuss the results obtained in previous works by closely following the approach taken in previous works like in Morisi:2009sc . Subsequently, we will generalize the analysis of previous works and discuss how the same results can be obtained using a more general setup and more detailed considerations.

We start from the charged lepton mass matrix obtained in Eq. 158. The correct charged lepton masses are reproduced if the vevs of the $A_{4}$ triplet fields $\Phi^{d}_{i}$ satisfy the alignment limit $v^{d}(1,\epsilon_{1},\epsilon_{2})$, where $v^{d}\gg\epsilon_{1},\epsilon_{2}$. Then, in similar notation and spirit as in Ref. Morisi:2009sc , we extract the three invariants of the Hermitian matrix $S=M_{l}M_{l}^{\dagger}$: $Det(S)$, $Tr(S)$ and $Tr(S)^{2}-Tr(S^{2})$. We then compute their values in the diagonal basis in terms of the charged lepton masses, $m_{e}$, $m_{\mu}$ and $m_{\tau}$. The equations are

[TABLE]

The expressions in Eq.159 can be readily solved in the vev alignment limit $v^{d}(1,\epsilon_{1},\epsilon_{2})$ discussed before. This amounts to the approximation

[TABLE]

or equivalently

[TABLE]

The solutions for $r$, $a_{l}$ and $b_{l}$ are given as

[TABLE]

Owing to the $A_{4}$ symmetry, the charged lepton mass matrix in Eq. 158 and down type quark mass matrix in Eq. 134 have the same structure. As a result the down quark mass matrix can also be decomposed using equations analogous to Eq. 159. For down type mass matrix of 134 we obtain

[TABLE]

Note that the parameters $\alpha$ and $r$ are common for both the charged lepton sector as well as in the down-type quark sector, as they are simply ratios between vevs of the fields $\Phi^{d}_{i}$. Thus comparing Eqs. 160 and 163 we obtain the following mass relation

[TABLE]

The procedure sketched above can be performed in a more general way by solving the equations numerically. The relevant equations for the case of charged leptons are

[TABLE]

Taking as input parameters the best fit values (at $M_{Z}$ scale) for the charged lepton masses Xing:2011aa and imposing $r>\frac{b}{a}$, there is a one-parameter family of solutions to these equations. These are related to the approximate solution described before. We build the functions $r(\alpha)$, $a_{l}(\alpha)$ and $b_{l}(\alpha)$ taking $\alpha$ as a free parameter. In the correct range for the parameter $\alpha$, the unique solution is found to be near the limit $r\gg\dfrac{b}{a}\gg 1$ and therefore it again leads to the mass relation in Eq. 1. Since the $(\alpha,r,a_{l},b_{l})$ are solutions of Eqs. III.1, the charged lepton masses are fitted exactly to their best-fit values. In order to underpin the relevant solution for down type quark masses, we also need to take into account not only the mass relation 1 and the charged lepton masses, but also the constraints for the experimental measurements (along with renormalization group evolution to $M_{Z}$ scale) of all the down-type quark masses Xing:2011aa . Here, we will impose the rather stringent 1-$\sigma$ bounds333Imposing 1-$\sigma$ is in fact rather stringent, and can easily be relaxed to a more conservative criterium e.g. 3-$\sigma$. We have deliberatively imposed the stringent 1-$\sigma$ bound in other to highlight the high precision obtained from our results. on the down-type quarks masses at $Z$ boson energy scale Xing:2011aa .

Then, for each valid $(\alpha,r,a_{l},b_{l},)$ we take $a_{q}$ and $b_{q}$ as

[TABLE]

where $\epsilon_{1}$ and $\epsilon_{2}$ are expected to be small.

Using this procedure gives sets of parameters $(\alpha,r,a_{l},b_{l},a_{q},b_{q})$ which give at the same time best fit values for charged lepton masses, down-type quarks inside the $1\sigma$ range and the mass relation in Eq. 1. In figure 2 we show the family-dependent bottom-tau mass prediction of our model for the s and d masses, along with their allowed 1-$\sigma$ ranges.

III.2 The charged piece of th lepton mixing matrix $U_{l}$

The charged lepton mass matrix Eq. 158 not only leads to correct lepton masses but also to non-trivial charged lepton rotation matrix $U_{l}$ as we discuss now. Just like the neutrino mass matrix, the charged lepton mass matrix can also be diagonalized by a bi-unitary transformation as

[TABLE]

The charged lepton mixing matrix $U_{l}$ in standard parameterization can be written as

[TABLE]

where $P$ and $P^{\prime}$ are diagonal matrix of phases and $U_{ij}$ is the usual complex rotation matrix appearing in the symmetrical parametrization of fermion mixing given in Schechter:1980gr , e.g.

[TABLE]

with an analogous definitions for $U_{13}$ and $U_{23}$. For the charged lepton mass matrix 158, we find that

[TABLE]

Where $\lambda_{C}\approx 0.22$ is the sine of the Cabbibo angle. In order to reproduce adequate values for the CKM matrix elements we may introduce a vector-like quark mixing with the up-type quarks, along the lines followed recently in Bonilla:2016sgx .

The diagonal phases in P and P’ are all are exactly [math] except for one which is $\pi$. Performing the numerical computation reconfirms the results obtained in Eq. 176 for the charged lepton mass matrix i.e. $\theta^{l}_{12}$ is finite and its value depends on the value of $\alpha$ in an inverse way, while $\theta^{l}_{13}$ and $\theta^{l}_{23}$ are both negligible (in particular, $\theta^{l}_{13}\sim 10^{-5}$ and $\theta^{l}_{23}\sim 10^{-7}$. Then, the charged lepton mixing matrix for our model is given as

[TABLE]

Thus in our model the lepton mixing matrix $U_{LM}=U^{\dagger}_{l}U_{\nu}$ receives significant charged lepton corrections which have interesting phenomenological consequences as we discuss in next section.

III.3 The lepton mixing matrix and neutrino mass ordering

As mentioned before in Section II, the light neutrino mass matrix in Eq. 82 leads to the neutrino mixing matrix $U_{\nu}$ which in standard parameterization Schechter:1980gr leads to

[TABLE]

As mentioned before, owing to the $A_{4}$ symmetry, we have that $\theta_{23}^{\nu}=45^{\circ}$ and $\delta^{\nu}_{CP}=90^{\circ}$ for both types of mass ordering: normal hierarchy (NH) or inverted hierarchy (IH). Since neutrinos in our model are Dirac fermions, the phases in the right in Eq. 178 i.e. $P^{\prime}$, are unphysical, while $\theta_{13}^{\nu}$ and $\theta_{12}^{\nu}$ are strongly correlated between each other. This result is completely general and follows from the $A_{4}$ symmetry, independently of the mass hierarchy, NH or IH. However, the behavior of the correlation between $\theta_{12}^{\nu}$ and $\theta_{13}^{\nu}$ does depend on the choice of NH or IH.

Taking into account the results in the previous sections, the lepton mixing matrix is

[TABLE]

One can regard the matrix $U_{l}$ as a correction to the neutrino mixing parameters obtained just by diagonalizing the neutrino mass matrix. For the NH case, the angle $\theta_{12}^{l}$ has to be big enough ($\sim>15^{\circ}$) so as to account for the correct mixing angles of the lepton mixing matrix, but at the same time it has to remain controlled ($<20^{\circ}$) otherwise the down-type quark masses cannot be fitted. This means that the parameter $\alpha$ has to be between $0.04$ and $0.08$. This lepton mixing matrix can fit the neutrino oscillation parameters within $3\sigma$ at the same time as the mass matrices fit the down-type quarks and the neutrino squared mass differences in the $1\sigma$ range and the charged lepton masses exactly. Once the lepton mixing matrix is written in the standard parametrization, two interesting features arise. On the one hand, $\theta_{23}>45^{\circ}$ and, on the other, a strong correlation appears between the atmospheric angle $\theta_{23}$ and $\delta_{CP}$, as shown in figure 3.

For IH, a different scenario arises. As in the case for NH, lepton corrections cannot be very big otherwise the down-type quark masses will not be fitted. However, the structure of the correlation between $\theta_{12}^{\nu}$ and $\theta_{13}^{\nu}$ implies that for allowed charged lepton corrections, the reactor angle $\theta_{13}$ is always outside the $3\sigma$ allowed range. Note that the model does not include any a priori theoretical bias in favour of normal hierarchy but it is a prediction of the model once one impose experimental constraints.

IV Discussion and Summary

We have proposed a $A_{4}\otimes Z_{4}\otimes Z_{2}$ flavor extension of the Standard Model with naturally small Dirac neutrino masses. Our lepton quarticity symmetry simultaneously forbids Majorana mass terms and provides dark matter stability. The flavor symmetry plays a multiple role, providing : (i) a generalized family-dependent bottom-tau mass relation, Eq. (1) and Fig. 2, (ii) a natural realization of the type-I seesaw mechanism for Dirac neutrino masses, as the tree level Dirac Yukawa term between left and right handed neutrinos is forbidden, (iii) a very predictive flavor structure to the lepton mixing matrix. The latter directly correlates the CP phase $\delta_{CP}$ and the atmospheric angle $\theta_{23}$, as shown in Fig. 3. This implies that (iv) CP must be significantly violated in neutrino oscillations, and the atmospheric angle $\theta_{23}$ lies in the second octant, (v) only the normal neutrino mass ordering is realized.

Our approach provides an adequate pattern of neutrino mass and mixing as well as a viable stable dark matter. This is achieved while providing testable predictions concerning the currently most relevant oscillation parameters, the atmospheric angle $\theta_{23}$ and the CP phase $\delta_{CP}$, as well as a successful family generalization of bottom-tau unification, despite the absence of an underlying Grand Unified Theory. Our lepton quarticity approach also leads to other interesting phenomena such as neutrinoless quadruple beta decay ($0\nu 4\beta$), which has recently been probed by the NEMO collaboration Arnold:2017bnh . The intimate connection between the Dirac nature of neutrinos and dark matter stability constitutes a key feature of our model. Other phenomenological implications will be taken up elsewhere.

Acknowledgements.

This research is supported by the Spanish grants FPA2014-58183-P, Multidark CSD2009-00064, SEV-2014-0398 (MINECO) and PROMETEOII/2014/084 (Generalitat Valenciana). The feynman diagram in Fig.1 was drawn using Jaxodraw Binosi:2003yf .