High-pT Signatures in Vector-Leptoquark Models

TL;DR

This paper analyzes collider signatures of TeV-scale vector bosons, including $Z'$, $U_1$ leptoquark, and $G'$, motivated by B-meson decay anomalies, deriving bounds on their masses and couplings.

Contribution

It provides a comprehensive analysis of collider signatures for three vector boson representations in models addressing B-meson decay anomalies, including interaction assumptions and bounds.

Findings

Derived bounds on masses and couplings of $Z'$, $U_1$, and $G'$

Constrained a wide range of new physics models

Analyzed interactions including right-handed and off-diagonal couplings

Abstract

We present a detailed analysis of the collider signatures of TeV-scale massive vector bosons motivated by the hints of lepton flavour non-universality observed in $B$-meson decays. We analyse three representations that necessarily appear together in a large class of ultraviolet-complete models: a colour-singlet ($Z'$), a colour-triplet (the $U_1$ leptoquark), and a colour octet ($G'$). Under general assumptions for the interactions of these exotic states with Standard Model fields, including in particular possible right-handed and flavour off-diagonal couplings for the $U_1$, we derive a series of stringent bounds on masses and couplings that constrain a wide range of explicit new-physics models.

Tables

Table 1. Table 1: Summary of the relevant experimental constraints. All searches have a centre of mass energy of 13 TeV.

Constrained BSM amplitude	Final state	Data set	Section	Reference
$U_1$ pair prod.	$b\overline{b}{\tau }^{-}{\tau }^{+}$	CMS, 35.9 fb-1	4.1	[41]
$U_1$ pair prod.	$t\overline{t}{\nu }_{\tau }\overline{{\nu }_{\tau }}$	CMS, 35.9 fb-1	4.1	[42]
$G^{\prime }$ pair prod.	$2b2\overline{b}$	CMS, 35.9 fb-1	4.1	[43]
$U_1$ [$t$ chan.] & $Z^{\prime }$ [$s$ chan.]	${\tau }_h^{+}{\tau }_h^{-}$	ATLAS, 36.1 fb-1	4.2	[44]
$U_1$ [$t$ chan.]	${\tau }_h\nu$	CMS, 35.9 fb-1	4.3	[45]
$U_1$ [$t$ chan.] & $Z^{\prime }$ [$s$ chan.]	${\tau }_h\mu$	ATLAS, 36.1 fb-1	4.4	[46]
$G^{\prime }$ [$s$ chan.] & $Z^{\prime }$ [$s$ chan.]	$t\overline{t}$	ATLAS, 36.1 fb-1	4.5	[47]

Table 2. Table 2 : Summary of the experimental constraints on pair produced leptoquarks in the b ​ b ¯ ​ τ − ​ τ + 𝑏 ¯ 𝑏 superscript 𝜏 superscript 𝜏 b\bar{b}\tau^{-}\tau^{+} [ 41 ] and t ​ t ¯ ​ ν τ ​ ν τ ¯ 𝑡 ¯ 𝑡 subscript 𝜈 𝜏 ¯ subscript 𝜈 𝜏 t\bar{t}\nu_{\tau}\bar{\nu_{\tau}} [ 42 ] final states, assuming the leptoquarks decay solely into third generation SM particles. When κ U = 1 subscript 𝜅 𝑈 1 \kappa_{U}=1 , QCD production processes become less important and lepton exchange (which depends on g U subscript 𝑔 𝑈 g_{U} ) is relevant. We thus show how the limit varies in the range g U ∈ [ 0 , 4 ] subscript 𝑔 𝑈 0 4 g_{U}\in[0,4] .

Parameters			$b\overline{b}{\tau }^{-}{\tau }^{+}$ final state			$t\overline{t}{\nu }_{\tau }\overline{{\nu }_{\tau }}$ final state
${\kappa }_U$	${\beta }_R^{33}$		$\text{BR}\left(U_1\to b{\tau }^{+}\right)$	Limit [TeV]		$\text{BR}\left(U_1\to t\overline{{\nu }_{\tau }}\right)$	Limit [TeV]
0	0		0.51	1.4		0.50	1.6
0	1		0.67	1.5		0.33	1.3
1	0		0.51	$1.1-1.3$		0.49	$1.1-1.2$
1	1		0.67	$1.2-1.4$		0.32	$1.0-1.2$

Table 3. Table 3 : Summary of the experimental cuts for the ATLAS τ h ​ τ h subscript 𝜏 ℎ subscript 𝜏 ℎ \tau_{h}\tau_{h} search  [ 44 ] . For the leading τ h subscript 𝜏 ℎ \tau_{h} we use the p T subscript 𝑝 𝑇 p_{T} cut p T τ h ​ 1 > 130 superscript subscript 𝑝 𝑇 subscript 𝜏 ℎ 1 130 p_{T}^{\tau_{h1}}>130  GeV as quoted in the HEPData entry for Ref.  [ 44 ] . Note that the corresponding cut was p T τ h ​ 1 > 85 superscript subscript 𝑝 𝑇 subscript 𝜏 ℎ 1 85 p_{T}^{\tau_{h1}}>85  GeV for 10 % percent 10 10\% of the data.

Particle selection	At least two ${\tau }_h$’s and no electrons or muons
Charge	${\tau }_{h1}$ ${\tau }_{h2}$ should be of opposite charge
${\tau }_h$ $p_T$	$p_T^{{\tau }_{h1}}>130$ GeV, $p_T^{{\tau }_{h2}}>65$ GeV
$\eta$	excluding
$\varphi$	$\left|\Delta \varphi ({\tau }_{h1},{\tau }_{h2})\right|>2.7$ rad

Table 4. Table 4 : Summary of the experimental cuts for the CMS τ h ​ ν subscript 𝜏 ℎ 𝜈 \tau_{h}\,\nu search  [ 45 ] .

Particle selection 1	No events with an electron ($p_T^e>20$ GeV, )
Particle selection 2	No events with a muon ($p_T^{\mu }>20$ GeV, )
${\tau }_h$ $p_T$	$p_T^{{\tau }_h}>80$ GeV
Missing energy	$E_T^{\mathrm{miss}}>200$
$p_T$ vs missing energy
$\varphi$	$\left|\Delta \varphi (p_T^{{\tau }_h},p_T^{\mathrm{miss}})\right|>2.4$ rad

Table 5. Table 5 : Summary of the experimental cuts for the ATLAS τ h ​ μ subscript 𝜏 ℎ 𝜇 \tau_{h}\,\mu search  [ 46 ] .

Particle selection	One single $\tau$ and $\mu$, no electrons
$p_T$	$p_T^{{\tau }_h}>65$ GeV, $p_T^{\mu }>65$ GeV
$\eta$	excluding ;
$\varphi$	$\left|\Delta \varphi ({\tau }_h,\mu )\right|>2.7$ rad
$\Delta R$	$\Delta R({\tau }_h,\mu )>0.4$

Equations

$${\mathcal{L}}_{\mathrm{NP}}\supset \frac{g_U}{\sqrt{2}}{U}_1^{\mu ,\alpha }\left[(J_U^L{)}_{\mu }^{\alpha }+(J_U^R{)}_{\mu }^{\alpha }\right]+h.c.\text{leavevmode},$$

$$(J_U^L{)}_{\mu }^{\alpha }={\beta }_L^{ij}{\bar{q}}_L^{i,\alpha }{\gamma }_{\mu }{\ell }_L^j\text{leavevmode},(J_U^R{)}_{\mu }^{\alpha }={\beta }_R^{ij}{\bar{d}}_R^{i,\alpha }{\gamma }_{\mu }{e}_R^j\text{leavevmode}.$$

$$\displaystyle\psi_{L}=\psi_{L}^{\rm SM}+\psi_{L}^{\rm exotic}\leavevmode\nobreak\ ,\qquad\psi^{\rm SM}_{L}=\left(\begin{array}[]{c}q^{\beta}_{L}\\ \ell_{L}\end{array}\right)\leavevmode\nobreak\ .$$

$$T^{\alpha}_{+}=\left(\begin{array}[]{cc}0&\delta_{\alpha\beta}\\ 0&0\end{array}\right)\leavevmode\nobreak\ .$$

$$T_{B-L}=\left(\begin{array}[]{cc}\frac{1}{3}\delta_{\beta\gamma}&0\\ 0&-1\end{array}\right)\leavevmode\nobreak\ ,\qquad\frac{1}{3}\sum_{\alpha,\delta=1}^{3}[T^{\alpha}_{+},T^{\delta}_{-}]=T_{B-L}\leavevmode\nobreak\ .$$

$$G_{\mathrm{NP}}^{\min }=SU(4)\times SU(2{)}_L\times U(1{)}_{T_R^3}\text{leavevmode},$$

$$(G_{\mathrm{NP}}^{\min }{)}^{\prime }=SU(4)\times SU(3{)}^{\prime }\times SU(2{)}_L\times U(1{)}_{T_R^3},$$

