{'fraction_non_alphanumeric': 0.038509348940798674, 'fraction_numerical': 0.02304863521002745, 'mean_word_length': 4.1044282881692, 'pattern_counts': {'":': 0, '<': 10, '<?xml version=': 0, '>': 9, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 2, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Quark lepton complementarity and renormalization group effects

We consider a scenario for the Quark-Lepton Complementarity relations between mixing angles in which the bi-maximal mixing follows from the neutrino mass matrix. According to this scenario in the lowest order the angle \theta_{12} is \sim 1\sigma (1.5 - 2^\circ) above the best fit point coinciding practically with the tri-bimaximal mixing prediction. Realization of this scenario in the context of the seesaw type-I mechanism with leptonic Dirac mass matrices approximately equal to the quark mass matrices is studied. We calculate the renormalization group corrections to \theta_{12} as well as to \theta_{13} in the standard model (SM) and minimal supersymmetric standard model (MSSM). We find that in large part of the parameter space corrections \Delta \theta_{12} are small or negligible. In the MSSM version of the scenario the correction \Delta \theta_{12} is in general positive. Small negative corrections appear in the case of an inverted mass hierarchy and opposite CP parities of \nu_1 and \nu_2 when leading contributions to \theta_{12} running are strongly suppressed. The corrections are negative in the SM version in a large part of the parameter space for values of the relative CP phase of \nu_1 and \nu_2: \phi>\pi/2.


I. INTRODUCTION
Implications of the observed pattern of neutrino mass and mixing (with two large angles) for fundamental physics are still an open question. This pattern has not yet led to a better understanding of the origins of the neutrino mass as well as fermion masses and mixing in general. In contrast, it made the situation more complicated and more intriguing [1]. In this connection, any hint from data and any empirical relation should be taken seriously and analyzed in details.
In fact, one feature has been realized recently that (if not accidental) may lead to a substantially different approach to the underlying physics. Namely, the sums of the mixing angles of quarks and leptons for the 1-2 and 2-3 generations agree with 45 • within 1σ. In other words, the quark and lepton mixings sum up to maximal mixing [2,3,4]: Here θ C is the Cabibbo angle, θ cb ≡ arcsin V cb , and V cb is the element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. According to eqs. (1) which are called the quarklepton complementarity (QLC) relations, the quark and lepton mixings are complementary to the maximal mixing. (A possibility that the lepton mixing responsible for the solar neutrino conversion equals maximal mixing minus θ C was first mentioned in [5], and corrections to the * Electronic address: michael.schmidt@ph.tum.de † Electronic address: smirnov@ictp.trieste.it bimaximal mixing [6] from the CKM type rotations were discussed in [7].) For various reasons it is difficult to expect exact equalities (1). However certain correlations clearly show up: • the 2-3 leptonic mixing is close to maximal because the 2-3 quark mixing, V cb , is very small; • the 1-2 leptonic mixing deviates from maximal substantially because the 1-2 quark mixing (i.e., the Cabibbo angle) is relatively large.
If not accidental coincidence, the QLC relations imply [3,4,8] (i) a kind of quark-lepton symmetry or quark-lepton unification which propagates the information about mixing from the quark sector to the lepton sector.
(ii) existence of some additional structure which produces maximal or bi-maximal mixing.
Even within this context one expects some deviations from exact quark-lepton complementarity due to -broken quark-lepton symmetry, -renormalization group (RG) effects.
There is a number of attempts to reproduce the QLC relations on the basis of already existing ideas about fermion mass matrices [3,8,9,10]. Usually they lead to too small deviations of θ 12 from π/4, and therefore require further corrections or deviations from the bimaximal mixing or from the Cabibbo mixing already in the lowest order. So, in the majority of the models proposed so far, an approximate QLC relation appears as a result of an interplay of different independent factors or as a sum of several independent contributions. In these cases the QLC relation seems to be accidental. There are few attempts to construct a consistent gauge model which explains the QLC relations. The simplest possibility is the SU (2) L × SU (2) R × SU (4) C model that implements the quark-lepton symmetry in the most straightforward way [8,10]. Phenomenology of schemes with QLC relations has been extensively studied [4,9,11,12].
The relation (1) is realized at some high energy scale, M F , of flavor physics and quark-lepton unification. Therefore one should take into account the renormalization group effects on the QLC relations when confronting them with the low energy data. In fact, it was marked in [4] that in the Minimal Supersymmetric Standard Model (MSSM) the corrections are typically positive but negative ∆θ 12 can be obtained from the RG effects in presence of non-zero 1-3 mixing. Also threshold corrections due to some new scale of physics, such as the low scale supersymmetry, can produce a negative shift of θ 12 , thus enhancing its deviation from π/4 [13].
The Cabibbo mixing can be transmitted to the lepton sector in a more complicated way (than via the quarklepton symmetry). In fact, sin θ C may turn out to be a generic parameter of the theory of fermion masses -the "quantum" of flavor physics. Therefore it may appear in various places: mass ratios, mixing angles. One can consider the Cabibbo angle as an expansion parameter for mixing matrices [4,7,14,15].
In this paper we study in details the RG effects in the QLC scenario where the bi-maximal mixing is generated by the neutrino mass matrix. We calculate corrections to the angles both in the Standard model (SM) and MSSM. We analyze the dependence of the corrections on various parameters and obtain bounds on the parameters from consistency condition with QLC. In particular, we find regions where the corrections are negative. The paper is organized as follows. In sec. 2 we formulate the scenario and comment on parameterization dependence of the QLC relations as well as confront the relations with experimental data. In sec. 3 we consider realization of the scenario in the seesaw type I mechanism. The RG effects in MSSM and SM are described in secs. 4 and 5 correspondingly. We consider the RG effects on 1-3 mixing and dependence of the effects on scale of new physics in sec. 6. Conclusions are formulated in sec. 7.

A. A scenario
A general scheme for the QLC relations is that "lepton mixing = bi−maximal mixing − CKM ′′ , (2) where the bi-maximal mixing matrix is [6] Here U m ij is the maximal mixing (π/4) rotation in the ij-plane.
We assume that the bi-maximal mixing is generated by the neutrino mass matrix. That is, the same mechanism which is responsible for the smallness of neutrino mass leads also to the large lepton mixing, and it is the seesaw mechanism [16] that plays the role of additional structure that generates the bi-maximal mixing. Therefore where Γ α ≡ diag(e iα1 , e iα2 , e iα3 ) is the phase matrix that can appear, in general, at diagonalization of the charged lepton or neutrino Dirac mass matrices. Similarity of the Dirac mass matrices in the lepton and quark sectors, related to the quark-lepton symmetry, is the origin of the CKM rotations in the lepton sector. Here, there are two possibilities: (i) In a certain ("symmetry") basis, where the theory of flavor is formulated, the neutrino mass matrix is of the bi-maximal form. So U ν = U bm , and the charged lepton mass matrix is diagonalized by the CKM rotation: The problem here is that the masses of charged leptons and down quarks are different: in particular, m e /m µ = 0.0047, whereas m d /m s = 0.04 − 0.06, and also m µ = m s at the grand unified theory (GUT) scale. Since m l = m d , the equality (5) implies particular structure of the mass matrices in which mixing weakly depends on eigenvalues.
(ii) In the "symmetry" basis both the bimaximal and CKM mixings come from the neutrino mass matrix, and the charged lepton mass matrix is diagonal. That is, the symmetry basis coincides with the flavor basis. In this case the Dirac mass matrix of neutrinos is the origin of the CKM rotation, whereas the Majorana mass matrix of the right-handed (RH) neutrinos is responsible for the bi-maximal mixing. Since the eigenvalues of the Dirac neutrino mass matrix are unknown we can assume an exact equality of the mass matrices as a consequence of the quark-lepton symmetry. The equality (6) propagates the CKM mixing from the quark to the lepton sector precisely. In this case, however, the Gatto-Sartori-Tonin relation between the Cabibbo angle and the ratio of down quark masses [17] turns out to be accidental. Furthermore, one needs to explain why in the symmetry basis both the charged lepton and down quark mass matrices are diagonal simultaneously in spite of difference of eigenvalues. These two cases have different theoretical implications, however the phenomenological consequences and the RG effects are the same.
In the scenario under consideration, the relation (1) is not realized precisely even for zero phases α i since the U CKM 12 rotation matrix should be permuted with U m 23 in (4) to reduce the mixing matrix to the standard parameterization form [4]. From (4) we obtain the following expressions for the leptonic mixing angles: U e2 ≡ cos θ 13 sin θ 12 = sin(π/4 − θ C ) where δ q is the quark CP-violating phase. This expression differs from the one derived in [4] by a factor cos θ 13 as well as by the last term, that turn out to be relevant at the level of accuracy we will consider here. The 1-3 mixing is large in this scenario [4,11,12]: and, hence, the Dirac CP phase δ is close to 180 • . So, for the 1-2 mixing we find from (7) and (8) sin and U e2 is given in (7). Expression for the 2-3 mixing reads U µ3 = cos θ 13 sin θ 12 The RG effect on V CKM is negligible.
The QLC relations are essentially parameterization independent. They can be expressed in terms of physical quantities (compare with [18]). Indeed, the moduli of elements of the mixing matrix, U αi , are physical quantities immediately related to observables and consequently, parameterization independent. In the standard parameterization and therefore | sin θ 12 | = 1 Notice that the presence of 1-3 mixing produces some ambiguity in formulation of the QLC relations. One can write the relations in terms of angles in the standard parameterization or in terms of matrix elements: arcsin (V us Both forms coincide in the limit U e3 → 0.  [19], , Bari [21]. The vertical dashed lines correspond to value (π/4−θC ) and tri-bimaximal mixing. The horizontal bar shows values of QLC predictions without RG corrections.

C. Experimental status
In fig. 1 we show results of determination of θ 12 by different groups. SNO collaboration did analysis of the data in terms of 2ν mixing [19], whereas in [20] and [21] complete 3ν analyses have been performed. Furthermore, in [21] a non-zero best fit value of 1-3 mixing has been obtained. Results of different analyses are in a very good agreement: Notice that the determination of θ 12 follows mainly from analysis of the solar neutrino data. In this analysis θ 12 and θ 13 correlate. In particular, the CC/NC ratio that gives the most important restriction on mixing is determined by P ∼ cos 4 θ 13 sin 2 θ 12 . The best fit values (14) are along with the trajectory P = constant.

D. QLC and tri-bimaximal mixing
The QLC prediction (16) is practically indistinguishable from the tri-bimaximal mixing [23] prediction sin 2 θ 12 = 1/3. So, it turns out that almost the same values of 1-2 mixing are obtained from two different and independent combinations of matrices There are two possible interpretations of this fact: The coincidence is accidental, which means that one of the two approaches (QLC or tri-bimaximal mixing) does not correspond to reality. To some extend that can be tested by measuring the 1-3 mixing. In the QLC-scenario one obtains (17), whereas the tri-bimaximal mixing implies sin 2 θ 13 = 0 unless some corrections are introduced.
The coincidence is not accidental, and therefore it implies a non-trivial expression for the Cabibbo angle. Indeed, from the equality sin θ QLC 12 = sin θ tbm 12 we obtain

III. SEESAW AND QLC
The seesaw mechanism [16] that provides a natural explanation of smallness of neutrino mass can also be the origin of the difference of the quark and lepton mixings, and in particular, the origin of bimaximal mixing. In the context of seesaw type-I this implies a particular structure of the RH neutrino mass matrix.
We assume that the type-I seesaw gives the dominant contribution to neutrino masses since it can provide the closest relation between the quark and lepton mass matrices, as required by the QLC. The relevant terms of the Lagrangian are where L ≡ (ν, l) is the leptonic doublet, N ≡ (ν R ) c , H d and H u are two different Higgs doublets in the MSSM and H u = iτ 2 H d in SM; Y e and Y ν are the charged lepton and neutrino Yukawa coupling matrices. We consider M R as a bare mass matrix of the RH neutrinos formed already at the GUT or even higher scale. It can be generated by some new interactions at higher scales. Decoupling of N leads to the low energy effective D = 5 operator After the electroweak symmetry breaking this operator generates the mass term for light neutrinos, ν T m ν ν with Let us consider the basis where the neutrino Dirac mass matrix is diagonal: Then the light neutrino mass matrix equals and m diag D ≡ Y diag ν H u . According to our assumption, the matrix (29) should generate the bimaximal rotation: where in general, Here is the phase matrix, is the diagonal matrix of the light neutrinos, and with ϕ i being the Majorana phases of light neutrinos.
According to our assumption, the CKM rotation follows from diagonalization of the charged lepton matrix and we will parameterize it as Here the diagonal matrix of the phase factors on the RH side has been absorbed in the charged lepton field redefinition; V CKM is the CKM matrix in the standard parameterization, θ q and δ q are the quark (CKM) mixing angles, and Thus, in general, there are three matrices of phases, Γ δ , Γ ϕ and Γ φ , relevant for relations between the mixing angles. Finally, from (36) and (31) we obtain and therefore in (4) The neutrino mass matrix in the flavor basis equals From (29) and (31) we find an expression for the RH neutrino mass matrix: Omitting the phase factor Γ(δ i ) (that can be absorbed in the definition of M R ) and including the CP phases Explicitly where (with phases ϕ i included). We can parameterize m diag with m t being the mass of top quark and ǫ ′ ≈ ǫ ∼ 3·10 −3 .
Using smallness of ǫ's it is easy to estimate the mass eigenvalues: Furthermore, the 1-2 and 2-3 mixing angles are of the order ǫ, whereas 1-3 mixing is of the order ǫ 2 .
In the case of normal mass hierarchy, in agreement with results of [24]. Notice a permutation character of these expressions: the masses of RH neutrinos 1, 2, 3 are determined by light masses 2, 3, 1. With (47) So, masses have a "quadratic" hierarchy.
In the case of inverted mass hierarchy, m 3 ≪ m 1 ≈ m 2 ≡ m A , and the same CP phases of ν 1 and ν 2 we obtain from (45) where m A ≡ |∆m 2 31 |. This leads again to a strong mass hierarchy. Notice that now the mass of the lightest RH neutrino is determined by the atmospheric mass scale. Thus, apart from special regions in the parameter space that correspond to level crossings (see sect. 5) the QLC implies generically a very strong ("quadratic") mass hierarchy of the RH neutrinos and very small mixing: Θ ij ∼ ǫ. As we will see, this determines substantially the size of the RG effects.
Let us introduce the unitary matrix, U R , which diagonalizes the right-handed neutrino mass matrix and the mixing matrix can be parameterized as where Θ ij and ∆ are the angles and CP-phase of the RH neutrino mixing matrix.
In what follows we will not elaborate further on the origin of particular structures of M R (42), just noticing that it can be related to the double (cascade) seesaw mechanism [25] with the "screening" of Dirac structure [26,27].

A. General consideration
The quark-lepton symmetry implied by the QLC relation means that physics responsible for this relations should be realized at some scale M F which is at the quark-lepton unification scale, M GUT , or even higher scales. An alternative possibility would be the quarklepton relation due to the Pati-Salam symmetry [28] broken below the GUT scale. Consequently, there are three different regions of RG running: (i) below the seesaw scales, µ < M 1 , where M 1 is the lightest RH neutrino mass. In this region all three neutrinos decouple and the D=5 operator (26) is formed; (ii) between the seesaw scales, The RG equation for the neutrino mass matrix is given by [29,30,31] 16π 2ṁ whereṁ ν ≡ µ dm ν /dµ, µ is the renormalization scale, and κm ν includes the gauge interaction terms that can influence the flavor structure in the SM case (see below); C e = −3/2, C ν = 1/2 in the SM and C e = C ν = 1/2 in the MSSM. From the evolution equation (51) for the mass matrix we can obtain the equations [42] for observables (masses and mixings) [32,33,34,35]. Above and below the seesaw scales the gauge interactions produce a flavor universal effect and all contributions (from all RH neutrinos) to the neutrino mass matrix have the same renormalization group equation.
In the limit of vanishing 1-3 mixing the evolution of θ 12 is described approximately by 32π 2θ 12 = Q 12 sin 2θ 12 P 11 − c 2 23 P 22 − s 2 23 P 33 + s 23 Re P 23 + +2 cos 2θ 12 (c 23 Re P 21 − s 23 Re P 31 )] + 4S 12 (c 23 Im P 21 − s 23 Im P 31 ) , (53) where s 23 ≡ sin θ 23 , c 23 ≡ cos θ 23 , etc., and Above the seesaw scale one needs to consider renormalization of couplings of the full Lagrangian (25). The evolution of the effective operator which gives masses to neutrinos after the electroweak symmetry breaking is determined by evolution of the neutrino Yukawa couplings Y ν and the mass terms of right-handed neutrinos. Below the seesaw scales, running is dominated by P 33 in the flavor basis which results in an increase of θ 12 in MSSM and a slight decrease in the SM due to different signs of C e : 32π 2θ 12 ≈ −Q 12 sin 2θ 12 s 2 23 P 33 .
Above the seesaw scales, the leading contribution is again given by P 33 , and the next-to-leading contribution comes from P 32 . This yields an increase of θ 12 when running to low scales both in the MSSM and in SM.
Explicitly the corresponding evolution equation can be written as and the phases φ i are determined in Eq. (36,37). Effect of running between the seesaw scales (about 10orders of magnitude in µ) is more complicated. In this range the Yukawa coupling matrix has two terms (contributions): (1) D=5 effective operators for the light neutrinos formed after decoupling of one or two RH neutrinos, Dirac type couplings and mass terms for undecoupled RH neutrinos given by the Lagrangian (25).
These terms are renormalized differently. In particular, for the terms of second type the neutrino Yukawa couplings are important.
The difference, however, cancels (between the seesaw scales) in the case of MSSM in which only the wave function renormalization takes place due to the non-renormalization theorem. In contrast, in the SM due to vertex corrections to the D=5 operators the difference does not cancel, and, as we will see, produces a significant effect. So, in the MSSM, the RG equations are the same for both contributions and eq. (53) is valid. In the SM they are not equal and eq. (53) can not be applied.
After the heaviest right-handed neutrino is integrated out, the right-handed neutrino mixing at the threshold influences running of θ 12 . In the second order of sin θ C , the expression forθ 12 reads: where Θ ij are the right-handed neutrino mixing angles at the scale at which the heaviest neutrino is integrated out. The unitary rotation of the right-handed neutrino fields is done at the threshold of the heaviest right-handed neutrino, and the exact definitions of the angles are given in Eq. (49, 50).

B. RG evolution and scales of flavor physics
We have performed running from the M F scale down to the electroweak scale and calculated ∆θ 12 ≡ θ 12 (M Z ) − θ 12 (M F ). For that we solved numerically a complete set of the RG equations including sub-leading effects due to the non-zero 1-3 mixing. In most of our calculations we take for definiteness M F = M GUT = 2 · 10 16 GeV. We consider separately dependence of the results on M F . Notice that the renormalization of θ 12 in the bimaximal scheme has been studied in [36].
The following free parameters determine the RG effects substantially: the absolute scale of light masses m 1 , the CP (Majorana) phases of light neutrinos, ϕ i , and the phases α i . We studied dependence of the RG effects on these parameter. For each set of the parameters we have calculated the RH mass matrix and running effects. The angles are fixed by the QLC relation at M F , and the mass squared differences are adjusted to lie in the experimentally allowed region at the electroweak scale. For the neutrino Yukawa couplings we take y 1 : y 2 : y 3 = ǫ 2 : ǫ : 1, (ǫ = ǫ ′ ) and ǫ = 3 · 10 −3 .
We consider RG evolution in the MSSM with a unique SUSY threshold 1 TeV. The RG effects depend on the absolute mass scale, m 1 , tan β, and the relative phase between the first and second mass eigenstates ϕ ≡ ϕ 2 − ϕ 1 . Dependence on other parameters (e.g., other phases) is rather weak. Still we will use explicitly the phase ϕ 2 keeping everywhere ϕ 1 = 0.
In what follows we will describe results of our numerical calculations. We give an interpretation of the results using approximate formulas presented in Secs. III and IV.1.

C. RG effect in MSSM with normal mass hierarchy
In fig. 2 we show some examples of the scale dependence of θ 12 for various values of parameters. With increase of m 1 two factors enhance the RG effects: (1) the largest mass M 3 decreases according to (46). Correspondingly, the region above the seesaw scale, M 3 − M GUT increases; (2) corrections to the mass matrix elements are proportional to values of elements: ∆m αβ ∝ m αβ and since with increase of m 1 the masses, m αβ , generically increase, the corrections increase correspondingly.
For relatively small tan β ∼ (3 − 10), the dominant contribution follows from region above the seesaw scales due to large (Y ν ) 33 . Evolution below M 3 is mainly due to the Yukawa couplings Y e which are relatively small. The effect increases fast with m 1 : Notice that M 3 ∝ 1/m 1 . Therefore for m 1 ∼ 10 −3 eV the running of θ 12 is mainly related to increase of region above the seesaw scale. For m 1 > 10 −2 eV the spectrum of light neutrinos becomes degenerate and ∆θ 12 ∝ Q 12 ∝ m 2 1 ( fig. 2(a)). For large tan β and small m 1 the dominant contribution to ∆θ 12 comes from the region below M 3 where ∆θ 12 ∝ tan 2 β (see fig. 2

(b)).
A combined dependence of corrections, ∆θ 12 , on m 1 and tan β is presented in fig. 3 where we show contours of constant ∆θ 12 in the (m 1 −tan β) plane. The change of behavior of contours at m 1 = 8·10 −4 eV is a consequence of our boundary condition: At m 1 < 8 · 10 −4 eV we have M 3 > M GUT , and therefore the region above seesaw scale disappears.
In fig. 4 we show the correction ∆θ 12 as functions of m 1 for different values of ϕ 2 . The dependence of ∆θ 12 on ϕ 2 , given essentially by the factor Q 12 , is weak for the hierarchical spectrum, m 1 ≪ 8 · 10 −3 eV, and very strong for the degenerate spectrum: ∆θ 12 ∝ (1 + cos ϕ). The corrections are strongly suppressed for the opposite CP-parities, ϕ 2 = 180 • , ( fig. 4) that agrees with the results of previous studies of corrections in the degenerate case [37,38,39,40].
Condition that the QLC prediction for θ 12 is within 1σ of the best fit experimental value requires ∆θ 12 < 0.5 • − 1 • . This, in turn, leads to bounds on parameters of neutrino spectrum and tan β. In particular, according to fig. 4 the degenerate neutrino spectrum is excluded for the same CP parities (ϕ 2 = 0). In the case of large tan β it requires strongly hierarchical spectrum: m 1 < 10 −3 eV that eliminates the running region above seesaw scale. However, a degenerate spectrum is allowed for ϕ ∼ 180 • .

D. MSSM and inverted mass hierarchy
In the case of inverted mass hierarchy, the states ν 1 and ν 2 associated to 1-2 mixing are strongly degenerate. Therefore, the RG effects are similar to those in the normal hierarchy case for m 1 = m A ∼ 5 · 10 −2 eV. So, the corrections, ∆θ 12 , are enhanced by the factor where subscripts NH and IH stand for normal and inverted mass hierarchy. This factor equals for the strong normal hierarchy and normal ordering (m 1 ≈ m 2 ) correspondingly.
In fig. 5 we show examples of running of θ 12 for different values of masses and phases. Dependences of θ 12 are well described by Q 12 , as in the case of normal mass hierarchy. Notice that now the heaviest RH neutrino mass is determined by m 3 , and two others by m A . With increase of m 3 (now the lightest neutrino mass) ( fig. 5(a)) the range above the seesaw scales, where the evolution of θ 12 is most strong, increases. The change of θ 12 below M 3 is slower being of the same size for different values of m 3 (until m 3 ≪ m A ). In this range the evolution is essentially due to Y e couplings, so that ∆θ 12 ∝ tan 2 β ( fig.  5(b)). The correction can be strongly suppressed for the opposite CP-parities of ν 1 and ν 2 : ∆θ 12 ∝ (1 + cos ϕ).
As in the case of normal hierarchy (ordering), in a large part of the parameter space the correction is positive, ∆θ 12 > 0, due to dominant effect of P 33 . For ϕ 2 = 0, consistency of the QLC prediction with data taken as ∆θ 12 < 2 • , implies tan β < 10 and m 3 < 8 · 10 −4  eV. For ϕ 2 ∼ π corrections can be strongly suppressed, so that a larger region of the parameter space becomes allowed. The corrections become negative for ϕ 2 = π (see figs. 6, 7) when the leading RG effects are strongly suppressed and the running is mainly due to sub-leading effect related to non-zero 1-3 mixing. This possibility has been mentioned in [4]. The sign of contribution due to non-zero θ 13 to the RG running of θ 12 due to non-zero θ 13 depends on the parameter (masses, phases) region.

V. RG EFFECTS IN THE STANDARD MODEL
In the standard model the evolution of θ 12 is more complicated. As we have already mentioned, apart from the Yukawa coupling contributions described by eq. (51) there are additional vertex diagrams that cancel in the SUSY case [41]. Furthermore, the vertex diagrams with the gauge bosons become important: their contribution to running between the seesaw scales influences the flavor structure of mass matrix and therefore changes θ 12 [33,34,35]. The point is that individual RH neutrinos, N i , have the flavor dependent couplings with the left handed components of neutrinos. Therefore the gauge boson corrections to the corresponding couplings will influence the flavor structure. Above the seesaw scales (where all RH neutrinos are operative) and below the seesaw scales (where all RH neutrinos decouple), flavor universality of the gauge interaction corrections is restored. There is no simple analytic formula for the θ 12 renormal- ization in the SM.
In fig. 8 we show examples of the θ 12 scale evolution. Above the seesaw scales the running is due to the Yukawa interactions, Y ν , and the effect is well described by the analytic results (53). Below the seesaw scales, µ < M 1 , the evolution is negligible: it is related to Y e couplings that are small in the SM. The main effect is collected between the seesaw scales. As we mentioned above, it is mainly due to the gauge vertex corrections: After decoupling of N 3 the D=5 operator is formed. The vertex diagram corrections to this operator due to the gauge (and also Yukawa) interactions produce running. Notice that for other RH neutrinos that do not decouple, the corresponding couplings produce box diagrams with propagators of the RH neutrinos. Those diagrams are finite and do not lead to logarithmic corrections. The gauge interaction effect dominates since N 3 with the largest Yukawa coupling is decoupled and Y e are small. The corrections increase with m 1 . The most interesting dependence of ∆θ 12 is the one on the CP-violation phase ϕ 2 ( fig. 8(b)). The corrections are positive, ∆θ 12 > 0, for ϕ 2 = 0. They are strongly suppressed for ϕ 2 ∼ π/2, in contrast to the SUSY case where suppression is realized for ϕ ∼ π. The corrections are negative for ϕ 2 > π/2. The angle of zero corrections, ϕ 2 (0), depends on m 1 and in general deviates from π/2. The deviation is due to the Yukawa interaction effects that produce the positive shift for strong Yukawa coupling hierarchy as we discussed before. The shift occurs both above and between the seesaw scales (see fig. 8(b)).
In fig. 9 we show contours of constant corrections in the m 1 − ϕ 2 plane, and in fig. 10 -an explicit dependence of ∆θ 12 on m 1 for different values of ϕ 2 . The line ∆θ 12 = 0, is close to ϕ 2 = π/2, 3π/2 for m 1 → 0, and it approaches π with increase of m 1 when spectrum becomes strongly degenerate. The pattern is nearly symmetric with respect to ϕ = π for small m 1 , the asymmetry appears for m 1 > 3 · 10 −3 eV.

A. Renormalization of 13 mixing
In the scenario discussed in this paper, the 1-3 mixing is non-zero and relatively large at the boundary (17). Notice that θ 13 (i) interferes with 1-2 mixing in the QLC relation as we discussed before; (ii) produces sub-leading effects in renormalization of θ 12 , (iii) can provide further bounds on the considered scenario if RG corrections are positive and large.
The running is suppressed by small mass ratio. Therefore only a small RG effect on 1-3 mixing appears for the hierarchical (normal as well as inverted) case. For instance, we find that for the parameter sets used in figure  2 (MSSM), the correction ∆θ 13 is always smaller than 0.2 • . In the SM, it is smaller than 0.3 • . For the degenerate spectrum, there can be a larger effect which strongly depends on the CP-phases. From   (68) Notice that for zero CP phases the cancellation occurs again. In the MSSM for m 1 = 0.03 eV and tan β = 50, we find ∆θ 13 ∼ 0.5 • . In contrast, for δ = ϕ 1 = π and ϕ 2 = 0 the two terms in (68) sum up and we obtain positive running: 64π 2θ 13 ≈ 4 sin 2θ 12 sin 2θ 23 . Consequently θ 13 becomes smaller at low energies.

B. Level crossing points
As we have established in sect. 3 the spectrum of the right-handed Majorana neutrinos is generically hierarchical. However, there are the level crossing points, where two of the RH neutrino masses become equal [24]. The case of degeneracy of two lighter RH neutrino states, M 1 ≈ M 2 , is of special interest from the point of view of generation of the baryon asymmetry in the Universe. In this case the resonance leptogenesis becomes possible which produces large enough asymmetry in spite of smallness of the masses and consequently, large wash out effect.
From (45) we find (Here the Majorana phases are included in m i ). It is easy to see that due to smallness of ǫ the condition M 1 ≈ M 2 can be satisfied only in the case of strong mass degeneracy |m 1 | ≈ |m 2 | ≈ m 0 when Then from the condition M 1 ≈ M 2 we find In this special case the mass is enhanced by factor 2m 0 / ∆m 2 21 ∼ 20 and the third mass is much smaller than in the hierarchical case: that is, smaller by factor m N H 1 /m 3 < 10 −3 . The level crossing condition (70) implies the opposite CP-violating phases; it coincides with the condition of strong suppression of the RG effects. It also implies smallness of the 11-element of m bm matrix. The condition for level crossing differs from that in [24] since here we require the neutrino Dirac matrix to be diagonal in the basis where the mass matrix of light neutrinos has exactly bimaximal form. If instead we use a generic matrix with non-maximal 1-2 mixing the level crossing condition can be realized for the hierarchical spectrum [24].
We find that in this crossing point the solar mass squared difference becomes large even if it is very small at the boundary. So, the 1-2 split has the radiative origin. The 1-3 split decreases by factor ∼ 2.

C. Evolution above the GUT scale
For M F > M GUT one should perform running also above the GUT scale. Restoration of the GUT sym- metry and unification of the gauge couplings does not prevent from different running of the Yukawa couplings, and therefore, from change of mixing angles. Renormalization of mixing angles would stop after possible unification of the Yukawa couplings which can be related, e.g., to restoration at M F a non-Abelian flavor symmetry. An alternative is the boundary at the string or Planck scale, where the Yukawa couplings are formed and their properties are determined immediately by some symmetry or/and string selection rules.
For illustration we performed the running in the MSSM up to the Planck scale (ignoring possible GUT effects). In fig. 12 we show the dependence of ∆θ 12 on m 1 for the same (QLC) initial conditions at the Planck scale: M F = M P l = 1.2 · 10 19 GeV. The RG effect becomes much larger. In particular the contribution from the region above the seesaw scale due to large Yukawa coupling Y ν increases substantially. It is enhanced in comparison with the case of running up to M GUT by the factor that can be as large as 3 -5 in some cases. Still for ϕ 2 = π or for small m 1 the RG effects are suppressed and can be consistent with the QLC relations.
Similar RG effects are expected in the SU(5) model with the RH singlet neutrinos. In fact, no new diagrams with large Y ν appear. Effect of the charged lepton couplings Y e is enhanced by factor 4 above M GUT due to the loop diagrams with down quarks (squarks) and H 1/3 charged Higgs bosons (Higgsinos).
The flavor-diagonal parts of the RG equations do influence the angles only indirectly through the change of the mass eigenvalues. Thus, the main effect of these interactions is due to the evolution of ∆m 2 12 .

VII. CONCLUSION
Experimental results on the 1-2 and 2-3 mixing in the quark and lepton sectors show certain correlations that can be interpreted as the quark-lepton complementarity. We considered the QLC scenario in which the bimaximal mixing follows from diagonalization of the neutrino mass matrix. In the lowest order of the perturbation theory, the value of angle θ 12 predicted in this scenario is about ∼ 1σ larger than the best fit experimental value. It coincides practically with the value given by the tribimaximal mixing. We commented on implications of this equality as well as on perspectives of future tests of the QLC relations.
In this paper, we assumed that the QLC relations are not accidental coincidences, but consequence of the quark-lepton symmetry and additional structure in the theory that produces the bimaximal mixing. Here the accidental coincidences would mean that values of mixing angles are the result of interplay of two or more independent contributions. In this connection, we proposed a realization of the QLC scenario that provides the closest relation between quarks and leptons. It is based on -the seesaw type-I mechanism that generates the bimaximal mixing due to specific structure of the RH neutrino mass matrix; -approximate equalities of the Dirac mass matrices: m u ≈ m D , m l ≈ m d that follow from the approximate quark-lepton symmetry or unification. A certain small violation of equalities of these matrices produces difference of mass hierarchies but does not affect substantially the mixing.
The only additional (and, in fact, unavoidable) factor that can affect the QLC relations is RG corrections.
One of the consequences of the proposed scenario is a very strong hierarchy of the RH neutrino masses, apart from several particular level crossing points. The latter are realized for strongly degenerate spectrum of light neutrinos, and particular values of the CP-violating phases. This determines substantially the RG effects.
We performed a systematic study of the RG effects in the SM and MSSM. We find that in the MSSM, the RG corrections to θ 12 are generically positive due to a dominant effect of the Yukawa coupling Y 33 . So, these corrections worsen agreement of the predicted θ 12 with data.
In the MSSM small negative corrections, |∆θ 12 | < 0.5 • , can appear for the opposite CP parities of ν 1 and ν 2 and inverted mass hierarchy, in which case the main terms in the RG equations are strongly suppressed and running is due to the sub-leading effects related to the non-zero 1-3 mixing. The correction increases with m 1 and strongly depends on the relative Majorana phase. For ϕ 2 = 0 the consistency of the QLC prediction for θ 12 with data implies strong mass hierarchy of light neutrinos and small tan β. For ϕ 2 = π corrections are suppressed and even the degenerate spectrum becomes allowed. For the inverted mass hierarchy the corrections are generically enhanced by larger values of masses of ν 1 and ν 2 .
The situation is qualitatively different in the SM. Here important contributions follow from the vertex corrections to the D=5 operator in the range between the seesaw scales. The Yukawa couplings (especially for small m 1 ) give sub-leading contribution. The RG corrections are negative in the interval ϕ 2 = π/2 − 3π/2 for small m 1 and the range of negative corrections becomes narrower, ϕ 2 = (0.9 − 1.2)π, for the degenerate neutrinos ν 1 and ν 2 .
Corrections depend substantially on the boundary scale M F . The value ∆θ 12 can be enhanced by factor 2-5 if M F increases from M GUT to M P l .
For the hierarchical mass spectrum renormalization of the 1-3 mixings is, in general, small: ∆θ 13 ∼ 0.2 • − 0.3 • . The correction can be large, ∆θ 13 ∼ θ 13 , for the degenerate spectrum. The sign of correction depends on values of CP-violating phases.
In conclusion, in a large part of the parameter space especially for the strong mass hierarchy and opposite CP phases of ν 1 and ν 2 , the RG corrections to the QLC relation are small. The corrections are positive in the MSSM apart from small region of parameter space that corresponds to the degenerate spectrum of light neutrinos and their opposite CP parities. The corrections are negative in the SM for ϕ 2 > π/2. In the considered QLC scenario the RG corrections allow one to reproduce the best fit experimental value of θ 12 exactly.