## Radiative and Seesaw Threshold Corrections to the S 3 Symmetric Neutrino Mass Matrix

Shivani Gupta a and C. S. Kim b

Department of Physics and IPAP, Yonsei University, Seoul 120-749, Korea

Pankaj Sharma c

Korea Institute of Advanced Study, Seoul 130-722, Korea

## Abstract

We systematically analyze the radiative corrections to the S 3 symmetric neutrino mass matrix at high energy scale, say the GUT scale, in the charged lepton basis. There are significant corrections to the neutrino parameters both in the Standard Model (SM) and Minimal Supersymmetric Standard Model (MSSM) with large tan β , when the renormalization group evolution (RGE) and seesaw threshold effects are taken into consideration. We find that in the SM all three mixing angles and atmospheric mass squared difference are simultaneously obtained in their current 3 σ ranges at the electroweak scale. However, the solar mass squared difference is found to be larger than its allowed 3 σ range at the low scale in this case. There are significant contributions to neutrino masses and mixing angles in the MSSM with large tan β from the RGEs even in the absence of seesaw threshold corrections. However, we find that the mass squared differences and the mixing angles are obtained in their current 3 σ ranges at low energy when the seesaw threshold effects are also taken into account in the MSSM with large tan β .

PACS numbers: 14.60.Pq, 14.60.St

a shivani@yonsei.ac.kr

b cskim@yonsei.ac.kr

c pankajs@kias.re.kr

## I. INTRODUCTION

The neutrino oscillation experiments have enriched our knowledge of masses and mixings of neutrinos and thus flavor structure of leptons. These developments aspire theorists to construct models for unraveling the symmetries of lepton mass matrices. For three neutrino flavors the solar and atmospheric mass squared differences are given as, ∆ m 2 12 ≈ 10 -5 eV 2 and | ∆ m 13 | 2 ≈ 10 -3 eV , where ∆ 2 m 2 ij = m 2 j -m 2 i and m i,j are the mass eigenvalues of neutrinos, respectively. With the evidence of nonzero value of reactor mixing angle θ 13 [1] we now have information of all three mixing angles contrary to earlier studies where only an upper bound on θ 13 existed. The lepton flavor mixing matrix comprising of three mixing angles and a CP violating phase is given as

$$\Lambda \, = \, \left ( \begin{array} { c c c c } c _ { 1 2 } c _ { 1 3 } & s _ { 1 2 } c _ { 1 3 } & s _ { 1 3 } e ^ { - i \delta _ { C P } } \\ - s _ { 1 2 } c _ { 2 3 } - c _ { 1 2 } s _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & c _ { 1 2 } c _ { 2 3 } - s _ { 1 2 } s _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & s _ { 2 3 } c _ { 1 3 } \\ s _ { 1 2 } s _ { 2 3 } - c _ { 1 2 } c _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & - c _ { 1 2 } s _ { 2 3 } - s _ { 1 2 } c _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & c _ { 2 3 } c _ { 1 3 } \end{array} \right ) \,.$$

Here c ij = cos θ ij , s ij = sin θ ij ; θ ij are the three mixing angles, δ CP is the Dirac CP phase. There are two additional CP phases if neutrinos are Majorana particles. The best fit values along with their 3 σ ranges of neutrino oscillation parameters [2] are shown in Table I. Dirac phase δ CP , which contributes to CP violation in the leptonic sector, is expected to be measured in the long baseline neutrino oscillation experiments. The strength of this leptonic CP violation is parameterized by Jarlskog rephasing invariant [3] J = c 12 s 12 c 23 s 23 c 2 13 s 13 sin δ CP . The two Majorana phases, however, contribute to the lepton number violating processes like neutrinoless double beta decay. The possible measurement of effective Majorana mass in the current and upcoming experiments like GERDA, EXO, CUORE, MAJORANA, SuperNEMO [4] will provide some additional constraints on the two Majorana phases and neutrino mass scale. The cosmological constraint on the sum of neutrino masses by the Planck Collaboration [5] is Σ m ν i &lt; 0.23 eV at 95% C.L.. Depending on the values chosen for the priors this sum can be in the range (0 23 . -0 933) eV. An additional bound from the Solar . Pole Telescope Collaboration [6] gives Σ m ν i = 0.32 ± 0.11 eV. However, these bounds are expected to be tested in forthcoming observations. There are some challenges left namely to determine the absolute mass scale, mass hierarchy of neutrinos and the CP violation in leptonic sector amongst others.

Evidence for nonzero θ 13 has lead to many studies for the deviation from the assumed

TABLE I. The experimental constraints on neutrino parameters taken from [2].

| Parameter                         | Best fit   | 3 σ              |
|-----------------------------------|------------|------------------|
| ∆ m 2 12 / 10 - 5 eV 2 (NH or IH) | 7.54       | 6.99 - 8.18      |
| ∆ m 2 13 / 10 - 3 eV 2 (NH)       | 2.43       | 2.19 - 2.62      |
| θ 12                              | 33.64 ◦    | 30.59 ◦ - 36.8 ◦ |
| θ 13                              | 8.93 ◦     | 7.47 ◦ -10.19 ◦  |
| θ 23                              | 37.34 ◦    | 35.1 ◦ - 52.95 ◦ |

symmetries that predict the vanishing θ 13 value. Among many possible discrete flavor symmetries to produce the current data, S 3 has been extensively studied in literature [7]. It is the smallest discrete non-Abelian group which is the permutation of three objects. Perturbations to S 3 symmetric leptonic mass matrices have been used to study the mass spectra of the leptons and predict well known democratic [8] and tri-bimaximal neutrino mixing scenario [9]. The effective light neutrino mass matrix, M , invariant under ν S 3 is given as

$$M _ { \nu } = p I + q D,$$

where

$$I = \left ( \begin{array} { c c c } 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{array} \right ), \quad D = \left ( \begin{array} { c c c } 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{array} \right ).$$

Here D is democratic matrix, p and q are in general complex parameters. There are quite a few studies about the breaking of S 3 symmetry in the leptonic sector [10-12] to produce the current observed neutrino oscillation data. It will be interesting to study order of breaking in the S 3 symmetric neutrino mass matrix that arise due to the RGE and the seesaw threshold corrections between the GUT scale Λ g and the electroweak scale Λ ew in the SM and MSSM. In the light of non zero θ 13 our main aim is to study the viability of producing the current neutrino oscillation data at Λ ew scale in an S 3 symmetric neutrino mass matrix M ν at Λ g when both radiative and seesaw threshold corrections are taken into account.