$$⟨0|{\bar{\chi }}_q^{\alpha }{\gamma }_{\mu }{\chi }_{\ell }|U_1(p,ϵ)⟩=F_U{ϵ}_{\mu }^{\alpha }\text{leavevmode},$$

$$⟨0|{\bar{\chi }}_q^{\alpha }{T}_{\alpha \beta }^a{\gamma }_{\mu }{\chi }_q^{\beta }|G(p,ϵ)⟩=F_{G^{\prime }}{ϵ}_{\mu }^a\text{leavevmode},⟨0|{\bar{\chi }}_q^{\alpha }{\gamma }_{\mu }{\chi }_q^{\alpha }|Z_q^{\prime }(p,ϵ)⟩=F_{Z_q^{\prime }}{ϵ}_{\mu }\text{leavevmode},⟨0|{\bar{\chi }}_{\ell }{\gamma }_{\mu }{\chi }_{\ell }|Z_{\ell }^{\prime }(p,ϵ)⟩=F_{Z_{\ell }^{\prime }}{ϵ}_{\mu }\text{leavevmode}.$$

$${\mathcal{L}}_{U_1}$$

$$-i{g}_Y\frac{2}{3}(1-{\tilde{\kappa }}_U)U_{1\mu }^{†}{U}_{1\nu }{B}^{\mu \nu }+\frac{g_U}{\sqrt{2}}[U_1^{\mu }({\beta }_L^{ij}{\bar{q}}_L^i{\gamma }_{\mu }{\ell }_L^j+{\beta }_R^{ij}{\bar{d}}_R^i{\gamma }_{\mu }{e}_R^j)+h.c.],$$

$${\mathcal{L}}_{Z^{\prime }}$$

$$+{\zeta }_d^{ij}{\bar{d}}_R^i{\gamma }_{\mu }{d}_R^j-3{\zeta }_{\ell }^{ij}{\bar{\ell }}_L^i{\gamma }_{\mu }{\ell }_L^j-3{\zeta }_e^{ij}{\bar{e}}_R^i{\gamma }_{\mu }{e}_R^j),$$

$${\mathcal{L}}_{G^{\prime }}$$

$$+g_{G^{\prime }}{G}^{\prime a\mu }({\kappa }_q^{ij}{\bar{q}}_L^i{T}^a{\gamma }_{\mu }{q}_L^j+{\kappa }_u^{ij}{\bar{u}}_R^i{T}^a{\gamma }_{\mu }{u}_R^j+{\kappa }_d^{ij}{\bar{d}}_R^i{T}^a{\gamma }_{\mu }{d}_R^j),$$

$$q_L^i=\left(\begin{array}{c}{V}_{ji}^{\ast }{u}_L^j\\ d_L^i\end{array}\right),{\ell }_L^i=\left(\begin{array}{c}{\nu }_L^i\\ e_L^i\end{array}\right),$$

$$\begin{array}{rlrlrl}{\beta }_L& =\left(\begin{array}{ccc}0& 0& {\beta }_L^{13}\\ 0& 0& {\beta }_L^{23}\\ 0& {\beta }_L^{32}& {\beta }_L^{33}\end{array}\right),& & & {\beta }_R& =\mathrm{diag}(0,0,{\beta }_R^{33}),\\ {\zeta }_{\ell }& =\left(\begin{array}{ccc}0& 0& 0\\ 0& {\zeta }_{\ell }^{22}& {\zeta }_{\ell }^{23}\\ 0& ({\zeta }_{\ell }^{23}{)}^{\ast }& {\zeta }_{\ell }^{33}\end{array}\right),& & & {\zeta }_e& =\mathrm{diag}(0,{\zeta }_e^{22},{\zeta }_e^{33}),\\ {\zeta }_Q& =\mathrm{diag}({\zeta }_Q^{ll},{\zeta }_Q^{ll},{\zeta }_Q^{33}),& & & {\kappa }_Q& =\mathrm{diag}({\kappa }_Q^{ll},{\kappa }_Q^{ll},{\kappa }_Q^{33}).\end{array}$$

