Spontaneous CP violation and μ-τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu -\tau $$\end{document} symmetry in two-Higgs-doublet models with flavour conservation

In multi-Higgs-doublet models, requiring simultaneously that (i) CP violation only arises spontaneously, (ii) tree level scalar flavour changing couplings are absent and (iii) the fermion mixing matrix is CP violating, can only be achieved in a very specific manner. A general approach with new clarifying insights on the question is presented. Considering the quark sector, that peculiar possibility is not viable on phenomenological grounds. We show that, considering the lepton sector, it is highly interesting and leads to viable models with μ-τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu -\tau $$\end{document} symmetric PMNS matrices. Phenomenological implications of the models, both for Dirac and Majorana (in a type I seesaw scenario) neutrinos, are analysed.


Introduction
An appealing motivation for the introduction of several Higgs doublets in SU (2) L ⊗U (1) Y models is the possibility to have spontaneous CP violation (SCPV) [1,2].One can indeed construct viable scenarios in which the sole origin of CP violation in the CKM matrix is the vacuum [3].On the other hand, a significant source of concern in multi-Higgs models is the possibility of tree level scalar flavour changing neutral couplings (SFCNC), severely constrained by experiment.In the context of Natural Flavour Conservation (NFC) à la Glashow and Weinberg [4], reference [5] addressed the implications of spontaneous CP violation; it was argued that even if CP is broken by the vacuum, in case one requires that SFCNC are absent, the CKM mixing matrix is nevertheless CP conserving for any number of quark families (notice that, although [5] considered NFC, this argument only involved the absence of SFCNC).Later on, it was pointed out in [6] and [7] that the previous general argument can be evaded, leading to CP violating CKM matrices with very specific properties together with flavour conservation (FC), i.e. absence of SFCNC; this exception, however, is outside the NFC paradigm, where quarks of a given charge receive their mass through the couplings of precisely one neutral Higgs boson.However, as discussed in [6], this very specific form of the CKM matrix was not phenomenologically viable.In this work, we address the question in the context of the lepton sector and present a simple model where the PMNS mixing matrix, as a consequence of the simultaneous requirement of flavour conservation and CP violation originated by the vacuum, has µ−τ symmetry [8].It is to be stressed that no additional symmetry is involved in order to enforce a PMNS mixing with that property.The paper is organised as follows.In section 2 we set our notation discussing general aspects of two-Higgs-doublet models (2HDMs).In section 3, we revisit the arguments in [5], [6] and [7] concerning FC and SCPV.A general approach to the question is addressed in section 4. Finally, in section 5 we present our minimal model with Dirac neutrinos, together with a straightforward extension to a type I seesaw scenario, and discuss some relevant phenomenological consequences.Additional aspects are addressed in the appendices: appendix A is a concise summary of the most relevant points concerning the scalar sector of a CP conserving 2HDM with SCPV; appendix B is devoted to an analysis of the stability of the model (and extensions of it) under renormalization group evolution of the Yukawa matrices; and appendix C recalls the basics of the diagonalization of O(3, R) matrices.

Flavour conservation in 2HDMs
The Yukawa couplings of 3 generations of quarks and 2 Higgs doublets have the following form where Q 0 L , u 0 R , d 0 R are 3-vectors in generation space (left-handed doublets with subscript L and right-handed singlets with subscript R), Y (u) j and Y (d) j are 3 × 3 Yukawa coupling matrices, and Φ j are the scalar (Higgs) doublets, with Φj = iσ 2 Φ * j (j = 1, 2).Electroweak symmetry is spontaneously broken by the vacuum expectation values It is helpful to consider the Higgs basis {H 1 , H 2 } [9-11] with In this basis, one has where θ ≡ θ 2 − θ 1 .Bidiagonalization of the mass matrices in Eq. ( 5) proceeds as follows: since ) are hermitian and positive definite with the same set of eigenvalues, there are unitary matrices U f L and U f R such that 7) with m f j real positive, and The matrices U f L and U f R are the unitary transformations of fields into the mass bases, namely 6).The off-diagonal elements of the matrices N u and N d encode SFCNC.Therefore, there is FC if the matrices N 0 u and N 0 d are bidiagonalized with the same where, in principle, If one imposes invariance of L Yq under standard CP transformations1 [12], the Yukawa matrices are real . Similarly, CP invariance requires real coefficients in the scalar potential.Finally, if θ ̸ = 0, π, one may have SCPV [13].

