Phenomenology of flavorful composite vector bosons in light of $B$ anomalies

We analyze the flavor structure of composite vector bosons arising in a model of vectorlike technicolor, often called hypercolor (HC), with eight flavors that form a one-family content of HC fermions. Dynamics of the composite vector bosons, referred to as HC rho in this paper, are formulated together with HC pions by the hidden local symmetry (HLS), in a way analogous to QCD vector mesons. Then coupling properties to the standard model (SM) fermions, which respect the HLS gauge symmetry, are described in a way that couplings of the HC rhos to the left-handed SM quarks and leptons are given by a well-defined setup as taking the flavor mixing structures into account. Under the present scenario, we discuss significant bounds on the model from electroweak precision tests, flavor physics, and collider physics. We also try to address B anomalies in processes such as B ->K(*) mu+ mu- and B ->D(*) tau nu, recently reported by LHCb, Belle, (ATLAS, and CMS in part.) Then we find that the present model can account for the anomaly in B ->K(*) mu+ mu- consistently with the other constraints while it predicts no significant deviations in B ->D(*) tau nu from the SM, which can be examined in the future Belle II experiment. The former is archived with the form C9 = -C10 of the Wilson coefficients for effective operators of b ->s mu+ mu-, which has been favored by the recent experimental data. We also investigate current and future experimental limits at the Large Hadron Collider (LHC) and see that possible collider signals come from dijet and ditau, or dimuon resonant searches for the present scenario with TeV mass range. To conclude, the present b ->s mu+ mu- anomaly is likely to imply discovery of new vector bosons in the ditau or dimuon channel in the context of the HC rho model. Our model can be considered as a UV completion of conventional U(1)' model.