Earlier studies have shown that there are significant RGE corrections to neutrino masses and mixing angles particularly for quasi-degenerate neutrino spectrum in the MSSM with large tan β . The SM is extended by three heavy right handed neutrinos at high energy scale to generate neutrino masses in Type I seesaw mechanism [13]. The seesaw threshold corrections arise due to subsequent decoupling of these heavy right handed Majorana neutrinos at their respective masses. The structure of the Dirac mass matrix M D is proportional to the neutrino Yukawa coupling matrix Y ν . We take a general Y ν and scan the parameter space to obtain the desired mixing pattern. The right handed Majorana mass matrix M is found by inverting the Type I seesaw formula at Λ . The charged lepton mass matrix g M l at Λ g is taken to be diagonal. We study the running behavior of neutrino masses and mixings from Λ g down to Λ ew using the RGE for the Yukawa couplings in the S 3 symmetric M ν both in the SM and MSSM. Above the heaviest seesaw scale ( M 3 ) there is a full theory and thus RGEs for Yukawa couplings Y e , Y ν and mass matrix M are considered. However, since our right handed neutrino mass matrix M is hierarchical ( M &lt; M 1 2 &lt; M 3 ), we also consider the seesaw threshold effects and thus the respective set of effective theories in between these scales, arising from the subsequent decoupling of heavy right handed fields at their respective masses. In the SM we find that the atmospheric mass squared difference (∆ m 2 13 ) along with the neutrino mixing angles are generated in their present 3 σ ranges at the low energy scale. However, the solar mass squared difference (∆ m 2 12 ) is greater than its allowed value ( ≈ 10 -4 ) in the SM. We find that it is possible to radiatively generate the current neutrino masses and mixing angles from the S 3 invariant neutrino mass matrix M ν in the charged lepton basis, when the seesaw threshold effects are taken into account in the MSSM with large tan β . In the MSSM with large tan β the solar mass squared difference ∆ m 2 12 can be produced in its current range along with the other neutrino oscillation parameters at Λ ew in the presence of these threshold corrections.

In section II we give the form of lepton mass matrices considered at Λ . In the subsequent g section, we give the RGE equations governing from Λ g to Λ ew , in presence of the seesaw threshold effects both in the SM and MSSM. In section IV we study the order of corrections to the neutrino mass matrix in the presence of seesaw threshold effects. Section V gives our numerical results for both cases under consideration. We conclude in the last section.

## II. FORM OF LEPTON MASS MATRICES AT THE GUT SCALE

We consider the basis where charged lepton mass matrix ( M l ) is diagonal and the effective light neutrino mass matrix ( M ν ) is S 3 symmetric as given in Eq. (2). The Yukawa coupling matrix for charged leptons is given as

$$Y _ { e } = \frac { 1 } { v } \left ( \begin{array} { c c c } m _ { e } & 0 & 0 \\ 0 & m _ { \mu } & 0 \\ 0 & 0 & m _ { \tau } \end{array} \right ),$$

where the Higgs vacuum expectation value (VEV) v is taken to be 246 GeV in the SM and 246 cos · β GeV in the MSSM. The Yukawa coupling matrix Y ν for the light neutrinos is taken of the form Y ν = y U D ν ν as given in [14] where D is the diagonal matrix Diag( r , r 1 2 , 1). The three parameters y ν , r 1 and r 2 are real, positive dimensionless parameters that characterize eigenvalues of Y ν . The unitary matrix U ν is the product of the three rotation matrices R 23 ( ϑ 2 ) · R 13 ( ϑ e 3 -iδ ) · R 12 ( ϑ 1 ) having one CP violating phase δ . Thus, Y ν has seven unknown parameters viz. three eigenvalues, three mixing angles and one CP phase. We vary the three hierarchy ( y ν , r 1 , r 2 ) parameters and, though they are completely arbitrary, but assumed to be &lt; O (1). Three angles ϑ 1 , ϑ 2 , ϑ 3 and δ are varied in the range of (0-2 π ).

The right handed mass matrix M is found by inverting the Type I seesaw formula as

$$M = - Y _ { \nu } M _ { \nu } ^ { - 1 } Y _ { \nu } ^ { T },$$

i.e. M in our study is found from Y ν and the S 3 symmetric neutrino mass matrix M ν at Λ g by inverting seesaw formula. The three right handed neutrino masses M 1 , M 2 and M 3 are obtained by diagonalizing the right handed Majorana mass matrix M . The light neutrino mass matrix M ν can be diagonalized by the unitary transformation R as R M R T ν . One of the possible form of R can be

$$R = U _ { T B M } = \left ( \begin{array} { c c c } \frac { 2 } { \sqrt { 6 } } & \frac { 1 } { \sqrt { 3 } } & 0 \\ - \frac { 1 } { \sqrt { 6 } } & \frac { 1 } { \sqrt { 3 } } & - \frac { 1 } { \sqrt { 2 } } \\ - \frac { 1 } { \sqrt { 6 } } & \frac { 1 } { \sqrt { 3 } } & \frac { 1 } { \sqrt { 2 } } \end{array} \right ).$$

The mass eigenvalues of M ν are p , p +3 q and p : corresponding to the light neutrino masses m 1 , m 2 and m 3 , respectively. Due to the degeneracy in the mass eigenvalues m 1 and m 3 , the diagonalizing matrix R is not unique. Degeneracy of masses implies that R is arbitrary up to

orthogonal transformation R 13 ( φ ), where φ is in 1-3 plane. Thus, most general diagonalizing matrix R is U TBM R 13 ( φ ), which implies the same physics as U TBM . In this work we set φ =0 without loss of generality [10, 12]. From the neutrino oscillation data we know ∆ m 2 12 ≈ 10 -5 and thus, there is small difference in the mass eigenvalues m 1 and m 2 which is a possible objection to this scenario of S 3 invariant approximation as here m 1 and m 3 are degenerate. The possible solution to this problem is given in [10, 12] where the complex numbers p and p +3 q are considered to have same magnitude but different directions. As shown in [10], q can be chosen completely imaginary and p is taken to be | p e | -i α 2 . The magnitudes of p and q can be written in terms of parameter x as

$$| p | = x \ s e c \frac { \alpha } { 2 }, \\ | q | = \frac { 2 } { 3 } x \ t a n \frac { \alpha } { 2 },$$

where x is a real free parameter and allowed range of α is 0 ≤ α &lt; π . The magnitude of p and p +3 q can be made equal by adjusting the phase α . The parameter x vanishes when α =180 ◦ and thus this value is disallowed. Substituting the values of p and q given in Eq. (7), the magnitudes of the mass eigenvalues are given as

$$| m _ { 1 } | = | m _ { 2 } | = | m _ { 3 } | = x \ s e c \frac { \alpha } { 2 }.$$

