Towards a complete A4 × SU(5) SUSY GUT

We propose a renormalisable model based on A4 family symmetry with an SU(5) grand unified theory (GUT) which leads to the minimal supersymmetric standard model (MSSM) with a ℤ9 × ℤ6 symmetry provides the fermion mass hierarchy in both the quark and lepton sectors, while ℤ4R symmetry is broken to ℤ2R, identified as usual R-parity. Proton decay is highly sup-pressed by these symmetries. The strong CP problem is solved in a similar way to the Nelson-Barr mechanism. We discuss both the A4 and SU(5) symmetry breaking sectors, including doublet-triplet splitting, Higgs mixing and the origin of the μ term. The model provides an excellent fit (better than one sigma) to all quark and lepton (including neu-trino) masses and mixing with spontaneous CP violation. With the A4 vacuum alignments, (0, 1, 1) and (1, 3, 1), the model predicts the entire PMNS mixing matrix with no free pa-rameters, up to a relative phase, selected to be 2π/3 from a choice of the nine complex roots of unity, which is identified as the leptogenesis phase. The model predicts a normal neutrino mass hierarchy with leptonic angles θ13ι ≈ 8.7∘, θ12ι ≈ 34∘, θ23ι ≈ 46∘ and an oscillation phase δι ≈ − 87∘.


Introduction
The Standard Model (SM), although highly successful, leaves many unanswered questions in its wake such as: what (if anything) stabilises the recently discovered Higgs boson mass? Are the three known gauge forces unified into a simple gauge group which also explains charge quantisation? What is the origin of the three families of quarks and leptons and their pattern of masses, mixing and CP violation? Why is CP so accurately conserved by the strong interactions? The answers to such questions may help resolve longstanding cosmological puzzles such as the nature of dark matter and the origin of matter-antimatter asymmetry, both of which are unexplained within the SM.
In this paper we propose a realistic and fairly complete model capable of addressing all the above questions unanswered by the SM. The basic ingredients of our model are Supersymmetry (SUSY) together with an SU(5) Grand Unified Theory (GUT), flavoured by an A 4 family symmetry (for a review see e.g. [1,2]). The model is minimal in the sense that SU(5) is the smallest GUT group and A 4 is the smallest family symmetry group that admits triplet representations. Also, below the GUT scale, the model yields the JHEP06(2015)141 minimal supersymmetric standard model (MSSM) supplemented by a minimal two righthanded neutrino seesaw mechanism. The model is realistic in the sense that it provides a successful (and natural) description of the fermion mass and mixing spectrum, including spontaneous CP violation, while resolving the strong CP problem. It is fairly complete in the sense that GUT and flavour symmetry breaking are addressed, including doublet-triplet splitting, Higgs mixing and the origin of the MSSM µ term.
The model also allows a WIMP dark matter candidate due to the conserved MSSM R-parity, and permits matter-antimatter asymmetry via leptogenesis involving the two right-handed neutrinos. We shall show that the leptogenesis phase is equal to the single phase appearing in the neutrino mass matrix, providing a direct link between neutrino oscillations and matter-antimatter asymmetry, although we shall not discuss cosmological aspects any further in this paper.
We emphasise the predictive nature of the model in the lepton sector, where the entire PMNS matrix is predicted without any free parameters, up to a discrete choice of a single phase. Large lepton mixing is accounted for by the seesaw mechanism [3][4][5][6] with constrained sequential dominance (CSD) [7][8][9][10][11][12]). With a diagonal two right-handed neutrino mass matrix M R , the dominant right-handed neutrino ν atm R mainly responsible for the atmospheric neutrino mass m 3 has couplings to (ν e , ν µ , ν τ ) proportional to (0, 1, 1), while the subdominant right-handed neutrino ν sol R giving the solar neutrino mass m 2 has couplings to (ν e , ν µ , ν τ ) proportional to (1,3,1). These couplings, corresponding to the so called CSD3 scheme [13,14], originate from A 4 vacuum alignment. 1 The model consequently predicts a normal neutrino mass hierarchy, m 3 > m 2 m 1 = 0. As mentioned above, the lepton sector is controlled by a relative phase which is selected to be 2π/3, chosen from the nine complex roots of unity arising from spontaneous CP violation of a Z 9 × Z 6 discrete symmetry, by a mechanism proposed in [18]. Such a spontaneous CP violating scenario had been proposed previously in order to account for the smallness of CP violation in the soft SUSY sector [19,20]. We also employ a Z R 4 discrete R-symmetry (as the origin of MSSM R-parity, as in [21]) and a missing partner (MP) mechanism [22,23] for doublet-triplet splitting as recently advocated for flavoured GUTs in [24]. The model predicts very sparse lepton and down-type quark Yukawa matrices, with five texture zeroes, and Yukawa elements involving simple SU(5) Clebsch-Gordan (CG) ratios of 4/9 and 9/2 for the first and second families, with m τ /m b = 1 for the third family, all in excellent agreement with their experimental values run up to the GUT scale [25,26].
Quark mixing originates predominantly from a non-diagonal and naturally hierarchical up-type Yukawa matrix, provided by the broken Z 9 discrete family symmetry. Quark CP violation, however, comes exclusively from a single off-diagonal element in the down Yukawa matrix. By contrast, to excellent approximation, all lepton mixing and CP violation originates from the neutrino mass matrix, whose structure is also controlled by the A 4 family symmetry and the Z 6 symmetry via the CSD3 type vacuum alignment as described above [13].
Although there have been many attempts in the literature based on A 4 flavoured SU(5) SUSY GUTs (for an incomplete list see e.g. [27][28][29][30][31][32][33][34][35][36]), we would argue that none are JHEP06(2015)141 as successful or complete as the present one. For example, many of the previous models predicted mixing very close to tri-bimaximal and are by now excluded. Indeed the present model is the first one based on CSD3 capable of predicting all the lepton mixing parameters consistent with current data on lepton mixing [13,14] (see also [15][16][17]). The full literature on flavoured SUSY GUTs, i.e. which involve a (discrete) family symmetry, is quite extensive (for an incomplete list see e.g. [37][38][39][40][41][42][43][44][45][46]). The goal of all these models is clear: to address the questions left unanswered by the SM. It will take some time and (experimental) effort to resolve all these models. However the most promising models are those that make testable predictions while being theoretically complete and consistent.
While there are many different chiral superfields in this model, indeed almost exactly a hundred, it is important to note that we are explicitly presenting a renormalisable model. Any "non-renormalisable terms" generated below the Planck scale are required to have a specific well defined realization through multiple renormalisable terms involving heavy messenger fields that can be integrated out around the GUT scale. The respective effective theory after they are integrated out is actually more predictive than otherwise, with a normal neutrino mass hierarchy, a zero lightest neutrino mass, and all lepton mixing angles and CP phases predicted. The model presented here is amongst the most predictive and complete SUSY GUTs of flavour, consistent with current data.
The layout of the remainder of the paper is as follows: in section 2 we describe the superfields directly related to the SM fermions and neutrinos, as well as their Yukawa structures as imposed by the GUT and family symmetries when certain A 4 breaking vacuum expectation values (VEVs) are applied; we also perform a global fit to the parameters of the model and present our predictions for the lepton sector. In section 3 we describe the superfields that are responsible for breaking the family symmetries and how the required A 4 breaking VEVs arise. In section 4 several aspects related to the GUT are discussed, particularly how to break SU(5) and Z R 4 down to the MSSM with R-parity in a viable way (i.e. addressing doublet-triplet splitting, the origin of the µ term and proton decay). We also discuss the resolution to the strong CP problem. In section 5 we discuss the link between leptogenesis and the oscillation phase in this model. Finally in section 6 we summarise our main results and conclude. Appendix A summarises the A 4 conventions used in this paper, in the basis of [47,48].