SCPV and SFCNC
Considering a CP invariant Lagrangian, that is real Yukawa matrices in Eq. ( 1), and SCPV, the following "CP conserving mixing" argument was presented in [5]: FC means that the matrices M 0 q and N 0 q (q = u, d) are simultaneously bidiagonalized, which is equivalent to Y (q) 1 and Y (q) 2 being bidiagonalized at the same time.Since Y (q) 1 and Y (q) 2 are real, the bidiagonalization is achieved with real orthogonal matrices, j2 , y rephasing matrices: the CKM matrix is thus essentially real, not CP violating.
Convincing as it is, the previous argument can be evaded [6,7]: even if the matrices are real, they can have complex eigenvalues and in that case they are not necessarily bidiagonalized simultaneously with real orthogonal matrices.In [6], a model with 2 Higgs doublets, 4 quark generations and the following Yukawa matrices, was presented: ), (10) where O is a real orthogonal matrix and (the Yukawa Lagrangian L Yq for quarks in Eq. ( 1) is CP invariant).The crucial ingredient in Eq. ( 10) are the 2 × 2 blocks in Y (u) j , namely They obey (no sum over j) B j B T j = B T j B j = (a2 j + b 2 j )1 2 (henceforward, 1 n is the n × n identity matrix), so that B j B T j and B T j B j have two degenerate real eigenvalues a 2 j + b 2 j while B j has two complex conjugate eigenvalues a j ± ib j : This implies that are simultaneously diagonalized unitarily: On the other hand, it is straightforward to check that the simultaneous real orthogonal bidiagonalization of Y (u) 1 and Y (u) 2 fails.In fact, one can bidiagonalize either B 1 or B 2 with real orthogonal matrices, but not both simultaneously if B 1 and B 2 are not proportional -and that is the case if the resulting quark masses are non-degenerate-.
The diagonalization of the mass matrices involves additional (diagonal) rephasings in order to have real mass terms, omitted in the resulting CKM matrix, which is From Eq. ( 15) it follows that √ 2, and thus the rephasing invariant relation This relation is not phenomenologically viable in a 4 quark generation model; a 3 quark generation implementation of the same idea is not feasible for the same essential reason: experimentally, two CKM elements in the same row/column cannot have the same modulus.However, if instead of the CKM matrix one considers the PMNS lepton mixing matrix, this possibility is open: before addressing it in section 5, the question is analysed in more generality in the following section.
4 Evading the "CP conserving mixing" argument: a general analysis As discussed in [7], a necessary ingredient to evade the "CP conserving mixing" argument is the presence of complex conjugate eigenvalues in the Yukawa matrices.In a fermion sector (f) with CP invariance at the Lagrangian level, this implies that the matrices have real degenerate eigenvalues, and thus the simultaneous real orthogonal bidiagonalization involved in the "CP conserving mixing" argument fails; nevertheless, a simultaneous unitary diagonalization of the Yukawa matrices which involves complex eigenvectors could exist and lead to a CP violating mixing matrix.
We address the question aiming to a more transparent understanding of the conclusions reached in [7], clarifying in the process how the simultaneous unitary diagonalization of the Yukawa matrices arises.We circumscribe our discussion to 3 fermion generations.Since the key ingredient is the degeneracy of the eigenvalues of we discuss separately the case of two degenerate eigenvalues and the case of all three eigenvalues degenerate.

