Flavor Non-universality Gauge Interactions and Anomalies in B-Meson Decays

TL;DR

This paper systematically constructs anomaly-free $U(1)'$ gauge theories addressing flavor non-universality and anomalies in B-meson decays, proposing models that induce flavor-changing interactions and are testable via collider and meson mixing experiments.

Contribution

It provides a general framework for building anomaly-free non-universal $U(1)'$ models with specific charge assignments and explores their implications for B-meson decay anomalies.

Findings

Derived a condition for anomaly cancellation with non-universal charges.

Identified five independent solutions for charge assignments.

Discussed experimental probes like meson mixing and collider searches.

Abstract

Motivated by flavor non-universality and anomalies in semi-leptonic B-meson decays, we present a general and systematic discussion about how to construct anomaly-free $U(1)'$ gauge theories based on an extended standard model with only three right-handed neutrinos. If all standard model fermions are vector-like under this new gauge symmetry, the most general family non-universal charge assignments, $(a,b,c)$ for three-generation quarks and $(d,e,f)$ for leptons, need satisfy just one condition to be anomaly-free, $3(a+b+c)=-(d+e+f)$. Any assignment can be linear combinations of five independent anomaly-free solutions. We also illustrate how such models can generally lead to flavor-changing interactions and easily resolve the anomalies in B-meson decays. Probes with $B_s-\bar{B}_s$ mixing, decay into $\tau^\pm$, dilepton and dijet searches at colliders are also discussed.

Figures

Tables

Table 1. Table 1 : An example with anomaly-free realization of U ​ ( 1 ) ′ 𝑈 superscript 1 ′ U(1)^{\prime} charges for SM chiral fermions. The last column shows the U ​ ( 1 ) ′ 𝑈 superscript 1 ′ U(1)^{\prime} charges for three generations, where a , b , c , d , e 𝑎 𝑏 𝑐 𝑑 𝑒 a,b,c,d,e and f 𝑓 f are arbitrary real numbers. The U ​ ( 1 ) ′ 𝑈 superscript 1 ′ U(1)^{\prime} gauge anomaly is canceled in the quark and lepton sectors separately if c = − ( a + b ) 𝑐 𝑎 𝑏 c=-(a+b) and f = − ( d + e ) 𝑓 𝑑 𝑒 f=-(d+e) , otherwise Eq. ( 11 ) has to hold.

	$SU{\left(3\right)}_c$	$SU{\left(2\right)}_L$	$U{\left(1\right)}_Y$	$U{\left(1\right)}^{\prime }$
$Q_L^i$	$3$	$2$	$+1/3$	$z_{Q^i}=(a,b,c)$
$u_R^{i\ast }$	$\overline{3}$	$1$	$-4/3$	$z_{u_i^{\ast }}=(-a,-b,-c)$
$d_R^{i\ast }$	$\overline{3}$	$1$	$+2/3$	$z_{d_i^{\ast }}=(-a,-b,-c)$
$L^i$	$1$	$2$	$-1$	$z_{L^i}=(d,e,f)$
$e_R^{i\ast }$	$1$	$1$	$+2$	$z_{e_i^{\ast }}=(-d,-e,-f)$
${\nu }_R^{i\ast }$	$1$	$1$	$0$	$z_{{\nu }_i^{\ast }}=(-d,-e,-f)$

Equations

$${\mathcal{R}}_K=\frac{Br\left(B^{+}\to K^{+}{\mu }^{+}{\mu }^{-}\right)}{Br\left(B^{+}\to K^{+}{e}^{+}{e}^{-}\right)},{\mathcal{R}}_{K^{\ast }}=\frac{Br\left(B\to K^{\ast }{\mu }^{+}{\mu }^{-}\right)}{Br\left(B\to K^{\ast }{e}^{+}{e}^{-}\right)}.$$

$${\mathcal{R}}_K$$

$${\mathcal{R}}_{K^{\ast }}$$

$$Q_{L}^{i}=\left(\begin{array}[]{c}u_{L}^{i}\\ d_{L}^{i}\end{array}\right),\;L^{i}=\left(\begin{array}[]{c}\nu_{L}^{i}\\ e_{L}^{i}\end{array}\right),\;\psi_{L,R}^{i}=P_{L,R}\psi^{i},$$