The Yukawa sector of the model
The model involves an A 4 × SU(5) CP invariant superpotential at the GUT scale, where all symmetries, including CP, are spontaneously broken along supersymmetric flat directions, as discussed in sections 3 and 4. As already noted, it involves a further Z 9 × Z 6 discrete family symmetry as well as a Z R 4 discrete R-symmetry. The purpose of this section is to describe those aspects of the model pertaining to the Yukawa sector, i.e. relevant for understanding quark and lepton masses, mixing and CP violation. The flavour sector of the model is very important in our approach, since we make a serious attempt to understand and, where possible, predict the experimentally observable fermion masses and mixing matrices.

JHEP06(2015)141
In table 1 we show the matter superfields F , T i that contain the quarks and leptons, as well as the right-handed neutrino superfields N c i and double seesaw superfield Γ, all of which carry unit Z R 4 charge. Apart from the A 4 × SU(5) assignments of F ∼ (3, 5), T i ∼ (1, 10), N c i ∼ (1, 1), under Z 9 they transform as F ∼ 0, T i ∼ (5, 7, 0), N c i ∼ (7,8). Unlike the rest of the quarks and leptons, the right-handed neutrinos are further charged under Z 6 (as are some of the symmetry breaking scalars).
In table 1 we also display the six Higgs superfields, generally denoted H (but also Λ) which serve to break the SU(5) gauge symmetry. The two light MSSM Higgs doublet superfields H u and H d will emerge from H 5 and a mixture of H5 and H 45 by a mechanism discussed later. The superfield ξ which breaks Z 9 is particularly central to this theory, as it is responsible for both right-handed neutrino masses and the up-type quark mass hierarchy. Finally we have the θ i superfields which break Z 6 and help to control Dirac neutrino masses, and nine A 4 breaking triplet flavons generally denoted φ, with various vacuum alignments, responsible for large lepton mixing.
With these assignments, only the top quark gets a mass from a renormalisable Yukawa coupling H 5 T 3 T 3 (which has Z R 4 charge 2 as required for an allowed superpotential term). All the other quark and lepton Yukawa couplings must arise through higher order terms. This provides the basic reason why most of the SM (or strictly MSSM) Yukawa couplings appear to be so small. The observed hierarchy of Yukawa couplings between the three families will be explained via a discrete Z 9 version of the Froggatt-Nielsen mechanism [49], with powers of the low VEV of ξ controlling the hierarchy in the up-type quark sector, and also, in part, the smallness of the down quark and electron.
In order to enhance predictivity we need the messengers listed in table 2, which is the price we pay for having a renormalisable theory at the GUT scale. We denote these superfields either as fermion messengers, X i , or scalar messengers, Σ i , depending on whether they carry similar quantum numbers to, respectively, the quarks and leptons (with odd Z R 4 charge) or the symmetry breaking scalars (with even Z R 4 charge). The fermion messengers X i carry similar quantum numbers to down-type quarks and charged leptons (and neutrinos). Scalar messengers Σ i have quantum numbers similar to H 5 (the superfield that gives the top quark a renormalisable mass term). The Σ i messengers do not get VEVs, which means we need not consider the effect of diagrams with Σ i superfields in external legs to the masses of SM fermions.
The messengers group themselves in pairs of two superfields with a renormalisable bare mass coupling which respects all the symmetries. Their masses are therefore expected to be at or around the GUT scale. Although there will be in general distinct masses for different pairs, for simplicity and because they are all expected to be at a similar mass scale, we take masses of all such pairs to be M and set it equal to the GUT scale in our numerical estimates. We emphasise that the successful predictions of the model in the lepton sector (namely predicting the PMNS matrix) is independent of the specific values of these mass parameters.