Introduction
The CERN Large Hadron Collider (LHC) has discovered [1,2] a Higgs boson, which is the last element to compensate the particle content predicted in the standard model (SM). Indeed, the SM explains particle phenomenologies in good agreement with experiments that include the Higgs coupling measurements. Though the structure of the electroweak symmetry breaking (EWSB) in the SM, i.e., the origin of the negative-mass squared for the Higgs field, is still mysterious, the present LHC data on the EW interactions and Higgs coupling properties are likely to imply that new physics beyond the SM would have no direct correlation with the dynamical structure of the EWSB. This might be why it is hard to see new particles in a TeV scale range. Nevertheless, it would still be worth investigating a bulk of new physics at TeV scale as long as it is reachable at the LHC and insensitive to EW precision variables. One interesting key to access such a new physics paradigm involves a vectorlike scenario, in which new particles are charged vector-likely under the SM gauges so that it is fairly insensitive to the EWSB structure. Among vectorlike models, one of viable candidates would be embedded in a strongly coupled sector, often called hypercolor (HC), or vectorlike-confinement, or chiral-symmetric technicolor [3][4][5][6]. This class of scenarios predicts variety of new phenomena in TeV scale physics such as presence of new composite particles which can be searched at the LHC [5] and affect the Higgs coupling property [4]. Dark matter candidates (HC composite scalars or fermions) are also involved in some models (see e.g., Refs. [6][7][8][9][10][11][12][13][14][15][16]) #1 .
Dynamics of composite vectors is formulated by use of the hidden local symmetry (HLS) [40][41][42][43][44] -the method established in addressing QCD vector mesons including offshell properties and couplings to the external gauge fields -based on an extension of nonlinear realization of "chiral" SU (8) L × SU (8) R symmetry for the HC fermions whose couplings to the SM gauge fields are vector-like. In this case, the SM gauge symmetries SU (3) c × SU (2) W × U (1) Y are realized as unbroken symmetries (diagonal subgroup), come down from a spontaneous breaking of the extended SU (3) c ×SU (2) W ×U (1) Y ×[SU (8)] HLS gauge symmetries. The spontaneously broken gauge-degree of freedom of SU (8) reflects the presence of 63 massive composite gauge bosons (HC ρ's), which in part arise as mixtures of the HLS and the SM gauge bosons. Then the couplings of the SU (8)-63 plet composite vectors to the SM fermions and SM gauge bosons are unambiguously determined by the HLS-gauge invariance. Since some of composite vectors can mix with the SM gauge bosons, the couplings to the SM fermions are limited by electroweak precision tests as fundamental requirements. As we will see, the present scenario is safe for flavor universal EW measurements such as the oblique parameters, while is severely constrained from flavor-dependent processes, in particular, from the forward-backward asymmetry of tau lepton A (0,τ ) FB and Z boson decay to bottom quark pair R b .
We then turn to flavor specific phenomena and investigate constraints from current experiments for flavor and collider physics on the parameter space of the HC ρ couplings to the SM fermions. The flavorful HC ρ couplings are severely constrained so that the third-generation couplings and their small mixings to the second-generation fermions are only viable. In particular, we will see whether the present scenario can accommodate recent anomalies which have been seen in several B meson decays such as B → K ( * ) µ + µ − and B → D ( * ) τν. A vast amount of works have been made in the last several years on the topics of anomalies in the b → sµ + µ − data (see, e.g., Refs. [45][46][47][48][49] for earlier works) #2 . Deviations of experimental data in the ratios R D ( * ) = B(B → D ( * ) τν)/B(B → D ( * ) ν) have also been studied in several specific models (see, e.g., Refs. [56][57][58][59][60][61][62][63][64][65] for earlier works.) We emphasize that the present HC scenario naturally includes various and flavorful new vector bosons that can address B physics, as a consequence of the presence of the strong HC sector at high energy scale. This is a crucial difference from other low-energy effective models (such as a U (1) model) in which new vector bosons are introduced without any dynamical reason. In particular, we will see that the present HC model can be considered as a UV completion of conventional U (1) models.
The implications to the collider physics at LHC are also addressed so that the HC ρ mesons with mass of TeV scale can be consistent with the current 13 TeV bounds from dijet, di-tau and di-muon searches. Finally, we discuss discovery and/or exclusion potential for the HC ρ mesons at future LHC experiments with higher luminosity and then see that the future LHC will be sensitive enough to discover or exclude those TeV mass HC ρ mesons as well as to confirm (in)consistency with the B anomalies. This paper is structured as follows: In Sec. 2 we start with introducing the one-family #2 Explanations by contributions of composite bosonic particles in the context of composite Higgs scenarios were also reported [50][51][52][53][54][55]. Table 1. The SM charge assignment for eight HC fermions F L/R = (Q, L) T L/R in the one-family model. model of the HC and then derive the low-energy effective Lagrangian including the HC pions and HC rho mesons with respect to the HLS formulation. Then we describe the HC rho couplings to the SM fermions and find that there are two types allowed by the HLS gauge invariance -one is of direct coupling type (generically dependent on fermion flavors) while another is of flavor-universal type (arising from the mixing between the HC rho mesons and the SM gauge bosons). After a brief introduction of the HC pion coupling forms as well, we summarize phenomenological features for the present HC rho model in Sec. 3 -mass splittings of the HC rho mesons, oblique corrections of the EW sector, and flavor-dependent corrections to A (0,τ ) FB and R b . These provide fundamental requirements for the present model. In Sec. 4, we discuss the flavor physics issues to give a list of flavor observables relevant for the present scenario of the model. The current flavor bounds are then placed on the HC rho mass and the coupling strength to obtain the allowed parameter space. Finally, in Sec. 5, we show the HC rho sensitivity at the LHC experiments by taking into account resonance searches from di-jet, di-tau and di-muon channels. Sec. 6 is devoted to summary of this work. We also provide a couple of Appendices to show details of computations leading to the results described in the text, and some details related to HC pion physics which are subdominant in the present study.  Table 1. Thus the HC sector possesses the approximate global "chiral" U (8) F L × U (8) F R symmetry (explicitly broken in part by the SM gauging and possible vectorlike fermion masses), among which the U (1) F A part is explicitly broken by the axial anomaly in the same manner as in the QCD case.
At the scale Λ HC , say Λ HC = O(1 − 10) TeV, the HC gauge interaction gets strong to develop the nonzero "chiral" condensate F A F B ∼ Λ 3 HC · δ AB (A and B being indices for SU (8) fundamental representations), which breaks the "chiral" symmetry of 8 HC fermions down to the vectorial one: According to the spontaneous breaking, the 63 Nambu-Goldstone (NG) bosons emerge, which will be pseudoscalars by the explicit breaking terms including the SM gauge interactions and possibly present vectorlike fermion mass terms like m 0 FF F , as discussed in Refs. [3][4][5][6]. By naively scaling the hadron spectroscopy in QCD, we may find 63 composite vectors (HC ρ mesons) as the next-to-lightest HC hadrons #3 #4 . Thus the low-energy effective theory of the HC sector would be constructed from the 63 HC pions (∼F A iγ 5 F B ) and also 63 HC rho mesons (∼F A γ µ F B ). Then the HC rho couplings to the SM particles arise indirectly from mixing between the SM gauge bosons ("indirect couplings"), and directly by an extended HC theory ("direct couplings"), which could be generated from extended (vector, or scalar) interactions communicating the HC and the SM fermion sectors (which would be like a generalized extended technicolor scenario). Note that the latter direct coupling can generically be flavor-dependent. Both types of couplings are unambiguously formulated by the HLS formalism, which is the main target of the subsections below.

HLS formulation
In this section we formulate the effective Lagrangian including the HC vectors along with the HC pions, arising from the one-family model of the HC introduced in the previous section.
The effective Lagrangian for those vectors can be formulated based on the HLS formalism, which has succeeded in QCD rho meson physics [40][41][42][43][44], where the ρ's are introduced as the gauge bosons of the HLS. Based on the nonlinear realization of the HLS and the "chiral" SU (8) F L × SU (8) F R symmetry, the Lagrangian is written as #5 The lightest f0(500) scalar meson in QCD may be a mixture ofqq and qqqq due to some feature specific to the three-flavor (u, d, s) structure, according to recent analyses, e.g., see Ref. [66]. #4 If one happens to take some special value of NHC, e.g. NHC = 3, then the one-family HC dynamics would turn from the ordinary QCD to be quite different, so-called, the walking gauge theory having the almost nonrunning gauge coupling (i.e. approximate scale invariance). Then the lightest hadron spectrum would include a dilaton with the mass as small as the HC pion mass, arising as the (approximate) spontaneous breaking of the scale invariance [67]. In the present study, for simplicity we will disregard the possibility of light dilaton. #5 We have imposed the C and P invariance as well, and assumed that the P invariance is violated only through the weak interactions when the HC pions and ρ mesons couple to the SM fermions through the external weak gauges (See also the ρ couplings to the SM fermions in Eq.(2.23), which take the left-handed form).
with the HLS gauge coupling g ρ , the HC pion decay constant f π , and the external gauge fields R µ and L µ that are associated by gauging the "chiral" symmetry. Ellipses include terms of higher derivative orders [44,68]. Under the HLS and the "chiral" symmetry, the transformation properties for basic variables -ξ L,R (nonlinear bases), ρ µ (HLS field), and α ⊥µ ,α ||µ (covariantized Maurer-Cartan one forms) -are described as The HLS is thus the gauge degree of freedom independent of the SM-external gauges and has spontaneously been broken together with the "chiral" symmetry in terms of the nonlinear realization: ξ L,R = 1. The nonlinear bases ξ L and ξ R can be parametrized by the NG bosons π for the "chiral" symmetry and P for the HLS. Hence, they are parametrized as ξ R L = e iP/f P · e ±iπ/fπ , where the HLS decay constant f P is related to the HC rho mass as m ρ = g ρ f P and then the Ps are eaten by the HLS gauge boson ρ µ due to the Higgs mechanism.
Note that, by construction, at the leading order of derivative expansion the HLS completely forbids triple vector vertices involving the external gauge fields along with HLS, such as G − G − ρ and G − ρ − ρ (G = L, R). When going beyond the leading order of derivative expansion, e.g., at O(p 4 ) one will encounter mixing terms such as tr[Ĝ µν ρ µν ] wherê G µν = ξ L/R G µν ξ † L/R [44,68]. The size of effects is on the order of one-loop, O(1/(4π) 2 ), and thus should potentially be small as long as the derivative expansion works as in the case of the chiral perturbation for pions.
Hereafter, we shall focus on the 63 composite ρ vectors, namely, the HC rho mesons that couple to the SM fermions. Possible interference effects from the HC pions on the vector meson physics will be addressed later on.

Particle assignment for HC ρ and HC π
The HC ρ fields are constructed from the underlying HC fermions as listed in Table 2 and expanded by the SU (8) group generators (with the Lorentz vector label µ suppressed) #6 : The way of embedding the ρ fields into the one-family 8 × 8 matrix form has been chosen to be on the basis of allowing to mix with the SM gauge bosons, so-called the Drell-Yan basis in light of the LHC production.
composite vector constituent color isospin Table 2. The HC rho mesons and their associated constituent HC quarks Q c = (U, D) c and leptons L = (N, E). Here λ a (a = 1, · · · , 8) are the Gell-Mann matrices, τ α SU (2) generators defined as τ α = σ α /2 (α = 1, 2, 3) with the Pauli matrices σ α , and the label c stands for the QCD-three colors, c = r, g, b. The numbers attached in lower scripts (1,3,8) correspond to the representations under the QCD color, i.e., singlet, triplet and octet for (1,3,8). Since the unbroken color symmetry acts as an internal degree of freedom, four color-octet (real), four color-triplet (complex), and seven color-singlet (real) states are observed as distinguishable particles. with T [2] (3)c = where τ α = σ α /2 (σ α : Pauli matrices), λ a and e c represent the Gell-Mann matrices and three-dimensional unit vectors in color space, respectively, and the generator T A is normalized as tr[T A T B ] = δ AB /2. For color-triplet components (leptoquarks), we define the following eigenforms which discriminate 3 and3 states of the SU (3) c gauge group, (3)c + iT [2] (3)c = 1 2 Also, the following relations hold: which are useful for rewriting the above decomposition in the eigenforms. The following normalization conditions of the 3 and3 states of the SU (3) c gauge group in the new bases are also useful: c ] = δ cc /2, (others) = 0. We can express the vector fields ρ of 8 × 8 matrix with block matrices as Thus, distinct particles consist of four color-octet (real), four color-triplet (complex), and seven color-singlet (real) states. The SU (2) W charge, in terms of which the HC rho field is decomposed as in Eq.(2.8), is identified with the one in the SM. The SM fermions will carry corresponding SU (2) W charges so that they couple to the HC rho once a relationship for the external SM gauge fields is given, as will be seen later. In a manner similar to the HC rho mesons in Eq.(2.9), the HC pions are embedded into the 8 × 8 matrix of SU (8) as , 2.4 Short sketch for masses of HC ρ and HC π As evident from Eq.(2.1), the HC ρ meson masses are degenerate due to the global SU (8) F V symmetry unless the external gauge fields are present. With the SM gauges tuned on, however, some of the HC ρ mesons charged under the SM gauge fields will get shifts to masses, due to the mixing with the SM gauge bosons, (see Eq.(2.26) in the later discussion.) Note that the massless property of the gluon and the photon is intact under the mass mixings because they are protected by the residual SU (3) c × U (1) em gauge symmetries after the Higgs mechanism works. Most sizable corrections will arise from mixing with the QCD gluon, which makes the mass of the color-octet isosinglet ρ 0 (8) lifted up by amount of O(4πα s /g ρ ) = O(0.1) for g ρ = O(10) and α s = g 2 s /(4π) 0.1. The HC pions listed in Eq.(2.10) get massive through the explicit breaking effect outside of the HC dynamics. One may make them massive by introducing some Dirac masses (m 0 F ) for HC fermions (like m 0 FF F ). To realize the HC pions as NG bosons for the spontaneous breaking of "chiral" symmetry, such an explicit HC fermion mass should be small so that M π O(Λ HC /(4π)) O(f π ). This implies M π ∼ O(100) GeV for Λ HC ∼ O(1) TeV. However, the present model can have another source which would allow some HC pions to get more massive: since the present HC theory consists of the one-family content with the number of HC fermions N F = 8, the masses of HC pions having the SM charges could be enhanced by the amplification of the explicit breaking effect, as discussed in [69] and references therein. The enhancement will then be most eminent for QCD colored pions, π (3) and π (8) due to the relatively large QCD coupling strength. Following [69], we evaluate the size of colored HC pion masses from the QCD gluon exchange contribution as , with C 2 = 4 3 (3) for color-triplet (octet) HC pions, where Λ UV denotes some ultraviolet high-energy scale up to which the HC theory is valid. Taking α s (M π ) ∼ 0.1 and Λ UV ∼ 10 16 GeV, for example, we thus estimate the π (3) and π (8) masses as M π (3) ∼ 3 TeV and M π (8) ∼ 4 TeV, respectively, for Λ HC ∼ 1 TeV.
In a similar way, the EW gauge interaction makes masses of EW-charged HC pions lifted up. This effect becomes operative for the π ±, 3 (1) and π ±,3 (1) pions to yield M π ±,3 (1),(1) ∼ 2 TeV for Λ HC ∼ 1 TeV and Λ UV ∼ 10 16 GeV as a benchmark. Hereafter, The indices '±' and '3' discriminate components of SU (2) W triplets. The index '0' emphasizes that the designated states are SU (2) W singlets. Thus, the sizes of the HC pion masses are roughly expected as for Λ HC ∼ 1 TeV and Λ UV ∼ 10 16 GeV. This is the significant feature for the HC pion in our model particularly when we discuss collider bound on the HC rho mesons. Hereafter we shall take the above HC pion spectroscopy as a benchmark in the present study on the HC rho meson physics. Most of HC pions thus become as heavy as HC rhos, which implies the loss of pseudo NG nature in the full quantum field theory, though at the classical level the "chiral" symmetry is conserved approximately enough due to the smallness of the gauge couplings. As was discussed in Ref. [69], this happens due to the amplification of explicit breaking effect, induced from non-perturbative feature of the QCD with many flavors (nearly conformal/scale-invariant dynamics). One might then think of integrating out heavy HC pions to get a reduced low-energy description characterized by smaller "chiral" manifold than the original class of SU (8). However, one cannot do it because heavy HC pions in the mass eigenstate basis involve all HC quark and leptons when they are formed, as seen from Eq.(2.10) (or in the same way of embedding as in Table 2 for HC rhos.). Thus one needs to construct the "chiral" SU (8) manifold to describe HC pions, which are manifestly NG bosons at the classical level of the "chiral" SU (8) dynamics, even though it becomes so heavy at the quantum level.