$$\displaystyle\frac{R_{D}-R_{D}^{\mathrm{SM}}}{R_{D}^{\mathrm{SM}}}\approx 0.2\times\left(\frac{g_{U}}{3}\right)^{2}\times\left[0.3+0.7\left(\frac{\beta_{L}^{23}}{0.1}\right)\right]\times\left\{\begin{array}[]{l l}\left(\frac{2.2\leavevmode\nobreak\ {\rm TeV}}{M_{U}}\right)^{2}\quad{\rm for}\quad\beta_{R}^{33}=0&\\[3.0pt] \left(\frac{4.2\leavevmode\nobreak\ {\rm TeV}}{M_{U}}\right)^{2}\quad{\rm for}\quad\beta_{R}^{33}=-1&\end{array}\right.\,.$$

$$m_T^{\mathrm{tot}}\equiv \sqrt{(p_T^{{\tau }_{h1}}+p_T^{{\tau }_{h2}}+E_T^{\mathrm{miss}}{)}^2-({\vec{p}}_T^{{\tau }_{h1}}+{\vec{p}}_T^{{\tau }_{h2}}+{\vec{p}}_T^{\mathrm{miss}}{)}^2},$$

$$m_T^{\mathrm{tot}}\equiv \sqrt{2p_T^{{\tau }_h}{p}_T^{\mathrm{miss}}[1-\cos \Delta \varphi ({\vec{p}}_T^{{\tau }_h},{\vec{p}}_T^{\mathrm{miss}})]},$$

$$\mathcal{B}({\rm Y}(2S)\to \tau \mu )\approx 2\cdot 10^{-6}\times {\left(\frac{g_U}{3}\right)}^4\times {\left(\frac{1\mathrm{TeV}}{M_U}\right)}^4\times |{\beta }_L^{32}{|}^2.$$

$$M_{G^{\prime }}$$

Full text

ZU-TH-19/05

High-$p_{T}$ Signatures in Vector–Leptoquark Models

Michael J. Baker, Javier Fuentes-Martín, Gino Isidori, Matthias König

Physik-Institut, Universität Zürich, CH-8057 Zürich, Switzerland

Abstract

We present a detailed analysis of the collider signatures of TeV-scale massive vector bosons motivated by the hints of lepton flavour non-universality observed in $B$-meson decays. We analyse three representations that necessarily appear together in a large class of ultraviolet-complete models: a colour-singlet ($Z^{\prime}$), a colour-triplet (the $U_{1}$ leptoquark), and a colour octet ($G^{\prime}$). Under general assumptions for the interactions of these exotic states with Standard Model fields, including in particular possible right-handed and flavour off-diagonal couplings for the $U_{1}$, we derive a series of stringent bounds on masses and couplings that constrain a wide range of explicit new-physics models.

1 Introduction

The hints of Lepton Flavour Universality (LFU) violation in semi-leptonic $B$ decays, namely the deviations from $\tau/\mu$ (and $\tau/e$) universality in $b\to c\ell\bar{\nu}$ decays [1, 2, 3, 4] and the deviations from $\mu/e$ universality in $b\to s\ell\bar{\ell}$ decays [5, 6], are among the most interesting departures from the Standard Model (SM) reported by experiments in the last few years. The attempt to find a single beyond-the-SM (BSM) explanation for the combined set of anomalies has triggered intense theoretical activity, whose interest goes beyond the initial phenomenological motivation. In fact, it has shed light on new classes of SM extensions that turn out to be very interesting per se and that have not been investigated in great detail so far, pointing to non-trivial dynamics at the TeV scale possibly linked to a solution of the SM flavour puzzle.

The initial efforts to address both sets of anomalies have been focused on Effective Field Theory (EFT) approaches via four-fermion effective operators (see [7, 8, 9, 10] for the early attempts). However, the importance of complementing EFT approaches with appropriate simplified models with new heavy mediators was soon realised [9, 11]. Given the relatively low scale of new physics hinted by the charged-current anomalies, the impact of considering a full model rather than an EFT on high-$p_{T}$ constraints are significant [12, 13, 14]. More recently, a further advancement has been achieved with the development of more complete (and more complex) models with a consistent ultraviolet (UV) behaviour (see in particular [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27]).

In early EFT attempts, it was realised that a particularly good mediator accounting for both sets of anomalies is a TeV-scale $U_{1}\sim(\mathbf{3},\mathbf{1},2/3)$ vector leptoquark, coupled mainly to third-generation fermions [8, 11]. The effectiveness of this state as a single mediator accounting for all available low-energy data has been clearly established in [28]. However, this state can not be the only TeV-scale vector particle in a realistic extension of the SM. Since it is a massive vector, the $U_{1}$ can be either a massive gauge boson of a spontaneously broken gauge symmetry $G_{\rm NP}\supset G_{\rm SM}$, as in the attempts proposed in [15, 16, 17, 18], or a vector resonance of some new strongly interacting dynamics, as e.g. in [19, 21]. As we show, in both cases the consistency of the theory requires additional vector states with similar masses. The purpose of this paper is to provide a comprehensive analysis of the high-$p_{T}$ constraints on the vector leptoquark $U_{1}$ and what can be considered its minimal set of vector companions, namely a colour octet $G^{\prime}\sim(\mathbf{8},\mathbf{1},0)$, which we will refer to as the coloron, and a colour singlet $Z^{\prime}\sim(\mathbf{1},\mathbf{1},0)$.

In our analysis we consider the most general chiral structure for the $U_{1}$ couplings to SM fermions. This is in contrast with many recent studies which considered only left-handed (LH) couplings. While this hypothesis is motivated by the absence of clear indications of right-handed (RH) currents in the present data and by the sake of minimality, it does not have a strong theoretical justification. Indeed, the quantum numbers of the $U_{1}$ allow for RH couplings, and in motivated UV completions such couplings naturally appear [18, 26]. We also analyse the impact of a non-vanishing mixing between the second and third family in high-$p_{T}$ searches, including in particular constraints from $pp\to\tau\mu$ and $pp\to\tau\nu$. As we show, the inclusion of right-handed couplings and/or a sizeable 2-3 family mixing yields significant modifications to the results found in the existing literature.

The structure of this paper is as follows: In section 2 we motivate our choice of TeV-scale vectors and in section 3 we introduce the phenomenological Lagrangian adopted to describe their high-$p_{T}$ signatures. We then present the results of the searches in section 4 and conclude with section 5.

2 The spectrum of vector states at the TeV scale

The bottom-up requirement for the class of models we are interested in is the following effective interaction of the $U_{1}$ field with SM fermions:

[TABLE]

Here $q_{L}\,(\ell_{L})$ denotes the left-handed quark (lepton) doublets, $d_{R}\,(e_{R})$ denotes the right-handed down-type quark (charge-lepton) singlets, $i\in\{1,2,3\}$ and $j\in\{1,2,3\}$ are flavour indices, $\alpha\in\{1,2,3\}$ is a $SU(3)_{c}$ index, and $\beta_{L,R}^{ij}$ are complex matrices in family-space.

The effective interaction in eq. 1 unambiguously identifies the representation of $U_{1}$ under $G_{\rm SM}=SU(3)_{c}\times SU(2)_{L}\times U(1)_{Y}$ to be $(\mathbf{3},\mathbf{1},2/3)$. There are two basic classes of well-defined UV theories where such interactions can occur:

1. i.

Gauge models. Here $U_{1}$ is the massive gauge boson of a spontaneously broken gauge symmetry $G_{\rm NP}\supset G_{\rm SM}$. The need for extra massive vectors follows from the size of the coset-space of $G_{\rm NP}/G_{\rm SM}$, that necessarily requires additional generators besides the six associated to $U_{1}$. 2. ii.

Strongly interacting models. Here $U_{1}$ appears as a massive resonance for a new strongly interacting sector. In this case the need of additional massive vectors is a consequence of the additional resonances formed by the same set of constituents leading to $U_{1}$.

Gauge models: the need for a $Z^{\prime}$

Within gauge models, let us start analysing the case of a single generation of SM fermions ($i=j=3$), and further assume that SM fermions belong to well-defined representations of $G_{\rm NP}$ (i.e. no mixing between SM-like and exotic fermions). Under these assumptions, $\beta_{L}$ is non-zero only if $q_{L}$ and $\ell_{L}$ belong to the same $G_{\rm NP}$ representation. We denote this representation $\psi_{L}$ and, without loss of generality, we decompose it as

[TABLE]

In this notation the left-handed current in eq. 1 can be written as $(J^{L}_{U})_{\mu}^{\alpha}=\bar{\psi}^{\rm SM}_{L}(T^{\alpha}_{+})\gamma_{\mu}\psi^{\rm SM}_{L}$ with the following explicit expression for the action of the $G_{\rm NP}$ generators on the SM projection of $\psi_{L}$:

[TABLE]

The closure of the algebra of the six generators $T^{\alpha}_{\pm}$ associated with the six components of $U_{1}$ implies the need of the following additional (colour-neutral) generator

[TABLE]

The same conclusion is reached by looking at the right-handed coupling in eq. 1. Moreover, since a possible mixing between SM and exotic fermions must occur in a $SU(3)_{c}$ invariant way, the decomposition in eq. 4 also holds for possible exotic fermions mixing with the SM ones. Hence the need of $T_{B-L}$ for the closure of the algebra is a general conclusion that holds independently of the possible mixing among fermion representations.

An equivalent way to deduce the need for an extra generator is the observation that the minimal group $G^{\rm min}_{\rm NP}\supset G_{\rm SM}$ containing generators associated to the representation $(\mathbf{3},\mathbf{1},2/3)$ is

[TABLE]

i.e. the subgroup of the Pati-Salam group $G_{\rm PS}=SU(4)\times SU(2)_{L}\times SU(2)_{R}$ [29]. $G_{\rm NP}^{\rm min}$ is obtained by considering the $U(1)$ subgroup of $SU(2)_{R}$ defined by its diagonal (electric-charge neutral) generator $T^{3}_{R}$. The coset $G^{\rm min}_{\rm NP}/G_{\rm SM}$ contains seven generators: the six $T^{\alpha}_{\pm}$ describing the coset $SU(4)/SU(3)_{c}\times U(1)_{B-L}$, and $T_{B-L}$.

In gauge models, the presence of an extra massive vector $Z^{\prime}\sim(\mathbf{1},\mathbf{1},0)$ associated with the breaking $U(1)_{B-L}\times U(1)_{T^{3}_{R}}\to U(1)_{Y}$ is thus unavoidable. Since the breaking of $U(1)_{B-L}$ necessarily implies a breaking of $SU(4)$, the breaking terms which lead to a non-vanishing $Z^{\prime}$ mass necessarily induce a mass term for the $U_{1}$ as well. Hence, the $Z^{\prime}$ state cannot be decoupled. The opposite is not true: since the $U_{1}$ generators are associated to the $SU(4)/SU(3)_{c}\times U(1)_{B-L}$ coset, mass terms for the $U_{1}$ do not necessarily contribute to the $Z^{\prime}$ mass.

Gauge models: the need for a $G^{\prime}$

While the minimal group in eq. 7 allows us to build a consistent model for a massive $U_{1}\sim(\mathbf{3},\mathbf{1},2/3)$, it does not leave us enough freedom to adjust $U_{1}$ and $Z^{\prime}$ couplings in order to comply with low- and high-energy data.

Under $G_{\rm NP}^{\rm min}$ the interaction strengths of both $U_{1}$ and $Z^{\prime}$ are unambiguously related to the QCD coupling ($g_{s}$) and to hypercharge, given that they all originate from the same $SU(4)$ group. In particular $g_{U}=g_{s}(M_{U_{1}})$, in a normalisation where $|\beta_{L,R}^{ij}|\leq 1$. Moreover, the couplings of the $Z^{\prime}$ to SM fermions are necessarily flavour universal.111This statement follows from the fact that the mixing of SM fermions among themselves (in flavour space) and with possible exotic representations necessarily involve states with the same $B-L$ charges. As a result, the mixing acts as a unitary rotation on the $Z^{\prime}$ couplings that remains proportional to the identity matrix in flavour space. A flavour-universal $Z^{\prime}$ is constrained by LHC dilepton searches to have M_{Z^{\prime}}\lower 3.01385pt\hbox{;\stackrel{{\scriptstyle\textstyle>}}{{\sim}};}5 TeV [30, 31]. Within $G_{\rm NP}^{\rm min}$, the $U_{1}$ should be necessarily close in mass [22] which, together with the low value of $g_{U}$, results in a negligible impact on $b\to c\ell\nu$ decays.

To avoid these constraints, $T^{\alpha}_{\pm}$, $T_{B-L}$, and the QCD generators $T^{a}$, should not be unified in a single $SU(4)$ group. Given the commutation rules between $T^{\alpha}_{\pm}$ and $T^{a}$, the next-to-minimal option is obtained with [15]

[TABLE]

where $SU(3)_{c}$ is the diagonal subgroup of $SU(4)\times SU(3)^{\prime}$ (see also [32, 33]). In this case we can achieve the two goals of 1) decoupling the overall coupling of $U_{1}$ from $g_{c}$, letting it reach the higher values needed to impact $B$-physics data with $M_{U_{1}}\sim$ few TeV; 2) having flavour non-universal couplings for both $U_{1}$ and $Z^{\prime}$. The latter can be achieved either via mixing with exotic fermions (as in [15]), and/or with a flavour-dependent assignment of the $SU(4)\times SU(3)^{\prime}$ quantum numbers (as in [18, 20]).

The enlargement of the coset space to $(G_{\rm NP}^{\rm min})^{\prime}/G_{\rm SM}$ directly requires a massive colour-octet vector (the “coloron” $G^{\prime}$) associated to the breaking $SU(3)_{[4]}\times SU(3)^{\prime}$, where $SU(3)_{[4]}$ is the “coloured” subgroup of $SU(4)$. Similarly to the case of the $Z^{\prime}$, breaking terms leading to a non-vanishing $G^{\prime}$ mass necessarily induces a mass term also for the $U_{1}$, while the opposite is not necessarily true.

