Revisiting the Universal Texture Zero of Flavour: a Markov Chain Monte Carlo Analysis

We revisit the phenomenological predictions of the Universal Texture Zero (UTZ) model of flavour originally presented in arXiv:1710.01741, and update them in light of both improved experimental constraints and numerical analysis techniques. In particular, we have developed an in-house Markov Chain Monte Carlo (MCMC) algorithm to exhaustively explore the UTZ's viable parameter space, considering both leading- and next-to-leading contributions in the model's effective operator product expansion. We also extract -- for the first time -- reliable UTZ predictions for the (poorly constrained) leptonic CP-violating phases, and ratio observables that characterize neutrino masses probed by (e.g.) oscillation, $\beta$-decay, and cosmological processes. We therefore dramatically improve on the proof-in-principle phenomenological analysis originally presented in arXiv:1710.01741, and ultimately show that the UTZ remains a minimal, viable, and appealing theory of flavour. Our results also further demonstrate the potential of robustly examining multi-parameter flavour models with MCMC routines.


Introduction
The bulk of the free, unexplained parameters in the Standard Model (SM) of particle physics originate in its flavour sector, thanks to the replication of SM fermion generations with distinct masses and quantum mixings.These parameters are technically natural, in that sending them to zero recovers a global U(3) 5 flavour symmetry of the Lagrangian [2,3].However, the Yukawa couplings of SM fermions to the Higgs boson break this symmetry in a deeply flavournon-universal manner, with a mass ratio of ∼ O (10 12 ) between (e.g.) neutrinos and the top quark.Furthermore the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix exhibits a hierarchical, approximately unit structure [4], while the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix is extremely non-hierarchical, with large mixings amongst generations [5].These highly disparate patterns of fermionic mass and mixing strongly hint that the origins of flavour in the SM may be dynamical, as opposed to a random, soft deviation from an accidental symmetric limit.
The flavour puzzle therefore remains a compelling motivation to search for physics Beyondthe-SM (BSM), as it can be solved dynamically via the breakdown of an ultraviolet (UV), BSM symmetry in specific directions of flavour space.This symmetry breaking typically occurs when exotic scalar familons develop special vacuum expectation values (vev) as determined by a family-symmetric scalar potential, although other alignment mechanisms are conceivable.When familons couple to SM fermions and the Higgs boson, their flavoured vevs shape the otherwise free Yukawa matrices of the SM, and therefore control their associated mass eigenvalues and mixing angles after electroweak symmetry breaking.These predictions can be compared to global flavour data sets to falsify the model, serving as an indirect probe of the new physics proposed.
While the predictions of flavour models -derived from either top-down or bottom-up considerations -are rich, they are also becoming increasingly difficult to falsify, given that experiment is rapidly resolving all SM flavour parameters to a high degree of precision, such that the models' predictions should actually be considered postdictions.Indeed, virtually all quark masses and CKM mixings are measured with exceptional accuracy, while only the PMNS angle θ l 23 ,1 the Dirac CP-violating phase δ l , absolute neutrino mass eigenvalues, and Majorana CP-violating phases (if relevant) are poorly constrained in the leptonic sector.While physical observables that depend on non-trivial combinations of these parameters, e.g.neutrinolessdouble-β decay rates (0νββ), single β-decay rates, and the sum of neutrino mass eigenvalues (as constrained by cosmology), offer additional independent probes of flavour models, it is conceivable that a believable BSM theory will also make falsifiable predictions for a subset of the aforementioned, unresolved constituent flavour parameters.Complicating matters further, many (most) BSM flavour models introduce a number of UV theory parameters that are difficult to numerically sample in a fully generic manner, and so extracting concrete predictions from said models is challenging in its own right.
In light of this experimental situation, and in response to the need for more robust analysis routines for exploring viable model parameter spaces, we will re-examine the Universal Texture Zero (UTZ) Model originally presented in [1].The UTZ is an effective theory (EFT) valid at mass scales above those characteristic of the SM, but below those of hypothetical (and potentially unfalsifiable), renormalizable UV completions, e.g.those incorporating ultra-heavy fermionic messenger fields V: Λ V > Λ UTZ > Λ SM .Its Yukawa sector is therefore generated only at the non-renormalizable level, with EFT expansion parameters in inverse powers of the messenger masses M i .The UTZ Lagrangian is symmetric under a ∆( 27) ≃ (Z 3 × Z 3 ) ⋊ Z 3 [6][7][8][9][10] non-Abelian discrete family symmetry and a further Z N discrete shaping symmetry, and is consistent with an underlying stage of SO (10) grand unification as all fermions and their conjugates -including right-handed (RH) gauge-singlet neutrinos -are assigned as triplets 3 under ∆ (27).Critically, the additional scalars introduced are charged such that a ∆ (27)-invariant scalar potential exists that drives family-symmetry breaking as mentioned above, yielding symmetric mass matrices with a characteristic texture zero in the (1,1) position for all family sectors.As shown in [1], this UTZ structure is capable of explaining quark and lepton flavour data with as few as nine infrared (IR) theory parameters, and therefore amounts to an appealing and predictive theory for the origin of SM flavour patterns.The UTZ stands as a continuation of similar solutions employing texture zeroes, explored already in e.g.[11,12].
However, the numerical exploration of the UTZ parameter space presented in [1] only achieved a 'proof-in-principle' fit demonstrating the model's phenomenological viability.It did not exhaustively explore the predictions of the UTZ Lagrangian at leading order (LO) in its EFT expansion parameters, nor did it consider the complete set of corrections generated by operators present at next-to-leading order (NLO) in 1/M i .Most importantly, the analysis in [1] did not present robust predictions for the aforementioned unresolved leptonic flavour parameters nor any other observables (e.g.β-decay rates) that depend on them, and hence it did not provide a reliable means of falsifying the UTZ model space as data continues to Fields ψ q,e,ν ψ c q,e,ν H Σ S θ 3 θ 23 θ 123 θ θ X ∆(27) 3 3 1 00 1 00 1 00 3 3 3 3 3 Table 1: The fields and ∆(27) × Z N family symmetry content of the UTZ flavour model.Note that θ X only appears in the scalar potential, and hence the only restriction on its Z N charge is that it does not contribute significantly to the fermionic mass matrices.Its Z N charge can therefore be left generic, as shown.
improve.In this paper we aim to remedy these shortcomings by applying a Markov Chain Monte Carlo (MCMC) fitting algorithm to the UTZ.Inspired by similar analyses [13,14] (also see [15][16][17] for advanced statistical analyses of flavour-texture models), the MCMC technology we employ allows for a robust exploration of multi-parameter models and their associated likelihoods.It also allows one to simultaneously extract predictions for poorly-constrained observables which are controlled (in part and in different combinations) by the same parameters that control exceptionally well-constrained observables, thereby accounting for the intricate correlations between UTZ theory parameters and their associated phenomenology.In this way we are capable of presenting predictions in experimentally-preferred regions of the UTZ parameter space, for both the LO and NLO UTZ Lagrangian.As we will show, the theory is phenomenologically viable at both orders in its operator product expansion, with the latter NLO terms yielding only minor corrections to the dominant LO predictions.The UTZ is therefore a stable, predictive, and minimal theory of flavour.The remainder of the paper develops as follows: in Section 2 we review the UTZ model as conceived in [1], including the field and symmetry content composing the (N)LO contributions to its operator product expansion, as well as its qualitative predictions in the quark and lepton sectors.Then in Section 3 we review the most up-to-date experiment that constrains its predictions in the Yukawa sector, and also discuss the uncertainties associated to renormalization group evolution (RGE) from the UV to the IR.In Section 4 we discuss the MCMC algorithm we have developed to explore the UTZ parameter space, and also present the results and analysis following from our scans.We conclude in Section 5.