$$D_{\alpha \beta \gamma }\equiv \text{tr}\left[\left\{T_{\alpha },T_{\beta }\right\}{T}_{\gamma }\right]$$

$$0$$

$$z_{Q_i}=-z_{u_i^{\ast }}=-z_{d_i^{\ast }},z_{L_i}=-z_{e_i^{\ast }}=-z_{{\nu }_i^{\ast }}.$$

$$3\sum _{i=1}^3{z}_{Q_i}=-\sum _{i=1}^3{z}_{L_i}.$$

$$z_Q=-z_L/3,$$

$$3\sum _{i=1}^3{z}_{Q_i}=-\sum _{i=1}^3{z}_{L_i},$$

$$3(a+b+c)=-(d+e+f).$$

$$L_e-L_{\mu },L_{\mu }-L_{\tau },B_u-B_c,B_c-B_t,B-L.$$

$$3(a+b+c)+(d+e+f)=0,\sum _j{X}_j=0,\sum _j{X}_j^3=0.$$

$$z_{Q_i}=(1,1,-2).$$

$${\mathcal{L}}_{H_1}=\sum _{i,j=1}^2\left(y_{ij}^u{\bar{Q}}_{L,i}{\tilde{H}}_1{u}_{R,j}+y_{ij}^d{\bar{Q}}_{L,i}{H}_1{d}_{R,j}\right)+y_{33}^u{\bar{Q}}_{L,3}{\tilde{H}}_1{u}_{R,3}+y_{33}^d{\bar{Q}}_{L,3}{H}_1{d}_{R,3}+h.c,$$

$$\mathcal{M}_{u,d}^{H_{1}}\sim\left(\begin{array}[]{ccc}\times&\times&0\\ \times&\times&0\\ 0&0&\times\end{array}\right).$$

$${\mathcal{L}}_{H_{2/3}}$$

$$\mathcal{M}_{u,d}^{H_{2/3}}\sim\left(\begin{array}[]{ccc}0&0&\times\\ 0&0&\times\\ \times&\times&0\end{array}\right).$$

$$\mathcal{M}_{u}^{H_{2}}\sim\left(\begin{array}[]{ccc}0&0&\times\\ 0&0&\times\\ 0&0&0\end{array}\right),\mathcal{M}_{d}^{H_{2}}\sim\left(\begin{array}[]{ccc}0&0&0\\ 0&0&0\\ \times&\times&0\end{array}\right).$$

$$⟨H_i⟩={\left(0v_i/\sqrt{2}\right)}^T,i=1,2,3;⟨S⟩=v_s/\sqrt{2}.$$

$${\mathcal{L}}_H=\sum _{i=1}^3{\left(D^{\mu }{H}_i\right)}^{†}\left(D_{\mu }{H}_i\right)+{\left(D^{\mu }S\right)}^{†}\left(D_{\mu }S\right),$$

$$J_{Z^{{}^{\prime }}}^{\mu }=\sum _{\psi }\sum _{i=1}^3{\bar{\psi }}_i{\gamma }^{\mu }\left[{ϵ}_i^{{\psi }_L}{P}_L+{ϵ}_i^{{\psi }_R}{P}_R\right]{\psi }_i,\psi =u,d,e,\nu .$$

$${\psi }_R^i$$

$$J_{Z^{{}^{\prime }}}^{\mu }=\sum _{\Psi =(U,D,\mathbf{e},\mathbf{\nu })}\sum _{i,j=1}^3{\bar{\Psi }}_i{\gamma }^{\mu }\left[{\left(V_{{\Psi }_L}^{†}{ϵ}^{\psi }{V}_{{\Psi }_L}\right)}_{ij}{P}_L+{\left(V_{{\Psi }_R}^{†}{ϵ}^{\psi }{V}_{{\Psi }_R}\right)}_{ij}{P}_R\right]{\Psi }_j.$$

$${ϵ}^{\psi }=z_{{\psi }_1}I+\mathrm{diag}\left(0,z_{{\psi }_2}-z_{{\psi }_1},z_{{\psi }_3}-z_{{\psi }_1}\right)\equiv z_{{\psi }_1}I+\delta {ϵ}^{\psi },$$

