This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

A Flavorful Composite Higgs Model :
Connecting the B anomalies with the hierarchy problem

Yi Chung [email protected] Center for Quantum Mathematics and Physics (QMAP), Department of Physics,
University of California, Davis, CA 95616, U.S.A.
Abstract

We present a model which connects the neutral current B anomalies with composite Higgs models. The model is based on the minimal fundamental composite Higgs model with SU(4)/Sp(4)SU(4)/Sp(4) coset. The strong dynamics spontaneously break the symmetry and introduce five Nambu-Goldstone bosons. Four of them become the Standard Model Higgs doublet and the last one, corresponding to the broken local U(1)U(1)^{\prime} symmetry, is eaten by the gauge boson. This leads to an additional TeV-scale ZZ^{\prime} boson, which can explain the recent B anomalies. The experimental constraints and allowed parameter space are discussed in detail.

I Introduction

The Standard Model (SM) of particle physics successfully describes all known elementary particles and their interactions. However, there are still a few puzzles that have yet to be understood. One of them is the well-known hierarchy problem. With the discovery of light Higgs bosons in 2012 [1, 2], the last missing piece of the SM seemed to be filled. However, SM does not address the UV-sensitive nature of scalar bosons. The Higgs mass-squared receives quadratically divergent radiative corrections from the interactions with SM fields, which require an extremely sensitive cancellation to get a 125125 GeV Higgs boson. To avoid the large quadratic corrections, the most natural way is to invoke some new symmetry such that the quadratic contributions cancel in the symmetric limit. This requires the presence of new particles related to SM particles by the new symmetry.

One appealing solution to the hierarchy problem is the composite Higgs model (CHM), where the Higgs doublet is the pseudo-Nambu-Goldstone boson (pNGB) of a spontaneously broken global symmetry of the underlying strong dynamics [3, 4]. Through the analogy to the chiral symmetry breaking in quantum chromodynamics (QCD), which naturally introduces light scalar fields, i.e., pions, we can construct models with light Higgs bosons in a similar way. In a CHM, an approximate global symmetry GG is spontaneously broken by some strong dynamics down to a subgroup HH at a symmetry breaking scale ff. The heavy resonances of the strong dynamics are expected to be around the compositeness scale 4πf\sim 4\pi f generically. The pNGBs of the symmetry breaking, on the other hand, can naturally be light with masses <f<f as they are protected by the shift symmetry.

Among all types of CHMs with different cosets, the CHMs with fundamental gauge dynamics featuring only fermionic matter fields are of interest in many studies [5, 6, 7, 8], which is known as the fundamental composite Higgs model (FCHM). In this type of CHMs, hyperfermions ψ\psi are introduced as the representation of hypercolor (HC) group GHCG_{HC}. Once the HC group becomes strongly coupled, hyperfermions form a condensate, which breaks the global symmetry. However, they always introduce more than four pNGBs, which means more light states are expected to be found. The minimal FCHM, which is based on the SU(4)/Sp(4)SU(4)/Sp(4) coset [9, 10, 11], contains five pNGBs. The four of them formes the SM Higgs doublet, and the fifth one, as a SM singlet, could be a light scalar boson (if the symmetry is global) or a TeV-scale ZZ^{\prime} boson (if the symmetry is local). No matter which, it should lead to some deviations in low energy phenomenology.

Although the direct searches by ATLAS and CMS haven’t got any evidence of new particles, LHCb, which does the precise measurement of B meson properties, shows interesting hints of new physics. There are discrepancies in several measurements of semileptonic B meson decays, especially the tests of lepton flavor universality (LFU), which are so-called the neutral current B anomalies [12, 13, 14, 15, 16, 17, 18]. Each anomaly is not statistically significant enough to reach the discovery level, but the combined analysis shows a consistent deviation from the SM prediction [19, 20, 21, 22, 23, 24]. These anomalies might be the deviation we are looking for.

One of the popular explanations is through a new ZZ^{\prime} vector boson which has flavor-dependent interactions with SM fermions. Many different types of ZZ^{\prime} models with diverse origins of U(1)U(1)^{\prime} gauge symmetry have been proposed [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52]. Depending on its couplings with fermions, the mass of the ZZ^{\prime} can range from sub-TeV to multi-TeV. For a ZZ^{\prime} boson at the TeV scale, it is natural to try to connect it with the hierarchy problem 111For our interest, we would like to mention some researches aiming at explaining the B anomalies within composite Higgs models. Different studies using different features of composite theory to address the problem, such as additional composite leptoquarks [54, 55, 57, 58] or composite vector resonances [59, 60, 61, 62, 63, 64, 65]. However, they are all different from this study, where we introduce a new ZZ^{\prime} boson..

In this paper, we realize this idea using a SU(4)/Sp(4)SU(4)/Sp(4) FCHM, where an U(1)U(1)^{\prime} subgroup within SU(4)SU(4) is gauged. The corresponding ZZ^{\prime} boson only couples to the third generation SM fermions F3F_{3} and the hyperfermions ψ\psi through the terms

int=gZZμ(F¯3γμF3+QHCψ¯γμψ),\displaystyle\mathcal{L}_{\text{int}}=g_{Z^{\prime}}Z^{\prime}_{\mu}\,(\,\bar{F}_{3}\gamma^{\mu}F_{3}+Q_{HC}\bar{\psi}\gamma^{\mu}\psi\,), (1)

where gZg_{Z^{\prime}} was normalized such that SM fermions F3F_{3} carry a unit charge and hyperfermions carry charge QHCQ_{HC}. When the hypercolor group becomes strongly coupled, the global symmetry SU(4)SU(4) and its gauged U(1)U(1)^{\prime} subgroup are broken. The 5th5^{th} pNGB is eaten by the U(1)U(1)^{\prime} gauge boson, which results in a TeV-scale ZZ^{\prime} boson. We will test the potential for this ZZ^{\prime} boson to explain the neutral current B anomalies. The parameter space allowed by different experimental constraints, mainly from neutral meson mixings and lepton flavor violation decays, will be discussed. The bounds on MZM_{Z^{\prime}} from the LHC direct searches are also shown.

