Probing New Physics Signals with Symmetry-Restored Yukawa Textures

TL;DR

This paper develops a symmetry-restored Yukawa texture framework within a two-Higgs-doublet model to connect neutrino masses, dark matter, and lepton-flavor universality violations, providing a unified explanation for several anomalies.

Contribution

It introduces a novel symmetry-restoration approach to Yukawa textures that naturally explains neutrino masses, dark matter, and flavor anomalies in a unified model.

Findings

Neutrino masses generated by low-scale seesaw mechanism.

Sterile neutrino as 7.1 keV dark matter candidate.

Model accounts for R_K(*) and R_D(*) anomalies.

Abstract

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.

Figures

Equations

$$R_{K^{(\ast )}}=\frac{\mathcal{B}(B\to K^{(\ast )}{\mu }^{+}{\mu }^{-})}{\mathcal{B}(B\to K^{(\ast )}{e}^{+}{e}^{-})}$$

$$R_{D^{(\ast )}}=\frac{\mathcal{B}(B\to D^{(\ast )}\tau \bar{\nu })}{\mathcal{B}(B\to D^{(\ast )}e(\mu )\bar{\nu })},$$

$$R_K=0.846_{-0.054}^{+0.060}(\text{stat}{)}_{-0.014}^{+0.016}(\text{syst}),$$

$$R_{K^{\ast }}=0.69_{-0.07}^{+0.11}(\mathrm{stat})\pm 0.05(\mathrm{syst}),$$

$$R_D=0.340\pm 0.027(\text{stat})\pm 0.013(\text{syst}),$$

$$R_{D^{\ast }}=0.295\pm 0.011(\text{stat})\pm 0.008(\text{syst}),$$

$$R_D^{\mathrm{SM}}=0.299\pm 0.003,R_{D^{\ast }}^{\mathrm{SM}}=0.258\pm \mathrm{0.005.}$$

$${\mathcal{L}}_Y$$

$${\mathcal{L}}_{H_1}$$

$$-{\bar{E}}_L{\tilde{Y}}_{\nu }{\tilde{H}}_1{N}_R+H.c.,$$

$${\mathcal{L}}_{H_2}$$

$$-{\bar{E}}_L{\mathcal{Y}}_{\nu }{\tilde{H}}_2{N}_R+H.c..$$

$$\displaystyle H_{1}=\left(\begin{array}[]{c}G^{+}\\ \frac{v+\phi_{1}+iG}{\sqrt{2}}\end{array}\right),\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ H_{2}=\left(\begin{array}[]{c}H^{+}\\ \frac{\phi_{2}+iA}{\sqrt{2}}\end{array}\right),$$

$$\tan (2\alpha )=\frac{2{\lambda }_6{v}^2}{\sqrt{(M_H^2-M_h^2{)}^2-4{\lambda }_6^2{v}^4}},$$

$$(u_L,d_L{)}_i$$

$$({\nu }_L,{\ell }_L{)}_i$$

$$y_i^f{y}_j^f{\delta }_{ij}v=m_i^f,$$

$${\mathcal{Y}}_{ij}^f=(y_i^f{)}^{{\tilde{n}}_i^f}\times (y_j^f{)}^{{\tilde{n}}_j^f},$$

$$M_f=v(y_i^f\times y_j^f{)}^{n_{ij}}.$$

