Trimaximal mixing with predicted \theta_13 from a new type of constrained sequential dominance

Following the recent T2K indication of a sizeable reactor angle, we present a class of models which fix \theta_{13} while preserving trimaximal solar mixing. The models are based on a new type of constrained sequential dominance involving new vacuum alignments, along the (1,2,0)^T or (1,0,2)^T directions in flavour space. We show that such alignments are easily achieved using orthogonality, and may replace the role of the subdominant flavon alignment (1,1,1)^T in constrained sequential dominance. In such models, with a normal hierarchical spectrum, the reactor angle is related to a ratio of neutrino masses by \theta_{13} = \frac{\sqrt{2}}{3}\frac{m^\nu_2}{m^\nu_3}, leading to \theta_{13} \sim 5^\circ - 6^\circ, while the atmospheric angle is given by the sum rule \theta_{23} \approx 45^\circ + \sqrt{2} \theta_{13} \cos \delta. We find that leptogenesis is unsuppressed due to the violation of form dominance and that the CP violating phase responsible for leptogenesis is precisely equal to the Dirac CP phase \delta, providing a direct link between leptogenesis and neutrino mixing in this class of models.


Introduction
Recently T2K have published evidence for a large non-zero reactor angle [1] which, when combined with data from MINOS and other experiments in a global fit, yields [2,3] 6 • θ 13 with a statistical significance of a non-zero reactor angle of about 3σ. If confirmed this would rule out the hypothesis of exact tri-bimaximal (TB) mixing [4], and a flurry of alternative proposals have recently been put forward [5]. For example, an attractive scheme based on trimaximal (TM) mixing remains viable [6]. TM mixing is defined to maintain the second column of the TB mixing matrix and hence preserves the solar mixing angle prediction sin θ 12 ≈ 1/ √ 3. However there is another variation of TM mixing which also preserves this good solar mixing angle prediction by maintaining the first column of the TB matrix, namely TM 1 mixing [7]: where 1 √ 3 sin ϑ = sin θ 13 , P ′ is a diagonal phase matrix required to put U TM 1 into the PDG convention [8], and P = diag(1, e i α 2 2 , e i α 3 2 ) contains the usual Majorana phases. In particular TM 1 mixing approximately preserves the successful TB mixing for the solar mixing angle θ 12 ≈ 35 • as the correction due to a non-zero but relatively small reactor angle is of second order. Although TM 1 mixing reduces to TB mixing in the limit that ϑ → 0, it is worth emphasising that in general TM 1 mixing involves an undetermined reactor angle θ 13 which could in principle be large or even maximal (e.g. 45 • ). The observed smallness of the reactor angle θ 13 compared to the atmospheric angle θ 23 ≈ 45 • and the solar angle θ 12 ≈ 34 • [2] is therefore not explained by the TM 1 hypothesis alone. Clearly the relative smallness of the reactor angle can only be explained with additional model dependent input. Although there are models of TM mixing which can account for the smallness of the reactor angle [9] so far there is no model in the literature for TM 1 mixing, let alone one which fixes the reactor angle.
In this paper we propose a model of TM 1 mixing where the magnitude of the reactor angle is fixed. The model we discuss is actually representative of a general strategy for obtaining TM 1 mixing using sequential dominance (SD) [10] and vacuum alignment. The strategy of combining SD with vacuum alignment is familiar from the constrained sequential dominance (CSD) approach to TB mixing [11] where a neutrino mass hierarchy is assumed and the dominant and subdominant flavons responsible for the atmospheric and solar neutrino masses are aligned in the directions of the third and second columns of the TB mixing matrix, namely φ ν 1 ∝ (0, 1, −1) T and φ ν 2 ∝ (1, 1, 1) T . The new idea here is to maintain the usual vacuum alignment for the dominant flavon, φ ν 1 ∝ (0, 1, −1) T as in CSD, but to replace the effect of the subdominant flavon vacuum alignment by a different one, namely either φ 120 ∝ (1, 2, 0) T or φ 102 ∝ (1, 0, 2) T , where such alignments may be naturally achieved from the standard ones using orthogonality arguments. We shall refer to this new approach as CSD2. We shall show that CSD2 leads to TM 1 mixing and a reactor angle which, at leading order, is predicted to be proportional to the ratio of the solar to the atmospheric neutrino masses, It is interesting to compare the predictions of CSD2 to another alternative to CSD that has been proposed to account for a reactor angle called partially constrained sequential dominance (PCSD) [12]. PCSD involves a vacuum misalignment of the dominant flavon alignment to (ε, 1, −1) T , with a subdominant flavon alignment (1, 1, 1) T , leading to tri-bimaximal-reactor (TBR) mixing [12] in which only the reactor angle is switched on, while the atmospheric and solar angles retain their TB values. However, in the case of PCSD, the value of the reactor angle is not predicted whereas CSD2 leads to the above relation.
The layout of the rest of the paper is as follows. In Section 2 we describe CSD2 including a discussion of vacuum alignment and an example of a model based on CSD2. In Section 3 we discuss the phenomenology of CSD2, first showing that it reproduces TM 1 mixing exactly, then comparing the second order analytic results to a numerical treatment. In Section 4 we show that leptogenesis is unsuppressed and moreover CSD2 leads to a link between the CP phase for leptogenesis and the Dirac CP phase δ. We conclude in Section 5.  [10]. The Yukawa couplings in A determine the atmospheric angle, while those in B determine the solar angle, with the reactor angle dependent on both [10]. For a strong hierarchy one can ignore the Yukawa couplings in C and the effect of the third right-handed neutrino.
TB mixing naturally emerges from CSD [11] where the neutrino Yukawa matrix Y ν and the right-handed neutrino mass matrix M R take the constrained form assuming m ν 1 = 0. In models with non-Abelian family symmetries this structure can be explained by two flavons pointing in the directions in flavour space defined by the columns of Y ν :