The Universal Texture Zero Model
The field and ∆(27) × Z N symmetry content of the UTZ [1] is given in Table 1.There one observes that all SM fermions ψ a are assigned as triplets 3 under the family symmetry, as are additional gauge singlet 'sterile' neutrinos that participate in a seesaw mechansim.Besides the fermionic content, we also have a set of BSM scalar familons θ i , charged as ∆ (27) antitriplets 3, a lepton-number-violating (LNV) anti-triplet familon θ necessary for describing the Majorana neutrino mass sector, and finally a triplet familon θ X necessary for successful vacuum alignment.All such familons are SM gauge singlets.There is also a ∆ (27) singlet sector composed of the Σ and Higgs H scalars, both associated to an underlying stage of grand unification consistent with the following symmetry-breaking chain: where the SO (10) breaking proceeds via an H vev and where The first term V 1 sets the scale of the scalar familon fields, and is sufficient to break the family symmetry spontaneously upon m 2 i being driven to negative values, perhaps via radiative corrections in the manner of [26].Then the second term V 2 aligns the θ 3,123 vevs in flavour space as a function of the sign of h i ; h 123 ≡ h i > 0 while h 3 ≡ h i < 0. The terms V 3,4,5,6 account for the final alignment of the θ 23,X and θ vevs, with V 3 sourcing the dominant coupling of θ X , V 4 selecting ⟨θ X ⟩ ∝ (2, −1, 1) out of the two degenerate vacua V 3 allows, and V 5 and V 6 respectively driving the final θ 23 and θ orientations upon minimization. 2 All subtleties considered, the potential in (3) aligns the scalar familon fields in special directions in flavour-space, Observe that in (2) we have only included terms that are consistent with a spontaneously broken, supersymmetric (SUSY) underlying theory with triplet mediators.Additional quartic terms may appear, but must be suppressed in order to preserve (4).All other aspects of the tree-level phenomenology of the UTZ model can be studied without reference to hypothetical UV completions, and we adopt this agnosticism to be as generic as possible in what follows.
where the parentheses on the first term indicate that both θ 3 and θ are aligned in the thirdfamily direction, and where we have included the generic phases α, β for completeness, although we will eventually set these to zero following the discussion in [27].We note that, of the above potential terms, many are not invariant under SU(3) F , and so the use of ∆( 27), a non-Abelian discrete subgroup of SU(3) F , was instrumental in the above discussion.

The Leading-Order Effective Yukawa Lagrangian
Upon demonstrating that a successful vacuum alignment is plausible upon family-symmetry breaking, a meaningful BSM Yukawa sector can be subsequently formed from the field and symmetry content of Table 1.This leads to the following LO UTZ effective Lagrangian in the Dirac sector of the theory: are free Wilson coefficients whose superscript denotes the mass dimension n of the operator, while M i,f represent the mass scales associated to heavy messenger fields that have been integrated out of the spectrum in forming the EFT, a lá the Froggatt-Nielsen mechanism [28].These messenger fields are associated to distinct UV completions and are typically taken to be vector-like fermions, although we do not wish to commit ourselves to any particular scenario.In what follows we will simply point out the implications and constraints on said UV messengers coming from the (falsifiable) IR spectrum associated to (4).
To that end, one quickly notices that a natural hierarchy for the third-family fermions is realized, thanks to the power suppression (assuming only mild hierarchies amongst messenger masses) of the second and third terms with respect to the first, which only contributes to the (3,3) entry of the Dirac mass matrices.While this helps realize an approximate SU(2) F symmetry of the quark mass matrices and associated CKM mixing matrix, it also implies that the ratio θ 3 /M 3,f is large [29], at least in the up sector.This is acceptable if θ 3 is the dominant contributor to the messenger mass, which we assume for all charged fermion sectors.For an alternative solution to this issue involving Higgs mediators, see [30].
Besides (4), the field and symmetry content of Table 1 also permits a Majorana mass Lagrangian, which at leading order in the OPE is of the following form: Here one notices that there are two insertions of the LNV scalar θ in each operator, as is consistent with our underlying SO(10) → SU(4)×SU(2) L ×SU(2) R GUT embedding, and also that the leading contributions in this effective Lagrangian are at dimension five and eight in the 1/M expansion of the EFT, as opposed to six and seven in the case of the Dirac Lagrangian given in (4).This results in an extremely dominant third-family hierarchy that has important phenomenological implications in the neutrino sector upon applying the seesaw, as mentioned below.Further discussion regarding the relative power suppression between Dirac and Majorana sectors will be given in Section 2.2.

Qualitative Charged Fermion Masses and Sum Rules
While the SM's quark and charged lepton flavour sector is exceptionally well-measured and therefore offers little opportunity for novel predictions, we do note that the UTZ Lagrangian above has been designed to realize successful charged fermion mass ratios, as well as two longstanding and successful phenomenological ansätze: the Georgi-Jarlskog mechanism and the Gatto-Sartori-Tonin sum rule.This is due to the UTZ structure of the Dirac mass matrices, given qualitatively by with f again indicating the family sector, f ∈ {u, d, e} and ϵ f associated small parameters.
Phenomenologically viable values are given by ϵ u ≈ 0.05 and ϵ d,e ≈ 0.15.This family splitting is accommodated via the UTZ relations which hold up to O(1) coefficients and signs.Here one sees that , and of course ϵ e and ϵ d can be equal given the symmetry breaking in (1) and the fact that both are T R,3 = −1/2 states which (in SUSY models) acquire their mass from the same Higgs boson (H d ).Note that we assume the messenger masses carry both lepton and quark quantum numbers, and so it is important that RH messengers associated to SU(2) R breaking in (1) dominate for ϵ u < ϵ d , as opposed to the LH messengers associated to SU(2) L , whose up and down masses are of course equal due to SU(2) L invariance.Also associated to this pattern of family suppression are the r f coefficients in (6), which are sourced from the Σ vev and therefore implement the Georgi-Jarlskog mechanism [31], resulting in at the scale of grand unification, as is consistent with RGE and threshold corrections [12].In addition to these characteristic mass relations, the quark sector realizations of (6) also implement the Gatto-Sartori-Tonin mixing sum rule [32], relating the Cabibbo angle θ c ≃ θ q 12 to mass ratios from the first-and second-generation quarks of both the up and down families.Setting δ ≈ π/2 and again accounting for RGE and threshold corrections [12] (cf.Section 3.1), both (8) and (9) remain successful predictions that are maintained in our UTZ framework. 3