Up quarks
Apart from the top quark mass, which originates from a lowest order Yukawa coupling, the remaining up-type quark Yukawa couplings appear from higher order terms that result from combining several renormalisable terms involving Σ i messengers and the GUT singlet superfield ξ. To be precise, the up-type quark Yukawa couplings arise from Σ i messenger tower diagrams shown in figure 1. For example, the most suppressed coupling arises from the first diagram in figure 1. Other less suppressed couplings arise from the diagrams where at the base one has the respective T i T j , with a shorter tower leading up to suppressed coupling, the renormalisable H 5 T 3 T 3 operator responsible for the top quark mass, is the last diagram in figure 1.
The effective superpotential responsible for the up-type Yukawa couplings is The resulting symmetric Yukawa matrix for up-type quarks is whereξ = ξ /M ∼ 0.1. The explicit form of Y u is given in eq. (2.18) and includes the coefficients u ij , which are O(1) and, by enforcing CP conservation at the GUT scale, necessarily real. Thus, the hierarchy of the up quark masses as well as the CKM mixing angles are given by powers ofξ. Due to the structure of this matrix, any phase introduced by ξ can be reabsorbed by appropriate redefinition of the three T i fields, so Y u does not contain a source of CP violation.

Down quarks, charged leptons and flavons
When considering the Yukawa structures of down quarks and charged leptons we must inevitably discuss A 4 triplet flavons. 2 The assignments of all the flavons under the family symmetries appear in table 1. Indeed, since the three families of F transform as a triplet of A 4 (see table 1), all T i H5F terms require a contraction with at least one A 4 triplet flavon to be invariant. From the diagrams shown in figure 2, integrating out the fermion messengers X, which acquire large masses as a result of either explicit mass terms or GUT scale Higgs VEVs, JHEP06(2015)141 we obtain effective operators of the form such that, apart from d 12 , the contraction appearing with T 1,2,3 isolates the respective F 1,2,3 family. This would lead to diagonal Yukawa structures if not for the additional term connecting T 1 (φ µ F ) (see figure 2d). The resulting effective Yukawa matrices are, schematically: where v Λ 24 and v H 24 are the respective VEVs of Λ 24 and H 24 (given in eq. (2.6)), and we include the subscripts LR to emphasise the role of the off-diagonal term to left-handed mixing from Y d . The off-diagonal term in Y e also provides a tiny contribution to left-handed charged lepton mixing θ e 12 ∼ m e /m µ which may safely be neglected. It also introduces CP violation to the CKM matrix via the phase of ξ .
Furthermore, because the underlying renormalisable theory is known, the diagrams in figure 2 are the only contributions for each family. The SU(5) contractions and associated CG coefficients appearing for each family are unique [24][25][26]. With the GUT scale symmetry breaking as discussed in section 4, each of the scalars here get a VEV with the group JHEP06(2015)141 structure: where the indices run a, b, c = 1, . . . , 5. This leads to the GUT scale prediction: The explicit forms of Y d and Y e , including CG and d ij coefficients, are given later in eq. (2.19) and eq. (2.20), respectively.

Neutrinos and CSD3
In order to obtain the CSD3 vacuum alignment in this model we couple the neutrinos to a set of flavons distinguished by the Z 6 symmetry. Of the superfields in table 1, only the righthanded neutrinos and some of the flavons are charged under this symmetry. For clarity, we relabel two of the flavon fields as φ atm ≡ φ 3 and φ sol ≡ φ 4 , to highlight their role in producing neutrino mixing. We also write N c atm ≡ N c 1 to denote the right-handed neutrino that dominantly leads to the atmospheric neutrino mass, and N c sol ≡ N c 2 as that which contributes mainly to the solar neutrino mass. The relevant terms in the superpotential giving neutrino masses are thus The flavons φ atm and φ atm gain VEVs according to the CSD3 alignment: where v atm and v sol are generally complex. Denoting the phases of VEVs as ρ i = arg(v i ), only the relative phase ρ atm − ρ sol between the VEVs is physically relevant, and is constrained to a discrete set of values, as discussed in section 3.3. The flavon ξ (already responsible for the up quark masses) is also acting as a Majoron by generating hierarchical right-handed neutrino masses. At the effective level, the Dirac terms result from coupling the neutrinos (and H 5 ) to φ atm and φ sol via the flavon θ 2 (an A 4 singlet carrying Z 6 charge). The corresponding diagrams with associated messengers appear in figure 3.
In turn, the Majorana mass term for N c atm is also non-renormalisable and we refer to the superfield Γ as the respective messenger. It couples only to N c atm and simply provides the non-renormalisable mass term for N c atm , suppressed relative to the mass of N c sol . As Γ has the quantum numbers of a third right-handed neutrino, one can also consider this field as mediating a double seesaw mechanism, responsible for the N c atm mass. The mixing term ξ 6 M 5 N c atm N c sol , though allowed by the symmetries, is absent as there is no combination of messengers able to produce it.
We write ξ = |v ξ |e iρ ξ , where ρ ξ is chosen from a discrete set of available phases, discussed in section 4.1 (see eq. (4.2)). This phase originates from the spontaneous CP violation of a discrete Abelian symmetry [18][19][20], in our case the Z 9 . We will now show that ρ ξ and ρ atm − ρ sol fix the relative phases within the effective neutrino mass matrix and consequently the leptonic mixing angles.
In a Supersymmetric (SUSY) model, the relevant terms in the superpotential giving neutrino masses, in the diagonal charged lepton basis, are where L i are three families of lepton doublets and the (CP conjugated) right-handed neutrinos N c atm and N c sol with real positive masses M atm and M sol do not mix. The structure of λ ν is determined by the vacuum alignments of φ atm and φ sol . The Dirac and Majorana matrices, derived by comparing the superpotential terms in eqs. (2.8) and (2.10) are For the see-saw mechanism we shall introduce a different convention for Yukawa and Majorana masses. The Yukawa matrices Y e , Y ν are defined in a LR convention by There is a simple dictionary between the seesaw basis and the SUSY basis, as follows: compared to the SUSY basis in eq. (2.8) used in leptogenesis calculations we see that Hence the neutrino matrices become, in the seesaw basis, (2.14) Seesaw produces the effective neutrino mass matrix and we have multiplied throughout by an overall phase which we subsequently drop, keeping only the (physical) relative phase η ≡ −ρ ξ + 2(ρ atm − ρ sol ), (2.16) where we recall the above definitions of phases, and that CP conservation at high energies ensures that y i and M are real. By arguments given in section 3.3 and section 4.1, we can restrict the physical phase η to a discrete choice, namely one of the nine complex roots of unity. The values η = ±2π/3 are preferred by CSD3 [13,14]. Note that the model predicts a normal neutrino mass hierarchy, namely m 3 > m 2 m 1 = 0, which will be tested in the near future. The sign of η has phenomenological significance, as it fixes the leptonic Dirac phase δ l . Specifically, a positive η uniquely leads to negative δ l , and vice versa [14]. As experimental data hints at δ l ∼ −π/2, the a posteriori preferred solution has positive η = +2π/3. The sign of η also has cosmological significance, as discussed in section 5. For example a positive η = +2π/3, together with the requirement that baryon asymmetry is positive, implies that the lightest right-handed neutrino should be N c should be somewhat heavier, which is the natural ordering in our model. 3 The conventions for Y ν,e and m ν differ, respectively, by overall Hermitian and overall complex conjugation compared to those used in the Mixing Parameter Tools package [50], which was used when performing global fits.