$$B_{ij}^{{\psi }_{L,R}}\equiv {\left(V_{{\psi }_{L,R}}^{†}{ϵ}^{\psi }{V}_{{\psi }_{L,R}}\right)}_{ij}=z_{{\psi }_1}{\delta }_{ij}+(V_{{\psi }_{L,R}}^{†}\delta {ϵ}^{\psi }{V}_{{\psi }_{L,R}}{)}_{ij}\equiv z_{{\psi }_1}{\delta }_{ij}+\delta B_{ij}^{{\psi }_{L,R}}.$$

$$z_U=(0,0,1),z_L=(0,q_{\mu },-3-q_{\mu }).$$

$$(B-L)-(B_u-B_c)-2(B_c-B_t)+(L_e-L_{\mu })+(q_{\mu }+2)(L_{\mu }-L_{\tau }),$$

Full text

I Introduction

Recently, some intriguing anomalies have been found in the branching ratios of semi-leptonic decays of B-mesons into electron and muon pairs,

[TABLE]

The LHCb collaboration presented the following values in Refs. Aaij:2014ora ; Aaij:2017vbb ; RKstar :

[TABLE]

where $q^{2}$ is the invariant mass for the final lepton pair. However, the standard model (SM) predicts $\mathcal{R}_{K}\approx 1\approx\mathcal{R}_{K^{*}}$ Hiller:2003js ; Bordone:2016gaq in the above kinematic region. It has been claimed that the overall deviation from the SM in B-physics is more than $4\sigma$ when global analyses are performed Capdevila:2017bsm ; DAmico:2017mtc ; Altmannshofer:2017yso by including other anomalies (the branching ratios of $B\rightarrow K^{(*)}\mu^{+}\mu^{-}$ Aaij:2014pli and $B_{s}\rightarrow\phi\mu^{+}\mu^{-}$ Aaij:2015esa , the angular distribution of decay rate of $B\rightarrow K^{(*)}\mu^{+}\mu^{-}$ and $P_{5}^{\prime}$ observables Aaij:2015oid ; Aaij:2014pli ; Aaij:2013qta ). If these anomalies are truly due to some physical effects, then lepton flavour universality (LFU) is violated and new physics beyond the SM is warranted.

We shall focus exclusively on anomalies in $\mathcal{R}_{K}$ and $\mathcal{R}_{K^{*}}$, since the theoretical uncertainty is expected to be small. Economical explanations are involved with four-fermion effective operators, such as $(\bar{s}\gamma_{\mu}P_{L}b)(\bar{\ell}\gamma^{\mu}\ell)$, see Ref. Alonso:2014csa for more systematic discussion. More concrete models have also been constructed to generate such an effective operator Alonso:2017uky ; Sala:2017ihs ; Bishara:2017pje ; Ellis:2017nrp ; Bonilla:2017lsq ; Feruglio:2017rjo ; Greljo:2017vvb ; Alonso:2017bff ; Wang:2017mrd ; Alok:2017sui ; Alok:2017jaf ; DiChiara:2017cjq ; Kamenik:2017tnu ; Cai:2017wry ; Ghosh:2017ber ; Becirevic:2017jtw ; Celis:2017doq . Various models Gauld:2013qja ; Buras:2013dea ; Altmannshofer:2014cfa ; Crivellin:2015mga ; Crivellin:2015lwa ; Celis:2015ara ; Greljo:2015mma ; Altmannshofer:2015mqa ; Niehoff:2015bfa ; Belanger:2015nma ; Falkowski:2015zwa ; Carmona:2015ena ; Chiang:2016qov ; Becirevic:2016zri ; Boucenna:2016wpr ; Megias:2016bde ; GarciaGarcia:2016nvr ; Ko:2017lzd ; Megias:2017ove ; Hiller:2014yaa ; Gripaios:2014tna ; Crivellin:2017zlb ; Alonso:2015sja ; Sahoo:2016pet ; Hu:2016gpe ; Chen:2017hir ; Ko:2017yrd ; Geng:2017svp ; Altmannshofer:2017poe in the literature, including extra $Z^{\prime}$, lepto-quark and loop-induced mechanisms, were proposed to address similar issues in the past.