Qualitative Neutrino Masses and Sum Rules
As discussed above, the bulk of the parameters left to be constrained experimentally are in the neutrino sector, and so it is worthwhile to discuss the qualitative, analytic predictions of the UTZ construction in this area.As with other flavour models, we employ the Type-I seesaw mechanism [34][35][36][37].In this framework the right-handed (RH) Majorana mass terms generated by (5) are parametrically heavier than the Dirac neutrino masses coming from (4).Integrating the RH neutrino fields out of the spectrum generates a left-handed (LH) Majorana neutrino mass term, where the phase η is defined from the predicted ratio of the heavy neutrino masses m 2 /m 3 (see the discussion in [1]).One notes that the relationships in (12) are similar to the renowned 'Tri-Bimaximal' (TBM) texture [42] (sin θ T BM 13 = 0, sin θ T BM ) that is often a starting point for neutrino mass model building.However, the salient difference with respect to prior models of this type -see e.g.[29,33] or the more recent [43][44][45] -is that the (1,1) texture zero of the mass matrix remains after application of the seesaw, such that our UTZ setup leads to a non-negligible departure from the TBM texture and naturally allows for a large(r) reactor mixing angle θ l 13 (which also receives corrections from the chargedlepton sector), in accord with data.Finally, we note that our ∆(27) family-symmetry breaking realizes the Z 2 × Z 2 [46] residual symmetry of the IR neutrino mass term only 'indirectly,'4 in that it is not a subgroup of ∆( 27) and appears only accidentally thanks to (4).

Higher-Order Contributions
The Lagrangians in (4)-( 5) represent the LO contributions to the UTZ operator product expansion.Higher-order terms in this series are suppressed by further powers of the relevant mediator masses, and should therefore represent small corrections to the qualitative structures and predictions discussed above.However, these corrections can a priori be non-negligible as noted in [1], and so it is important that we consider them robustly as we revisit the UTZ.
In the Dirac sector, the NLO ∆( 27) × Z N invariant terms composed of the same field content as in Table 1 arise at mass-dimension eight, i.e. with four powers of mediator mass suppression, from which once can readily conclude that the HO contributions ∝ S 2 in ( 13) are also parametrically smaller than those ∝ ΣS: In [1] we used (15) to justify ignoring the S 2 contribution to the Dirac mass matrix entirely.However, we will now include both terms in (13) for completeness.The UTZ's operator product expansion is of course infinite-dimensional in the absence of an explicit UV completion.Hence further, next-to-next-to-leading order contributions can also be written down.However, these operators will have at least three additional insertions of ∆( 27) triplets, and are therefore highly suppressed.We neglect their contributions as a result.We also note that the NLO contributions to the Dirac Lagrangian given in (13) enter at the same mass-dimension as LO contributions to the Majorana Lagrangian given in (5), and so we do not consider any corrections to (5) to be consistent in our power counting.

Complete Effective Mass Matrices in the Ultraviolet
The discussions in the Subsections above lead to the LO and NLO Lagrangians of ( 4), ( 5) and ( 13).After family-and electroweak-symmetry breaking, these Lagrangians generate the following Dirac and Majorana fermion UTZ mass matrices: y e i(α+β+ρ) (x e −iρ + 2y e −iϕ ) e i(2α+ρ+ϕ) x e i(α+ϕ) Here the matrices have been normalized such that where M θ and s are the overall scale-setting parameters of ( 17)-( 18) which, along with the relative-scalesetting parameters {a, b, c, d, x, y}, are defined in terms of scalar vevs and other coefficients: with r u,d,e,ν = (1, 1, −3, −3)/3 and ⟨Θ 23,123 ⟩ ≡ ⟨θ k θ k θ k 23,123 ⟩, i.e. the vev of the singlet contractions with k superscript in (5).The relationship between primed and unprimed parameters, along with associated complex phases, is then given by and it is clear that c and d account for the HO Dirac corrections in (13).Given ( 17)-( 18), the values of the 'physical' fermionic mass, mixing, and CP-violating parameters can be extracted numerically as described in [1] or analytically, using flavourinvariant theory as described in (e.g.) [49][50][51].Then, given ( 19)- (20), one can compare the number of IR theory parameters vs. IR physical parameters as a measure of the predictivity of the UTZ.At LO, there are a priori two coefficients (a,b) and two phases (γ, δ) for each charged fermion sector, plus the additional two family-universal phases (α, β) from vacuum alignment.However, following [27], we can set all but two of these phases to zero without loss of generality.Assuming the GUT embedding discussed above to relate the down quarks to charged leptons taking into account the Georgi-Jarlskog factors, one then has (2 + 2) • 3 + 2 − 4 − 4 = 6 UTZ model parameters (including two phases) to describe three CKM mixing angles, one CKM Dirac phase, four quark mass ratios and two charged lepton mass ratios, totalling 10 physical parameters.The neutrino sector's predictivity is even more striking, in the sequentially dominant limit of (12).There, only three parameters, including a phase and an overall mass scale, are necessary to reproduce the neutrino mass differences, which when combined with the aforementioned charged lepton parameters also generate PMNS angles and phases.In total, we see that only nine theory parameters are required to reproduce 18 physical parameters at LO in the UTZ OPE.This is to be compared to the SM where, before allowing for weak basis transformations or rephasing freedoms, the same physical parameters are controlled by three, 3 × 3 complex matrices (charged fermion Yukawas) and unspecified neutrino mass operator(s) (taking the Weinberg operator [52] as an infrared limit of the seesaw mechanism, at least one additional complex symmetric 3 × 3 matrix must be introduced).