Full parameter fit
The structure of the Yukawa matrices and neutrino mass matrix is set by the theory, up to O(1) coefficients. The VEVs of the fields ξ, Λ 24 and H 24 are at or near the GUT scale, but otherwise undetermined. This freedom coincides with the choice of coefficients in the Yukawa matrices, providing no extra degrees of freedom in the determination of the Yukawas other than to provide the appropriate scale. The same is true for the flavon fields φ e , φ µ and φ τ , which provide the necessary hierarchy in the down-quark and charged lepton Yukawa sector. The neutrino matrix m ν is given in eq. (2.15). Letting v f represent the VEV of a field f , the Yukawa matrices are as follows: (2.20) As already remarked, the phases in Y u from powers of ξ = |v ξ |e iρ ξ can be removed by field redefinition. Without loss of generality we have rephased fields such that the only phase appearing in Y d and Y e is the phase ζ as shown in eqs. (2.19), (2.20), so all quark CP violation originates from the single phase ζ appearing in Y d 12 . In turn, ζ is determined by a combination of phases coming from various field VEVs; more precisely, ζ = ρ ξ −2ρ H 24 −ρ Λ 24 . As long as it is reasonably far from zero, it can produce the necessary CP violation. Different choices of ζ do not affect the goodness-of-fit, corresponding simply to different but equally valid choices of O(1) coefficients. For our fit we choose ζ = π/3. Note that the corresponding phase in Y e 21 does not contribute to leptonic CP violation, since this term does not affect left-handed mixing, to an accuracy of O(m e /m µ ).
To fit the real coefficients u ij , d ij , m a and m b , we propose a function χ 2 that relates the N physical predictions P i ({x}) for a given set of input parameters {x} to their current best-fit values µ i and their associated 1σ errors, denoted σ i , by The errors σ i are equivalent to the standard deviation of the experimental fits to a Gaussian distribution. For most parameters, this is essentially the case, with the exception of the (lepton) atmospheric angle θ l 23 . For a normal hierarchy, the distribution is roughly centered on θ l 23 = 45 • , while the best fit value is given by θ l 23 = 42.3 • . So as to not overstate the error for θ l 23 , we approximate its distribution by a Gaussian about 42.3 • , setting σ θ l 23 = 1.6 • . We now wish to minimise χ 2 . In this analysis, N = 18, corresponding to six mixing angles θ l ij (neutrinos) and θ q ij (quarks), the CKM phase δ q , nine Yukawa eigenvalues for the quarks and charged leptons, and two neutrino mass-squared differences ∆m 2 21 and ∆m 2 31 . In the lepton sector, we use the PDG parametrisation of the PMNS matrix [51] U PMNS = R l 23 U l 13 R l 12 P PDG in terms of s ij = sin θ l ij , c ij = cos θ l ij , the Dirac CP violating phase δ l and further Majorana phases contained in P PDG = diag(1, e i α 21 2 , e i α 31 2 ). Experimentally, the leptonic phase δ l is poorly constrained at 1σ (and completely unconstrained at 3σ), so is not fit, and left as a pure prediction of the model, as are the (completely unconstrained) Majorana phases α 21 and α 31 . As this model predicts only two massive left-handed neutrinos, i.e. m 1 = 0, one Majorana phase is zero, which we take to be α 31 = 0.
The running of best-fit and error values to the GUT scale are generally dependent on SUSY parameters, notably tan β, as well as contributions from SUSY threshold corrections. We extract the GUT scale CKM parameters and all Yukawa couplings (with associated errors) from [52] for judicious choices of tan β. In further reference to [52], we choose for the parameterη b parametrising the threshold corrections a valueη b = −0.24375; a non-zero value is required primarily to produce a necessary (small) difference in b and τ Yukawa couplings. Experimental neutrino parameters are extracted from [53]. All data is reproduced in tables 3 and 4.
Minimisation by differential evolution was performed in Mathematica, producing the set of O(1) input coefficients and the corresponding physical parameters seen in table 5, with an associated χ 2 = 7.98 (for tan β = 5) and χ 2 = 7.84 (for tan β = 10).
In this fit, the VEVs of ξ, Λ 24 , H 24 and the three φ e,µ,τ are fixed by hand in terms of the scale M , which is taken to be the GUT scale, i.e. M ≈ 3 × 10 16 GeV. Similarly, the Higgs doublet VEV enters only implicitly through m a and m b , but is understood to take (at the GUT scale) the value v H = 174 GeV. We set The value of v ξ is chosen to accommodate not only the fit to Y u parameters but also to control the µ-term, as discussed in section 4.2. Meanwhile the factor ∼3 split between v Λ 24 and v H 24 assists in establishing a hierarchy between the e and µ families. With the above numerical values for the VEVs, we get the following Yukawa matrices in terms only of O(1) coefficients and the complex phase ζ:   Table 3. Best fit values for quark and charged lepton parameters when run to the GUT scale as calculated in [52], with the SUSY breaking scale set at 1 TeV. We have included an overall contribution from threshold corrections corresponding toη b = −0.24375 which affects primarily the b quark Yukawa coupling y b . µ i represents the best-fit value and σ i the error, as defined in eq. (2.21).
It is worth reiterating that the neutrino mass matrix phase η can be forced to admit only phases coming from the nine complex roots of unity, essentially due to spontaneous CP violation with the Z 9 symmetry, where we select η = ±2π/3, both of which yield equally good χ 2 fits, with only minor adjustments to O(1) coefficients. The primary effect is in the prediction of δ l ; as previously stated, positive η corresponds to negative δ l . As this is preferred by experiment, the results presented in table 5 are for η = +2π/3. The fit also predicts the Majorana phases α 21 = 72 • and α 31 = 0.
In order to understand the significance of the χ 2 fit, and assess the strength of the model overall, it is prudent to enumerate the parameters and predictions of the model. The nominal parameter count at the GUT scale is very large, owing to the diverse field content. However, at the scale where we are able to make predictions, many of these parameters combine to give a constrained set of free parameters that need to be determined. Notably, the VEVs of Higgs and flavon fields such as those given in eq. (2.22) do not constitute true degrees of freedom, as they can be absorbed by redefining other parameters.