CSD2 and vacuum alignment
As our starting point, we consider the vacuum alignment sector of any SD flavour model where the flavons, which furnish triplet representations under a family symmetry G F , are typically aligned in the three orthogonal directions Applying these flavons to build up the neutrino Yukawa matrix, one would obtain TB neutrino mixing and a form diagonalisable neutrino mass matrix. In many models three additional orthogonal flavons are present, The third flavon φ e 3 is generically introduced to generate hierarchical quark and charged lepton mass matrices, and typically governs the masses of the third generation of charged particles.
In this paper we augment this set of CSD flavons by another flavon pointing in the direction These alignments can be obtained by requiring the orthogonality conditions φ 120 · φ ν 3 = φ 120 · φ e 3 = 0 or alternatively φ 102 · φ ν 3 = φ 102 · φ e 2 = 0, where the "·" denotes the usual SO(3) inner product. Implicitly we assume real triplet representations, but it is straightforward to extend this mechanism to complex representations. The orthogonality conditions can be realised easily with "Lagrange multiplier" superfields D 1 and D 2 , which are singlets under the family symmetry G F . They give the following terms in the superpotential The F -term conditions |F D 1 | = |F D 2 | = 0 are equivalent to the orthogonality conditions and therefore yield the desired alignments. The shaping symmetries of a specific model select the flavons which couple to the RH neutrinos and thus determine the structure of the neutrino Yukawa matrix. CSD2 corresponds to the case where the dominant flavon is φ ν 1 as in CSD, and the subdominant flavon is taken to be either φ 120 or φ 102 . In the following subsection, we sketch a model in which this is achieved.

