Neutrinoless double beta decay and lepton flavour violation in broken $\mu-\tau$ symmetric neutrino mass models

TL;DR

This paper investigates neutrinoless double beta decay and lepton flavor violation within broken $$ symmetric neutrino mass models in a left-right symmetric framework, analyzing different mixing schemes and their experimental implications.

Contribution

It introduces a comprehensive analysis of neutrinoless double beta decay and lepton flavor violation in various $$ symmetric neutrino models, considering new physics effects at TeV scales.

Findings

New physics contributions can significantly affect neutrinoless double beta decay half-life.

Correlations between lepton flavor violation processes and neutrino parameters are established.

Predictions are compared with current experimental limits, constraining model parameters.

Abstract

We have studied neutrinoless double beta decay and charged lepton flavour violation in broken $\mu-\tau$ symmetric neutrino masses in a generic left-right symmetric model (LRSM). The leading order $\mu-\tau$ symmetric mass matrix originates from the type I (II) seesaw mechanism, whereas the perturbations to $\mu-\tau$ symmetry in order for generation of non-zero reactor mixing angle $\theta_{13}$, as required by latest neutrino oscillation data, originates from the type II (I) seesaw mechanism. In our work, we considered four different realizations of $\mu-\tau$ symmetry, viz. Tribimaximal Mixing (TBM), Bimaximal Mixing (BM), Hexagonal Mixing (HM) and Golden Ratio Mixing (GRM). We then studied the new physics contributions to neutrinoless double beta decay (NDBD) ignoring the left-right gauge boson mixing and the heavy-light neutrino mixing within the framework of LRSM. We have considered the mass of the gauge bosons and scalars to be around TeV and studied the effects of the new physics contributions on the effective mass and the NDBD half life and compared with the current experimental limit imposed by KamLAND-Zen. We further extended our analysis by correlating the lepton flavour violation of the decay processes, $\left(\mu\rightarrow 3e\right)$ and $\left(\mu\rightarrow e\gamma\right)$ with the lightest neutrino mass and atmospheric mixing angle $\theta_{23}$ respectively.

Figures

Tables

Table 1. Table 1: Global fit 3 σ 𝜎 \sigma values of ν 𝜈 \nu oscillation parameters sigma

PARAMETERS	$3\sigma$ RANGES	BEST FIT$\pm 1\sigma$
$\Delta m_{21}^2\left[{10}^{-}5e{V}^2\right]$	7.11-8.18	${7.60}_{-0.18}^{+0.19}$
$\Delta m_{31}^2\left[{10}^{-}3e{V}^2\right]$(NH)	2.30-2.65	${2.48}_{-0.07}^{+0.05}$
$\Delta m_{23}^2\left[{10}^{-}3e{V}^2\right]$(IH)	2.26-2.48	${2.38}_{-0.06}^{+0.05}$
${\sin }^2{\theta }_{12}$	0.278-0.375	0.323$\pm 0.016$
${\sin }^2{\theta }_{23}$(NH)	0.392-0.643	${0.567}_{-0.128}^{+0.032}$
(IH)	$0.403-0.640$	${0.573}_{-0.043}^{+0.025}$
${\sin }^2{\theta }_{13}$(NH)	0.177-0.294	$0.234\pm 0.020$
(IH)	$0.183-0.297$	0.240$\pm$ 0.019
$\delta$	0-2$\pi$(NH)	${254}^0$
	0-2$\pi$(IH)	${266}^0$

Table 2. Table 2: The values of G 0 ​ ν ​ ( Q , Z ) superscript G 0 𝜈 Q Z \rm G^{0\nu}(Q,Z) , Q of the initial isotope for several NDBD processes of experimental interest.