JHEP06(2015)141
Parameter (from [53]) µ i ± 1σ 3σ range Relevant parameters that require consideration include: six u ij , four d ij , masses m a and m b , phases η and ζ, the threshold factorη b , and tan β, for a total of N I = 16 input parameters. However three of these parameters, namely tan β, η and ζ, are fixed prior to the fit, with the latter two phases restricted to discrete choices, as discussed previously. Finally, the factorη b affects only the coupling y b and is fitted by hand. As mentioned earlier, the model fits N = 18 observables, including nine Yukawa eigenvalues, two neutrino mass squared differences, six mixing angles and the quark CP phase. In addition the model predicts the leptonic CP phase δ l , two Majorana phases (one of which is zero) and a massless physical neutrino.

A 4 symmetry breaking and the flavon vacuum
In order to address A 4 symmetry breaking we need to address three aspects of the flavon vacuum: what drives some flavons to have VEVs a few orders of magnitude below the GUT scale, what determines their vacuum alignment, and what fixes the relative vacuum phase ρ atm − ρ sol ≡ arg(v atm v * sol ) (and consequently the physical phase η). In this section we consider each of these issues in turn.

Driving the flavon vacuum expectation values
The flavon φ VEVs are driven by radiative breaking [54] (see e.g. [55] for a recent review). The soft squared mass terms appear as m 2 and become negative at the required scales Λ i < M , lifting the respective flat direction to a few orders of magnitude below the GUT scale for each flavon: φ i φ † i ∼ Λ 2 i . Hierarchies of VEVs are thus naturally expected due to the logarithmic nature of this mechanism. The precise symmetry breaking scale Λ i for each field φ i depends on otherwise undetermined parameters in the model which are different for each flavon, such as superpotential terms involving messengers. 4   Table 5. Fitted input quark Yukawa coefficients u ij and d ij (arranged by their position in the Y u and Y d matrices, respectively) and neutrino mass parameters m a and m b , and associated physical parameters produced by minimising the function χ 2 . We choose ζ = π/3. With tan β = 5, the fit gives χ 2 = 7.98, while tan β = 10 gives χ 2 = 7.84, both very good fits. The largest single contribution to χ 2 is from the fit to the atmospheric angle θ l 23 . These results are given with η = +2π/3. The non-zero Majorana phase is also predicted to be α 21 = 72 • , and is insensitive to tan β, as indeed are all the mixing angles and phases. Therefore a hierarchy for such flavon VEVs is generated and remains stable due to radiative breaking [56][57][58][59][60][61]. On the other hand those for v e , v µ , v τ in eq. (2.22) arise from F -terms as discussed later.

Flavon vacuum alignment
Thus far we have assumed that the A 4 triplet VEVs are aligned in special directions. In this section we describe how these directions are obtained by the superpotential terms allowed by the symmetries. In doing this, the role of Z 6 becomes clearer. The driving sector, a set JHEP06(2015)141 of superfields A i and O ij with Z R 4 charge 2, is listed fully in table 7. 5 The inclusion of the Z 6 symmetry is necessary because the driving superpotential, responsible for aligning the flavons, needs to have each driving field isolated, as shown below (see also [13]): The additional Z 6 charges ensure each term is separated from all others, leading to an array of vanishing F -term conditions that force mutual orthogonality conditions between many of the vacuum alignments. Since this was fully discussed in [13], we need only state the results here, namely that eq. (3.1) leads to the following vacuum alignment patterns: The role of the VEVs (containing two zero entries) of the superfields φ e,µ,τ appearing in eq. (2.3) was already discussed in section 2.2. Meanwhile, the role of the VEVs of the flavons φ 3,4 (redubbed φ atm,sol ) was described in section 2.3. It is the special structure of these vacuum alignments, combined with the phase of η in the neutrino mass matrix, that leads to the very successful prediction of the leptonic mixing angles (as described in section 2.4). The remaining VEVs are not directly relevant to the masses and mixings of SM fermions, but help shape the VEVs of φ atm and φ sol .

Flavon vacuum phases
With the direction of the A 4 triplet flavons φ fixed, we turn now to a discussion of how to fix the relative phase ρ atm − ρ sol ≡ arg(v atm v * sol ) to a discrete choice. We present a mechanism which does this by adding a number of fields P that are A 4 and SU(5) singlets, also given in table 7, and which resemble the O fields except they do not force orthogonality between the flavons φ.