$$m_{i_0}^f(y_{i_0}^f)=0\Rightarrow \{\begin{array}{l}(V_L^f{)}_{i_{0}j}=(V_L^f{)}_{j{i}_0}=0,j\ne i_0,\\ (V_L^f{)}_{i_0{i}_0}=1,\end{array}.$$

$$M_{\nu }\simeq -M_D{M}_R^{-1}{M}_D^T,$$

$$({\mathcal{Y}}_f{)}_{ij}=\{\begin{array}{l}{\left(\frac{m_i\times m_j}{{\Lambda }^2}\right)}^n{e}^{-i{\theta }_{ij}},f=u,d,\ell ,\\ {\left(\frac{m_i\times M_i}{{\Lambda }_j^2}\right)}^n{e}^{-i{\theta }_{ij}},f=\nu ,\end{array}$$

$$L_{N_{R1}^{\prime }}=0,L_{N_{R2}^{\prime }}=-L_{N_{R3}^{\prime }}=1,L_{E_{Li}^{\prime }}=L_{e_{Ri}^{\prime }}=1.$$

$$\Gamma (N_1\to \nu \gamma )=\frac{9\sqrt{2}{\alpha }_{\mathrm{em}}{G}_F}{1024{\pi }^4}\sum _{i=1}^3|({\tilde{Y}}_{\nu }{)}_{i1}{|}^2{M}_1^3,$$

$${\theta }_{\nu }^2\simeq \frac{\sqrt{2}}{4G_F{M}_1^2}\sum _{i=1}^3|({\tilde{Y}}_{\nu }{)}_{i1}{|}^2.$$

$$\text{NO}:0\simeq m_1\ll m_2<m_3,\text{IO}:0\simeq m_3\ll m_1<m_2.$$

$${\Omega }_{N_1}{h}^2$$

$$\times \left(M_A^{-1}+M_H^{-1}+2M_{H^{\pm }}^{-1}\right),$$

$$:m_1\simeq 10^{-6}\text{leavevmode}\mathrm{eV},\text{leavevmode}{m}_2=0.009\text{leavevmode}\mathrm{eV},m_3=0.05\text{leavevmode}\mathrm{eV},$$

$$:m_3\simeq 10^{-6}\text{leavevmode}\mathrm{eV},\text{leavevmode}{m}_1=0.049\text{leavevmode}\mathrm{eV},\text{leavevmode}{m}_2=0.05\text{leavevmode}\mathrm{eV}.$$

$$C_{9\ell }^{\mathrm{NP}}$$

$$=\frac{G_F^{-1}{\alpha }_{\mathrm{em}}^{-1}}{64\sqrt{2}\pi M_W^2}\sum _{i=1}^3|({\mathcal{Y}}_{\nu }{)}_{\ell i}{|}^2|({\mathcal{Y}}_u{)}_{33}{|}^2\mathcal{G}({\lambda }_{\pm },{\lambda }_i,{\lambda }_t),$$

$${\mathcal{O}}_{9(10)}=\frac{{\alpha }_{\mathrm{em}}}{4\pi }\left(\bar{s}{\gamma }_{\mu }{P}_{L}b\right)\left(\bar{\ell }{\gamma }^{\mu }({\gamma }_5)\ell \right).$$

$$\mathcal{G}(x,y,z)$$

$$-\frac{x}{(x-y)(x-z)}.$$

Full text

Introduction.—The neutrino mass and dark matter (DM) are two well-known signatures that require new physics (NP) beyond the Standard Model (SM). Towards solving these two momentous issues in particle physics, some tantalizing NP signals observed in B-meson decays, such as the ratios of the branching fractions

[TABLE]

and

[TABLE]

can provide complementary incentives in pinning down the underlying NP theory Bifani et al. (2019); Ciezarek et al. (2017).

With the updated LHCb measurement in the dilepton invariant mass squared range $q^{2}=[1.1,6.0]\leavevmode\nobreak\ {\rm GeV}^{2}$ Aaij et al. (2019),

[TABLE]

together with the previous $R_{K^{*}}$ data Aaij et al. (2017),

[TABLE]

discrepancies at the level of $\sim 2.5\sigma$ are found with respect to the corresponding SM predictions (both being equal to one up to a few percent corrections Bordone et al. (2016)). The new Belle result of $R_{K^{*}}$ Abdesselam et al. (2019a), given its sizable uncertainty, is also compatible with the LHCb measurement Aaij et al. (2017). For the ratios $R_{D^{(*)}}$, the latest combination of the BaBar Lees et al. (2012, 2013), Belle Huschle et al. (2015); Hirose et al. (2017, 2018); Abdesselam et al. (2019b), and LHCb Aaij et al. (2015, 2018a, 2018b) measurements performed by the Heavy Flavor Averaging Group gives Amhis et al. (2019a)

[TABLE]

which, after taking into account their correlation of $-0.38$, exhibit a $3.1\sigma$ deviation from the SM predictions,

[TABLE]

See Ref. Amhis et al. (2019a) 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 Crivellin et al. (2013) has been investigated intensively. This framework can explain the $R_{D^{(\ast)}}$ anomalies with sizable top-charm couplings Crivellin et al. (2016); Iguro and Tobe (2017); Li et al. (2018) and, when embedded with three right-handed neutrinos, address the neutrino mass and the $R_{K^{(\ast)}}$ anomalies simultaneously Li et al. (2018), or the $7.1\leavevmode\nobreak\ {\rm keV}$ sterile neutrino dark matter Adulpravitchai and Schmidt (2015a); Baumholzer et al. (2018). 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}\times 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 mass-powered parametrization of these non-diagonal Yukawa matrices, which enables us to incorporate correlatively the neutrino mass, DM, as well as the $R_{K^{(\ast)}}$ and $R_{D^{(\ast)}}$ anomalies within a single framework.

The model.—In many model buildings concerned with the Yukawa structures, the usual strategy is to implement symmetry invariance into the Yukawa sector by invoking some heavy dynamical fields, as adopted e.g. in the Froggatt-Nielsen mechanism Froggatt and Nielsen (1979) and the minimal flavor violation hypothesis Buras et al. (2001); D’Ambrosio et al. (2002). 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

[TABLE]

Here, $Q_{L}\equiv(V^{\dagger}u_{L},d_{L})^{T}$ and $E_{L}\equiv(U_{\nu}^{*}\nu_{L},e_{L})^{T}$, with $V$ and $U_{\nu}$ 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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures), the terms that are either proportional to the active neutrino mass or further suppressed by the light-heavy neutrino mixing parameter $R^{*}\simeq M_{D}M_{R}^{-1}=v\tilde{Y}_{\nu}M_{R}^{-1}/\sqrt{2}$ Ibarra et al. (2010). Note that the matrices $\hat{Y}_{u,d,\ell}=\sqrt{2}m_{u,d,\ell}/v$ are already in a diagonal form, but $\tilde{Y}_{\nu}$ is not. In addition, the matrices $\mathcal{Y}$ remain non-diagonal 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 Branco et al. (2012)