Two degenerate eigenvalues Let us consider first that
is real and symmetric, there exists a 3 This does not fix O L completely, since Eq. ( 17) is invariant under If one writes explicitly and thus Comparing Eq. ( 17) with Eq. ( 21), and assuming µ j1 , µ j2 ̸ = 0, one obtains ⃗ w ′T j = (0, 0) and, most importantly, has two degenerate eigenvalues is thus of the form with The appearance of an orthogonal submatrix is most important.On the one hand, those having det They have complex conjugate eigenvalues e iα and e −iα , with corresponding complex conjugate eigenvectors which is the central ingredient of the counterexample presented in [6] to the "CP conserving mixing" argument discussed in section 3: it is sufficient to realize that, introducing λ j , then a j ± ib j = λ j e ±iφ j and the B j matrix in Eq. ( 11) is simply Notice in addition that U 0 in Eq. ( 24) does not depend on α: all SO(2, R) matrices are simultaneously diagonalizable following Eq.( 24).On the other hand, if the orthogonal submatrix of interest satisfies det O 2 = −1 instead, it has the form and presents real eigenvalues ±1: although the degeneracy in is still present, there is no simultaneous irreducibly complex unitary diagonalization such as Eq. ( 24).

Three degenerate eigenvalues
Having stressed the role of orthogonal matrices, the case of having all eigenvalues degenerate becomes almost straightforward: one necessarily has matrices have two complex conjugate eigenvalues with corresponding complex conjugate eigenvectors, and a real ±1 eigenvalue corresponding to a real eigenvector.We emphasize two important points (for the sake of clarity, a short reminder of the unitary diagonalization of O(3, R) matrices is provided in appendix C): (i) the eigenvectors are orthonormal and thus a unitary diagonalization of Y (f) j exists, (ii) the eigenvectors do not depend on the eigenvalues and thus there are different O(3, R) matrices that can be unitarily diagonalized simultaneously.
One can then construct counterexamples to the "CP conserving mixing" argument extending the case analysed in [6] in a transparent manner via , with all O j simultaneously diagonalizable unitarily.In [6,7] there is no mention of the intimate connection between degeneracies, complex eigenvalues, a simultaneous unitary diagonalization, and the role of orthogonal matrices: the previous discussion makes it clear (without invalidating, of course, the results in these references).
Concerning the resulting mixing matrices, the crucial common point in the previous two cases (Y with two or three degenerate eigenvalues) is the following.Consider the unitary matrix U that diagonalizes simultaneously the Y (f) j for all j; its columns, as vectors, are {⃗ r, ⃗ v + , ⃗ v − }, which are orthonormal and, crucially, | as in section 3 (starting with a different ordering of the columns of U, the equality of moduli of V elements would apply to the correspondingly permuted rows).If the Y (f) j are such that one has instead V = O T 3 U, the equality of moduli is in terms of columns rather than in terms of rows.One can even consider a mixing matrix V = U † U ′ with both U and U ′ complex unitary matrices with orthonormal columns In this case, rather than equality of moduli in two full different rows/columns, one has − : such a possibility is not phenomenologically viable neither for CKM, nor for PMNS; besides the equality of moduli, the resulting mixing matrix is CP conserving [7] (one can readily check that the imaginary part of a rephasing invariant quartet vanishes).It is thus clear how, in general, the basic ingredient underlying the particular example of section 3 can yield a mixing matrix with the same kind of specific property: equality of moduli in pairs of rows/columns.As mentioned, this kind of property is ruled out in the quark sector, but not in the lepton sector, to which we now turn our attention.

Lepton models
This kind of scenario with FC and a spontaneous origin of CP violation giving an irreducibly complex mixing matrix is not viable in the quark sector.Consequently, we do not discuss the quark sector further and only mention that one would need, in order to have SCPV and a complex CKM, to ensure that SFCNC are under control, as done for example in [3].Nevertheless, the idea can still work in the lepton sector.In subsection 5.1 we present a simple implementation with Dirac neutrinos; then, in subsection 5.2, we extend it to a seesaw scenario with Majorana neutrinos, and finally, in subsection 5.3 we discuss some relevant phenomenological consequences concerning neutrinos and charged leptons interactions.

