Flavour bounds on the ﬂavon of a minimal and a non-minimal Z 2 × Z N symmetry

We investigate ﬂavour bounds on the Z 2 × Z 5 and Z 2 × Z 9 ﬂavour symmetries. These ﬂavour symmetries are a minimal and a non-minimal forms of the Z 2 × Z N ﬂavour symmetry, that can provide a simple set-up for the Froggatt-Nielsen mechanism. The Z 2 × Z 5 and Z 2 ×Z 9 ﬂavour symmetries are capable of explaining the fermionic masses and mixing pattern of the standard model including that of the neutrinos. The bounds on the parameter space of the ﬂavon ﬁeld of the Z 2 × Z 5 and Z 2 × Z 9 ﬂavour symmetries are derived using the current quark and lepton ﬂavour physics data and future projected sensitivities of quark and lepton ﬂavour effects. The strongest bounds on the ﬂavon of the Z 2 × Z 5 symmetry come from the D 0 − ¯ D 0 mixing. The bounds on the Z 2 × Z 9 ﬂavour symmetry are stronger than that of the minimal Z 2 × Z 5 symmetry. The ratio R µµ provides rather robust bounds on the ﬂavon parameters in the future phase-I and phase-II of the LHCb by leaving only a very small region in the allowed parameter space of the models.


Introduction
The Z 2 × Z N flavour symmetry [1] provides a new framework for the celebrated Froggatt-Nielsen (FN) mechanism that eventually furnishes an elegant solution to the flavour problem of the standard model (SM) [2].The flavour problem of the SM comprises a set of fundamental questions, including the origin of the mass pattern of fermions of the SM, an explanation for the observed quark-mixing, and the source of neutrino masses and oscillations.There are various approaches to address this problem in literature.For instance, it can have a solution through the hierarchy of vacuum-expectation-values (VEVs) in a technicolour-framework where VEVs are sequential chiral condensates of an extended dark-technicolour sector [3,4,5].A possible explanation can be obtained using Abelian flavor symmetries [2,6,7,8,9,10], creating loop-suppressed couplings to the Higgs [11], through wave-function localization [12] or via compositeness [13].
The central idea of the FN mechanism mechanism is based on an Abelian flavour symmetry U (1) F , which can distinguish different flavours of fermions among and within the fermionic generations in the SM.This is achieved by introducing a flavon field χ in such a way that only top quark gets mass from a renormalized SM interaction, and masses for other fermions are obtained from appropriate non-renormalized higher dimensional operators, which are constructed using the flavon field χ.For instance, if under the U (1) F symmetry the fermions ψ c i and ψ j have charges θ i and θ j respectively and the charge of the the SM Higgs field is zero then the Yukawa Lagrangian of the SM is forbidden by the U (1) F symmetry.In this scenario, the masses of the SM fermions can be recovered by the non-renormalizable operators of the form, = y (θ i +θ j ) ψϕψ = Y ψϕψ where y is the dimensionless coupling constant, Λ is the scale at which these operators are renormalized, = χ Λ , Y = y (θ i +θ j ) is the effective Yukawa coupling, and the gauge singlet flavon scalar field χ transforms under the SU (3) c × SU (2) L × U (1) Y symmetry of the SM as, The U (1) F flavour symmetry is broken spontaneously when the flavon field χ acquires a VEV.The scale Λ is not provided by the theory, and it can be anywhere between the weak and the Planck scale.
We only require that the flavour symmetry should be broken weakly which means the ratio χ Λ should be less than unity.The flavon exchange effects in the SM phenomenology will be highly suppressed if the scale of new physics Λ is much larger than the weak scale.However, if the flavour symmetry is broken close to the weak scale, we can hope to see observable effects on the direct or indirect experimentally measured SM observables such as mixing and CP -violation in mesons.Therefore, we need to ask how low the flavour scale could be such that it respects the bounds on flavour-changing and CP-violating processes.Moreover, the nature of the flavour symmetry also plays an important role in the investigation of the flavour scale.For example, if the flavour symmetry is a continuous U (1) F , then we should ask whether it is a gauged or a global symmetry.In a gauged U (1) F scenario, the phenomenology of the flavon field will be affected by the exchange of the corresponding gauge boson.
If the continuous U (1) F is global, then a massless Goldstone boson must exist.The Z 2 × Z N flavour symmetry, unlike the conventional continuous U (1) flavour symmetry that is employed to achieve the FN mechanism, is a product of two discrete symmetries which can implement the FN mechanism in a unique way such that the flavour structure of the SM including neutrino masses and mixing parameters can be parametrized in terms of a small parameter which is the ratio of the VEV of the flavon field and the flavour scale Λ [1].We notice that the origin of the Z 2 ×Z N flavour symmetry may be traced to an underlying Abelian or non-Abelian continuous symmetry or their products.For instance, the Z 2 × Z N symmetry may be a by-product of a spontaneous breaking of U (1) × U (1) continuous product symmetry.
We note that the discrete Z 2 symmetry is extensively used in studying the different versions of the two-Higgs-doublet model (2HDM) and the minimal supersymmetric SM (MSSM).In particular, in the Z 2 × Z N flavour symmetry, the discrete Z 2 symmetry exactly behaves like the one used in the type-II 2HDM [14].Therefore, the Z 2 × Z N flavour symmetry may also be used to implement the FN mechanism in the type-II 2HDM and the MSSM.Moreover, the discrete Z 2 symmetry is also found to be useful in model building, for instance, see references [3,5,15,16,17,18].
In this work, we investigate flavour bounds on the dynamics of the flavon field of a minimal and a non-minimal form of the Z 2 × Z N flavour symmetry that provides a simple set-up for the FN mechanism.We do not consider any ultraviolet completion of the Z 2 × Z N based FN mechanism, and present our results in a model-independent manner.The phenomenological investigations of the flavon field of the FN mechanism in the framework of a continuous U (1) symmetry and its extensions are dedicatedly performed in literature, for instance, flavour bounds are investigated in reference [19], the LHC phenomenology is explored in references [20,21,22,23,24,25,26,27], a low flavour breaking scale is studied in reference [28], a study for a future high energy collider is presented in reference [29], flavon exchange effects in the dark matter interactions are studied in reference [30], the texture based investigation of the FN mechanism can be found in reference [31].
We shall present our phenomenological analysis along the following line: In section 2, we investigate a minimal form of the Z 2 × Z N flavour symmetry that can implement the FN mechanism.An explanation to neutrino masses and mixing parameters is discussed in section 2.4.1.A non-minimal form of the Z 2 × Z N flavour symmetry that implements the FN mechanism is discussed in section 3. The scalar potential of our model is discussed in section 4. Phenomenological bounds based on the quark flavour physics on the parameter space of a minimal and a non-minimal form of the Z 2 × Z N flavour symmetry are derived in section 5. Leptonic flavour constraints are investigated in section 6.A summary of the work is presented in section 7.