Experimental Constraints
The core experimental constraints on the UTZ model presented in Section 2 are of course the fermionic mass eigenvalues and CKM/PMNS mixings extracted from a host of low-and high-energy flavour experiments.Regarding the charged fermion sector, this information is regularly collated in the PDG review [4], which reports bounds on fermion masses and mixing angles.We have reported these IR bounds for the mass sector in Table 3, translating the uncertainties on individual masses into uncertainties on mass ratios, given that the UTZ only predicts the charged fermion mass spectrum up to a common scale.On the other hand, uncertainties on mixing angles and the Dirac CP phase can be extracted from global fits to the CKM matrix and Jarlskog invariant given by [4] where the left equality defines the CKM as the overlap of the matrices U u,d diagonalizing the up / down Yukawa couplings.The translation of these bounds to the θ q ij and δ q basis is given in Table 2.
Leptonic mass and mixing constraints are of course deeply sensitive to ongoing neutrino oscillation, cosmology, and β-decay experiments.The authors of [5] have compiled a global fit to the available oscillation data, finding (e.g.) where the LHS again gives the standard definition of the PMNS matrix as it appears in the charged-current interactions in terms of constituent charged-lepton and neutrino mixing matrices U l,ν , and the 3σ confidence bounds on the RHS further assume a unitary V PMNS and include Super-Kamiokande atmospheric data -see [5] for details.As seen in [5] and also in Table 2, current 3σ oscillation constraints do not yet fully determine the quadrant of the atmospheric mixing angle θ l 23 and, at least in the normal ordering scenario, have only excluded ∼ 43% of the available domain of the leptonic CP-violating phase δ l , i.e. δ l is only constrained within a ∼ 200 • arc.This is reduced to an exclusion of only ∼ 20% of the phase domain when not including SK data.
The authors of [5] have also obtained global constraints on the differences of squared neutrino mass eigenvalues, finding at the 3σ confidence level  [4], while the UV bounds are estimated by considering various input RGE/threshold correction parameter choices from [12], and accounting for the propagated IR experimental uncertainties.In the PMNS sector we take 3σ global bounds from NuFit, in the normal ordering scenario and incorporating Super-Kamiokande atmospheric data.
in the normal-mass-ordering scenario relevant to the UTZ construction, and again including Super-Kamiokande atmospheric data. 5We have translated this to a bound on the ratio ξ ≡ ∆m 2 sol /∆m 2 atm in Table 3.A second class of neutrino mass constraints comes directly from cosmological probes.For example, assuming the ΛCDM model and using data from the Cosmic Microwave Background's (CMB) angular spectra, the Planck experiment has put an upper bound on the sum of cosmologically stable neutrino masses m Σ of [54] m Σ ≡ i m ν i < 0.26 eV (24) at the 95 % confidence level.When also including data from Baryon Acoustic Oscillations this bound is reduced to m Σ < 0.12 eV [54], which can be yet further reduced to m Σ < 0.09 eV when including Type Ia supernova luminosity distances and growth rate parameter determinations [57] -see [58] for a recent review.Finally, additional constraints in the neutrino mass and mixing sector originate in the effective mass terms controlling electron-neutrino 0νββ decay and single β decays, where V ei is the matrix element of the first row and i-th column of the PMNS matrix defined in (22), and m i is the corresponding neutrino mass eigenvalue.Robust bounds for these quantities Table 3: The same as Table 2, but for fermion masses.We estimate the neutrino mass squared difference in the UV from [53], and recall that the ξ ratio only differs from the IR when tan β is large and/or the neutrino mass spectrum is partially degenerate (see text), hence the two UV bounds for ξ in the last row, with the left (right) cell corresponding to the low (high) tan β scenario.IR bounds for m β(β) and m Σ are taken from [54][55][56], and their corresponding UV bounds are given by conservatively setting s = 1.4 (see text).
are provided by dedicated terrestrial experiments.In (25) we have cited the KamLAND-Zen collaboration [55], while the limit in (26) is the 90% confidence-level bound from KATRIN [56].KATRIN's future sensitivity is expected to reach m β < 0.2 eV [59]. 6t this point we should clarify that, while in practice we must specify a numerical value for the scale M θ when applying the seesaw formula (10), 7 and must therefore count this as a relevant IR model parameter in our theory (unlike the charged fermion case), this choice completely determines the overall neutrino mass scale.As this is a free UTZ model parameter, we can vary it between sensible UV seesaw scales to accommodate (e.g) ( 23)-( 26), and probe said variation's effect on the relative-scale-setting parameters {a ν , b ν , x, y}.We have done so between M θ ∈ 10 10−12 GeV, and observe in Figure 1 that, at least qualitatively, larger values of M θ are preferred.Regardless, as with the charged fermion spectrum, we only consider UTZ predictions for ratios of observables that depend on the overall neutrino mass scale, where M θ cancels, as truly meaningful.For this reason we will present ξ, ⟨m β(β) ⟩/m Σ , and ⟨m ββ ⟩/m β as predictions in Section 4, but not their individual constituent parameters m ν i or m β(β) , although we can report these based on the various M θ 's identified in the MCMC evolution, and indeed m β(β) , ∆m 2  sol,atm , and m Σ should be seen as giving reliable constraints on the MCMC algorithm, in that once a value for M θ has been settled on, their associated values must still be consistent with observation.

Renormalization Group Evolution Uncertainties
A critical uncertainty for any prediction of the UTZ model comes from the fact that ( 17)-( 18) are the textures associated to the UV theory.Any comparison with data must be made at the scale where said experimental constraints are obtained, which in the case of the UTZ model is orders of magnitude below where ( 17)-( 18) hold.Thankfully the Renormalization Group (RG) evolution required to account for this scale separation is well-studied in the context of a background SM or minimally-SUSY SM (MSSM) spectrum (consistent with our vacuum alignment mechanism discussed above), for both the charged fermion [12,62,63] and neutrino sectors [26,53,[64][65][66][67].For example, assuming a large flavour-breaking scale and a background SUSY spectrum allowing for high-scale gauge-coupling unification, the phenomenological charged-fermion structures discussed between ( 6)-( 9) are already consistent with RGE and threshold corrections from the IR to UV -see e.g.[12] -up to uncertainties regarding (e.g.) the underlying SUSY breaking scale and parameter spaces (in particular the ratio of Higgses, tan β), which can drive some mass and mixing-angle splittings.As we do not specify tan β or other parameters and/or fields beyond those of the Yukawa sector of the EFT presented in Section 2.1, we have used [12] to estimate the overall uncertainty associated to UV quark mass and mixing parameters, accounting for the broad range of possible theory parameters studied therein, and of course propagating updated IR experimental uncertainties from [4] to the UV.These estimates are reported in Table 2-3.
Moving to the neutrino sector and assuming a Type-I seesaw mechanism, detailed RGE and threshold correction analyses can be found in [53,64] and references therein.There one concludes that, in the absence of a conspiracy between special alignments of phases, large tan β, and/or a (partially-)degenerate 8 light neutrino spectrum, radiative corrections to PMNS mixing angles and phases between disparate scales is generally minimal, for the lowest seesaw scale of M 1 , and with √ ξ} for {θ l 12 , θ l 13 , θ l 23 , δ l , ϕ i }, respectively.Taking rough order-of-magnitude estimates for Γ enh and allowing for M 1 as large as M GUT ∼ 10 16 GeV and tan β as large as 50, one sees that typically ∆θ l ij ≲ O(10 −2 ), which is largely insignificant in comparison to the experimental uncertainties on mixing parameters given in Table 2, except for possible corrections to θ l 13 .Given that we predict a normallyordered, hierarchical mass spectrum as a result of the sequential dominance condition of (11),9  we can take the 3σ bounds of (22) to hold in the UV, implying that the experimental uncertainties are large in comparison to radiative effects.
On the other hand, the light neutrino mass eigenvalues are far more sensitive to RGE than are the PMNS parameters, even in a hierarchical system.Assuming small tan β, neutrino mass eigenvalues generally evolve with a common scaling, m ν,i (µ) ≈ s(µ, µ 0 ) m ν,i (µ 0 ) with (e.g.) s ≈ 1.1 − 1.2 for tan β ≈ 10 or s ≈ 1.35 − 1.4 for SM-like running.This obviously leads to a UV enhancement of the neutrino mass differences in (23) ∝ s 2 , but this effect cancels in the ratio ξ.On the other hand, large tan β can drive UV flavour splittings amongst the neutrino mass eigenvalues, evolving both ∆m 2 sol,atm and ξ, an effect which is especially enhanced in the case of a (partially-)degenerate spectrum, and which is considerably uncertain when allowing for generic phase configurations.We have used [53] to estimate the effect on ξ in this regime in Table 2, where one sees that an uncertainty greater than an order of magnitude in principle exists, although this is quite conservative given the neutrino mass domain considered in [53], and the fact that we can constrain our MCMC scan to prefer a hierarchical mass spectrum, i.e. ∆m 2 sol /m 2 ν 1 ≫ 1. 10 Finally, we note that RGE discussed above also impacts the UV values of ( 24)- (26), which serve as constraints on the MCMC system.In Table 3 we have estimated these in the (conservative) SM-like scenario, with s = 1.4 for all neutrino species.
In summary, we will apply the UV bounds in Tables 2-3 to account for a rather generic class of RGE and threshold corrections to fermionic mass and mixing in the UTZ.They will allow us to robustly explore the UTZ's predictions without introducing unnecessary assumptions about the background field content and/or non-flavour parameter spaces that are irrelevant to the EFT construction at hand, which is designed to be as model-independent as possible.