In this paper we focus on $Z^{\prime}$ models with an extra $U(1)$ gauge symmetry. In the literature, usually just a specific charge assignment is chosen, without noting that many other options could be equally possible. Here, we provide a systematic investigation of general family $U(1)$ gauge symmetry and illustrate how to choose charges consistently to get anomaly-free models without introducing new fermions, except for three right-handed neutrinos. For family universal models, there is only one non-trivial charge assignment, the well-known $B-L$ symmetry. For family non-universal models, however, infinitely many solutions exist as linear combinations of five independent anomaly-free bases. We also show how some models can provide explanations for the anomalies in B-meson decays.

This paper is organized as follows. In Section II, we discuss the consistent conditions for $U(1)^{\prime}$ charge assignment first, then give an example to show how a realistic model can be constructed to match the observed fermion masses and mixings. In Section III, we exemplify one charge assignment in the context of anomalies in B-meson decays. Finally, we give our conclusion.

II Anomaly free Family-Nonuniversal $U(1)^{\prime}$ Models

In this section, we give some general discussion about anomaly-free conditions for $U(1)^{\prime}$ models, without introducing extra chiral fermions other than three right-handed neutrinos. We denote the weak doublets and singlets as follows, $\psi=u,d,e,\nu$,

[TABLE]

with $P_{L}=(1-\gamma_{5})/2,P_{R}=(1+\gamma_{5})/2$ and $i=1,2,3$ as the family/generation index. The anomaly is proportional to the completely symmetric constant factor,

[TABLE]

$T_{\alpha}$ is the representation of the gauge algebra on the set of all left-handed fermion and anti-fermion fields, and $``\textrm{tr}"$ stands for summing over those fermion and anti-fermion species. Note that the $T$s above may or may not be the same since they depend on the referred gauge groups and also the chiral fermions running in the loop of the triangle anomaly-diagram.

The anomaly free conditions for the theory are given by

[TABLE]

So far, the discussion has been standard and the solution space of the above equations is expected to be large since we have more variables than equations. Interestingly, one can easily check that the first four equations are satisfied automatically if fermions are vector-like under the new $U(1)^{\prime}$ gauge symmetry, namely

[TABLE]

With vector-like charge assignment, we only need take care of the last two linear equations, which are actually reduced to just one,

[TABLE]

This equation is much easier to solve, but could have multiple solutions. For example,

Family universal model:

[TABLE]

which is the unique non-trivial solution, the well-known $B-L$ gauge symmetry. 2. 2.

Family non-universal models:

[TABLE]

where $z_{i}$ are not identical. Since we have six variables but just one constraint, infinitely many solutions exist. For example, we are free to choose just one generation to be charged, the other two as singlets, or any assignments for quark sector with a proper choice of charges for leptons. Some models have been discussed in Refs. Liu:2011dh ; Xing:2015fdg ; Kownacki:2016pmx ; Asai:2017ryy . In general, we can have the charge assignment as in Table 1, where $a,b,c,d,e$ and $f$ are arbitrary real numbers but satisfy

[TABLE]

As a special case, we could also imagine that anomalies are canceled separately in the quark and leptons sectors, namely $\sum z_{Q_{i}}=0=\sum z_{L_{i}}$ if $c=-(a+b)$ and $f=-(d+e)$. Such a parametrization includes some well-studied models, such as $a=b=c=0$ and $d=0,e=-f\neq 0$ corresponds to $L_{\mu}-L_{\tau}$, $d=-e\neq 0$ and $f=0$ for $L_{e}-L_{\mu}$, and so on. Note that Eq. 8 is linear, so any linear combinations of anomaly-free realizations would also satisfy this equation, like $x(B-L)+y(L_{\mu}-L_{\tau})+z(L_{e}-L_{\mu})+...$ . The solution space for Eq. (11) is five-dimensional, so we can choose the following five independent solutions as the bases,

[TABLE]

As emphasized above, we are restricting ourselves to extended models with only three additional right-handed neutrinos. If more particles are to be introduced, requirements on the charge assignment should change correspondingly. For example, one could also introduce more SM-singlet Weyl fermions $\chi_{j}$ with $U(1)^{\prime}$ charge $X_{j}$, in cases where SM fermions are vector-like in $U(1)^{\prime}$, giving

[TABLE]

