Deviations to Tri-Bi-Maximal mixing in the limit of $\mu-\tau$ symmetry

TL;DR

This paper investigates how small deviations from $$ symmetry affect Tri-Bi-Maximal neutrino mixing, predicting observable effects on neutrino parameters and neutrinoless double beta decay, especially for quasi-degenerate neutrino masses.

Contribution

It introduces two ansatzes for corrected mixing patterns under $$ symmetry, linking Majorana $CP$ phases to measurable neutrino parameters and decay observables.

Findings

Constraints on Majorana $CP$ phases from mixing parameters.

Predictions for neutrinoless double beta decay in quasi-degenerate neutrino mass regime.

Testable deviations in upcoming neutrino experiments.

Abstract

In the limit of an approximate $\mu-\tau$ symmetry in the neutrino mass matrix, we explore deviations to the Tri-Bi-Maximal mixing pattern in the neutrino sector. We consider two different ansatzes for the corrected pattern to predict the current values of neutrino mixing parameters. We show that it is possible to constrain the Majorana $CP$ phases by studying their correlation to the mixing parameters and we study their effects on neutrinoless double beta decay observables. These predictions are sharp for the quasi-degenerate ordering and can be tested in upcoming experiments.

Figures

Equations

$$U_{\mathrm{PMNS}}$$

$$U_{\mu-\tau}=\left(\begin{array}[]{ccc}c_{12}&s_{12}&0\\ \frac{-s_{12}}{\sqrt{2}}&\frac{c_{12}}{\sqrt{2}}&\frac{-1}{\sqrt{2}}\\ \frac{-s_{12}}{\sqrt{2}}&\frac{c_{12}}{\sqrt{2}}&\frac{1}{\sqrt{2}}\end{array}\right),$$

$$U_{\mathrm{PMNS}}=U_{TBM}{U}_{\mathrm{Corr}}$$

$$U_{\mathrm{Corr}}=U_{ij}(\varphi ,\sigma )\mathrm{Diag}\left(1,{\mathrm{e}}^{-\mathrm{i}\frac{{\alpha }_1}{2}},{\mathrm{e}}^{-\mathrm{i}\frac{{\alpha }_2}{2}}\right),$$

$$U_{13}=\left(\begin{array}[]{ccc}\cos\phi&0&\sin\phi~{}e^{-i\sigma}\\ 0&1&0\\ -\sin\phi~{}e^{i\sigma}&0&\cos\phi\end{array}\right),$$

$$U_{23}=\left(\begin{array}[]{ccc}1&0&0\\ 0&\cos\phi&\sin\phi~{}e^{-i\sigma}\\ 0&-\sin\phi~{}e^{i\sigma}&\cos\phi\end{array}\right),$$

$${\sin }^2{\theta }_{12}=\frac{\mid U_{e2}{\mid }^2}{1-\mid U_{e3}{\mid }^2},{\sin }^2{\theta }_{23}=\frac{\mid U_{\mu 3}{\mid }^2}{1-\mid U_{e3}{\mid }^2},{\sin }^2{\theta }_{13}=\mid U_{e3}{\mid }^2,$$

$$J_{CP}$$

$${\beta }_1={\alpha }_1,{\beta }_2={\alpha }_2+2(\sigma -{\delta }_{CP}).$$

$${\sin }^2{\theta }_{12}$$

$${\sin }^2{\theta }_{23}$$

$${\sin }^2{\theta }_{13}$$

$$J_{CP}=-\frac{1}{6\sqrt{3}}\sin 2\varphi \sin \sigma ,\sin {\delta }_{CP}=-\frac{(2+\cos 2\varphi )\sin \sigma }{{\left[(2+\cos 2\varphi {)}^2-3{\sin }^{2}2\varphi {\cos }^2\sigma \right]}^{1/2}},$$

$${\sin }^2{\theta }_{12}$$

$${\sin }^2{\theta }_{23}$$

$${\sin }^2{\theta }_{13}$$

