Fully constrained Majorana neutrino mass matrices using Σ(72×3)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\varSigma (72\times 3)}$$\end{document}

In 2002, two neutrino mixing ansatze having trimaximally mixed middle (ν2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _2$$\end{document}) columns, namely tri-chi-maximal mixing (TχM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\chi \text {M}$$\end{document}) and tri-phi-maximal mixing (TϕM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\phi \text {M}$$\end{document}), were proposed. In 2012, it was shown that TχM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\chi \text {M}$$\end{document} with χ=±π16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi =\pm \,\frac{\pi }{16}$$\end{document} as well as TϕM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\phi \text {M}$$\end{document} with ϕ=±π16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi = \pm \,\frac{\pi }{16}$$\end{document} leads to the solution, sin2θ13=23sin2π16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sin ^2 \theta _{13} = \frac{2}{3} \sin ^2 \frac{\pi }{16}$$\end{document}, consistent with the latest measurements of the reactor mixing angle, θ13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{13}$$\end{document}. To obtain TχM(χ=±π16)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\chi \text {M}_{(\chi =\pm \,\frac{\pi }{16})}$$\end{document} and TϕM(ϕ=±π16)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {T}\phi \text {M}_{(\phi =\pm \,\frac{\pi }{16})}$$\end{document}, the type I see-saw framework with fully constrained Majorana neutrino mass matrices was utilised. These mass matrices also resulted in the neutrino mass ratios, m1:m2:m3=2+21+2(2+2):1:2+2-1+2(2+2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_1:m_2:m_3=\frac{\left( 2+\sqrt{2}\right) }{1+\sqrt{2(2+\sqrt{2})}}:1:\frac{\left( 2+\sqrt{2}\right) }{-1+\sqrt{2(2+\sqrt{2})}}$$\end{document}. In this paper we construct a flavour model based on the discrete group Σ(72×3)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma (72\times 3)$$\end{document} and obtain the aforementioned results. A Majorana neutrino mass matrix (a symmetric 3×3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3\times 3$$\end{document} matrix with six complex degrees of freedom) is conveniently mapped into a flavon field transforming as the complex six-dimensional representation of Σ(72×3)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma (72\times 3)$$\end{document}. Specific vacuum alignments of the flavons are used to arrive at the desired mass matrices.

In this paper we construct a flavour model based on the discrete group Σ(72 × 3) and obtain the aforementioned results. A Majorana neutrino mass matrix (a symmetric 3 × 3 matrix with six complex degrees of freedom) is conveniently mapped into a flavon field transforming as the complex sixdimensional representation of Σ (72 × 3). Specific vacuum alignments of the flavons are used to arrive at the desired mass matrices.

Introduction
The neutrino mixing information is encapsulated in the unitary PMNS mixing matrix which, in the standard PDG parameterisation [1], is given by where s i j = sin θ i j , c i j = cos θ i j . The three mixing angles θ 12 (solar angle), θ 23 (atmospheric angle) and θ 13 (reactor angle) along with the C P-violating complex phases (the Dirac phase, δ, and the two Majorana phases, α 21 and α 31 ) parameterise U PMNS . In comparison to the small mixing angles observed in the quark sector, the neutrino mixing angles are found to be relatively large [2]: sin 2 θ 13 = 0.01934 → 0.02392.
The values of the complex phases are unknown at present. Besides measuring the mixing angles, the neutrino oscillation experiments also proved that neutrinos are massive particles. These experiments measure the mass-squared differences of the neutrinos and currently their values are known to be [2], Several mixing ansatze with a trimaximally mixed second column for U PMNS , i.e. |U e2 | = |U μ2 | = |U τ 2 | = 1 √ 3 , were proposed during the early 2000s [3][4][5][6][7]. Here we briefly revisit two of those, the tri-chi-maximal mixing (Tχ M) and the tri-phi-maximal mixing (TφM), 1 which are relevant to Table 1 The standard PDG observables θ 13 , θ 12 , θ 23 and δ in terms of the parameters χ and φ. Note that the range of χ as well as φ is − π 2 to + π 2 . In TχM (TφM), the parameter χ (φ) being in the first and the fourth quadrant correspond to δ equal to + π 2 (0) and − π 2 (π ), respectively sin 2 θ 13 sin 2 θ 12 sin 2 θ 23 δ TχM Both Tχ M and TφM have one free parameter each (χ and φ) which directly corresponds to the reactor mixing angle, θ 13 , through the U e3 elements of the mixing matrices. The three mixing angles and the Dirac C P phase obtained by relating Eq. (1) with Eqs. (7), (8) are shown in Table 1.
In Tχ M, since δ = ± π 2 , C P violation is maximal for a given set of mixing angles. The Jarlskog C P-violating invariant [10][11][12][13][14] in the context of Tχ M [5] is given by On the other hand, TφM is C P conserving, i.e. δ = 0, π, and thus J = 0. Since the reactor angle was discovered to be non-zero at the Daya Bay reactor experiment in 2012 [15], there has been a resurgence of interest [16][17][18][19][20][21][22][23][24][25][26][27] in Tχ M and TφM and their equivalent forms. For any C P-conserving (δ = 0, π) mixing matrix with non-zero θ 13 and trimaximally mixed ν 2 column, we can have an equivalent parameterisation realised using the TφM matrix. Here the "equivalence" is with respect to the neutrino oscillation experiments. The oscillation scenario is completely determined by the three mixing angles and the Dirac phase (Majorana phases are not observable in neutrino oscillations), i.e. we have a total of four degrees of freedom in the mixing matrix. If we assume C P conservation and also maximal mixing [4]. By this notation, both TχM and TφM fall under the category of T M 2 . To be more specific, T M 2 , which breaks C P maximally, is TχM and T M 2 , which conserves C P is TφM. assume that the ν 2 column is trimaximally mixed, then there is only one degree of freedom left. It is exactly this degree of freedom which is parameterised using φ in TφM mixing. Similarly any mixing matrix with δ = ± π 2 , θ 13 = 0 and trimaximal ν 2 column is equivalent to Tχ M mixing.
In 2012 [22], shortly after the discovery of the non-zero reactor mixing angle, it was shown that Tχ M (χ =± π 16 ) as well as TφM (φ=± π 16 ) results in a reactor mixing angle, sin 2 θ 13 = 2 3 sin 2 π 16 = 0.025, (10) consistent with the experimental data. The model was constructed in the Type-1 see-saw framework [28][29][30][31]. Four cases of Majorana mass matrices were discussed: where M Maj is the coupling among the right-handed neutrino fields, i.e. (ν R ) c M Maj ν R . In Ref. [22], the mixing matrix was modelled in the form with ω = e i 2π 3 andω = e -i 2π 3 , in which the 3 × 3 trimaximal contribution came from the charged-lepton sector. U ν , on the other hand, was the contribution from the neutrino sector. The four U ν s vis à vis the four Majorana neutrino mass matrices given in Eqs. (11) and (12), gave rise to Tχ M (χ =± π 16 ) and TφM (φ=± π 16 ) , respectively. All the four mass matrices, Eqs. (11), (12) . Due to the see-saw mechanism, the neutrino masses become inversely proportional to the eigenvalues of the Majorana mass matrices, resulting in the neutrino mass ratios . (14) Using these ratios and the experimentally measured masssquared differences, the light neutrino mass was predicted to be around 25 meV. In this paper we use the discrete group Σ(72 × 3) to construct a flavon model that essentially reproduces the above results. Unlike in Ref. [22] where the neutrino mass matrix was decomposed into a symmetric product of two matrices, here a single sextet representation of the flavour group is used to build the neutrino mass matrix. A brief discussion of the group Σ(72×3) and its representations is provided in Sect. 2 of this paper. Appendix A contains more details such as the tensor product expansions of its various irreducible representations (irreps) and the corresponding Clebsch-Gordan (C-G) coefficients. In Sect. 3, we describe the model with its fermion and flavon field content in relation to these irreps. Besides the aforementioned sextet flavon, we also introduce triplet flavons in the model to build the charged-lepton mass matrix. The flavons are assigned specific vacuum expectation values (VEVs) to obtain the required neutrino and charged-lepton mass matrices. A detailed description of how the charged-lepton mass matrix attains its hierarchical structure is deferred to Appendix B. In Sect. 4, we obtain the phenomenological predictions and compare them with the current experimental data along with the possibility of further validation from future experiments. Finally, the results are summarised in Sect. 5. The construction of suitable flavon potentials which generate the set of VEVs used in our model is demonstrated in Appendix C.

The group Σ(7× 3) and its representations
Discrete groups have been used extensively in the description of flavour symmetries. Historically, the study of discrete groups can be traced back to the study of symmetries of geometrical objects. Tetrahedron, cube, octahedron, dodecahedron and icosahedron, which are the famous Platonic solids, were known to the ancient Greeks. These objects are the only regular polyhedra with congruent regular polygonal faces. Interestingly, the symmetry groups of the platonic solids are the most studied in the context of flavour symmetries too -A 4 (tetrahedron), S 4 (cube and its dual octahedron) and A 5 (dodecahedron and its dual icosahedron). These polyhedra live in the three-dimensional Euclidean space. In the context of flavour physics, it might be rewarding to study similar polyhedra that live in three-dimensional complex Hilbert space. In fact, five such complex polyhedra that correspond to the five Platonic solids exist as shown by Coxeter [32]. They are 3{3}3{3}3, 2{3}2{4} p, p{4}2{3}2, 2{4}3{3}3, 3{3}3{4}2 where we have used the generalised schlafli symbols [32] to represent the polyhedra. The polyhedron 3{3}3{3}3 known as the Hessian polyhedron can be thought of as the tetrahedron in the complex space. Its full symmetry group has 648 elements and is called Σ(216 × 3). Like the other discrete groups relevant in flavour symmetry, Σ(216 × 3) is also a subgroup of the continuous group U (3).
The principal series of Σ(216 × 3) [33] is given by We find that, in the context of flavour physics and model building, Σ(72 × 3) has an appealing feature: it is the smallest group containing a complex three-dimensional representation whose tensor product with itself results in a complex six-dimensional representation, 2 i.e.
With a suitably chosen basis for 6 we get where (a 1 , a 2 , a 3 ) T and (b 1 , b 2 , b 3 ) T represent the triplets appearing in the LHS of Eq. (16). All the symmetric components of the tensor product together form the representation 6 and the antisymmetric components form3. For the SU (3) group it is well known that the tensor product of two 3s gives rise to a symmetric 6 and an antisymmetric3. Σ(72×3) being a subgroup of SU (3), of course, has its 6 and3 embedded in the 6 and3 of SU (3).  k  1  1  1  24  9  9  9  18  18  18  18  18  18  18  18  18  or d(C k )  1  3  3  3  2  6  6  4  12  12  4  12  12  4  12  12   1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 1 Consider the complex conjugation of Eq. (16), i.e.3 ⊗ 3 =6 ⊕ 3. Let the right-handed neutrinos form a triplet, ν R = (ν R1 , ν R2 , ν R3 ) T , which transforms as a3. Symmetric (and also Lorentz invariant) combination of two such triplets leads to a conjugate sextet,S ν , which transforms as a6, where ν Ri .ν R j is the Lorentz invariant product of the righthanded neutrino Weyl spinors. We may coupleS ν to a flavon field which transforms as a 6 to construct the invariant term In general, the 3 × 3 Majorana mass matrix is symmetric and has six complex degrees of freedom. Therefore, using Eq. (20), any required mass matrix can be obtained through a suitably chosen vacuum expectation value (VEV) for the flavon field.
To describe the representation theory of Σ(72 × 3) we largely follow Ref. [33]. Σ(72 × 3) can be constructed using four generators, namely C, E, V and X [33]. For the threedimensional representation, we have The characters of the irreducible representations of Σ(72×3) are given in Table 2. Tensor product expansions of various representations relevant to our model are given in Appendix A. There we also provide the corresponding C-G coefficients and the generator matrices.

The model
In this paper we construct our model in the Standard Model framework with the addition of heavy right-handed neutrinos. Through the type I see-saw mechanism, light Majorana neutrinos are produced. The fermion and flavon content of the model, together with the representations to which they belong, are given in Table 3. In addition to Σ(72 × 3), Table 3 The flavour structure of the model. The three families of the left-handed-weak-isospin lepton doublets form the triplet L and the three right-handed heavy neutrinos form the triplet ν R . The flavons φ α , φ β and ξ , are scalar fields and are gauge invariants. On the other hand, they transform non-trivially under the flavour groups we have introduced a flavour group C 4 = {1, −1, i, −i} for obtaining the observed mass hierarchy for the charged leptons. The Standard Model Higgs field is assigned to the trivial (singlet) representation of the flavour groups.
For the charged leptons, we obtain the mass term where H is the Standard Model Higgs, Λ is the cut-off scale, y τ and y μ are the coupling constants for the τ -sector and the μ-sector, respectively.Ā βα is the conjugate triplet obtained from φ β and φ α , constructed in the same way as the second part of Eq. (17), where φ α = (φ α1 , φ α2 , φ α3 ) T and φ β = (φ β1 , φ β2 , φ β3 ) T . L † τ R transforms as 3 × i under the flavour group, Σ(72 × 3) × C 4 . The flavonφ β transforms as3 × −i and hence it couples to L † τ R as shown in Eq. (22). No other coupling involving τ R , μ R or e R with eitherφ β orφ α is allowed, given the C 4 assignments in Table 3. However, L † μ R andĀ βα , which transform as 3 × 1 and3 × 1, respectively, can couple, Eq. (22). Note thatĀ βα is a second-order product of φ β and φ α and it is antisymmetric. No other second-order product transforming as3 exists, since the antisymmetric product of φ β with itself or φ α with itself vanishes. H.T . represents all the higher-order terms, i.e. the terms consisting of higherorder products of the flavons, coupling to e R , μ R and τ R . It can be shown that, for obtaining a flavon term coupling to the e R , we require at least quartic order. 3 The VEV of the Higgs, (0, h o ), breaks the weak gauge symmetry. For the flavonsφ α andφ β , we assign the vacuum alignments 4 3 Refer to Appendix B for an analysis of the higher order products of φ α and φ β . 4 Refer to Appendix C for the details of the flavon potential that leads to these VEVs.
where V is one of the generators of Σ(72 × 3) given in Eq. (21) and is proportional to the 3 × 3 trimaximal matrix. The constant m has dimensions of mass. Substituting these vacuum alignments in Eq. (22) leads to the following charged-lepton mass term: Now, we write the Dirac mass term for the neutrinos: whereH is the conjugate Higgs and y w is the coupling constant. With the help of Eq. (20), we also write the Majorana mass term for the neutrinos: where y m is the coupling constant. Let ξ be the VEV acquired by the sextet flavon ξ , and let ξ be the corresponding 3 × 3 symmetric matrix of the form given in Eq. (20). Combining the mass terms, Eqs. (26) and (27), and using the VEVs of the Higgs and the flavon, we obtain the Dirac-Majorana mass matrix: The 6 × 6 mass matrix M, forms the coupling where ν L = (ν e , ν μ , ν τ ) T are the left-handed neutrino flavour eigenstates.
Since y w h o is at the electroweak scale and y m ξ is at the high energy flavon scale (> 10 10 GeV), small neutrino masses are generated through the see-saw mechanism. The resulting effective see-saw mass matrix is of the form From Eq. (30), it is clear that the see-saw mechanism makes the light neutrino masses inversely proportional to the eigenvalues of the matrix ξ . We now proceed to construct the four cases of the mass matrices, Eqs. (11), (12), all of which result in the neutrino mass ratios, Eq. (14). To achieve this we choose suitable vacuum alignments 6 for the sextet flavon ξ .
Here we assign the vacuum alignment Using the symmetric matrix form of the sextet given in Eq. (20), we obtain Diagonalising the corresponding effective see-saw mass matrix M ss , Eq. (30), we get , 1, leading to the neutrino mass ratios, Eq. (14). The unitary matrix U ν is given by The product of the contribution from the charged-lepton sector i.e. V from Eqs. (25), (21) and the contribution from the neutrino sector i.e. U ν from Eqs. (33), (34) results in the Tχ M (χ =+ π 16 ) mixing: 6 Refer to Appendix C for the details of the flavon potentials that lead to these VEVs.

Tχ M
Here we assign the vacuum alignment resulting in the symmetric matrix In this case, the diagonalising matrix is and the corresponding mixing matrix is with χ = − π 16 .

)
Here we assign the vacuum alignment resulting in the symmetric matrix In this case, the diagonalising matrix is and the corresponding mixing matrix is with φ = + π 16 .

)
Here we assign the vacuum alignment resulting in the symmetric matrix In this case, the diagonalising matrix is and the corresponding mixing matrix is with φ = − π 16 . As stated earlier, the four cases, Eqs.  11)) are composed of real numbers implying they remain invariant under complex conjugation. Therefore, they do not contribute to C P violation. In our model, , the resulting leptonic mixing, V U ν , is also maximally C P-violating (Tχ M). Note that U Tχ M , Eq. (7), is symmetric under the conjugation and the exchange of μ and τ rows. This generalised C P symmetry under the combined operations of μ-τ exchange and complex conjugation is referred to as μ-τ reflection symmetry in previous publications [5,16,[41][42][43]. The conjugation symmetry in the neutrino VEVs together with maximal C P violation from the charged-lepton sector produces the μ-τ reflection symmetry of U PMNS .
Consider the exchange of the first and the third rows as well as the columns of the mass matrix, Eq. (20). This is equivalent to the exchange of the first and the third elements and the fourth and the sixth elements of the sextet flavon, Eq. (19). In Σ(72 × 3), this exchange can be realised using the group transformation by the unitary matrix E.V.V , with E and V given in Eqs. (21), (67). By the group transformation we imply left and right multiplication of the mass matrix using the 3 × 3 unitary matrix and its transpose or equivalently left multiplication of the sextet flavon using the 6 × 6 unitary matrix. The mass matrices, Eqs. (12), and the corresponding flavon VEVs, Eqs. (40), (44), are invariant under the transformation by E.V.V together with the conjugation. The VEVs break Σ(72×3) almost completely except for E.V.V with conjugation which remains as their residual symmetry. 7 The resulting mixing matrix, U PMNS = V U ν , is TφM, which is real and C P conserving. E.V.V -conjugation symmetry in the neutrino VEVs together with maximal C P violation from the charged-lepton sector produces the C P symmetry of U PMNS . Both Tχ M and TφM have a trimaximal second column. This feature of the mixing matrix was linked to the "magic" symmetry of the mass matrix [42,[44][45][46]. In our model, the charged-leptonic contribution, V , is trimaximal. Because of the vanishing of the fourth and the sixth elements of the sextet VEVs, Eqs. (31), (36), (40), (44) which correspond to the offdiagonal zeros present in the mass matrices, Eqs. (11), (12), the trimaximality of V carries over to U PMNS = V U ν .
Consider the unitary matrix, Group transformation by A is the multiplication of the offdiagonal Majorana matrix elements by −1. Invariance under A, implies these elements vanish and ensures trimaximality.
In the literature, small groups like the Klein group [16,[47][48][49][50][51] are often used to implement symmetries like the generalised C P and the trimaximality as the residual symmetries of the mass matrix. However, A is not a group member of Σ(72 × 3). In our model, the vanishing mass matrix elements arise as a consequence of the specific choice of the flavon potential, Eq. (142), rather than the result of a residual symmetry under Σ(72 × 3). The presence of a simple set of numbers in the VEVs (and the mass matrices) is suggestive of additional symmetry transformations (like the one generated by A, Eq. (49)) which are not a part of Σ(72 × 3). The present model only serves as a template for constructing any fully constrained Majorana mass matrix using Σ(72 × 3). We impose additional symmetries on the mass matrix by using flavon potentials with a carefully chosen set of parameters, Table 8. Realising these symmetries naturally by incorporating more group transformations along with Σ(72 × 3) in an expanded flavour group requires further investigation.

Predicted observables
For comparing our model with the neutrino oscillation experimental data, we use the global analysis done by the NuFIT group and their latest results reproduced in Eqs. (2)-(6). They are a leading group doing a comprehensive statistical data analysis based on essentially all currently available neutrino oscillation experiments. Their results are updated regularly and published on the NuFIT website [2]. The value sin 2 θ 13 = 2 3 sin 2 π 16 = 0.02537, 8 is slightly more than the upper limit of the 3σ range, 0.02392. We provide a solution to this discrepancy in the following discussion.
In our previous analysis in Sects. 3.1-3.4, we used the relation U PMNS = V U ν where V is the left-diagonalising matrix for the charged-lepton mass matrix, Eq. (25). However, the diagonalisation achieved by V is only an approximation. In Eq. (25), the presence of the O( 4 ) element in the e L -μ R offdiagonal position in relation to the O( 4 ) electron mass and O( 2 ) muon mass produces an O( 2 ) correction to the diagonalisation, i.e. a more accurate left-diagonalisation matrix is The resulting correction in the e3 element of U PMNS is Since (U PMNS ) e3 = sin θ 13 e −iδ and ≈ m μ The above correction is sufficient to reduce 9 our prediction for sin 2 θ 13 to within the 3σ range. For the solar angle, using the formula given in Table 1, we get This is within 3σ errors of the experimental values, although there is a small tension towards the upper limit. For the atmospheric angle, Tχ M predicts maximal mixing: which is also within 3σ errors. The NuFIT data as well as other global fits [52,53] are showing a preference for nonmaximal atmospheric mixing. As a result there has been a lot of interest in the problem of octant degeneracy of θ 23 [54][55][56][57][58][59][60][61]. TφM predicts this non-maximal scenario of atmospheric mixing. TφM (φ= π 16 ) and TφM (φ=− π 16 ) correspond to the first and the second octant solutions, respectively. Using the formula for θ 23 given in Table 1, we get TφM (φ=+ π 16 ) : TφM (φ=− π 16 ) : The Dirac C P phase, δ, has not been measured yet. The discovery that the reactor mixing angle is not very small has raised the possibility of a relatively earlier measurement of δ [62][63][64]. Tχ M having δ = ± π 2 should lead to large observable C P-violating effects. Substituting χ = ± π 16 in Eq. (9), our model gives which is about 40% of the maximum value of the theoretical range, − 1 . On the other hand, TφM, with δ = 0, π and J = 0, is C P conserving.
The neutrino mixing angles are fully determined by the model, Eqs. (10), (53), (54), (55), (56). Hence, we simply compared the individual mixing angles with the experimental data in the earlier part of this section. Regarding the neutrino masses, the model predicts their ratios, Eq. (14). To compare this result with the experimental data, which gives the masssquared differences, Eqs. (5), (6), we utilise a χ 2 analysis, We report that the predicted neutrino mass ratios are consistent with the experimental mass-squared differences. Using the χ 2 analysis we obtain, are the mass ratios The best fit values correspond to χ 2 min = 0.03 and the error ranges correspond to Δχ 2 = 1, where Δχ 2 = χ 2 − χ 2 min . The results from our analysis are also shown in Fig. 1.
Note that the mass ratios Eq. (14), are incompatible with the inverted mass hierarchy. Considerable experimental studies are being conducted to determine the mass hierarchy [63,[65][66][67][68][69][70][71] and we may expect a resolution in the nottoo-distant future. Observation of the inverted hierarchy will obviously rule out the model.
Cosmological observations can provide limits on the sum of the neutrino masses. The strongest such limit has been set recently by the data collected using the Planck satellite [72,73]:  (61) which is not far below the current cosmological limit. Improvements in the cosmological bounds from Planck data are expected. Future ground-based CMB polarisation experiments such as Polarbear-2 [74] and Square Kilometer Array-2 [75], could lower the cosmological limit to below 100 meV and could also determine the mass hierarchy. Such results may support or rule out our model. Neutrinoless double beta decay experiments seek to determine the nature of the neutrinos as Majorana or not. These experiments have so far set limits on the effective electron neutrino mass [76] |m ββ |, where with U representing U PMNS . In all the four mixing scenarios predicted by the model, Eqs. (35), (39), (43), (47), we have and |U e3 | = 2 3 sin π 16 . Also, all of them result in the phases: 10 Therefore the model predicts Substituting the neutrino masses from Eqs. (59) in Eq. (64) we get m ββ = 23.47 +0.16 −0.14 meV.

Renormalisation effects on the observables
The see-saw mechanism requires the existence of a heavy Majorana mass term coupling the right-handed neuntrinos together. Our model, combined with the observed neutrino mass-squared differences, predicts that the neutrino masses are a few tens of meVs. This places the see-saw scale (also the flavon scale) at around 10 12 GeV. As such, this is the scale at which the fully constrained mass matrices, as proposed in our model, are generated. In order to accurately compare the model with the observed masses and mixing parameters, it is necessary to calculate its renormalisation group (RG) evolution from the high energy scale down to the electroweak scale.
We use the Mathematica package, REAP [87], to numerically study the RG evolution of the masses and the mixing observables. The Mathematica code for calculating the RG evolution relevant to the model is given below: In the above code, the initial values of the mixing observables and the masses are set at 10 12 GeV. The mixing observables are chosen such that they correspond to Tχ M (χ =+ π 16 ) . We set the masses to be 30.55, 32.33 and 69.24 meV. These specific values are chosen such that they are consistent with Eq. (14) (at 10 12 GeV) and give the best fit to the observed mass-squared differences when renormalised to the electroweak scale (100 GeV). MSNParameters in the code gives the renormalised parameters at 100 GeV as its output and thus we get θ 12 = 33.78 • , θ 23 = 45.00 • , θ 13 = 9.165 • , δ = 90.00 • , m 1 = 24.63 meV, m 2 = 26.07 meV, m 3 = 56.16 meV 11 . From these values we conclude that, under the conditions of our model, renormalisation has virtually no effect on the mixing parameters. On the other hand, it affects our predictions for the masses, Eqs. (59), (61), (65), by a few percentage points.
Analysis of RG equations [87][88][89][90][91][92] show that, even though the neutrino masses (the fermion masses in general) evolve appreciably, their ratios evolve slowly. This behaviour is sometimes referred to as "universal scaling". For our model, the light neutrino masses evolve by around 20%, while their ratios by less than 1%. This ensures that the mass ratios, Eq. (14), theorised at the high energy scale remain practically valid at the electroweak scale as well.

Summary
In this paper we utilise the group Σ(72×3) to construct fully constrained Majorana mass matrices for the neutrinos. These mass matrices reproduce the results obtained in Ref. [22] i.e. Tχ M (χ =± π 16 ) and TφM (φ=± π 16 ) mixings along with the neutrino mass ratios, Eq. (14). The mixing observables as well as the neutrino mass ratios are shown to be consistent with the experimental data. Tχ M (χ =± π 16 ) and TφM (φ=± π 16 ) predict the Dirac C P-violating effect to be maximal (at fixed θ 13 ) and null, respectively. Using the neutrino mass ratios in conjunction with the experimentally observed neutrino mass-squared differences, we calculate the individual neutrino masses. We note that our predicted neutrino mass ratios are incompatible with the inverted mass hierarchy. We also predict the effective electron neutrino mass for the neutrinoless double beta decay, |m ββ |. We briefly discuss the current status and future prospects of determining experimentally the neutrino observables leading to the confirmation or the falsification of our model. In the context of model building, we carry out an in-depth analysis of the representations of Σ(72 × 3) and develop the necessary groundwork to construct the flavon potentials satisfying the Σ(72 × 3) flavour symmetry. In the charged-lepton sector, we use two triplet flavons with a suitably chosen set of VEVs which provide a 3 × 3 trimaximal contribution towards the PMNS mixing matrix. It also explains the hierarchical structure of the charged-lepton masses. In the neutrino sector, we discuss four cases of Majorana mass matrices. The Σ(72 × 3) sextet acts as the most general placeholder for a fully constrained Majorana mass matrix. The intended mass matrices are obtained by assigning appropriate VEVs to the sextet flavon. It should be noted that we need additional symmetries to 'explain' any specific texture in the mass matrix. This work was supported by the UK Science and Technology Facilities Council (STFC). Two of us (RK and PFH) acknowledge the hospitality of the Centre for Fundamental Physics (CfFP) at the Rutherford Appleton Laboratory. RK acknowledges the support from the University of Warwick. RK thanks the management of the School of the Good Shep-herd, Thiruvananthapuram, for providing a convenient and flexible working arrangement conducive to research.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix A: Irreps of Σ(72 × 3) and their tensor product expansions
The generator matrices for the triplet representation are provided in Eq. (21). We define the basis for the sextet representation using Eqs. (17). The resulting generator matrices are With (a 1 , a 2 , a 3 ) T and (b 1 ,b 2 ,b 3 ) T transforming as 3 and 3, the tensor product expansion, Eq. (68), is given by In this basis, the generator matrices of the octet representation are The octet is a real representation.
We define the basis for the doublet representation in such a way that 6 is simply the Kronecker product of 2 and3, i.e.
where (a 1 , a 2 ) T and (b 1 ,b 2 ,b 3 ) T represent 2 and3, respectively. In such a basis, the generator matrices for the doublet are The singlets 1 ( p,q) transform as In terms of the tensor product expansion, Eq. (74), these singlets are given by 1 ≡ a T u b, where a and b represent the doublets in Eq. (74) and u, u 1 , u ω and uω are unitary matrices, with σ 2 being the second Pauli matrix and V , X being the generators of the doublet representation, Eq (73).
The C-G coefficients for the above tensor product expansion are given by where (a 1 , a 2 , a 3 , a 4 , a 5 , a 6 ) T and (b 1 , b 2 , b 3 ) T represent the sextet and the triplet appearing in the LHS of Eq. (78).

Appendix B: Hierarchical structure of the charged-lepton mass matrix
The triplet flavons, φ α and φ β , transform as 3 × −i and Table 3. In this appendix, we analyse the cubic and the quartic products which give rise O( 3 ) and O( 4 ) mass matrix elements, respectively in Eq. (25). We neglect the products beyond quartic order.
Cubic products The above expansions, Eqs. (101)-(103), do not contribute to any coupling, since3 does not appear in their RHS.
This tensor product corresponds to the conjugation of Eq. (107). The conjugate expansion will have four3s 13 in the RHS. The product3 ⊗3 ⊗3 ⊗3 can be obtained in terms of the flavons φ α and φ β in several different ways. All these are listed in Table 5. For each combination of flavons, we provide the corresponding C 4 representation. Under φ α ∝ (1, 0, 0) and φ β ∝ (0, 0, 1), we calculate the vacuum alignments of the above-mentioned four3s and list them in the table. ( This tensor product corresponds to the conjugation of Eq. (109). The conjugate expansion will have four3s 14 in the RHS. All the products of φ α and φ β in the form of 3 ⊗ 3 ⊗3 ⊗ 3 are listed in Table 6, along with their respective C 4 representations. Under φ α ∝ (1, 0, 0) and φ β ∝ Table 6 Quartic products of φ α and φ β of the form 3 ⊗ 3 ⊗3 ⊗ 3 leading to3s

Appendix C: Flavon potentials
Here we discuss the flavon potentials that lead to the vacuum alignments assumed in our model. It should be noted that even though our construction results in the required VEVs, we are not doing an exhaustive analysis of the most general flavon potentials involving all the possible invariant terms. However, the content we include is sufficient to realise our VEVs.

The sextet flavon: ξ
We studied the invariants that can be constructed using the sextet ξ up to the quartic order (renormalisable) and found that these terms are insufficient to obtain a potential devoid of continuous symmetries (SU (3) and its continuous subgroups). Therefore, as in the previous subsection, we introduce extra flavons to break the continuous symmetries and to ensure that the potential has a discrete set of minima. The extra flavons introduced here are a doublet η and two triplets φ a , φ b . The flavons used in the charged-lepton sector (φ α , φ β , η α , η β ) are kept distinct from the flavons used in the neutrino sector (ξ , η, φ a , φ b ) in order to avoid unwanted couplings between the two sectors. Table 7 provides the complete list of flavons in the model along with the fermions.
Our first step is to write the potential terms for η, φ a and φ b , similar to Eq. (126), The individual invariant terms in Eq. (127) are T η = (|η| 2 − m 2 ) 2 + Re 2 (η T u 1 η) + Re 2 (ω η T u ω η) + Re 2 (ω η T uω η), where S a , K a and S b , K b are defined similar to S α , K α in Eqs. (112), (120) having φ α , η α replaced with φ a , η and φ b , η, respectively. As described earlier, it is straightforward to show that each term in Eqs. (128)-(133) vanishes, if we assign the following VEVs: S a and S b are the sextets constructed from φ a and φ b , respectively. We may also construct a sextet combining φ a and φ b together, Using the VEVs, Eqs.
If the second column of U PM N S is trimaximally mixed, then the VEV, ξ , as well as the resulting X have non-zero elements only in the first, second, third and the fifth positions. As shown in Eqs. (138), (139), (140), S a , S b and S ab have non-zero elements only in the first, third and the fifth position, respectively. Therefore, a linear combination of ξ , X ξ , S a , S b and S ab can be constructed which fully vanishes. With this information in hand, we construct the potential term, T ξ = mξ +c 1X +c 2Sa +c 3Sb +c 4Sab T × (m ξ + c 1 X + c 2 S a + c 3 S b + c 4 S ab ) , where c 1 , c 2 , c 3 and c 4 are constants. T ξ couples the sextet flavon, ξ with the triplet flavons, φ a and φ b . Any neutrino mass matrix which leads to a trimaximally mixed column can be obtained using a potential of the form, Eq. (142). The values of the constants resulting in the four VEVs, Eqs. (31), (36), (40), (44), are given in Table 8.
Using an appropriate choice of the constants, c 1 , c 2 , c 3 , and c 4 , we may obtain any mixing scheme within the constraint of a trimaximal column. It can be shown that, having the symmetry of c 2 and c 3 being real (invariant under complex conjugation) leads to Tχ M. In the case of TφM, the symmetry is the simultaneous conjugation and interchange of c 2 and c 3 . Additionally, the fact that c 1 , c 2 , c 3 and c 4 are related by simple ratios points to the presence of more symmetries, the study of which is beyond the scope of this paper.
With the help of first-and second-order partial derivatives of a given potential, its minima can be calculated, as was done in previous work, e.g. in Ref. [93]. Using such a procedure, Table 8 The values of constants appearing in the potential, Eq. (142), for the sextet flavon, ξ , corresponding to the four cases. We have t = tan( π along with numerical analysis, we have verified that every potential discussed here has a discrete set of minima and that the quoted VEVs are included among those minima in each case.