Couplings to SM particles
The SM fermion fields are written as an eight-dimensional vector on the base of the fundamental representation of SU (8), where q and l are SU (2) F L ,F R doublets for the quark and lepton fields. The SM-covariant derivatives that act on the f -fermion multiplets are then expressed as the 8 × 8 matrix forms: where G µ , W µ and B µ are the SU (3) c × SU (2) W × U (1) Y gauge fields along with the gauge couplings g s , g W and g Y , respectively; and Q q,l em is the electromagnetic (EM) charge defined as As clearly seen from Eqs.(2.13) and (2.14), the global SU (8) F L × SU (8) F R symmetry is explicitly broken in the SM fermion sector by the SM gauges, hence of course the SU (8) is not a good symmetry for the SM fermions. Just for convenience to address SM fermion couplings to the HC rhos which form the SU (8) adjoint representation, we have introduced the SU (8) F L /F R -multiplet form for the SM fermions as in Eq.(2.12). The covariant derivatives for the HC fermions can also be written in terms of the 8 × 8 matrix form. We may relate the charges of the HC fermions with those of the SM quark and lepton charges, involving the HC-quark and -lepton numbers. Then the nonlinear bases ξ L,R in Eq.
It is useful to expandα ||µ andα ⊥µ in Eq.(2.2) in powers of the HC pion fields π with the unitary gauge for the HLS (P ≡ 0): with V µ = (R µ + L µ )/2 and A µ = (R µ − L µ )/2. Then the 1-forms in Eqs. (2.18) and (2.19) are represented asα (2.20) We may define the dressed fields for the left-handed SM fermions, which transform as These transformations allow us to write down the HC ρ couplings to the left-handed SM fermions in the HLS-invariant way as where i and j label the generations of the SM fermions (i, j = 1, 2, 3). Using Eqs. (2.13) and (2.20), one can thus extract the HC ρ and V SM (SM gauge boson) couplings to the left-handed SM fermions. As a result, we have where ρ µ QQ , ρ µ LL , and ρ µ QL are combinations of the HC ρ mesons as defined in Eq.(2.9) and g ij L = (g 1L + 2g 2L + g 3L ) ij . Note that the V SM -f L -f L term in Eq.(2.24) is not the normal SM interactions but additional contributions in this model.
The HLS invariance actually allows one to write down vector couplings other than those in Eq.(2.24), which would take the form like with the generation-dependent coupling h ij L . As seen from Eq.(2.20), however, the 1-form α µ⊥ goes to vanish in the unitary gauge of the HLS; ξ L/R → 1 up to HC pion terms ∂ µ π/f π + · · · . Since HC rhos and pions are generically composed of identical HC fermion bilinears, they are expected to develop the couplings to the SM fermions in the same form, which could be generated from an extended HC theory. Explicit modeling of the extended HC sector is beyond the scope of current analysis. We will briefly address possible effects from those HC pion couplings in the later section.

ρ -V SM mixing structures and induced-indirect couplings to SM fermions
In addition to the direct interactions of Eq.(2.24), the HC ρ mesons also have interactions induced by the mixing with the SM gauge bosons. The mixing term is involved in the mass matrix of the vector boson, which is written by where we used the relations in Eqs.(2.6), (2.7), (2.13) and the normalization of the SU (8) generators as tr T A T B = δ AB /2. Note that the mixing form is manifestly custodialsymmetric. At first, we illustrate the case for charged bosons, The corresponding matrix element is then written by where v VEV is the Higgs VEV which is the same one as in the SM. Expanding the above mass matrix with respect to r g , one can obtain the mass eigenvalues as with r m = m W /m ρ and then the mass eigenstates W , ρ (1) , and ρ (1) as ('±' is omitted above), to the nontrivial leading order in expansion for r g 1 and r m 1. Next, the neutral vector fields have the following mass matrix: (2.33) Here, our notation on the electroweak sector is as follows: One immediately finds that the above mass matrix has zero determinant, which ensures the existence of the massless photon. The mass eigenvalues are then obtained and expanded as follows: The Z mass (M Z ) multiplied by cos θ W in the above appears to differ from the W mass (M W ) in Eq.(2.28) at the order of O(r 2 g ). The deviation can actually be absorbed by redefinition of the Weinberg angle, so that no sizable correction to the ρ parameter is present, consistently with the custodial symmetric mixing form in Eq. (2.26). (This will also be clearly seen when the ρ parameter is explicitly evaluated below.) The mass eigenstate A and Z boson fields arise from the gauge eigenstate fields to the nontrivial leading order of r g and r m . One can easily see that, to the leading order of the (r g , r m )-perturbation, the weak mixing structure for the neutral A and Z gauge bosons are precisely the same as in the SM. As will be shown later, the vector mixing in the present model does not give significant corrections to SM interactions of the EW sector so that oblique corrections are negligible up to the order of r 2 g . Finally, we examine the colored part written by 1 2 which holds zero determinant and thus guarantees the massless property of the gluon. The mass eigenstatesG a andρ (8)a are given as Here, the ratio r gs is defined as g s /g ρ . The indirect couplings of the HC ρ mesons to SM fermions thus arise from the above flavor-universal V SM -ρ mixings in the mass eigenstates. As seen from the expressions of the mixings, such flavor-universal couplings are suppressed for r g 1, namely, for large g ρ , which is also required for the oblique corrections to be negligible. This is, indeed, inferred from the QCD case. (See the later sections.) On the other hand, flavor-specific couplings of the HC ρ mesons to SM fermions are also given with the form g ρ g ij L as in Eq.(2.24) and then it can significantly contribute to variety of flavor processes, as we will see in the next section.

Couplings including HC π
From the chiral Lagrangian in Eq.(2.1) with the concrete form of the covariantized Maurer-Cartan one forms in Eqs. (2.18) and (2.19), we find that the following types of HC pion coupling terms emerge after the expansion (up to the quartic order in fields): ρ-π-π, V-π-π, V-V-π-π and π-π-π-π. Their interaction forms easily read The specific choice, a = 2, turns out to make the V-π-π term vanishes at the leading order. This is referred to as the vector dominance in which the chiral perturbation theory reproduces experimental results regarding QCD. For the present study, we may therefore assume the vector dominance scenario and then see that the ρ-π-π coupling g ρππ is completely set by the g ρ : g ρππ = g ρ . As a reference number, we take #7 g ρ | vector dominance = g ρππ | QCD 6 . (2.44) Note that this reference size of g ρ is consistent with the perturbative expansion in the charged and neutral gauge boson sectors as shown above. (The full set of concrete forms of the ρ-π-π couplings are provided in appendix A.) We also note that taking into account of chiral anomalies results in the presence of (covariantized) Wess-Zumino-Witten (WZW) terms [72,73], which contain three-point interactions of ρ-ρ-π, ρ-V-π and V-V-π. The HC pion couplings to the SM fermions would be provided by an extended-HC theory, as written down in Eq.(2.25). Phenomenology regarding those couplings will briefly be mentioned later when the collider signatures are discussed. #7 The dependence of the ρ-π-π coupling gρππ on the number of constituent fermions (NF ) would not be significant, although the overall coefficient of ρ-π-π vertex scales with 1/ NF /2 as discussed in Refs. [70,71].