An example of a model with CSD2
In the following we briefly outline how to implement the discussed alignment into a concrete model of lepton masses and mixing angles. As the model is mainly for the purpose of illustration, we do not discuss the quark sector and the charged lepton mass matrix is diagonal by construction.
We start with a short discussion on how to achieve the alignments in Eqs. (2.4) and (2.5). This is rather standard and we follow closely the alignment as discussed in [13]. Although we do not implement spontaneous CP violation we stick to Z 4 symmetries to make the analogy more obvious.
The vacuum alignment of φ e 1,2,3 is generated by the renormalisable superpotential where "⋆" denotes the symmetric cross product in A 4 , which is defined as (x ⋆ y) i = s ijk x j y k , with s ijk = +1 on all permutations {i, j, k} ∈ {1, 2, 3} and zero otherwise, see also [16]. Here and in the rest of this section we drop coupling constants in the superpotential to increase the  4 . We also assume a U (1) R symmetry under which the matter fields carry a charge of one, the flavons are neutral and the driving fields are doubly charged.
For the alignment of the φ ν flavons we employ the superpotential where we use a similar notation as before and the charges under the A 4 and Z 4 symmetries are distributed in the same way. Note that we have introduced an additional auxiliary singlet flavon field ξ with a non-zero vev to align φ ν 2 . Having obtained the alignment φ ν 2 from the A ν 2 driving field, the remaining neutrino flavons φ ν 1 and φ ν 3 are aligned by the driving fields O νe 11 and O ν ij . Eventually, the new flavon alignments are achieved by adding the two singlet driving fields D 1 and D 2 of Eq. (2.7) as discussed previously. D 1 is charged under Z for φ 102/120 . With the symmetries and the field content as given in Tab. 1 we end up with the following Yukawa superpotential which give the Yukawa couplings after the flavons develop their vevs. Λ is a generic messenger mass scale. The charged lepton Yukawa matrix is diagonal due to the alignment of the φ e flavons: Y e = diag(y e , y µ , y τ ) , (2.11) and the neutrino Yukawa couplings depending on the choice of the subdominant flavon, either φ 120 or φ 102 . The parameters a and b can be determined from the parameters in the superpotential. Later on we will see, that a relative phase difference arg(a/b) = 45 • , 135 • , 225 • or 315 • , which translates into a Dirac CP phase δ = 90 • or 270 • , is preferred by experimental data and that this would also maximise the generated baryon asymmetry. Such a phase difference can be easily obtained in the context of spontaneous CP violation from discrete symmetries as discussed in [13], which could be applied here straightforwardly. Due to the Z 4 symmetries the RH neutrinos have no mass terms at the renormalisable level, but they become massive after the flavons develop their vevs due to the following terms in the superpotential From the symmetries alone also terms like φ ν 1 · φ 120/102 N 1 N 2 would be allowed, but we assume, that the messenger fields mediating such operators are absent. Under this common assumption the RH neutrino mass matrix is diagonal (2.14) 3 The phenomenology of CSD2