Vector spectrum in strongly interacting models

In strongly interacting models, the leptoquark $U_{1}$ is a composite state with two elementary fermions charged under the new confining group $G_{\rm strong}$ as constituents. These fermions are necessarily charged under $SU(3)_{c}$ in order to generate a colour-triplet state.

The simplest option is the case of a vector triplet ($\chi_{q}^{\alpha}$) and a vector singlet ($\chi_{\ell}$), both in the fundamental of $G_{\rm strong}$, such that

[TABLE]

where we have not explicitly indicated the $G_{\rm strong}$ indices. With these basic constituents one expects also one $G^{\prime}$ and two $Z^{\prime}$:

[TABLE]

The masses of these states are not precisely related to that of the $U_{1}$ as in the case of gauge models, but they are expected to be of similar size since they originate from the same dynamics. In principle one can enlarge the multiplicity of the constituents, e.g. the colour triplet $U_{1}$ can be achieved by combining ${\bf 3}$ and $\bf{8}$ of $SU(3)_{c}$, but this can only increase the number of extra coloured vectors. A further exotic option is to consider $U_{1}$ as a fermion bilinear in an antisymmetric combination of $G_{\rm strong}$, as allowed e.g. in $SU(2)$. However, beside this peculiar case where symmetric $G_{\rm strong}$ combinations are forbidden (or much heavier in mass), this does not prevent the presence of at least one $G^{\prime}$ and one $Z^{\prime}$ with masses comparable to the $U_{1}$.

3 Phenomenological Lagrangian

Having motivated the minimal set $\{G^{\prime},Z^{\prime},U_{1}\}$ of massive vectors for a meaningful description of TeV scale dynamics, we proceed to set up a versatile framework for analysing the high-$p_{T}$ signatures of these states in a general way. In our analysis we restrict our attention to the interactions of these vectors with SM fermions and gauge bosons. We neglect possible Higgs couplings to the $Z^{\prime}$ since they are severely constrained by electroweak precision data (see e.g. [28]) and are typically very small in the model realisations we are interested in. We also ignore any possible interactions of the extra vectors among themselves and to any other particles related to the UV completion of the model (either scalars or fermions). While some of the high-$p_{T}$ signatures related to these interactions can be quite interesting [22], they are highly dependent on the details of the UV completion. Here we only consider their possible indirect effects on the widths of the vectors, which we treat as an additional free parameter.222We will consistently assume that right-handed neutrinos, if present, are heavy enough so that they effectively decouple and do not play any relevant role. Models with light $\nu_{R}$ in connection to the $B$-anomalies can be found in [34, 35, 36, 37], and in connection to the vector leptoquark in [38].

We define the general Lagrangian for these vectors as follows:

[TABLE]

where $T^{a}=\lambda^{a}/2$, with $\lambda^{a}$ ($a=1,\dots,8)$ the Gell-Mann matrices. In both $\mathcal{L}_{U_{1}}$ and $\mathcal{L}_{G^{\prime}}$ we include possible non-minimal interactions with SM gauge fields, which play a role in the pair production of the heavy vectors at the LHC. In gauge models these couplings vanish, $\kappa_{U}=\tilde{\kappa}_{U}=\kappa_{G^{\prime}}=\tilde{\kappa}_{G^{\prime}}=0$. However, this is not necessarily the case in strongly interacting models. The so-called minimal-coupling scenario for the leptoquark corresponds to $\kappa_{U}=\tilde{\kappa}_{U}=1$. Since a triple coupling of the type $GGG^{\prime}$ would lead to a huge enhancement of the colour-octet production at LHC and with that, to very strong constraints from high-energy data, in what follows we will take $\kappa_{G^{\prime}}=0$.

Without loss of generality, we choose the flavour basis of the $SU(2)_{L}$ fermion doublets to be aligned to the down-quark sector, i.e.

[TABLE]

where $V_{ji}$ denote the CKM matrix elements. We assume that the new vectors are coupled dominantly to third generation fermions. The couplings to light quarks are assumed to respect a $\mathrm{U(2)}_{q}$ flavour symmetry broken only in the leptoquark sector by the same leading spurion controlling the $3\to q$ mixing in the CKM matrix [39]. We parameterise the strength of this spurion by $\beta_{L}^{23}$. In the lepton sector we assume vanishing couplings to electrons. These assumptions are phenomenologically motivated by the tight constraints from low-energy observables, in particular $\Delta F=2$ amplitudes and lepton flavour violation in charged leptons (see e.g. [28, 26]). More precisely, we take the following textures for the vector couplings ($Q=q,u,d$):

[TABLE]

As shown in [28], the assumption of a single $\mathrm{U(2)}_{q}$ breaking spurion in both leptoquark and SM Yukawa couplings implies the relation $\beta_{L}^{13}=V_{td}^{*}/V_{ts}^{*}\;\beta_{L}^{23}$. More generally, from $\mathrm{U(2)}$ symmetries acting on both quark and lepton sectors, we expect the following hierarchy: $|\beta_{L}^{31}|\ll|\beta_{L}^{23}|,\leavevmode\nobreak\ |\beta_{L}^{32}|\ll|\beta_{R}^{33}|,\leavevmode\nobreak\ |\beta_{L}^{33}|=\mathcal{O}(1)$, and analogously for the $\zeta_{\ell,e,Q}^{ij}$ and $\kappa_{Q}^{ij}$ couplings.

Since our main motivation for analysing the high-$p_{T}$ phenomenology of the $U_{1}$ is its success in addressing $B$-physics anomalies, it is worth recalling for which values of its couplings this can happen. Detailed analyses of this question in specific UV models can be found for instance in [15, 26, 22, 40]. To remain sufficiently generic, we focus here on the contribution to $R_{D}=\Gamma(B\to D\tau\nu)/\Gamma(B\to D\mu\nu)$. This is the most interesting low-energy observable to constrain the $U_{1}$ couplings relevant at high-$p_{T}$ in the presence of right-handed currents. Setting $\beta_{L}^{33}=-1$ one has [40]

[TABLE]

This expression illustrates well the dependence on the relevant couplings. For $g_{U}$ within the pertubative regime $(g_{U}<\sqrt{4\pi})$, the $U_{1}$ mass cannot exceed a few TeV. As a reference benchmark, a $20\%$ enhancement in $R_{D}$ (in good agreement with present data) is obtained for $M_{U}=2.2$ TeV if $\{g_{U},\beta_{L}^{23}\}=\{3,0.1\}$ and $\beta_{R}^{33}=0$, or $M_{U}=4.2$ TeV with the same values of $\{g_{U},\beta_{L}^{23}\}$ but setting $\beta_{R}^{33}=-1$.

4 Results

We consider a variety of high-$p_{T}$ searches at the LHC which place limits on the model discussed above. The most constraining ones, which we discuss in detail below, are shown in table 1. In some cases the searches are optimised for the BSM processes we are interested in, allowing a simple translation of the reported limits in terms model parameters. In most cases however, a reinterpretation of the reported limits is necessary.

A relatively simple case is that of the leptoquark pair production. The differential and total cross-sections for these processes are well-known [48]. Here we use the recent CMS analyses dedicated to pair-produced (scalar) leptoquarks decaying primarily to third generation SM fermions [41, 42]. Since the leptoquarks are predominantly produced via their strong couplings to gluons, the limits only depend on the branching ratios to the relevant final states. Bounds on the coloron mass are extracted from a search for pair-produced resonances decaying to quark pairs, reported by CMS in the same way [43].

The case of the $\tau^{+}\tau^{-}$ final state, which constrains both the $Z^{\prime}$ ($s$ channel production) as well as the $U_{1}$ ($t$ channel exchange), is significantly more involved. Here we re-interpret the limits on resonances decaying into tau-lepton pairs, with hadronically decaying taus, reported by ATLAS [44]333We do not consider the corresponding analysis by CMS [49], which focuses on heavy Higgs bosons. (bounds from leptonic tau decays turn out to be significantly weaker at large ditau invariant masses). We first consider the bounds placed on the $U_{1}$ and on the $Z^{\prime}$ in isolation, for various choices of couplings and widths, and then in combination. As we emphasise below, it is essential to include all relevant experimental information when deriving limits in this case.

We extract further bounds on $U_{1}$ by recasting CMS searches for $pp\to\tau\nu$ [45] and limits on both $Z^{\prime}$ and $U_{1}$ from the $pp\to\tau\mu$ search by ATLAS [46]. In both cases the 2-3 family mixing of the leptoquark plays a key role. As far as other dilepton final states are concerned, we explicitly checked that constraints from $pp\to\mu\mu$ (see e.g. [31]) do not significantly constrain the parameter space relevant to our model.

