The framework for a common origin of $\delta_{\rm CKM}$ and $\delta_{\rm PMNS}$

We analyse a possible connection between CP violations in the quark and lepton sectors, parametrised by the CKM and PMNS phases. If one assumes that CP breaking arises from complex Yukawa couplings, both in the quark and lepton sectors, the above connection is not possible in general, since Yukawa couplings in the two sectors have independent flavour structures. We show that both the CKM and PMNS phases can instead be generated by a vacuum phase in a class of two Higgs doublet models, and in this case a connection may be established. This scenario requires the presence of scalar FCNC at tree level, both in the quark and lepton sectors. The appearance of these FCNC is an obstacle and a blessing. An obstacle since one has to analyse which models are able to conform to the strict experimental limits on FCNC, both in the quark and lepton sectors. A blessing, because this class of models is falsifiable since FCNC arise at a level which can be probed experimentally in the near future, specially in the processes $h\to e^\pm\tau^\mp$ and $t\to h c$. The connection between CP violations in CKM and PMNS is explicitely illustrated in models with Minimal Flavour Violation.


Introduction
In this paper we study the possibility of having a relation between CP violation in the quark and lepton sectors, parametrized by the phases δ CKM and δ PMNS , respectively, of the CKM [1,2] and PMNS [3,4] mixing matrices. This question is specially important in view of the efforts to detect leptonic CP violation in neutrino experiments such as Dune, T2K and Nova [5,6,7]. At present, there is solid experimental evidence that CKM is complex, even if one allows for the presence of New Physics contributing to CP violation [8,9]. The fact that CKM is complex does not imply that CP violation is violated at the Lagrangian level through complex Yukawa couplings. Indeed it has been pointed out [10] that one may have a vacuum induced CP violation generating a complex CKM matrix in agreement with experiment. If one considers an extension of the Standard Model (SM) with non-vanishing neutrino masses and assumes that the origin of CP violation is the presence of complex Yukawa couplings, then there is no relation between δ CKM and δ PMNS . This just reflects the fact that in the SM the lepton and the quark Yukawa couplings are completely independent quantities. If one wants to obtain a relation between δ CKM and δ PMNS , an interesting possibility is to assume that CP is spontaneously broken, with the vacuum phase generating both the phase in the quark and in the lepton sectors. In this paper, we consider the extension to the leptonic sector of a previously proposed viable minimal model where the Lagrangian respects CP invariance but the vacuum is CP violating, with a complex phase which generates a complex CKM matrix. The model consists of a generalised Branco-Grimus-Lavoura (BGL) model [11,12,13,14] with a flavoured Z 2 symmetry in the context of two Higgs doublet models (2HDMs) [15,16]. Two of the quark families are even under the Z 2 symmetry while the remaining one is odd. It was shown [10] that there is a profound connection between the possibility of generating a complex CKM matrix and the existence of tree level Scalar Flavour Changing Neutral Couplings (SFCNC) both in the up and down quark sectors. The same correlation also occurs when one extends the model to the leptonic sector. In order for δ PMNS = 0 to be generated one has to have SFCNC both in the charged lepton and neutrino sectors. The fact that SFCNC necessarily arise in the class of models which we are considering is both an obstacle and a blessing. It is an obstacle because one has to analyse which models are able to conform to the strict experimental limits on SFCNC in the quark and lepton sectors. It is a blessing because the class of models we are considering are falsifiable in the sense that they imply SFCNC at a level which can be probed at different experiments. The paper is organised as follows. In the next section we present the model, including the quark and lepton sectors. In section 3 we analyse how δ CKM and δ PMNS are generated. In section 4 we analyse how CP violation requires the presence of SFCNC. Section 5 is dedicated to a study of a general relation between δ CKM and δ PMNS . In section 6 we present simplified models incorporating Minimal Flavour Violation (MFV) [17,18] and the connection to SFCNC; a specific model is analysed in detail, including the prediction of δ PMNS arising in this subclass of models. Finally in the last section, we present our conclusions. Additional details on different aspects are provided in the appendices.