Phenomenological features of HC rho mesons
As explained in the previous section, the HC composite ρ mesons arise from the SU (N HC ) HC gauge theory with 8 new HC fermions at the strong scale of O(1 − 10) TeV. The total 63 (or 15 distinct) massive composite vectors, ρ µ , listed in Table 2 have identical masses m ρ in the gauge basis. Those 63 HC ρ meson properties are characterized as follows: • the HC ρ mesons are categorized into three types for SU (3) c : octet (ρ α,0 (8) ), triplet (ρ α,0 (3) ), and singlet (ρ α,0 • ρ (3) couples to quark-lepton pairs, referred to as leptoquark [74], • no direct interaction to gluon even though ρ (8) and The new interactions of V µ (= ρ µ , V µ SM ) to the SM fermions are described in (2.24). There also exists the ρ µ -V µ SM mixing in the mass matrix. In the following part, we will discuss phenomenological features of the present HC ρ model, derived from the above fundamental interaction properties.

Mass splitting of HC ρ mesons
As seen in the previous section, one of the charged HC ρ meson, ρ ± (1) , gets the correction on the mass due to the ρ -V SM mixing, while another one ρ ± (1) does not change its mass. Applying the same analysis, the masses of all the HC ρ mesons in the mass eigen-basis are obtained as where r g = g W /g ρ and r gs = g s /g ρ . Recall that we now take g ρ = 6 [in Eq.(2.44)] as a reference number by following the discussion in the previous section. Thus, one finds that the masses of the HC ρ mesons are almost degenerated even in the mass eigenbasis.

Flavor universal corrections to the EW sector
The mass mixing forms in Eq.(2.27) and (2.33) imply the existence of non standard EW couplings to SM fermions. In the case of the large g ρ , as in the reference point inferred from the QCD case (See Eq.(2.44)), the EW gauge fields (W ±,3 , B) can be treated as external gauge fields. In that case, one can directly read off the current correlators (Π: transversely polarized gauge boson amplitudes) composed of SM fermions from the mass matrices in Eqs. (2.27) and (2.33). Taking into account new physics contributions from the ρ ±,3 (1) and ρ 0 (1) exchanges, one thus finds where Q 2 ≡ −q 2 is the space-like external momentum squared, and terms of O(r 2 g ) corrections have only been included. With these correlators at hand, one can extract the flavor-universal EW (oblique) parameters, calledŜ,T , W and Y [75]: which are converted into the Peskin-Takeuchi S and T parameters [76] as [77] where α em = e 2 /(4π). From Eq.(3.2) we obtain the oblique parameters aŝ where r m = m W /m ρ and r g = g W /g ρ . Therefore, the S and T parameters in the present HC ρ model are suppressed by a factor r 2 m · r 2 g . Up to the order of O(r 2 The result implies that the HC ρ mesons are insensitive to the flavor-universal EW sector.

Key contributions of
Besides the (negligible) flavor-universal interactions of the HC ρ, there are lots of flavorspecific interactions in the ρ -f L -f L terms as described in Eq.(2.24). Here we summarize and extract significant terms relevant for flavor and collider physics.
Decomposing the Lagrangian of Eq.(2.24), we have for neutral HC ρ and for charged HC ρ. We also have terms from the indirect interaction of Eq.(2.26), but contributions are suppressed due to large g ρ (= 6 as in Eq.(2.44)). Combining the two contributions, it is typically given such as For almost all case, we simply neglect the suppressed terms of ρ (1,8) -f -f ( ) in the following study and will get back to this point for As for leptoquark-type HC ρ, we have 14) Here, we suppressed the indices of the SU (3) c fundamental representations. The leptoquark-type HC ρ has no indirect term from the beginning.
Finally, ρ (8) couples to quarks with the following form.
Note that all the above ρ − f − f couplings accompany with g ij L , the size of which will actually be constrained by flavor-dependent EW precision tests as seen below.

Deviations of flavor-dependent
The flavor-dependent couplings with the strength of g ij L in Eq.(2.24) give rise to corrections to the SM gauge-fermion-fermion coupling properties. The effective vertex for such a V SMf -f interaction consists of two sources: the direct term of V SM -f -f and the term of ρ-f -f through the ρ-V SM mixing. It is typically described by at the low energy scale (q 2 = −Q 2 m 2 V SM ). This is the salient feature of the HLS formulation: the composite vector bosons as gauge bosons in HLS are correlated with the SM gauge bosons V SM (= G, W, B), through the nonlinear realization.
As a consequence, the corrections in the EW sector are obtained as with the usual notations of I ± = τ 1 ± iτ 2 , and I 3 = τ 3 , where the effective corrections ∆ W/Z are written by .
Note again that m V is the mass in the gauge basis whereas M V in the mass-eigen basis.
For the final form of the above result, we keep the term up to O(r 2 x ) for x = m, g. A key point of this model is that the coupling g ij L in the shift parameter of Eq.(3.22) share with the HC ρ couplings to the SM fermions [in Eq.(2.24)] as a consequence of the HLS formulation, which allows to introduce the SM gauges as a remnant of the spontaneous breaking of the (gauged)"chiral" and hidden local symmetries. This crucial feature puts severe constraints on both diagonal/off-diagonal components of g ij L in the gauge basis, as we will see later.

Flavor-dependent constraints from the EW sector
In this subsection we show constraints on the flavorful coupling strength ∆ ij W/Z (g ij L ) in Eq.(3.21), required from flavor-dependent EW precision tests. We also discuss a reasonable setup for the parameters, followed by flavor and collider limits in the next sections.
As can be seen in Eq.(2.24), our model involves lots of new interactions at the tree level, most of which are obviously already disfavored. In particular, it is easily expected that couplings to the first and second generations are severely constrained. To avoid such matters as well as to address the flavor anomalies in B decays, the reasonable setup may be given as #8 (3.23) Even when we put the above assumption, constraints on g 33 L from precision measurements of the electroweak sector must be concerned, which is rephrased in a way that the Z-boson couplings to the SM fermions were measured very precisely and a sizable deviation from #8 The flavor-blind choice such as g ij L = diag[1, 1, 1], where flavor-changing interactions emerge in the leptoquark part in mass eigenbasis, is disfavored mainly due to severe constraints from the LHC direct searches of new resonances. Note again that the third-generation-philic form is shared with both HC rhos and pions. Its realization by an extended HC sector could be interesting as a future study. the SM is immediately disfavored. In the present model, the form of the Z-f -f couplings is given as  [78]. In the present model, left-handed and righthanded couplings of the Z boson are read off from Eq.(3.25) as where we ignored the mixing effect of the lepton fields by rotating gauge to mass eigenbases. The present experimental value A (0,τ ) FB | exp = 0.0188 ± 0.0017 and the SM prediction A (0,l) FB | SM = 0.01622 ± 0.00009 [79] still allow the small discrepancy, within the uncertainty. This is then compared with the HC-rho model value where we take s 2 W = 0.2336 in the numerical analysis #9 . As a result, the allowed range is obtained as The shifts for the Z-boson onshell couplings also contribute to the Z → hadrons decay amplitudes. The most precise measurement has been done by R b = Γ(Z → bb)/Γ(Z → hadrons). #9 The Weinberg angle (so-called the Z mass shell value) is derived from the relation with the input variables GF = 1.1663787 × 10 −5 GeV −2 , MZ = 91.1876 GeV, αem(MZ )| −1 MS = 127.950 [79].

It is described as R b
(|g bL | 2 + |g bR | 2 )/ q=u,d,s,c,b |g qL | 2 + |g qR | 2 neglecting the tiny phase-space difference, where the Z boson couplings to quarks (g q L ,q R ) in the present model are read off from Eqs. (3.25): where we again ignored the mixing effect of the quark fields via the flavor rotation in the mass-eigen basis, (which is indeed the case for our scenario as seen in the next section.) The present experimental value of R b | exp = 0.21629 ± 0.00066 and the SM prediction R b | SM = 0.21579 ± 0.00003 [79] allow the tiny discrepancy, Comparing this with the HC-rho model value Combing the above constraints, the allowed range for g 33 L is obtained as (3.36) at 95% C.L.. As we will discuss later, the 1σ range points to a negative value for g 33 L which excludes a low-mass range for the HC rho mass scale to be consistent with flavor limits.
To address flavor-changing processes, we may introduce the mixing structure between the second and third generations of the down-type quarks and leptons as will be seen in Eq.(4.2). From the flavor issues discussed in the next section, the down-quark mixing angle is required to be as small as 10 −3 − 10 −2 and thus R b is the most relevant EW limit. As for the lepton mixing angle, on the other hand, we will consider the two cases θ L π/4 and θ L ∼ π/2. In the former case, A  We will get back to this result later.

Contribution to flavor observables
The HC ρ mesons provide rich phenomenologies addressing flavor physics as will be discussed in this section. One easily finds that the interaction terms of ρ-f -f in Eq.(2.24) and Eqs.(3.7)-(3.19) produce flavor changing neutral currents (FCNC) that are severely constrained.
With the assumption g 33 L = 0 (and others = 0) as in Eq.(3.23), there is no FCNC term in the gauge basis. This setup, however, still causes FCNCs in transforming from the gauge basis to the mass basis: where U , D, and L are three-by-three unitary matrices and the spinors with the prime symbol denote the fermions in the mass basis #10 . The capital latin indices I, J identify the mass eigenstates. The CKM matrix element is then given by According to the literature [80], in order to address several flavor anomalies recently reported in measurements ofB → Kµ + µ − (andD ( * ) τν) as well as to avoid severe constraints of FCNCs in the first and second generations, the mixing structures of D and L are reasonably parametrized by Through these flavor mixings, we will see significant contributions to flavor phenomenologies. The following factors are useful, 3)

7)
#10 Ambiguities can remain in the transformation of the right-handed neutrinos when (active) neutrinos are massive. Here, we consider massless neutrinos, where no ambiguity remains.
with their charge conjugations On the other hand, the indirect couplings derived from Eq.(2.26) have no contribution to FCNC phenomena at a tree-level #11 since these flavor transformations do not produce FCNCs. Note that charged flavor transitions such as b → c can be affected and we will see this point soon later.

Effective four-fermion operators from the HC ρ and π meson contributions
We see from Eq.(2.24) that the octet HC mesons ρ α,0 along with the Wilson coefficients having the form with an appropriate factor, derived from Eqs.(2.9) and (2.24), for every HC ρ meson. When we recall the degenerated ρ meson masses, we have where ∆ ij;kl = 1/2 for two identical currents (i = k, j = l) and 1 for the others. Note that the O ijkl 2q2 (3) term is canceled out in this model with the degenerated masses. This term can be nonzero only if . The HLS-gauge invariance in the present model, however, does not allow such a sizable mass-split between them. On the other hand, scalar-type contributions to four fermion processes are given by HC pion exchanges in principle. From the interaction of Eq.(2.25), corresponding contributions involve fermion masses (derived from derivative couplings of HC pions to external fermions), which is the NG boson nature of HC pions. Therefore, a typical size of the contribution is written as where h ij L are the coefficients of the operators shown in Eq.(2.25), m π represents a typical scale of the HC pion mass, and m f, max indicates the largest mass of the external fermions which describes a typical scale of the momentum in meson-associated flavor physics.
This is compared with the one from the HC rho, g 2 ρ g ij L g kl L /m 2 ρ . The couplings could be comparable such as h ij L ∼ g ρ g ij L whereas m π becomes comparable with m ρ as discussed in Sec. 2.5 and Sec. 2.4, respectively. Thus, we expect that the scalar contribution is suppressed by the factor m 2 f,max /f 2 π , which is maximally of O(m b /O(10 2 GeV)) 2 (1/20) 2 . Thereby, it is totally natural to neglect HC pion contributions to flavor physics and then we focus only on the HC rho contributions. #12

Effect on flavor observables
When we transform the left-handed SM fermion fields to those in the mass basis by U , D, and L under the assumptions of Eqs.(3.23) and (4.2), the Lagrangian of Eq. (4.14) generates nonzero FCNC contributions in the b -s system and the lepton flavor violating system. To be specific, intriguing processes are summarized as follows: • O 2q2 (n) :B → D ( * ) τν (only for n = 3),B → K ( * ) µ + µ − ,B → K ( * ) νν, and τ → φµ, Then Eq.(4.14) implies that the net effect of the HC ρ mesons (which arises from the direct ρ-f -f term of Eq.(2.24)) onB → D ( * ) τν is exactly zero at the tree level, even though each color-triplet and color-singlet HC ρ meson does contribute.
In the following main body of the manuscript, we skip to show the prime symbol for fermion fields to emphasize them being mass eigenstates.

4.2.1B → K ( * ) + −
The effective Hamiltonian for b → s + − in the present model is described by