[TABLE]

with the vacuum expectation value $v\simeq 246\leavevmode\nobreak\ {\rm GeV}$. The neutral scalars $\phi_{1,2}$ are the superposition of the two mass eigenstates $H$ and $h$ via $\phi_{1(2)}=\cos\alpha(-\sin\alpha)H+\sin\alpha(\cos\alpha)h$, with the mixing angle determined by

[TABLE]

where $\lambda_{6}$ is the quartic mixing coupling of the term $[(H_{1}^{\dagger}H_{1})(H_{1}^{\dagger}H_{2})+\rm H.c.]$ in the Higgs-basis scalar potential Branco et al. (2012). Let us consider both $\mathbb{Z}_{2}$ and $CP$ symmetries that are conserved by the Higgs-basis scalar potential but violated by the Yukawa sector. In this case, $\lambda_{6}=0$ and hence the two neutral Higgs bosons ($H$ and $h$) decouple from each other. Furthermore, choosing one of the solutions, $\alpha=\pi/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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures) breaks a global $U_{Q}(1)^{3}\times U_{L}(1)^{3}$ symmetry, with the following transformation rules:

[TABLE]

where $\mathcal{Q}_{i}$, $\mathcal{U}_{i}$, $\mathcal{D}_{i}$, $\mathcal{L}_{i}$, $\mathcal{E}_{i}$ and $\mathcal{N}_{i}$ denote the corresponding $U(1)$ charges of the $i$-th fermion generation. 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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures) exhibit a complete breaking of the $U_{Q}(1)^{3}\times U_{L}(1)^{3}$ symmetry.

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 $\mathcal{L}_{H_{1}}$, one can see that the $U_{Q}(1)^{i}$ symmetry would be restored when the $i$-th quark masses hypothetically vanish, $m_{i}^{u,d}=0$, since these terms are all determined by the Yukawa matrices $\hat{Y}_{u,d}$. Based on this observation, we can conjecture the quark masses, or equivalently the dimensionless Yukawa eigenvalues $y^{f}_{i}$ defined by