Some fermion $\chi_{k}$ actually could be a dark matter (DM) candidate. For instance, a Majorana mass term $\bar{\chi}_{k}^{c}\chi_{k}$ would be induced after $U(1)^{\prime}$ symmetry breaking by a SM-singlet scalar $S$ with $U(1)^{\prime}$ charge $2X_{k}$, since interactions like $\bar{\chi}_{k}^{c}\chi_{k}S^{\dagger}$ are allowed. Vector-like $\chi_{k}$ is another popular scenario for DM where the Dirac mass term $\bar{\chi}_{k}\chi_{k}$ is allowed. In both cases, $Z_{2}$ symmetry can protect the stability of DM.

To build realistic models with correct SM fermion masses and mixings, we need to introduce some scalar fields $H_{i}$ to spontaneously break gauge symmetries. The scalar contents would be highly dependent on the charge assignments for these chiral fermions. In the most general cases, for the quark sector we can introduce several Higgs doublets with hypercharge $Y=-1$ and $U(1)^{\prime}$ charges, $a-b,a-c$ and $b-c$, to make renormalizable Yukawa interactions, giving the desired quark masses and CKM mixing matrix. In the lepton sector, Higgs doublets with $U(1)^{\prime}$ charges, $d-e,d-f$ and $e-f$, suffice to give lepton masses and neutrino mixing.

Below, we shall give an example with explicit charge assignment to illustrate how consistent models can be constructed Liu:2011dh . Let us focus on the quark sector first. We shall use the following setup:

[TABLE]

The above symmetry can be regarded as $3(B_{u}-B_{c})+6(B_{c}-B_{t})$, expanded in the five bases of Eq. 12. Some phenomenologies have been studied first in Ref. Liu:2011dh , and later in Ref. Crivellin:2015lwa along with $L_{\mu}-L_{\tau}$ symmetry in the lepton sector. Here, this model is introduced just for illustration and will be referred to in comparison with the model for the B-anomaly in Section III.

With the above $U(1)^{\prime}$ charges, a SM Higgs doublet $H_{1}$ with zero $U(1)^{\prime}$ charge can cause spontaneous electroweak symmetry breaking to generate the masses of all the SM particles, but not correct flavor mixing. To see what happens in the quark sector, we can write the gauge-invariant Yukawa terms as

[TABLE]

where $y_{ij}^{u,d}$ are the Yukawa couplings. After $H_{1}$ gets a vacuum expectation value (VEV), the resulting mass matrices for $u$ and $d$ have the following form:

[TABLE]

This kind of mass matrix cannot give the correct CKM matrix, since the third generation will not mix with the other two. Now if we have two more Higgs doublets, $H_{2}$ with $U(1)^{\prime}$ charge $-3$ and $H_{3}$ with $+3$, the following Yukawa term are allowed:

[TABLE]

When both $H_{2/3}$ get VEVs, these terms contribute to the mass matrices with

[TABLE]

Now diagonalizing the total mass matrices, $\mathcal{M}_{u,d}^{H_{1}}+\mathcal{M}_{u,d}^{H_{2/3}}$, would result in three-flavor mixing. Note that we cannot replace $\tilde{H}_{3}(H_{3})$ with $H_{2}(\tilde{H}_{2})$ in Eq. (16) because the $U(1)_{Y}$ symmetry would forbid that, although only one of them is necessary to give three-flavor mixing. In the case of no $H_{3}$ or $H_{3}$ not getting a VEV, the mass matrices are:

[TABLE]

Three-flavor mixing can still arise after diagonalization of $\mathcal{M}_{u,d}^{H_{1}}+\mathcal{M}_{u,d}^{H_{2}}$. One can easily discuss leptons, since similar physics appears. For example if $z_{L_{i}}=(0,1,-1)$, extra Higgs doublets with charges $\pm 1$ and/or $\pm 2$ would be able to achieve the required lepton masses and mixing.

Gauge bosons will get their masses through the Higgs mechanism. When $H_{2}$ and $H_{3}$ get VEVs, the $U(1)^{\prime}$ gauge symmetry is also broken. If the $U(1)^{\prime}$ gauge coupling is comparable to the electroweak coupling, the $Z^{\prime}$ boson is expected to have a mass around the electroweak scale, which is highly constrained. To get a heavy $Z^{\prime}$ boson, an electroweak singlet scalar $S$ with $U(1)^{\prime}$ charge $z_{s}$ can be introduced. Then the following vacuum configuration would break the gauge symmetries to $U(1)_{em}$,

[TABLE]

The kinetic terms for scalars are