$$J_{CP}=-\frac{1}{6\sqrt{6}}\sin 2\varphi \sin \sigma ,\sin {\delta }_{CP}=-\frac{(5+\cos 2\varphi )\sin \sigma }{{\left[(5+\cos 2\varphi {)}^2-24{\sin }^{2}2\varphi {\cos }^2\sigma \right]}^{1/2}}.$$

$$|{\sin }^2{\theta }_{23}-\frac{1}{2}|$$

$$M_{\nu }=M_{\mu -\tau }+\delta M\left(\hat{\delta },\widehat{ϵ}\right).$$

$$\hat{\delta }$$

$$\widehat{ϵ}$$

$$|m_2|$$

$$|m_1|$$

$$0.273<{\sin }^2{\theta }_{12}<0.379,0.0199<{\sin }^2{\theta }_{13}<0.0244,0.453<{\sin }^2{\theta }_{23}<0.598$$

Full text

I Introduction

Neutrino mixing angles have been determined with unprecedent precision in recent years Tanabashi:2018oca ; Capozzi:2013csa ; deSalas:2017kay . These angles define the structure of the lepton mixing matrix kwown as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix Pontecorvo:1957cp ; Maki:1962mu , which can be written in the standard form

[TABLE]

Here, $c_{ij}$ and $s_{ij}$ stand for $\cos\theta_{ij}$ and $\sin\theta_{ij}$, respectively, and $\theta_{ij}$ denotes the mixing angles $\theta_{12}$, $\theta_{13}$, and $\theta_{23}$ . $\delta_{CP}$ is the Dirac-$CP$ violating phase, whereas $\beta_{1}$ and $\beta_{2}$ are two additional phases which account for the Majorana nature of neutrinos. Despite the success in determining the mixing angles Tanabashi:2018oca ; Capozzi:2013csa ; deSalas:2017kay , a precise determination of the $CP$ violating phases is still missing. The next-to-next generation of neutrino experiments could finally help to determine the Dirac $CP$ violating phase ($CPVP$) Abi:2018dnh . However, given the lack of experimental observables directly related to the Majorana phases, the needed of indirect, but complementary determinations of their values becomes an interesting task to be explored.

Oscillation experiments confirm that the reactor $\theta_{13}$ angle is small but not zero, whereas the atmospheric $\theta_{23}$ angle its close to its maximal value, $\pi/4$. A direct consequence of choosing these critical values ($\theta_{13}=0$ and $\theta_{23}=\pi/4$) is that the mixing matrix in Eq. (5), without including Majorana pases, takes the form

[TABLE]

with the mixing angle $\theta_{12}$ as the only free parameter. The subindex $\mu-\tau$ in Eq. (6) refers to the so-called $\mu-\tau$ symmetry Fukuyama:1997ky ; Harrison:2002et ; Babu:2002ex ; Ohlsson:2002rb ; Mohapatra:2004mf ; Ghosal:2004qb ; Ma:2004zv ; Mohapatra:2006un ; Joshipura:2005vy ; Fuki:2006xw ; Riazuddin:2007aa ; Luhn:2007sy ; Koide:2008zza ; Honda:2008rs ; Ishimori:2008gp ; Ge:2010js ; Ge:2011ih ; He:2011kn ; Ge:2011qn ; Talbert:2014bda ; He:2012yt since it satisfies the relations $|U_{\mu i}|=|U_{\tau i}|$. Among the symmetry approaches based on the $\mu-\tau$ symmetry, it has been of great interest the Tri-Bi-Maximal (TBM) mixing pattern Vissani:1997pa where $\sin^{2}\theta_{12}=1/3$. The TBM pattern has been the starting point of many theoretical works since it can be generated from larger flavor symmetries (see for instance Altarelli:2012ss ; Xing:2015fdg and references therein). Nonetheless, the predicted mixing angles within this symmetric approach are in conflict with current experimental determinations.