This results in equal magnitude of all three mass eigenvalues and thus a degenerate spectrum of neutrinos to begin with at Λ . g As pointed out earlier in [10] that the phase α affects the rate of neutrinoless double beta decay but will not affect neutrino oscillation parameters. Thus, this phase is of Majorana type. When we run these masses from Λ g to highest seesaw scale M 3 , the degeneracy of the mass eigenvalues is lifted by RGE corrections. We consider the normal hierarchical spectrum of masses where m 1 is the lowest mass. The other two masses are given as m 2 = √ m 2 1 +∆ m 2 12 and m 3 = √ m 2 1 +∆ m 2 13 . Since the three mass eigenvalues at Λ g have equal magnitude the two mass squared differences are vanishing to begin with. Once the degeneracy of three mass eigenvalues is lifted, their non zero values are generated. In subsequent sections we will explore the generation of the solar and atmospheric mass squared differences, together with the three mixing angles in their current 3 σ limit at Λ ew through the radiative corrections from S 3 invariant neutrino mass matrix at Λ . g

## III. RGE EQUATIONS WITH SEESAW THRESHOLD EFFECTS

In Type I seesaw the SM is extended by introducing three heavy right handed neutrinos and keeping the Lagrangian of electroweak interactions invariant under SU (2) L × U (1) Y gauge transformation. In this case, the leptonic Yukawa terms of the Lagrangian are written as

$$- \, \mathcal { L } _ { \nu } = \bar { l } _ { e } \phi Y _ { e } e _ { R } + \bar { l } _ { e } \tilde { \phi } Y _ { \nu } \nu _ { R } + \frac { 1 } { 2 } \bar { \nu } _ { R } ^ { c } M \nu _ { R } + h. c..$$

Here φ is the SM Higgs doublet and ˜ = φ iσ φ 2 ∗ . e R and ν R are right handed charged lepton and neutrino singlets, Y e and Y ν are the Yukawa coupling matrices for charged leptons and Dirac neutrinos, respectively. The last term of Eq. (9) is the Majorana mass term for the right handed neutrinos.

Quite intensive studies have been done in literatures [15, 16] regarding the general features of RGE of neutrino parameters. At the energy scale below seesaw threshold i.e. when all the heavy particles are integrated out, the RGE of neutrino masses and mixing angles is described by the effective theory which is same for various seesaw models. But above the seesaw scale full theory has to be considered and thus, there can be significant RGE effects due to the interplay of heavy and light sector. The RGE equations and subsequent decoupling of heavy fields at their respective scales are elegantly given in Refs. [17-19]. The comprehensive study of the RGE and seesaw threshold corrections to various mixing scenarios is recently done in [20].

The effective neutrino mass matrix M ν above M 3 is given as

$$M _ { \nu } ( \mu ) = - \frac { v ^ { 2 } } { 2 } Y _ { \nu } ^ { T } ( \mu ) M ^ { - 1 } ( \mu ) Y _ { \nu } ( \mu ),$$

where v = 246 sin · β GeV in the MSSM and µ is the renormalization scale. Y ν and M are µ dependent. Since we study the evolution of leptonic mixing parameters from Λ g to Λ ew scale in a generic seesaw model we need to take care of the series of effective theories that arise by subsequent decoupling of the heavy right handed fields M i ( i =1,2,3) at their respective mass thresholds. The Yukawa couplings Y ν and M are dependent on the energy scale Λ. At the GUT scale we consider the full theory and the one loop RGEs for Y e , Y ν and M are given as

$$\dot { Y } _ { e } = \frac { 1 } { 1 6 \pi ^ { 2 } } Y _ { e } \left [ \alpha _ { e } + C _ { 1 } H _ { e } + C _ { 2 } H _ { \nu } \right ],$$

$$\dot { Y } _ { \nu } & = \frac { 1 } { 1 6 \pi ^ { 2 } } Y _ { \nu } \left [ \alpha _ { \nu } + C _ { 3 } H _ { e } + C _ { 4 } H _ { \nu } \right ], \\ \dot { M } & = \frac { 1 } { 1 6 \pi ^ { 2 } } C _ { 5 } \left [ ( Y _ { \nu } Y _ { \nu } ^ { \dagger } ) M + M ( Y _ { \nu } Y _ { \nu } ^ { \dagger } ) ^ { T } \right ],$$

where ˙ Y i = dY i dt ( = i e, ν ), t= ln µ/µ ( 0 ) with µ µ ( 0 ) being the running (fixed) scale, and H i = Y † i Y i ( i = e, ν ). The coefficients are C 1 = 3 2 , C 2 = -3 2 , C 3 = -3 2 , C 4 = 3 2 , C 5 = 1 in the SM and C 1 = 3 , C 2 = 1 , C 3 = 1 , C 4 = 3 , C 5 = 2 in the MSSM, respectively. The expressions for α e and α ν in the SM and MSSM are explicitly given as

$$\text{nd} \, \alpha _ { \nu } \text{ in the SM and MSSM are explicitly given as} \\ \alpha _ { e ( S M ) } & = T r ( 3 H _ { u } + 3 H _ { d } + H _ { e } + H _ { \nu } ) - ( \frac { 9 } { 4 } g _ { 1 } ^ { 2 } + \frac { 9 } { 4 } g _ { 2 } ^ { 2 } ), \\ \alpha _ { \nu ( S M ) } & = T r ( 3 H _ { u } + 3 H _ { d } + H _ { e } + H _ { \nu } ) - ( \frac { 9 } { 2 0 } g _ { 1 } ^ { 2 } + \frac { 9 } { 4 } g _ { 2 } ^ { 2 } ), \\ \alpha _ { e ( M S S M ) } & = T r ( 3 H _ { d } + H _ { e } ) - ( \frac { 9 } { 5 } g _ { 1 } ^ { 2 } + 3 g _ { 2 } ^ { 2 } ), \\ \alpha _ { \nu ( M S S M ) } & = T r ( 3 H _ { u } + H _ { \nu } ) - ( \frac { 3 } { 5 } g _ { 1 } ^ { 2 } + 3 g _ { 2 } ^ { 2 } ), \\ \text{re} \, q _ { 1 } \, \text{ are the } U ( 1 ) _ { V } \text{ and } S U ( 2 ) _ { I } \text{ gauge coupling constants.}$$

where g 1 2 , are the U (1) Y and SU (2) L gauge coupling constants.

The heavy right handed mass matrix M obtained from Eq. (5) is non diagonal and thus is diagonalized by the unitary transformation U R as

$$U _ { R } ^ { T } M U _ { R } = D i a g ( M _ { 1 }, M _ { 2 }, M _ { 3 } ).$$