Dirac neutrinos
For Dirac neutrinos we introduce 3 right-handed neutrino fields ν 0 R , which are singlets under the gauge group, giving rise to the following leptonic Yukawa sector: where the Yukawa matrices are given by with λ j , φ j , y jk real and O a real orthogonal 3 × 3 matrix.Thus, L Yℓ is CP invariant.The procedure is analogous to the steps recalled for the quark sector in section 2: transformation of the scalars into the Higgs basis, identification of mass matrices and additional Yukawa matrices and diagonalization of mass matrices; with the first step, Electroweak and CP invariance are spontaneously broken by the vacuum in Eq. ( 2).Then, the mass matrices read while the additional Yukawa couplings are The diagonalization of the charged lepton mass matrix M 0 ℓ is implemented as where and the rephasing R ℓ R to obtain real diagonal elements is5 Although we analyse N ℓ in subsection 5.3, the important point to notice meanwhile is that N ℓ is diagonal -there is FC as intended-but not necessarily real, and this has potential implications for CP violation to be addressed later.The neutrino mass matrix M 0 ν , and also N 0 ν , are already diagonal -only a rephasing is needed to have real eigenvalues in M ν = U † ν L M 0 ν U ν R -and thus the PMNS matrix (up to rephasings) is Then, leading to the rephasing invariant property The PMNS matrix in Eq. (37) has "µ − τ symmetry" [8].One salient consequence of this "µ − τ symmetry" is that, in the standard PDG parametrization [14,15], the angle θ 23 = π/4 and the Dirac CP violating phase is δ = ±π/2.Furthermore, with U 1j = O 1j , the parametric freedom in O allows to fix the values of the remaining parameters in agreement with experiment [16,17].Some comments are in order.
• It is remarkable, and to some extent counterintuitive, that CP violation in PMNS is independent of the specific value of θ ̸ = 0, π: that is, no matter the specific value of θ (as long as θ ̸ = 0, π), the Dirac CP phase in the PMNS matrix is δ = ±π/2.
• On general grounds, if two charged leptons were degenerate, the PMNS matrix would be CP conserving: following Eq.(32), one thus needs, in terms of the parameters of the model, In addition to having two relatively complex vevs, i.e. c β s β ̸ = 0 and θ ̸ = 0, π, nondegeneracy of the charged leptons requires the presence of 2 × 2 Eq. ( 25)-like blocks in both and 2 , and furthermore these blocks have to be not proportional.• Notice the peculiar fact that in the "basic" Jarlskog CP-violation sensitive invariant [11,18,19], non-fulfilment of Eq. ( 38) would yield det[M ℓ M † ℓ , M ν M † ν ] = 0 because of the (squared) mass differences factors.In this model, checking CP violation at the level of weak basis invariants is crucial [20].
• Eq. ( 28) is the only viable implementation of the basic idea to obtain a complex PMNS matrix with SCPV together with flavour conserving j : relations like Eq. (36) for other rows, or instead for columns if one exchanges j , are in conflict with with the current knowledge of the PMNS matrix [16,17].One can also consider couplings like Y (ℓ) j in Eq. ( 28) both for neutrinos and charged leptons: in that case, equality of moduli in two rows/columns of the PMNS matrix would not extend to the whole row/column, but that would nevertheless remain not viable (and would lead to a CP conserving PMNS matrix, contrary to the intended goal).
• It is to be noticed that the diagonalization of the neutrino mass matrix has had almost no role in the obtention of the µ − τ symmetric PMNS matrix (we have indeed chosen diagonal Y to start with).
• It is to be stressed that the µ−τ symmetry of the PMNS matrix does not derive from an additional fundamental symmetry of the model; as a consequence, as discussed in appendix B, the minimal implementation of the basic idea -to combine SCPV with FC and CP violating mixings-, is not stable under one loop renormalization group evolution (one loop stability can be nevertheless obtained within non-minimal implementations).
We have thus achieved the main goal of this work: Eq. ( 28), together with SCPV, gives a model without tree level SFCNC where the PMNS matrix is not only CP violating, it is phenomenologically viable and has µ − τ symmetry.In appendix B we discuss the stability of this scenario under renormalization group evolution (RGE).