Deviations to the TBM scenario have been investigated in order to restore the compatibility with latest neutrino data Petcov:2014laa ; Garg:2013xwa ; Garg:2017mjk ; Shimizu:2014ria ; Sruthilaya2015 ; Sruthilaya2016 . Some parametrizations written in the form

[TABLE]

have been considered, where $U_{\rm Corr}$ is a correction matrix which encodes the deviation from the TBM mixing matrix $U_{TBM}$. This correction is usually written in terms of up to two rotation angles in the orthogonal case Garg:2013xwa ; Garg:2017mjk , and including one complex phase in the unitary case Shimizu:2014ria ; Sruthilaya2015 ; Sruthilaya2016 . These approaches have been able to reproduce the pattern of mixing parameters in the limit of small rotation angles and, in some cases, have provided some predictions for the Dirac $CP$ phase. However, this type of ansatz could not give any hint about the Majorana phases. One possibility of exploring the Majorana case is to include additional phases into the TBM corrected matrix as in Petcov:2014laa ; Chen:2018eou ; Chen:2018zbq .

On the other hand, the possibility of still having an approximate $\mu-\tau$ symmetry in the neutrino mass matrix has been explored in Gupta:2013it , where the viability of such scenario has been confronted with neutrino data. In particular, the connection between an approximate $\mu-\tau$ symmetry in the neutrino mass matrix and both Dirac and Majorana phases have been also explored without considering a particular ansatz of the mixing matrix Rivera2 .

In this paper, we investigate different ansatz which correct the TBM mixing and study their effect on the neutrino mass matrix, which, to the best of our knowledge, has not been explored. The novelty of our approach is that, by defining an approximate $\mu-\tau$ symmetry in the mass matrix, the Majorana phases could be bounded and affect neutrino observables as the neutrinoless double beta decay amplitude ($|m_{ee}|$). The remainder of this work is organized as follows. First, we present the different scenarios which correct the TBM pattern. Then, we describe the $\mu-\tau$ symmetric limit of the neutrino mass matrix and its connection with the correction parameters. Third, we show the main results of our analysis and the phenomenological implications over the $CP$ parameters and $|m_{ee}|$. Finally, we give our final comments and conclusions.

II Deviations to TBM pattern

For the sake of simplicity, we will consider in our forthcoming analysis only deviations to the TBM pattern coming from the neutrino sector, such that the corrected mixing matrix takes the form of Eq. (7). Let us now consider the following ansatz for the correction matrix

[TABLE]

where $U_{ij}$ is a unitary matrix which depends on the rotation angle $\phi$, and the complex phase $\sigma$. Here $i,j=1,2,3$, and $i\neq j$. The last term in Eq. (8) is a diagonal matrix with two complex phases, $\alpha_{1}$ and $\alpha_{2}$. This approach incorporates in a similar fashion to the standard parametrization in Eq. (5) two new complex phases, which are expected to be related to the physical Majorana $CP$ phases. These type of parametrizations, which can be generated in a model dependent way from specific flavor symmetries, has been discussed elsewhere Petcov:2014laa .

Among the three possibilities in writing $U_{ij}$, it is direct to see from Eqs. (7) and (8) that only the rotations

[TABLE]

and

[TABLE]

lead to a mixing matrix able to accommodate $\theta_{13}\neq 0$, as we will show later. We should notice that for recovering the TBM matrix it is enough to consider the limit where $\phi=\alpha_{1}=\alpha_{2}=0$. Clearly, in this case, $U_{\rm Corr}$ corresponds to the identity matrix as it is expected. In contrast, however, a $U_{12}$ rotation directly leads to a complex $\mu-\tau$ symmetric matrix predicting $\theta_{13}=0$ and $\theta_{23}=\pi/4$, with an unphysical Dirac phase. For these reasons, we will not consider the rotation $1-2$ in our following discussion444In fact, a rotation of the $1-2$ neutrino sector may come from higher order perturbative corrections, which introduces a correction to the zeroth order TBM prediction $\sin^{2}\theta_{12}=1/3$. In such case, the predicted mixing matrix still preserves a $\mu-\tau$ symmetric structure..