A minimal Z × Z N flavour symmetry
We now discuss the question of a minimal form of the Z 2 × Z N flavour symmetry that can provide a simple set-up of the FN mechanism.Our guiding principle for this purpose is the observation that a minimal suppression of the effective Yukawa couplings will require a minimal form of the Z 2 × Z N flavour symmetry.For instance, we assume that the mass of the top quark originates from the tree level SM Yukawa operator, then, following the principle of minimum suppression (PMS), the mass of the bottom quark is obtained from the operator having the suppression of the order y , the mass of the charm quark from the operator having the suppression of the order y 2 , the mass of the strange quark from the operator having the suppression of the order y 3 , and the mass of the up and down quarks from the operators having at least the suppression of the order y 4 .
Additionally, we need to count the number of hierarchical energy scales needed to account for the fermionic mass hierarchy in the SM.For instance, for the quark sector, we need three energy scales to explain the mass hierarchy among the three fermionic families.We note that only the second and the third quark families have intra-generational mass hierarchies, which require only two hierarchical energy scales to achieve an explanation for the mass hierarchy within the second and third quark families.These hierarchical energy scales are created by the different non-renormalizable operators of the flavon fields as given in equation 1.Since the mass of the top quark is generated by the renormalized SM Yuakwa operators, at least four energy scales are required to be created through the operators of the form given in equation 1 for providing an explanation for the hierarchical quark mass pattern.
The symmetry Z N in the Z 2 × Z N flavour symmetry is responsible for providing such operators.The symmetry Z 2 distinguishes between the up-type and the down-type quarks, which makes sure that identical non-renormalizable operators of the flavon fields, as given in equation 1, do not appear in the up and down-type quark mass matrices.For creating four energy scales, the required Z N symmetry, therefore, should have at least four non-trivial charges.Therefore, the size of a minimal symmetry will be determined by this requirement and through the application of the PMS.
Finally, we must note that a minimal form of the Z 2 ×Z N flavour symmetry should not only produce correct pattern of the charged fermion masses, it should be capable of explaining the quark mixing pattern, neutrino masses, and more importantly, it should predict correct pattern of the neutrino mixing angles.
After taking into account above considerations, Z 2 × Z N flavour symmetry allows us to write the following generic Lagrangian which provides masses to the charged fermions of the SM, where χ or χ † may appear in the numerator of the term inside the square brackets.We note that in the above Lagrangian i and j are family indices, ψ q L , ψ L are quark and leptonic doublets, ψ u R , ψ d R , ψ R are right-handed up, down type singlet quarks and leptons, ϕ and φ = −iσ 2 ϕ * are the SM Higgs field and its conjugate and σ 2 is the second Pauli matrix.The effective Yukawa couplings Y ij are defined in terms of the expansion parameter

Z 2 × Z 2 flavour symmetry
The simplest choice is the Z 2 × Z 2 flavour symmetry, which turns out to be a trivial selection since the only charges of the Z 2 symmetry are ±1, which are too trivial to provide four energy scales or equivalently non-trivial operators of the form given in equation 1.Hence, we conclude that this symmetry cannot be used to create a simple FN mechanism.