JHEP06(2015)141
These fields and their respective charge assignments result in the following invariant superpotential terms: where each term technically has an associated real coupling λ which is O(1) and may be made positive by field redefinitions. We omit these for simplicity as they have no effect on the general argument presented here, with one caveat: the two superfields P 1,2 33 have exactly the same quantum numbers but different λ couplings to flavons. Due to this duplication there are two independent relations between the flavon VEVs involving different λ couplings which leads to an additional constraint on the phases of the respective VEVs. Exact values of these λ are not specified; it suffices that they are not equal.
Furthermore, the primary role of the SU(5) adjoint fields Z 2 and Z 3 which couple to P µµ is in the GUT breaking mechanism (as discussed in section 4.1). Their phases are fixed separately by other superpotential terms.
We begin the analysis of these terms by noting they do not affect the alignments of the flavons φ. The corresponding F -terms for each field P ij produces a set of coupled equations that admit a solution where none of the A, O, and P fields but all the flavons obtain a VEV. Omitting the (real, positive, O(1)) λ coefficients, these VEVs have the structure: Regarding the magnitudes of the VEVs, two comments are in order. We assumed above that M sets the scale of the VEV of φ e , which is in contradiction with our previous assumption that it be O(10 −3 )M . This violates our simplifying assumption that all mass scales are equal, and demonstrates that some spectrum of mass scales is in fact required in this model. As for the VEV v 3 , as discussed earlier, it is driven to a specific scale Λ 3 radiatively [56][57][58][59][60][61]. Writing ρ i ≡ arg(v i ), this VEV structure gives (up to multiples of π) the phase relation where n is an integer, and similar relations for the other flavons as linear combinations of ρ 3 , (ρ Z 2 + ρ Z 3 ) and multiples of 2π/3. This is an important equation since it fixes the relative phase ρ 3 − ρ 4 = ρ atm − ρ sol in terms of 1 6 (ρ Z 2 + ρ Z 3 ). As discussed in section 4.1,

GUT scale symmetry breaking, proton decay and the strong CP problem
In this section we discuss the aspects of the model related to grand unification, starting with how the R-symmetry and the GUT gauge group are spontaneously broken. We refer to the superfields involved as the scalar sector; they are listed in tables 1 and 6. We then describe the details of the MP mechanism, and finish this section with an analysis that justifies the absence of dangerous proton decay operators in the model. (5) and Z R 4 breaking As previously discussed, the Υ messengers form pairs; their mass scale, unprotected by any symmetry, is near the highest scale of the theory, which we represent generically as M . The GUT breaking superpotential with non-renormalisable terms is then 6

SU
(4.1) We have five GUT adjoint superfields, three of which (the Z fields) are charged by 2 and two (Λ 24 and H 24 ) by 0 under the R-symmetry. Also appearing in W GUT is the Majoron ξ, the GUT singlet field which we have seen is involved in giving mass to several SM fermions and whose VEV breaks lepton number by giving the right-handed neutrinos their Majorana mass. The supressions of the non-renormalisable terms in eq. (4.1) come precisely from the mass of the Υ messengers, as displayed in figure 4. W GUT has a non-trivial minimum: where all the adjoint scalars get a VEV of the form Φ 24 = v Φ 24 diag(2, 2, 2, −3, −3). By themselves, the F -terms associated with W GUT also allow a trivial minimum where the 6 A renormalisable term of the form Z2H24Π2, allowed by the symmetries, has been dropped to make the discussion more transparent. This term mixes the VEVs of the GUT breaking scalars with the ones in the MP mechanism so they should be naturally around the same scale (M ∼ MGUT). Beyond this, its practical effect is minimal: the fields obtain VEVs with or without this term. Since the VEVs get very complicated when this "mixing" term is included, we ignore it for simplicity, simply bearing in mind that VEVs from both sets of fields are related.

Field
Representation  magnitude of each VEV vanishes. But after SUSY is broken and we consider the effects of the small contribution from radiative breaking [54] to the scalar components of the GUT breaking superfields (as we did in section 3 for the A 4 breaking flavons), the stationary point with vanishing magnitudes is no longer a minimum due to the radiatively induced negative squared mass term. To a very good approximation the true minima are given by the magnitudes in eq. (4.2), which are now a lower energy state than the trivial F -term solution.
We conclude that eq. (4.1) can generate GUT and R-symmetry breaking at high scale, with Z R 4 broken to Z R 2 (standard R-parity preserved) by the Z i VEVs. Because Z R 4 is broken at a high scale, the no-go theorem from [62] does not apply to our model and we verified that all the components of the SU(5) adjoints acquire GUT scale masses. Figure 4. Diagrams for non-renormalisable potential terms.
A slightly unappealing issue with W GUT is that the minimum requires some non-O(1) choice of λ parameters if we are to obtain a hierarchy between the VEVs of H 24 and Λ 24 , and an appropriate value for ξ /M as shown in eq. (2.22). These requirements come from the successful fit to up and down quark and charged lepton masses, as discussed in sections 2.1, 2.2 and 2.4, and partly also for the µ term, as will be discussed shortly. However, since the messengers will in general have different masses (recall we set them all equal to M only for simplicity), the λ parameters need not be as hierarchical as eq. (4.2) appears to indicate. For example, if the masses of messengers Σ are slighly larger than the GUT scale masses of messengers Υ, this would allow all λ to be O(1).
We note also that, although we are considering a situation where the superpotential parameters (the M and λ couplings) are real due to CP conservation, the VEVs of the GUT breaking scalars may be complex, since they depend on n th order roots of real numbers. As has been shown in sections 2.3 and 3.3, the phases of the fields ξ, Z 2 and Z 3 are relevant for establishing the physical phase η in the neutrino mass matrix, which controls neutrino masses and mixing. We see immediately that ρ ξ = 2πk 9 , for integer k, i.e. one of nine roots of unity. While ρ Z 2 and ρ Z 3 individually can be any of six roots, originating in the factor λ 1/6 5 , their product Z 2 Z 3 cancels this factor such that the largest root is a third, giving ρ Z 2 + ρ Z 3 = 2πk 3 , for integer k .