We can see from the standard parametrization of the PMNS matrix in Eq. (5) that the mixing angles can be defined in terms of the elements of the neutrino mixing matrix Tanabashi:2018oca :

[TABLE]

while the $CP$ violation parameters can be obtained from

[TABLE]

By comparing Eq. (5) with the corrected $\mu-\tau$ matrix of Eq. (7), we obtain that the Majorana $CP$ phases are related to the new phases via

[TABLE]

As we can see, the relations in Eqs. (11-13) hold independently of the rotation adopted in $U_{ij}$, but may differ when they are written in terms of the correction parameters involved in Eq. (8).

As a first case, let us consider the $1-3$ rotation (Case I) given in Eq. (9). From Eq. (11), the experimental mixing angles are then linked to the correction parameters in Eq. (7) by means of Shimizu:2014ria ; Sruthilaya2015

[TABLE]

The Jarlskog invariant and the Dirac $CP$ can be obtained from

[TABLE]

and the Majorana phases from Eq. (13).

On the other hand, in the second case we take the $2-3$ rotation (Case II) of Eq. (10). In this case, the mixing parameters are given by Shimizu:2014ria ; Sruthilaya2015

[TABLE]

while

[TABLE]

As in the previous case, the Majorana phases can be obtained from Eq. (13).

A direct inspection of Eqs. (II) and (II) shows that current determinations of the $\theta_{13}$ mixing angle forbid null values of the $\phi$ angle, but are in favor of small values of this parameter. Moreover, this parameter allows determining the size of the departure from the maximal value of the atmospheric via the reactor angle, where we obtain the approximate sum rules

[TABLE]

Thus, this shows that deviations of $\theta_{23}$ from the maximal mixing are correlated to deviations in $\theta_{13}$ from zero through the phase $\sigma$, which also correlates these mixings with the $CP$ phases, $\delta_{CP}$ and $\beta_{2}$, as we can see from Eqs. (13), (15) and (17). From the theoretical point of view, this could be of great interest since it point towards the possibility of explaining the observed mixings with a common physical origin in favor of a well-defined flavor symmetry Rivera1 ; Rivera2 .

III $\mu-\tau$ symmetryc limit in the mass matrix

From our previous discussion, we can see that it is mandatory to adopt departures from the TBM mixing pattern in order to explain the observed mixings. Nevertheless, it is still missing a description of the effects of such deviations in the neutrino mass matrix. Before doing this, let us first discuss the $\mu-\tau$ symmetry in the mass matrix.

In the basis where the charged lepton are diagonal, the neutrino mass matrix is obtained from $M_{\nu}=U_{\nu}{\rm diag}(m_{1},m_{2},m_{3})U_{\nu}^{\rm T}$, where $U_{\nu}\equiv U_{\rm PMNS}$. In the $\mu-\tau$ symmetric limit we can write $U_{\nu}=U_{\mu-\tau}$ (as in Eq. (6)). In this case, the resulting mass matrix reflects an exchange symmetry between the $\mu$ and $\tau$ entries, $|m_{e\mu}|=|m_{e\tau}|~{}~{}\text{and}~{}~{}m_{\mu\mu}=m_{\tau\tau}$. It is straightforward to show that the corrections adopted in the mixing matrix will change in consequence this symmetric structure of the neutrino mass matrix.

Following the approach in Rivera1 ; Gupta:2013it , the breaking of the $\mu-\tau$ symmetry in the neutrino mass matrix can be accommodated in two parameters which encode the strength of such breaking. In this way, the mass matrix takes the form

[TABLE]

Here, the matrix $M_{\mu-\tau}$ does posses the characteristics of a $\mu-\tau$ symmetric matrix, whereas $\delta M$ is defined by only two nonzero breaking parameters, $\hat{\delta}$ and $\hat{\epsilon}$. In terms of the matrix elements, these parameters are given by

[TABLE]

It is worth noticing that an approximate $\mu-\tau$ symmetric mass matrix Gupta:2013it is obtained when $|\hat{\delta}|,|\hat{\epsilon}|\ll 1$.