An MCMC Scan of Parameter Space
A proof-in-principle numerical analysis of the UTZ predictions derived from ( 17)-( 18) was originally performed in [1], in order to show that the model was consistent with available mass and mixing data at the time.This semi-analytic study, while successful, relied on a largely heuristic contour analysis to identify a viable region of the UTZ parameter space.However, the analysis was incomplete in many ways, in that it did not 1.exhaustively explore the available UTZ model space, robustly accounting for all theory correlations amongst its Lagrangian parameters and therefore conclusively determine whether the LO UTZ effective Lagrangian adequately describes nature; 2. explore the complete set of corrections coming from NLO effective operators as discussed in Section 2.2.Only the largest corrections identified in the Dirac Lagrangian were briefly considered in [1], and only in the down-quark sector (the corrections parameterized by d d and ψ d ).
3. identify sufficiently generic predictions for (e.g.) the CP-violating phases δ l and ϕ 1,2 or PMNS atmospheric angle θ l 23 , when all other (well-measured) flavour parameters were simultaneously resolved by the UTZ; 4. consider in any way the experimental constraints from, nor predictions for, neutrinosector observables like 0νββ, single β-decay rates, or the sum of neutrino masses m Σ .
Furthermore, the experimental datasets available for theory comparison have been updated since the original publication of [1].All of these considerations motivate us to revisit our phenomenological analysis of the UTZ in order to better determine its viability and identify means of falsifying it.However, given the number of free parameters introduced by ( 4), (5), and even ( 13), numerical techniques more sophisticated than those applied in [1] will be necessary to achieve 1-4.To that end, in this work we consider a Markov Chain Monte Carlo (MCMC) algorithm for exploring the UTZ.

The Generic MCMC Algorithm
Our numerical analysis will rely on a Metropolis-Hasting MCMC algorithm.The purpose of this approach is to find the posterior distribution of the model after applying relevant experimental constraints, thereby obtaining viable, high-likelihood UTZ parameter regions.
The MCMC technique has proven to be very powerful when applied to the exploration of high-dimensional parameter spaces, with physics applications originating in phenomenological studies of SUSY extensions of the SM [68][69][70], cosmology [71], and the determination of parton distribution functions [72].More recently, two publications have used this approach to study the viability of a flavoured SUSY SU(5) model [13] and a scotogenic model for loop-induced neutrino masses [14], from which we will follow most of the methodology.The algorithm is based on an iterative process where every new proposed parameter point is selected in an area near to the previous one, and its estimated viability drives its acceptance in the chain.To be more explicit, every Markov chain starts on a randomly selected point within the parameter interval ranges.Then, on every iteration, a new point with parameters ⃗ θ n+1 is proposed in the vicinity of the previous point with parameters ⃗ θ n .In our study, the new proposed parameter value is computed according to where Π {a, b} is a Gaussian distribution with mean value a and standard deviation b.The parameter κ parametrizes the allowed jump length between two iterations, and its value is chosen empirically in order to maximize the efficiency of the algorithm.If the calculated value exceeds the limits of the defined intervals for the model parameters, the point is rejected.
We then compute the global likelihood associated with the proposal point L n+1 , which is accepted with a probability p = min 1, which enforces the acceptance of points with higher likelihood and conditions the acceptance to the viability of the proposal point with respect to the previously accepted one.Since our objective is to evaluate the global posterior distribution for the parameters, accepting points with lower likelihoods and using a large number of chains is important both in attempting to avoid an MCMC evolution where the chain gets trapped in local maxima that might be fine-tuned, and also in helping to enforce a better distribution across the full parameter space.
For simplicity, we assume here that our experimental constraints are not correlated and that the global likelihood is simply the products of the individual likelihoods, i.e where ⃗ O is the set of experimental observables used as constraints.Furthermore, we assume a Gaussian shape for all the constraint likelihoods where uncertainties are given in Table 2-3, except for constraints that only correspond to upper or lower bounds.In these latter cases we apply a step function whose likelihood is assigned to 1 if the bound is satisfied, and which otherwise employs a Gaussian 'corrective' factor that diminishes the likelihood assigned to the phase-space point as a function of the extent to which the bound is violated.Within this numerical setup, the chain will converge to high-likelihood domains whose area represents the viability of the models according to the applied uncertainties on the constraints.
Additionally, in order to speed up the convergence process, we modify our jump parameter κ to include a memory of proposal tries where t is the number of tries before accepting a new point in the chain.This becomes extremely helpful as the chain converges since some parameters might have a very thin datacompatible range.As soon as one point is accepted, t is set to 0 again, maintaining the chain's ability to jump to another parameter region. 11inally, we focus solely on points within a chain for which the convergence already occurred (i.e.where maximums of likelihoods are reached).Therefore, we set up a 'burning-length' parameter which automatically removes the first N burn points of each chains.This parameter is once again chosen empirically during the pre-runs by studying multiple likelihood evolution plots.
For clarification and to summarize, we emphasize that the MCMC strategy presented above is associated to Bayesian statistical methods.Our goal is to obtain the posterior distribution for UTZ model parameters while taking into account available experimental constraints.During this process, our algorithm will successively compare parameter points while keeping the best ones, within a probability associated to their respective likelihood ratios.In this approach, prior knowledge is represented by the parameter distributions in the absence of any constraints.We have performed such 'constraint-free' scans and found that every parameter in the scan follows a uniform prior distribution, which demonstrates that the final posterior distributions we present in the following Sections are solely explained by the use of experimental constraints (also see Footnote 12).The scan ranges of UTZ model parameters, along with the value of the model parameter in the global best-fit dataset, for both LO and HO fits.Recall that only two charged-fermion phase parameters are non-redundant at LO [27], and so we have chosen {γ d , δ d }, as in [1].Graphical representations of the MCMC evolution of these parameters are given in Figure 1.Also recall that there are no relevant HO Majorana corrections, that we have kept the HO corrections real, and that the global best-fit values identified for the phases are not terribly meaningful, as we do not observe very strong MCMC preferences for any phase values in our scans (they are all relatively evenly distributed across [0, 2π]).