The leading bound on the $G^{\prime}$ is extracted by the unfolded $t\bar{t}$ invariant mass spectrum provided by ATLAS [47]. In principle, the $U_{1}$ and the $Z^{\prime}$ could be constrained by dijet searches. However, in our setup resonances tend to be very wide, with a width-over-mass $\sim 25\%$. As a result, the limits reported in the literature on narrow dijet peaks over a data driven background spectrum [50, 51, 52] are not directly applicable. Furthermore, dijet signatures are mostly produced for light quarks and gluons, which couple only weakly to $Z^{\prime}$ and $G^{\prime}$ in our setup.444Searches for $b$-tagged dijet signatures would remedy this, but they tend to be sensitive to relatively low mass ranges and narrow widths [53, 54]. Indeed, dedicated recasts of dijet searches performed in a setup similar to ours have shown that these constraints are less significant than those from the $t\bar{t}$ final state [22]. Although one can envision scenarios where current dijet searches are more constraining than $t\bar{t}$ searches, such as when third-generation couplings are suppressed or when light-generation couplings are large, these limits are less relevant for the class of models which fit the flavour anomalies and so we do not consider dijet searches.

To perform recasts of these searches we implement the model described in section 3 in FeynRules 2.3.32 [55] and generate the corresponding UFO model file. The FeynRules model files as well as the corresponding UFO model are available at https://feynrules.irmp.ucl.ac.be/wiki/LeptoQuark. In our Feynrules implementation and in all our results throughout this paper, we include only tree-level effects. While some NLO QCD corrections are available for the vector leptoquark case [56], in specific models these are expected to be supplemented by additional NLO contributions that can be of similar (or even larger) size. Hence we opt not to include them and we add a systematic error in our signal to (partially) account for them. Other Feynrules implementations for the vector leptoquark (but with interactions to third-generation left-handed fields only) are available [57]. We have cross-checked our leptoquark implementation (with $\beta_{R}^{33}=0$) against the one in [57], finding a perfect agreement between the two.

4.1 Limits from resonance pair production

We first briefly discuss limits on the leptoquark coming from their pair production. For a large fraction of the parameter space, the dominant production modes are governed by QCD and the relevant couplings are the strong gauge coupling and $\kappa_{U}$, see section 3. However, for $\kappa_{U}\sim 1$ the QCD-induced production cross-section is smaller and pair-production via lepton exchange becomes relevant for large values of $g_{U}$. The most constraining searches in our scenario are those for the $b\bar{b}\tau^{-}\tau^{+}$ [41] and $t\bar{t}\nu_{\tau}\bar{\nu_{\tau}}$ [42] final states.

In table 2 we report the limits for various values of $\kappa_{U}$ and $\beta_{R}^{33}$, which determines the branching ratios (the branching ratios deviate slightly from the expected 1/2, 1/3, 2/3 due to phase space effects). We assume that the leptoquark decays only into third generation SM particles and find that the limits range from 1 TeV to 1.6 TeV. Similar limits have also been obtained in the literature, see e.g. [42, 58, 33], although using lower luminosity in the $b\bar{b}\tau^{-}\tau^{+}$ channel. Whenever it is possible to compare, we find good agreement between our results and those in the aforementioned references. With $\kappa_{U}=0$ there is an extra coupling to the gluon field strength tensor boosting the production cross-section and strengthening the limit. As $\beta_{R}^{33}$ increases, the branching ratio to $b\tau^{+}$ increases while the branching ratio to $t\bar{\nu}_{\tau}$ decreases, which is reflected in a strengthening and weakening of the limits, respectively. For illustration, we include the strongest bound from pair-production, i.e the limit $M_{U}>1.6$ TeV, in figs. 1, 6 and 7.

In a similar fashion, bounds on the coloron mass can be extracted from a search for pair-produced resonances decaying to quark pairs, performed by the CMS collaboration [43]. The search excludes a coloron in the whole mass range considered, from 80 GeV to 1500 GeV, so provides an upper bound of $M_{G^{\prime}}>1.5$ TeV. However, a stronger upper bound can be estimated by extrapolating the production cross-section and exclusion limit to higher energies, where bounds of 1.7 TeV and 2.1 TeV for $\tilde{\kappa}_{G^{\prime}}=0$ and $\tilde{\kappa}_{G^{\prime}}=1$ are obtained. The stronger bound in the latter case can be understood from the fact that the corresponding operator in section 3 adds significantly to the $gg\to G^{\prime}G^{\prime}$ amplitude. The estimated limits are practically independent of the choices of the couplings to quarks, because the production cross section is dominated by the gluon-initiated processes. In setting these limits, we fix the coloron gauge coupling to $g_{G^{\prime}}=3$, $\kappa_{G^{\prime}}=0$ and $\kappa_{q,u,d}^{33}=1$.

4.2 $pp\to\tau\tau$ search

The ATLAS collaboration has performed a search of heavy resonances in the ditau final state using $36.1\leavevmode\nobreak\ \textrm{fb}^{-1}$ of 13 TeV data [44]. In this section we recast this search to set limits on the $U_{1}$ and $Z^{\prime}$ masses for different choices of the couplings. In section 4.2.3 and section 4.2.2 we consider separate limits for the $Z^{\prime}$ and $U_{1}$ assuming that one of the two has fully decoupled. The interplay of the two resonances in this search is considered at the end, in section 4.2.4.

4.2.1 Search strategy

We focus on the analysis with $\tau_{h}\tau_{h}$ since this channel presents the highest sensitivity to high-mass resonances. The contributions to the $pp\to\tau^{+}\tau^{-}$ process from new heavy resonances, including the interference with the SM, are computed using Madgraph5_aMC@NLO v2.6.3.2 [59], with the NNPDF23_lo_as_0119_qed PDF set [60]. Hadronization of the $\tau$ final states is performed with Pythia 8.2 [61] with the A14 set of tuned parameters [62]. Detector simulation is done using Delphes 3.4.1 [63]. The ATLAS Delphes card has been modified to satisfy the object reconstruction and identification requirements. In particular we include the $\tau$-tagging efficiencies quoted in the experimental search [44]. After showering and detector simulation, we apply selection cuts using MadAnalysis 5 v1.6.33 [64] (see table 3 for details on the applied cuts). We have validated our results by generating the SM Drell-Yan $pp\to\tau\tau$ background and comparing our results with the one quoted by ATLAS. A good agreement is found between the two samples (we find a discrepancy with the quoted central values of less than 20%, well within the given $1\sigma$ region).

After passing through selection cuts, the resulting events are binned according to their total transverse mass,

[TABLE]

where $p_{T}^{\tau_{h1,2}}$ are the transverse momenta of the visible decay products for the leading and sub-leading taus, respectively, and $E_{T}^{\rm\,miss}$ and $\vec{p}_{T}^{\;\rm miss}$ are the total missing transverse energy and missing momentum in the reconstructed event. We compare our binned events with the histogram in Fig. 3b of the supplementary material of [44], which contains the corresponding $m_{T}^{\rm tot}$ histograms for the SM background and the experimental data, with $b$-tag inclusive event selection. For the statistical analysis we use the modified frequentist $\textrm{CL}_{s}$ method [65]. We compute the $\textrm{CL}_{s}$ using the ROOT [66] package Tlimit [67] and exclude model parameter values with $\textrm{CL}_{s}<0.05$. In our statistical analysis we include all the bins and SM backgrounds errors, provided by the ATLAS collaboration in the corresponding HEPData entry [44].555The power of each of the bins in excluding a signal is shown in fig. 5, where we plot the 95% CL exclusion limit in the leptoquark mass, as a function of the number of the bins included in the statistical analysis. When combining the bins, we ignore any possible correlation in the bin errors, since they are not provided by the experimental collaboration. We also include a systematic uncertainty of $20\%$ for the signal to account for possible uncertainties related to the PDF, tau hadronization, detector simulation and unaccounted NLO corrections.

4.2.2 Limits on the $U_{1}$ leptoquark

In this section we decouple the $Z^{\prime}$ and concentrate on the limits arising exclusively from the leptoquark exchange. In our search we take maximal values for $\beta_{L}^{33}$ (i.e. $\beta_{L}^{33}=1$) and consider three benchmarks for the right-handed coupling: $|\beta_{R}^{33}|=\{0.0,\,0.5,\,1.0\}$. Note that the search is not sensitive to the relative sign choice between $\beta_{R}^{33}$ and $\beta_{L}^{33}$ but only to their magnitudes. The reason for this is that the New Physics (NP) amplitudes of different chiralities do not interfere with each other and the amplitude proportional to $\beta_{R}^{33}\,\beta_{L}^{33}$ does not interfere with the SM ones. We further fix the leptoquark width to its natural value. The leptoquark width only mildly affects the results of this search, contrary to the $Z^{\prime}$ case discussed in the next section, since the NP contribution is generated via a $t$ channel exchange.