Z 2 × Z 3 flavour symmetry
The first non-trivial form of the Z 2 × Z N flavour symmetry, which may provide an implementation of the FN mechanism, is Z 2 × Z 3 .The symmetry Z 3 has two non-trivial charges characterized by ω and ω 2 , where ω is the cube root of unity.In the first scenario, following the PMS, we assign the charges to the SM and flavon fields as given in table 1. Fields Table 1: The charges of left and right-handed fermions of three families of the SM, Higgs and the flavon field under the Z 2 and Z 3 product symmetry, where ω is the cube root of unity.
We observe that the masses of s and b quarks can be recovered from this charge assignment, for instance, mass of the s quark is of the order 3 , and that of the b quark is of the order .However, the mass of the u and d quarks are produced by the operators y( Any other charge assignment also does not reproduce masses of every quark.
As an additional check, we may assume that exactly identical diagonal operators for the u and d quarks in their mass matrices as given in table 2. This charge assignment is against the original sprite of the Z 2 × Z N flavour symmetry, where the Z 2 is exactly like the symmetry used in the type II 2HDM.It turns out that even in this case, the mass of the u quark is of the order .This conclusion does not change even if we provide different non-trivial charge assignments to the fermions and flavon fields under the Z 2 × Z 3 flavour symmetry.
Table 2: The charges of left and right-handed fermions of three families of the SM, Higgs and the flavon fields under the Z 2 and Z 3 product symmetry, where ω is the cube root of unity.

Z 2 × Z 4 flavour symmetry
The charges of the symmetry Z 2 × Z 4 are characterized by the fourth roots of unity, which are 1, ω, ω2 and ω 3 where ω = i, ω 3 = ω * and ω 2 = −1.We particularly note that minimally suppressed diagonal operator of the form y( χ Λ ) 4  ψL ϕu R is not the dominant operator no matter what charge we assign to the flavon and fermionic fields.This is because the tree-level SM Yukawa operator y ψL ϕu R is allowed for any charge assignment for the y( χ Λ ) 4  ψL ϕu R operator under the Z 2 × Z 4 flavour symmetry.
Table 3: The charges of left and right-handed fermions of three families of the SM, Higgs and the flavon fields under the Z 2 and Z 4 product symmetry, where ω is the fourth root of unity.
Therefore, to produce the mass of the u-quark, we either choose a non-trivial transformation of the u R -quark or the first family of the quarks under the Z 2 symmetry, in addition to the next to the minimal suppressed operator of the order 5 .One such charge assignment is given in table 3.In this case, the masses of the down-type quarks are produced correctly through the operators with minimal suppression.However, in the case of up-type quarks, still non-diagonal tree-level SM operators dominate the mass of the u-quark.Other alternative charge assignments also do not work for creating an FN mechanism through the Z 2 × Z 4 symmetry.

Z 2 × Z 5 flavour symmetry
We now impose the next flavour symmetry, that is, the Z 2 × Z 5 symmetry on the SM in a way that the various fields of the SM transform under this symmetry as given in table 4 [1].As discussed earlier, we need to create at least four hierarchical energy scales as an origin of the quark mass spectrum.This means, for creating these energy scales, a non-trivial and a minimal Z N symmetry should have at least four non-trivial charges.Thus, the symmetry Z 5 could be such a symmetry.Moreover, we note that the transformation of fields under the Z 2 × Z 5 symmetry is chosen such that the symmetry Z 2 exactly acts like the way used in the type-II 2HDM 1 .
Table 4: The charges of left and right-handed fermions of three families of the SM, right-handed neutrinos, Higgs, and singlet scalar fields under Z 2 and Z 5 symmetries, where ω is the fifth root of unity.
The Z 2 × Z 5 flavour symmetry allows us to write the following Lagrangian which provides masses to the charged fermions of the SM, The mass matrices for up-and down-type quarks and charged leptons can be written now in terms of the expansion parameter , (4) The masses of quarks and charged leptons approximately are [32] The mixing angles of quarks are found to be [32], We notice that the sin θ 12 and sin θ 23 have the same order.The similar result is also reported in reference [7].
We present a fit of the experimental data to the masses of fermions in appendix.It turns out that some of the couplings are not order one.We discuss a theoretical scenario for such couplings in the appendix.This kind of not order one couplings are also reported in references [7,19].

Neutrino masses and mixing
The neutrino masses are obtained by adding three right-handed neutrinos as shown in table 4. The Lagrangian for the tree-level Majorana mass is, where i, j are flavour indices.The Majorana mass matrices M R is, where The Dirac mass Lagranginan for neutrinos can be written as, The Dirac mass matrix is given by, The mass matrix of neutrinos after including the Majorana mass terms can be written as, Since v << f , we ignore the contribution of the mass matrix M L to the neutrino masses. 2 Now, we can use the type-I seesaw mechanism to determine the neutrino masses by assuming M D << M R [33].Thus the light neutrino mass matrix is, where = v f 3 , and we have assumed each and every coupling exactly one except those appearing in above equation.
We obtain two degenerate neutrino masses.The masses approximately are given by, This kind of approximate degenerate neutrino masses are well studied in literature, for instance, see references [34,35,36,37].The leptonic mixing angles can be written as, The remarkable observation is the pattern of the neutrino mixing angles.The mixing angle θ 12 and θ 13 are of the same order of magnitude, where θ 13 is closer to the Cabibbo angle, and the mixing angle θ 23 is completely unsuppressed. Fields Table 5: The charges of left and right-handed fermions of three families of the SM, right-handed neutrinos, Higgs, and singlet scalar field under Z 2 and Z 9 symmetries, where ω is the ninth root of unity.

A non-minimal Z 2 × Z 9 flavour symmetry
We note from the previous section that some of the Yukawa couplings for the minimal model based on the Z 2 × Z 5 flavour symmetry are not order one, which is a preferred choice in literature.However, so far purpose has been to introduce the Z 2 × Z N flavour paradigm.In this section, we show a nonminimal model based on the Z 2 × Z N flavour paradigm where the Yukawa couplings turn out to be order one, and are given in the appendix.We adopt a non-minimal Z 2 × Z 9 flavour symmetry, and assign the charges to different fields as shown in table 5.
The Z 2 × Z 9 flavour symmetry allows us to write the following Lagrangian which provides masses to the charged fermions of the SM, The mass matrices for up and down-type quarks and charged leptons turn out to be, (18) The masses of charged fermions are approximately given by [32], Similarly, the mixing angles of quarks read [32],

Neutrino masses and mixing
The masses and mixing of neutrinos in the non-minimal model is identical to that of the minimal model.Thus, we write the Majorana Lagrangian for right-handed neutrinos as, where i, j are flavour indices.The Majorana mass matrices M R can be written as, where The Dirac mass Lagranginan for neutrinos is, The Dirac mass matrix for neutrinos now reads, The mass matrix of neutrinos after including the Majorana mass term is, The light neutrino mass matrix is,

where
= v f , and we have again assumed each and every coupling exactly one except those appearing in above equation.
The neutrino masses approximately are,

The scalar potential
The scalar potential of the model can be written in the following form, where we have introduced a soft breaking of the Z 5 symmetry in the fifth term.We are assuming λ ϕχ = 0, i.e., no Higgs-flavon mixing [29].If this term is non-zero, the phenomenology of the flavon field will be different, for instance, see reference [21].The only parameter which can have a phase in the scalar potential is ρ.However, this phase can be removed by a phase rotation of the flavon field χ leading to a real value of the VEV of the field χ.
We can parametrize the flavon field by excitations around its VEV, In a similar manner, the Higgs field can be written as, The minimization conditions can be written in terms of the scalar and pseudo-scalar components having the following masses: We observe that the mass of the pseudoscalar component of the flavon field depends on the softbreaking parameter ρ.Therefore, it is a free parameter of the model.Now using equation 33, we can write The couplings of the scalar and pseudoscalar components of the flavon field are obtained from equation 26 by writing the effective Yukawa couplings in the following form: where f = u, d, , and n f ij is the power of the parameter appearing in the mass matrices M f .We note that for our phenomenological investigation, we have only retained the terms linear in the flavon field components s and a in equation 37.The terms which are higher than the linear terms are not interesting in the present work.The couplings of the Higgs boson field h to the charged fermions are real and diagonal since the mass matrices M f can be diagonalized resulting in real and positive masses of the charged fermions.However, the couplings of the scalar and pseudoscalar components s and a of the flavon field are given by n f ij M f .This product cannot be diagonalized exactly and as a consequence, the couplings of s and a cannot be made real and diagonal.This in turn gives rise to the flavour-changing and CP -violating interactions of the flavon field.
The couplings of a field with fermions for minimal Z 2 × Z 5 symmetry are now given by, In the similar way, the couplings of a field with fermions for non-minimal Z 2 × Z 9 symmetry are given by For the pseudoscalar component of flavon field, the following notation is used: The couplings of a to fermions are identical to that of s except with a relative phase factor i. This factor becomes trivial in the squared amplitude of a Feynman diagram mediated by a. Thus, it does not give rise to CP -violating interactions to the order investigated in this work.
5 Quark flavour physics in the minimal and non-minimal Z 2 ×Z N flavour symmetry The quark flavour physics places stronger bounds on the parameter space of our model, which is parametrized by the VEV of the flavon field f , the mass of the pseudoscalar flavon m a , and the quartic coupling λ χ .In particular, measurement of the loop-induced processes in the SM such as neutral meson mixing and rare mesonic decays constrain the parameter space of the model.The numerical inputs used in this work are given in table 6.

Neutral meson mixing
The non-diagonal couplings of the flavon to fermions introduce the FCNC interactions at tree-level.Therefore, they are expected to be highly suppressed from neutral meson-antimeson mixing.These interactions of the flavon can be parametrized by writing the ∆F = 2 effective Hamiltonian as follows, where q R,L = 1±γ 5 2 q and the colour indices are omitted for simplicity.Table 6: Values of the experimental and theoretical quantities used as input parameters.
The tree-level contribution to neutral meson mixing due to the flavon exchange gives rise to the following Wilson coefficients [47,48], where m s and m a are the masses of scalar and pseudoscalar component of flavon field, respectively.
The Wilson coefficients C i are computed at a scale Λ, where heavier new degrees of freedom are integrated out.They need to be evolved down to the hadronic scales 4.6 GeV for bottom mesons, 2.8 GeV for charmed mesons, and 2 GeV for kaons.These particular scales are used in the lattice computations of the corresponding matrix elements [40,44,45].In this work, renormalization group running of the matrix elements is implemented as discussed in reference [45] and matrix elements are taken from reference [40,44].Thus, the new physics contribution to the B q − Bq mixing amplitudes due to the Wilson coefficients C i at a scale Λ can be written as [45], where q = d, s, α s is the strong coupling constant, η = α s (Λ)/α s (m t ), and a j , b are the so-called the magic numbers which are taken from reference [49].We can write a similar formula for D 0 − D0 mixing with magic numbers given in reference [45].For K 0 − K0 mixing, the formula becomes [45], where R r are the ratio of the matrix elements of NP operators over that of SM [50] and their numerical values are directly taken from reference [45] for our analysis.The magic numbers for K 0 − K0 mixing are taken from reference [40].
The mixing observables of the K 0 − K0 mixing can be used now to constrain the flavon mass and VEV by employing their experimental measurements.These are [45], where numbers are given at 95% C.L., H ∆F =2 eff contains the SM and flavon contributions, and H ∆F =2 SM represents only the SM contribution.
The mixing observables for the B q − Bq mixing are, q where q = s, d for B s and B d mixing respectively.The following measurements at 95 % CL limits are used in this work [45], The new physics contributions to neutral meson mixing can be written as, We assume the minimum flavour violation scenario that corresponds to σ d,s,K = 0. We adopt the following future sensitivity phases in this work [51] Figure 1: The allowed parameter space by flavour observables C K and C ∆m K in the m a − f plane for the minimal (Z 2 × Z 5 ) flavour symmetry.On the left panel in figure 1a, the allowed bounds for λ χ = 2 with current limits for C ∆m K and C K are shown by red and violet boundaries, respectively.Also, the allowed bound with projected limits of C ∆m K is shown with yellow boundary.The effect of the variation of the quartic coupling λ χ on the observable C K is shown on the right panel in figure 1b.
In figure 1, we show the bounds on the VEV of the flavon and the mass of the pseudo-scalar flavon arising due to the neutral kaon mixing observables C K and C ∆m K for the minimal model based on the Z 2 × Z 5 flavour symmetry.On the left, the allowed region by the observables C K and C ∆m K is shown for the quartic coupling λ χ = 2.There is a sudden dip in the allowed parameter space given by C K which appears due to a cancellation in the Wilson coefficients C ij 2 and Cij 2 when masses of scalar and pseudoscalar flavon become identical.For C ∆m K , this dip is not visible in this plot, and is excluded for the quartic coupling for λ χ = 2.The region bounded by the yellow curve is the allowed parameter space by the future projected sensitivity as shown in table 7. On the right panel, we show the allowed regions of the parameter space by the observable C K for λ χ = 0.5, 2, 4π.It is observed that the allowed region shrinks as λ χ approaches smaller values.
Similar results for non-minimal model based on the Z 2 × Z 9 flavour symmetry are shown in figure 2. The constraints on the allowed parameter space by the neutral kaon mixing observables C K and C ∆m K , in this case, are more stringent in comparison to that of the minimal (Z 2 × Z 5 ) flavour symmetry.In particular, we observe that the region with the sudden dip is excluded for the non-minimal model.Moreover, there is no allowed parameter space for the future projected sensitivity of the observable C ∆m K for the bench-mark values of the couplings given in the appendix.Therefore, we show these bounds for different values ( Figure 2: The allowed parameter space by flavour observables C K and C ∆m K in the m a − f plane for the non-minimal Z 2 × Z 9 flavour symmetry.On the left panel in figure 2a, the allowed bounds for λ χ = 2 with current limits for C ∆m K and C K are shown by green and magenta boundaries, respectively.Also, the allowed bound with projected limits of C ∆m K is shown with olive coloured boundary.The effect of the variation of the quartic coupling λ χ on the observable C K is shown in the right panel in figure 2b.
We show the allowed regions of parameter space by the B s − Bs mixing observables C Bs and φ Bs for λ χ = 2 in the m a − f plane for the minimal model based on the Z 2 × Z 5 flavour symmetry and the non-minimal model based on the Z 2 × Z 9 flavour symmetry in figure 3.In the left panel in figure 3a, the red and yellow coloured boundaries are representing allowed flavon contribution for current values of C Bs and φ Bs , respectively in the minimal model, while that for the non-minimal model are shown in the right panel in figure 3b, by magenta and green coloured boundaries, respectively.We note that the region in the left panel, bounded by the blue curve, represents the allowed bounds for the observable C Bs with projected limits of LHCb Phase-II for the minimal model, while the same for the non-minimal model is shown by the region surrounded by the olive coloured curve in the right panel.For the LHCb phase-I, the bounds are shown by the pink coloured curves for the minimal as well as for the non-minimal models, which are not appreciably different than that of the LHCb phase-II.We also notice an isolated allowed strip of the parameter space for the non-minimal model below the green boundary in the right panel.The effects of the observables C Bs and φ Bs are relatively strong in the non-minimal model based on the Z 2 × Z 9 flavour symmetry, which is obvious from the figure itself.The SM contribution to D 0 − D0 mixing is marred by large hadronic uncertainties.Therefore, for constraining the parameter space of our model, we keep only the flavon contribution to D 0 − D0 mixing such that it always lies within the 2σ experimental bound [53].
The bound in the m a − f panel arising from the D 0 − D0 mixing is shown in figure 5 for the minimal model based on the Z 2 × Z 5 flavour symmetry and the non-minimal model based on the Z 2 × Z 9 flavour symmetry.The first remarkable observation is the allowed parameter space for the minimal model is much smaller than that of the non-minimal model.This is because the D 0 − D0 mixing has an enhancement of the order 2 (see the coupling 2y u 21 2 in equation 38) in the minimal model based on the Z 2 × Z 5 flavour symmetry, which is not present in the case of the non-minimal model.Therefore, the bound derived by the D 0 − D0 mixing is the most stringent bound among the bounds given by the mixing observables in the minimal model.However, this is not the case for the non-minimal model.
The one-loop contribution to this mixing from the box diagram depends on the relatively large y ct and y tc couplings of the flavon to fermions.In the minimal model, y tc is zero.Thus, this contribution is proportional to 4 /(4π 2 f 2 ) in the minimal model and proportional to 16 /(4π 2 f 2 ) in the non-minimal model.Therefore, this contribution is highly suppressed with respect to the tree-level contribution used in deriving the bounds in figure 5.

Leptonic decays of mesons
The effective Hamiltonian for flavon mediated decays of neutral mesons into two charged leptons can be written as, The branching ratio of a meson decaying to two charged leptons reads, where β(x) = √ 1 − 4x 2 with x = m /m M .The Wilson coefficients having tree-level contribution of the flavon are given as [47,48], In the SM, processes of mesons decaying to two charged leptons are induced by one-loop contribution, and for the B s meson it is given by [48], where Inami-Lim function Y (x) is given by [54], where η QCD = 1.0113 includes NLO QCD effects [55].For B d meson, the SM predictions are obtained by a simple replacement of indices in equation 51.
The average of the branching fraction of B s → µ + µ − from HFLAV group is [56], The latest measurement of the branching fraction of As observed in reference [59], due to sizeable width difference, of the B s meson, theoretical branching ratio can be converted to experimental branching ratio by multiplying (1 − y s ) −1 , where y s = 0.088 ± 0.014 [60].This correction is negligible in the case of the B d meson.
We have also taken HFLAV measurement average into account for effective lifetime τ µµ , which is [56], with BR(B s → µ + µ − ) given in equation 53.
The current and future sensitivities of these observables are summarized in table 8.The effective lifetime can be written in the following form [63], where we have assumed the SM value of the final state dependent observable Figure 6: The allowed parameter space by BR(B d → µ + µ − ) in the left panel and that for the BR(B s → µ + µ − ) in the right panel for λ χ = 2 with current experimental bounds for the minimal and the nonminimal models.
In figure 6, the bounds coming from the branching ratios BR(B d → µ + µ − ) and the BR(B s → µ + µ − ) for λ χ = 2 in the m a − f plane are shown.The branching ratios BR(B d → µ + µ − ) and BR(B s → µ + µ − ) place weaker constraints on the parameter space of the minimal model.However, for the non-minimal model, the bounds, particularly from the branching ratios BR(B s → µ + µ − ), are quite stronger.For the projected sensitivities of the LHCb Phase-I and II and of the ATLAS for the BR(B s → µ + µ − ), we do not obtain any appreciable improvements in our bounds.Therefore, we do not show them in figure 6.
Figure 7: The allowed parameter space by τ µµ for λ χ = 2 in the m a − f plane for the recent measurement and the future projected sensitivity of the LHCb for the minimal (Z 2 × Z 5 ) model on the left, and for the non-minimal (Z 2 × Z 9 ) model on the right.
Figure 7 shows the allowed bounds derived from the HFLAV average of the effective lifetime τ µµ and the branching-ratio BR(B s → µ + µ − ) by the green coloured curve for the minimal model in figure 7a and through the olive coloured boundary for the non-minimal model in figure 7b.The bounds from the current measurement by CMS, and from the future projected sensitivity of the LHCb phase-II for the minimal model based on Z 2 × Z 5 , are shown by red and yellow coloured boundaries respectively, in figure 7a.The same bounds for the non-minimal model based on the Z 2 × Z 9 flavour symmetry are shown in figure 7b by green and magenta coloured curves respectively.We observe that the bounds are stronger for the non-minimal model in comparison to that of the minimal model.For the projected sensitivity of the LHCb Phase-I, we do not find any appreciable improvement in the allowed bounds over that from the current measurement.Therefore, we do not show it in figure 7.
Figure 8: The left panel represents the allowed parameter space by R µµ for the minimal (Z 2 × Z 5 ) model with λ χ = 2 with the current measurement and for the future projected sensitivity of the LHCb Phase-I and LHCb Phase-II.The same allowed parameter space for the non-minimal (Z 2 × Z 9 ) model is shown in the right panel.
One of the most important observables for our models is the ratio R µµ whose future projected sensitivity will be crucial to constrain our models.On the left panel in figure 8a, we show the bounds arising from the ratio R µµ for the minimal model.These bounds are weaker for the current measurement as well as for the future projected sensitivity of the LHCb Phase-I shown by the red and yellow boundaries, respectively.However, the future projected sensitivity of the LHCb Phase-II dramatically changes this scenario and provide extremely stringent constraints on the parameter space of the minimal model shown by the purple coloured strip.Therefore, the LHCb Phase-II will be decisive for the minimal model based on the Z 2 × Z 5 flavour symmetry.We also do not find any improvement over the bounds given by the LHCb Phase-I using the future projected sensitivity of the CMS experiment.Therefore, we do not show it in figure 8a.
On the other hand, for the non-minimal model based on the Z 2 × Z 9 symmetry, the bounds from the ratio R µµ are shown in figure 8b in the right panel.The bounds from the current measurement are shown by the green boundary while the bounds from the LHCb Phase-I are surrounded by the purple coloured curve.Moreover, we also have bounds from the future projected sensitivity of the CMS experiment surrounded by olive coloured boundary.These are more stringent than that of the projected sensitivity of the LHCb Phase-I .We note that similar to the minimal model, the bounds for the projected sensitivity of the LHCb Phase-II are highly stringent depicted by the orange coloured strip.This result does not change even if we deviate from the bench-mark values of the Yukawa couplings used for these bounds.Therefore, these bounds are robust, and crucial to test the parameter space of the non-minimal model based on the Z 2 × Z 9 symmetry in the future projected sensitivity of the LHCb Phase-II.
For the K L → µ + µ − decay, we have only reliable estimate of the so-called short distance (SD) part of K L → µ + µ − decay [47].We use the SM prediction obtained in reference [47] and given by where at NNLO Y NNL = λ 4 P c (Y ), λ = |V us | and P c (Y ) = 0.113 ± 0.017 [64].The short distance contribution can be extracted from the experimental measurement and has an upper limit [48], For the case of D → µ + µ − decay, the SM contribution is plagued by large non-perturbative effects.Therefore, we only require that the flavon contribution does not generate more than the experimental upper bound on the branching ratio that is given at 90% C.L. [65], BR(D → µ + µ − ) < 6.2 × 10 −9 .
(62)  We show bounds arising from the BR(K L → µ + µ − ) SD for λ χ = 2 for the minimal and the nonminimal models in the m a − f plane in figure 9.It turns out that the constraints on the parameter space of the minimal model are weaker.However, the bounds on the parameter space of the nonminimal model are more stringent.The constraints from the D → µ + µ − decay are much weaker than that from the BR(K L → µ + µ − ) SD .Therefore, we do not show them in this work.Belle II [69] Table 9: Experimental upper limits on various Leptonic flavour violation (LFV) processes.

Radiative leptonic decays
The effective Lagrangian for the radiative leptonic decays can be written as, The radiative leptonic decays are mediated by dipole operators and their branching ratio reads, The one-loop contribution to the radiative leptonic decays is shown in figure 10.The corresponding Wilson coefficients are [28],  The radiative leptonic decays µ → eγ places weak constraints on the parameter space of the minimal model as shown in figure 11a.This does not change much even for the future projected sensitivities of the MEG-II experiment.For the non-minimal model, the bounds from the current measurement of the MEG experiment, shown in figure 11b by the green boundary, are quite stringent relative to that of the minimal model.Moreover, for the future projected sensitivities of the MEGII experiment, the constraints are further stronger depicted by the magenta coloured boundary.The constraints coming from the decays τ → eγ and τ → µγ are weaker than that of the decay µ → eγ.Therefore, they are not shown in this work.

A µ → A e conversion
The effective Lagrangian describing A µ → A e conversion can be written as, Moreover, there is additional contribution to A µ → A e conversion from the dipole operators given in equation (63).The Feynman diagram for A µ → A e conversion is shown in figure 12.The corresponding Wilson coefficients arise due to the diagram on the left in figure 12 [28], The nuclear effects which include effects of quarks inside the nucleons as well as the contribution of the Feynman diagram on the right side of figure 12 are absorbed in the nucleon-level Wilson where the quark content of the proton is accounted by vector and scalar couplings f p Vq , f p q , and f p heavy = 2/27 1 − f p u − f p d − f p s [80].For right-handed operators, analogous expressions are obtained by replacing L with R, and for the neutron p is replaced by n.The vector operators contribute extremely less than the scalar operators and can be neglected [28].The numerical values of vector and scalar couplings are taken from references [81,82], which are based on the lattice average given in reference [83], The A µ → A e conversion rate including nuclear effects can be written as [28],  13a for the current measurement and for the different future projected sensitivities, and the same for the nonminimal model are shown in figure 13b.The strongest bounds among them arise from the projected future sensitivities of the PRISM/PRIME experiment for the minimal as well as non-minimal models.

µ → 3e and τ → 3µ decays
The three body flavour violating leptonic decays µ → 3e and τ → 3 where = e, µ provide additional tests of the dipole operators given in equation (63).Their decay width can be written as [28], where the tree-level contribution is ignored due to the strong chiral-suppression which is dominated by the logarithmic enhancement of the dipole operators [28].Other contributions, such as Z-mediated penguin are strongly suppressed and ignored [85].In figure 14, we show the allowed parameter space by BR(µ → 3e) for λ χ = 2 in the m a − f plane.For the minimal model, the bounds are shown in figure 14a.They are weaker for the current measurement and for the future projected sensitivity as well.In figure 14b, the same bounds on the allowed parameter space are shown for the non-minimal model.For the current measurement, the bound is given by the green boundary, and for the future projected sensitivity, it is depicted by the magenta coloured curve.As obvious from figure 14b, the bounds for the non-minimal model are quite stringent relative to the minimal model.The bounds from the decay τ → 3 are much weaker and are not shown in this work.Finally, we show a summary of the most relevant and stringent bounds on the prameter space of sensitivities of the experiments such as MEG II, Mu3e, DeeMe, COMET, Mu2e and PRISM/PRIME will be able to constrain the parameter space of the non-minimal model based on the Z 2 × Z 9 flavour symmetry.On the quark side, the B d − Bd mixing in the future phase-I and II of the LHCb will be able to eliminate a sufficient region of the parameter space.In future phase-II of the LHCb, the ratio R µµ will be crucial in ruling out the major part of the flavon parameter space.Thus, future projected sensitivities of the LHCb phase-I and II will play a defining role in determining the fate of the flavon models discussed in this work.

Outline of a possible ultraviolet completion of the Z 2 × Z N model
We present an outline of the underlying renormalizable theory, which could be suitable to the models discussed in this work using the idea discussed in reference [5].Let us suppose that the underlying theory is a technicolour (TC) theory containing two technicolour symmetries.The SM Higgs field comes from the conventional TC group and the flavon field χ is derived from a different dark technicolour symmetry (DTC).
The TC chiral condensate which play the role of the SM Higgs VEV can be parametrized as, where ∆χ TC is the chirality of the operator on the left of above equation, Λ TC is the scale of the underlying gauge TC theory, and k TC is a constant.In a similar manner, we can write DTC chiral condensate which represents the flavon VEV, Now the couplings y ij are given by the following equation: where the scale Λ ETC corresponds to the extended TC theory in which the SM, TC and DTC fermions are embedded.
The masses of the fermions can be written as, We observe that if the function f (Λ TC , Λ DTC , Λ ETC ) is always generated by a tree-level exchange of the underlying theory, all the couplings y ij will have the same order of magnitude (which could be order one).However, depending of the structure of the underlying theory, some of the y ij could be loop-induced and some could come from tree-level contributions.Therefore, this will result in some couplings being suppressed compared to the tree-level contributions.This scenario will lead to the numerical couplings for the minimal model based on the Z 2 × Z 5 flavour symmetry.

|y d 12 | = 1 and |y d 21 |
= π) of the couplings which are allowed by a more relaxed fit.

2 × Z 5 )Figure 3 :
Figure 3: The parameter space allowed by flavour observables C Bs and φ Bs for λ χ = 2 in the m a − f plane for the minimal (Z 2 × Z 5 ) model in the left panel and for the non-minimal (Z 2 × Z 9 ) model in the right panel.