The model
The Yukawa sector of the most general 2HDM, in the case of Dirac neutrinos, can be written as We extend the viable model presented in reference [10] to the leptonic sector. The model is enforced by a Z 2 symmetry, extended now to the leptonic sector. The fields Φ 2 , Q 0 L 3 and L 0 L 3 are odd under Z 2 , the rest of the fields are even under Z 2 . This symmetry gives rise to the so-called generalized BGL (gBGL) textures [14] for the original Yukawa coupling matrices Γ with × denoting generic entries.
Note that Γ where P 3 is the projector Spontaneous electroweak symmetry breaking is achieved with the vacuum As usual, we introduce v 2 ≡ v 2 1 + v 2 2 and β such that v 1 = vc β , v 2 = vs β (here and in the following, c β = cos β, s β = sin β, t β = tan β, t −1 are introduced; {H 1 , H 2 } is the "Higgs basis" [19,20], that is One can rewrite the Yukawa couplings in eq. (1) as Note that even if N 0 f are proportional to M 0 f , this proportionality involves a diagonal matrix different from the identity. This means that in general it will not be possible to bi-diagonalize both matrices simultaneously. The matrices N 0 f control the scalar mediated flavour changing neutral couplings. Concerning the scalar sector, the scalar potential has the following form with µ 2 jk ∈ R and λ k ∈ R. For a scalar potential with an exact Z 2 symmetry, µ 2 12 = 0 and even if λ 5 could be complex, there is no CP violation. In eq. (13), µ 2 12 = 0 softly breaks the Z 2 symmetry and even though µ 2 12 , λ 5 ∈ R give a CP invariant V (Φ 1 , Φ 2 ), the possibility to have CP violation arising from the vacuum is open. The expansion of the fields around the vacuum in the Higgs basis reads The neutral mass eigenstates are {h, H, A}: with R a real orthogonal 3 × 3 matrix. h is usually assumed to be the SM-like Higgs with mass m h = 125 GeV; the alignment limit, in which its couplings are SM-like, corresponds to R 11 → 1. For further details concerning the scalar sector we refer to [10].

Generation of CP violating CKM and PMNS matrices
It is clear that the global phase θ 1 in the previous expressions can be rephased away by redefining the phases of the H i fields or of the right handed fermion fields, and thus we set θ 1 = 0 from now on without loss of generality. Note also that invariance under CP of the entire Lagrangian implies that As shown in reference [10], taking into account the reality condition in eq. (16) and the position of the irremovable phase θ in the textures in eqs.
(2)-(3), the mass matrices can be factorized as where M 0 f are arbitrary real mass matrices and σ d = σ e = θ, σ u = σ ν = −θ. Note that the diagonal matrix of phases ϕ 3 (σ f ) can be written as Correspondingly, the CKM mixing matrix V = U † u L U d L and the PMNS mixing matrix U = U † e L U ν L , defined with the usual conventions in the charged currents W ± interactions, are: Since O f L are arbitrary real rotations, it is evident that there is enough freedom in eq. (24) to obtain arbitrary V and U , except for the fact that any CP violating observable in the quark sector and any CP violating observable in the lepton sector, must vanish with θ → 0. It is thus interesting to scrutinize in detail the relation that must exist among the CP violating phases in V and U , δ CKM and δ PMNS respectively. Anticipating the discussion in section 5, δ CKM and δ PMNS will simply correspond to the CP phases in a standard parametrization; notice that, a priori, the change of sign in θ entering V and U in eq. (24), does not imply in general that δ CKM = −δ PMNS .

CP violation and the presence of SFCNC
It was shown in [10] that in this class of 2HDMs with spontaneous CP violation, there is a deep connection between the complexity of the CKM matrix and the presence of SFCNC. Since there is no evidence yet of SFCNC beyond the SM, the simplest approach in the analysis of these models would be to impose that SFCNC are absent. As discussed in [10], this leads to a real CKM, contrary to evidence, and thus SFCNC are necessary. We recall here the essence of this connection. The appearance of SFCNC is encoded in the N 0 f matrices in eq. (11), which control the Yukawa couplings of H 2 ; in the fermion mass bases, N 0 f → N f , and eq. (12) gives where we have introduced the projection operators In eqs. (25)- (26), SFCNC are controlled by the real projectors P f 3 , in particular the off-diagonal entries of P f 3 , which are controlled by the O f L matrices, which also give the CKM and PMNS mixing matrices. It is important to notice that, by construction, Equation (27) means that SFCNC in the up and down quark sectors are not independent, they are related through the CKM matrix. For example, if one fixes SFCNC in the up quark sector, SFCNC in the down quark sector are completely determined; this fact will be particularly relevant in order to address appropriately the count and the election of the independent parameters in the model. The situation in the lepton sector is analogous. The elements of the matrices P f 3 are for a fixed k, i.e. P f 3 = P k for that given f . Consider, for example, the absence of SFCNC in the neutrino sector, that is P ν 3 = P k . Using eq. (18) in eq. (24), Then the PMNS matrix U is written as a real rotation times a diagonal matrix of phases (with e −2iθ in position k and the rest of them 1), and thus there is no CP violation. Similarly, if one starts with the absence of SFCNC in the charged lepton sector, P e 3 = P k and giving again a CP conserving mixing matrix. Therefore, in this model, in order to have a non-vanishing CP violating phase in the CKM matrix, there must be tree level SFCNC both in the up and in the down quark sectors and, mutatis mutandis, in order to have a non-vanishing CP violating phase in the PMNS matrix, there must be tree level SFCNC both in the neutrino and in the charged lepton sectors.

The general relation between δ CKM and δ PMNS
From the discussion in the previous sections, it is important to recall that The CKM and PMNS matrices can be parametrized, up to rephasing of fields, in terms of 4 quantities each; in the usual PDG parametrizations [21], the parameters are {θ q 12 , θ q 13 , θ q 23 , δ q } and {θ 12 , θ 13 , θ 23 , δ } respectively. In the quark sector, experimental information allows the extraction of θ q 12 , θ q 13 , θ q 23 , and of the CP violating phase δ q , which is neatly different from zero. In the lepton sector, experimental information allows to extract θ 12 , θ 13 , θ 23 , but the phase δ remains 'the last frontier', where some sensitivity is emerging in current analyses [22,23], but still far from a neat determination. Although the form of the CKM and PMNS matrices in eq. (24) is different from the PDG parametrization, one can impose in an invariant manner that V and U agree with the experimental information encoded in {θ q 12 , θ q 13 , θ q 23 , δ q } and {θ 12 , θ 13 , θ 23 , δ } (see appendix B for a detailed explanation). In particular, δ CKM , the CP violating phase in V , is the model prediction for δ q and similarly δ PMNS , the CP violating phase in U , is the model prediction for δ . We already know that if θ = 0, then δ CKM = δ PMNS = 0, but θ = 0 δ CKM = 0, δ PMNS = 0. At this point, we need to discuss the independent parameters in the model. We start with the quark sector, keeping in mind that an important goal is to analyse how information on CP violation in the quark sector, i.e. the constraint δ q = δ CKM , can translate into some prediction on δ PMNS . The CKM matrix in eq. (24) involves 7 real parameters: 3 in the rotation O u L , 3 in the rotation O d L , and θ. It should match the four independent quantities measured in the CKM matrix, equivalent to {θ q 12 , θ q 13 , θ q 23 , δ q }. Besides the CKM matrix, the (28)) also control the SFCNC which, we recall, must be present: they involve 2 of the parameters in O u L , O d L . This means that measurements of CKM and SFCNC processes (e.g. h → uc,ūc, t → hc, t → hu, h → bs,bs, etc), that is measurements of {θ q 12 , θ q 13 , θ q 23 , δ q } and of {r [u]1 ,r [u]2 } (or, equivalently, of 2 independentr [u]j ,r [d]k ), can provide sufficient constraints to fix the parameters of the model. Although there are in principle 7 parameters, one can reduce this number to 6 in a simple manner, as we now discuss. Each rotation can be written as the product of three two-dimensional rotations controlled by one parameter each, where each R j (x) can be one of the following (with and c x = cos x, s c = sin x.
If one chooses a parametrization with the leftmost rotation only acts in the 1-2 plane), for example it is clear (note the ϕ 3 (2θ) diagonal structure in eqs. (17)- (18)) that, rather than p u 1 and p d 1 separately, only p u 1 − p d 1 enters V and we can eliminate one parameter at once (we simply set p d 1 = 0 without loss of generality). In summary, the experimental information constrains {θ q 12 , θ q 13 , θ q 23 , δ q ,r [u]1 ,r [u]2 }, and could fix the model parame- θ} (a full analysis along these lines was presented in [10]). The most important aspect here, in which we are interested, is the fact that, ideally, one can fix θ with this procedure, since CP violation is well established in the quark sector. One can address the lepton sector similarly with Again, one can set p ν 1 = 0 without loss of generality, ending up with {p e 1 , p e 2 , p e 3 , p ν 2 , p ν 3 , θ} as parameters in the lepton flavour sector. The experimental information on PMNS strongly constrains {θ 12 , θ 13 , θ 23 }; additional information from SFCNC sensitive processes like h → i¯ j , i = j, is needed in order to constrain or fix the parameters in eq. (37). The crucial point is that θ can be a priori fixed in the quark sector and thus, with one less experimental input in the leptonic sector, one could in principle predict the value of the CP violating Dirac phase δ in PMNS prior to its measurement. It is in this sense that we can ideally relate the PMNS phase to the CKM phase in this class of generalized BGL-2HDM with spontaneous CP violation.