16)
#12 An exceptional case is top-quark-associated processes, as we see around Eq.(5.5). with where the conventional coefficients C SM 9 0.95 and C SM 10 −1.00 indicate the loop functions from the SM contributions, at the m b scale, (see, e.g., [81]).
Experimental measurements regarding b → s + − have been developed in recent years to search for new physics and then the net results have suggested significant deviations from the SM predictions: • Since the decayB → K * µ + µ − includes detectable final state particles (µ + , µ − , K * → Kπ), it has variety of angular observables. LHCb and Belle have then measured the optimized observable, so-called P 5 [82], and found a sizable discrepancy between their results and the SM prediction, as in Refs. [83][84][85][86]. Recently, ATLAS and CMS have also reported (preliminary) results for several observables including P 5 in Refs. [87,88], which are in agreement with the LHCb and Belle results.
• The tests of lepton-flavor-universality violation (LFUV) have been examined by measuring the ratio LHCb. The result, again, suggests a deviation from the SM prediction in R K [91] and very recently in R K * [92]. Besides R K ( * ) , Belle has recently shown their result on a new measurement of the LFUV effect, P µ 5 − P e 5 [93], which points to a presence of LFUV although the result is not yet statistically significant [90].
Immediately after the very recent report on R K * [92], global fit analyses for new physics have been investigated by several theory groups, see Refs. [55,[94][95][96][97][98][99][100][101][102]. The results are usually shown in terms of the conventional coefficients C ( ) X for X = 9, 10, (where the primed coefficients are those for the operators replacing P L → P R in Eq.(4.16).) According to the literatures, favored solutions to accommodate the present data are obtained in the following three cases (among single degree of freedom of new physics contributions); C µµ 9 = 0, C µµ 9 = −C µµ 10 = 0, and C µµ 9 = −C µµ 9 = 0 #13 . The present model corresponds to the second case C µµ 9 = −C µµ 10 and then the favored region is given as at the 2σ level, whereas the best fit point is −0.61. We have quoted the favored region given in Ref. [94], where all available associated data from LHCb, Belle, ATLAS and CMS were combined.
#13 To be precise, fit to the LFUV data (RK and RK * ) prefers C µµ 9 = −C µµ 10 = 0 while that to all the data favors C µµ 9 = 0.

4.2.2B → K ( * ) νν
The effective Hamiltonian forB → K ( * ) νν is written by where the SM loop function is C SM L −6.36. The experimental upper bounds on the branching ratios ofB → K ( * ) νν [103,104] put constraints on new physics contributions [105]. The most severe constraint is then obtained as [80] −13 at 90% C.L., where

τ → φµ
Our model produces lepton flavor violating decays such as τ → φµ and τ → 3µ at the tree level while they are much suppressed in the SM. The effective Hamiltonian for τ → φµ is written by We take the constraint obtained in [80] from the 90% C.L. upper limit of B(τ → µφ) < 8.4 × 10 −8 [106], which leads (4.23)

τ → 3µ
The lepton flavor violating decay τ → 3µ plays an important role in this model. The effective Hamiltonian is and then the branching ratio is obtained by where the factor 0.94 came from the phase space suppression for the decay [80]. The following experimental upper bound at 90% C.L. is available [107]: This is compared with the experimental measurement [108] ∆M exp. Bs = (17.757 ± 0.021) ps −1 . (4.29) Note that a theoretical uncertainty comes from the input parameters of V tb V * ts and f 2 BsB Bs , which are much more dominant than the above experimental uncertainty e.g., ∆M SM Bs = (17.4 ± 2.6) ps −1 . For a conservative choice, we take ±1σ range for the theoretical uncertainty.

4.2.6B → D ( * ) τν
The semi-tauonic B meson decays ofB → D ( * ) τν were measured [109][110][111] and then it has turned out that the experimental data deviate from the SM predictions. To be specific, with respect to the ratios (for = e or µ), we see the deviations of the combined experimental data [109][110][111] from the SM of 1.7σ and 3.3σ, respectively. On the other hand, the recent update of the analysis for R D * at Belle reported in Ref. [112] shows consistency with the SM. Combining this new result with the others #14 , one finds where the deviations are now 1.7σ and 3.0σ, respectively. Although they are still significant, the SM predictions are consistent within 3σ range. The corresponding effective Hamiltonian that contributes toB → D ( * ) ν can be written as The recent Belle analysis in Ref. [112] was performed using different decay modes of the tag-side B mesons from the previous one in Ref. [110] so that the results in Refs. [110,112] are independent with each other. Hence the two results can be statistically combined.
As we pointed out in Sec. 4.1, the net effect of the charged HC ρ mesons [from Eq.(2.24)] on d → u ν is exactly zero, namely, for the degenerated HC ρ masses. This is actually the consequence of the global SU (8) invariance for HC rhos: the cancellation turns out to take place between ρ α (1) -ρ α (1) and ρ α

Allowed region in parameter space
Now we investigate the parameter space of the model consistent with the flavor measurements. Theoretical input for our evaluation is summarized in Table 3. In Fig. 1 Table 3. One easily sees that the constraint from the B s mixing (magenta region) is much stringent |θ D | 1, while the other constraints are consistent with the b → sµ + µ − anomaly (blue region) in some limited regions.
To more precisely see the allowed region near θ D ∼ 0, in Fig. 2 we show the close-up version focused on the θ D ∼ 0 region, for various values of g ρ g 33 L with m ρ = 1 TeV fixed, where we have taken the significant constraints, namely from ∆M s , B(τ → 3µ), and the b → sµ + µ − global fit. In the close-up plot, we have taken the parameter range favored up to ±3σ level for the b → sµ + µ − anomaly. For |g ρ g 33 L | 1, it turns out that the b → sµ + µ − anomaly is not consistent with the constraints from ∆M s and B(τ → 3µ) in the present model. As for the range 0.35 |g ρ g 33 L | 1, several comments are in order: • We found that there are two isolated regions, θ L π/4 ("left-side") and θ L ∼ π/2 ("right-side"), where all the constraints are (marginally) satisfied.   • The left-side spot is barely viable when the b → sµ + µ − anomaly is 3σ above in terms of the coefficient C µµ 9 from the best-fit point (C µµ 9 | best = −0.61 [94]), where C µµ 9 | +3σ = −0.23, which means that the deviation from the SM has to be rather small.
• On the other hand, the right-side spot fairly satisfies all the constraints. In particular, the point of θ L = π/2 can accommodate the best fit point for the b → sµ + µ − anomaly. Note that θ L = π/2 is the point in which τ in the gauge basis is exactly equivalent to µ in the mass basis. In this case, the lepton sector only has the connection to the HC ρ from the µ-µ-ρ term.
For |g ρ g 33 L | 0.35 we also found that the two allowed spots (which were divided by the τ → 3µ constraint) are merged into a single spot and then the region near θ L ∼ π/2 is only allowed. In Fig. 3, we survey the allowed range of |g ρ g 33 L | for the case θ L = π/2. The result implies that the present model for θ L = π/2 requires |g ρ g 33 L | 0.1 in order to explain the b → sµ + µ − anomaly consistently with the bound from ∆M s .
To summarize, we investigated the allowed regions in the parameter space of θ L , θ D , and |g ρ g 33 L |, which satisfy all the flavor constraints. The situation is then divided by two cases; θ L π/4 ("left-side") and θ L ∼ π/2 ("right-side"). As a result, the allowed region  L , m ρ ) plane fixing g ρ = 6 for θ L π/4 (left) and θ L ∼ π/2 (right). The (darker) gray region is allowed by the EW precision measurements at 95% (68%) C.L. The regions in red (blue) line boundaries are favored by the b → sµ + µ − data that satisfies all the other flavor constraints, for θ L π/4 (θ L ∼ π/2).  The left-side spot is considered as the τ -dominant case where the τ -τ -ρ coupling is relatively larger than the other lepton couplings to HC ρ (including LFV.) On the other hand, the right-side spot only involves the µ-µ-ρ coupling. Indeed, we have to pay attention to this difference when we consider collider limit at LHC as will be discussed in the next section. The combined plots with the constraint from the EW precision measurements are shown in Fig. 4 on the (g ρ g 33 L , m ρ ) plane for the θ L π/4 and θ L ∼ π/2 cases. The regions in red (blue) line boundaries are favored by the b → sµ + µ − data, along with the other constraints, for the left-side (θ L π/4) and right-side (θ L ∼ π/2) spots. The shaded regions show 68% and 95% C.L. constraints from the EW measurements as obtained in Sec. 3.5. (Note that, for the θ L ∼ π/2 case, the EW limit is obtained by combining the R b and A (0,µ) FB constraints.) In the figure, the reference number g ρ = 6 is taken. We can see that the favored regions from the b → sµ + µ − data are consistent with the 95% C.L. EW precision measurements. One also finds that, in the 1σ range (68% C.L.), the EW precision measurements exclude the HC rho mass m ρ 1.8 TeV and the positive value of g ρ g 33 L for the left-side spot case with θ L π/4. The allowed parameter space can be examined by direct searches at the LHC. It will be discussed below.

Collider-related issues
In this section, we discuss constraints from the latest null results in the new physics searches and future prospects at the 13 TeV LHC. The HC ρ's as well as the HC π's will be resonantly, or non-resonantly produced at the hadron collision machinery, to be constrained by the present experimental data.