The Yukawa coupling Y ν is accordingly transformed as Y U ν ∗ R . At the highest seesaw scale M 3 , the effective operator κ (3) is given by the matching condition as

$$\kappa _ { ( 3 ) } = 2 Y _ { \nu } ^ { T } M _ { 3 } ^ { - 1 } Y _ { \nu },$$

in the basis where M is diagonal. The Yukawa coupling Y ν above is a 3 × 3 matrix and all the variables are set to the scale M 3 . At the scale lower than M 3 ( µ &lt; M 3 ) the effective neutrino mass matrix M ν is given as

$$M _ { \nu } = - \frac { v ^ { 2 } } { 4 } \{ \kappa _ { ( 3 ) } + 2 Y _ { \nu ( 3 ) } ^ { T } M _ { ( 3 ) } ^ { - 1 } Y _ { \nu ( 3 ) } \}.$$

As can be seen that M ν is the sum of κ (3) given in Eq. (14) and the seesaw factor which is obtained after decoupling M 3 . Thus, Y ν (3) is 2 × 3 and M (3) is 2 × 2 mass matrices. RGE between the scale M 3 and M 2 is governed by the running of κ (3) , Y ν (3) and M (3) . The running of κ (3) is given as

$$\dot { \kappa } _ { ( 3 ) ( S M ) } = \frac { 1 } { 1 6 \pi ^ { 2 } } \left [ ( C _ { 3 } H _ { e } + C _ { 6 } H _ { \nu ( 3 ) } ) ^ { T } \kappa _ { ( 3 ) } + \kappa _ { ( 3 ) } ( C _ { 3 } H _ { e } + C _ { 6 } H _ { \nu ( 3 ) } ) + \alpha _ { ( 3 ) } \kappa _ { ( 3 ) } \right ], \quad ( 1 6 )$$

$$\dot { \kappa } _ { ( 3 ) ( M S S M ) } = \frac { 1 } { 1 6 \pi ^ { 2 } } \left [ ( H _ { e } + C _ { 6 } H _ { \nu ( 3 ) } ) ^ { T } \kappa _ { ( 3 ) } + \kappa _ { ( 3 ) } ( H _ { e } + C _ { 6 } H _ { \nu ( 3 ) } ) + \alpha _ { ( 3 ) } \kappa _ { ( 3 ) } \right ],$$

where C 6 = 1 2 and H ν (3) = Y † ν (3) Y ν (3) . α (3) in the SM and MSSM is explicitily given as

$$\alpha _ { ( 3 ) ( S M ) } & = 2 \ T r ( 3 H _ { d } + 3 H _ { u } + H _ { e } + H _ { \nu ( 3 ) } ) - 3 g _ { 2 } ^ { 2 } + \lambda, \\ \alpha _ { ( 3 ) ( M S S M ) } & = 2 \ T r ( 3 H _ { u } + H _ { \nu ( 3 ) } ) - \frac { 6 } { 5 } g _ { 1 } ^ { 2 } - 6 g _ { 2 } ^ { 2 }.$$

The effective operator κ (2) at M 2 is given by the matching condition as

$$\kappa _ { ( 2 ) } = \kappa _ { ( 3 ) } + 2 Y _ { \nu ( 2 ) } ^ { T } M _ { ( 2 ) } ^ { - 1 } Y _ { \nu ( 2 ) }.$$

where all the parameters are set to seesaw scale M 2 . M ν at the scale below M 2 has the RGE equation given in Eq. (16). In their RGE expression we have Y ν (2) which is 1 × 3 and M (2) which is 1 × 1. The low energy effective theory operator κ (1) is obtained after integrating out all three heavy right handed fields. The one loop RGE for κ (1) from lowest seesaw scale M 1 down to Λ e scale is given as

$$\dot { \kappa } _ { ( 1 ) } = ( C _ { 3 } H _ { e } ^ { T } ) \kappa _ { ( 1 ) } + \kappa _ { ( 1 ) } ( C _ { 3 } H _ { e } ) + \alpha \kappa _ { ( 1 ) },$$

where

$$\begin{matrix} 3 g _ { 2 } ^ { z } + \lambda & \text{in the SM} \\ \\ 2 & \text{in the MSSM}. \end{matrix}$$

$$\alpha = 2 \ T r ( 3 H _ { u } + 3 H _ { d } + H _ { e } ) - 3 g _ { 2 } ^ { 2 } + \lambda \ \text{ in the SM},$$

$$\alpha = 2 \ T r ( 3 H _ { u } ) - \frac { 6 } { 5 } g _ { 1 } ^ { 2 } - 6 g _ { 2 } ^ { 2 } \text{ in } t$$

When the Higgs field gets VEV, the light neutrino mass matrix is obtained from κ (1) as M ν = κ (1) v 2 4 . We diagonalize M ν to obtain neutrino masses, mixing angles and CP phases.

## IV. NEUTRINO MASSES AND MIXINGS

The RGE above the highest seesaw scale M 3 depends on more parameters than below the lowest seesaw scale M 1 due to presence of the neutrino Yukawa couplings ( Y ν ). The RGE equations consist of H e , M and H ν out of which latter can be large. In the basis where charged lepton mass matrix is diagonal, M ν at two different energy scales Λ ew and Λ g are homogeneously related as [18, 21, 22]

$$M _ { \nu } ^ { \Lambda _ { e w } } = I _ { K } \cdot I ^ { T } \cdot M _ { \nu } ^ { \Lambda _ { g } } \cdot I.$$

Here I K is flavor independent factor arising from gauge interactions and fermion antifermion loops. It does not influence the mixing angles. The matrix I has the form

$$I & = D i a g ( e ^ { - \Delta _ { e } }, e ^ { - \Delta _ { \mu } }, e ^ { - \Delta _ { \tau } } ), \\ & \succeq D i a g ( 1 - \Delta _ { e }, 1 - \Delta _ { \mu }, 1 - \Delta _ { \tau } ) + \mathcal { O } ( \Delta _ { e, \mu, \tau } ^ { 2 } ),$$

where

$$\Delta _ { j } = \frac { 1 } { 1 6 \pi ^ { 2 } } \int [ 3 ( H _ { j } ) - ( H _ { \nu _ { j } } ) ] d t,$$

where j = e, µ, τ . Numerically, ∆ SM τ can be of the order of 10 -3 when Y τ ∼ . 01 and Y ν = 0 2 and the scale . µ and µ 0 are 10 12 and 10 , 2 respectively. In the MSSM, ∆ MSSM τ ∼ 10 -3 (1 + tan 2 β ) for the same values of Y ν and Y τ . In the absence of seesaw threshold effects, ∆ τ is small ≈ 10 -5 in the SM for the above mentioned values of Y τ .

