Probing New Physics Signals with Symmetry-Restored Yukawa Textures

We implement a symmetry violation guideline into a two-Higgs-doublet model embedded with three right-handed neutrinos, and exploit the generic Yukawa structures of the model via a hypothetical symmetry restoration of a global $U_Q(1)^3\times U_L(1)^3$ symmetry. We then apply a mass-powered parametrization to construct the phenomenologically motivated Yukawa interactions, which enables us to incorporate correlatively the neutrino mass, dark matter, as well as the lepton-flavor universality violations in $R_{K^{(*)}}$ and $R_{D^{(*)}}$. Specifically, two atmospheric-scale neutrino masses are generated by a low-scale seesaw mechanism, while the much lighter one, being of $\mathcal{O}(10^{-6})~{\rm eV}$, is fixed by a $7.1~{\rm keV}$ sterile neutrino dark matter produced primordially by the freeze-in mechanism. On the other hand, the neutrino and the charged-lepton mass hierarchies encoded in the mass-powered textures can naturally account for the $R_{K^{(*)}}$ and the $R_{D^{(*)}}$ anomalies, respectively. As a further application, a milder discrepancy of the muon $g-2$ is observed, which has also been implied by the recently refined lattice results.

See Ref. [16] for a collection of refined SM predictions. Benefiting from cancellations of large parts of the hadronic uncertainties, the R K ( * ) and R D ( * ) anomalies, if confirmed, would hint at lepton-flavor universality violation (LFUV) that is not accountable within the SM.
In recent years, a two-Higgs-doublet model (2HDM) with generic Yukawa structures [17] has been investigated intensively. This framework can explain the R D ( * ) anomalies with sizable top-charm couplings [18][19][20] and, when embedded with three right-handed neutrinos, address the neutrino mass and the R K ( * ) anomalies simultaneously [20], or the 7.1 keV sterile neutrino dark matter [21,22]. One might, therefore, expect the general 2HDM with three right-handed neutrinos as a promising candidate to unify these NP signals. However, the model invokes generically overabundant (unknown) Yukawa parameters and hence limits its capabilities of theory prediction and correlation. Thus, understanding the origin of the phenomenologically motivated Yukawa structures becomes especially crucial to further exploit the model.
In this paper, we shall be concerned with the possible correlations among these generic Yukawa couplings in a 2HDM with three right-handed neutrinos (2HDM+3N R ). To this end, we consider here a hypothetical symmetry restoration of a global U Q (1) 3 × U L (1) 3 symmetry that should have been broken by the Yukawa couplings in the mass-eigenstate basis. By attributing the symmetry breaking sources to the known fermion masses (or equivalently the dimensionless Yukawa eigenvalues), we then find that the unknown Yukawa matrices of the model can be built out of these known fermion masses. To visualize this idea phenomenologically, we propose a simple masspowered parametrization of these non-diagonal Yukawa matrices, which enables us to incorporate correlatively the neutrino mass, DM, as well as the R K ( * ) and R D ( * ) anomalies within a single framework.
The model.-In many model buildings concerned with the Yukawa structures, the usual strategy is to imple-ment symmetry invariance into the Yukawa sector by invoking some heavy dynamical fields, as adopted e.g. in the Froggatt-Nielsen mechanism [23] and the minimal flavor violation hypothesis [24,25]. Such considerations are strongly supported by the low-energy flavor observations. Nevertheless, in the low-energy regime, there seems to be no exact but rather broken symmetries, or approximate symmetries with some random perturbations. Assuming that the broken symmetries are not recovered by any dynamical fields, we exploit here an interesting scenario that does not receive much attention. Explicitly, we shall assume that the Yukawa interactions have indeed completely and explicitly broken some family symmetries in the mass-eigenstate basis, but with limited and correlative perturbations, or symmetry breaking sources (SBS). Furthermore, the SBS are located only in the Yukawa sector such that, when the SBS vanish, the corresponding symmetries can be restored universally and completely throughout the Yukawa Lagrangian.
It should be emphasized that, if there are various unrelated SBS, we must tune these SBS synchronously to vanish, in order to make the symmetry restoration universal throughout the whole Lagrangian. This is technically feasible but not quite natural. In this context, we shall assume further that the various SBS have some common origins, so that the naturally synchronous vanishing of the SBS is triggered by these common origins. This prescription allows the Yukawa structures to be generated only by the common SBS, and hence reduces dramatically the number of free parameters. Obviously, such a setup for constructing the Yukawa textures is not based on the symmetry invariance principle, but rather on the symmetry violation guideline, which is concerned with how the flavor symmetries are broken completely by the Yukawa interactions in the mass-eigenstate basis, and can be introduced via the following criteria: (i). The Yukawa interactions in the mass-eigenstate basis have broken explicitly some family symmetries but only via limited and correlative SBS.
(ii). The Yukawa textures generated by these SBS are subject to a hypothetical symmetry restoration. This indicates that the corresponding symmetry will be recovered universally and completely throughout the whole Lagrangian under vanishing SBS.
(iii). There are no other adjustable free parameters beyond those SBS, so that the symmetry restoration is triggered only by vanishing SBS.
The universal symmetry restoration in criteria (ii) corresponds to the natural setup in which all the symmetry perturbations can vanish synchronously when the common SBS origin is tuned to zero, while the criteria (iii) forces the number of free parameters to be further reduced during the construction of the Yukawa matrices.
In the rest of our work, we shall apply these criteria in the 2HDM+3N R framework to construct the Yukawa interactions, which are found to be able to account for the NP signals considered. In the mass-eigenstate basis, the Yukawa interaction Lagrangian can be parametrized in a compact form as Here, Q L ≡ (V † u L , d L ) T and E L ≡ (U * ν ν L , e L ) T , with V and U ν corresponding respectively to the quark and neutrino mixing matrices observed in experiments, no longer form the SU L (2) doublets. More importantly, we have particularly neglected, in Eq. (8), the terms that are either proportional to the active neutrino mass or further suppressed by the light-heavy neutrino mixing parameter [26]. Note that the ma-tricesŶ u,d, = √ 2m u,d, /v are already in a diagonal form, butỸ ν is not. In addition, the matrices Y remain nondiagonal and encode the unknown Yukawa interactions arising in the 2HDM+3N R framework.
The two scalar doublets H 1,2 are given in the Higgs basis as [27] with the vacuum expectation value v 246 GeV. The neutral scalars φ 1,2 are the superposition of the two mass eigenstates H and h via φ 1(2) = cos α(− sin α)H + sin α(cos α)h, with the mixing angle determined by where λ 6 is the quartic mixing coupling of the term in the Higgs-basis scalar potential [27]. Let us consider both Z 2 and CP symmetries that are conserved by the Higgs-basis scalar potential but violated by the Yukawa sector. In this case, λ 6 = 0 and hence the two neutral Higgs bosons (H and h) decouple from each other. Furthermore, choosing one of the solutions, α = π/2, h would return to the SM Higgs boson (the so-called alignment limit). As a consequence, there exist no tree-level flavor-changing neutral currents (FCNC) involving h, as is observed in the SM.
It can be readily seen that, in the mass-eigenstate basis, Eq. (8) breaks a global U Q (1) 3 × U L (1) 3 symmetry, with the following transformation rules: The specific charge assignments are irrelevant, provided that no particular relations among these charges would change the broken symmetry. This makes it possible for us to consider the most general case where the transformation rules specified by Eq. (11) exhibit a complete breaking of In the following, we shall seek for the common SBS that satisfy the criteria given above. Let us take the neutral currents for a start. Considering the quark neutral currents (as well as the quark mass terms) of L H1 , one can see that the U Q (1) i symmetry would be restored when the i-th quark masses hypothetically vanish, m u,d i = 0, since these terms are all determined by the Yukawa ma-tricesŶ u,d . Based on this observation, we can conjecture the quark masses, or equivalently the dimensionless Yukawa eigenvalues y f i defined by as the common SBS. Turning now to the quark neutral currents of L H2 , the non-diagonal Yukawa matrices Y u,d are also expected to be generated by the common SBS. Following the guideline specified by the three criteria given above, we can construct Y u,d in the following form: where f = u, d, and the factorsñ f i denote the powers to which the generic Yukawa matrices Y f are built out of the eigenvalues y f i . The reason for such a power realization is that these power factors, while being free, do not play the role of SBS. In fact, these powers can be readily understood by noting that the Yukawa matrices are built out of the eigenvalue products themselves, as required by the criterion (iii). It is then readily seen that, when the i-th quark masses hypothetically vanish, or y u,d i = 0, the corresponding U Q (1) i symmetry would be restored throughout the neutral currents of Eq. (8).
Next, let us consider the quark charged currents. In this part, due to the presence of quark mixing, the SBS role is now played by both the mixing matrix and the quark masses. To still retain the masses as the common SBS, some correlation should be established between the quark masses and the mixing matrix, such that the mixing would synchronously vanish when quarks are hypothetically massless. Since the general mass matrix M f is constructed in terms of the mass eigenvalues through a bi-unitary rotation, M f = V f † LM f V f R , and the combination of up-and down-quark rotations gives the physical CKM mixing matrix, V = V u L V d † L , the same guideline allows us to build the general mass matrix with the following texture: Again, the factor n ij denotes the power to which the general mass matrix M f is expanded in terms of the Yukawa eigenvalues. Note that the general mass matrix cannot exhibit a form like Eq. (13), otherwise the mass matrix would be of rank one and unacceptably render only one fermion generation massive, while this is not the case for Y f as it is not responsible for fermion masses in the considered mass-eigenstate basis. Eq. (14) indicates that, if the general mass matrices in a non-physical basis were also constructed in terms of the Yukawa eigenvalues, the flavor spectra would be induced by the random powers (non-SBS parameters) n ij without affecting the broken symmetry pattern after changing to the physical basis. In addition, the synchronous vanishing of the physical mixing matrix triggered by fermion masses can now be seen as follows: In this way, the hypothetically vanishing of quark masses as the common SBS is able to prompt a universal U Q (1) 3 symmetry restoration throughout the Yukawa Lagrangian specified by Eq. (8). In fact, the same restoration is also valid in the W ± -mediated charged currents due to Eq. (15).
The above analysis can be directly applied to the charged-lepton sector. However, it is non-trivial to see a hypothetical symmetry restoration in the neutrino sector. This can be seen in several ways. Firstly, unlike the Dirac fermions, the left-handed neutrinos do not combine with the corresponding right-handed counterparts to form the physical masses. Secondly, as the generation of active neutrino masses relies on the seesaw mechanism which entails a non-singular right-handed neutrino mass matrix M R , turning hypothetically the right-handed Majorana neutrinos massless would cause conceptual issues for the active neutrino masses. Finally, as the seesaw mass formula, is a leading-order result, and the light-heavy neutrino mixing is also given at the leading order, taking into account the higher-order terms would sophisticate the exploitation of common SBS in the neutrino sector, although these higher-order corrections do not significantly affect the observables considered. Given these observa-tions and in order to simplify the analysis, we shall assume that the heavy neutrinos do not carry the hypothetical U (1) charges in the physical mass-eigenstate basis, and set N = 0 in the transformation rules defined by Eq. (11). In this case, the same criteria lead us to build the neutrino Yukawa matrices with only the rows specified, (Ỹ ν ) ij , (Y ν ) ij ∝ y ν i , where y ν i denote the dimensionless Yukawa eigenvalues of the effective neutrino mass matrix. The j-columns of (Ỹ ν ) ij and (Y ν ) ij need not be specified, because the right-handed counterparts do not participate in the symmetry violation. Nevertheless, the effective active neutrino mass matrix in the non-physical basis can still exhibit a texture like Eq. (14) due to the seesaw formula specified by Eq. (16), which is necessary to guarantee the synchronous vanishing of neutrino mixing via Eq. (15). In this way, the hypothetically vanishing of charged lepton and active neutrino masses (or equivalently their respective Yukawa eigenvalues) renders a universal U L (1) 3 symmetry restoration throughout the Yukawa Lagrangian (Eq. (8)), as well as in the W ± -mediated charged currents.
In the current work, as we are mainly concerned with the effects of additional Yukawa interactions on explaining the NP signals considered, we shall focus only on Y f but without delving into the flavor structures of Eq. (14). Eq. (13), on the other hand, can stir up even richer phenomenologies, as such a structure indicates that, for the additional Yukawa interactions, the flavor-specific couplings y f i can be either enhanced or suppressed by the non-universal powersñ i . Here we are interested in how large the tree-level FCNC couplings of the scalars to the lighter flavors can be allowed when those couplings to the heavier flavors are responsible for the NP signals considered. To this end, we propose a simple mass-powered parametrization based on Eq. (13) by scaling the SBS with a common power and a dimensional normalization factor, where Λ denotes the dimension normalization factor, n is the power, and θ represents the possible CP -violating phase. Note that, as Y ν is not specified in the j-columns, we have phenomenologically introduced the right-handed neutrino masses M i , with the index collocation being different from the Dirac fermions. In addition, the dimension normalization factors in the neutrino sector Λ j can be traced back to two origins of heavy neutrinos: TeVscale heavy neutrinos and keV-scale sterile neutrino DM, as will be shown in the subsequent discussions. It should also be mentioned that the index collocation in Y ν is motivated by the explanation of R K ( * ) anomalies under the i → j γ constraints, as will be clarified later. Finally, Eq. (17) indicates that the FCNC mediated by the neu-tral (pseudo)scalars (H, A) would be controlled by the fermion mass hierarchies. This can provide, therefore, a compelling suppression of FCNC in the lighter-flavor sector while allowing for sizable FCNC to explain the B-meson anomalies.
Phenomenology.-We now apply the mass-powered Yukawa textures to the neutrino mass, DM, as well as the R K ( * ) and R D ( * ) anomalies. In addition, we calculate the NP contribution to the muon g − 2 as a further application. For numerical analyses, we set Λ u = m t , Λ = m τ and Λ d m b to allow sizable top-and tau-associated Yukawa couplings, which are found to be essential for the LFUV explanations and, meanwhile, to suppress the FCNC in the down-quark sector. On the other hand, possible CP -violating phases in Eq. (17) will be neglected unless stated otherwise. The scalar masses are considered around M S O(500) GeV (S = H ± , H, A), to sufficiently alleviate the scalar effects on electroweak precision observables [28].
Neutrino mass and DM : As the neutrino mass spectrum is currently unknown, we cannot apply Eq. (17) directly to the neutrino sector. It is known that, to have a stable neutrino mass generation via a low-scale seesaw mechanism, a lepton-like U (1) symmetry is usually considered in the flavor basis. In Ref. [20], the lepton-like charges in the (primed) flavor basis are chosen as There, the Yukawa couplings involving N 1 and N 3 break the lepton-like U (1) symmetry, and N 2 and N 3 are nearly degenerate. After rotating to the physical masseigenstate basis in which the neutrino mixing information is encoded in the Dirac Yukawa matrixỸ ν , as has been constructed in Eq. (8), one can find that, if the lightest sterile neutrino (denoted as N 1 ) plays the role of keV DM, all the elements of the first column ofỸ ν will be strongly suppressed by the cosmological X-ray observation (see Ref. [29] for an updated review), via the W ± -mediated N 1 → νγ decay [30][31][32] where G F and α em are the Fermi and the fine-structure constants, respectively. Note that, we have expressed the decay width in terms of the Dirac Yukawa matrix, so that the usual constraint on the light-heavy neutrino mixing angle θ ν can be translated onto (Ỹ ν ) i1 by the relation In this case, only two neutrino masses will be generated at the atmospheric scale ∆m atm 0.05 eV [33], and the resulting mass hierarchy assumes either the normal ordering (NO) or the inverted ordering (IO), As the DM candidate, the lightest sterile neutrino N 1 would be a Feebly Interacting Massive Particle [34], since the associated Yukawa couplings that break the leptonlike U (1) symmetry are expected to be suppressed. On the other hand, as the corresponding Dirac Yukawa cou-plingỸ ν is severely constrained, the relic density of N 1 will then be accumulated by the freeze-in mechanism [34] via the decays S → N 1 + l (l = or ν), when the scalars are in thermal equilibrium with the SM bath. The relic abundance is then given by [34] On the other hand, the 7.1 keV N 1 populates the warm DM region [36,37], and the explanation of the 3.5 keV X-ray line via the radiative decay N 1 → νγ can also be consistent with the Lyman-α observation [38]. In Fig. 1, the current DM relic density, Ω N1 h 2 = 0.12 [39], is fitted with a numerical hierarchy Λ ν1 = 10 5 Λ ν2(3) MeV (blue line). Here the choice of the hierarchy between Λ ν1 and Λ ν2(3) , as well as the approximation Λ ν2 Λ ν3 is sensible, as Λ ν1 corresponds to the DM interactions and the small difference Λ ν2 −Λ ν3 can be treated as the breaking source of the lepton-like symmetry. In addition, we have taken M 2 > M S so as to open the decay channel N 2,3 → S +l, which can sufficiently decrease the lifetime of N 2,3 . The resulting impact of N 2,3 on the effective relativistic degrees of freedom at the Big Bang Nucleosynthesis epoch can be, therefore, neglected safely [40].  [41][42][43][44][45], with the black line corresponding to the best-fit point. The yellow and green regions are excluded by the µ → eγ and τ → µγ constraints, respectively. The light blue region is excluded by Z → µ + µ − and the red one by R W µe . Finally, the gray region is excluded by the perturbative unitarity requirement Y > √ 4π. The blue line corresponds to the DM relic density ΩDM h 2 = 0.12 [39] with a numerical hierarchy Λν1 = 10 5 Λ ν2(3) MeV.
effective Wilson coefficients: corresponding to the semi-leptonic operators Here, It is seen that, due to the neutrino mass hierarchies encoded in the mass-powered parametrization (see Eq. (17)), the necessary flavor non-universal couplings appear naturally, and the LFUV observables R K ( * ) can be explained by the fact that M 2(3) M 1 and m 2 > m 1 . Nevertheless, such a muon-neutrino specific Yukawa coupling may give significant effects on muon-associated observables. To clarify that a successful explanation of the R K ( * ) anomalies is possible with such a coupling, we further take into account the primary constraints from µ → eγ, τ → µγ, Z → µ + µ − , as well as the ratio R W µe ≡ Γ(W → µν)/Γ(W → eν) that tests the leptonflavor universality [39]. Finally, the perturbative unitarity requirement Y > √ 4π is also imposed. From Fig. 1, it is clearly seen that the updated model-independent analyses [41][42][43][44] that prefer C NP 9µ = −C NP 10µ < 0 with a significance at the level of ∼ 4σ can be well reproduced here, as shown by the shaded band.
Here we have taken the IO pattern of the active neutrino masses (see Eq. (23)). Note that, for the sake of plotting, we have compiled the possibly negative values of Y ν to that of Λ ν2 , since the considered observables depend on |Y ν | 2 . With such a numerical setup, the resolution of the R K ( * ) anomalies exists under all the constraints considered. However, we find that if the active neutrino masses have the NO pattern, the parameter region allowed by the R K ( * ) explanation will be completely ruled out by the τ → µγ constraint.
Furthermore, we have also numerically confirmed that, if Y ν exhibits a texture similar to that of the Dirac fermions, rather than the one proposed in Eq. (17), the resolution of the R K ( * ) anomalies would be excluded either by µ → eγ or by τ → µγ, no matter which hierarchies of the active neutrino masses are taken. Therefore, being in association with the lightest heavy neutrino as a 7.1 keV DM, the mass-powered texture of Y ν in Eq. (17) serves a twofold role, on the one hand accounting for the LFUV in R K ( * ) and on the other hand predicting an IO pattern of the active neutrino masses. R D ( * ) : The NP effect on b → c i ν transitions arises from the tree-level charged-Higgs contribution. Neglecting the suppressed down-quark Yukawa couplings, the effective Hamiltonian is given by It can be readily seen that, the lepton-flavor nonuniversality is induced by mass hierarchy encoded in Y . For numerical analysis, we implement the Wilson coefficient to the updated formulae in Ref. [46]. Besides, to verify that our model can provide a feasible resolution of the R D ( * ) anomalies, we consider the bounds from the branching ratio B(B c → τ ν) derived from the B c lifetime [47][48][49], and the D * longitudinal polarization fraction [50], which differs from its SM prediction [51], by ∼ 1.6σ. Following Ref. [52], we also consider the mass difference ∆M s in the B s −B s system and the inclusive B → X s γ branching ratio, both of which receive sizable NP effects from the top-associated Yukawa couplings. As shown in Fig. 2, a 1σ-level explanation of R D ( * ) anomalies is feasible under the constraint B(B c → τ ν) < 60% [46,53] [55,56]. On the other hand, by taking n 1, the chargedlepton mass hierarchies encoded in Y would induce significant effects on the τ channel but suppress those in the µ/e modes. This ensures negligible effects on B(B → D ( * ) µν)/B(B → D ( * ) eν) [57] and, at the same time, provides a natural explanation of the R D ( * ) anomalies.
Concerning the constraint from H ± → τ ν searches at the LHC, it is found that the decay width of H ± is now dominated by B(H ± → τ ν) (30%) and B(H ± → tb) (70%), and the strength of cross section times branching ratio, σ(pp → H ± ) × B(H ± → τ ν), will be reduced compared to the case with B(H ± → τ ν) 90%. Therefore, the compatibility between the R D ( * ) explanation and the LHC constraint from H ± → τ ν searches [58] can be realized, as the latter becomes weaker by at least a factor of three. ∆a µ : As a further application of the mass-powered parametrization, let us now consider the longstanding puzzle observed in the muon g − 2, a µ = (g µ − 2)/2. The current value [59], exhibits a 3.7σ discrepancy between theory and experiment. However, the present SM prediction is still plagued by large hadronic uncertainties, leading to a possible range of 0.7−4.2σ deviations [60]. Furthermore, a milder discrepancy, including the no-NP solution, has also been implied by the recently refined lattice calculations of the hadronic contributions to the muon g − 2 [61][62][63].  [59] are shown in green. The region below the black, blue, and red curves are excluded by τ → 3µ, µ → eγ and τ → µγ, respectively. θ u 33 = π is chosen to obtain a maximal NP contribution to ∆aµ.
Here, the NP contributions arise from both the oneloop diagrams involving the muon-tau Yukawa couplings, as well as the two-loop Barr-Zee diagrams involving the top and tau Yukawa couplings even in the degenerate mass regime of the scalar bosons [52]. To sum over these contributions, we have taken the formulae [52] by rescaling the Yukawa couplings with the mass-powered parametrization given in Eq. (17). In this regime, we show, in Fig. 3, the parameter regions required to explain the current g − 2 data, as well as the constraints from lepton-flavor changing (LFV) processes, with the primary one resulting from τ → 3µ. It is seen that the mass-powered parametrization cannot explain the current 3.7σ discrepancy of the muon g − 2 at 1σ level. Instead, the 3σ-level accommodation suggests a milder discrepancy or even no-NP solution of the muon g − 2, which has also been hinted by the latest refined lattice results [61][62][63].
Conclusions.-We have considered an interesting scenario where the Yukawa interactions are constructed by a symmetry violation guideline, rather than by the usually adopted symmetry invariance principle. Implementing such a guideline into the 2HDM+3N R framework, we found that the model becomes very constrained but its predictive power is enhanced dramatically. By attributing the SBS of a hypothetical U Q (1) 3 × U L (1) 3 symmetry to the known fermion masses (or equivalently the dimensionless Yukawa eigenvalues), we found that the unknown Yukawa matrices of the model can be built out of the known fermion masses. As a phenomenological application, we have investigated the interplay between the severely constrained light-flavor FCNCs and the explanations of NP signals from heavy-flavor FCNCs by proposing a simple mass-powered parametrization of the additional Yukawa matrices.
The mass-powered texture renders a correlative incorporation of the neutrino mass, DM, as well as the LFUV observed in R K ( * ) and R D ( * ) . Let us summarize the main results. Two atmospheric-scale neutrinos are generated by a U (1)-protected low-scale seesaw mechanism, while the remaining much lighter one is fixed by a 7.1 keV sterile neutrino DM. Such a DM can explain the 3.5 keV X-ray line, and its relic density is produced by the freezein mechanism via the thermalized scalar decays. In light of the mass-powered parametrization, the IO pattern of the active neutrino masses, 0 m 3 m 1 < m 2 , together with a strong hierarchical heavy neutrino masses, M 1 M 2 M 3 , accounts for the R K ( * ) anomalies, while the charged-lepton mass hierarchies explain the R D ( * ) data. Finally, a milder discrepancy of the muon g − 2 suggested here can be served as a future probe for the particular parametrization.