Simplified models incorporating MFV and their connections to SFCNC
In the full analysis of the quark sector presented in [10], it was shown that the model was viable after surmounting a large set of constraints related to flavour transitions, Higgs signals, electroweak precision and requirements on the scalar potential. In particular, the right amount of CP violation in the CKM matrix could be accommodated, together with the presence of SFCNC. Surprisingly, this could be achieved with significant freedom in the values of θ and SFCNC. Therefore, direct generalization of the full analysis to include the lepton flavour sector does not appear to be the most promising avenue to explore the connection among CP violation in the CKM and PMNS mixing matrices in this kind of model, since this connection is blurred by this remaining freedom in θ and SFCNC. An interesting possibility is to restrict the model by making simplifying assumptions about SFCNC, with these assumptions guided by -and thus compatible with -experimental data. In the following we focus on that kind of restricted scenario. We will first address the quark sector and then discuss the extension to the lepton sector.

Quark sector
As explained in section 4, eliminating SFCNC either in the up or in the down quark sector, definitely simplifies SFCNC but eliminates altogether a CP violating CKM matrix. The next level of simplification would be to impose the absence of some SFCNC. Considering the general structure of SFCNC, proportional tor [f ]jr[f ]k , the only simplified alternative between the most general possibility and the absence of SFCNC is to assume only one vanishing componentr [f ]i = 0 for a given i = 1, 2 or 3; in that case only one SFCNC transition j k is present (i, j, k all different). Notice that, a priori, one can make this kind of assumption simultaneously in both the up and the down sectors: that is, for given i and l, one can in principle imposê r [u]i = 0 andr [d]l = 0. According to the discussion in section 5, we have 6 parameters that need to satisfy the 4 constraints to obtain a realistic CKM matrix, to which we are now adding these 2 new requirements on SFCNC. One can indeed implement these 2 SFCNC conditions directly with a reduction from 6 to 4 parameters in our CKM matrix in eq. (34). Notice, in that case, that the models incorporate the MFV ansatz [17,18], since they have exactly as many parameters in the flavour sector as the CKM matrix requires. Furthermore, the only non-vanishing SFCNC coupling in each sector will be controlled by one of the mixing angles of CKM. As discussed in appendix C, one can consider, in principle, 81 different models of this type, out of which only one appears to satisfy the requirements to be viable: this model is the specific example that we consider now, step by step, in order to illustrate the central idea of the paper.
1. We start imposing the following form of ther [f ] vectors controlling SFCNC: with parameters p u 2 , p d 2 .
The assumption in eq. (40) is that t c and b d SFCNC are present while u c, u t, d s, b s SFCNC are absent.