[TABLE]

where $D_{\mu}$ is the covariant derivative. From this Lagrangian, the $W^{\pm}$ mass can be simply read out, $g_{2}\sqrt{v_{1}^{2}+v_{2}^{2}+v_{3}^{2}}/2$. Neutral gauge bosons, on the other hand, are generally mixed, but it is possible to make $Z^{\prime}$ heavy when $v_{s}\gg v_{i}$ such that experimental constraints from $Z-Z^{\prime}$ mixing can be safely evaded, since the mixing is proportional to $v^{2}_{i}/v^{2}_{s}$; see Ref. Langacker:2008yv for a general review.

The interaction for $\bar{\psi}\psi Z^{\prime}$ can be obtained from $gZ_{\mu}^{{}^{\prime}}J_{Z^{{}^{\prime}}}^{\mu}$, where $g$ is the gauge coupling constant of $U(1)^{\prime}$ and the current $J_{Z^{{}^{\prime}}}^{\mu}$ in the gauge eigenstates is given by

[TABLE]

The above $\epsilon_{i}^{\psi_{L/R}}$s are the $U(1)^{\prime}$ charges $z_{\psi_{i}}$ for fermions $\psi^{L/R}_{i}$. Rotating the fermion fields with unitary transformations such that their mass matrices are diagonalized, we get

[TABLE]

where $\Psi=U,D,\bm{e},\bm{\nu}$ are the mass eigenstates. The CKM matrix is given by $V_{\text{CKM}}=V_{U_{L}}^{\dagger}V_{D_{L}}$ and the neutrino mixing matrix by $V_{\text{PMNS}}=V_{\bm{e}_{L}}^{\dagger}V_{\bm{\nu}_{L}}$. The rotation of fermion fields in Eq.(19) leads to

[TABLE]

We have used $\epsilon^{\psi}\equiv\epsilon^{\psi_{L}}=\epsilon^{\psi_{R}}$, since we are considering the vector-like charge assignment. One can immediately notice that generally $V^{\dagger}\epsilon V\not\propto I$ if $\epsilon\not\propto I$, namely family non-universal gauge interactions. In our previous examples, we have $\epsilon^{\psi}\propto\mathrm{diag}\left(1,1,-2\right)$ or $\mathrm{diag}\left(0,1,-1\right)$, and we expect flavor-changing effects to arise. Since only $V_{\text{CKM}}$ or $V_{\text{PMNS}}$ is experimentally measured, the individual matrix $V_{\psi_{L,R}}$ is unknown. Thus the resulting products $V_{\Psi_{L,R}}^{\dagger}\epsilon^{\psi}V_{\Psi_{L,R}}$ are also unknown.

III Phenomenologies and anomalies in B-meson decays

In this section, we discuss how the above framework can address recent anomalies in B physics. Since left-handed fermions have the same charges as the right-handed ones, we can reparametrize

[TABLE]

where $z_{U}=(a,b,c),z_{L}=(d,e,f)$, and

[TABLE]

Flavor changing processes can happen when $\delta\epsilon^{\psi}\neq 0$ or $\delta B^{\psi_{L,R}}_{ij}\neq 0$. Note that elements in the matrix $\delta B^{\psi_{L,R}}$ are not necessarily smaller than $z_{\psi_{1}}$ for a general setup, since $z_{\psi_{1}}$ can be zero if fermions in the first generation are $U(1)^{\prime}$ singlets.

To illustrate how it affects B meson decay, we exemplify the following anomaly-free charge assignment

[TABLE]

This assignment can be expanded by the bases in Eq. (12),

[TABLE]

which is a nice example in the sense that it involves all five anomaly-free bases. If $q_{\mu}=-3/2$, the lepton sector has some kind of $L_{\mu}+L_{\tau}$ symmetry. If $|q_{\mu}|\ll 1$, only the third generation is effectively $U(1)^{\prime}$-charged. We should emphasize again that it is free to change the above assignment by adding any linear combinations of other anomaly-free solutions. For example, we could use $z^{\prime}_{U}=(1,1,-1)$ which is just the sum of the above charges with $(1,1,-2)$ mentioned earlier. However, these two models give different signal strengths in experiments, such as LHC dijet events, therefore they are subject to different constraints.