This paper is organized as follows. In section II, we introduce the SU(4)/Sp(4)SU(4)/Sp(4) FCHM. The calculations of the gauge sector, including SM gauge group and U(1)U(1)^{\prime} gauge symmetry, are presented. To study its phenomenology, we specify the transformation between flavor basis and mass basis in section III. The resulting low energy phenomenology is discussed in section IV, including the B anomalies and other experimental constraints. Section V focuses on the direct searches, which play an important role in constraining a TeV-scale ZZ^{\prime} boson. Section VI and Section VII contains our discussions and conclusions.

II The SU(4)/Sp(4)SU(4)/Sp(4) FCHM

In fundamental composite Higgs models, additional hyperfermions ψ\psi are added to generate composite Higgs. The hyperfermions are representations of hypercolor group GHCG_{HC}, whose coupling becomes strong around the TeV scale. The hyperfermions then form a condensate, which breaks the global symmetry and results in the pNGBs as the Higgs doublet. In this paper, we study the minimal fundamental composite Higgs model based on the global symmetry breaking SU(4)Sp(4)SU(4)\to Sp(4). The fermionic UV completion of a SU(4)/Sp(4)SU(4)/Sp(4) FCHM only require four Weyl fermions in the fundamental representation of the SU(2)=Sp(2)SU(2)=Sp(2) hypercolor group [7, 8]. The four Weyl fermions transform under GSM=SU(3)C×SU(2)L×U(1)YG_{SM}=SU(3)_{C}\times SU(2)_{L}\times U(1)_{Y} as

ψL=(UL,DL)=(1,2,0),\displaystyle\psi_{L}=(U_{L},D_{L})=(1,2,0),
UR=(1,1,1/2),DR=(1,1,1/2).\displaystyle U_{R}=(1,1,1/2),\quad D_{R}=(1,1,-1/2). (2)

Next, we rewrite the two right-handed hyperfermions as U~L=iσ2CU¯RT\tilde{U}_{L}=-i\sigma^{2}C\bar{U}^{T}_{R} and D~L=iσ2CD¯RT\tilde{D}_{L}=-i\sigma^{2}C\bar{D}^{T}_{R}. Since all the four Weyl fermions are according to the same representation of the hypercolor group, we can recast them together as

ψ=(UL,DL,U~L,D~L)T,\psi=(U_{L},D_{L},\tilde{U}_{L},\tilde{D}_{L})^{T}~{}, (3)

which has a SU(4)SU(4) global symmetry (partially gauged). The hypercolor group becomes strongly coupled at the TeV scale, which forms a non-perturbative vacuum and breaks the SU(4)SU(4) down to Sp(4)Sp(4). In CHMs, the condensate ψψΣ0\langle\psi\psi\rangle\propto\Sigma_{0} is chosen such that electroweak symmetry is preserved. It will be broken after the Higgs interactions and loop-induced potentials are taken into account. However, we will only focus on some key ingredients here and leave the complete analysis to the future.

II.1 Basics of SU(4)/Sp(4)SU(4)/Sp(4)

To study the SU(4)/Sp(4)SU(4)/Sp(4) symmetry breaking, we can parametrize it by a non-linear sigma model. Consider a sigma field Σ\Sigma, which transforms as an anti-symmetric tensor representation 𝟔\mathbf{6} of SU(4)SU(4). The transformation can be expressed as ΣgΣgT\Sigma\to g\,\Sigma\,g^{T} with gSU(4)g\in SU(4). The scalar field Σ\Sigma has an anti-symmetric VEV Σ\langle\Sigma\rangle, where

Σ=Σ0=(iσ200iσ2).\langle\Sigma\rangle=\Sigma_{0}=\begin{pmatrix}i\sigma_{2}&0\\ 0&i\sigma_{2}\\ \end{pmatrix}. (4)

The Σ\Sigma VEV breaks SU(4)SU(4) down to Sp(4)Sp(4), producing five Nambu-Goldstone bosons.

The 15 SU(4)SU(4) generators can be divided into the unbroken ones and broken ones with each type satisfying