Figure 4
Figure 4 shows the allowed parameter space by flavour observables C B d and φ B d of the B d − Bd mixing for λ χ = 2 in the m a − f plane.This is shown for the minimal model on the left in figure 4a and for the non-minimal model on the right in figure 4b.In the left panel in figure 4a , the red and yellow coloured boundaries are representing allowed flavon contribution for current values of C B d and φ B d , respectively for the minimal model based on the Z 2 × Z 5 symmetry while that for non-minimal model is shown in the right panel in figure 4b by green and purple coloured boundaries, respectively.Moreover, the region surrounded by the blue curve in the left panel shows the allowed parameter space for the observable C B d with projected limits of LHCb Phase-II for the minimal model while the same for the non-minimal model is shown by olive coloured boundary in the right panel.The bounds for the LHCb phase-I are shown by the pink coloured boundaries for the minimal as well as for the non-minimal model, and are very similar to that of the LHCb phase-II.

Figure 4 :
Figure 4: The allowed parameter space by flavour observables C B d and φ B d for λ χ = 2 in the m a − f plane for the minimal (Z 2 × Z 5 ) model in the left panel and for the non-minimal (Z 2 × Z 9 ) model in the right panel.

Figure 5 :
Figure 5: The allowed parameter space by |M D 12 | for λ χ = 2 in the m a − f plane for the minimal (Z 2 ×Z 5 ) model is shown by the magenta coloured boundary while that for the non-minimal (Z 2 ×Z 9 ) model is shown by the red coloured boundary.