At above seesaw scale appreciable deviations may occur only for large values of Y τ or Y ν . In the absence of seesaw threshold effects i.e. when there is no H ν term in Eq. (22) the radiative corrections are governed by ∆ τ term as ∆ e,µ is too small. On the other hand, due to the presence of large H ν , the seesaw threshold effects play a crucial role in the running of mixing angles. Below the seesaw scales, deviations are obtained from Y τ ∼ √ 2 m /v τ ≈ O (10 -2 ) in the SM, and ∼ √ 2 m τ v cos β in the MSSM. There can be significant deviations in the MSSM with large tan β which enhances Y τ . In the presence of seesaw threshold, appreciable deviations are possible as it is natural to have Y ν ∼ O (1). The analytic expressions for the running neutrino mixing angles, masses and CP phases are quite long and have been earlier derived in literatures [19, 21, 23]. As shown earlier in [16, 24], the mass squared differences in the denominator of the RGEs play a significant role in the evolution of neutrino mixing angles and phases, if the left-handed neutrinos are nearly degenerate.

## V. NUMERICAL RESULTS AND DISCUSSIONS

At the GUT scale Λ g we have seven free parameters in Y ν and two free parameters α and x in M ν . The three mixing angles in ϑ 1 , ϑ 2 and ϑ 3 and phase δ are allowed to take the values in the range (0-2 π ). The hierarchy parameters y ν , r 1 , r 2 and x are randomly varied and are expected to be &lt; O (1). The physical range of phase α is from (0π ). The mass spectrum at the high scale is degenerate and thus we have vanishing solar and atmospheric mass squared difference to begin with. The parameter space at the Λ g with which the low energy neutrino

TABLE II. Numerical values of input and output parameters that are radiatively generated via the RGE and seesaw threshold effects both in the SM and MSSM. The input parameters are taken at Λ g = 2 × 10 16 GeV and tan β =55 in the MSSM.

| Parameters        | SM             | MSSM              |
|-------------------|----------------|-------------------|
| Input             |                |                   |
| r 1               | 7 × 10 - 3     | 3 × 10 - 4        |
| r 2               | 0.33           | 0.34              |
| δ                 | 158.1 ◦        | 96.3 ◦            |
| y ν               | 0.72           | 0.49              |
| θ 1               | 225.2 ◦        | 188.5 ◦           |
| θ 2               | 244.6 ◦        | 245.8 ◦           |
| θ 3               | 345.5 ◦        | 241.7 ◦           |
| x ( eV )          | 7.46 × 10 - 2  | 1 . 9 × 10 - 3    |
| α                 | 102.5 ◦        | 170.2 ◦           |
| Output            |                |                   |
| m 1 ( eV )        | 9 . 4 × 10 - 2 | 1 . 47 × 10 - 2   |
| θ 12              | 34.3 ◦         | 35.2 ◦            |
| θ 13              | 7.88 ◦         | 9.85 ◦            |
| θ 23              | 47.7 ◦         | 47.4 ◦            |
| ∆ m 2 12 ( eV 2 ) | 4.18 × 10 - 4  | 7 . 14 × 10 - 5   |
| ∆ m 2 13 ( eV 2 ) | 2.48 × 10 - 3  | 2 . 35 × 10 - 3   |
| M 1 (GeV)         | 3.1 × 10 5     | 2 . 86 × 10 3     |
| M 2 (GeV)         | 4.43 × 10 8    | 9 × 10 8          |
| M 3 (GeV)         | 2.83 × 10 9    | 4 . 7 × 10 9      |
| J                 | -2.99 × 10 - 2 | - 3 . 87 × 10 - 2 |
| | m ee | ( eV )   | 9.2 × 10 - 2   | 1 . 03 × 10 - 2   |

data is obtained at the Λ ew is illustrated in Table II. The set of input parameters in that particular parameter space are also given in the table. We choose the set of parameters in parameter space at the high scale for which maximum value of θ 13 is obtained and the other mixing angles and mass squared differences are simultaneously obtained in current 3 σ range at Λ ew . However, the parameter space under consideration is only for illustration and not unique. Search for complete parameter space is an elaborate study and thus independent future work.

## A. Radiative Threshold corrections in the SM

Study of radiative corrections to S 3 symmetric neutrino mass matrix can be divided into three regions that are governed by different RGE equations. The first region is above the highest seesaw scale M 3 to Λ , g where there can be considerable contribution of Y ν . The second region is the region between the three seesaw scales and the third is below the lowest seesaw scale M 1 , where all heavy fields are decoupled. The solar mixing angle θ 12 can have large RGE corrections as the running is enhanced by the factor proportional to m 2 ∆ m 2 12 at the leading order, which can be large for degenerate spectrum. The RGE is comparatively small for other two mixing angles θ 23 and θ 13 , where the RGE is proportional to m 2 ∆ m 2 13 . However, for the degenerate neutrino mass spectrum there can be considerable corrections for these mixing angles, too. Below the seesaw scale the RGE corrections to the mixing angles in the SM are negligible as they get contributions only from Y τ . In Fig 1, we show the RGE corrections to the mixing angles and masses in the SM for the set of input parameters given in second column of Table II. Fig. 1 shows that below M 1 scale there is no significant corrections to the mixing angles. Below the lowest seesaw scale the running of the mass eigenvalues is significant even in the SM for degenerate as well as hierarchical neutrinos [16] due to the factor α given in Eq. (19), which is much larger than Y 2 τ . The running of masses is given by a common scaling of the mass eigenvalues [25]. Clearly, the RGE of each mass eigenvalue is proportional to the mass eigenvalue itself. The running of masses in Fig. 1 can be seen to start from the degenerate values of masses at Λ g and there are significant corrections to the masses below M 1 . Earlier analyses [24] studied the successful generation of mass squared differences and mixing angles for degenerate neutrinos in the SM. In their analysis it is shown that the generation of mass squared differences is very sensitive

to the value of sin 2 θ 12 , which should be greater than 0.99 to fit the mass squared differences simultaneously with the consistent angles. This limit is, however, completely ruled out by the current oscillation data. The degeneracy of three mass eigenvalues is lifted by the RGE running from Λ g to M 3 . Potentially significant breaking of neutrino mass degeneracy is provided by the RGE effects. The seesaw threshold effects in addition increase the mass splitting between the masses m 2 and m 3 required to fit the masses with the current data in

50

0.12

40

30

20

10

0

0.115

0.11

0.105

0.1

0.095

0.09

(deg)

ij

θ

0

2

4

6

8

10

θ

θ

θ

12

14

16

0

2

4

6

8

Log10(

µ

/GeV)