Concerning CKM, sincer
One can add a factor R 12 (p d 1 ) to the left of R 13 (p d 2 ) in O d L but, as discussed previously, this would just amount to a redefinition p u The main point is that the CKM matrix is In order to relate δ CKM and δ PMNS it is specially relevant that the quark sector fixes θ. with and

Lepton sector
We address now the lepton sector, applying analogous requirements on SFCNC to the ones considered for the quark sector in the previous subsection. In the lepton sector, the most stringent constraint comes from bounds on µ → e + γ. If we only allowed µ e SFCNC, the coupling would be controlled by |U ei U µi | 2 with i = 1, 2 or 3. Since PMNS is rather non-hierarchical, one can estimate [26] that avoiding the current bound Br(µ → e + γ) < 4.2 × 10 −13 [27] requires a cancellation or finetuning at the 10 −4 -10 −5 level among scalar and pseudoscalar contributions in 2-loop Barr-Zee contributions [28,29].
1. It is mandatory to eliminate µ e SFCNC: this can be achieved either witĥ r [e]1 = 0 orr [e]2 = 0. In the neutrino sector all three choicesr [ν]k = 0, k = 1, 2, 3 are in principle available. However, out of these 6 combined possibilities which automatically evade µ → e + γ bounds, the only one which is allowed experimentally (see appendix C for details) iŝ and where the permutation P 23 interchanges the third and second rows of R 12 (p ν 2 ): (Including P 23 does not depart from the general form in eq. (32); as explained in appendix D). As discussed previously, we do not include a left factor R 12 (p ν 1 ) in O ν L , since it amounts to a redefinition of p e 1 → p e 1 − p ν 1 .