Predictive trimaximal mixing from CSD2
With the charged lepton mass matrix being diagonal, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing originates solely from the neutrino sector. As discussed in the previous section we introduce two RH neutrinos N i (i = 1, 2) which entails one massless light neutrino. The RH neutrino mass matrix M R is assumed to be diagonal and each N i couples to its own flavon. Adopting φ ν 1 for the dominant and φ 120 for the subdominant term, the resulting neutrino Yukawa matrix is Y (120) ν , see Eq. (2.12). 2 Then the (type-I) seesaw formula leads to a simple effective 2 We comment on the case where the subdominant flavon is taken to be φ102 below Eqs. (3.6-3.13). Note that this resembles the case with two-texture zeros in the neutrino Yukawa matrix, whose phenomenology was extensively discussed in [14,15]. Here we go beyond those papers by giving an explicit vacuum alignment mechanism for such a texture, where all mixing angles and phases depend on a single complex parameter.
light neutrino mass matrix, given by can in general be complex. We assume |m b | ≪ |m a |, which can originate from a hierarchy in M A and M B , in the parameters a and b, or a combination of both. This is nothing but sequential dominance.
The scale of family symmetry breaking, the messenger scale Λ and the scale of the righthanded neutrino masses is undetermined and can be chosen appropriately. This is due to the fact, that from which this freedom of choice is obvious. Clearly, the unitary matrix that diagonalises M ν depends on only one complex parameter An overall Majorana phase related to the parameter m a can be set to zero without loss of generality. With |m b | ≪ |m a | we can use ǫ as an expansion parameter for the resulting neutrino masses and mixing angles. We have computed our results to second order in ǫ. General formulas can be found in [10]. Before determining the PMNS matrix explicitly, it is worth to show why the neutrino mass matrix M ν in Eq. (3.1) implies the TM 1 mixing form in Eq. (1.2) where the first column is proportional to (−2, 1, 1) T . The reason is simply that φ ν 3 ∝ (−2, 1, 1) T is an eigenvector of M ν in Eq. (3.1) with a zero eigenvalue corresponding to the first neutrino mass being zero. The reason for this is that M ν in Eq. (3.1) is a sum of two terms, the first being proportional to φ ν 1 φ ν 1 T and the second being proportional to φ 120 φ T 120 . Since φ ν 3 ∝ (−2, 1, 1) T is orthogonal to both φ ν 1 and φ 120 it is then clearly annihilated by the neutrino mass matrix, i.e. it is an eigenvector with zero eigenvalue. Therefore we immediately expect M ν in Eq. (3.1) to be diagonalised by the TM 1 mixing matrix where the first column is proportional to φ ν 3 ∝ (−2, 1, 1) T . Although the remainder of this subsection gives a perturbative diagonalisation of M ν in Eq. (3.1), we already know it must lead to TM 1 mixing exactly to all orders according to this general argument.
Mindful of the PDG phase conventions for the PMNS mixing matrix we write where c ij = cos θ ij and s ij = sin θ ij . Furthermore, P = diag(1, e i α 2 2 , e i α 3 2 ) is the Majorana phase matrix and P ′ = diag(e iδe , e iδµ , e iδτ ) is an unphysical phase matrix which is required to bring the PMNS matrix into PDG form.
In order to determine U PMNS we first calculate the eigenvalues of (M † ν M ν )/m 2 a which will be functions of ǫ and α. Requiring orthogonality, the corresponding eigenvectors can be obtained analytically. These eigenvectors, normalised to the unit length, comprise the unitary matrix which diagonalises M † ν M ν and thus also M ν . It can therefore be identified with P ′ U PMNS once the Majorana phases have been adjusted to give real masses in Eq. (3.4). Using ǫ as our expansion parameter we obtain the second order result, 3 sin α (only up to order ǫ) , (3.11) Note that the PMNS matrix has only one non-trivial Majorana phase as one of the neutrinos is exactly massless. These results are only slightly modified if we choose the (1, 0, 2) T alignment for the subdominant neutrino term: θ 23 → π 2 − θ 23 , δ → π + δ, δ e → π + δ e , and δ µ ↔ δ τ . All observables in the neutrino sector can be expressed in terms of m a , ǫ and α. Excluding Majorana phases (and the mass of the massless neutrino), this means that the model class makes three predictions which should be testable in future oscillation experiments since θ 13 is comparatively large.
It is useful to compare the above predictions to a general leading order parametrisation of the PMNS mixing matrix in the PDG convention in terms of deviations from TB mixing [17], where the deviation parameters s, a, r are defined as [17], (3.15) 3 We remark that the CP violating Dirac phase δ is only determined to order ǫ as it always appears together with sin θ13 which in turn is already of order ǫ. For completeness we also list the results for the unphysical phases δe = − ǫ 2 sin α(1 − 5ǫ cos α), δµ = − ǫ 2 sin α(3 − 7ǫ cos α), and δτ = π + ǫ 2 sin α(1 − ǫ cos α).
At leading order the above predictions can be expressed by a = r cos δ , where the predicted reactor angle may be compared to Eq. (1.1). 4 We emphasise that these predictions hold true for both the (1, 2, 0) T as well as the (1, 0, 2) T alignment. In both cases, with a suitable choice of phase convention, the leading order mixing matrix can be written in the form, where Eq. (3.18) corresponds to a small angle expansion of TM 1 mixing in Eq. (1.2). However, from the general argument given earlier in this subsection, we expect TM 1 mixing in Eq. (1.2) to be valid to all orders beyond the small angle approximation.

Numerical results
In this section we present the numerical results for the two CSD2 cases defined by Eqs. (2.12) and (2.14) in Figs. 1 and 2. We used random values for m a , m b and α as input and calculated with the Mixing Parameter Tools provided with the REAP package [19] the resulting neutrino masses, mixing angles and CP violating phases. We used then the recent global fit results from [2] for the solar and atmospheric neutrino mass squared differences and mixing angles deduced from the new reactor fluxes as a constraint. We have also checked that the numerical results agree well with the second order analytic results in Eqs. (3.6-3.13).
There are some interesting features of the plots. First of all, note that θ 13 can go up to more than 7 • in the 3σ interval. It is, however, more interesting to look at the 1σ regions. The atmospheric mixing angle θ 23 has an upper 1σ bound of 45 • , which is very restrictive for the (1, 2, 0) T alignment. Indeed, by this bound, θ 13 > 5 • is disfavoured in the (1, 2, 0) T case, while for the (1, 0, 2) T alignment values up to 6.4 • are still allowed, see upper panels of Fig. 1. This is due to the fact that the deviations from θ 23 = 45 • have opposite signs for both cases. Turning to the solar mixing angle, the 1σ region for θ 12 induces a lower bound on θ 13 of approximately 4.5 • , which is identical in both cases, see lower panels of Fig. 1.
It is also interesting to look at the phases in Fig. 2. In the (1, 2, 0) T alignment case a phase difference α, cf. Eq.  Fig. 2. Actually, the maximally CP violating cases δ = ±90 • are in both cases at the edge of the preferred regions. This is due to the fact that the corrections to a maximal atmospheric mixing are very small for α ≈ δ = ±90 • , see Eq. (3.8). Such a phase can emerge naturally in models with spontaneous CP violation from discrete symmetries [13].