Typical constraints on HC π
Even though the details of the HC pion sector is out of our major interests, we briefly comment on possible constraints on this part. Here, we focus on the color-singlet isospinsinglet HC pion π 0 (1) as a typical signature. Two types of interactions can be derived for the π 0 (1) . The first one is from the global chiral anomalies of hypercolor fermions, which are represented by the covariantized WZW terms in the present non-Abelian [116,117] (based on the discussion on Ref. [118]), where we focus on V-V-π interactions. M 4 means the four-dimensional Minkowski manifold, U denotes the chiral nonlinear basis U ≡ e 2iπ/fπ = ξ † L ξ R , and the differential one-forms are defined as One-form of the external vector gauge boson includes the gluon as Then we encounter the trace, tr ( where we used the relation dG a dG a = d 4 x G a µν G a ρσ ε µνρσ /4. This effective action describes the anomaly-induced G-G-π 0 (1) interaction. (Note that G indicates gluon in this article.) Another possible origin of the G-G-π 0 (1) coupling comes from a top-quark loop contribution constructed from the t-t-π 0 (1) interaction, which could be, in the present model, induced through an extended HC as given in Eq.(2.25). As discussed in Ref. [119] the amplitude corresponding to such a top-loop contribution is generically written down with the factor g A = α s /(2πv VEV ), the t-t-π 0 (1) coupling strength r y being real and defined as a ratio to the SM like top-Higgs case (m t γ 5 /v VEV ), and the gluon polarization vectors µ,ν . Here, we adopt the Feynman rule of the fundamental pseudoscalar to parametrize the coupling. The loop function f (τ ) is a function of the parameter τ ≡ 4m 2 t /M 2 with η ± = 1 ± √ 1 − τ . Combined with the above two sources, the squared amplitudes are computed as with the factor Through the well known formula for cross section with a spin-J resonance that arises from a proton-proton collision with gluonic initial state [120] σ(GG → π 0 (1) → γγ) = (where s is the center of mass energy and C GG | 13 TeV = 2137 [120] denotes the luminosity coefficient for a pair of gluons as initial partons,) we can immediately calculate the diphoton cross section. First, we shall consider the simplest case with r y = 0, (namely, no coupling to top quark pair.) In this case, we estimate the diphoton cross section (5.14) Note that, in the case of r y = 0, B(π 0 (1) → γγ) is completely free from the π 0 (1) mass dependence because π (1) decays only to the massless final states, 2γ and 2G. Hence the diphoton cross section in Eq.(5.14) is controlled only by the ratio (M π 0 (1) /f π ). To survey a generic parameter space in the present model, we shall momentarily take the value of M π 0 (1) in a range from O(100) GeV [low mass] up to O(TeV) [high mass] #15 , and discuss the phenomenological constraints from the diphoton cross section of Eq.(5.14).
In the high mass case (M π 0 (1) ∼ TeV), the estimated diphoton cross section in Eq.(5.14) is compared with the 95% C.L. upper bound on fiducial cross sections as σ spin 0 ggF (γγ) ≡ σ(GG → [spin-0 resonance] → γγ) 1 fb for narrow-width resonance with mass ∼ 1 TeV reported in Refs. [121,122]. Thus we can naively say that the prediction in the HC theory has a tight tension with the present experimental data unless M π 0 (1) /f π 3. To be consistent with the above bound for M π 0 (1) ∼ 1 TeV implies that the HC pion decay constant f π should be (at least) as large as around v VEV . However, several nonperturbative estimates with some approximation [123,124], as well as recent lattice simulations in QCD with many flavors [67,125], suggest the mass relation in magnitude of m ρ ∼ O(10)f π (for the m ρ /f π ratio in QCD with 8 flavors, see also Ref. [69].) Therefore, the diphoton constraint (M π 0 (1) /f π 3) along with the estimated mass relation (m ρ ∼ O(10)f π ) indicates the bound on the HC rho mass scale as m ρ 3 TeV. As for the low mass case, the ATLAS 8 TeV bound is available in the range 65 -200 GeV. In particular, the fluctuating 95% C.L. bounds around 65 -100 GeV is ∼ 50 fb as a crude average [126], which corresponds to ∼ 250 fb at 13 TeV [120]. Thereby the diphoton bound looks not serious for M π 0 , based on the nonperturbative observation.
The condition M π 0 (1) /f π 3 is violated, for example, if f π ∼ 100 GeV and M π 0 (1) 300 GeV. In such a case, we may still evade the 95% C.L. upper bound at 13 TeV by taking the additional contribution from the coupling to top quark (i.e. r y = 0) to make the partial width Γ(π 0 (1) → GG) reduced. Such a reduction can originate from the interference effect in Eq.(5.12). We illustrate this case in the following. In Fig. 5, we show how Γ(π 0 (1) → #15 As listed in Eq.(2.11), a typical size of the M π 0 (1) is expected to be O(100) GeV. However, the TeV mass range might be achieved when one consider possible effects from extended HC sector, which could be enhanced in the case of many flavor QCD (nearly conformal/walking gauge theory), in a way similar to extended technicolor scenarios. = 300 GeV (σ spin 0 ggF (γγ) 4 fb) is looser than that for M π 0 (1) = 400 GeV (σ spin 0 ggF (γγ) 2 fb). Thus, the interference contribution helps us to revive the possibility of π 0 (1) with the mass around 400 GeV even when f π ∼ 100 GeV, while no cancellation may be required for 300 GeV M π 0 (1) 400 GeV. This consequence is followed by discussion for a LHC bound on the HC ρ meson mass.
Another constraint comes from the scalar leptoquark search derived from a nonresonant pair production of π 0 (3) in our model, as has been reported by ATLAS in Ref. [127] for first and second generation leptoquarks. The result says that the corresponding mass scale of π 0 (3) should be 1050 − 1100 GeV at 95% C.L. assuming 100% branching fraction for a π 0 (3) decay. We note that final states are less ambiguous since possible decay branches are limited, but still there is a parameter dependence on a in general as shown in Eq.(2.41). Here, we provide a simple conclusion such that the bound from the pair production is harmless if a typical mass scale of HC pions is sufficiently greater than 1 TeV. As shown in Eq.(2.11), vector leptoquarks indeed obtain ∼ 3 TeV masses for Λ HC ∼ 1 TeV, which come most dominantly through QCD gluon exchange corrections amplified due to the characteristic feature of many flavor QCD. Thus our heavy HC pion scenario can still avoid the current scalar leptoquark bound. #16 For NHC ≥ 4, the cancellation still works for greater ry values to realize the complete cancellation in Γ(π 0 (1) → GG).

Resonant productions of HC ρ
In this part, we discuss constraints from resonance searches at the 13 TeV LHC, which provide stringent bounds on mass scales of HC rho mesons directly.

Basic backgrounds
As we have pointed out, although the V SM -ρ mixing (in the covariantized HLS formulation of this model) generates mass splittings among HC ρ mesons in physical eigenstates, only a few-percent splittings are allowed because of the large coupling g ρ ∼ 6 as in Eq.(2.44), which would be supported from the QCD-like vector dominance. Thereby, it is a good approximation to consider all of the HC ρ components to be degenerated in the common mass scale m ρ , and interactions induced through such mass mixings to be negligible. We should also recall that the mixing angles θ D and θ L defined in Eq.(4.2) were already restricted from the EW (in Sec. 3.5) and flavor (in Sec. 4) observables. To address the b → sµ + µ − anomaly consistently with the other constraints, we found that θ D needs to be much tiny such as θ D ∼ 10 −2 -10 −3 . Thus we can take the limit of θ D → 0 for all the calculations on collider phenomena (though the minuscule values are mandatory for addressing the b → sµ + µ − anomaly.) On the other hand, it has turned out that the favored regions for θ L are categorized by two spots, θ L π/4 and θ L ∼ π/2 depending on the magnitude of the coupling g ρL ≡ g ρ g 33 L . This angle θ L determines the relative coupling strength of HC ρ mesons to 2τ and 2µ, which are cos 2 θ L and sin 2 θ L , respectively, hence the coupling to 2τ becomes dominant for θ L π/4, while the coupling to 2µ does for θ L ∼ π/2. Thus, significant collider bounds would be derived from the 2τ and 2µ channel searches, depending on the θ L . To be concrete and conservative, we shall hereafter take θ L = 0 (for the 2τ channel) and θ L = π/2 (for the 2µ channel) as the reference values corresponding to the two cases.
To calculate resonant processes, values of total decay widths are important. Possible decay branches are shared by a pair of the SM fermions and a pair of HC pions [see Appendix A]. Here, we assume that the latter case is kinematically blocked (m ρ < 2m π ), which turns out to be reasonable. To see how it works, first recall that some HC pion masses should be sufficiently as heavy as O(1) TeV, as listed in Eq.(2.11). More precisely, we see from Eq.(2.11) that when f π ∼ O(100) GeV and Λ HC ∼ 1 TeV, the mass of the colorsinglet isospin-triplet mesons is ∼ 2 TeV, while those of the colored mesons are ∼ 3 TeV (∼ 4 TeV) for color triplets (color octets). To make sure how the decay channels to HC pion pairs open, we list up the total values of the final-state particle masses (m ππ ): (c.f., Appendix A),  Figure 6. The ratio Γ ρ /m ρ in Eq.(5.15) as a function of g ρL .
Then, additional contributions to the HC rho's decay branches appear when m ρ 3 TeV at the present benchmark point, f π ∼ O(100) GeV and Λ HC ∼ 1 TeV. On the other hand, when f π is somewhat greater than ∼ 100 GeV (with a sizable explicit breaking scale m 0 F ), the HC pions becomes heavier and we may block the HC rho's decays to the HC pions consistently, keeping the relation m ρ ∼ O(10)f π intact.
Thus our assumption m ρ < 2m π may be justified even in the range m ρ 3 TeV, so that we may be able to ignore decays to HC pion pairs. Of interest enough is then that all of the physical HC ρ components have the common value in the total width as where we simply ignored tiny contributions through mixing effects. Details of partial widths are provided in appendix C. The curve of the ratio Γ ρ /m ρ as a function of g ρL is illustrated in Fig. 6.

Forms of resonant cross sections
We summarize the forms of differential production (on the solid angle Ω in the center-ofmass flame) cross section at the LHC. As we pointed out, we set the mixing angles θ D = 0 and θ L = 0 or π/2 and consider that all of the HC ρ mesons are completely degenerated. We note that in the limit of θ D = 0, the possible initial state is bb only. Due to this mass degeneracy, we should take all of the HC ρ contributions simultaneously.
First, we look at the dijet final state that originates from b (or anti-b) quark where ρ 0 (8) , ρ 3 (8) , ρ 3 (1) , ρ 3 (1) , and ρ 0 (1) contribute as intermediate states (bb → ρ's → bb). The differential cross section forms for color-singlets and color-octets are summarized as where we take all of the quarks are massless. Now, both of color-octets and color-singlets contribute in s-channels, and no interference term appears between the singlet amplitudes and those of octets. Here, the following relations hold in the hatted Mandelstam variables in the parton system where the angle θ is defined in the center-of-mass frame. The factors f bb singlet and f bb octet represent summations of all possible combinations of the couplings plus the overall color factor, Next, we go for the ditau and dimuon final states, where ρ 3 (3) , ρ 0 (3) , ρ 3 (1) , ρ 3 (1) , ρ 0 (1) are possible intermediate states (bb → ρ's → ττ /µμ). We note that in the present limit (θ D = 0 and θ L = 0 or π/2), the cross section formula for the ditau channel is exactly the same as that for the dimuon channel. In the present case, color-triplets/-singlets contribute in t/schannels, and interference terms are observed between the color-triplets and color-singlets. The differential cross section forms for color-singlets, color-triplets, and interferences are summarized as (5.26)

Results
The convolution with the (anti-)bottom quark parton distribution function (PDF) inside the proton f b/p (x, µ F ) (fb /p (x, µ F )) with the Bjorken x and the PDF factorization scale µ F is formulated as Here, the total energy squared of the present LHC s(= 13 2 TeV 2 ) is related to the total energy squared of the focused parton system s as s = x 1 x 2 s ≡ τ s. Kinematically, the minimal of the fraction τ is estimated as [130] where µ F is set as m ρ . We adopt the CTEQ6L1 PDF [131] in calculations in Mathematica with the help of a PDF parser package, ManeParse 2.0 [132], and set τ 0 as 10 −6 of the minimal value of the Bjorken x in the CTEQ6L1 PDF set, which is not far from the values shown in Eq.(5.28).
In Fig. 7, we summarize the cross sections of To calculate the numerical integration including the PDF convolution in the Divonne method [133,134], we use the CUBA package [135] with the Mathlink protocol in Mathematica. The values of cross sections were cross-checked with MadGraph5 aMC@NLO [136,137], where the UFO-style model file [138] was generated by the FeynRules package [139,140]. For estimating the acceptance A of dijet events, we generated parton-level events in MadGraph5 aMC@NLO and analyzed them in the ROOT framework [141] with the help of ExRootAnalysis, which is a part of the integrated package of MadGraph5 aMC@NLO. We obtained A 0.69 in our case, which a bit deviates from the isotropic case (A ≈ 0.6) shown in Ref. [128]. The 95% C.L. upper bounds at √ s = 13 TeV were extracted from Refs. [128] (dijet, based on 12.9 fb −1 CMS data), #17 [129] (ditau, based on 2.2 fb −1 CMS data), and [143] (dimuon, based on 36.1 fb −1 ATLAS data). The expectations for the limits after 300 fb −1 data accumulation were simply calculated by rescaling from the present bounds [128,129,143].
From Fig. 7, we see the constraints on g ρL and m ρ for the θ L = 0 case. The dijet bound in the left panel shows that no constraint is imposed on our flavor-specific HC ρ mesons if the value of g ρL is maximized so as to be consistent with Eq.(4.33), which is the combined constraints from all the appreciable flavor observables, (the corresponding maximal value is taken for each of m ρ from Eq.(4.33).) On the other hand, the ditau channel excludes a part of possibilities to have the maximized g ρL , where m ρ is greater than ∼ 1500 GeV, whereas it excludes the HC rho mass scale as m ρ 900 (2100) GeV for the fixed coupling value g ρL = 1 (2), as shown in the right panel at 95% C.L.s. Therefore, the ditau channel plays a significant role in probing this scenario at the LHC #18 .
In Fig. 8, we show the collider bound for the θ L = π/2 case from the dimuon searches. It indicates that our scenario was already tested and excluded up to m ρ = 4 (1.5) TeV at 95% C.L. when g ρL = 1 (0.5). We should keep in mind, however, that the present scenario with g ρL = 1 for the mass range 4 TeV m ρ 10 TeV can still accommodate the b → sµ + µ − anomaly as seen in Eq. (4.34). Note that the dijet bound for the θ L = π/2 case is not significant when g ρL < 1.

#17
The latest ATLAS result was reported in Ref. [142] after [128], where the constraints on W and W * scenarios do not overwhelm the bound of [128] in the range of the invariant mass, q 2 2.5 TeV. #18 We observed that the largeness of the t-channel effective coupling shown in Eqs.(5.24)-(5.26) results in the situation that the t-channel contribution becomes a major part to the cross section in the ditau and dimuon production. Thus some deviations from the present acceptance times efficiency may be expected when a dedicated collider simulation is performed. We do not take into account of this point in our ballpark estimations of current constraints and future prospects.

Summary
In this paper, we analyzed the flavorful composite physics of composite vector bosons arising from the one-family model of the vectorlike HC theory, which we dubbed the HC rho model formulated based on the HLS construction. The model structure has been unambiguously fixed by the "chiral" symmetry for the HC fermions and the HLS gauge invariance. The model involves the 63 HC rho mesons and 63 HC pions with phenomenologically rich structure, which couple to the SM gauge bosons as the essential consequence of the HLS gauge invariance. Coupling properties to the SM fermions are restricted due to the HLS gauge invariance of the vectorlike theory. The flavor-dependent couplings are required to be the third-generation fermion-philic by the flavor-dependent EW precision measurements. The flavor-universal couplings, on the other hand, hardly get limited by the EW sector constraints (oblique corrections), which is indeed due to the vectorlike model construction.
To be specific, we have found that the forward-backward asymmetry of tau lepton A (0,τ ) FB and the Z boson decay to bottom quark pair (R b ) are fairly sensitive to the flavor-dependent couplings of the HC rho mesons.
In turn, we surveyed the allowed coupling parameter space of the HC rho mesons from the relevant flavor observables by taking into account the flavor mixing structures between the second and third generations involved in the left-handed down-quarks and leptons, parametrized by the angle θ D and θ L as in Eq.(4.2). The most stringent bound come from the B 0 s -B 0 s mixing. In conjunction with other flavor observables, it has turned out that the down-quark mixing θ D has to be much tiny (but non-zero) while the lepton mixing θ L has wider allowed range depending on the mass scale of HC rho mesons and the coupling strength as seen in Fig. 1. In particular, the viable scenarios are classified into two spots with respect to θ L : one class is the case that the HC rho mesons predominantly couple to third-generation leptons (θ L π/4), and another is to the second-generation (θ L ∼ π/2), as obtained in Figs. 1 -3. The exclusion plots combined with the flavor-dependent EW precision tests for the above two cases are shown in Fig. 4. Of interest for both two scenarios is that the HC rho mesons hardly contribute to B → D ( * ) ν, which in turn implies that the ratios R D ( * ) do not significantly deviate from the SM predictions. This is essentially because of the almost complete degeneracy in the HC rho mass spectra due to the large flavor-universal coupling (g ρ ∼ 6). In contrast, the HC rho mesons with mass of TeV scale can give significant contributions to b → sµ + µ − . Hence the present scenario can achieve the large values of the Wilson coefficients for the effective operators of b → sµ + µ − with the V − A form (C 9 = −C 10 ), which can account for the present anomalies in the experimental data from LHCb, Belle, ATLAS, and CMS.
We then discussed the implications to the collider physics at LHC and found that the HC rho mesons with mass of TeV scale can be consistent with the current 13 TeV LHC data: In the case with θ L π/4, the most stringent limit comes from dijet and ditau channels, which turns out to exclude the mass up to ∼ 1 − 2 TeV for the flavored HC rho coupling g ρL = 1 − 2 [ Fig. 7]. The HC rhos in the other case with θ ∼ π/2 is, on the other hand, more severely constrained by the dimuon channel, and the mass has already been excluded up to ∼ 4 (1.5) TeV for g ρL = 1 (0.5).
Thus, our HC rho model has interesting correlations sensitive to the EW precision measurements, the B-decay anomalies and the LHC collider signatures. Through the present anatomy of our model, we can reach a definite conclusion: if the current R D ( * ) anomaly goes away, but the R K ( * ) deficit further grows to be explained by the flavorful HC rhos on TeV mass scale, then those HC rhos will show up also in the future LHC data on the ditau, dijet and dimuon channels with higher luminosity. In particular, the case with θ L ∼ π/2 is the most intriguing even when viewed from any point of EW precision, flavor and collider physics; (i) the experimental bound for b → sµ + µ − anomaly can easily be satisfied if and only if a tiny θ D is at hand, while the case with θ L π/4 is barely allowed within the 3σ range even for the tiny θ D if g ρ g 33 L (= g ρL ) is close to unity [ Fig. 2]; (ii) the EW precision tests give the severer constraint on the HC rho mass when the 68% C.L. bound is taken into account in the case with θ L π/4, which is much milder in the case with θ L ∼ π/2 [ Fig. 4]; (iii) the dimuon channel search at the LHC generically has higher sensitivity than the ditau channel, so does the case with θ L ∼ π/2 [Figs. 7 and 8].
It is also interesting to note that the net effects on the flavor observables in the present model are similar to those in a low-energy effective U (1) model, in spite of the fact that our model includes 63 components of new vector bosons. This can be checked by looking at the analytic formulae for the observables presented in this paper. To be specific, the U (1) model can be described, in terms of the HC rho model, by replacing the coefficient of the four-fermion operators (q L γ µ q L )(q L γ µ q L ), (q L γ µ q L )(¯ L γ µ L ), (¯ L γ µ L )(¯ L γ µ L ): [for (a, b) = (q, q), (q, ), ( , )], see e.g., Ref. [144] and c.f., Eqs.(4.17), (4.21), (4.23), (4.25), (4.28). Therefore, the HC rho model can be considered as one of the UV completed models that realize the extended U (1) gauge symmetry to the SM at the low energy scale. Unlike the U (1) model, on the other hand, vector leptoquark bosons are involved in the present model. Therefore, individual searches for signals of the vector leptoquark bosons would be of great importance to distinguish our scenario with the low energy U (1) model that only contains Z [145]. Besides the flavorful HC rhos and HC pions, the one-family model of the HC predicts a number of other HC hadrons like composite scalars and baryons. All those HC hadrons are expected to have the same order of the mass as the HC rhos (on the order of TeV scale), and thus are potentially sensitive enough to be detected at the LHC. The lightest HC baryon might be a candidate of dark matter, because of the stability by the HC baryon number conservation. Those interesting issues are to be discussed in another publication.
In closing, we briefly sketch how to discriminate our scenario based on the vectorlike confinement and other proposals based on composite Higgs scenarios [50][51][52][53][54][55]. Two general guiding principles can be proposed. One is to measure the properties of the observed 125 GeV Higgs boson precisely. At the leading order no deviation is expected in Higgs couplings in our scenario with the SU (2) W doublet fundamental Higgs boson, while deviations are expected to be observed in Higgs couplings in composite Higgs scenarios (see e.g., [146][147][148].) The other is to clarify the species of vector-ρ mesons, which depends on how global symmetries break down, at the LHC. For example in the composite Higgs scenario discussed in Ref. [54] (which follows discussions in Refs. [149,150]), the original global symmetry is SU (4) × SU (2) L × SU (2) R × U (1) X while that of ours is SU (8). A possible difficulty is that expected spectra of these vector-ρ mesons may be fairly degenerated to evade the bounds from the electroweak precision measurements as discussed in section 3.5. Thereby, more dedicated discussions are required to declare how relevant this way is for discriminating models with (hidden) strong dynamics.

Note added:
The new LHCb measurement of R D * has been reported as in Ref. [151] after when we submitted the first version of this work on the arXiv.org. Its result is consistent with the SM prediction within ∼ 1σ, but a naive private combination with the other results leads R exp D * /R SM D * ≈ 1.21 ± 0.07. Although the deviation is unchanged, the central value becomes smaller, which is good direction for the HC ρ model.

Acknowledgments
This work was supported in part by the JSPS Grant-in-Aid for Young Scientists (B) #15K17645 (S.M.). K.N. thanks Dong Woo Kang, Pyungwon Ko, Takaaki Nomura, Tomoki Nosaka, Seong Chan Park, Yasuhiro Yamamoto and Naoki Yamatsu for useful discussions, and especially Takuya Morozumi for various fruitful comments and pointing out a typo. K.N. is also grateful to David London and the Group of Particle Physics of Université de Montréal for the kind hospitality at the middle stage of this work. R.W. thanks Ryuichiro Kitano for a discussion at very early stage of this work.
A The ρ − π − π coupling terms In this Appendix we give the explicit formulae for the ρ-π-π coupling terms derived from Eq. ( where we have used ξ R,L = e ±iπ/fπ in the unitary gauge of the HLS. From Eq.(A.1) [or Eq.(2.40) equivalently], we derive concrete forms of the ρ-π-π couplings.
Here, αβγ is the SU (2) W antisymmetric tensor, f abc (antisymmetric) and d abc (symmetric) describe the SU (3) c group structure.

B Decomposition of four-fermion currents
In this Appendix we derive the four-fermion operators induced from the HC rho meson exchanges, which are all relevant to flavor physics study in the text. By integrating out the HC rho mesons coupled to the SM fermions with the g ij L in Eq.(2.24) and using Eq.(2.9), the following effective Lagrangian forms are obtained: with T a ≡ λ a /2. Here, the factor ∆ IK;JL describes combinatorics factor, which is defined as #19 ∆ ik;jl = ∆ ki;jl = ∆ ik;lj = ∆ ki;lj = 1/2 for i = k and j = l, 1 for others.

(B.5)
With the help of the following relations in [80,152], #19 In general, the effective current-current interaction by the exchange of vector particle should be defined where MV means the mass of the vector boson and Jµ = i,j εijf i R/L γµf j R/L with charges εi,j (ignoring the SU (3)C and SU (2)W generators to be shown). Here, the overall minus sign originate from the vector-boson propagator, and the factor two is introduced to compensate the combinatoric factor of two from expanding the quadratic form of JµJ µ or deriving corresponding Feynman rules. In any case, no additional factor of two is found at the corresponding amplitudes. The factor ∆ ik;jl should be introduced to derive correct combinatoric factor in the case that the two fermion bi-linear forms are the same, where no factor of two emerges at the stage of expanding the quadrature, while the factor appears at the stage of deriving Feynman rules. = C [3] q i q j q k q l (q i L γ µ σ α q j L )(q k L γ µ σ α q l L ) + C [3] l i l j l k l l (l i L γ µ σ α l j L )(l k L γ µ σ α l l L ) + C [3] q i q j l k l l (q i L γ µ σ α q j L )(l k L γ µ σ α l l L ) + C [1] q i q j q k q l (q i L γ µ q j L )(q k L γ µ q l L ) + C [1] l i l j l k l l (l i L γ µ l j L )(l k L γ µ l l L ) + C [1] q i q j l k l l (q i L γ µ q j L )(l k L γ µ l l L ), (B.9) with C [3] q i q j q k q l = ∆ ik;jl 1 2 1 2 α il;kj − 1 6 α ij;kl 1 (M ρ α where the two factors α ij;kl and β ij;kl are defined as α ij;kl ≡ g 2 ρ g ij L g kl L , (B.16) Here, the overall minus sign has a clear nature of vector-boson exchanges, namely, originating from the vector propagator. However, to avoid clumsy negative signs in the contexts of the Wilson coefficients, we do not include these signs in the definition of the C's. Note that all of the coefficients have mass dimensions of two. The definition, (q i L γ µ σ α q j L )(l k L γ µ σ α l l L ) = 2(u i L γ µ d j L )(e k L γ µ ν l L ) + 2(d i L γ µ u j L )(ν k L γ µ e l L ) (l i L γ µ l j L )(l k L γ µ l l L ) = (ν i L γ µ ν j L )(ν k L γ µ ν l L ) + (ν i L γ µ ν j L )(e k L γ µ e l L ) + (e i L γ µ e j L )(ν k L γ µ ν l L ) + (e i L γ µ e j L )(e k L γ µ e l L ), (B.23) (l i L γ µ σ α l j L )(l k L γ µ σ α l l L ) = 2(ν i L γ µ e j L )(e k L γ µ ν l L ) + 2(e i L γ µ ν j L )(ν k L γ µ e l L ) + (ν i L γ µ ν j L )(ν k L γ µ ν l L ) + (e i L γ µ e j L )(e k L γ µ e l L ) − (ν i L γ µ ν j L )(e k L γ µ e l L ) − (e i L γ µ e j L )(ν k L γ µ ν l L ), (B.24) where u, d, ν, e are up-type, down-type quarks, neutrinos, charged leptons, respectively.
Here we shall extract operators relevant for our discussion in the next stage, −L eff ⊃ C [1] q i q j l k l l + C [3] q i q j l k l l (d i L γ µ d j L )(e k L γ µ e l L ) + C [1] q i q j l k l l − C [3] q i q j l k l l (d i L γ µ d j L )(ν k L γ µ ν l L ) + 2C [3] q i q j l k l l (u i L γ µ d j L )(e k L γ µ ν l L ) + (d i L γ µ u j L )(ν k L γ µ e l L ) + C [1] q i q j q k q l + C [3] q i q j q k q l (d i L γ µ d j L )(d k L γ µ d l L ) + C [1] l i l j l k l l + C [3] l i l j l k l l (e i L γ µ e j L )(e k L γ µ e l L ). (B.25) When (g L ) ij takes unsuppressed generic form, as widely known, such possibilities are immediately discarded by severe constraints on flavor-changing neutral current (FCNC) processes, especially the K 0 -K 0 mixing. We adopt the safety setup in Eq.(3.23) adopted in e.g., [80], which leads to the relations, where m ρ represents the universal vector boson mass and the coupling g 33 L is a real quantity. It would be noted that the vanishing condition for C [3] qqll is, concretely speaking, M ρ α (1) = M ρ α (1) and M ρ α .
In the current setup where only 2 ↔ 3 flavor transition exists, it suffices to consider the effective Hamiltonians for b → s + − , b → sνν, b → cτ − ν, τ → µss, B 0 s (= sb) ↔ B 0 s (= bs), τ − → µ − µ + µ − . The ones for the first three categories are (B.28) where α and G F are the QED fine structure constant and the Fermi constant, respectively; the Wilson coefficients include both of the SM and new physics (NP) contributions: C X = C X (SM) + C X (NP); the coefficients C IJ 9 , C IJ 10 , C IJ L , C IJ V have zero mass dimensions. Let us remind that, as defined in Eq.(4.1), the fermion fields with the prime symbol represent the mass eigenstates, where we skipped to show the symbol throughout almost all the part of the main text for clarity. We note that when only the left-type vector interactions exist as in the case of the present model, the following relation is realized, Eq.(B.25), thus brings us to the concrete results, qqll + C [3] qqll X 23 dd X IJ ll = −C IJ 10 (NP) , (B.32) qqll − C [3] qqll X 23 dd X IJ ll , (B.33) where the coefficients C IJ V (NP) become zero when the HC vector bosons are completely degenerated [c.f. Eq.(B.27)]. The effective Lagrangians which describe the remaining three processes are also discussed in Ref. [80] as L eff (τ → µss)(NP) = − C [1] qqll + C [3] qqll X 22 dd X 23 ll (µ L γ µ τ L )(s L γ µ s L ), (B.35) llll + C [3] llll X 23 ll X 22 ll (µ L γ µ τ L )(µ L γ µ µ L ) = + C [1] llll + C [3] llll sin 3 θ L cos θ L (µ L γ µ τ L )(µ L γ µ µ L ). (B.37)

C Decay widths of ρ mesons
In this Appendix we summarize the partial decay widths of the ρ mesons. Here we ignore the mixing effect of the ρ mesons and the SM vector bosons as have been done throughout the present paper [see Eq.(3.1)]. We also assume the kinematic relation of m ρ < 2m π , where no decay branch into a pair of pions is possible. This situation could be realized in the case of HC with many flavors such as the present one-family model (See Eq.(2.11)). We use the short-hand notation Γ 0 = g 2 ρL m ρ /8π and we claim that all of the ρ mesons take the universal mass m ρ by ignoring small mass splittings due to the ρ-V SM mixing. We treat all of the final-state fermions as massless particles. where we note that (ρ + (3) ) * = (ρ − (3) ). The electromagnetic charges of ρ + (3) , ρ − (3) , ρ 3 (3) , ρ 0 (3) are +5/3, −1/3 2/3, 2/3, respectively.