The PMNS matrix is
It is fully fixed by 3 mixing angles and the CP violating phase θ already obtained in the quark sector.
4. We can fit now eq. (53) to the experimental information on PMNS encoded in {θ 12 , θ 13 , θ 23 } (see appendix B). In this fit, θ is already set to the value obtained from the fit to the CKM matrix in eq. (43). Although different PMNS analyses [22,23] show some sensitivity to the phase δ , we do not include that information in the fit since we are precisely interested in its prediction. The fit gives the following two solutions,   ) show the values that correspond to δ PMNS attending to the ∆χ 2 profiles for δ obtained for normal and inverted neutrino mass orderings in [22]. We stress that using the information on CP violation in the quark sector, we have been able to predict the phase in PMNS using the connection that SCPV provides in this model; in particular, Solution 1 has δ PMNS = 1.629π, which is in good agreement with the most likely values in PMNS analyses. 6. With the values of e τ SFCNC in eq. (56), we have predictions (see appendix A) such as (59) Using eq. (45), we have the sharp range which should be seen or disproved in the near future since the current bound is 6 Br(h → eτ +ēτ ) Exp ≤ 4.7 × 10 −3 [31,32,33] (although there is some freedom in Γ(h)/Γ(h SM ), it does not modify substantially this conclusion).

Conclusions
We have discussed the possibility of having a framework where there is a connection between the CP violations in the quark and the lepton sectors, parametrised by the two phases δ CKM and δ PMNS . In general, it is not possible to establish this connection, since the Yukawa couplings Y f -with f = u, d, , ν -generating the quark and lepton masses are complex matrices, with no relation between Y u,d and Y ,ν . In this paper we have investigated this connection in a framework where CP violation in the quark and lepton sectors have a common origin, being both generated by a complex vacuum phase. We have pointed out that in order to construct experimentally viable models of this class, some conditions have to be satisfied, namely the vacuum phase has to be able to generate complex CKM and PMNS matrices. This is not an easy task since it is assumed that CP is spontaneously broken, so the Yukawa couplings are real. We have shown that in order to generate a complex CKM matrix, one has to have SFCNC both in the up and down quark sectors. An entirely analogous requirement applies to the lepton sector, where the generation of a complex PMNS matrix also requires the presence of scalar leptonic FCNC. Since there are stringent bounds on these SFCNC, we consider generalised BGL 2HDMs where there is natural suppression of some of the most dangerous SFCNC. We have shown that within some of the generalised BGL models, there is indeed a connection between δ CKM and δ PMNS . The interplay among CPV and the existence of SFCNC makes that these relations are quite involved implying connections or predictions for processes mediated by SFCNC in all the sectors: up, down quarks and charged leptons 7 . To clarify all these relations we have worked a subclass of our most general model of SCPV, but guided by two important hints: (i) an experimental fact: the absence of any convincing evidence of the presence of SFCNC at the actual level of precision; (ii) a theoretical discovery: the necessity of having at least one type of SFCNC in each sector: quarks up and down, charged leptons and even in the neutrino sector. Therefore we have worked with models that have the minimal amount of SFCNC needed to keep SCPV. These simplified models verify the MFV ansatz. Because they are controlled by the four unit vectorsr [u] ,r [d] ,r [ν] ,r [e] having a zero in some entry, there are 3 4 = 81 possible models of this type. In the model discussed in section 6, the connection between δ CKM and δ PMNS gives a prediction for δ PMNS in agreement with recent PMNS analyses. In the appendices we explain how the remaining models confront experimental data requirements.