LO UTZ Model
Given the above discussion, we present the different MCMC hyper-parameters that we use in our setup : N burn , κ 0 , ϵ, the number of chains launched N chains , and the length of the chains L chain which determine the target number of accepted points for every chain.All these parameters are chosen empirically depending on the model and the final statistics desired for the distributions.As a final comment, we note that it is usually better to allow for more chains, rather than longer chains, as this ensures a more reliable parameter exploration.

UTZ Specifics
Following the algorithm above, we now specify the constraints that will guide our likelihood evolution in the MCMC, and also the hyper-and model-parameter choices that control our statistics.Regarding the former, we have identified / implemented the following set of MCMC  constraints and predictions: where R f i f 3 corresponds to the ratio of the ith generation over third-generation mass (R i3 ≡ m f i /m f 3 ) for the corresponding family f , and where 'n.h.' corresponds to the constraint ∆m 2 sol /m 2 ν 1 ≫ 1, which enforces a strictly-hierarchical normal-ordered light neutrino spectrum.The associated numerical constraints correspond to the UV bounds from Tables 2-3.Hence there are N cons = 21 constraints to guide the MCMC likelihood evolution, and N obs = 7 additional predictions that depend on correlated theory parameters, but which do not impact MCMC likelihoods.Observe that sin δ q,l and ξ are listed as quasi-predictions because, as discussed above, the UV bounds associated to them are extremely large, either due to IR experimental uncertainties (sin δ l ) or due to theory uncertainties associated to radiative corrections (sin δ q , ξ).We will therefore use Tables 2-3 as (weak) MCMC constraints, but will also present these results as novel predictions of the UTZ framework, along with those already listed as such in (32).
Given (32), we then set the values of the MCMC hyper-parameters we have employed to 12   N chains = 2500, L chains = 500, N burn = 40, κ 0 = 0.01, ϵ = 0.00005 , while Table 4 gives the ranges scanned over for the actual UTZ model parameters outlined in Section 2. The ranges listed for both were identified from successful preliminary MCMC runs with broader model-parameter ranges, coarser hyper-parameter specifications and, most importantly, general physics considerations from Section 2, which we now discuss.
Considering the Majorana sector, we heuristically observe that establishing the sequential dominance condition in (11) with M 3 /M 2 ∼ 10 n GeV requires max(x, y) ∼ O(10 −n ), and this is largely independent of the scale M θ .We have required n ≥ 3, to truly establish the third-family Majorana dominance implied by (5).Requiring M 2 ≫ M 1 of course requires further suppression between x and y, such that M 2 /M 1 ∼ 10 n GeV (roughly) corresponds to min(x, y) ∼ max(x, y) • 10 −n/2 , and we recall that the qualitative physics leading to (12) does indeed imply said additional hierarchy.However, we also notice from (5) that the two coefficients are sourced from Lagrangian terms that enter at the same power counting (suppressed 12 We have explored the stability of our results under variations of each of these parameters.We already mentioned in Footnote 11 that setting ϵ → 0 does not qualitatively impact our conclusions, a fact that we have also confirmed with respect to κ 0 , by varying it by a factor of three.Furthermore, N burn was determined empirically by conservatively analyzing the likelihood evolution of numerous individual chains; decreasing it will only serve to include lower-likelihood regions of parameter space.Finally, preliminary scans of the UTZ implementing (32) were done using much lower statistics than implied by {N, L} chains , finding results in qualitative agreement with the high-statistics run generated by (33).Hence we believe these choices are quite conservative.
Finally, as mentioned above, we performed a preliminary scan with the same statistics as in (33), but without applying the constraints of (32) (all model parameter configurations generate a likelihood of 1), in order to confirm that the UTZ does not exhibit built-in preferred regions.Taking the generic scan range −5 • 10 −1 ≤ {a f , b f , x, y} ≤ 5 • 10 −1 , we find that all parameter distributions are flat.In other words, the shapes of the distributions presented in Fig. 1 are truly driven by (32) Phase Combos   2-6, to readily compare with [1].
by M 4 ).While the combinations of vevs and coefficients can lead to additional suppression, in Table 4 we have kept the scan range for y on the same generic order of magnitude as x, which of course allows for the additional suppression, but does not enforce it.Scan ranges for the LO Dirac parameters {a, b} f are determined by observing that the LO Dirac Lagrangian in (4) only exhibits one order of messenger mass suppression w.r.t. the leading third-generation scale-setting terms.Allowing for a broad range of Wilson coefficients and flavon vevs, we consider −5 • 10 −1 < {a, b} f < 5 • 10 −1 as a reasonable first constraint on preliminary MCMC scans, which we then iteratively refine given observed preferential domains.Following this procedure, we have noticed that the up-family parameters prefer to be (roughly) symmetrically distributed around zero, and extend to ±O(10 −4 ) (O(10 −3 )) for the a u (b u ) terms.The down-quark and charged-lepton parameters are also symmetric about zero, but with centers around O(10 −3 ) (O(10 −2 )) for a d (b d ).In Table 4 and Figure 1 we have only considered the positive branch of these parameters.Finally, the Dirac neutrino parameters are also distributed in a roughly symmetric way about zero, with both a ν and b ν peaked around O(10 −1 ).
Upon identifying the final LO scan ranges as above, we then consider the HO Dirac parameters {c, d} f , which we recall from (13) contribute at 1/M 4 i,f in the UTZ OPE, i.e. at one order higher than the leading terms.Consistent with our power-counting philosophy at LO, we require these terms be at least one order of magnitude smaller than their LO counterparts.We then consider the analytic hierarchy between the HO operators identified in (15) suggesting the c f correction ∝ S 2 be yet further suppressed w.r.t.d f .To this end, if we have identified a scan range of |min{a, b} f | < O(10 −n ), we require |d f | < O(10 −n−1 ) and |c f | < O(10 −n−2 ).While this of course does not forbid the possibility that |d f | ∼ |c f | (as is also in principle allowed given slightly non-universal messenger masses and/or hierarchical Wilson coefficients), it is sufficiently generic for our purposes and, in any event, we observe that these HO corrections do not converge on highly-preferential domains regardless, which is clearly visible in the last two rows of Figure 1.Note that, for simplicity, we have kept these HO corrections real.
Finally, we note that in all of the above considerations we have allowed for generic LO phase configurations in the neutrino sector, 13 but have chosen the two non-redundant LO phases in the charged fermion sector as in [1], i.e. we allow for non-zero γ d and δ d .This allows us to readily compare the physical conclusions of our analysis with those of [1].In Table 5 we show that this choice is amongst the higher-likelihood configurations given the six possible pairings, having considered other configurations in preliminary MCMC scans with limited statistics.Up to this choice, we have otherwise allowed for arbitrary phases in our MCMC scans; all are constrained to [0, 2π], and we observe that there is typically no strong MCMC preference for said phase domains.For this reason we do not show their MCMC histogram distributions in Fig. 1.