The $b\rightarrow s$ transitions are usually analyzed in terms of the following effective Hamiltonian

[TABLE]

Here $V$ is the CKM matrix and $\alpha=1/137$ is the fine-structure constant. Note that the coefficients $C_{i}$ and $C_{i}^{\prime}$ are scale-dependent, governed by the renormalization group equation. They are first calculated at high scales and then run to a lower scale, which is usually taken as the bottom quark mass $m_{b}$ for decay processes. We just list some relevant operators for our later discussions:

[TABLE]

In general, all the above operators can be generated. Since anomalies are closely related to $O_{9}$, we calculate the induced coefficient for $O_{9}$ by $Z^{\prime}$-mediated new physics

[TABLE]

To resolve the anomalies, $C_{9}^{\textrm{NP}}$ should be around $\simeq-1.1$ Capdevila:2017bsm ; DAmico:2017mtc , which can be translated into

[TABLE]

The above formula is generally applicable to any non-trivial charge assignment. In some cases, we can simplify it further. For instance, since $B^{\bm{e}_{L/R}}_{\mu\mu}$ are elements in the diagonal, we could expect $B^{\bm{e}_{L}}_{\mu\mu}\sim q_{\mu}$ if $|q_{\mu}|\gg|3+q_{\mu}|$, or no rotation in the charged lepton sector ($V_{\textrm{PMNS}}=V_{\bm{\nu}_{L}}$), and Eq. (27) can be approximated as

[TABLE]

Now with the charge assignment as in Eq. (23), we explicitly have $\delta B^{D_{L}}_{sb}=(V_{D_{L}}^{\dagger})_{23}(V_{D_{L}})_{33}$. If the CKM matrix comes solely from the rotation of down quarks, we would have $\delta B^{D_{L}}_{sb}=V_{tb}V_{ts}^{*}$ and

[TABLE]

Other coefficients can also be calculated similarly. Also, if $q_{\mu}\neq-3$, we would expect new physics effects to show up in $B\rightarrow K^{(*)}\tau^{+}\tau^{-}$. Since we mainly focus on $O_{9}$-related anomalies in B-meson decays to muons, we shall neglect other operators as long as the setup does not violate current limits. For example, we can freely choose $B^{\bm{e}_{L}}_{\mu\mu}=B^{\bm{e}_{R}}_{\mu\mu}$, which results in $C_{10}^{\textrm{NP}}=0=C_{10}^{{}^{\prime}\textrm{NP}}$.

$Z^{\prime}$ may also mediate $B_{s}-\bar{B}_{s}$ mixing in the above scenario, since the operator $(\bar{s}\gamma_{\mu}P_{L}b)^{2}$ is inevitably induced, which actually gives the most stringent limit at the moment. Current bounds Arnan:2016cpy can be put on the following quantity:

[TABLE]

Comparing with Eq. (29), we can safely evade this constraint for $|q_{\mu}|\gtrsim 3$ and resolve B anomalies at the same time.

In Fig. 1, motivated by B physics anomalies, we plot several contours with $C_{9}^{\textrm{NP}}\simeq-1.1$ in the $M_{Z^{\prime}}/g$ and $\delta B^{D_{L}}_{sb}$ plane for $|q_{\mu}|=1,3$ and $5$. They are shown by the dashed purple, dot-dashed red and dotted blue lines, respectively. The region above the black line is excluded by $B_{s}-\bar{B}_{s}$ mixing. As expected from Eq. (28), increasing $|q_{\mu}|$ would allow a larger parameter space. For small $|q_{\mu}|$, $Z^{\prime}$ could then be tested by other means.

Since $Z^{\prime}$ couples to both quarks and leptons, dilepton and dijet searches for heavy resonances at colliders can probe $Z^{\prime}$. The expected signal strength depends on

[TABLE]