Log10(

µ

10

/GeV)

FIG. 1. The RGE of the mixing angles, masses and mass squared differences between the GUT scale Λ g and the electroweak scale Λ ew in the SM. The initial values of the parameters are given in the second column of Table II. The boundaries of three grey shaded areas, i.e. dark, medium and light denote the points when heavy right handed singlets M 3 , M 2 and M 1 are integrated out respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0172v1-figure-9482aff0-2841-48b7-872f-4caca71b4c3e.jpg)

13

23

12

mi (eV)

12

m1

m3

m2

14

16

terms of mass squared differences.

Running between and above the seesaw scales is modified by the contribution of Yukawa couplings Y ν . The contribution of Y ν in the RGE can result all the three mixing angles in their current allowed ranges at the EW scale in the SM. Thus, there are significant corrections to mixing angles even in the SM in the presence of the seesaw threshold effects when there is exactly the degenerate mass spectrum to begin with. As can be seen from Fig. 1, the mixing angles at the GUT scale are θ 23 = 45 , ◦ θ 12 = 35 3 . ◦ and θ 13 = 0 . For the set of ◦ parameters given in the second column of Table II we get the mixing angles in the allowed range at the electroweak scale in the SM. θ 23 is found to have values below maximality. The presence of seesaw threshold corrections and thus Y ν in the RGE equations make this possible even in the SM, as can be seen in Fig. 1. The grey shaded area in Fig. 1 and Fig. 2 illustrates the ranges of effective theories that emerge when we integrate out heavy right handed singlets. At each seesaw scales, i.e. M 1 , M 2 and M 3 , one heavy singlet is integrated out and thus (n-1) × 3 sub-matrix of Y ν remains. Therefore, the running behavior between these scales can be different from running behavior below or above these scales. Between these scales the neutrino mass matrix comprises of two terms κ ( n ) and 2 Y T ν n ( ) M -1 R n ( ) Y ν n ( ) , as given in Eq. (15). It is shown in [16] that in the SM these two terms between the thresholds are quite different which can give dominant contribution to the running of mixing angles in this region. Both θ 12 and θ 13 in Fig. 1 get large corrections between three seesaw scales. θ 23 gets the deviation of ≈ 2.7 ◦ in the upper direction.

Three neutrino mixing angles can be generated in their current allowed ranges starting from S 3 symmetric neutrino mass matrix at the Λ g in the SM as there can be large corrections when we consider the seesaw threshold effects. In the SM, atmospheric mass squared difference ∆ m 2 13 is generated within the current oscillation data limit ( ≈ 2 48 . × 10 -3 ) starting from the vanishing value at the Λ g since all masses are degenerate, as seen from Fig. 1. The solar mass squared difference ∆ m 2 12 of the order of ≈ 10 -4 eV 2 is simultaneously generated at the Λ ew , as shown in Fig.1 which is larger than its present value. The byproduct of this analysis is the masses of right-handed neutrinos that are determined from Eqs. (5) and (13) and are not free parameters. The values of effective Majorana mass | M ee | and Jarlskog rephasing invariant ( J ) at the Λ ew are also calculated for particular set of parameters given in Table II. We conclude that in the SM the RGE in addition to seesaw threshold effects cannot simultaneously produce the solar mass squared difference in its allowed range at low

scale along with the other neutrino oscillation parameters.

## B. Radiative Threshold corrections in the MSSM

As stated earlier we divide the radiative corrections to S 3 symmetric neutrino mass matrix into three regions governed by different RGE equations in the MSSM. In the region below

FIG. 2. The RGE of the mixing angles, masses and mass squared differences between the GUT scale Λ g and the electroweak scale Λ ew in the MSSM with tan β =55. The initial values of the parameters are given in the third column of Table II. The boundaries of three grey shaded areas, i.e. dark, medium and light denote the points when heavy right handed singlets M 3 , M 2 and M 1 are integrated out respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0172v1-figure-9dc0756c-9adb-4b90-b5ab-d37afbaccc29.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0172v1-figure-c2f5bf30-c9c2-4ed3-bdb7-7c2dd745033a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0172v1-figure-58a67539-e944-429a-87b5-efcfa51bb3b6.jpg)

the lightest seesaw scale, for the SM and MSSM with small tan β , Yukawa coupling Y τ ∼ 0.01 is small, and thus there cannot be large RGE corrections in these cases. For large tan β these corrections can be larger due to the presence of factor Y 2 τ (1+tan 2 β ) in the MSSM. We show the RGE of mixing angles and masses in the MSSM with tan β =55 in Fig. 2 for the set of input parameters given in third column of Table II. In the region above the energy scale M 1 , we get contributions from another Yukawa coupling Y ν which brings in more free parameters in the analysis. In the region for the MSSM with large tan β the presence of seesaw threshold effects can enhance the RGE of the mixing angles significantly. As can be seen from Fig. 2, we can have all three mixing angles simultaneously in the current limit at the EW scale starting from S 3 symmetric neutrino mass matrix at the Λ . g Fig. 2 shows that there are large corrections to θ 13 ( ∼ 9.85 ) between the Λ ◦ g and M 3 scale due to the presence of Y ν . As mentioned earlier, the running of neutrino mass matrix between the seesaw thresholds gets contributions from two terms κ ( n ) and 2 Y T ν n ( ) M -1 R n ( ) Y ν n ( ) given in Eq. (15). As shown in the SM, these two terms are different and thus results in the enhanced running between the different energy regimes. On the other hand, in the MSSM, as can be seen from Fig. 2, there are not much deviations in the mixing angles between the energy scales. It is because in the MSSM the two contributions are almost identical and thus cancel each other resulting in minimum deviation in those regions. The only significant corrections occurs in the region above the highest seesaw scale M 3 due to relatively large Y ν .

The mixing angle θ 12 does not have much corrections and θ 23 receives the correction of 2.5 ◦ in the upper direction and is thus above maximal. The running of masses in the MSSM (Fig. 2) is much larger than the SM due to the presence of tan β which in our case is large. The dominant effect, however, is the corrections in the range M 3 ≤ µ ≤ Λ g where the flavor dependent terms ( Y l and Y ν ) can be large. The interesting dependence of α ν (MSSM) and tan β on the running contributions of flavor dependent terms is given in [16]. For large tan β the contribution of Y e and Y ν become important. We also show the radiative corrections to the solar and atmospheric mass squared difference from the Λ g to the Λ ew for degenerate masses at the GUT scale in the MSSM with tan β =55. To begin with both the mass squared differences are zero at the GUT scale. From the mass squared differences shown in Fig 2, we see that the RGE in combination with seesaw threshold corrections can result both mass squared differences in their current 3 σ ranges at the low scale. The hierarchy parameters y ν , r 1 and r 2 though arbitrary but are found to be &lt; O (1). For given set of input parameters