Exclusion limits in the $(g_{U},M_{U})$ plane, setting $\beta^{23}_{L}=0$, are shown in fig. 1 (left). Similar recasts for the case with $\beta_{R}^{33}=0$ can be found in the literature [58, 68]. We obtain slightly stronger limits than those in the previous references. As we show in fig. 5, this difference can be understood from the fact that we consider the full $m_{T}^{\rm tot}$ distribution and not only the highest bin. The lower bins are important since a $t$ channel exchange gives rise to a broad tail in the spectrum. Exclusion limits for the scenario where $\beta_{R}^{33}\neq 0$ have not been discussed in the literature. We find that the additional chirality significantly enhances the cross section, yielding limits that are about $70\%$ stronger than in the case when $\beta_{R}^{33}=0$.

Finally, we also study the limits on $M_{U}$ for non-zero values of $\beta_{L}^{23}$, fig. 1 (right). Here we fix $g_{U}=3$ and $\beta_{L}^{33}=1$ and plot the corresponding exclusion limits for the three benchmark values of $\beta_{R}^{33}$ discussed above. As can be seen, only a mild increase of the limits is found for $\beta_{L}^{23}\lesssim 0.4$. For larger values of $\beta_{L}^{23}$, the PDF enhancement is enough to make $s\bar{s}\to\tau^{+}\tau^{-}$ the dominant partonic channel and the limits start growing linearly with $\beta_{L}^{23}$.

4.2.3 Limits on the $Z^{\prime}$ resonance

We now proceed to the limits set on the $Z^{\prime}$, decoupling the leptoquark. Throughout this section we fix $\zeta_{q,u,d}^{33}=\zeta_{\ell,e}^{33}=1$ and focus on the impact of varying the overall $Z^{\prime}$ coupling $g_{Z^{\prime}}$, varying the coupling to left-handed light quarks $\zeta_{q}^{ll}$, and varying the width of the $Z^{\prime}$.

In the left panel of fig. 2 we set $\zeta_{q}^{ll}=0$ and show the exclusion in the $(g_{Z^{\prime}},M_{Z^{\prime}})$ plane. For small couplings, $g_{Z^{\prime}}<0.5$, the $Z^{\prime}$ is not excluded above 1 TeV as the production cross section is too small. In the range $0.5<g_{Z^{\prime}}<1.0$ the limit increases from 1 TeV to 2 TeV and it approaches a regime where it increases linearly with the coupling. This can be understood by the fact that, having set $\zeta_{q}^{ll}=0$, the $Z^{\prime}$ is dominantly produced from $b$-quarks, which carry only low momentum fractions of the protons. As a result, even for relatively low masses the effective cross-section scales like a contact interaction $\sigma_{Z}^{\prime}\sim g_{Z}^{\prime 4}/M_{Z}^{\prime 4}$.

Finally, we also show the impact of varying the width. As can be noted, doubling the width (dashed line in fig. 2) has a relatively minor impact. This is consistent with the observation that the limits does not come from the on-shell production of the $Z^{\prime}$, but rather from its tail (that scales like a contact interaction).

In fig. 2 (right) we fix $g_{Z^{\prime}}=3$ and vary the couplings to left-handed light quarks $\zeta_{q}^{ll}$. Since the light quarks have less PDF suppression than the third-generation quarks, the limit increases rapidly. For $\zeta_{q}^{ll}\ll 1$, the width is not affected by increasing $\zeta_{q}^{ll}$, while for larger values of $\zeta_{q}^{ll}$ the width starts to be affected leading to a change of slope.

We again show that doubling the natural width decreases the limit by around 10 %. We also show the impact of changing the relative sign between the light quark couplings and the third-generation coupling. With opposite signs the interference term contributes constructively, strengthening the limit, whereas when the signs are the same the interference term contributes destructively, weakening the limit.

In fig. 3 we fix $g_{Z^{\prime}}=3$ and vary the width for $\zeta_{q}^{ll}\in\{0.0,\,0.5,\,1.0\}$. As noted above, we see that the limit depends only weakly on the width. For all values of $\zeta_{q}^{ll}$, a doubling of the width from 25% to 50% decreases the limit by around 10%. The grey area show values of the width which are below the corresponding natural width.

In summary, the $Z^{\prime}$ mass limit of the ditau search depends weakly on the universal coupling $g_{Z^{\prime}}$, is very sensitive to the light-quark couplings (it is excluded below 5 TeV for $\zeta_{q}^{ll}\approx 1$), and is only weakly relaxed by an increase of the total width of the $Z^{\prime}$.

4.2.4 Combined limits for the $Z^{\prime}$ and the $U_{1}$ leptoquark

We now consider the limits when both the $Z^{\prime}$ and the leptoquark are present. For the $Z^{\prime}$ we set $\zeta_{q,u,d}^{33}=\zeta_{\ell,e}^{33}=1$ and $\zeta_{q}^{ll}=0$. For the leptoquark we set $\beta_{L}^{33}=\beta_{R}^{33}=1$ and $\beta_{L}^{23}=0$. In both cases we assume natural widths.

In fig. 4 we show the exclusion limit on the $(M_{U},\leavevmode\nobreak\ M_{Z^{\prime}})$ plane for a variety of overall coupling strengths, $g_{U}=g_{Z^{\prime}}\in\{2.5,\,3.0,\,3.5\}$. The increase of the limits with growing coupling in each step is relatively small for the $Z^{\prime}$ ($\sim 200{\,\rm GeV}$), while it is larger for the leptoquark ($\sim 600{\,\rm GeV}$). We see that the decoupling regimes considered in the previous two sections hold when the $Z^{\prime}$ is heavier than (roughly) $3$ TeV, and when the leptoquark is heavier than $5-6$ TeV.

Below the decoupling regime, the limits on both particles strengthen by a few hundred GeV, since they both contribute to the $m_{\text{T}}^{\text{tot}}$ distribution.

We now highlight the importance of including more than just the highest bin in $m_{\text{T}}^{\text{tot}}$ in setting the mass limit. In fig. 5 (left) we plot the $m_{\text{T}}^{\text{tot}}$ distribution of the data and background from [44], along with our simulated leptoquark and $Z^{\prime}$ contributions. We show the distributions for $g_{U}=3$ and $g_{Z^{\prime}}=3$, for masses at the 95% C.L. limit. After a peak, the background steadily falls with increasing $m_{\text{T}}^{\text{tot}}$. The final bin has a larger number of events than the preceding bin as this bin is wider and as it includes overflow events. As such, the final three bins each contain a similar number of background events. Since tau pair production via a $Z^{\prime}$ proceeds through an $s$ channel, it is more peaked in $m_{\text{T}}^{\text{tot}}$ and the events from a multi-TeV $Z^{\prime}$ cluster in the highest energy bin. However, tau pair production via a leptoquark proceeds through $t$ channel process, so there is no clear peak in the invariant mass distribution. This leads the distribution in $m_{\text{T}}^{\text{tot}}$ to extend to lower values. We see in fig. 5 (right) the impact of including only the $N$ highest bins in the CLs calculation. For the $Z^{\prime}$, the limit obtained with only the highest bin is almost $200{\,\rm GeV}$ lower than the limit including all bins. For the leptoquark, when only the highest bin is included, the 95% C.L. limit is around $400{\,\rm GeV}$ weaker than when all bins are included. When the highest two bins are included the difference reduces to around $100{\,\rm GeV}$, and slowly improves as more bins are added. We see that it is crucial to include more than the highest bin in $m_{\text{T}}^{\text{tot}}$ to produce an accurate estimate of the leptoquark exclusion limit. However, it should be noted that we have not been able to account for possible correlations between the bin errors, which could impact the derived exclusion limits.

4.3 $pp\to\tau\nu$ search

The ATLAS and CMS collaborations have performed searches for heavy resonances decaying to $\tau\nu$ (with hadronically decaying $\tau$) using $36.1\leavevmode\nobreak\ \textrm{fb}^{-1}$ [69] and $35.9\leavevmode\nobreak\ \textrm{fb}^{-1}$ [45] of 13 TeV data, respectively. In this section we reinterpret this search in the context of the model in section 3 to set limits on the vector leptoquark mass as a function of $\beta_{23}^{L}$. In our limits we use the CMS data. Since ATLAS data presents a (small) upper fluctuation with respect to the SM background, a combination of ATLAS and CMS data yields slightly weaker limits than CMS data alone (see e.g. [14]).

4.3.1 Search strategy

We compute the NP contribution to the $pp\to\tau_{h}\nu$ process, including the interference with the SM, using Madgraph5_aMC@NLO v2.6.3.2 [59], with the NNPDF23_lo_as_0119_qed PDF set [60]. Hadronization of the $\tau$ final state is done with Pythia 8.2 [61] using the CUETP8M1 set of tuned parameters [70]. The detector response is simulated using Delphes 3.4.1 [63]. The CMS Delphes card has been modified to satisfy the object reconstruction and identification requirements, in particular we include the $\tau$-tagging efficiencies quoted in the experimental search [45].

After showering and detector simulation, we apply the selection cuts specified in table 4 using MadAnalysis 5 v1.6.33 [64]. As a cross-check we have generated the SM Drell-Yan $pp\to\tau\nu$ background and compared our results to the one quoted by CMS. A good agreement is found between the two samples, within $20\%$ of the quoted central values.