[TABLE]

as the common SBS. Turning now to the quark neutral currents of $\mathcal{L}_{H_{2}}$, the non-diagonal Yukawa matrices $\mathcal{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 $\mathcal{Y}_{u,d}$ in the following form:

[TABLE]

where $f=u,d$, and the factors $\tilde{n}^{f}_{i}$ denote the powers to which the generic Yukawa matrices $\mathcal{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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures).

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_{L}^{f\dagger}\hat{M}_{f}V_{R}^{f}$, and the combination of up- and down-quark rotations gives the physical CKM mixing matrix, $V=V_{L}^{u}\,V_{L}^{d\dagger}$, the same guideline allows us to build the general mass matrix with the following texture:

[TABLE]

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 $\mathcal{Y}_{f}$ as it is not responsible for fermion masses in the considered mass-eigenstate basis.

Eq. (18) 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:

[TABLE]

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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures). In fact, the same restoration is also valid in the $W^{\pm}$-mediated charged currents due to Eq. (19).

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,

[TABLE]

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 observations 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 $\mathcal{N}=0$ in the transformation rules defined by Eq. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures). In this case, the same criteria lead us to build the neutrino Yukawa matrices with only the rows specified, $(\tilde{Y}_{\nu})_{ij},(\mathcal{Y}_{\nu})_{ij}\propto y^{\nu}_{i}$, where $y^{\nu}_{i}$ denote the dimensionless Yukawa eigenvalues of the effective neutrino mass matrix. The $j$-columns of $(\tilde{Y}_{\nu})_{ij}$ and $(\mathcal{Y}_{\nu})_{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. (18) due to the seesaw formula specified by Eq. (20), which is necessary to guarantee the synchronous vanishing of neutrino mixing via Eq. (19). 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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures)), as well as in the $W^{\pm}$-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 $\mathcal{Y}_{f}$ but without delving into the flavor structures of Eq. (18). Eq. (17), 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 $\tilde{n}_{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. (17) by scaling the SBS with a common power and a dimensional normalization factor,

[TABLE]

where $\Lambda$ denotes the dimension normalization factor, $n$ is the power, and $\theta$ represents the possible $CP$-violating phase. Note that, as $\mathcal{Y}_{\nu}$ 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 $\Lambda_{j}$ can be traced back to two origins of heavy neutrinos: TeV-scale 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 $\mathcal{Y}_{\nu}$ is motivated by the explanation of $R_{K^{(*)}}$ anomalies under the $\ell_{i}\to\ell_{j}\gamma$ constraints, as will be clarified later. Finally, Eq. (21) indicates that the FCNC mediated by the neutral (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 $\Lambda_{u}=m_{t}$, $\Lambda_{\ell}=m_{\tau}$ and $\Lambda_{d}\gg 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. (21) will be neglected unless stated otherwise. The scalar masses are considered around $M_{S}\simeq\mathcal{O}(500)\leavevmode\nobreak\ {\rm GeV}$ ($S=H^{\pm},H,A$), to sufficiently alleviate the scalar effects on electroweak precision observables Haller et al. (2018).

Neutrino mass and DM: As the neutrino mass spectrum is currently unknown, we cannot apply Eq. (21) 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. Li et al. (2018), the lepton-like charges in the (primed) flavor basis are chosen as

[TABLE]

There, the Yukawa couplings involving $N^{\prime}_{1}$ and $N^{\prime}_{3}$ break the lepton-like $U(1)$ symmetry, and $N^{\prime}_{2}$ and $N^{\prime}_{3}$ are nearly degenerate. After rotating to the physical mass-eigenstate basis in which the neutrino mixing information is encoded in the Dirac Yukawa matrix $\tilde{Y}_{\nu}$, as has been constructed in Eq. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures), 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 $\tilde{Y}_{\nu}$ will be strongly suppressed by the cosmological X-ray observation (see Ref. Boyarsky et al. (2019) for an updated review), via the $W^{\pm}$-mediated $N_{1}\to\nu\gamma$ decay Bulbul et al. (2014); Boyarsky et al. (2014, 2015)

[TABLE]

where $G_{F}$ and $\alpha_{\rm 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 $\theta_{\nu}$ can be translated onto $(\tilde{Y}_{\nu})_{i1}$ by the relation

[TABLE]

In this case, only two neutrino masses will be generated at the atmospheric scale $\Delta m_{{\rm{atm}}}\simeq 0.05\leavevmode\nobreak\ {\rm eV}$ Esteban et al. (2017), and the resulting mass hierarchy assumes either the normal ordering (NO) or the inverted ordering (IO),

[TABLE]

As the DM candidate, the lightest sterile neutrino $N_{1}$ would be a Feebly Interacting Massive Particle Hall et al. (2010), since the associated Yukawa couplings that break the lepton-like $U(1)$ symmetry are expected to be suppressed. On the other hand, as the corresponding Dirac Yukawa coupling $\tilde{Y}_{\nu}$ is severely constrained, the relic density of $N_{1}$ will then be accumulated by the freeze-in mechanism Hall et al. (2010) via the decays $S\to N_{1}+l$ ($l=\ell$ or $\nu$), when the scalars are in thermal equilibrium with the SM bath. The relic abundance is then given by Hall et al. (2010)

[TABLE]

where $g^{*}=110.75$ denotes the effective number of relativistic degrees of freedom in the thermal bath around freeze-in temperature $T\sim M_{S}$, after taking into account the additional Higgs bosons besides those in the SM (the heavy Majorana neutrinos are not included and the active neutrinos are of Majorana type).

If $N_{1}$ is the source of $3.5\leavevmode\nobreak\ {\rm keV}$ X-ray line Kusenko (2009), then $M_{1}=7.1\leavevmode\nobreak\ {\rm keV}$ is fixed. In this case, the active neutrino mass spectrum is determined by

[TABLE]

On the other hand, the $7.1$ keV $N_{1}$ populates the warm DM region Merle et al. (2014); Adulpravitchai and Schmidt (2015b), and the explanation of the $3.5\leavevmode\nobreak\ {\rm keV}$ X-ray line via the radiative decay $N_{1}\to\nu\gamma$ can also be consistent with the Lyman-$\alpha$ observation Merle and Schneider (2015).

In Fig. 1, the current DM relic density, $\Omega_{N_{1}}h^{2}=0.12$ Tanabashi et al. (2018), is fitted with a numerical hierarchy $\Lambda_{\nu 1}=10^{5}\Lambda_{\nu 2(3)}\leavevmode\nobreak\ {\rm MeV}$ (blue line). Here the choice of the hierarchy between $\Lambda_{\nu 1}$ and $\Lambda_{\nu 2(3)}$, as well as the approximation $\Lambda_{\nu 2}\simeq\Lambda_{\nu 3}$ is sensible, as $\Lambda_{\nu 1}$ corresponds to the DM interactions and the small difference $\Lambda_{\nu 2}-\Lambda_{\nu 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}\to 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 Hernandez et al. (2014).

$R_{K^{(*)}}$: The dominant NP contribution to $b\to s\ell^{+}\ell^{-}$ stems from the box diagrams mediated by charged Higgs boson, top quark and heavy neutrino, which result in the effective Wilson coefficients:

[TABLE]

corresponding to the semi-leptonic operators

[TABLE]

Here, $\lambda_{\pm}=M_{H^{\pm}}^{2}/M_{W}^{2}$, $\lambda_{i}=M_{i}^{2}/M_{W}^{2}$, $\lambda_{t}=m_{t}^{2}/M_{W}^{2}$, and the loop function

[TABLE]

It is seen that, due to the neutrino mass hierarchies encoded in the mass-powered parametrization (see Eq. (21)), the necessary flavor non-universal couplings appear naturally, and the LFUV observables $R_{K^{(*)}}$ can be explained by the fact that $M_{2(3)}\gg 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 $\mu\to e\gamma$, $\tau\to\mu\gamma$, $Z\to\mu^{+}\mu^{-}$, as well as the ratio $R^{W}_{\mu e}\equiv\Gamma(W\to\mu\nu)/\Gamma(W\to e\nu)$ that tests the lepton-flavor universality Tanabashi et al. (2018). Finally, the perturbative unitarity requirement $\mathcal{Y}>\sqrt{4\pi}$ is also imposed.

From Fig. 1, it is clearly seen that the updated model-independent analyses Aebischer et al. (2019); Ciuchini et al. (2019); Algueró et al. (2019); Arbey et al. (2019) that prefer $C_{9\mu}^{\rm NP}=-C_{10\mu}^{\rm NP}<0$ with a significance at the level of $\sim 4\sigma$ 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. (Probing New Physics Signals with Symmetry-Restored Yukawa Textures)). Note that, for the sake of plotting, we have compiled the possibly negative values of $\mathcal{Y}_{\nu}$ to that of $\Lambda_{\nu 2}$, since the considered observables depend on $|\mathcal{Y}_{\nu}|^{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 $\tau\to\mu\gamma$ constraint.

Furthermore, we have also numerically confirmed that, if $\mathcal{Y}_{\nu}$ exhibits a texture similar to that of the Dirac fermions,

[TABLE]

rather than the one proposed in Eq. (21), the resolution of the $R_{K^{(*)}}$ anomalies would be excluded either by $\mu\to e\gamma$ or by $\tau\to\mu\gamma$, no matter which hierarchies of the active neutrino masses are taken. Therefore, being in association with the lightest heavy neutrino as a $7.1\leavevmode\nobreak\ {\rm keV}$ DM, the mass-powered texture of $\mathcal{Y}_{\nu}$ in Eq. (21) 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\to c\ell_{i}\nu$ transitions arises from the tree-level charged-Higgs contribution. Neglecting the suppressed down-quark Yukawa couplings, the effective Hamiltonian is given by

[TABLE]

It can be readily seen that, the lepton-flavor non-universality is induced by mass hierarchy encoded in $\mathcal{Y}_{\ell}$. For numerical analysis, we implement the Wilson coefficient to the updated formulae in Ref. Blanke et al. (2019a). Besides, to verify that our model can provide a feasible resolution of the $R_{D^{(*)}}$ anomalies, we consider the bounds from the branching ratio $\mathcal{B}(B_{c}\to\tau\nu)$ derived from the $B_{c}$ lifetime Alonso et al. (2017); Celis et al. (2017); Li et al. (2016), and the $D^{*}$ longitudinal polarization fraction Abdesselam et al. (2019c),

[TABLE]

which differs from its SM prediction Hu et al. (2019),

[TABLE]

by $\sim 1.6\sigma$. Following Ref. Li et al. (2019), we also consider the mass difference $\Delta M_{s}$ in the $B_{s}-\bar{B}_{s}$ system and the inclusive $B\to X_{s}\gamma$ branching ratio, both of which receive sizable NP effects from the top-associated Yukawa couplings.

As shown in Fig. 2, a $1\sigma$-level explanation of $R_{D^{(*)}}$ anomalies is feasible under the constraint $\mathcal{B}(B_{c}\to\tau\nu)<60\%$ Blanke et al. (2019a, b), while a simultaneous explanation can be realized only at $2\sigma$ level if $\mathcal{B}(B_{c}\to\tau\nu)<30\%$ Alonso et al. (2017) is imposed. Here the $2\sigma$ range allowed by $\Delta M_{s}^{\rm exp}/\Delta M_{s}^{\rm SM}$ (green region) is obtained with $\Delta M_{s}^{\rm exp}=(17.757\pm 0.021)\leavevmode\nobreak\ \rm ps^{-1}$ Amhis et al. (2019b) and $\Delta M_{s}^{\rm SM}=(20.01\pm 1.25)\leavevmode\nobreak\ \rm ps^{-1}$ Di Luzio et al. (2018); Aoki et al. (2019). On the other hand, by taking $n_{\ell}\simeq 1$, the charged-lepton mass hierarchies encoded in $\mathcal{Y}_{\ell}$ would induce significant effects on the $\tau$ channel but suppress those in the $\mu/e$ modes. This ensures negligible effects on $\mathcal{B}(B\rightarrow D^{(\ast)}\mu\bar{\nu})/\mathcal{B}(B\rightarrow D^{(\ast)}e\bar{\nu})$ Jung and Straub (2019) and, at the same time, provides a natural explanation of the $R_{D^{(*)}}$ anomalies.

Concerning the constraint from $H^{\pm}\to\tau\nu$ searches at the LHC, it is found that the decay width of $H^{\pm}$ is now dominated by $\mathcal{B}(H^{\pm}\to\tau\nu)$ ($30\%$) and $\mathcal{B}(H^{\pm}\to tb)$ ($70\%$), and the strength of cross section times branching ratio, $\sigma(pp\to H^{\pm})\times\mathcal{B}(H^{\pm}\to\tau\nu)$, will be reduced compared to the case with $\mathcal{B}(H^{\pm}\to\tau\nu)\simeq 90\%$. Therefore, the compatibility between the $R_{D^{(*)}}$ explanation and the LHC constraint from $H^{\pm}\to\tau\nu$ searches Iguro et al. (2019) can be realized, as the latter becomes weaker by at least a factor of three.

$\Delta a_{\mu}$: As a further application of the mass-powered parametrization, let us now consider the longstanding puzzle observed in the muon $g-2$, $a_{\mu}=(g_{\mu}-2)/2$. The current value Keshavarzi et al. (2018),

[TABLE]

exhibits a $3.7\sigma$ 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\sigma$ deviations Jegerlehner and Nyffeler (2009). 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$ Borsanyi et al. (2018); Blum et al. (2018); Davier et al. (2019).

Here, the NP contributions arise from both the one-loop 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 Li et al. (2019). To sum over these contributions, we have taken the formulae Li et al. (2019) by rescaling the Yukawa couplings with the mass-powered parametrization given in Eq. (21). 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 $\tau\to 3\mu$. It is seen that the mass-powered parametrization cannot explain the current $3.7\sigma$ discrepancy of the muon $g-2$ at $1\sigma$ level. Instead, the $3\sigma$-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 Borsanyi et al. (2018); Blum et al. (2018); Davier et al. (2019).

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}\times 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\leavevmode\nobreak\ {\rm keV}$ sterile neutrino DM. Such a DM can explain the $3.5\leavevmode\nobreak\ {\rm keV}$ X-ray line, and its relic density is produced by the freeze-in mechanism via the thermalized scalar decays. In light of the mass-powered parametrization, the IO pattern of the active neutrino masses, $0\simeq m_{3}\ll m_{1}<m_{2}$, together with a strong hierarchical heavy neutrino masses, $M_{1}\ll M_{2}\simeq 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.

This work is supported by the National Natural Science Foundation of China under Grant Nos. 11675061 and 11435003, as well as by the Fundamental Research Funds for the Central Universities under Grant Nos. CCNU18TS029 and 2019YBZZ079.