PMNS-leptogenesis link
As has been noticed in [20][21][22], in models where TB mixing is realised via flavons which are orthogonal to each other, as in CSD or more generally in scenarios which satisfy the conditions of form dominance [22,23], the CP asymmetries for leptogenesis vanish.
On the contrary, in models with (1, 2, 0) T or (1, 0, 2) T vacuum alignment -since the flavon vevs of the model are now no longer orthogonal -the asymmetry does not vanish, rendering models of this type attractive for cosmology.
Furthermore, the two zero textures in Y ν imply a direct link between the CP violation for leptogenesis and the Dirac CP phase δ, as has been discussed for models with sequential dominance and hierarchical RH neutrinos in [15,20,24]. The produced baryon asymmetry Y B from leptogenesis in models with (1, 2, 0) T or (1, 0, 2) T vacuum alignment satisfies Since the baryon asymmetry Y B is positive, it follows that, in models with a fixed alignment and RH neutrino masses, leptogenesis requires δ in a specific range. In models where the "+" sign applies in Eq. (4.1) only the region around δ = 90 • generates the correct positive Y B , while in models where the "−" sign holds only the region around δ = 270 • is valid. Current global fits [3] favour sin δ being negative and this suggests that the negative sign is favoured in Eq. (4.1). According to the above discussion this suggests either the (1, 0, 2) T alignment with M B ≪ M A or the (1, 2, 0) T alignment with M A ≪ M B . The latter possibility corresponds to so called "light sequential dominance" which plays a special role in leptogenesis within the framework of two right-handed neutrino models as recently discussed in [25].

Summary and conclusions
Recently T2K have published evidence for a large non-zero reactor angle which, if confirmed, would exclude the tri-bimaximal mixing pattern. In this paper we have presented a model which fixes the reactor angle while preserving trimaximal solar mixing. In particular we have shown how a variant of trimaximal mixing, called TM 1 mixing in Eq. (1.2) with the solar angle given by sin θ 12 ≈ 1/ √ 3, results from an extension of constrained sequential dominance involving new vacuum alignments along the (1, 2, 0) T or (1, 0, 2) T directions in flavour space. We have shown that such alignments are naturally achieved using orthogonality, and may replace the role of the subdominant flavon alignment (1, 1, 1) T in constrained sequential dominance. We have proposed the first model in the literature of this kind leading to TM 1 mixing where the reactor angle is related to the ratio of the solar to the atmospheric neutrino masses, θ 13 = . We emphasise that the considered model is merely representative of a general strategy based on CSD2 for obtaining TM 1 mixing together with the above prediction for the reactor angle.
We have studied the phenomenological consequences of CSD2 both analytically and numerically. The analytic treatment confirms that TM 1 mixing results, at leading order, in a reactor angle which is predicted to be proportional to the ratio of the solar to the atmospheric neutrino masses, yielding θ 13 ∼ 5 • − 6 • , while the atmospheric angle is given by the sum rule θ 23 ≈ 45 • + √ 2θ 13 cos δ, where the leptonic Dirac CP phase δ is undetermined by CSD2, but experimentally preferred to lie in a range of 90 • − 130 • or 230 • − 270 • . The numerical results agree well with the second order analytic results, and demonstrate the full range of neutrino mixing parameters possible with CSD2, although this range could be extended in models which contain additional contributions from charged lepton mixing.
Finally we have seen that in CSD2 leptogenesis is unsuppressed due to the violation of form dominance and that the decay asymmetries feature a direct link between the CP phase for leptogenesis and the Dirac CP phase δ, with the produced baryon asymmetry Y B ∝ ± sin δ.
In conclusion, CSD2 leads to a highly predictive form of leptonic mixing, with the solar angle tightly constrained to its trimaximal value, the reactor angle predicted to be within the range of recent global fits, and the atmospheric angle correlated with the Dirac CP phase δ which is precisely equal to the leptogenesis phase. The large reactor angle indicated by T2K therefore opens up the exciting possibility of an early measurement of the low energy CP violating phase δ which is also responsible for the matter-antimatter asymmetry of the Universe within this class of models.