Majorana neutrinos -type I seesaw
Attending to the fact commented above that the diagonalization of the neutrino mass matrix has a secondary role in the previous scenario with Dirac neutrinos, it is natural and simple to extend this scenario with the inclusion of Majorana mass terms for the right-handed neutrinos (and thus incorporate, as usual, a rationale for the smallness of neutrino masses).Aiming for simplicity, we introduce diagonal Majorana mass terms

Neutrinos
In the Dirac neutrino case, both N 0 ν and N ν are diagonal and there are no SFCNC.Furthermore, if M 0 ν in Eq. (30) does not involve large cancellations, one reasonably has (N ν ) jj ∼ m ν j and thus the flavour conserving Yukawa interactions of neutrinos appear to be safely negligible with a mν j v ∼ 10 −13 suppression.Concerning the Majorana neutrinos of the seesaw scenario, there are some differences with respect to the Dirac neutrino case.First, since the Yukawa interactions arise from L 0 L N 0 ν H2 ν 0 R + H.c. terms in Eq. ( 29), and ν 0 L (ν 0 R ) are mostly light (heavy) neutrinos, the Yukawa couplings involving only light neutrinos are suppressed by sin α j in Eq. ( 44).In addition, since one has a scenario with three independent pairs with one light and one heavy neutrino each (mixings and couplings between different pairs are absent), there are in fact no SFCNC in the light-light sector.There are, nevertheless, SFCNC involving one light and one heavy neutrino (one non-vanishing coupling for each independent pair): since our concern lies with effects not at the large scale M R j , we do not discuss them further.
As a final comment on the neutrino sector, while the original motivation -obtaining a CP violating PMNS mixing matrix when one requires that (i) the only source of CP violation is the vacuum state, and that (ii) there are no SFCNC-fully applies in the Dirac neutrino case, this motivation is weakened in the second scenario with Majorana neutrinos, since neutrino masses arise differently (they involve a different diagonalization).However, it is remarkable that the most relevant features of the Dirac neutrino case in subsection 5.1, that is a µ − τ symmetric PMNS and no SFCNC for charged leptons, are preserved in the simple seesaw Majorana scenario of subsection 5.2.