in Table II, the value of effective Majorana mass | M ee | is ≈ 10 -2 eV and that of J is ≈ 3.87 × 10 -2 . Thus, we find that it is possible to simultaneously obtain the neutrino oscillation parameters at the electroweak scale for S 3 mass matrix at the Λ g in the MSSM with large tan β .

## VI. CONCLUSIONS

We studied the RGE corrections to the S 3 symmetric neutrino mass matrix in the presence of seesaw threshold corrections both in the SM and MSSM. In the absence of seesaw threshold effects there are negligible corrections to the mixing angles in the SM and MSSM with low tan β . However, large corrections are possible in neutrino parameters once the seesaw threshold effects are taken into consideration both in the SM and MSSM. Above all seesaw scales there is the contribution of the Yukawa coupling, Y ν in addition to Y l and thus the RGE depends on more parameters than below the seesaw scales. Thus, we can say the corrections depend on the form of Y ν which has free parameters. In the SM we found that the mixing angles can be obtained in their current 3 σ range at the electroweak scale when we begin with S 3 neutrino mass matrix at the GUT scale Λ . The significant running occurs g between and above the seesaw threshold scales. Below the lowest seesaw scale there are no significant corrections as the only contribution comes from Y τ which is small. However, in this case of exactly equal magnitude of mass eigenvalues, the two mass squared difference are not simultaneously generated in the current range at the electroweak scale. The solar mass squared difference is found to be large ( O (10 -4 )) in comparison to its allowed value at the electroweak scale.

There can be large radiative corrections in the MSSM with tan β =55 when threshold effects are taken into consideration. The large corrections to the mixing angles occur at the scale above the seesaw threshold where the Yukawa coupling, Y ν , is present and has large free parameters which can enhance running for large tan β . All the mixing angles and the mass squared differences are obtained in their current 3 σ range at the low scale in this case. The free input parameters y ν , r 1 and r 2 are expected to be small &lt; O (1). Thus, for the SM it is not possible to radiatively generate the solar mass squared difference at the electroweak scale in its current range. In the MSSM with large tan β we can generate all the masses and mixing angles in the current allowed range at the electroweak scale, starting from exactly

degenerate mass spectrum at the GUT scale.

## VII. ACKNOWLEDGEMENTS

The work of C. S. K and S. G. is supported by the National Research Foundation of Korea (NRF) grant funded by Korea government of the Ministry of Education, Science and Technology (MEST) (Grant No. 2011-0017430) and (Grant No. 2011-0020333).

- [1] K. Abe et al. [T2K Collaboration], Phys. Rev. Lett. 107 , 041801 (2011) [arXiv:1106.2822 [hepex]]; Y. Abe et al. [DOUBLE-CHOOZ Collaboration], Phys. Rev. Lett. 108 , 131801 (2012) [arXiv:1112.6353 [hep-ex]]; F. P. An et al. [DAYA-BAY Collaboration], Phys. Rev. Lett. 108 , 171803 (2012) [arXiv:1203.1669 [hep-ex]]; J. K. Ahn et al. [RENO Collaboration], Phys. Rev. Lett. 108 , 191802 (2012) [arXiv:1204.0626 [hep-ex]].
- [2] G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A. M. Rotunno, Phys. Rev. D 86 (2012) 013012 [ arXiv:hep-ph/1205.5254 ].
- [3] C. Jarlskog, Phys. Rev. Lett. 55 , 1039 (1985).
- [4] A. A. Smolnikov [GERDA Collaboration], arXiv: 0812.4194 [nucl-ex]; M. Auger et al. [EXO Collaboration], Phys. Rev. Lett. 109 , 032505 (2012); http://www.mpi-hd.mpg.de/gerda/; P. Gorla et al. [CUORE Collaboration], J. Phys. Conf. Ser. 375 , 042013 (2012) ; Wilkerson J. F et al. [MAJORANA Collaboration], J. Phys. Conf. Ser. 375 , 042010 (2012); Barabash A. S et al. [SuperNEMO Collaboration], J. Phys. Conf. Ser. 375 , 04201 (2012).
- [5] P. A. R. Ade et al. [Planck Collaboration], Planck 2013 results. XVI. Cosmological parameters , [ arXiv:astro-ph.CO/1303.5076 ].
- [6] Z. Hou, C. L. Reichardt, K. T. Story, B. Follin, R. Keisler et al. , arXiv:1212.6267 [astroph.CO].
- [7] S. Pakwasa, H. Sugawara, Phys. Lett. B 73 , 61 (1978); E. Derman, Phys. Rev. D 19, 317 (1979); E. Ma, Phys. Rev. D 61 , 033012 (2000) [hep-ph/9909249]; M. Tanimoto, Phys. Lett. B 483 , 417 (2000) [hep-ph/0001306]; P. F. Harrison and W. G. Scott, Phys. Lett. B 557 , 76 (2003) [hep-ph/0302025]; S. -L. Chen, M. Frigerio and E. Ma, Phys. Rev. D 70 , 073008 (2004) [Erratum-ibid. D 70 , 079905 (2004)] [hep-ph/0404084]; F. Caravaglios and S. Morisi,

hep-ph/0503234; W. Grimus and L. Lavoura, JHEP 0508 , 013 (2005) [hep-ph/0504153]; R. N. Mohapatra, S. Nasri and H. -B. Yu, Phys. Lett. B 639 , 318 (2006) [hep-ph/0605020]; N. Haba and K. Yoshioka, Nucl. Phys. B 739 , 254 (2006) [hep-ph/0511108]; Y. Koide, Eur. Phys. J. C 50 , 809 (2007) [hep-ph/0612058]; C. -Y. Chen and L. Wolfenstein, Phys. Rev. D 77 , 093009 (2008) [arXiv:0709.3767 [hep-ph]]; F. Feruglio and Y. Lin, Nucl. Phys. B 800 , 77 (2008) [arXiv:0712.1528 [hep-ph]]; Y. H. Ahn, S. K. Kang, C. S. Kim and T. P. Nguyen, Phys. Rev. D 82 , 093005 (2010) [arXiv:1004.3469 [hep-ph]]; D. A. Dicus, S. -F. Ge and W. W. Repko, Phys. Rev. D 82 , 033005 (2010) [arXiv:1004.3266 [hep-ph]]; Z. -z. Xing, D. Yang and S. Zhou, Phys. Lett. B 690 , 304 (2010) [arXiv:1004.4234 [hep-ph]]; D. Meloni, S. Morisi and E. Peinado, J. Phys. G 38 , 015003 (2011) [arXiv:1005.3482 [hep-ph]]; R. Jora, J. Schechter and M. N. Shahid, Phys. Rev. D 82 , 053006 (2010) [arXiv:1006.3307 [hep-ph]].