Doublet-triplet splitting, Higgs mixing and the µ term
Given that we have a number of GUT representations containing (SM gauge group) SU (2) doublets and triplets we turn now to a brief discussion of how doublet-triplet splitting is achieved in this model. Although one could alternatively introduce the double MP mechanism [63] as demonstrated in [24], here we limit ourselves to describing how the MP mechanism is implemented with the fields listed in table 6.
We have a superpotential  With eq. (4.4), the MP mechanism proceeds through the superpotential: (4.5) The very high order non-renormalisable term at the end arises through the Σ messengers already listed in table 2 (half of which participate in the up quark Yukawa terms, as discussed in section 2.1 and illustrated in figure 1). Strictly speaking this term does not participate in splitting the masses of SU(2) doublets and triplets, rather it is the source of the µ term in our model, as shown below.
The terms in W MP generate mixing between the 45s and 5s of SU (5). The mass matrix for the SU(2) triplets contained in the 45s and 5s is: (4.8) The usual MSSM term µH d H u comes from this mechanism with: where v Π 1 provides the necessary Z R 4 breaking. Using the fit from eq. (2.22) we see that ξ 8 ∼ 1.6×10 −10 M GUT . If we choose the couplings at the vertices of the tower that generates the ξ 8 term to be ∼ 0.5 we may get a term µ ∼ O(10 2 − 10 3 ) GeV without any fine-tuning.

Proton decay
A classic problem in GUT theories, and in particular those based on SU (5), is the prediction of excessively fast proton decay. The most dangerous processes come from the "dimension 5" operators T T T F (for a discussion of dimension 6 operators we refer the reader to [24]). The "dimension 5" operators are forbidden by the symmetries of the model, but related higher order operators of the following form are allowed: where the extra superfields shown are needed for such terms to be invariant under the symmetries. Since we are working with the renormalisable theory, in order for this type of effective term to be present at the GUT scale, there must be messengers allowing them. In this case, an analysis of the SU(5) index structure revels there should either be messengers that are SU(5) 10, or SU(5) 5 that are also charged under Z R 4 . As one can confirm from table 2, our model has neither: 10 messengers were not used, and the 5 messengers are all neutral under Z R 4 . We conclude therefore that our symmetry content, together with the existing set of messengers, do not allow any such GUT scale suppressed operators that would lead to excessively fast proton decay to be generated. The operators in eq. (4.10) may in principle be generated by physics at the Planck scale, with the scale M replaced by the Planck mass, leading to highly suppressed proton decay.

Strong CP problem and the Nelson-Barr resolution
We first recall the strong CP problem, namely that the physical angle θ = θ QCD − θ q , where θ QCD multiplies the topological gluon term (g 2 s /32π 2 )GG and θ q = arg det(Y u Y d ), is limited to be θ < 10 −10 by the non-observation of the neutron EDM [51,64,65]. The origin of such a small number, θ < 10 −10 , is commonly called the strong CP problem. It is interesting to compare this with the CP violation related to the weak interaction in the quark sector; the relevant quantity is the Jarlskog invariant

JHEP06(2015)141
which, when compared to data, is required to be non-vanishing, and indeed in the standard parameterisation, requires a large phase angle δ q ∼ 1.
It turns out that our model resolves the strong CP problem without relying on the introduction of axions. Unlike the axion solution, which requires a global U(1) symmetry with a colour anomaly, we shall rely on the fact that the high energy theory conserves CP, ensuring that θ QCD = 0. CP is then spontaneously broken in such a way as to yield δ q ∼ 1 while maintaining θ < 10 −10 and in particular θ q < 10 −10 . How it achieves this feat can be seen from eq. (2.18) where Y u is real, while the structure of Y d in eq. (2.19) gives it a real determinant. This is due to the lack of a Yukawa term Y d 21 , meaning the coupling Y d 12 (which is the only complex Yukawa coupling) does not appear in the determinant of Y d . Therefore there are no contributions to θ q even after spontaneous CP breaking. This is similar to the Nelson-Barr mechanism [66][67][68][69], where the triangular form of Yukawa matrices was proposed, although in our model θ q vanishes due to the triangular form of Y d only, with Y u being non-triangular and real.
For a successful resolution of the strong CP problem, we must ensure that no higher order corrections to the Yukawa matrices arise which would violate the bound θ q < 10 −10 . The main focus is on the Yukawa coupling Y d 21 which is zero at leading order but which may in principle receive higher order corrections, violating the triangular structure. However in our model such higher order corrections are absent at the field theory level with the specified messenger sector. This entry in the Yukawa matrix would arise from the coupling of the bilinears T 2 H5 ,45 to the bilinear φ e F . Since these terms are non-renormalisable, we require messengers to form them. The messengers that could produce such terms are the X i fields in table 2. With these messengers, the only allowed connection to φ e F is T 1 H5 (contributing to Y d 11 ), thus forbidding the Y d 21 term, even when allowing for all higher-order corrections. Therefore the specified model has no strong CP violation arising from Y d 21 since the required operators are not generated at the field theory level.
It is also important to consider the effect of higher order corrections arising from the Planck scale, since such operators only have to respect the symmetries of the model, and do not require the specified messenger sector to generate them. The biggest contribution would come from the term 7 With a general choice of phase, such a term would lead to θ q ∼ 10 −4 which is far too big. However the contribution to θ q may be avoided by a judicious choice of GUT breaking phases. As stated in section 2.4, the physical phase in the down quark Yukawa matrix is ζ = ρ ξ −2ρ H 24 −ρ Λ 24 . The new Planck suppressed term has a phase ζ = −ρ ξ +2ρ Λ 24 . Choosing a relation between phases 2ρ H 24 = ρ Λ 24 , then ζ = −ζ and this way the contribution to θ q vanishes. This happens for one in three cases. The next biggest contribution comes from a term