It is direct to show that the breaking parameters in Eq. (III) can be written in terms of the correction parameters ($\phi,~{}\sigma$, $\alpha_{1}$, and $\alpha_{2}$) by using Eqs. (7) and (8). The complete expressions of the breaking parameters depend on the parametrization adopted in the mixing matrix and are rather large. Also, they depend on the mass ordering selected and will not be displayed here. The absolute masses $|m_{1,2,3}|$ can be expressed in terms of the lightest neutrino mass $m_{0}$ as

[TABLE]

Here, the square mass difference $\Delta m_{21}^{2}=m_{2}^{2}-m_{1}^{2}$ is also known as the solar mass scale, and $\Delta m_{31}^{2}=m_{3}^{2}-m_{1}^{2}$ as the atmospheric mass scale, where $m_{0}$ becomes $|m_{1}|$ for the normal mass hierarchy (NH), and $m_{3}$ for the inverted mass hierarchy (IH). Thus, the expressions in Eq. (III) give us a direct relation between the correction parameters of the mixing matrix and the breaking parameters of the mass matrix. In addition, small departures from the $\mu-\tau$ symmetric limit in the mass matrix restrict the parameter space of $|\hat{\delta}|$ and $|\hat{\epsilon}|$ to small values. As a consequence, it could be possible to restrict even more the values of the correction parameters and hence the predicted neutrino mixings.

IV Results

For the numerical analysis, we use the results of the latest global fit for the various neutrino oscillation experiments deSalas:2017kay . For the sake of simplicity, we will work hereafter in the inverted hierarchy (IH), however, as we will show later, our results are valid for both approaches. The $3\sigma$ intervals of the mixing angles are

[TABLE]

and the squared mass differences $\Delta m_{21}^{2}=7.55^{+0.20}_{-0.16}\times 10^{-5}$ eV2 and $|\Delta m_{31}^{2}|=2.42^{+0.03}_{-0.04}\times 10^{-}3$ eV2. Let us divide our discussion into the two cases of interest. First, we analyse the $1-3$ rotation, then, the $2-3$ rotation case.

IV.1 Case I: 1-3 rotation

Based on Eqs. (13), (II), and (15), we will analyse the impact of the correction parameters ($\phi,\sigma$ and $\alpha_{1,2}$) over the mixing parameters when an approximate $\mu-\tau$ mass matrix is required by imposing $|\hat{\delta}|,|\hat{\epsilon}|<0.25$. We will also analyse the case when this condition is omitted. In any case, we will select those values of the correction parameters which predicts mixing angles within the allowed $3\sigma$ region of the experimental values, and show several relations among the predicted mixings. $\phi$ will be scattered in the ($0,\pi/2$) interval, and the phases $\sigma$, $\alpha_{1}$, and $\alpha_{2}$ within the ($-\pi,\pi$) range.

Some particular cases can be analysed by selecting specific values of $\alpha_{1}$ and $\alpha_{2}$. Let us first discuss the case $\alpha_{2}=0$. From Eq. (13), this limit implies that $\beta_{1}=\alpha_{1}$ and $\beta_{2}=2(\sigma-\delta_{CP})$. In Figs. 1 and 2, we can observe that the correlations between the $\theta_{23}$ mixing angle and the $CP$ parameters are narrow when both conditions, small symmetry breaking and $\alpha_{2}=0$, are imposed (blue region). In this case, the deviation of $\theta_{23}$ from its maximal value ($\theta_{23}=\pi/4$) is rather large, i.e. the predicted value lies in a very narrow range near the 3$\sigma$ current limit. A most precise experimental determination of such observable could, in principle, rule out this scenario in the near future. The allowed $3\sigma$ region of $\delta_{CP}$ deSalas:2017kay is also shown in Fig. 1, which shows that this scenario is also disfavored by the last indications of maximal $CP$ violation ($\delta=-\pi/2$) Tanabashi:2018oca ; Capozzi:2013csa ; deSalas:2017kay . However, as we can see from Figs. 1 and 2, there is a strong correlation between $\beta_{2}$ and $\delta_{CP}$ with the atmospheric angle even when the small breaking condition is omitted (gray region), as expected from Eq. (13). Thus, in this approach it could be possible to extract some indirect information about $CP$ violating phases through a precise determination of the atmospheric angle. Moreover, we can identify that the restricted blue region of $\delta_{CP}$ in Fig. 1 is due only to the small symmetry breaking conditions since $\delta_{CP}$ is independent of the correction phases $\alpha_{1,2}$, as we can verify from Eq. (15).