- [8] H. Fritzsch and Z. -Z. Xing, Phys. Lett. B 372 , 265 (1996) [hep-ph/9509389]; M. Fukugita, M. Tanimoto and T. Yanagida, Phys. Rev. D 57 , 4429 (1998) [hep-ph/9709388]; H. Fritzsch and Z. -z. Xing, Phys. Rev. D 61 , 073016 (2000) [hep-ph/9909304].
- [9] P. F. Harrison and W. G. Scott, Phys. Lett. B 535 , 163 (2002) [hep-ph/0203209]; P. F. Harrison and W. G. Scott, Phys. Lett. B 557 , 76 (2003) [hep-ph/0302025]; X. -G. He and A. Zee, Phys. Rev. D 68 , 037302 (2003) [hep-ph/0302201]; Prog. Theor. Phys. 109 , 795 (2003) [Erratumibid. 114 , 287 (2005)] [hep-ph/0302196].
- [10] S. Dev, S. Gupta and R. R. Gautam, Phys. Lett. B 702 , 28 (2011) [arXiv:1106.3873 [hep-ph]]; S. Dev, R. R. Gautam and L. Singh, Phys. Lett. B 708 , 284 (2012) [arXiv:1201.3755 [hep-ph]]. [11] S. Zhou, Phys. Lett. B 704 , 291 (2011) [arXiv:1106.4808 [hep-ph]].
- [12] R. Jora, S. Nasri and J. Schechter, Int. J. Mod. Phys. A 21 , 5875 (2006) [hep-ph/0605069]; R. Jora, J. Schechter and M. Naeem Shahid, Phys. Rev. D 80 , 093007 (2009) [Erratum-ibid. D 82 , 079902 (2010)] [arXiv:0909.4414 [hep-ph]]; R. Jora, J. Schechter and M. N. Shahid, Int. J. Mod. Phys. A 28 , 1350028 (2013) [arXiv:1210.6755 [hep-ph]].
- [13] P. Minkowski, Phys. Lett. B 67 , 421 (1977); T. Yanagida, (O.Sawada and A. Sugamoto, eds.), KEK, Tsukuba, Japan, 1979, p. 95; M. Gell-Mann, P. Ramond, and R. Slansky, (P. Van Nieuwenhuizen and D. Z. Freedman, eds.), North Holland, Amsterdam, 1979, p.315; R. N. Mohapatra and G. Senjanovic Phys. Rev. Lett. 44 , 912 (1980).
- [14] J. -w. Mei and Z. -z. Xing, Phys. Rev. D 70 , 053002 (2004) [hep-ph/0404081].

- [15] K. S. Babu, C. N. Leung and J. T. Pantaleone, Phys. Lett. B 319 , 191 (1993) [arXiv:9309223 [hep-ph]]; S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 519 , 238 (2001) [hep-ph/0108005]; S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 525 , 130 (2002) [hep-ph/0110366]; S. Antusch, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 538 , 87 (2002) [hep-ph/0203233]; A. S. Joshipura, S. D. Rindani and N. N. Singh, Nucl. Phys. B 660 , 362 (2003) [hep-ph/0211378]; T. Ohlsson, H. Zhang and S. Zhou, Phys. Rev. D 87 , 013012 (2013) [arXiv:1211.3153 [hep-ph]].
- [16] S. Antusch, J. Kersten, M. Lindner and M. Ratz, Nucl. Phys. B 674 , 401 (2003) [hepph/0305273].
- [17] P. H. Chankowski and Z. Pluciennik, Phys. Lett. B 316 , 312 (1993) [hep-ph/9306333]; S. Antusch, J. Kersten, M. Lindner, M. Ratz and M. A. Schmidt, JHEP 0503 , 024 (2005) [hepph/0501272]; T. Ohlsson and S. Zhou, arXiv:1311.3846 [hep-ph].
- [18] J. R. Ellis and S. Lola, Phys. Lett. B 458 , 310 (1999) [hep-ph/9904279].
- [19] J. -w. Mei, Phys. Rev. D 71 , 073012 (2005) [hep-ph/0502015].
- [20] S. Gupta, S. K. Kang and C. S. Kim, arXiv:1406.7476 [hep-ph].
- [21] S. Goswami, S. Khan and S. Mishra, arXiv:1310.1468 [hep-ph].
- [22] A. Dighe, S. Goswami and W. Rodejohann, Phys. Rev. D 75 , 073023 (2007) [hep-ph/0612328]; A. Dighe, S. Goswami and P. Roy, Phys. Rev. D 76 , 096005 (2007) [arXiv:0704.3735 [hep-ph]].
- [23] J. Bergstrom, M. Malinsky, T. Ohlsson and H. Zhang, Phys. Rev. D 81 , 116006 (2010) [arXiv:1004.4628 [hep-ph]]; J. Bergstrom, T. Ohlsson and H. Zhang, Phys. Lett. B 698 , 297 (2011) [arXiv:1009.2762 [hep-ph]].
- [24] J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl. Phys. B 556 , 3 (1999) [hepph/9904395]; J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl. Phys. B 569 , 82 (2000) [hep-ph/9905381].
- [25] P. H. Chankowski and S. Pokorski, Int. J. Mod. Phys. A 17 , 575 (2002) [hep-ph/0110249].
Shivani Gupta C. S. Kim Pankaj Sharma
Radiative and Seesaw Threshold Corrections to the $S_3$ Symmetric Neutrino Mass Matrix
We systematically analyze the radiative corrections to the $S_3$ symmetric
neutrino mass matrix at high energy scale, say the GUT scale, in the charged
lepton basis. There are significant corrections to the neutrino parameters both
in the Standard Model (SM) and Minimal Supersymmetric Standard Model (MSSM)
with large tan$\beta$, when the renormalization group evolution (RGE) and
seesaw threshold effects are taken into consideration. We find that in the SM
all three mixing angles and atmospheric mass squared difference are
simultaneously obtained in their current 3$\sigma$ ranges at the electroweak
scale. However, the solar mass squared difference is found to be larger than
its allowed 3$\sigma$ range at the low scale in this case. There are
significant contributions to neutrino masses and mixing angles in the MSSM with
large tan$\beta$ from the RGEs even in the absence of seesaw threshold
corrections. However, we find that the mass squared differences and the mixing
angles are obtained in their current 3$\sigma$ ranges at low energy when the
seesaw threshold effects are also taken into account in the MSSM with large
tan$\beta$.