JHEP06(2015)141
giving θ q ∼ 10 −14 which is several orders of magnitude below the current experimental bound. Any other Planck suppressed terms allowed by the symmetries are further suppressed so we need not consider them. Therefore the model may resolve the strong CP problem even in the presence of Planck scale operators controlled only by symmetry. Finally, extra contributions may come from SUSY breaking terms. If we assume that there is no extra CP violation in this sector, which is controlled by the spontaneously CP violating flavons, the SUSY flavour problem is under control and such contributions to θ are also expected to be negligible [70].

The leptogenesis link
The link between leptogenesis and the PMNS matrix was first studied for sequential dominance in [71]. In the original form of CSD, the columns of the Dirac mass matrix in the flavour basis were orthogonal to each other and consequently the CP asymmetries for cosmological leptogenesis [72,73] vanished [74][75][76][77][78]. In this model, leptogenesis does not vanish since the columns of the Dirac mass matrix in the flavour basis are not orthogonal. Interestingly, since the seesaw mechanism in this model with two right-handed neutrinos only involves a single phase η = 2π/3, both the leptogenesis asymmetries and the neutrino oscillation phase must necessarily originate from this phase, providing a direct link between the two CP violating phenomena in this model.
Following the arguments in [74], the produced baryon asymmetry Y B from leptogenesis in the seesaw model in eq. (2.11) satisfies where the "+" sign applies to the case M atm M sol and the "−" sign holds for the case M sol M atm . Since the observed baryon asymmetry Y B is positive, it follows that, for M atm M sol , we must have sin η to be positive, while for M sol M atm we must have sin η to be negative. We have seen that positive η is associated with negative δ l and vice versa. Although the global fits do not distinguish the sign of η, the present hint that δ l ∼ −π/2 would require positive η, then in order to achieve positive Y B we require M atm M sol , which is natural in our model, corresponding to "light sequential dominance", where successful leptogenesis may be achieved in the two right-handed neutrino model as discussed in [79].

Conclusion
We have presented here a fairly complete realisation of an SU(5) GUT flavoured with A 4 , which leads to the MSSM plus two right-handed neutrinos below the GUT scale. The A 4 family symmetry unifies the three families of 5-plets F and its vacuum alignment determines the Yukawa matrices. In addition a Z 9 × Z 6 symmetry provides the mass hierarchy and controls spontaneous CP violation in both the quark and lepton sectors while a Z R 4 symmetry is broken to Z R 2 , identified as the usual R-parity. Proton decay is under control in this model, with the symmetries forbidding dangerous dimension-5 operators,

JHEP06(2015)141
and similar (but higher order) operators being very suppressed. The strong CP problem is resolved in a similar way to the Nelson-Barr mechanism. The model is highly predictive and satisfies many distinct (and non-trivial) phenomenological requirements.
Imposing CP at the high scale is an important feature of the model. If we do not impose CP then all couplings become complex, leading to all VEVs having undetermined phases. In particular the phase η, present in both neutrino mixing and leptogenesis, is no longer restricted to a discrete choice. However the link between leptogenesis and low energy phenomenology remains. On the other hand we would no longer solve the strong CP problem.
We highlight the ubiquitous nature of the flavon field ξ across all the sectors of the model: it triggers spontaneous CP violation in both the quark and lepton sectors, generates up-type quark mass hierarchies and CKM mixing, explains the smallness of down quark and electron masses and breaks lepton number, providing the hierarchy between solar and atmospheric right-handed neutrino masses. In addition, the ξ field is responsible for generating the small (complex) µ term. The phase of the ξ field VEV also contributes to the relative phase η appearing in the neutrino mass matrix. This phase, the only one appearing in the neutrino mass matrix and in the formula for the baryon asymmetry of the universe, provides a direct link between the PMNS phase δ l and leptogenesis.
We emphasise that the entire PMNS matrix is predicted without any free parameters, up to nine choices for η, where we select η = +2π/3, since it is preferred by comparing CSD3 to data. The required vacuum alignments are provided from A 4 . The model effectively serves to yield the CSD3 scheme (with two right-handed neutrinos) within a fully working and viable SUSY GUT of flavour in which all quark and lepton masses and mixings are successfully described. Indeed the model provides an excellent fit (better than one sigma) to all quark and lepton (including neutrino) masses, mixing and CP violation. All fermion mass hierarchies are understood in the sense that purportedly O(1) couplings indeed contain no strong hierarchies. However the most immediate predictions of the model are those of CSD3 with two right-handed neutrinos and η = +2π/3, namely a normal neutrino mass hierarchy with m 1 = 0, a reactor angle of θ l 13 ≈ 8.7 • , a solar angle θ l 12 ≈ 34 • , close to maximal atmospheric mixing θ l 23 ≈ 46 • and almost maximal leptonic CP violation, with an oscillation phase δ l ≈ −87 • consistent with the current hint δ l ≈ −π/2.
The reason why the field content is so large is that the model is fairly complete. In particular it is renormalisable at the GUT scale, which requires a large explicit field content including many heavy messenger superfields. It also addresses many aspects relevant both to a GUT and to family symmetry models (stopping short of discussing the details of SUSY breaking and its string theory completion). In particular, the A 4 , SU(5) and R-symmetry symmetry breaking sectors all require large field content. For example, the GUT symmetry is broken by an explicit superpotential at the GUT scale, including doublet-triplet splitting via a missing partner mechanism (leaving no light exotic degrees of freedom at the low scale), Higgs mixing and the origin of the MSSM µ term of the right order of magnitude.
Despite the many successes of the model, there are inevitably several important issues that lie beyond the scope of this paper. To take one example, we do not discuss GUT scale threshold corrections, which will be important in maintaining successful gauge coupling JHEP06(2015)141 unification in the presence of many fields, including colour triplets, at the GUT scale. In fact all the additional superfields in non-trivial representations of the gauge group may have masses at or above the GUT scale. Another important issue is that of the low energy superpartner spectrum in this model. While we expect SUSY induced flavour changing to be under control for the 5-plets, which are unified into an A 4 triplet, this is not the case for the 10-plets T i which are singlets of A 4 , leading to flavour violation in the super-CKM basis. It would be interesting to study this in the future.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.