where $\sigma$ is the cross section for $Z^{\prime}$ production, $f$ and $f^{\prime}$ are SM fermions, and $Br$ denotes the decay branching ratio. For hadron colliders, we shall integrate the above quantity over the quark parton distribution functions (PDFs) (throughout our calculations, we have used MMHT2014 Harland-Lang:2014zoa PDFs). In the case of charge assignment for quarks, $(0,0,1)$, hadron colliders such as LHC with energy $\sqrt{s}=13\ \mathrm{TeV}$ have less discovery potential for $Z^{\prime}$, since $Z^{\prime}$ would couple weakly to the first two generations through quark mixing only, but strongly to the third generation, which has small PDFs. A future $100\ \mathrm{TeV}$ hadron collider has a better chance because the production rate is increased thanks to the enhancement of the PDFs of bottom and top quarks. In Fig. 2(a), we give the ratio of $Z^{\prime}$ production from bottom and top quarks in our model to that from light quarks if $Z^{\prime}$ also couples to $u$ and $d$. We have normalized the cross section to a $3\ \mathrm{TeV}$ $Z^{\prime}$ at LHC with $\sqrt{s}=13\ \mathrm{TeV}$. As shown, $Z^{\prime}$ from the bottom channel is reduced by a factor of $\mathcal{O}(10^{3})$ at $\sqrt{s}=13\ \mathrm{TeV}$ and $\mathcal{O}(10^{2})$ at $\sqrt{s}=100\ \mathrm{TeV}$. Because of that, the limits from hadron collider searches are relaxed dramatically and $M_{Z^{\prime}}\lesssim 1\ \mathrm{TeV}$ would still be allowed, which can be inferred from Fig. 2(b), where we show dilepton searches for a Sequential SM (SSM) $Z^{\prime}$ (SSM $Z^{\prime}$ is identical to SM $Z$ except for the mass) as the dashed black line. The region above the solid red curve is excluded by dilepton searches Aaboud:2016cth . However, if the signal strength is reduced by $10$ or $100$, the exclusion limit would be shifted to $\simeq 2.4\ \mathrm{TeV}$ (dashed blue) and $1.2\ \mathrm{TeV}$ (dot-dashed purple), respectively. Since in the model discussed, the cross section is lowered by $\mathcal{O}(10^{3})$, taking the branching ratio into account would give $M_{Z^{\prime}}\gtrsim\mathcal{O}(600\ \mathrm{GeV})$, with some dependences on $q_{\mu}$. Similarly, constraints from dijets are also weakened.

In comparison, the charge assignment $(1,1,-1)$, which is the linear combination of $(1,1,-2)$ in Section II and $(0,0,1)$, will give different results. In such a case, $Z^{\prime}$ can couple to light quarks and the cross section for production can be sizable. If $g$ is at the same order as the weak coupling, the limit would be similar to the dilepton search at LHC with $\sqrt{s}=13\ \mathrm{TeV}$ for SSM $Z^{\prime}$ $M_{Z^{\prime}}\gtrsim 3.4\ \mathrm{TeV}$ Aaboud:2016cth ; Khachatryan:2016zqb and dijet channel $M_{Z^{\prime}}\gtrsim 3.4\ \mathrm{TeV}$ Aaboud:2017yvp for $g=0.5$. These limits might fluctuate, since the values of the decay branching ratio of $Z^{\prime}$ would be different from those in $Z^{\prime}_{SSM}$. In general, direct searches at colliders are complementary to the bound from $B_{s}-\bar{B}_{s}$ mixing, Eq. (30).

IV Conclusions

Motivated by the anomalies in semi-leptonic B-meson decays, we have discussed an explanation in models with general family $U(1)^{\prime}$ gauge symmetry. We have presented a systematic investigation on how to consistently assign charges to SM chiral fermions with three right-handed neutrinos. If fermions in the standard model are vector-like under this new $U(1)^{\prime}$ symmetry, their charges in Table 1 have to, and only need to, satisfy the condition given in Eq. (11). Generally, infinitely many anomaly-free family non-universal models exist, as linear combinations of five independent anomaly-free bases. If both bottom quark and muon couple to this new $U(1)$, typically the anomalies in B-meson decays can be explained.

We have also discussed several other experimental searches for such models, including $Z^{\prime}$-mediated effects in $B_{s}-\bar{B}_{s}$ mixing, dilepton and dijet searches for heavy resonances at colliders. Some viable parameter space has already been probed by these searches. Future searches in colliders and other B-meson decay modes should be able to provide more powerful information on the physical parameters and test different scenarios for $Z^{\prime}$ charge assignment.

Acknowledgement

YT is grateful to Koichi Hamaguchi, Chengcheng Han and Kazunori Nakayama for enlightening conversations. The work of YT was supported by the Grant-in-Aid for Innovative Areas No.16H06490.