Results and Analysis
We implement the MCMC scan as described above on a computing cluster.The output is a data-set composed of UTZ model parameters, associated values for the constraints and predictions from (32), and the corresponding likelihood of said predictions for each saved MCMC iteration.We denote the corresponding data-set Θ j i , where our notation implies that the ith data-set has j entries corresponding to the model / physical / likelihood parameter(s).Hence i ∈ {1, 2, ..., N Θ }, where i.e the overall number of data-sets, each of which has j ∈ {1, 2, ..., L Θ } constituents, where with the additional unit in L Θ 's counting coming from the standard likelihood function L for a set of model predictions compared to experiment.Given (33), N Θ = 1.15 • 10 6 , and we now examine the physical and model parameters embedded therein.

Fermion Flavour Mass Ratios
We first examine the MCMC results for the UTZ's predictions in the fermion mass sector. Figure 2 illustrates these for mass ratios in the up quark, down quark, charged lepton, and neutrino families.Both Figure 2 (and upcoming figures) and Table 6 give results for the LO and, when indicated, HO MCMC scans, with the former given in blue and the latter in red.
Note that these figures represent density plots, in that darker regions correspond to parameter domains where more Markov chains evolved.Also, the black 'target' marker in Figures 2-6 corresponds to the location of the overall (global) best-fit data-set Θ i , which is also given numerically in Table 6.
The gray bands correspond to the global data available from the PDG (NuFit) collaborations for the charged fermion (neutrino) masses, corrected to the UV according to the discussion in Section 3.1.Comparing these to the blue and red regions, we see that the UTZ is capable of successfully resolving the entire charged fermion mass spectrum, for both quarks and leptons, up to the RGE and threshold correction uncertainties.Furthermore, the UTZ predictions for (currently unmeasured) neutrino mass ratios are shown in the bottom-right panel; given the model parameter ranges explored in Table 4, the ratio R i3 ≡ m ν i /m ν 3 is densely populated within 2.5 • 10 −2 < R 23 < 2 • 10 −1 for the heavier generations while the smaller mass ratio is densely populated between 1.5 • 10 −3 < R 13 < 2 • 10 −2 .However, we see that, albeit less frequent, much larger neutrino mass hierarchies are also resolved, with R 13 (R 23 ) falling below 10 −6 (5 • 10 −3 ).    4. The blue (red) regions correspond to the LO (HO) MCMC scan results, with darker regions corresponding to places of higher density.Gray regions represent the UV bounds for the mass ratios as presented in Table 3, and the black target markers correspond to the global best-fit values shown numerically in Table 6.
Finally, we notice from the up-family plot that the inclusion of the red HO corrections sourced from (13) do not qualitatively change the physics conclusions of the blue LO regions.We have in fact observed this quite generically across family and observable sectors, and hence for visual clarity we only display the dominant LO results in what follows, unless otherwise specified.

Fermion Mixings and CP-Violation
In Figure 3 we have presented the MCMC UTZ predictions for the CKM (top two plots) and PMNS (bottom two plots) mixing angles θ q,l ij as well as the associated Dirac CP-violating phases δ q,l .Here we again compare to (radiatively corrected) data from the PDG and NuFit given in gray, and notice that the blue (red) LO (HO) UTZ Lagrangian is again highly successful at resolving these parameters.Indeed, while we observe that the overlap with PMNS uncertainty bands is perhaps qualitatively more successful than that of the quarks, the regions overlap Note that this conclusion differs from the naive analysis in [1], which found elements in the third row and column of the CKM to be outside of the UV uncertainty bands considering only the LO UTZ Lagrangian, a deviation sourced by the θ q 23 mixing angle.While we observe that the bulk of the MCMC sample points for sin θ q 23 are indeed lower than the allowed uncertainty region, a significant number of LO points do overlap successfully.Studying [12], one concludes that lower values of θ q 23 tend to correspond to higher tan β RGE scenarios.Hence independent evidence that a background spectrum imitating this UV MSSM structure14 is not physical would in principle also disfavor the UTZ theory of flavor, up to the extent the bounds on sin θ q 23 drive our current MCMC likelihoods.
In addition to θ q 23 , the bottom-right panel of Figure 3 suggests that resolving values of sin δ l ∼ 0 simultaneously with θ l 23 in the UTZ is disfavored compared to larger | sin δ l | in the UV.Hence the {θ l 23 , δ l } sector of the PMNS represents an exceptional opportunity to constrain significant portions of the UTZ parameter space, as information on δ l from neutrino oscillations continues to improve.
To fully present the CP-violating sector of the UTZ, we have also presented our MCMC   results for sin δ q,l side-by-side in Figure 4, along with results for the Majorana CP-violating phases ϕ 1,2 of the PMNS matrix.Reliable experimental constraints on ϕ 1,2 are presently nonexistent, and so they also represent opportunities to falsify / further constrain our model space.However, one observes that a broad range of Majorana phases are predicted in the UTZ.We have contextualized this observation by including the MCMC histograms for these phases, analogous to the model parameters presented in Fig. 1, in the bottom two panels of Fig. 4.These histograms reveal that, while it is true that virtually all values of sin ϕ 1,2 are acceptable, a huge number of Markov chains evolved to | sin ϕ 1,2 | ∼ 1. Hence improving data (and therefore more rigid constraints in (32)) could allow us to resolve more precise predictions for these Majorana phases.
Finally, as with the masses presented in Figure 2, we have noticed that HO corrections (as shown in red the top left panel of Fig. 3, and in the histograms of 4) do not qualitatively alter the LO physical conclusions we present above.This is because we have been conservative in Table 4 regarding the relative magnitude of said HO corrections w.r.t.LO parameters, which is of course motivated by the relative suppression of these Lagrangian terms. 15

β-Decay and Cosmological Probes
We now focus on the sector of observables sensitive to the absolute neutrino mass scale and the Majorana (vs.Dirac) nature of the neutrino field, i.e. m Σ and m β(β) .As mentioned above, we only consider ratios of these observables as highly meaningful predictions in the UTZ, and we present their regions in Fig. 5, where it is clear that, at least for the parameter-space we have explored, the UTZ largely prefers values for the ratios m β(β) /m Σ and m ββ /m β given by 5.9 (36) In the event that a positive signal for m β(β) is ever observed, (36) will serve as an excellent probe of the UTZ construction.One also notices in the right panel of Fig. 5 that the MCMC has evolved such that relatively small values of the neutrino-mass-squared difference ratio ξ are preferred, with respect to the possible UV upper bound in Table 3.However, the observed region is still consistent with both the low-tan β / SM-like RGE and high-tan β-like RGE scenarios discussed in Section 3.1. 16f course, as mentioned, we can also report the actual values of the constituent functions m Σ and m β(β) , despite them being less meaningful due to their sensitivity to M θ .For completeness we do so in Figure 6.As expected due to their use as constraint in (32), we observe that the UTZ readily evades available bounds from (e.g.) KATRIN, Planck, and KamLAND-Zen.However we emphasize that this statement effectively amounts to a consistency check on the MCMC framework implemented.