After passing through selection cuts, the resulting events are binned according to their total transverse mass,

[TABLE]

with $p_{T}^{\tau_{h}}$ and $p_{T}^{\rm miss}$ being, respectively, the transverse momenta of the visible decay products of the $\tau$ and the missing transverse momentum in the reconstructed event. We compare our binned events with the data and background estimates presented in fig. 3 (left) of [45]. As we did in section 4.2, we consider all the available bins in the $m_{T}^{\rm tot}$ distribution, treating their errors as uncorrelated. For the statistical analysis we use the modified frequentist $\textrm{CL}_{s}$ method [65] computed with the ROOT [66] package Tlimit [67]. In the determination of the limit, we include a systematic uncertainty of $20\%$ in the NP signal to account for possible uncertainties related to the PDF, tau hadronization, detector simulation and unaccounted NLO corrections.

4.3.2 Limits on the $U_{1}$ leptoquark

For this search, we fix $\beta_{L}^{33}=1$ and consider two different benchmarks for the right-handed coupling, $|\beta_{R}^{33}|=0,\,1$. In this case the relative sign between $\beta_{R}^{33}$ and $\beta_{L}^{33}$ is not observable in this channel. Since the leptoquark width plays a marginal role, we fix it to its natural value. We furthermore set $g_{U}=3$.

We compute exclusion limits for the vector leptoquark in the $(\beta_{L}^{23},M_{U})$ plane, see fig. 6. For comparison, we overlay the corresponding limits from $pp\to\tau\tau$ (see fig. 1 (right)) and pair-production limits. As can be seen, these limits give complementary information to those presented in section 4.2.2, offering more stringent limits only when the $\beta_{L}^{23}$ coupling becomes large. Analogous limits for the case $\beta_{R}^{33}=0$ have already been derived in the past literature [14]; we find good agreement between these limits and the ones quoted here. Interestingly, and as happens in the $pp\to\tau\tau$ search, the exclusion bounds get significantly affected by non-zero values of $\beta_{R}^{33}$. The different shapes in the exclusion bands can be understood from the fact that, for $|\beta_{R}^{33}|=1$, the dominant partonic process is $bc\to\tau\nu$, whose cross section scales as $\sigma_{bc\to\tau\nu}\sim|\beta_{L}^{23}|^{2}/M_{U}^{4}$ in the EFT limit. On the contrary, for $\beta_{R}^{33}=0$, the relative contribution from $sc$ production, for which $\sigma_{sc\to\tau\nu}\sim(|\beta_{L}^{23}|/M_{U})^{4}$, is important and even becomes dominant for medium-size values of $\beta_{L}^{23}$.

4.4 $pp\to\tau\mu$ search

The ATLAS collaboration has published a search for heavy particles decaying into different-flavour dilepton pairs using $36.1\leavevmode\nobreak\ \textrm{fb}^{-1}$ [46] of 13 TeV data. In this section we recast the ATLAS data and reinterpret the collider bounds in terms of the model in section 3 to set limits on $\beta_{32}^{L}$ and $\zeta^{23}_{\ell}$, as a function of the leptoquark and $Z^{\prime}$ masses, respectively.

4.4.1 Search strategy

We use Madgraph5_aMC@NLO v2.6.3.2 [59] with the NNPDF23_lo_as_0119_qed PDF set [60] to compute the NP contribution to the $pp\to\tau\mu$ process. The output is passed to Pythia 8.2 [61] for tau hadronization and the detector effects are simulated with Delphes 3.4.1 [63]. The ATLAS Delphes card has been adjusted to satisfy the object reconstruction and identification criteria in the search. In particular we have modified the muon efficiency and momentum resolution to match the High-$\mathit{p_{T}}$ muon operating point, and adjusted the missing energy reconstruction to account for muon effects. We have further included the $\tau$-tagging efficiencies quoted in the experimental search [46].

After showering and detector simulation, we apply the selection cuts specified in table 5 using MadAnalysis 5 v1.6.33 [64]. The resulting events are binned according to their dilepton invariant mass. Following the approach described by ATLAS [46], the tau momentum is reconstructed from the magnitude of the missing energy and the momentum direction of the visible tau decay products. This approach relies on the fact that the momentum of the visible tau decay products and the neutrino momentum are nearly collinear.

In order to validate our procedure, we have simulated the $Z^{\prime}$ signal quoted in the experimental search [46], finding good agreement between our signal and the one by ATLAS.

We compared our results with the binned invariant mass distribution in [46]. Since the error correlations are not provided, we treat the bin errors as uncorrelated. We use the modified frequentist $\textrm{CL}_{s}$ method [65] to obtain $95\%$ CL limits. These limits are computed using the ROOT [66] package Tlimit [67]. In the determination of those limits, we include a systematic uncertainty of $20\%$ for the NP signal to account for possible uncertainties related to the PDF, tau hadronization, detector simulation and unaccounted NLO corrections.

4.4.2 Limits on the $U_{1}$ leptoquark

Following a similar strategy as for the other channels, we fix $g_{U}=3$ and $\beta_{L}^{33}=1$, and take two benchmark values for the right-handed coupling $|\beta_{R}^{33}|=0,1$ (different sign choices for this parameter do not have an impact on the high-$p_{T}$ signal). Varying the leptoquark width only yields a subleading effect so we keep it fixed to its natural value.

We decouple the $Z^{\prime}$ and compute the exclusion limits for the vector leptoquark mass as a function of $\beta_{L}^{32}$, see fig. 7 (left). As in previously analysed channels, the exclusion limits vary significantly for different values of $\beta_{R}^{33}$. We additionally overlay the corresponding exclusion limit obtained from the $pp\to\tau\tau$ search and searches for pair-production. The limits from $pp\to\tau\mu$ become stronger than those obtained from $pp\to\tau\tau$ only for large values of the $\beta_{L}^{32}$ parameter, especially in the case when $|\beta_{R}^{33}|=1$.

The limits presented here offer complementary constraints to the ones obtained from low-energy flavour observables. Indeed, one can establish a one-to-one correspondence between this search and the experimental limits from $\Upsilon(nS)\to\tau\mu$ decays. Using the expression in [24] we get

[TABLE]

This is to be compared to the current experimental limit: $\mathcal{B}(\Upsilon(2S)\to\tau\mu)_{\rm exp}<3.3\cdot 10^{-6}$ ($90\%$ CL). Interestingly, we find the current bounds from high-$p_{T}$ data (see fig. 7 left) to be much more constraining than those from its low-energy counterpart.666Current limits from loop-mediated transitions, such as $\tau\to\mu\gamma$, offer stronger bounds in certain UV completions [26]. However, these bounds are more sensitive to the details of the UV completion and are therefore less robust. Future improvements on $pp\to\tau\mu$ searches can serve as a valuable probe of the leptoquark flavour structure.

4.4.3 Limits on the $Z^{\prime}$

We finally comment on the limits on the $Z^{\prime}$, decoupling the leptoquark. We fix $g_{Z^{\prime}}=3$, $\zeta_{q,u,d}^{33}=\zeta_{\ell,e}^{33}=1$ and set the $Z^{\prime}$ width to its natural value. Limits on the $Z^{\prime}$ mass as a function of $\zeta_{\ell}^{32}$ are shown in fig. 7 (right). As in the leptoquark case, we overlay the corresponding limits on the $Z^{\prime}$ mass extracted from $pp\to\tau\tau$. As can be seen, these limits are always stronger than those from the present search, irrespective of the value of $\zeta_{\ell}^{32}$.

4.5 $pp\to t\bar{t}$ search

We finally turn our attention to searches in the ditop final state, which is subject to NP effects from $s$ channel colorons and $Z^{\prime}$ bosons. We focus our analysis on the coloron since the bounds from this channel on the $Z^{\prime}$ are weaker than the ones reported in section 4.2.3.

4.5.1 Search strategy

We perform a recast of the ATLAS study [47], using 36 fb*-1* of collected data. Since the data was unfolded in this work, we can compute parton-level predictions and directly compare them to the unfolded distributions provided in the reference study.

We choose to derive the constraints from the normalised parton-level differential cross-sections as a function of the $t\bar{t}$-invariant mass, shown in fig. 14(b) of [47]. As in the other searches, we do not include possible error correlations between the bins in the invariant-mass distribution since they are not provided. Our signal predictions are derived by integrating the leading-order SM QCD partonic cross-sections $q\bar{q}\to t\bar{t}$ and $gg\to t\bar{t}$ and the NP contributions from coloron and $Z^{\prime}$ over the parton distribution functions, employing the NNPDF30_nlo_as_0119 PDF set [60] and fixing the factorisation and renormalization scale to the center of the corresponding $t\bar{t}$-invariant mass bin. We use the running strong coupling constant as provided by the PDF set. The only cut applied is on transverse momentum of either top quark: $p_{T}^{t}>500\leavevmode\nobreak\ \mathrm{GeV}$. Note that our reference study places the cuts as $p_{T}^{t,1}>500\leavevmode\nobreak\ \mathrm{GeV}$ on the leading top, and $p_{T}^{t,2}>350\leavevmode\nobreak\ \mathrm{GeV}$ on the subleading one. For a fully exclusive, partonic $t\bar{t}$ final state, $p_{T}^{t,1}=p_{T}^{t,2}$ and hence the second cut does not influence our calculation. However, the unfolded distributions are derived from data which employ this slightly milder cut, leading to slight deviations in bins of lower invariant mass. We therefore drop the bins $m_{t\bar{t}}<1.2\leavevmode\nobreak\ \mathrm{TeV}$ and then find excellent agreement with the SM predictions presented in the ATLAS study. While the analysis also provides unfolded spectra differential in $p_{T}$ and various other kinematic observables, we find the invariant mass spectrum to be the most constraining distribution. We therefore focus solely on the invariant mass spectrum and do not consider searches in the angular spectra.