Charged leptons
The most important sources of potential phenomenological concern are the flavour conserving couplings of the charged leptons in N ℓ .From Eqs. (30)-(32), with φ = φ 2 − φ 1 , we have (48) The charged lepton masses appear explicitely in Eq. (48) because, from Eq. (32), we have used 2 The couplings of the charged leptons with the physical neutral scalars (see appendix A for further details) are with j = 1, 2, 3, corresponding to {S 1 , S 2 , S 3 } = {h, H, A}, and where R jk is the scalar mixing in Eq. ( 62).
The observables that can be affected and deserve some attention are listed below.
• Similarly, scalar resonances in pp → S → ℓ + ℓ − for q 2 ∼ m 2 S (with q 2 the invariant mass of the lepton pair); one can also consider production of the charged scalar H ± and decays H + → ℓ + ν, H − → ℓ − ν, but they are controlled by the same couplings and the experimental sensitivity is poorer.
• Contributions to precision observables sensitive to loops involving the new scalars and the couplings in Eq. ( 48), such as the g − 2 of the electron and the muon (g is the gyromagnetic ratio), and the electric dipole moment of the electron (eEDM).
One could address them through a full numerical exploration of the parameter space of the model (which in principle has enough parametric freedom to satisfy those constraints), but that is beyond the scope of this work -establishing that the model can be viableand would require, in addition, that the quark sector is also specified.Fortunately, references [21][22][23] can provide a satisfactory answer.In these works, a general flavour conserving 2HDM is considered in order to address deviations from SM expectations in g − 2 of both the electron and the muon.In the full numerical analyses presented in [21][22][23], signal strengths of the SM Higgs-like and bounds on resonant production of the new scalars are relevant constraints.These constraints are satisfied despite the fact that the diagonal entries of N ℓ required to obtain sizable contributions to (g − 2) e,µ are much larger than the corresponding charged lepton masses (in addition, one important ingredient is that the scalar sector has to be close to the scalar alignment limit).Since here we do not necessarily require such values of the couplings, no problem is to be expected from these observables.The only exception is the eEDM d e , since CP conservation is assumed in [21][22][23] that is As a guidance to how problematic the constraint could be, consider the contribution of a scalar with mass m S ≃ 500 GeV and couplings a e S = b e S ≃ 10m e ; with these values, It is thus clear that, even considering a scalar S not too heavy together with enhancements in both the scalar and pseudoscalar couplings a e S and b e S with respect to m e , these contributions are not expected to be in conflict with the experimental constraint.One might worry then about Barr-Zee two loop contributions, where the key point is that one of the "small" electron Yukawa couplings is replaced by the Yukawa coupling of an internal fermion (which can compensate other loop suppression factors).The contribution of the illustrative diagram in Figure 1b to d e is where N f c and Q f are the number of colours and the electric charge of fermion f , a S f and b S f are the scalar and pseudoscalar couplings of f (as in Eq. ( 50 | ∼ 3 × 10 −13 e • GeV −1 , which is incompatible with the experimental bound.From this rough comparison it is clear that the eEDM is indeed a relevant constraint and these illustrative values (considered simultaneously) are too naive or generic: although not "automatically granted", compatibility with experiment can be regained by a combination of heavier scalars, smaller Yukawa couplings of the electron, 6 smaller Yukawa couplings of the top quark (unless the bottom quark contribution is much enhanced by the couplings, the dominant contribution is the top quark one), and even cancellations from the contributions of the different scalars.As an example consider the joint contribution of two scalars S 1 and S 2 with masses m S 1 = m S 2 = 1 TeV, couplings to the electron given by Eq. ( 51) with n e = e iπ/ 4  15 m e , couplings to the top quark analogous to Eq. ( 51) with ℓ → t, R 3j → −R 3j (this change of sign is simply due to the opposite weak isospin of up quarks and charged leptons) and n t = e iπ/ 4  15 m t , and finally scalar mixings  One might also worry about potentially dangerous contributions to µ → eγ transitions similar to Figure 1.At one loop, rather than the diagram 1a with a neutral scalar and an electron running in the loop (with the external photon attached to the electron) one has an internal charged scalar and a neutrino (the external photon now attached to the scalar).At two loops, diagram 1b is modified replacing the neutral scalars S with the charged scalar H ± and the internal photon with a W ∓ .For Dirac neutrinos, the smallness of neutrino masses and couplings (N ν ) jj ∼ m ν j safely suppress µ → eγ transitions.For Majorana neutrinos in the scenario of section 5.2, one might worry about the contributions with internal "mostly heavy" ν 0 R .One expects again suppressed contributions from the large heavy neutrino masses in the one loop diagrams, and (less transparently, see [24]) from the heavy neutrino in the initial-to-final fermion line, together with the heavier W ∓ in the internal vector boson line.In any case, although plausible, this does not guarantee automatic agreement with experimental constraints on µ → eγ transitions, and, as in the case of the electron EDM, they can play a relevant role in shaping the available parameter space in more detailed phenomenological studies.

Conclusions
In the context of multi-Higgs-doublet models where CP invariance holds at the Lagrangian level, the simultaneous requirement of (i) absence of SFCNC and (ii) spontaneous CP violation, is analysed.In [5], it was argued that these two conditions necessarily lead to a CP conserving mixing matrix.As shown in [6], the Yukawa couplings of a very particular model can reconcile both requirements with a CP violating mixing matrix.Further aspects of the question were addressed in [7].There are two central ingredients.First, the combinations must have (real) degenerate eigenvalues.Second, Y (f) j must have complex conjugate eigenvalues with corresponding complex conjugate eigenvectors leading to a simultaneous unitary diagonalization of all Y (f) j (j = 1, 2 . ..), which guarantees the absence of SFCNC.A general approach to the model presented in [6] and the analysis presented in [7] has been addressed in this work, stressing the important role of orthogonal matrices and their intimate connection with degeneracies, complex eigenvalues and a simultaneous unitary diagonalization of the Yukawa couplings.We present the only viable implementation of this idea in the lepton sector, considering both Dirac and Majorana neutrinos (the latter in a type I seesaw scenario).In all cases, the resulting PMNS mixing matrix, as a consequence of the simultaneous requirement of SFCNC absence and SCPV, is experimentally viable and µ − τ symmetric.Phenomenological implications of the model are discussed, with particular emphasis on the role of the electric dipole moment of the electron in shaping the available parameter space.one can identify the would-be Goldstone bosons G 0 , G ± , and read the mass terms for the new neutral scalars (for H ± one can readily obtain m 2 with M 2 0 real and symmetric.The neutral scalars mass matrix M 2 0 is diagonalized with a real 3 which gives the physical neutral scalars {h, H, A}.The important point to recall here is that, since are in principle non-zero, R "mixes" all neutral scalars and we do not have, as usual in other contexts, scalars h, H and pseudoscalar A (as clearly seen Eqs. ( 50)-( 51)).

B Renormalization group evolution
As commented in precedence, the µ − τ symmetric character of the PMNS matrix does not derive from an imposed additional symmetry.It is thus interesting to analyse the one loop renormalization group evolution of this kind of scenario, paying special attention to two aspects: (i) despite being absent at tree level, how SFCNC may arise at one loop; (ii) is it possible to depart from the minimal implementation of the model (while maintaining the essential ingredients) and still obtain absence of SFCNC at one loop?
The one loop renormalization group evolution (RGE) of the lepton Yukawa couplings is It is rather clear from Eqs. ( 64)-(65) that in principle SFCNC will arise through one loop RGE.If FC is imposed at some scale µ FC , one would thus expect SFCNC at a scale µ proportional to ln(µ/µ FC ), which might be sufficient to keep them under control.In both Eqs. ( 64) and (65), the terms that can spoil FC are the ones with products of 3 Yukawa coupling matrices within square brackets.Half of these terms are nevertheless harmless.With Y a , are diagonal.From the point of view of FC, the problematic terms which give rise to SFCNC through RGE are Y a .If one thinks of extending the simple model of section 5.1 in order to preserve FC at one loop, it is by no means impossible to consider scenarios where this is achieved, as we discuss for illustration in the following.First, one can observe that if more than two Higgs doublets are considered and if no Higgs doublet couples simultaneously to ℓ 0 R and ν 0 R , then all terms Y where, for convenience, a parametrization of the real unit vector (n 1 , n2 , n3 ) with sphericallike coordinates is used.The set {⃗ v 0 , ⃗ v + , ⃗ v − } is orthonormal, 9 and thus we have the unitary diagonalization For O = exp(αA) it is immediate that the eigenvalues are 1, e iα and e −iα , the corresponding eigenvectors are ⃗ v 0 , ⃗ v + and ⃗ v − above, and its unitary diagonalization reads This construction makes clear that the eigenvectors of O do not depend on α; this implies, in particular, that O 1 = exp(α 1 A) and O 2 = exp(α 2 A) with α 1 ̸ = α 2 and the same A in Eq. (70), can be unitarily diagonalized simultaneously.

8π 2 m e v 2 j=h
. Let us analyse the eEDM in some detail.The experimental constraint is |d e | < 1.1 × 10 −29 e • cm = 5.6 × 10 −16 e • GeV −1 .With the couplings in Eq. (51), the contribution of one loop diagrams like Figure 1a to d e is |d e | = e )), and F and G are loop functions.With the previous values m S ≃ 500 GeV, a e S = b e S ≃ 10m e , and an internal top quark f = t with couplings a t S = b t S ≃ m t , we have |d (S,t) e

S e γ 5 (
sin θ 23 with θ 23 = 0.2.Considering these illustrative values, which have similar size scalar and pseudoscalar couplings both for e and t, one has |d (S 1 ,t) e +d (S 2 ,t) e | ≃ 3.7×10 −16 e•GeV −1 , in agreement with the experimental bound.The main point, to close the discussion, is that the model remains viable but the eEDM plays an important role in shaping the available parameter space.b) Two loops (Barr-Zee).
where c ℓ and c ν are linear combinations of squared SU (3) c , SU (2) L and U (1) Y gauge couplings, D ≡ 16π 2 d d ln µ with µ the energy scale, and n d is the number of Higgs doublets.