Figure 9 :
Figure 9: The parameter space allowed by BR(K L → µ + µ − ) SD with λ χ = 2 in the m a − f plane for the minimal (Z 2 × Z 5 ) and non-minimal (Z 2 × Z 9 ) model are shown with red and magenta coloured boundaries, respectively.

2 × Z 9 )Figure 11 :
Figure 11: In the left panel, the red and yellow coloured boundaries represent the allowed parameter space by the current and projected sensitivities of the MEG experiment on the BR(µ → eγ) for the minimal (Z 2 × Z 5 ) model with λ χ = 2 in the m a − f plane.The same bounds for the non-minimal (Z 2 × Z 9 ) model are shown in the right panel with green and magenta coloured curves, respectively .

Figure 13 :
Figure 13: The left panel represents the allowed parameter space by BR(µ → e) conversion in the m a − f plane with different experimental limits of various target nuclei in table 9 with λ χ = 2, for the minimal (Z 2 × Z 5 ) model.Similar allowed parameter space for the non-minimal (Z 2 × Z 9 ) model is shown in the right panel.

Figure 14 :
Figure 14: The allowed parameter space by BR(µ → 3e) for the minimal (Z 2 × Z 5 ) model on the left and for the non-minimal (Z 2 × Z 9 ) model on the right with λ χ = 2.

RFigure 15 :
Figure 15: Summarized significant bounds on the minimal Z 2 × Z 5 and the non-minimal Z 2 × Z 9 model in the left and the right panel respectively.

Table 7 :
Future projected sensitivity of the neutral meson mixing.2. Phase II which is 300f b −1 LHCb and 250ab −1 Belle II (late 2030s).The expected sensitivities to C ∆m K and C Bq in future phase I and II of LHCb and Belle II can be obtained from table 7.

6 Leptonic flavour physics in the minimal and non-minimal
Z 2 × Z N flavour symmetry Charged lepton flavour violation (CLFV) has a substantial potential to provide bounds on flavon physics in future upcoming experiments.The quark flavour constraints currently dominate our flavon model.However, it is expected that the future projected sensitivities of CLFV, which are shown in table9, may significantly improve the bounds from quark flavour physics.

Table 10 :
6s −1 ] Numerical values of the dimensionless coefficients D, S p,n , V p,n and the muon capture rate for different nuclei.