Now, we consider the case $\alpha_{1}=0$. From Eq. (13), it is straightforward to show that this case directly implies that $\beta_{1}=0$, and that $\beta_{2}$ is unbounded since $\alpha_{2}$ remains as a free parameter. Without considering small symmetry breaking requirements, this case does not give useful information about Majorana phases. Moreover, this case is totally ruled out when the requirement of a small symmetry breaking is imposed.

In the trivial case $\alpha_{1}=\alpha_{2}=0$, we also obtain $\beta_{1}=0$, while the correlation of $\delta_{CP}$ ($\beta_{2}$) with $\theta_{23}$ is given again by the gray region in Fig. 1 (2) when the symmetric limit is not included. Again, this case is completely disfavored when the requirement of a small symmetry breaking is added.

Finally, when we leave $\alpha_{1,2}$ as free parameters it is direct to show that $\beta_{1}$ remains free. The predicted region of $\delta_{CP}$ in this scenario is the same for both cases, when we include or omit the approximate symmetry condition, and coincides with the gray region in Fig. 1. The correlation between $\beta_{2}$ and $\sin^{2}\theta_{23}$ is shown in Fig. 3. Here, we can observe that the parameter space of $\beta_{2}$ is reduced when the small braking conditions are imposed (blue region). It is worth stressing that this scenario is the most favored by the experimental data when the small breaking requirements are included.

Concerning to our Case I, we should stress that we have identified two possible scenarios predicting nice correlations between $CP$ violating parameters and the atmospheric mixing angle: the case $\alpha_{2}=0$, and the general case where $\alpha_{1}$ and $\alpha_{2}$ remain free. Through these correlations, it would be possible to extract valuable information about $CP$ violation in the lepton sector in the advent of precise determinations of the mixing angles in forthcoming neutrino experiments.

Now, let us discuss some phenomenological implications of the different scenarios over the effective mass $|m_{ee}|$ describing neutrinoless double beta decay. In Fig. 4, we show the predicted regions of $|m_{ee}|$ in the case of a 1-3 rotation for the inverted hierarchy. We plot separately the cases where $\alpha_{1,2}\neq 0$ (right) and $\alpha_{2}=0$ (left). In each case, we can observe how the region is reduced when the approximate symmetry requirement is included. As we have previously discussed, these are the only cases where at least one of the Majorana phases can be bounded. For comparison, we also plot the allowed region of the mass element $|m_{ee}|$ in the general case when both of the Majorana phases are free (yellow region).

For $\alpha_{2}=0$ (left plot of Fig. 4), and without including restrictions on the breaking parameters (gray boxes), the predicted region of $|m_{ee}|$ is slightly narrow compared to the case of totally free Majorana phases (yellow band). Owing to the correlation between $\beta_{2}$ and the atmospheric angle, this region could be narrowed even more with a more precise determination of $\theta_{23}$, as can be seen from Fig. 2. The implementation of small breakings in the mass matrix restricts the values of the lightest neutrino mass to lie above the 0.03 MeV, showing a marked preference for the quasi-degenerate hierarchy as was noted in Gupta:2013it ; Rivera1 ; Rivera2 . This scenario is disfavored given its preference for maximal deviations in the atmospheric angle.