4.5.2 Limits on the coloron

We are now ready to present the constraints on the various parameters related to the coloron. Throughout this section we set $\kappa_{q}^{33}=1$ and $\kappa_{q}^{ll}=\kappa_{u}^{ll}=\kappa_{d}^{ll}=-(g_{s}/g_{G^{\prime}})^{2}$, and we fix $g_{G^{\prime}}=3$, unless otherwise stated.

Before discussing each plot by itself, a few general comments are in order. Since the search is in a normalised spectrum, it is mostly sensitive to signals that create a shape sufficiently distinct to the background. A distinct shape in this case means a change in the spectrum that is unaffected by normalisation, meaning a peak or a change in the overall slope. Uniform shifts in the spectra (originating from resonances with both very high and very low masses, or large widths777A very wide resonance also leads to a suppression in the overall signal cross-section, further decreasing the constraining power of the search.) are washed out by the overall normalisation. As a result, the strongest bounds are obtained when the coloron can be produced nearly on-shell and the width is moderate. For example, in fig. 8 we show the signal of two parameter points with $\kappa_{q}^{ll}=0$, $M_{G^{\prime}}=2.5$ TeV and different choices of the coloron width. We see that the narrower coloron results in a larger change in the slope.

In fig. 9, we show exclusion regions for the coloron with its natural width and with a width enhanced by a factor of two. In the left panel, the exclusion limits in the $(g_{G^{\prime}},M_{G^{\prime}})$ plane are shown for the natural width and twice this value. An interesting feature of these exclusion regions is that the boundaries bend towards smaller masses for larger couplings. This can be understood by the fact that while the cross section grows with the coupling, so does the width. For the reasons discussed above, the search then loses sensitivity to the resulting signal.

In the right panel of fig. 9, exclusions are shown for varying values of the coupling to left-handed light quarks $\kappa_{q}^{ll}$, keeping the right-handed couplings $\kappa_{u,d}^{ll}$ fixed. With larger couplings to the light quarks, the production of the coloron from valence quarks of the proton increases drastically. Since the valence quarks tend to carry more of the protons’ momenta, they can produce the coloron closer to its mass shell, leading to a signal that the search can more easily discriminate from the background. If we were to set $\kappa_{u}^{ll}=\kappa_{d}^{ll}=\kappa_{q}^{ll}=0$, we would found only very mild bounds, in which case the pair-production search discussed in section 4.1 outperforms this one. If $\kappa_{q}^{ll}$ is chosen to be positive, the bounds tend to be weaker due to interference between the NP and the SM contributions.

Finally, fig. 10 shows exclusion limits with varying widths of the coloron. The different curves (solid, dashed, dotted) show various different choices of relations between the couplings to left- and right-handed light quarks. As expected, limits get weaker with increasing width of the resonance. When the sign of $\kappa_{q}^{ll}$ is chosen to be opposite of $\kappa_{u,d}^{ll}$, the bounds also become weaker for the same reason as discussed above. The grey bands denote the regions in which the floating width parameter is below the partial width to quarks. Note that for $\kappa_{G^{\prime}}\neq 0$ the coloron can decay to two gluons, in which case the actual width would become significantly larger than the partial width to quarks alone.888In this case, the production cross section of the coloron would also be drastically increased, leading to much stronger bounds on its mass from this search.

5 Conclusions

The high-$p_{T}$ phenomenology of models predicting a TeV-scale $SU(2)_{L}$ singlet vector leptoquark which is able to account for the hints of LFU violations observed in $B$-meson decays is quite rich. This is both because this exotic mediator can manifest itself in different final states accessible at the LHC, and also because this state cannot be the only TeV-scale exotic vector. As we have shown, the minimal consistent set of massive vectors comprising a $U_{1}$ also includes a coloron and a $Z^{\prime}$. In this paper we have presented a comprehensive analysis of the high-$p_{T}$ signatures of this set of exotic TeV states, deriving a series of bounds on their masses and couplings.

The results have been discussed in detail in the previous section and will not be repeated here. Here we limit ourself to summarise a few key messages, emphasising the novelties of our analysis compared to the results in the existing literature:

• •

In most of the relevant parameter space the most stringent bound on the leptoquark is obtained by the $pp\to\tau\tau$ process. In this channel a possible $O(1)$ right-handed coupling ($\beta_{R}^{33}$) has a very large impact, as shown in fig. 1.

• •

A non-vanishing off-diagonal coupling of the $U_{1}$ to quarks has a modest impact in $pp\to\tau\tau$, provided |\beta_{L}^{23}|\lower 3.01385pt\hbox{;\stackrel{{\scriptstyle\textstyle<}}{{\sim}};}0.2 (as expected from a natural flavour structure), but a significantly larger impact in $pp\to\tau\nu$. However, the latter search remains subleading compared to $pp\to\tau\tau$ up to |\beta_{L}^{23}|\lower 3.01385pt\hbox{;\stackrel{{\scriptstyle\textstyle<}}{{\sim}};}0.8 for for $|\beta_{R}^{33}|=1$ (or up to |\beta_{L}^{23}|\lower 3.01385pt\hbox{;\stackrel{{\scriptstyle\textstyle<}}{{\sim}};}0.6 for $|\beta_{R}^{33}|=0$).

• •

For large non-vanishing off-diagonal coupling to leptons, a potentially interesting channel is $pp\to\tau\mu$. In the pure left-handed case, the bound from this channel is stronger than the one from $pp\to\tau\tau$ if $|\beta_{L}^{32}|\geq 0.5$ (see fig. 7).

• •

Taking $g_{U}=g_{Z^{\prime}}$ and assuming dominant third-generation coupling to fermions and small couplings to the light families, the constraints on the $Z^{\prime}$ mass are significantly weaker than those on the $U_{1}$ (see fig. 4). The combination of $U_{1}$ and $Z^{\prime}$ signals in $pp\to\tau\tau$ leads to a modest increase on the corresponding bounds, confined to a relatively narrow region of the parameter space.

• •

The bound on the coloron from $pp\to t\bar{t}$ is quite sensitive to the width of this state, and to the possible coupling to light quarks. Due to the increase of the width, the bounds become weaker at large couplings (see fig. 10 (left)).

The bounds we have obtained are very general and can be applied to a large class of models. One of the advantages of having analysed the three states together is the possibility of performing a direct comparison of the bounds obtained from the different mediators (via different processes) on the same model parameter space. As an illustration of this fact, in fig. 11 we show a comparison of $U_{1}$ and coloron bounds in the $(g_{U},M_{U})$ plane, assuming the following relation between their masses and couplings

[TABLE]

This relation follows from the gauge symmetry in eq. 8 assuming two breaking terms transforming as $\omega_{3}\sim(4,\bar{3})$ and $\omega_{1}\sim(4,1)$ under $SU(4)\times SU(3)^{\prime}$ [15, 18]. As can be seen, there is an interesting interplay between the two types of bounds, which changes according to the (model-dependent) ratio $\omega_{1}/\omega_{3}$. Once more, it is worth stressing the importance of the possible right-handed coupling of the $U_{1}$ (neglected in previous analyses): while the coloron sets the most stringent bounds on most of the parameter space for $\beta_{R}^{33}=0$, this is no-longer true for $\beta_{R}^{33}=O(1)$. This fact has relevant phenomenological consequences. For instance, the benchmark point corresponding to $20\%$ correction in $R(D)$ discussed after Eq.(18) is excluded if $\beta_{R}^{33}=0$, while it is not yet excluded for $\beta_{R}^{33}=-1$. More generally, given the minor role of the coloron bounds when $\beta_{R}^{33}=O(1)$, we can state that there is a more direct connection between high-$p_{T}$ physics and $B$-physics anomalies in models with a large right-handed leptoquark coupling.

Acknowldgements

We thank Claudia Cornella and Admir Greljo for useful comments on the manuscript. JFM is grateful to the Mainz Institute for Theoretical Physics (MITP) for its hospitality and its partial support while this work was being finalised. We would like to thank Hubert Spiesberger for allowing us to use the THEP Cluster in Mainz for parts of this study. This research was supported in part by the Swiss National Science Foundation (SNF) under contract 200021-159720.