Summary Comments
Before concluding, we summarize the results in the above Sections: • The LO UTZ Lagrangian in (4) and ( 5) is sufficient to describe all available data on fermionic mass and mixing, as well as data constraining the overall scale of neutrino masses.This result is novel, and represents a substantial improvement on the phenomenological findings of [1], which found that HO corrections were necessary to describe the third row and column of the CKM matrix (due to θ q 23 ).This illustrates the power of our MCMC algorithm to robustly explore viable UTZ parameter spaces, in comparison to less sophisticated methods.However, θ q 23 still represents an excellent parameter to exclude UTZ parameter spaces in the future.
• The MCMC algorithm also allows us to present robust predictions for observables that are not well constrained by data -e.g.leptonic CP-violating phases sin{δ l , ϕ 1 , ϕ 2 } and neutrino mass ratios m ν i /m ν 3 , ξ, m β(β) /m Σ , and m ββ /m β -despite the fact that said observables depend sensitively on theory parameters that are highly-correlated to other, well-constrained observable sectors.These findings provide excellent opportunities for the falsification / exclusion of UTZ parameter spaces.We have presented these predictions using the (N)LO UTZ Lagrangian in Figs.2-5.
• The HO corrections generated by the operators in (13) do not qualitatively change the physics conclusions driven by the dominant operators in (4)-( 5).This is due to our (natural) assumption that said HO parameters are suppressed with respect to LO parameters, a constraint that we did not impose as rigorously in [1].As a result, the UTZ's predictions are dominated by as few as nine IR theory parameters.Hence the UTZ is realized as a well-defined, stable, and predictive effective theory of flavour.
• The results we have presented are of course sensitive to the hyper-and model-parameter ranges we have explored, which are presented in (33) and Table 4.While we have taken care in identifying these ranges, and have indeed demonstrated that they are successful, they are not necessarily unique.Exploring alternative parameter spaces, possibly with even more statistics than implied by (33), will be especially motivated in the event data fully excludes the predictions presented in Figures 2-5, and/or a specific renormalizable completion (with exact RGE / threshold behavior) of the UTZ is identified.

Summary and Outlook
We have re-examined the Universal Texture Zero (UTZ) model of flavour presented in [1] in light of updated experimental constraints and in the context of a novel Markov Chain Monte Carlo (MCMC) analysis routine.We have considered the UTZ's predictions at both leading-and next-to-leading orders in its effective theory operator product expansion, and the associated phenomenological pre-and post-dictions are given in Figures 2-5.There we observe that the UTZ is capable of fully resolving the fermionic mass and mixing spectrum as constrained by global data sets, for both the quark and lepton sectors, up to uncertainties regarding radiative corrections to/from the ultraviolet.We have also presented a host of novel, robust predictions for poorly-constrained leptonic observables, in particular the PMNS CPviolating phases δ l , ϕ 1,2 and neutrino mass-sector ratio observables m β(β) /m Σ and m ββ /m β .These latter results offer a route to UTZ falsification and/or parameter-space exclusions.
Our analysis therefore greatly improves on the proof-in-principle phenomenology pursued in the original UTZ publication [1], which was incapable of yielding a robust prediction for even (e.g.) δ l , and which did not consider observables like m β(β) and m Σ .Indeed, our updated MCMC analysis revises the claim from [1] that the leading UTZ Lagrangian is insufficient to account for all fermionic data, before considering next-to-leading corrections.However, as discussed at the end of Section 4, there is still room for improvement, as a yet more exhaustive scan of the UTZ parameter space is in principle possible.We also note that, while our MCMC algorithm fully accounts for theory correlations amongst UTZ model parameters, we have not accounted for experimental correlations, beyond what is already accounted for in the global fits presented in Section 3.While we do not expect such correlations to qualitatively change our conclusions, pursuing such an analysis in the future could be interesting.
Besides these future technical/phenomenological improvements, we also note that significant progress has recently been made in rigorously connecting theories of flavour controlled by non-Abelian discrete symmetries (and additional shaping symmetries) to string theories with toroidal orbifold compactifications -see e.g.[73][74][75][76].It would be interesting to determine whether the UTZ (or a close cousin) could be formally embedded into one of these structures, thereby providing a UV origin for the field and symmetry content of Table 1, and an unambiguous background spectrum that would minimize the radiative uncertainties that we have considered agnostically in our effective field theory setup.After all, the absence of ∆( 27) contractions with non-trivial singlets in (4), (5), and ( 13) is already consistent with the stringy models examined in [77].
We leave these questions to future work, and simply conclude that Figures 2-6 indicate that the UTZ represents an appealing, minimal, and phenomenologically viable model of flavour physics, and therefore provides some support for the idea that observed flavour patterns are the result of yet-discovered Beyond-the-Standard Model dynamics, rather than (e.g.) random chance.

Figure 1 :
Figure 1: Histograms demonstrating the distribution of MCMC iterations for the Dirac (Majorana) scale-setting UTZ parameters {a, b, c, d} d,u,ν ({x, y, M }), in the LO (blue) and HO (red) scans.We have distributed our results across 100 horizontal bins, while the sum of all vertical histogram values in a given plot is equal to N Θ .By and large, phases, like the {c, d} f shown, do not exhibit strong preferential values.

Figure 2 :
Figure 2: MCMC density plots for UTZ quark and lepton flavoured mass ratio predictions.Plots are generated with the hyper-parameter choices in (33) with model-parameter variations as given in Table4.The blue (red) regions correspond to the LO (HO) MCMC scan results, with darker regions corresponding to places of higher density.Gray regions represent the UV bounds for the mass ratios as presented in Table3, and the black target markers correspond to the global best-fit values shown numerically in Table6.

Figure 3 :
Figure 3: The same as Figure 2, but for the CKM and PMNS mixing angles and Dirac CP-violating phases.

Figure 4 :
Figure 4: The same as Figure 3, but now comparing the quark and lepton Dirac CP-violating phases, and presenting the novel predictions for the Majorana CP-violating phases ϕ 1,2 .

Figure 5 :
Figure 5: The same as Figure 2, but now comparing the neutrino mass ratio observables predicted by the UTZ.

Figure 6 :
Figure 6: MCMC scatter results for β-decay and absolute neutrino mass scale observables.Gray regions correspond to IR bounds from KamLAND-Zen (m ββ ), KATRIN (m β ), and Planck (m Σ ), corrected to the UV (if stated) with the conservative s evolution factor for the heaviest generation.As noted in the text, these results are presented as consistency tests of our approach.

Table 5 :
The maximum likelihoods returned from preliminary LO MCMC scans with N chains = 500, L chain = 250, and L burn = 20, upon choosing different charged fermion phase configurations.Note that we have shown the results for the {γ d , δ d } configuration in Figures

Table 6 :
The UTZ global best-fits for fermionic mass (top) and mixing (bottom) parameters in the UV, as determined from the MCMC scan described in Section 4. The upper lines correspond to the fit allowing for only LO UTZ Lagrangian parameters, while the lower lines also account for HO parameters, both of whose MCMC distributions are given in Table4.Figures2-3show the total spread of MCMC predictions in this sector, and also highlight the LO global fits presented in this table with a black target marker.