On the other hand, in the general case $\alpha_{1},\alpha_{2}\neq 0$, and when the small breaking condition is adopted (blue points), a strong preference for the quasi-degenerate hierarchy is also obtained (right plot of Fig. 4). Hence, such scenario could be tested with improved limits of cosmology and/or neutrinoless double beta decay experiments. When the symmetry breaking condition is removed, the predicted region (gray boxes) is extended over the full range of $m_{0}$, as it is expected since both of the Majorana phases remains unbounded in this case.

IV.2 Case II: 2-3 rotation

Let us now move our attention to the case of a 2-3 rotation matrix. The selection criteria is similar to the one used in the previous case by allowing the predicted mixing angles to lie within the 3$\sigma$ experimental range. In Fig. 5 we show the predicted regions of $\delta_{CP}$ versus $\sin^{2}\theta_{23}$ when we include (blue region) and omit (gray region) the requirements for small symmetry breaking. It should be noticed that blue and gray regions are overlapped, which means that the breaking parameters does not affect the parameters $\phi$ and $\sigma$ involved in predicting $\delta_{CP}$ (see Eq. 17).

In Fig. 6, we can see that the $\alpha_{2}=0$ scenario predicts well-defined regions of both Majorana $CP$ phases when the condition of an approximate symmetry is included (blue region). The omission of this last condition leaves $\beta_{1}$ undetermined (gray region in left plot) while $\beta_{2}$ remains unchanged (right plot), as it is expected due to the fact that $\beta_{2}$ does only depend on $\phi$ and $\sigma$ in this scenario.

Other combinations of the phases $\alpha_{1}$ and $\alpha_{2}$ leave Majorana phases undetermined, or fixed to zero, which give no relevant information even if we ask for small symmetry breakings in the neutrino mass matrix.

The predicted region of $|m_{ee}|$ is shown in Fig. 7, which corresponds to the case $\alpha_{2}=0$. The scattered gray boxes (blue points) show the predicted values of $|m_{ee}|$ when the small symmetry breaking condition is omitted (included). As we have discussed, the case $\alpha_{2}=0$ with the small breaking requirement is the only one where both Majorana phases can be bounded. This is reflected in a narrow region of $|m_{ee}|$ (blue points in Fig. 7), which is disfavored by the current cosmological limit of $m_{0}$.

V Summary

It is well known that the TBM mixing pattern predicts maximal value for the atmospheric angle ($\theta_{23}=\pi/4$) and null $\theta_{13}$. The latter value is disfavored by recent experimental data while the former is still allowed. Despite the confirmation of a sizable reactor angle, the TBM pattern remains as a simple alternative to explain neutrino mixings if some modifications are adopted.

In this study, we have parametrized the neutrino mixing matrix by considering deviations from the TBM pattern through a correction matrix in the neutrino sector, which we wrote in terms of four correction parameters. Given its comparison with the standard parametrization, we obtained analytic expressions for the mixing angles, the Dirac $\delta_{CP}$, and the Majorana phases, in terms of these parameters. We have analyzed, in addition, two breaking parameters which helped to define an approximate $\mu-\tau$ symmetry in the mass matrix. These breaking parameters could also be written in terms of the correction parameters and played an important role in bounding the correction phases. In total, we have considered five restrictions (three from the experimental values of the mixing angles, and two from the symmetric limit in the mass matrix) in order to bound the correction parameters, and hence to predict the allowed regions for Majorana and Dirac $CP$ phases in this context.

Concerning our scheme, the predicted $CP$ phases could be related with the atmospheric angle in different scenarios. We have also analyzed the phenomenological implications of the predicted Majorana phases over the effective mass $|m_{ee}|$ describing neutrinoless double beta decay. In this line of thought, our results showed marked preference for the quasi-degenerate hierarchy when the approximate $\mu-\tau$ symmetry requirement was included. Future improvements in the determination of the mixing angles and the Dirac $CP$ phase, joint to the effective mass $|m_{ee}|$, will play a fundamental role in testing such scenarios.

Acknowledgements

SLT acknowledge Universidad Santiago de Cali for the hospitality.