A Higgs SFCNC
The Yukawa couplings of the Higgs-like scalar h, following section 2, read The sign ± appearing with R 3s is + for f = d, and − for f = u, ν. One can immediately read the SFCNC terms relevant for h → eτ (neglecting m e terms) In the SM, the coupling of h SM toτ τ is simply − mτ v h SMτ τ and, with |R 21 − iR 31 | 2 = 1 − R 2 11 , one can obtain straightforwardly eq. (59). Similarly, for t → hc decays (neglecting m c terms) the interaction term is
To illustrate how one fits the CKM and PMNS matrices, in the models under consideration, to the experimental information condensed in eqs. (67) and (68), let us consider the CKM case. In terms of the parameters of the particular model under consideration, a CKM matrix V is computed. Then, one computes (for simplicity, we use here 1, 2, 3 indices rather than u, c, t, d, s, b) and δ = arg These are the values, obtained in a rephasing invariant manner, to be compared with eq. (67) (in the actual fit, this comparison gives the usual likelihood function of the model parameters, which is maximized). For the PMNS matrix, we follow the same procedure with the values in eq. (68), and include no constraint on the phase δ , as discussed in section 6.

C MFV models
In section 6 we have introduced a class of models in which SFCNC are only present in a single transition in the up and down quark sectors: that is the minimal assumption which can produce a complex CKM matrix. This assumption is then extended to the lepton sector. One can consider, in principle, 3 4 = 81 different models in this class (3 choices for the vanishing component ofr [f ] for f = u, d, ν, e). In the lepton sector, however, following the considerations on µ e, we directly discard models withr [e]3 = 0. In the following we analyse how there is only one model which satisfies some basic requirements, the model considered in section 6. We use the following notation to identify the different models: model (ud) = (jk), j, k = 1, 2, 3, is the model in whichr [u]j = 0 andr [d]k = 0: j identifies which generation does not have SFCNC in the up quark sector and similarly for k in the down quark sector. For the quark sector, the model in section 6 withr [u] and r [d] in eq. (40) is model (ud) = (12). The extension of the notation to the lepton sector is straightforward: withr [ν] ,r [e] in eq. (49), the complete label of the model is (ud, ν ) = (12,32). In principle, one would need to perform an analysis of each model similar to the one in [10], extended to the lepton sector, a task which is beyond the scope of this work. It is nevertheless possible to consider simpler analyses and a few requirements to understand how all but one model have to be discarded. The first requirement in the quark sector is that a CKM matrix in agreement with data can be obtained. All 9 models (ud) = (jk), j, k = 1, 2, 3, can give a good CKM matrix. The next requirement concerns SFCNC: if one sorts the absolute value of the components ofr [f ] , |r  [35]), which cannot be satisfied within models (ud) = (31) and (32). The only model left in the quark sector is (ud) = (12). The problem of exploring the viable models is then reduced to analyse the requirements on the lepton sector of the 6 models (ud, ν ) = (12, jk), j = 1, 2, 3, k = 1, 2. Models (ν ) = (11), (21), (31), can produce a realistic PMNS matrix, but then |r [e]2r[e]3 | is close to its maximal value, 1/2, yielding too large Br(h → µτ ) predictions. On the other hand, models (ν ) = (12), (22), cannot produce a realistic PMNS matrix; model (ν ) = (32), on the contrary, can produce a good PMNS matrix. We are left with (ud, ν ) = (12, 32), the model considered in section 6, as the only viable candidate.

D Parametrizations
In section 6.2, eq. (51), a permutation P 23 is introduced in the parametrization of O ν L which deviates, apparently, from the general parametrization introduced in eq. (32). While R 13 (p) and R 23 (p) naturally have one vanishing component in the third row (the one controlling SFCNC) in position 2 and 1 respectively, parametrizing the MFV models introduced in section 6 in which the vanishing component is in position 3 requires an additional consideration. This permutation P 23 is introduced in order to properly implement that case. One can indeed rewrite that is, P 23 can be simply viewed as a product of a fixed rephasing and a 2-3 rotation with fixed angle π/2 and thus there is nothing essentially different with respect to the other cases. It is to be mentioned that while the PDG parametrization of CKM and PMNS in appendix B is such that one can reduce the ranges of θ ij to [0; π/2] (while δ ∈ [0; 2π[) through rephasings, in our case the ranges of the parameters p f 1 , p f 2 in the different orthogonal matrices O f L entering CKM and PMNS, require some care since they cannot be completely reduced to those ranges (the form of V and U in eqs. (24), (34)-(37), is different from eq. (66)).