$\beta \beta -decay$	$G^{0\nu }[{10}^{-14}{y}^{-1}$	Q$\left[KeV\right]$	Experiments
${48}_{Ca}\to {48}_{Ti}$	6.3	4273.7	CANDLES
${76}_{Ge}\to {76}_{Se}$	0.63	2039.1	GERDA, Majorana
${82}_{Se}\to {82}_{Kr}$	2.7	2995.5	SuperNEMO, Lucifer
${100}_{Mo}\to {100}_{Ru}$	4.4	3035.0	MOON, AMoRe
${116}_{Cd}\to {116}_{Sn}$	4.6	2809	Cobra
${130}_{Te}\to {130}_{Xe}$	4.1	2530.3	CUORE
${136}_{Xe}\to {136}_{Ba}$	4.3	2461.9	EXO, KamLAND-Zen, NEXT, XMASS
${150}_{Nd}\to {150}_{Sm}$	19.2	3367.3	SNO+, DCBA/MTD

Table 3. Table 3: Experimental limits on LFV muon decays.

DECAY CHANNEL	EXPERIMENT	BRANCHING RATIO LIMIT
$\mu \to e\gamma$	MEG	muegamma
$\mu \to$eee	SINDRUM	SINDRUM
$\mu Au\to$ e Au	SINDRUM II	sindrum2

Equations

$$\mathrm{N}\left(\mathrm{A},\mathrm{Z}\right)\to \mathrm{N}\left(\mathrm{A},\mathrm{Z}+2\right)+{\mathrm{e}}^{-}+{\mathrm{e}}^{-},$$

$${\mathrm{\Gamma }}^{0\nu }=\frac{1}{{{\mathrm{T}}_{\frac{1}{2}}}^{0\nu }}={\mathrm{G}}^{0\nu }(\mathrm{Q},\mathrm{Z}){\left|{\mathrm{M}}^{0\nu }\right|}^2\frac{{\left|{\mathrm{m}}_{\beta \beta }\right|}^2}{{{\mathrm{m}}_{\mathrm{e}}}^2}.$$

$${\mathrm{m}}_{\beta \beta }=\sum _{\mathrm{i}}\mathrm{U}{{}_{\mathrm{ej}}}^2{\mathrm{m}}_{\mathrm{j}},\mathrm{j}=1,2,3,$$

$$\rm U_{PMNS}=\left[\begin{array}[]{ccc}c_{12}c_{13}&s_{12}c_{13}&s_{13}e^{-i\delta}\\ -c_{23}s_{12}-s_{23}s_{13}c_{12}e^{i\delta}&-c_{23}c_{12}-s_{23}s_{13}s_{12}e^{i\delta}&s_{23}c_{13}\\ s_{23}s_{12}-c_{23}s_{13}c_{12}e^{i\delta}&-s_{23}c_{12}-c_{23}s_{13}s_{12}e^{i\delta}&c_{23}c_{13}\end{array}\right]U_{Maj}.$$

$${\mathrm{m}}_{\beta \beta }={\mathrm{m}}_1{\mathrm{c}}_{12}^2{\mathrm{c}}_{13}^2+{\mathrm{m}}_2{\mathrm{s}}_{12}^2{\mathrm{c}}_{13}^2{\mathrm{e}}^{2i\alpha }+{\mathrm{m}}_3{\mathrm{s}}_{13}^2{\mathrm{e}}^{2i\beta }.$$

$$\mathrm{Q}={\mathrm{T}}_{3L}+{\mathrm{T}}_{3R}+\frac{\mathrm{B}-\mathrm{L}}{2}={\mathrm{T}}_{3L}+\mathrm{Y},$$

$$\rm q_{L}=\left[\begin{array}[]{c}u_{L}\\ d_{L}\end{array}\right]\equiv\left(3,2,1,\frac{1}{3}\right),q_{R}=\left[\begin{array}[]{c}u_{R}\\ d_{R}\end{array}\right]\equiv\left(3,1,2,\frac{1}{3}\right)$$

$$l_{L}=\left[\begin{array}[]{c}\nu_{L}\\ e_{L}\end{array}\right]\equiv\left(1,2,1,-1\right),l_{R}=\left[\begin{array}[]{c}\nu_{R}\\ e_{R}\end{array}\right]\equiv\left(1,1,2,-1\right)$$

$$\rm\phi=\left[\begin{array}[]{cc}\phi_{1}^{0}&\phi_{1}^{+}\\ \phi_{2}^{-}&\phi_{2}^{0}\end{array}\right]\equiv\left(\phi_{1},\widetilde{\phi_{2}}\right),$$

$$\rm\Delta_{L,R}=\left[\begin{array}[]{cc}{\delta_{\frac{L,R}{\sqrt{2}}}}^{+}&\delta_{L,R}^{++}\\ \delta_{L,R}^{0}&-{\delta_{\frac{L,R}{\sqrt{2}}}}^{+}.\end{array}\right]$$

$$<\phi>=\left[\begin{array}[]{cc}\frac{k_{1}}{\sqrt{2}}&0\\ 0&\frac{k_{2}}{\sqrt{2}}\end{array}\right]$$

$$<\Delta_{L,R}>=\left[\begin{array}[]{cc}0&0\\ \frac{v_{L,R}}{\sqrt{2}}&0\end{array}\right].$$

$$<{\Delta }_L>=v_L=\frac{\gamma k^2}{v_R}.$$

$$\mathcal{L}={\mathrm{h}}_{\mathrm{ij}}{\overline{\mathrm{\Psi }}}_{\mathrm{L},\mathrm{i}}\varphi {\mathrm{\Psi }}_{\mathrm{R},\mathrm{j}}+\widetilde{{\mathrm{h}}_{\mathrm{ij}}}{\overline{\mathrm{\Psi }}}_{\mathrm{L},\mathrm{i}}\tilde{\varphi }{\mathrm{\Psi }}_{\mathrm{R},\mathrm{j}}+{\mathrm{f}}_{\mathrm{L},\mathrm{ij}}{{\mathrm{\Psi }}_{\mathrm{L},\mathrm{i}}}^{\mathrm{T}}\mathrm{Ci}{\sigma }_2{\mathrm{\Delta }}_{\mathrm{L}}{\mathrm{\Psi }}_{\mathrm{L},\mathrm{j}}+{\mathrm{f}}_{\mathrm{R},\mathrm{ij}}{{\mathrm{\Psi }}_{\mathrm{R},\mathrm{i}}}^{\mathrm{T}}\mathrm{Ci}{\sigma }_2{\mathrm{\Delta }}_{\mathrm{R}}{\mathrm{\Psi }}_{\mathrm{R},\mathrm{j}}+h.c.$$

$$\rm M_{\nu}=\left[\begin{array}[]{cc}M_{LL}&M_{D}\\ {M_{D}}^{T}&M_{RR}\end{array}\right],$$

$${\mathrm{M}}_{\mathrm{D}}=\frac{1}{\sqrt{2}}({\mathrm{k}}_1\mathrm{h}+{\mathrm{k}}_2\tilde{\mathrm{h}}),{\mathrm{M}}_{\mathrm{LL}}=\sqrt{2}{\mathrm{v}}_{\mathrm{L}}{\mathrm{f}}_{\mathrm{L}},{\mathrm{M}}_{\mathrm{RR}}=\sqrt{2}{\mathrm{v}}_{\mathrm{R}}{\mathrm{f}}_{\mathrm{R}}.$$

$${\mathrm{M}}_{\nu }={{\mathrm{M}}_{\nu }}^{\mathrm{I}}+{{\mathrm{M}}_{\nu }}^{\mathrm{II}},$$

$${\mathrm{M}}_{\nu }={\mathrm{M}}_{\mathrm{LL}}+{\mathrm{M}}_{\mathrm{D}}{{\mathrm{M}}_{\mathrm{RR}}}^{-1}{{\mathrm{M}}_{\mathrm{D}}}^{\mathrm{T}}=\sqrt{2}{\mathrm{v}}_{\mathrm{L}}{\mathrm{f}}_{\mathrm{L}}+\frac{{\mathrm{k}}^2}{\sqrt{2}{\mathrm{v}}_{\mathrm{R}}}{\mathrm{h}}_{\mathrm{D}}{{\mathrm{f}}_{\mathrm{R}}}^{-1}{{\mathrm{h}}_{\mathrm{D}}}^{\mathrm{T}}.$$

$${\mathrm{h}}_{\mathrm{D}}=\frac{({\mathrm{k}}_1\mathrm{h}+{\mathrm{k}}_2\tilde{\mathrm{h}})}{\sqrt{2}\mathrm{k}},\mathrm{k}=\sqrt{{\left|{\mathrm{k}}_1\right|}^2+{\left|{\mathrm{k}}_2\right|}^2}$$

$${\mathrm{v}}_{\mathrm{L}}=\frac{\gamma {{\mathrm{M}}_{\mathrm{W}}}^2}{{\mathrm{v}}_{\mathrm{R}}}.$$

$${\mathrm{M}}_{\nu }=\gamma {\left(\frac{{\mathrm{M}}_{\mathrm{W}}}{{\mathrm{v}}_{\mathrm{R}}}\right)}^2{\mathrm{M}}_{\mathrm{RR}}+{\mathrm{M}}_{\mathrm{D}}{{\mathrm{M}}_{\mathrm{RR}}}^{-1}{{\mathrm{M}}_{\mathrm{D}}}^{\mathrm{T}}.$$

$$\gamma =\frac{{\beta }_1{\mathrm{k}}_1{\mathrm{k}}_2+{\beta }_2{{\mathrm{k}}_1}^2+{\beta }_3{{\mathrm{k}}_2}^2}{(2{\rho }_1-{\rho }_3){\mathrm{k}}^2}.$$

$${{\mathrm{A}}_{\nu }}^{\mathrm{LL}}\cong \frac{1}{{{\mathrm{M}}_{{\mathrm{W}}_{\mathrm{L}}}}^4}\sum \frac{{{\mathrm{U}}_{{\mathrm{L}}_{{\mathrm{e}}_{\mathrm{i}}}}}^2{\mathrm{m}}_{\mathrm{i}}}{{\mathrm{p}}^2}.$$

$${{\mathrm{A}}_{\mathrm{N}}}^{\mathrm{RR}}\propto \frac{1}{{{\mathrm{M}}_{{\mathrm{W}}_{\mathrm{L}}}}^4}\frac{{{\mathrm{U}}_{{\mathrm{R}}_{{\mathrm{e}}_{\mathrm{i}}}}^{\ast }}^2}{{\mathrm{M}}_{\mathrm{i}}}.$$

$${{\mathrm{A}}_{{\mathrm{\Delta }}_{\mathrm{R}}}}^{\mathrm{RR}}\propto \frac{1}{{{\mathrm{M}}_{{\mathrm{W}}_{\mathrm{R}}}}^4}\frac{1}{{{\mathrm{M}}_{{\mathrm{\Delta }}_{\mathrm{R}}}}^2}{\mathrm{f}}_{\mathrm{R}}{\mathrm{v}}_{\mathrm{R}}.$$

$${\mathrm{\Gamma }}_{\mu }\equiv \mathrm{\Gamma }\left({\mu }^{-}\to {\mathrm{e}}^{-}{\nu }_{\mu }\overline{{\nu }_{\mathrm{e}}}\right),{\mathrm{\Gamma }}_{\mathrm{capt}}^{\mathrm{Z}}\equiv \mathrm{\Gamma }\left({\mu }^{-}+\mathrm{A}\left(\mathrm{Z},\mathrm{N}\right)\to {\nu }_{\mu }+\mathrm{A}\left(\mathrm{Z}-1,\mathrm{N}+1\right)\right).$$

$$\mathrm{B}{\mathrm{R}}_{\mu \to \mathrm{e}\gamma }\equiv \frac{\mathrm{\Gamma }\left({\mu }^{+}\to {\mathrm{e}}^{+}\gamma \right)}{{\mathrm{\Gamma }}_{\mu }},$$

$$\mathrm{B}{\mathrm{R}}_{\mu \to \mathrm{e}}^{\mathrm{Z}}\equiv \frac{\mathrm{\Gamma }\left({\mu }^{-}+\mathrm{A}\left(\mathrm{N},\mathrm{Z}\right)\to {\mathrm{e}}^{-}+\mathrm{A}\left(\mathrm{N},\mathrm{Z}\right)\right)}{{\mathrm{\Gamma }}_{\mathrm{capt}}^{\mathrm{Z}}},$$

$$\mathrm{B}{\mathrm{R}}_{\mu \to 3e}\equiv \frac{\mathrm{\Gamma }\left({\mu }^{+}\to {\mathrm{e}}^{+}{\mathrm{e}}^{-}{\mathrm{e}}^{+}\right)}{{\mathrm{\Gamma }}_{\nu }}.$$

$$\mathrm{BR}\left(\mu \to 3e\right)=\frac{1}{2}{\left|{\mathrm{h}}_{\mu \mathrm{e}}{\mathrm{h}}_{\mathrm{ee}}^{\ast }\right|}^2\left(\frac{{{\mathrm{m}}_{{\mathrm{W}}_{\mathrm{L}}}}^4}{{{\mathrm{M}}_{{\mathrm{\Delta }}_{\mathrm{L}}}^{++}}^4}+\frac{{{\mathrm{m}}_{{\mathrm{W}}_{\mathrm{R}}}}^4}{{{\mathrm{M}}_{{\mathrm{\Delta }}_{\mathrm{R}}}^{++}}^4}\right).$$

Full text

I INTRODUCTION

The milestone discovery of neutrino oscillation and the corresponding realization that neutrinos are massive particles has been one of the compelling revelation which suggests physics beyond the Standard Model (SM). The recent neutrino experiments MINOS minos , T2K t2k , Double Chooz doublechooz , Daya Bay dayabay and RENO reno have not only confirmed the earlier observations but also measured the neutrino parameters more accurately. The 3$\sigma$ global fit values of the neutrino oscillation parameters according to recent analysis are shown in the table 1.

Notwithstanding, the absolute neutrino mass scale is still unperceived. However, the Planck experiment has given an upper bound on the sum of the light neutrino mass to be $\rm\sum_{i}\left|m_{i}\right|<$0.23 eV planck in 2012 and recently the bound has been constrained to $\rm\sum_{i}\left|m_{i}\right|<$0.17 eV planck1 . The simplest hypothesis (way) to account for a neutrino mass is to introduce atleast two right handed (RH) neutrino in the Standard Model (SM). This will allow a Dirac coupling with the Higgs, like other fermions in the SM. However, corresponding Yukawa coupling has to be fine tuned around $\rm 10^{-12}$ which is quite unnatural. This kind of fine tuning can be avoided to explain the neutrino masses in the seesaw mechanism, a mechanism beyond SM (BSM) physics which is categorised into type I type1 , type II type2 , type III type3 , inverse inverse seesaw mechanism. The BSM physics also unveils various phenomenon like Baryon Asymmetry of the Universe (BAU), Lepton Number Violation (LNV), Lepton Flavour Violation (LFV), existence of dark matter etc. One of the theoretical framework to make the first three processes observable is the left-right symmetric model (LRSM) LRSM which is considered to be an appealing candidate for physics BSM. Here, the gauge group is very simple extension of the SM gauge group. It provides a natural framework to understand the spontaneous breaking of parity and origin of small neutrino mass via seesaw mechanism.

Furthermore, the physics community worldwide is embarking on the next challenging problem in finding out the nature of the neutrinos, whether they are four component Dirac particles possesing a conserved lepton number or two component Majorana particles, along with the absolute scale of neutrino mass. This problem is directly related to the issue of LN conservation, which is one of the most obscure sides of the SM not supported by an underlying principle. One of such process of fundamental importance in particle physics which pops up almost in any extension of the SM is neutrinoless double beta decay (NDBD) NDBD . It is defined as a second order, slow radioactive process that transforms a nuclide of atomic number Z into its isobar with atomic number Z+2,

[TABLE]

thereby violating the total LN conservation. Its existence is directly linked to that of the Majorana neutrinos maj (i.e., identical to its own anti particle).

The general expression for the total decay width of $0\nu\beta\beta$ taking into account the coulomb interaction of the electrons and the final nucleus is given by,

[TABLE]

The numerical values of $\rm G^{0\nu}(Q,Z)$, Q and the natural abundance of several nuclei of experimental interest are given in the table 2 which are adopted from reference table .

The main aim of the experiment on the search for 0$\rm\nu\beta\beta$ decay is the measurement of the effective Majorana neutrino mass, which is a combination of the neutrino mass eigenstates and neutrino mixing matrix terms, given by,

[TABLE]

where, $\rm U_{ej}$ are the elements of the first row of the neutrino mixing matrix, $\rm U_{PMNS}$ (dependent on the known parameters $\rm\theta_{13},\theta_{12}$ and the unknown Majorana phases $\rm\alpha$ and $\rm\beta$ alpha ). $\rm U_{PMNS}$ is the diagonalizing matrix of the light neutrino mass matrix, $m_{\nu}$ given by equation 4.

[TABLE]

The abbreviations used are $c_{ij}$= $\cos\theta_{ij}$, $s_{ij}$=$\sin\theta_{ij}$, $\delta$ is the Dirac CP phase while the diagonal phase matrix, $\rm U_{Maj}$ is $\rm diag(1,e^{i\alpha},e^{i(\beta+\delta)})$ P contains the Majorana phases $\rm\alpha$ and $\rm\beta$. The Majorana phases $\rm\alpha$ and $\rm\beta$ have an effect in the process, which are allowed only if massive neutrinos are Majorana particles and are characterized by a violation of total LN, such as NDBD. In the standard parameterization of the mixing matrix, $\rm m_{\beta\beta}$ is given by,

[TABLE]

A huge amount of experimental and theoretical activity is persued in order to detect and predict the decay process. Although no convincing experimental evidence of the decay exists till date, but new generation of experiments that are already running or about to run assures to expedite the current limits exploring the degenerate-hierarchy region of neutrino masses. In addition, from the life time of this decay combined with sufficient knowledge of the nuclear matrix elements (NME), one can set a constraint involving the neutrino masses. Moreover, if one incorporates the recent results of neutrino oscillation experiments, one can set a stringent limit on the neutrino mass scale. The latest experiments ndbd that have improved the lower bound of the half life of the decay process include KamLAND-Zen kamland and GERDA gerda which uses Xenon-136 and Germanium-76 nuclei respectively. Incorporating the results from first and second phase of the experiment, KamLAND-Zen imposes the best lower limit on the decay half life using Xe-136 as $\rm T_{1/2}^{0\nu}>1.07\times 10^{26}$ yr at $90\%$ CL and the corresponding upper limit of effective Majorana mass in the range (0.061-0.165)eV.

Again one of the most important BSM framework to understand the origin of neutrino mass and large leptonic mixing is to identify the possible underlying symmetries. Symmetries can relate two or more free parameters of the model or make them vanish, making the model more predictive. The widely studied $\rm\mu-\tau$ symmetric mutau neutrino mass matrix giving zero $\rm\theta_{13}$ is one such scenerio where discrete flavor symmetries can relate two or more terms in the neutrino mass matrix. The neutrino oscillation data before the discovery of non zero $\rm\theta_{13}$ were in perfect agreement with $\rm\mu-\tau$ symmetric neutrino mass matrix. The four different realizations of neutrino mixing pattern generally found in literature which can generate from $\rm\mu$-$\tau$ symmetric mass matrices are tribimaximal mixing (TBM), bimaximal mixing (BM), hexagonal mixing (HM), golden ratio mixing (GRM) matrices. But, after discovery of non zero $\rm\theta_{13}$, one needs to go beyond these $\rm\mu$-$\tau$ symmetric framework. Since the experimental value of $\theta_{13}$ is still much smaller than the other two mixing angles, $\rm\mu$-$\tau$ symmetry can still be a valid approximation and the non zero $\rm\theta_{13}$ nonzero1 can be accounted for by incorporating the presence of small perturbation to $\rm\mu$-$\tau$ symmetry.

The discovery of neutrino oscillation has provided clear evidence of the fact that neutrinos are massive as well as the violation of the lepton flavour lfv1 during the propagation of the neutrinos. Lepton flavour is consequently a broken symmetry and the SM has to be adapted to incorporate massive neutrinos and thus we can also hope that lepton flavour violation (LFV) will be visible in the charged lepton sector lfv5 . The exact mechanism of LFV being unknown, its study is of large interest as it is linked to neutrino mass generation, CP violation and new physics BSM. The LFV effects from new particles at TeV scale are naturally generated in many models and therefore considered to be a prominent signature for new physics. In LRSM, where electroweak symmetry is broken dynamically, an experimentally accessible amount of LFV is predicted in a large region of parameter space. In a wide range of models for physics BSM, highest sensitivity in terms of BR is expected for $\rm\mu\rightarrow 3e$ and $\rm\mu\rightarrow e\gamma$ decay processes.

To study these phenonomenon theoretically or phenomenologically, many works have been performed in LRSM based framework ls1 . In most of these works, authors mostly considered the TBM like neutrino mass as leading order contribution and arising from type I seesaw and using the type II seesaw as a perturbation to generate non zero $\rm\theta_{13}$ perturbations . More recently, the authors of ndbddb2 ndbddb1 studied the new physics contribution to NDBD with prominent type I and type II as well as equally dominating type I and type II seesaw. Again, many works have been done in charged lepton flavour violation sector in literature considering type I and type II dominant cases as well as equally dominant type I and type II in the TeV scale LRSM framework which is within the presently accessible reach of the colliders and implements the two seesaw mechanisms naturally lfv5 .

In this context, we present a phenomenological study of different $\rm\mu-\tau$ symmetric mutau neutrino mass models to check their consistency with the stringent constraints from cosmology, with various processes like LNV, LFV etc. We have taken the leading order mass matrices obeying $\rm\mu$-$\tau$ symmetry originating from type I (II) seesaw then incorporating type II (I) seesaw as perturbations to generate non zero $\theta_{13}$. Then we studied the LFV in the LRSM framework and further correlated the LFV of the processes $\left(\mu\rightarrow e\gamma\right)$ and $\left(\rm\mu\rightarrow 3e\right)$ with lightest neutrino mass and atmospheric mixing angle, $\rm\theta_{23}$ in different neutrino mass models favouring $\mu-\tau$ symmetry. In NDBD, we discuss the different contributions ndbddb2 from right handed (RH) neutrinos and RH gauge bosons, triplet Higgs higgs as well as light heavy neutrino mixing that can contribute to the effective mass governing the process and identify the significant ones. In this work, we have considered only the dominant new physics contribution as coming from the diagrams containing purely RH current mediated by the heavy gauge boson, $\rm W_{R}$ by the exchange of heavy right handed neutrino, $\rm N_{R}$ and another from the charged Higgs scalar $\rm\Delta_{R}$ mediated by the heavy gauge boson $\rm W_{R}$ ndbddb1 . We have ignored the contributions coming from the left-right gauge boson mixing and heavy light neutrino mixing.

This paper is structured as follows. In section II we briefly discuss the left-right symmetric model framework and the origin of neutrino mass and summarize the NDBD process in this framework in section III. We also discuss the different feynmann diagrams contributing to the amplitude of the decay process (the new physics contribution) in this section. In section IV, we briefly discuss lepton flavor violating processes, mainly $\rm\left(\mu\rightarrow 3e\right)$ and $\rm\left(\mu\rightarrow e\gamma\right)$. In section V, we present our numerical analysis and results and then in section VI, we conclude by giving a brief overview of our work.

II LEFT RIGHT SYMMETRIC MODEL(LRSM) AND NEUTRINO MASS

The explaination of the smallness of neutrino mass and the profile of its mixing as required by the recent experiments has been taken as a great puzzle in particle physics. The fact that neutrino has mass implies the requirement of new physics beyond the $\rm SU(3)_{c}\times SU(2)_{L}\times U(1)_{Y}$ SM electroweak . One possibility to introduce neutrino mass is the so called seesaw mechanism wherein we introduce right handed heavy singlet neutrino, $\rm\nu_{R}$ (type I seesaw), scalar Higgs triplet (type II seesaw) and hypercharge-less fermion triplets (type III seesaw). Left-right symmetric model (LRSM) LRSM can be considered to be very appealing model for Physics beyond the Standard model. The seesaw mechanisms can be realized in the context of left-right symmetric model or GUTs where seesaw scale might be related to other physical scales.

In LRSM, the gauge group is a very simple extension of the standard model gauge group, $\rm SU(3)_{c}\times SU(2)_{L}\times U(1)_{Y}$. Most of the problems like parity violation of weak interaction, masssless neutrinos, CP problems, hierarchy problems etc can be explained in the framework of LRSM, based on the gauge group, $\rm SU(3)_{c}\times SU(2)_{L}\times SU(2)_{R}\times U(1)_{B-L}$ genericlrsm LRSM . In this model, the electric charge is related to the generators of the group as

[TABLE]

where $\rm T_{3L}$ and $\rm T_{3R}$ are the 3rd components of isospin under $\rm SU(2)_{L}$ and $\rm SU(2)_{R}$. In LRSM, the left and right handed components of the fields are treated on the same footing. If the Higgs sector of the model is choosen so that RH symmetry is spontaneously broken by triplets, the model gives rise to tiny neutrino masses naturally via seesaw mechanism. Herein, there are 2 sources of lepton number violation, the Majorana masses of neutrinos and Yukawa interaction of triplet Higgs. The Quarks and leptons transform under the L-R symmetric gauge group as,

[TABLE]

[TABLE]

We consider the general class of left-right symmetric model which are invariant under $\rm SU(3)_{c}\times SU(2)_{L}\times SU(2)_{R}\times U(1)_{B-L}$ symmetry with the Higgs content, $\rm\phi(1,2,2,0)$, $\Delta_{L}(1,2,1,-1)$, $\Delta_{R}(1,1,2,-1)$. A convenient representation of fields is given by 2$\times$ 2 matrices for the Higgs bidoublets and the $\rm SU(2)_{L,R}$ triplets as,

[TABLE]

[TABLE]

The neutral Higgs fields $\rm\delta_{L,R}^{0}$, $\rm\phi_{1}^{0},\phi_{2}^{0}$ can potentially acquire VEVS $\rm v_{R},v_{L},k_{1},k_{2}$ respectively.

[TABLE]

[TABLE]

The VEV $\rm v_{R}$ breaks the $\rm SU(2)_{R}$ symmetry and sets the mass scale for the extra gauge bosons $\rm(W_{R}$ and Z$\rm{}^{\prime})$ and for right handed neutrino field $\rm(\nu_{R})$. The VEVs $\rm k_{1}$ and $\rm k_{2}$ serves the twin purpose of breaking the remaining the $\rm SU(2)_{L}\times U(1)_{B-L}$ symmetry down to $\rm U(1)_{em}$, thereby setting the mass scales for the observed $\rm W_{L}$ and Z bosons and providing Dirac masses for the quarks and leptons. Clearly, $\rm v_{R}$ must be significantly larger than $\rm k_{1}$ and $\rm k_{2}$ in order for $\rm W_{R}$ and Z ′ to have greater masses than the $W_{L}$ and Z bosons. $v_{L}$ is the VEV of $\Delta_{L}$, it plays a significant role in the seesaw relation which is the characteristics of the LR model and can be written as,

[TABLE]

The acceptable breaking pattern is, $\rm SU(2)_{L}\times SU(2)_{R}\times U(1)_{B-L}\xrightarrow{<\Delta_{R}>}SU(2)_{L}\times U(1)_{Y}\xrightarrow{<\phi>}U(1)_{em}$.

The Yukawa lagrangian in the lepton sector is given by,

[TABLE]

Where the family indices i, j are summed over, the indices i, j=1, 2, 3 represents the three generations of fermions. $\rm C=i\gamma_{2}\gamma_{0}$ is the charge conjugation operator, $\rm\widetilde{\phi}=\tau_{2}\phi^{*}\tau_{2}$ and $\gamma_{\mu}$ are the Dirac matrices. Considering discrete parity symmetry, the Majorana Yukawa couplings $f_{L}=f_{R}$ (for left-right symmetry) gives rises to Majorana neutrino mass after electroweak symmetry breaking when the triplet Higgs $\rm\Delta_{L}$ and $\rm\Delta_{R}$ acquires non zero vacuum expectation value. Then equation (14) leads to $6\times 6$ neutrino mass matrix as shown in reference 2 of ls1

[TABLE]

where

[TABLE]

Where $M_{D}$, $M_{LL}$ and $M_{RR}$ are the Dirac neutrino mass matrix, left handed and right handed mass matrix respectively. Assuming $M_{L}\ll M_{D}\ll M_{R}$, the light neutrino mass, generated within a type I+II seesaw can be written as,

[TABLE]

[TABLE]

Where the first and second terms in equation (18) corresponds to type II seesaw and type I seesaw mediated by RH neutrino respectively. Here,

[TABLE]

In the context of LRSM both type I and type II seesaw terms can be written in terms of $M_{RR}$ which arises naturally at a high energy scale as a result of spontaneous parity breaking. In LRSM the Majorana Yukawa couplings $f_{L}$ and $f_{R}$ are same (i.e, $f_{L}=f_{R}$) and the VEV for left handed triplet $v_{L}$ can be written as,

[TABLE]

Thus equation (18) can be written as ,

[TABLE]

In literature, (reference breaking ndbddb2 ) author define the dimensionless parameter $\gamma$ as,

[TABLE]

Here the terms $\beta$, $\rho$ are the dimensionless parameters that appears in the expression of the Higgs potential.

III $0\nu\beta\beta$ decay in LRSM

Many theoretical groups has studied NDBD in connection with LRSM dbd1 . In the context of LRSM, there are several contributions to NDBD in addition to the standard contribution via light Majorana neutrino exchange owing to the presence several heavy additional scaler, vector and fermionic fields . Many of the earlier works have explained it in details with the corresponding feynmann diagrams (see ref. ndbddb2 ). The various contributions to $\rm 0\nu\beta\beta$ decay transition rate in LRSM are briefly summarized below.

• •

Standard Model contribution to NDBD where the intermediate particles are the $\rm{W_{L}}^{-}$ bosons and light neutrinos. The amplitude of this process depends upon the leptonic mixing matrix elements elements and light neutrino masses.

• •

Heavy right handed neutrino contribution to NDBD in which the mediator particles are the $\rm{W_{L}}^{-}$ bosons. The amplitude of this process depends upon the mixing between light and heavy neutrinos as well the mass of the heavy neutrino, $\rm N_{i}$.

• •

Light neutrino contribution to NDBD in which the intermediate particles are $\rm{W_{R}}^{-}$ bosons. The amplitude of this process depends upon the mixing between light and heavy neutrinos as well as the mass of the right handed gauge boson, $\rm{W_{R}}^{-}$ boson.

• •

Heavy right handed neutrino contribution to NDBD in which the mediator paticles are the $\rm{W_{R}}^{-}$ bosons. The amplitude of this process depends upon the elements of the right handed leptonic mixing matrix and the mass of the right handed gauge boson, $\rm{W_{R}}^{-}$ boson as well as the mass of the heavy right handed Majorana neutrino, $\rm N_{i}$.

• •

Light neutrino contribution from the Feynman diagram mediated by both $\rm{W_{L}}^{-}$ and $\rm{W_{R}}^{-}$. The amplitude of this process depends upon the mixing between light and heavy neutrinos, leptonic mixing matrix elements, light neutrino masses and the mass of the gauge bosons, $\rm{W_{L}}^{-}$ and $\rm{W_{R}}^{-}$.

• •

Heavy neutrino contribution from the Feynman diagram mediated by both $\rm{W_{L}}^{-}$ and $\rm{W_{R}}^{-}$. The amplitude of the process depends upon the right handed leptonic mixing matrix elements, mixing between the light and heavy neutrinos as well as the mass of the gauge bosons, $\rm{W_{L}}^{-}$ and $\rm{W_{R}}^{-}$ and the mass of the heavy right handed neutrino, $\rm M_{i}$.

• •

Triplet Higgs $\rm\Delta_{L}$ contribution to NDBD in which the mediator particles are $\rm{W_{L}}^{-}$ bosons. The amplitudes for the process depends upon the masses of the $\rm{W_{L}}^{-}$ bosons, left handed triplet Higgs, $\rm\Delta_{L}$ as well as their coupling to leptons, $\rm f_{L}$.

• •

Right handed triplet Higgs $\rm\Delta_{R}$ contribution to NDBD in which the mediator particles are ${W_{R}}^{-}$ bosons. The amplitude for the process depends upon the masses of the $\rm{W_{R}}^{-}$ bosons, right handed triplet Higgs, $\rm\Delta_{R}$ as well as their coupling to leptons, $f_{R}$.

However in our present work, we have considered only three of the above mentioned contributions to NDBD. One from the standard light neutrino contribution through exchange of $\rm{W_{L}}^{-}$ as shown in figure 1(a) and the other two are the new physics contributions to NDBD which corresponds to figures 1(b) and 1(c), that is the ones mediated by ${W_{R}}^{-}$ and $\Delta_{R}$ respectively. The amplitudes of the contributions are given in several earlier works like ndbddb2 . For simple approximations, an assumption of similar mass scales for the heavy particles has been made in the LRSM, where, $\rm M_{R}\approx M_{W_{R}}\approx M_{\Delta_{L}^{++}}\approx M_{\Delta_{R}^{++}}\approx TeV$, at a scale accessible at the LHC. Under these assumptions, the amplitude for the light-heavy mixing contribution which is proportional to $\rm\frac{{m_{D}}^{2}}{M_{R}}$ remains very small (since $\rm m_{\nu}\approx\frac{{m_{D}}^{2}}{M_{R}}\approx(0.01-0.1)eV$, $m_{D}\approx(10^{5}-10^{6})$ eV which implies $\rm\frac{m_{D}}{M_{R}}\approx(10^{-7}-10^{-6})$ eV) Thus, we ignore the contributions involving the light and heavy neutrino mixings. For a simplified approach, we have also ignored the mixing between $\rm{W_{L}}^{-}$ and $\rm{W_{R}}^{-}$ bosons owing to the above mentioned assumptions, which would cause a further supression in the amplitude of the process (for reference see ndbddb1 ). Again, the contribution from $\rm{\Delta_{L}}^{-}$, $\rm{W_{L}}^{-}$ is suppressed by the type II seesaw contribution to light neutrino mass and hence neglected here.

Considering these contributions we have studied the NDBD. Different neutrino mass satisfying the mixing criteria namely, TBM, BM, HM and GRM are considered as leading contribution in either type I or type II seesaw. The perturbation is added for generation of non zero $\theta_{13}$ nonzero in either of the seesaw terms.

The amplitude of the corresponding processes which we have considered in our work are given by,

• •

Standard light neutrino contribution,

[TABLE]

where, $\left|p\right|\sim$ 100 MeV momentum is the typical momentum transfer at the leptonic vertex, $\rm m_{ee}=\sum{U_{L_{e_{i}}}}^{2}m_{i}$ is the effective neutrino mass. $\rm U_{L_{e_{i}}}$ represents the elements of the first row of the neutrino mixing matrix, $\rm U_{PMNS}$.

• •

Heavy RH neutrino contribution,

[TABLE]

• •

Scalar triplet contribution,

[TABLE]

.

Here, $\rm\rm U^{*}_{R_{e_{i}}}$ denotes the first row of the unitary matrix diagonalizing the right handed neutrino mass matrix, $\rm M_{RR}$ with mass eigen values, $\rm M_{i}$.

IV Lepton Flavour Violation (LFV)

There have been various attempts to observe and predict theoretically the manifestation of LFV involving various modes of muon decay since long. The most promising LFV low energy channels are probably $\mu\rightarrow e\gamma$, $\mu\rightarrow 3e$, $\mu\rightarrow e$ conversion in nuclei which occur in rates accessible in recent experiments. Defining the decay rates (from reference lfv1 )as,

[TABLE]

The relevant branching ratios (BR) for the processes are,

[TABLE]

[TABLE]

[TABLE]

The selected limits for lepton flavour violating muon decays and muon to electron conversion experiments are shown in table 3

In the SM seesaw, the LFV decay rates induced by neutrino mixing are suppressed by tiny neutrino masses,$\left(\frac{{\Delta m_{A}}^{2}}{{M_{W}}^{2}}\right)\sim 10^{-50}$ and hence are well below the current experimental limits and even the distant future sensitivities. New physics beyond the standard model is required to make the process observable, there are several theoretical frameworks BSM that could provide the necessary operators. One of those theories is the LRSM in which several new contributions appear due to the additional RH current interactions, which could lead to sizeable LFV rates for TeV scale $v_{R}$ that occur at rates observable in current experiments. LFV in the LRSM has been studied in many previous works. There are various LFV processes providing constraints on the masses of the right handed neutrinos and doubly charged scalars. It turns out that the process $\mu\rightarrow 3e$ induced by doubly charged bosons $\Delta_{L}^{++}$ and $\Delta_{R}^{++}$ and $\mu\rightarrow e\gamma$ provides the most relevant constraint. In our present work, we consider these processes in minimal left-right symmetric model (MLRSM). The limit of branching ratio of the process $\mu\rightarrow 3e$ as shown in table 3 is $<1.0\times 10^{-12}$ at $90\%$ CL was obtained at the Paul Scherrer institute (PSI) over 20 years ago by the SINDRUM experiment SINDRUM . Presently the Mu3e collaboration has submitted a letter of intent to PSI to perform a new improved search for the decay $\mu\rightarrow 3e$ with a sensitivity of $10^{-16}$ at $95\%$ CL sindrum2 which corresponds to an improvement by four orders of magnitude compared to the former SINDRUM experiment. Whereas the new upper limit for BR of the process $\mu\rightarrow e\gamma$ is established to be $<4.2\times 10^{-13}$ at $90\%$ CL by the MEG collaboration. Taking into account the contributions from heavy righthanded neutrinos and Higgs scalars, the expected branching ratios and conversion rates of the above processes have been calculated in the LRSM in the work (first reference in lfv ).

The BR for the process $(\mu\rightarrow 3e)$ is given by,

[TABLE]

Where $h_{ij}$ describes the lepton Higgs coupling in LRSM and is given by,

[TABLE]

For $\mu\rightarrow e\gamma$, the relevant BR is given by, lfv1

[TABLE]

where, $g_{lfv}$ is defined as,

[TABLE]

The sum is over the heavy neutrinos only. $\rm M_{\Delta_{L,R}}^{++}$ are the masses of the doubly charged bosons, $\rm{\Delta_{L,R}}^{++}$, V is the mixing matrix of the right handed neutrinos with the electrons and muons. $\rm M_{n}(n=1,2,3)$ are the right handed neutrino masses.

V NUMERICAL ANALYSIS AND RESULTS

In our present work we have studied LNV (NDBD) for standard as well as non standard contributions for the effective mass as well as the half life governing the decay process in the framework of LRSM. We have also correlated the LFV of the process, $\rm\mu\rightarrow 3e$ and $\rm\mu\rightarrow e\gamma$ with the lightest neutrino mass and atmospheric mixing angle, $\rm\theta_{23}$ respectively for both normal and inverted mass hierarchies. In this section we present a detailed analysis of our work and we have divided it into different subsections, firstly the standard light neutrino contribution to NDBD and then the new physics contribution to NDBD considering perturbation in type II and then type I seesaw. Lastly we have shown the analysis of correlating LFV with $\rm m_{lightest}$ and $\rm\theta_{23}$.

V.1 Standard light neutrino contribution

For NDBD mediated by the light Majorana neutrinos, the half life of the decay process is given by equation (1) and the effective mass governing the process is as given in equation (5). In our present work, we first evaluated the effective light neutrino mass within the standard mechanism using the formula (3) where, $\rm U_{ej}$ are the elements of the first row of the neutrino mixing matrix, $\rm U_{PMNS}$ (dependent on the known parameters $\rm\theta_{13},\theta_{12}$ and the unknown Majorana phases $\rm\alpha$ and $\rm\beta$). $\rm U_{PMNS}$ is the diagonalizing matrix of the light neutrino mass matrix, $\rm m_{\nu}$, such that

[TABLE]

where $\rm{M_{\nu}}^{(diag)}=diag(m_{1},m_{2},m_{3})$. In the case of 3 neutrino mixing, 2 $\nu$ mass spectra are possible,

• •

Normal Hierachy $\rm(NH)$ which corresponds to $m_{1}<m_{2}\ll m_{3}$ ; $\Delta{m_{12}}^{2}\ll\Delta{m_{23}}^{2}$.

• •

Inverted Hierarchy $\rm(IH)$ which corresponds to $m_{3}\ll m_{1}\sim m_{2}$ ; $\Delta{m_{12}}^{2}\ll\left|\Delta{m_{13}}^{2}\right|$.

In both the spectra, $\rm\Delta{m_{12}}^{2}=\Delta{m_{solar}}^{2}$. For NH, $\rm\Delta{m_{23}}^{2}=\Delta{m_{atm}}^{2}$ and for IH, $\left|\Delta{m_{13}}^{2}\right|=\Delta{m_{atm}}^{2}$. In the case of NH, the neutrino masses $m_{2}$ and $m_{3}$ are connected with the lightest mass $m_{1}$ by the relation,

[TABLE]

In IH, $\rm m_{3}$ is the lightest mass and we have,

[TABLE]

For both the normal and inverted hierarchies, equation (5) can be written in terms of lightest neutrino mass as,

for NH,

[TABLE]

for IH,

[TABLE]

The $3\sigma$ ranges of the mass squared differences and mixing angles from global analysis of oscillation data are outlined as in the table 1. Using the best fit values of the mass squared differences and the $\rm 3\sigma$ ranges of the three mixing angles from a global analysis of oscillation data (as shown in table 1), we have shown the variation of the effective Majorana mass as a function of the lightest neutrino mass $\rm m_{1}$ (for NH) and $\rm m_{3}$ (for IH). During our calculation, we have varied the Majorana phase $\rm\alpha$ and $\rm\beta$ from 0 to $\rm 2\pi$. The effective mass assumes different values depending on whether the neutrino mass states follows normal hierarchy (NH) or inverted hierarchy (IH). We have used equations (37)and (38) in evaluating the effective mass in terms of the lightest neurino mass. The variation is shown in figure 2. It is seen from the figure that the light neutrino contribution to neutrinoless double beta decay $(0\nu\beta\beta)$ can saturate the bound imposed by KamLAND-ZEN $\rm(\leq 0.061-0.165eV)$ [reference (kamland )] only for the higher values of lightest neutrino masses which is disallowed by the Planck data (lightest mass for NH $\sim$ 0.07 and lightest mass for IH $\sim$ 0.065).

Again, we have evaluated the effective majorana mass for different leptonic mixing patterns possessing $\rm\mu-\tau$ symmetry, namely, tribimaximal, golden ratio and hexagonal mixing mutau using equation (37)and (38). In all the different $\rm\mu-\tau$ symmetric mixing patterns which we have considered, i.e., TBM, HM, GRM, the reactor mixing angle $\theta_{13}$ is 0 and $\theta_{23}$ is $45^{0}$. Whereas $\rm\theta_{12}$=$35.5^{0}$,(for TBM), $\theta_{12}=30^{0}$(for HM), $\theta_{12}$=$31.71^{0}$(for GRM). Since, $\rm\theta_{12}$=$45^{0}$, i.e, BM has been ruled out by experiments, we have ignored this case for the standard contributions. Again, it is to be noted that there are two values of $\rm\theta_{12}$ for GRM, which are, $31.7^{0}$ and $35.96^{0}$ GRM . In our present study, we have considered the first value which is allowed as mentioned in reference GRM GRM1

The variations of $\rm m_{\nu}^{eff}$ for the different mixing patterns for NH and IH in terms of lightest neutrino mass are shown as in figure 3.

V.2 New physics contribution to NDBD considering perturbation in type II seesaw.

For the new Physics contribution, we have considered the contributions of $\rm 0\nu\beta\beta$ from the right handed current and from the triplet Higgs $\rm(\Delta_{R})$. The contributions from the left handed Higgs triplet, $\rm\Delta_{L}$ is suppressed by the light neutrino mass. Also we consider the mixing between LH and RH sector to be so small that their contributions to $0\nu\beta\beta$ can be neglected. The total effective mass is thus given by the formula, (as in ndbddb1 )

[TABLE]

Here, $\rm<p^{2}>=m_{e}m_{p}\frac{M_{N}}{M_{\nu}}$ is the typical momentum exchange of the process, where $\rm m_{p}$ and $\rm m_{e}$ are the mass of the proton and electron respectively and $\rm M_{N}$ is the NME corresponding to the RH neutrino exchange. We know that TeV scale LRSM plays an important role in 0$\nu\beta\beta$ decay. We have considered the values $\rm M_{W_{R}}$ = 3.5 TeV, $\rm M_{W_{L}}$ = 80 GeV, $\rm M_{\Delta_{R}}\approx$3TeV, the heavy RH neutrino $\approx$ TeV which are within the recent collider limits collider . The allowed value of p (the virtuality of the exchanged neutrino) is in the range $\sim$ (100-200) MeV. In our analysis, we have taken p$\simeq$180 MeV ndbddb2 .

Thus,

[TABLE]

However, equation (39) is valid only in the limit $\rm{M_{i}}^{2}\gg\left|<p^{2}>\right|$ and $\rm{M_{\Delta}}^{2}\gg\left|<p^{2}>\right|$.

The formula for light $\nu$ masses in the presence of both type I and type II seesaw can be written as,

[TABLE]

[TABLE]

where, $\rm U_{PMNS}$ and $U_{(\mu-\tau)}$ represents the diagonalizing matrix of $\rm M_{\nu}$ and $\rm{M_{\nu}}^{I}$. The Majorana phases have been taken in the type I seesaw term ndbddb1 . From equation(42) we can evaluate $\rm{M_{\nu}}^{II}$. We have considered the case when $\rm{M_{\nu}}^{I}$ possess $\rm\mu-\tau$ symmetry, with the various choices for mixing matrices such as TBM, BM, HM, GRM, with uniquely predicting $\rm\theta_{13}=0$. We have considered $\rm{M_{\nu}}^{I(diag)}=X{M_{\nu}}^{(diag)}$, where we have introduced the parameter X to describe the reltive strength of the type I and II seesaw terms. The parameter X can take any numerical value provided the two seesaw terms gives rise to correct light neutrino mass matrix. In our case, we have considered X=0.5, i.e., equal contributions from both the seesaw terms. The required correction to $\mu-\tau$ type $\nu$ mass matrix for generation of non zero reactor mixing angle ($\rm\theta_{13}$) can be obtained from the perturbation matrix, $\rm{M_{\nu}}^{II}$ mass matrix. $\rm{M_{\nu}}^{II}$ can be constructed as,

[TABLE]

It can be derived using equation (41). The type II seesaw mass matrix is evaluated in terms of light neutrino mass matrix, constructed using the best fit neutrino data and $\mu-\tau$ symmetric type I mass matrices (TBM, BM, HM, GRM). The elements are shown in appendix.

To evaluate $\rm{m_{N+\Delta_{R}}}^{eff}$, we need the diagonalizing matrix of the heavy right handed Majorana mass matrix $\rm M_{RR}$, $U_{Rei}$ and its mass eigenvalues, $\rm M_{i}$.

$\rm M_{RR}$ can be written in the form(from reference newphysics ) and is evident from equation (21)

[TABLE]

[TABLE]

In the above equation, $\rm U_{(\mu-\tau)}$ represents $\rm U_{TBM}$, $\rm U_{BM}$, $\rm U_{HM}$, $\rm U_{GRM}$ mutau , i.e, the diagonalizing matrices of the TBM, BM, HM and GRM mass matrices. For TeV scale type I + type II seesaw, we have fine tuned the dimensionless parameter, $\gamma\sim 10^{-10}$, we have considered $v_{R}\sim TeV$. Thus after obtaining $\rm M_{RR}$, we diagonalized it and obtained the eigenvalues, $\rm M_{i}$ and its diagonalizing matrix in terms of the lightest neutrino mass ( $\rm m_{1}$ or $\rm m_{3}$)for (NH/IH) and the Majorana phases ($\rm\alpha$ and $\beta$). We have varied the Majorana phases $\rm\alpha$ and $\rm\beta$ from 0 to $2\pi$ and evaluated the effective mass for new physics contribution using formula (39) in terms of lightest neutrino mass. This is shown in figure (4).We have imposed the KamLAND-Zen bound on the new physics contribution to effective mass and the Planck bound on the sum of the absolute neutrino mass.

V.3 New physics contribution to NDBD considering perturbation in type I seesaw.

Alternatively, we have again considered the type II seesaw to give rise to $\rm\mu-\tau$ type neutrino mass matrix and the necessary correction to obtain non-zero $\rm\theta_{13}$ is obtained from the type I seesaw term. Thus, $\rm{M_{\nu}}^{II}$ in equation (44) can be written as,

[TABLE]

where, $\rm U_{(\mu-\tau)}$ represents $\rm U_{TBM},U_{BM},U_{HM},U_{GRM}$.

[TABLE]

[TABLE]

Like in the previous case, we have again evaluated the right handed Majorana mass matrix using equation (44). We have fine tuned the dimensionless parameter $\gamma$ and then by diagonalizing the right handed Majorana mass matrix $\rm M_{RR}$, we have obtained $\rm U_{Rei}$ and the eigenvalues, $\rm M_{i}$ $\rm(i.e.{M_{RR}}^{(diag)})$ where,

[TABLE]

We then evaluated the effective Majorana mass, $\rm{m_{N+\Delta_{R}}}^{eff}$ using equation (39) as a function of the lightest left handed neutrino mass. This is shown in figure (5). When we consider the type II seesaw term to be $\rm\mu-\tau$ symmetric and the perturbation from the type I seesaw term, the type I seesaw mass matrix can be derived as in the previous case and is shown in appendix.

For the new physics contribution in which the type II term acts as the perturbation, we have also evaluated the half life of the $0\nu\beta\beta$ decay process using equation (2), where

[TABLE]

By substituting the values of the phase factors(${G_{0}}^{\nu}$) phasefactor ndbdg , nuclear matrix element(NME) nme ndbdg and mass of electron, we have obtained the half life as a function of the lightest mass in the different mixing patterns for both NH and IH, as shown in figure (6). In the similar process, we have also computed the half life for new physics contribution to NDBD in which the type I term acts as the perturbation, for generation of non zero $\theta_{13}$. It is shown in figure (7).

V.4 Correlating LFV with lightest neutrino mass and $\rm\theta_{23}$

To correlate LFV with neutrino mass in our analysis, we have considered the LFV processes, $\rm\mu\rightarrow 3e$ and $\rm\mu\rightarrow e\gamma$. The BR for both the processes have a strong flavour dependence on the RH mixing matrix. Since the process $\rm\mu\rightarrow 3e$ is controlled by $\rm h_{\mu e}h_{ee}^{*}$ whereas $\rm\mu\rightarrow e\gamma$ is controlled by the factor $\rm\left[M_{R}{M_{R}}^{*}\right]_{\mu e}$ , the later is independent of the Majorana CP phases and the lightest neutrino mass, $\rm m_{j}$. We have correlated the BR of the process $\rm\mu\rightarrow 3e$ with the lightest neutrino mass $\rm(m_{1}/m_{3})$ for (NH/IH). The BR of the process $\rm\mu\rightarrow e\gamma$ is correlated with the atmospheric mixing angle, $\rm\theta_{23}$, since the other two mixing angles $\rm\theta_{12}$ and $\rm\theta_{13}$ are measured precisely. For calculating the BR, we used the expression given in equation (30) and (32). The lepton Higgs coupling $\rm h_{ij}$ in (31) can be computed explicitly for a given RH neutrino mass matrix as shown in equation (44) by diagonalizing the RH neutrino mass matrix and obtaining the mixing matrix element, $V_{i}$ and the eigenvalues $M_{i}$. For evaluating $M_{RR}$, we need to know $\rm{M_{\nu}}^{II}$, as evident from equation (44). We computed ${M_{\nu}}^{II}$ from equation (42). For determining the BR for $\mu\rightarrow 3e$, we imposed the best fit values of the parameters, $\rm{\Delta m_{sol}}^{2}$, $\rm{\Delta m_{atm}}^{2}$, $\rm\delta$, $\rm\theta_{13}$, $\rm\theta_{23}$, $\rm\theta_{12}$ in $\rm M_{\nu}$ . The numerical values of $\rm{M_{\nu}}^{I}$ can be computed as before for different mixing patterns, TBM, BM, HM, GRM. Thus, we get $\rm{M_{\nu}}^{II}$ as a function of the parameters $\rm\alpha,\beta$ and $\rm m_{lightest}$. Then varying both the Majorana phases, $\rm\alpha,\beta$ from 0 to 2$\pi$, we obtained $\rm{M_{\nu}}^{II}$ as a function of $m_{lightest}$. Similarly, for $\mu\rightarrow e\gamma$ we substituted the values of the lightest mass (m1/m3)for(NH/IH) as (0.07eV/0.065eV) and best fit values for the parameters ${\Delta m_{sol}}^{2}$, ${\Delta m_{atm}}^{2}$, $\rm\delta$, $\rm\theta_{13}$, while varying both the Majorana phases, $\alpha,\beta$ from 0 to 2$\pi$ and thus obtained $\rm{M_{\nu}}^{II}$ and hence $\rm M_{RR}$ as a function of the atmospheric mixing angle $\theta_{23}$. Thus BR can be obtained as a function of $\sin^{2}\theta_{23}$ from equation (32). We have varied the value of $\rm\sin^{2}\theta_{23}$ in its 3$\sigma$ range as in table 1 and the lightest neutrino mass from $10^{-3}$ to $10^{-1}$ and obtained the values of BR for different mixing patterns, TBM, BM, HM, GRM. The variation is shown in figure (8),(9),(10) and (11) for both NH and IH.

VI CONCLUSION

The quest for NDBD and its interrelation with neutrino mass makes it a very interesting and enthralling topic of research at present time. Its existence would not only confirm the intrinsic nature of the neutrinos but also would provide a stringent limit on the absolute scale of the neutrino mass. In this paper, we contemplated the implications of NDBD and LFV in LRSM framework. Owing to the presence of new scalars and gauge bosons in this model, various additional sources would give rise to contributions to NDBD process, which involves RH neutrinos, RH gauge bosons, scalar Higgs triplets as well as the mixed LH-RH contributions. For a simplified analysis we have ignored the left-right gauge boson mixing and heavy light neutrino mixing. We have considered the extra gauge bosons and scalars to be of the order of TeV. Again the existence of non zero $\rm\theta_{13}$ has many implications in neutrino sector beyond SM. A simple way to accomodate non zero $\rm\theta_{13}$ is by adding a perturbation matrix to the neutrino mass matrix. A well known neutrino mass mixing pattern is the one obeying $\rm\mu-\tau$ symmetry. In our present analysis, we have considered the different realizations of the $\mu-\tau$ symmetric mass matrices, namely, TBM, BM, HM, GRM matrices. The perturbation to this matrices to generate non zero $\theta_{13}$ is obtained from either of the seesaw terms, type I and type II. We have considered two different approaches, type I giving $\rm\mu-\tau$ symmetry and type II as perturbation, type II giving $\rm\mu-\tau$ symmetry and type I as perturbation, for generation of non zero $\rm\theta_{13}$. We analysed the standard as well as new physics contribution to the effective mass $\rm m_{eff}$ governing NDBD as well as the half life considering both type I and type II seesaw. We have shown the variations of the effective mass as well as the half life with the lightest neutrino mass which corresponds to the standard as well as the non standard contributions. We have seen from our analysis that both the approaches yields different consequences in NDBD. The various parameters we have chosen for our numerical analysis are consistent with constraints from $\nu$ oscillation experiments. We have also discussed the impacts of the lightest neutrino mass and not so precisely known atmospheric mixing angle, $\theta_{23}$ on the behaviour of LFV of the decay process, $\rm\mu\rightarrow 3e$ and $\rm\mu\rightarrow e\gamma$ respectively. Based on our observations, the following conclusions could be arrived at,

• •

In the standard light neutrino contribution to NDBD, it is observed that all the mass patterns (TBM, HM, GRM) yields almost similar results for NH mass spectrum. The effective mass governing NDBD is found to be of the order of $\rm 10^{-3}$ eV and are within and much below the current experimental limit kamland . Whereas in case of IH mass spectrum, for TBM, HM and GRM, the values of effective mass are found to be within and close to the experimental limit and are of the order of $10^{-2}$ eV. However, in all the cases, the light neutrino contribution can saturate the experimental limit for lightest neutrino mass ($\rm m_{1}/m_{3}$) for (NH/IH) of around 0.1 eV.

• •

In new physics contribution considering perturbation in type II seesaw, for IH, TBM, HM and GRM shows results within the recent experimental bound for lightest mass varying from (0.001-0.1) eV. Whereas, for NH the effective mass lies within experimental limit for lightest mass in the range (0.01-0.1) eV. In case of half life also, except BM mass pattern, TBM, HM and GRM schemes shows better results. In all the cases, both NH and IH seems to be more compatible with the experimental results.

• •

In new physics contribution considering perturbation in type I seesaw, the values that are consistent with experimental bound imposed by KamLAND-Zen are found for lightest mass (0.001-0.1) eV for TBM and about (0.01-0.1 eV) for all other cases. Whereas for half life, TBM shows better results. In all other mixing patterns, half life lies within experimental bound for values of lightest mass lying from (0.005-0.1) eV for IH.

• •

It is observed from our analysis that the BR for the process $\rm\mu\rightarrow 3e$ in the LRSM remains consistent with the experimental bound for a wide range of light neutrino mass. However, it depends on the neutrino mass spectrum as evident from fig 8 and 9. In case of IH, the BR is spread over a wide range and lies even in the range of the recently proposed limit with a sensitivity of $10^{-16}$. For the process, $\rm\mu\rightarrow e\gamma$, the results for BR are found to be consistent with the experimental limit for all the mixing patterns, except for HM and BM (NH) in the $3\sigma$ range of $\theta_{23}$ . In this case, the dependence of LFV on the neutrino mass spectrum is not much significant as seen in fig 10 and 11.

The effective neutrino mass depends on the character of the neutrino mass spectrum. In most of our analysis in case of NDBD as well as LFV, we have observed that both the hierachial patterns shows almost equal dominance. However, it is easier to observe the process if we consider the leading order mass matrices obeying $\rm\mu$-$\tau$ symmetry originating from type I seesaw and using type II seesaw as perturbations to generate non zero $\rm\theta_{13}$. Nevertheless, a more detailed analysis considering the presence of all the mechanisms which can generate the process in the LRSM framework should be persued to give a general conclusion.

VII APPENDIX

Elements of the type II Seesaw mass matrix (case B) and type I Seesaw mass matrix (Case C):

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Where, $\rm{c_{ij}^{\mu\tau}}=\cos\theta_{ij}^{\mu\tau}$, $\rm{s_{ij}^{\mu\tau}}=\sin\theta_{ij}^{\mu\tau}$ represents the mixing angles for $\rm\mu-\tau$ symmetric neutrino mass matrix (TBM, BM, HM, GRM).