{unbroken generatorsTa:TaΣ0+Σ0TaT=0,broken generatorsXa:XaΣ0Σ0XaT=0.\begin{cases}\text{unbroken generators}&T_{a}:T_{a}\Sigma_{0}+\Sigma_{0}T_{a}^{T}=0~{},\\ \text{broken generators}&X_{a}:X_{a}\Sigma_{0}-\Sigma_{0}X_{a}^{T}=0~{}.\end{cases} (5)

The Nambu-Goldstone fields can be written as a matrix with the broken generator:

ξ(x)eiπa(x)Xa2f.\xi(x)\equiv e^{\frac{i\pi_{a}(x)X_{a}}{2f}}. (6)

Under SU(4)SU(4), the ξ\xi field transforms as ξgξh\xi\to g\,\xi\,h^{\dagger} where gSU(4)g\in SU(4) and hSp(4)h\in Sp(4). The relation between ξ\xi and Σ\Sigma field is given by

Σ(x)=ξΣ0ξT=eiπa(x)XafΣ0.\Sigma(x)=\xi\,\Sigma_{0}\,\xi^{T}=e^{\frac{i\pi_{a}(x)X_{a}}{f}}\Sigma_{0}~{}. (7)

The broken generators and the corresponding fields in the matrix can be organized as follows:

iπaXa=\displaystyle i\pi_{a}X_{a}= (ia𝕀2(H~H)2(H~H)ia𝕀)\displaystyle\begin{pmatrix}{ia}\,\mathbb{I}&\sqrt{2}\left(\tilde{H}H\right)\\ -\sqrt{2}\left(\tilde{H}H\right)^{\dagger}&-ia\,\mathbb{I}\\ \end{pmatrix} (8)

In this matrix, there are five independent fields. The four of them form the Higgs (complex) doublet HH. Besides, there is one more singlet aa, which will turn out to be the longitudinal part of the ZZ^{\prime} boson. By these matrices, we can construct the low energy effective Lagrangian for these pNGB fields.

II.2 The SM gauge sector

The SM electroweak gauge group SU(2)W×U(1)YSU(2)_{W}\times U(1)_{Y} is embedded in SU(4)×U(1)XSU(4)\times U(1)_{X} with generators given by

SU(2)W:12(σa000),U(1)Y:12(0000000000100001)+X𝐈.SU(2)_{W}:\frac{1}{2}\begin{pmatrix}\sigma^{a}&0\\ 0&0\\ \end{pmatrix},\quad U(1)_{Y}:\frac{1}{2}\begin{pmatrix}0&0&0&0\\ 0&0&0&0\\ 0&0&-1&0\\ 0&0&0&1\\ \end{pmatrix}+X\mathbf{I}~{}. (9)

The extra U(1)XU(1)_{X} factor accounts for the different hypercharges of the fermion representations but is not relevant for the bosonic fields. These generators belong to Sp(4)×U(1)XSp(4)\times U(1)_{X} and are not broken by Σ0\Sigma_{0}. Using the Σ\Sigma field, the Lagrangian for kinetic terms of Higgs boson comes from

h=f28tr[(DμΣ)(DμΣ)]+,\mathcal{L}_{h}=\frac{f^{2}}{8}\text{tr}\left[(D_{\mu}\Sigma)(D^{\mu}\Sigma)^{\dagger}\right]+\cdots, (10)

where DμD_{\mu} is the electroweak covariant derivative. Expanding this, we get

h=12(μh)2+f28gW2sin2(hf)[2Wμ+Wμ+ZμZμcosθW].\mathcal{L}_{h}=\frac{1}{2}(\partial_{\mu}h)^{2}+\frac{f^{2}}{8}g_{W}^{2}\,\text{sin}^{2}\left(\frac{h}{f}\right)\left[2W^{+}_{\mu}W^{-\mu}+\frac{Z_{\mu}Z^{\mu}}{\text{cos}\,\theta_{W}}\right]. (11)

The non-linear behavior of the Higgs boson in the CHM is apparent from the dependence of trigonometric functions. When hh obtains a nonzero VEV h=V\langle h\rangle=V, the WW boson acquires a mass of

mW2=f24gW2sin2(Vf)=14gW2v2,m_{W}^{2}=\frac{f^{2}}{4}g_{W}^{2}\,\text{sin}^{2}\left(\frac{V}{f}\right)=\frac{1}{4}g_{W}^{2}v^{2}, (12)

where vfsin(V/f)Vv\equiv f\,\text{sin}(V/f)\approx V. The non-linearity of the CHM is parametrized by

ξv2f2=sin2(Vf).\xi\equiv\frac{v^{2}}{f^{2}}=\sin^{2}\left(\frac{V}{f}\right)~{}. (13)

The Higgs boson couplings to SM fields in the SU(4)/Sp(4)SU(4)/Sp(4) CHM are modified by the non-linear effect due to the pNGB nature of the Higgs boson. For example, the deviation of the Higgs coupling to vector bosons is parameterized by

κVghVVghVVSM=cos(Vf)=1ξ1ξ2.\kappa_{V}\equiv\frac{g_{hVV}}{g^{SM}_{hVV}}=\text{cos}\left(\frac{V}{f}\right)=\sqrt{1-\xi}\approx 1-\frac{\xi}{2}~{}. (14)

To decide the bound on the parameter ξ\xi, we also need to determine the Yukawa coupling in the model, which is beyond the scope of the present work. The most conservative bound requires ξ0.06\xi\lesssim 0.06 [66, 67], which implies the symmetry breaking scale f1f\gtrsim 1 TeV.

II.3 U(1)U(1)^{\prime} gauge symmetry

Besides the SM gauge symmetry, we also gauge the U(1)U(1)^{\prime} subgroup of SU(4)SU(4) with the generator given by

U(1):QHC(𝕀00𝕀).U(1)^{\prime}:~{}Q_{HC}\begin{pmatrix}\mathbb{I}&0\\ 0&-\mathbb{I}\\ \end{pmatrix}. (15)

The U(1)U(1)^{\prime} behaves like the lepton number of hyperfermions, where a hyperfermion carry charge QHCQ_{HC} and an anti-hyperfermion carry charge QHC-Q_{HC}. To explain the neutral current B anomalies without violating the experimental constraints, we assume SM fermions (but only the third generation) also carry a nonzero, universal charge, which is set to 1 for simplicity as mentioned in eq. (1). To make the U(1)U(1)^{\prime} gauge symmetry anomaly-free, we need to take QHC=2Q_{HC}=-2 in the minimal FCHM. Now the U(1)U(1)^{\prime} gauge symmetry becomes the difference between the third generation SM number and the hyperfermion number, or written as SM3HFSM_{3}-HF, which is like the hyper version of anomaly-free BLB-L symmetry.

When SU(4)SU(4) global symmetry is broken down by the Σ\Sigma VEV to Sp(4)Sp(4) at the symmetry breaking scale, the U(1)U(1)^{\prime} subgroup is also broken down. It results in a massive ZZ^{\prime} gauge boson with

MZ=gZ(2|QHC|f)gZf,M_{Z^{\prime}}=g_{Z^{\prime}}\left(2\,|Q_{HC}|f\right)\equiv g_{Z^{\prime}}f^{\prime}, (16)

where we define the scale

f2|QHC|f=4f,f^{\prime}\equiv 2\,|Q_{HC}|f=4f, (17)

which is relevant in the study of ZZ^{\prime} phenomenology.

To sum up, in this flavorful SU(4)/Sp(4)SU(4)/Sp(4) FCHM, five pNGBs are generated below the compositeness scale. The four of them become the SM Higgs doublet we observed but with non-linear nature, which will be tested in the future Higgs measurements. The 5th one is eaten by the U(1)U(1)^{\prime} gauge boson and results in a heavy ZZ^{\prime} boson around the TeV scale. Other model construction issues and phenomenology of SU(4)/Sp(4)SU(4)/Sp(4) CHM have been studied comprehensively in [7, 8]. In the following sections, we will focus on the ZZ^{\prime} phenomenology and the connection with the B anomalies.

III Specify the mixing matrices for phenomenology

To discuss the phenomenology, we need to first rewrite the ZZ^{\prime} interaction terms in eq. (1) to cover all generations and separate different chirality as

int=gZZμ(F¯LfγμQFLfFLf+F¯RfγμQFRfFRf),\displaystyle\mathcal{L}_{\text{int}}=g_{Z^{\prime}}Z^{\prime}_{\mu}\,(\,\bar{F}_{L}^{f}\gamma^{\mu}Q_{F_{L}}^{f}F_{L}^{f}+\bar{F}_{R}^{f}\gamma^{\mu}Q_{F_{R}}^{f}F_{R}^{f}\,), (18)

where F=(F1,F2,F3)F=(F_{1},F_{2},F_{3}) includes SM fermions of all the three generations with superscript ff for flavor basis. The 3×33\times 3 charge matrices in the flavor basis look like

QFL/Rf=(000000001).Q_{F_{L/R}}^{f}=\begin{pmatrix}0&0&0\\ 0&0&0\\ 0&0&1\\ \end{pmatrix}. (19)

However, to study phenomenology, we need to transform them to the mass basis FL/RmF_{L/R}^{m} through the mixing matrices as FL/Rf=UFL/RFL/RmF^{f}_{L/R}=U_{F_{L/R}}F^{m}_{L/R}. After the transformation, we get

int=gZZμ(F¯LmγμQFLmFLm+F¯RmγμQFRmFRm),\displaystyle\mathcal{L}_{\text{int}}=g_{Z^{\prime}}Z^{\prime}_{\mu}\,(\,\bar{F}_{L}^{m}\gamma^{\mu}Q_{F_{L}}^{m}F_{L}^{m}+\bar{F}_{R}^{m}\gamma^{\mu}Q_{F_{R}}^{m}F_{R}^{m}\,), (20)

where the charge matrices becomes

QFL/Rm=UFL/RQFL/RfUFL/R.Q_{F_{L/R}}^{m}=U_{F_{L/R}}^{\dagger}Q_{F_{L/R}}^{f}U_{F_{L/R}}. (21)

Therefore, we need to know all the UFL/RU_{F_{L/R}} to determine the magnitude of each interaction. However, The only information about these unitary transformation matrices is the CKM matrix for quarks and PMNS matrix for leptons. The two relations that need to be satisfied are

VCKMUuLUdLandVPMNSUνLUeL,V_{CKM}\equiv U_{u_{L}}^{\dagger}U_{d_{L}}\quad\text{and}\quad V_{PMNS}\equiv U_{\nu_{L}}^{\dagger}U_{e_{L}}, (22)

which only tells us about the left-handed part with no information about the right-handed part. Even with these two constraints, they only give the difference between two unitary transformations, but not the individual one. Therefore, we need to make some assumptions about the matrices so there won’t be too many parameters.

To simplify the analysis, we assume all the UFRU_{F_{R}} are identity matrices. Therefore, for right-handed fermions, only the third generation joins in the interaction with no flavor changing at all. The couplings are the same for all the right-handed fermions it couples to with coupling strength gZg_{Z^{\prime}}.

For the left-handed side, due to the observation of VCKMV_{CKM} and VPMNSV_{PMNS}, there is a guarantee minimal transformation for UFLU_{F_{L}}. Because we only care about the transition between the second and third generation down-type quarks and charged leptons, we will only specify the rotation θ23\theta_{23} between the second and third generation of UdLU_{d_{L}} and UeLU_{e_{L}} as

UFL=(1000cosθFsinθF0sinθFcosθF)U_{F_{L}}=\begin{pmatrix}1&0&0\\ 0&\text{cos}~{}\theta_{F}&\text{sin}~{}\theta_{F}\\ 0&-\text{sin}~{}\theta_{F}&\text{cos}~{}\theta_{F}\\ \end{pmatrix} (23)

where F=d,eF=d,\,e. Keeping only the angle θ23\theta_{23} is a strong assumption but a good example case for phenomenological study because it avoids some of the most stringent flavor constraints from light fermions and leaves a simple parameter space for analysis. Following this assumption, the rest of the matrices are fixed as UuL=VCKMUdLU_{u_{L}}=V_{CKM}^{\dagger}U_{d_{L}} and UνL=VPMNSUeLU_{\nu_{L}}=V_{PMNS}^{\dagger}U_{e_{L}}. Notice that, although they looks similar, the magnitude we expect for the two angles are quite different. For θd\theta_{d}, we expect it to be CKM-like, i.e. sin θd𝒪(0.01)\theta_{d}\sim\mathcal{O}(0.01). However, for θe\theta_{e}, it could be as large as sin θe1\theta_{e}\sim 1.

We can then calculate the charge matrices as

QFLm=(0000sin2θF12sin2θF012sin2θFcos2θF),Q_{F_{L}}^{m}=\begin{pmatrix}0&0&0\\ 0&\text{sin}^{2}~{}\theta_{F}&-\frac{1}{2}\,\text{sin}~{}2\theta_{F}\\ 0&-\frac{1}{2}\,\text{sin}~{}2\theta_{F}&\text{cos}^{2}~{}\theta_{F}\\ \end{pmatrix}, (24)

where F=d,eF=d,\,e, and write down all the coupling for left-handed fermions. To study the B anomalies, two of them, gsbg_{sb} and gμμg_{\mu\mu}, are especially important, so we further define

gsbgZϵsb\displaystyle g_{sb}\equiv-g_{Z^{\prime}}\epsilon_{sb} withϵsb=12sin 2θd,\displaystyle\quad\text{with}\quad\epsilon_{sb}=\frac{1}{2}\,\text{sin}\,2\theta_{d}, (25)
gμμgZϵμμ\displaystyle\quad g_{\mu\mu}\equiv g_{Z^{\prime}}\epsilon_{\mu\mu} withϵμμ=sin2θe.\displaystyle\quad\text{with}\quad\epsilon_{\mu\mu}=\text{sin}^{2}\,\theta_{e}. (26)

We will see later that constraints will be put on the three key parameters: the scale ff^{\prime}, the mixings ϵsb\epsilon_{sb}, and ϵμμ\epsilon_{\mu\mu}.

IV Low Energy Phenomenology

With the specified mixing matrices, we can then discuss the parameter space allowed to explain the B anomalies. Also, the constraints from other low energy experiments are presented in this section.

IV.1 Neutral Current B Anomalies

To explain the observed neutral current B anomalies, an additional negative contribution on bsμ+μb\to s\mu^{+}\mu^{-} is required. Based on the assumption we make, after integrating out the ZZ^{\prime} boson, we can get the operator

Δ=4GF2VtbVtse216π2CLL(s¯LγρbL)(μ¯LγρμL)\Delta\mathcal{L}=\frac{4G_{F}}{\sqrt{2}}V_{tb}V_{ts}^{*}\frac{e^{2}}{16\pi^{2}}C_{LL}(\bar{s}_{L}\gamma^{\rho}b_{L})(\bar{\mu}_{L}\gamma_{\rho}\mu_{L}) (27)

in the low energy effective Lagrangian with coefficient

CLL=gsbgμμMZ2(35TeV)2=ϵsbϵμμf2(35TeV)2.C_{LL}=\frac{g_{sb}g_{\mu\mu}}{M_{Z^{\prime}}^{2}}~{}(35~{}\text{TeV})^{2}=-\frac{\epsilon_{sb}\epsilon_{\mu\mu}}{f^{\prime 2}}~{}(35~{}\text{TeV})^{2}. (28)

The global fit value for the Wilson coefficient, considering all rare B decays [19], gives

CLL=0.82±0.14,C_{LL}=-0.82\pm 0.14~{}, (29)

which requires

ϵsbϵμμf2=1(39TeV)2fϵsbϵμμ(39TeV).\frac{\epsilon_{sb}\epsilon_{\mu\mu}}{f^{\prime 2}}=\frac{1}{(39~{}\text{TeV})^{2}}\implies f^{\prime}\sim\sqrt{\epsilon_{sb}\epsilon_{\mu\mu}}~{}(39~{}\text{TeV}). (30)

The generic scale with large mixing angles is f40f^{\prime}\sim 40 TeV. However, as we mentioned, the value ϵsb𝒪(0.01)\epsilon_{sb}\sim\mathcal{O}(0.01), which will bring it down to the TeV scale.

IV.2 Neutral Meson Mixing

The measurement of neutral meson mixing put strong constraints on the ZZ^{\prime} solution. Based on our specified mixing matrices, which have suppressed mixings between the first two generations, the BsB¯sB_{s}-\bar{B}_{s} mixing turns out to be the strongest constraint. The measurement of mixing parameter [68] compared with SM prediction by recent lattice data [69] gives the bound on the s¯bZ\bar{s}bZ^{\prime} vertex as

gZMZϵsb1194TeVfϵsb194(TeV).\frac{g_{Z^{\prime}}}{M_{Z^{\prime}}}\epsilon_{sb}\leq\frac{1}{194~{}\text{TeV}}\implies f^{\prime}\geq\epsilon_{sb}\cdot 194~{}(\text{TeV}). (31)

Combining with the requirement from eq. (30), we can rewrite the constraint as

fϵμμ7.7(TeV).f^{\prime}\leq\epsilon_{\mu\mu}\cdot 7.7~{}(\text{TeV})~{}. (32)

The constraint can be understood as that, in the bsμ+μb\to s\mu^{+}\mu^{-} process, the bsbs side, which is constrained by the BsB¯sB_{s}-\bar{B}_{s} mixing measurement, should be extremely suppressed. Therefore, the μμ\mu\mu side needs to be large enough to generate the observed B anomalies. We can also find a hierarchy ϵμμ/ϵsb25\epsilon_{\mu\mu}/\epsilon_{sb}\geq 25, which leads to the bound ϵsb0.04\epsilon_{sb}\leq 0.04, which is consistent with what we expected.

IV.3 Lepton Flavor Violation Decay

In the lepton sector, there is also a strong constraint from the flavor changing neutral currents (FCNCs). The off-diagonal term in the charge matrix of charged lepton will introduce lepton flavor violation decay, in particular, τ3μ\tau\to 3\mu, from the effective term

LFV=gZ2MZ2se3ce(τ¯LγρμL)(μ¯LγρμL),\mathcal{L}_{LFV}=\frac{g_{Z^{\prime}}^{2}}{M_{Z^{\prime}}^{2}}s_{e}^{3}c_{e}(\bar{\tau}_{L}\gamma^{\rho}\mu_{L})(\bar{\mu}_{L}\gamma_{\rho}\mu_{L})~{}, (33)

where se=sinθes_{e}=\text{sin}\,\theta_{e} and ce=cosθec_{e}=\text{cos}\,\theta_{e}. The resulting branching ratio can be expressed as

BR(\displaystyle BR( τ3μ)=2mτ51536π3Γτ(gZ2MZ2se3ce)2\displaystyle\tau\to 3\mu)=\frac{2m_{\tau}^{5}}{1536\pi^{3}\Gamma_{\tau}}\left(\frac{g_{Z^{\prime}}^{2}}{M_{Z^{\prime}}^{2}}s_{e}^{3}c_{e}\right)^{2}
=3.28×104(1 TeVf)4ϵμμ3(1ϵμμ).\displaystyle=3.28\times 10^{-4}\,\left(\frac{1\text{ TeV}}{f^{\prime}}\right)^{4}\epsilon_{\mu\mu}^{3}(1-\epsilon_{\mu\mu})~{}. (34)

The value should be <2.1×108<2.1\times 10^{-8} at 90%90\% CL by the measurement [70]. It also puts a strong constraint on the available parameter space. The exclusion plot combining the constraint from BsB¯sB_{s}-\bar{B}_{s} mixing on the parameter space ff^{\prime} v.s. ϵμμ\epsilon_{\mu\mu} is shown in Fig. 1.

Refer to caption
Figure 1: The viable parameter space from the experimental constraints. The shaded region is excluded by the corresponding measurements. The bright blue line labels the upper edge of the available parameter space.

The small ϵμμ\epsilon_{\mu\mu} region is excluded, which give a minimal value ϵμμ0.82\epsilon_{\mu\mu}\geq 0.82. It implies the angle θe\theta_{e} is quite large. The value of ff^{\prime} is bounded from above as shown in eq. (32) but not from below as it could be small in the ϵμμ=1\epsilon_{\mu\mu}=1 limit. However, due to the connection with symmetry breaking scale f1f\gtrsim 1 TeV, we are interested in f4f^{\prime}\gtrsim 4 TeV, which corresponds to the upper region of the parameter space. In this region, the ZZ^{\prime} contributions to neutrino trident production [71, 72] and muon (g2)(g-2) [73, 74] are negligible, so we will only focus on the experimental constraints we mention in this section.

V Direct ZZ^{\prime} Searches

The measurements from flavor physics in the last section can only put the constraints on the mixings and the scale f=MZ/gZf^{\prime}=M_{Z^{\prime}}/g_{Z^{\prime}}. The direct searches, on the other hand, can give the lower bound on the mass of MZM_{Z^{\prime}} directly. A general ZZ^{\prime} collider search has been discussed in [75]. In this section, we will focus on the scenario determined by our model.

V.1 Decay width and branching ratios

The partial width of the ZZ^{\prime} boson decaying into Weyl fermion pairs fi¯fj\bar{f_{i}}f_{j} is

Γij=C24πgij2MZ,\Gamma_{ij}=\frac{C}{24\pi}g_{ij}^{2}M_{Z^{\prime}}, (35)

where gijg_{ij} is the coupling of fi¯fjZ\bar{f_{i}}f_{j}Z^{\prime} vertex and CC counts the color degree of freedom. In the limit that all mfm_{f} are negligible, we get the total relative width as

ΓZMZ=1624πgZ20.2gZ2.\frac{\Gamma_{Z^{\prime}}}{M_{Z^{\prime}}}=\frac{16}{24\pi}g_{Z^{\prime}}^{2}\sim 0.2~{}g_{Z^{\prime}}^{2}~{}. (36)

The value is important when we try to pick up the bound from the LHC searches.

The dominant decay channels are the diquarks channel of the third generation quarks as

Br(tt¯)Br(bb¯)37.5%.Br(t\bar{t})\sim Br(b\bar{b})\sim 37.5\%. (37)

Decays to the light quarks and exotic decays like tctc and bsbs are also allowed but strongly suppressed due to the small rotational angles.

The main constraint is expected to come from the clear dilepton channels. Based on the specified mixing matrices we gave, the branching ratios are

Br(ττ)6.25(1+(1ϵμμ)2)%,\displaystyle Br(\tau\tau)\sim 6.25~{}(1+(1-\epsilon_{\mu\mu})^{2})~{}\%, (38)
Br(τμ)12.5ϵμμ(1ϵμμ)%,\displaystyle Br(\tau\mu)\sim 12.5~{}\epsilon_{\mu\mu}(1-\epsilon_{\mu\mu})~{}\%, (39)
Br(μμ)6.25ϵμμ2%.\displaystyle Br(\mu\mu)\sim 6.25~{}\epsilon_{\mu\mu}^{2}~{}\%. (40)

We already get ϵμμ0.82\epsilon_{\mu\mu}\geq 0.82 from the flavor constraints, which implies Br(μμ)4.2%Br(\mu\mu)\geq 4.2\%. Therefore, the μμ\mu\mu final state is the most promising channel but also puts the stringent constraint on the MZM_{Z^{\prime}}.

V.2 Production cross section

In the model, the ZZ^{\prime} boson only couples to the third generation quarks in the flavor basis. Even after rotating to the mass basis, the couplings to the first and second generation quarks are still suppressed due to the small mixing angles. Therefore, the dominant production come from the process bb¯Zb\bar{b}\to Z^{\prime}. In the following discussion, we will ignore all the other production processes and the small mixing angle θd\theta_{d}. In this way, the cross section can be written as

σ(bb¯Z)gZ2σbb(MZ){\sigma}(b\bar{b}\to Z^{\prime})\equiv g_{Z^{\prime}}^{2}\cdot\sigma_{bb}(M_{Z^{\prime}}) (41)

where the coupling dependence is taken out. The σbb\sigma_{bb} is determined by the bottom-quark parton distribution functions [76, 77], which is a function of MZM_{Z^{\prime}}.

V.3 The μμ\mu\mu channel search

From the branching ratios and the production cross section we got, we can calculate the cross section for dimuon final state

σμμσ×Br(μμ)=116σbbgZ2ϵμμ2.\sigma_{\mu\mu}\equiv\sigma\times Br(\mu\mu)=\frac{1}{16}\,\sigma_{bb}\cdot g_{Z^{\prime}}^{2}~{}\epsilon_{\mu\mu}^{2}. (42)

Moreover, from the BsB¯sB_{s}-\bar{B}_{s} constraint, we get the lower bound on ϵμμ\epsilon_{\mu\mu} as a function of ff^{\prime} in (32), which gives

σμμ116σbbgZ2(f7.7 TeV)2=σbb(MZ31 TeV)2.\sigma_{\mu\mu}\geq\frac{1}{16}\,\sigma_{bb}\cdot g_{Z^{\prime}}^{2}\left(\frac{f^{\prime}}{7.7\text{ TeV}}\right)^{2}=\sigma_{bb}\left(\frac{M_{Z^{\prime}}}{31\text{ TeV}}\right)^{2}. (43)

The equality holds when ϵμμ=f/7.7\epsilon_{\mu\mu}={f^{\prime}}/{7.7} TeV, which corresponds to the blue line in Fig. 1. It gives the minimal cross section as a function of MZM_{Z^{\prime}} that allows us to compare with the experimental results. The current best search comes from the ATLAS [78] with an integrated luminosity of 139 fb-1. The result is shown in Fig. 2.

Refer to caption
Figure 2: Upper limits at 95% CL on the cross section times branching ratio σμμ\sigma_{\mu\mu} as a function of MZM_{Z^{\prime}} for 10% (red) and 0.5% (black) relative width signals for the dimuon channel. Observed limits are shown as a solid line and expected limits as a dashed line. Also shown are theoretical predictions of the minimal cross section for ZZ^{\prime} in the model (blue) assuming CLL=0.82C_{LL}=-0.82 (solid) and 0.68-0.68 (dotted).

Notice that, the bound by collider searches depends on the width. In Fig. 2, we show relative width of 10%(red) and 0.5%(black). The wider one gives a weaker bound. However, it require a larger gZ0.7g_{Z^{\prime}}\sim 0.7 and thus a smaller f1.7f^{\prime}\sim 1.7 TeV, which is excluded as shown in Fig. 1. The bright blue segment in Fig. 1 is the available parameter space with the minimal cross section. In this region, the value f7f^{\prime}\sim 7 TeV, which implies a smaller gZ0.17g_{Z^{\prime}}\sim 0.17. Therefore, we should use the black line with 0.5% width in the plot, which requires MZ1200M_{Z^{\prime}}\gtrsim 1200 GeV. If we relax the best-fit value in the eq. (29) to one sigma region, we get a weaker bound as MZ900M_{Z^{\prime}}\gtrsim 900 GeV.

V.4 Other decay channels

To looks for other decay channels, we need to first set up benchmark points. From the previous discussion, we choose the value MZ=1.4M_{Z^{\prime}}=1.4 TeV, which is right above the current bound. For simplicity, we set ϵμμ=1\epsilon_{\mu\mu}=1, which makes σττ=σμμ\sigma_{\tau\tau}=\sigma_{\mu\mu} and στμ=0\sigma_{\tau\mu}=0. Once we pick up a value for ff^{\prime}, other parameters are automatically set. We can then calculate all the cross sections we are interested in. The results are listed in table 1. For a fixed MZM_{Z^{\prime}}, a larger ff^{\prime} implies a smaller gZg_{Z^{\prime}} and thus smaller cross sections. We can check that the σμμ\sigma_{\mu\mu} for these benchmark points are still below the bound. Other channels, even with a larger cross section, are well below the observed limits but will be tested during the HL-LHC runs.

ff^{\prime}(TeV) gZg_{Z^{\prime}} σtot\sigma_{tot}(fb) σtt/bb\sigma_{tt/bb}(fb) σττ/μμ\sigma_{\tau\tau/\mu\mu}(fb)
5.0 0.28 11.21 4.20 0.70
6.0 0.23 7.79 2.92 0.49
7.0 0.20 5.72 2.15 0.36
Table 1: The cross sections for each decay channel based on MZ=1.4M_{Z^{\prime}}=1.4 TeV with different choice of ff^{\prime}.

We only show the flavor conserving final states so far, but the ZZ^{\prime} boson can also have flavor violating decays. However, their cross sections are already constrained by the absence of FCNCs. In the quark sector, the mixings are strongly constrained and thus the branching ratios for these decays are suppressed. However, in the lepton sector, a larger mixing is allowed and the search for flavor violating decays like ZμτZ^{\prime}\to\mu\tau might be viable.

Although other channels are unlikely to be the discovery channel, once the ZZ^{\prime} boson is discovered, the next thing to do will be to look for the same resonance in other channels. Through the searches, we can decide the partial widths and figure out the couplings of the ZZ^{\prime} boson to other fields. The structure of couplings can help us distinguish between different ZZ^{\prime} models. For example, the ZZ^{\prime} boson in our model couples universally to all the third generation SM fermions in the flavor basis. Even considering the transformation to the mass basis, it still has a unique partial width ratio

Γtt:Γbb:Γ:Γνν3:3:1:1,\Gamma_{tt}:\Gamma_{bb}:\Gamma_{\ell\ell}:\Gamma_{\nu\nu}\sim 3:3:1:1, (44)

where Γ\Gamma_{\ell\ell} is the sum of all the charged lepton partial widths. The measurement will allow us to probe the nature of the ZZ^{\prime} boson and the underlying U(1)U(1)^{\prime} symmetry.

VI Discussions

In this study, we are interested in the value of ff^{\prime}, which is related to the breaking scale ff, and the bound on MZM_{Z^{\prime}}, which is important for the collider searches. In the last section, we found that a certain straight line (such as the blue line) in Fig. 1 corresponding to a predicted cross section σμμ(f0)\sigma_{\mu\mu}(f^{\prime}_{0}), which is given by

Line: ϵμμ=ff0σμμ(f0)=σbb(MZ4×f0)2,\text{Line: }\epsilon_{\mu\mu}=\frac{f^{\prime}}{f^{\prime}_{0}}\implies\sigma_{\mu\mu}(f^{\prime}_{0})=\sigma_{bb}\left(\frac{M_{Z^{\prime}}}{4\times f^{\prime}_{0}}\right)^{2}, (45)

where f0f^{\prime}_{0} represents the slope of the line, e.g. for the blue line in Fig. 1, f0=7.7f^{\prime}_{0}=7.7 TeV. Using this relation, we can calculate the cross section σμμ\sigma_{\mu\mu} for each point in the parameter space in Fig. 1 with a certain value of MZM_{Z^{\prime}}. It allows us to combine “the constraints in the parameter space in ff^{\prime} v.s. ϵμμ\epsilon_{\mu\mu} plot” (as shown in Fig. 1) with “the direct μμ\mu\mu channel search results from the ATLAS [78]” into “the viable parameter space in ff^{\prime} v.s. MZM_{Z^{\prime}} plot” as shown in Fig. 3.

Refer to caption
Figure 3: Constraints on ff^{\prime} v.s. MZM_{Z^{\prime}} plot for MZM_{Z^{\prime}} below 33 TeV. The white region is currently allowed, where ϵμμ\epsilon_{\mu\mu} and ϵsb\epsilon_{sb} are chosen to satisfy (28) from the requirement of the B anomalies. The shaded regions are excluded by the corresponding constraints from Fig. 1 combining with the direct searches, where we use the ATLAS 139 fb-1 dimuon searches. The three straight lines represent different values of gZg_{Z^{\prime}}.

The blue region is excluded by the BsB¯sB_{s}-\bar{B}_{s} meson mixing, which gives the lower bound MZ1.2M_{Z^{\prime}}\gtrsim 1.2 TeV. The bright blue line corresponds to the same parameter space as in Fig. 1 with MZ1.2M_{Z^{\prime}}\sim 1.2 TeV. The yellow region, also excluded by the BsB¯sB_{s}-\bar{B}_{s} meson mixing, sets the maximum value for ff^{\prime} as shown in eq. (32), which can also be found directly in Fig. 1. Once the stronger constraint from BsB¯sB_{s}-\bar{B}_{s} meson mixing is placed, the yellow line will move downward and the blue line will move rightward. The red region, which is excluded by τ3μ\tau\to 3\mu, restricts the parameter space from below. It places the lower bound on ff^{\prime}, which will be pushed upward if the constraint becomes stronger. We can also see the data fluctuations in dimuon search become the fluctuations on the red curve. The strength of the coupling gZg_{Z^{\prime}} with three different values is also labeled as the black straight line in the plot.

There are two regions worth noticed in the plot: (1) The region with the light ZZ^{\prime} that corresponds to a small gZg_{Z^{\prime}} but a large ff^{\prime} region, i.e. (gZ,f)(0.2,7 TeV)(g_{Z^{\prime}},f^{\prime})\sim{(0.2,7\text{ TeV})}. (2) For a natural CHM without a large fine-tuning, a smaller ff (and thus f=4ff^{\prime}=4f) is preferred, which corresponds to a larger gZg_{Z^{\prime}} region, such as (gZ,f)(0.5,4 TeV)(g_{Z^{\prime}},f^{\prime})\sim{(0.5,4\text{ TeV})} with a heavier ZZ^{\prime}. Both regions are around the boundary. The direct searches will extend both blue and red exclusion regions rightward, so both points we mentioned will be probed soon. The lower bound on MZM_{Z^{\prime}} will be pushed to 2 TeV and most of the interesting parameter space will be explored during the HL-LHC era [79, 80].

VII Conclusions

In this paper, we presented a new ZZ^{\prime} solution to the B anomalies, whose scale is related to the symmetry breaking scale of the underlying strong dynamics. We found that the anomaly-free U(1)U(1)^{\prime} symmetry can arise from SM3HFSM_{3}-HF, the difference between the third generation SM fermion number and the hyperfermion number. This type of U(1)U(1)^{\prime} is naturally broken at the TeV scale in many fundamental composite Higgs models, which allow us to connect it with the hierarchy problem. We constructed a concrete model based on SU(4)/Sp(4)SU(4)/Sp(4) minimal FCHM. The relation f=2|QHC|f=4ff^{\prime}=2\,|Q_{HC}|f=4f connects the flavor anomalies scale ff^{\prime} with the symmetry breaking scale ff in the FCHM.

The potential for the ZZ^{\prime} boson to explain the B anomalies is discussed in detail. Other flavor physics measurements, like neutral meson mixings and lepton flavor violation decays, put constraints on the allowed parameter space as shown in Fig. 1. The direct searches also give the bound on the mass of ZZ^{\prime} as MZ1.2M_{Z^{\prime}}\gtrsim 1.2 TeV. The combined constraints on the scale ff^{\prime} v.s. mass MZM_{Z^{\prime}} are shown in Fig. 3, which gives a clear picture about how the parameter space will be probed in the future. Some attractive regions are still viable and will be tested during the HL-LHC era.

Acknowledgements.
I thank Hsin-Chia Cheng for many useful discussions. I am also grateful to Ben Allanach and Wolfgang Altmannshofer for reading the previous version and giving many helpful suggestions. This work is supported by the Department of Energy Grant number DE-SC-0009999.

References