Experimental studies of superconducting gap structure and quantum fluctuations in novel superconductors and heavy fermion compounds
Kartik Panda

of Physics
Mission Vivekananda Educational and Research Institute


CTPU-PTC-23-38

Current status on pair-produced muon-philic vectorlike leptons in multilepton
channels at the LHC

Junichiro Kawamura${}^{\ a}$<EMAIL_ADDRESS>Seodong Shin${}^{\ a,b}$
<EMAIL_ADDRESS>

a Particle Theory and Cosmology Group, Center for Theoretical Physics of the
Universe, Institute for Basic Science (IBS), Daejeon, 34126, Korea

b Laboratory for Symmetry and Structure of the Universe, Department of
Physics, Jeonbuk National University, Jeonju, Jeonbuk 54896, Korea

In this work, we obtain the current limits on the pair production of
vectorlike leptons decaying to a Standard Model gauge boson and a lepton in
the second generation using the Run-2 data at the LHC. Since there is no
dedicated search, we recast the ATLAS analyses searching for the type-III
seesaw heavy leptons in the multi-lepton channels. There is no limit for the
$SU(2)_{L}$ singlet vectorlike lepton beyond about 100 GeV, while the limit is
about 780 GeV for the doublet one. Thus, dedicated searches for the vectorlike
leptons are necessary, especially for the singlet one. We also study the
general cases of the vectorlike lepton decays and future sensitivities the HL-
LHC.

##  1 Introduction

Vectorlike (VL) fermions are key ingredients in many new physics models beyond
the Standard Model (SM) adopted to resolve both theoretical and experimental
issues. Since chiral fermions in the fourth family are excluded experimentally
[1, 2], these are considered to be vectorlike and their masses are given
independently to the Higgs mechanism in the SM. The VL fermions are introduced
in, for instance, supersymmetric models [3, 4, 5, 6, 7, 8, 9, 10], gauge
mediated supersymmetry breaking scenario [11, 12, 13, 14, 15, 16], composite
Higgs models [17, 18, 19], KSVZ axion models [20, 21], axion-like particle
models [22, 23, 24, 25], alternative solutions of the strong CP problem [26,
27], two Higgs doublet model augmented by VL fermions [28, 29, 30, 31, 32, 33,
34, 35, 36, 37, 38], and models for gauge coupling unification [39, 40].

Among the fourth family fermions, VL leptons (VLLs) with nonzero lepton number
play unique roles in constructing lepton-philic dark matter (DM) models [41,
42, 43, 44, 45, 46], mirror sector models [47, 48, 49], and explanations for
the muon anomalies [50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
64, 65] 111 The latest experimental result [66] confirms the previous results
[67, 68] which might deviate from the SM prediction [69]. . Interestingly, the
lightest VLL is expected to be in the reach of the Large Hadron Collider (LHC)
or High-Luminosity (HL)-LHC. Both ATLAS and CMS collaborations search for the
pair production of VLLs each of which dominantly decays to a SM boson and a
tau lepton [70, 71, 72]. For the doublet VLL 222Throughout this work, doublet
(singlet) for VLL means iso-doublet (iso-singlet) under the $SU(2)_{L}$ gauge
symmetry., the ATLAS search excludes the mass range of
$130<m_{\mathrm{VLL}}<900~{}\mathrm{GeV}$ [72], and the CMS search excludes
the mass up to 1045 GeV [71]. The singlet VLL is less constrained and the
limit is less than $150~{}\mathrm{GeV}$ [71]. Prospects of such VLLs at the
future colliders are discussed in Ref. [73]. The pair productions of the VLLs
decaying to a SM boson and a muon (neutrino) are studied by the ATLAS [74] and
the theorists [75, 76] using the Run-I data. The limit is obtained for $m_{\rm
VLL}\lesssim 500$ GeV when the neutral component of the lightest doublet VLL
dominantly decays to a $W$ boson and a muon [75]. There are also studies for
the VLL produced from cascade decays of extra neutral Higgs bosons [28, 29,
30, 31, 35, 36], and signals from the VLLs decaying to a DM particle [77] or
$Z^{\prime}$ boson [78].

In this paper, we study pair-productions of the VLLs, through the Drell-Yan
process, decaying to the second generation lepton, namely muon-philic VLL.
Such VLL is well motivated to explain the experimental anomalies in the muon
$g-2$ [79, 69, 80] and the semi-leptonic $B$ decays [81]. 333 The LHCb has
recently announced the new result of $R_{K^{(*)}}$ consistent with the SM
expectation [82] which requires efforts such as separate measurements of the
branching ratio at Belle-II [83]. In this work, we obtain the current limits
using the Run-2 data at $\sqrt{s}=13$ TeV by simply recasting the ATLAS
analyses searching for the triplet lepton in the type-III seesaw [84, 85]. We
then estimate expected sensitivities at the HL-LHC using the same channels.

This paper is organized as follows. We briefly explain the VLLs in Sec. 2 and
discuss the analysis strategy in Sec. 3. Our main results are shown in Sec. 4.
Finally, we summarize the paper in Sec. 5.

##  2 Vectorlike lepton scenario

Figure 1:  The NLO production cross sections at $\sqrt{s}=13$ TeV.

In this work, we consider the two kinds of VLLs whose gauge quantum numbers
are given by

$\displaystyle(E,\overline{E})=(\mathbf{1}_{-1},\mathbf{1}_{1}),\quad(L,\overline{L})=(\mathbf{2}_{-1/2},\mathbf{2}_{1/2}),$
(1)

where $\mathbf{r}_{Y}$ is a representation $\mathbf{r}$ under $SU(2)_{L}$ and
the index $Y$ is the hypercharge under $U(1)_{Y}$. Here the electric charge is
given by $Q=T_{3}+Y$, with $T_{3}$ being isospin. The doublet $L$ has two
components denoted by $(L^{0},L^{-})$, where the superscript is the electric
charge. These VLLs can be considered as the fourth generation leptons since
$L$ and $E$ have the same quantum numbers as the SM leptons.

For simplicity, we assume that one of the VLLs is heavy enough so that the
lightest charged VLL is almost either singlet or doublet 444 See, e.g. Ref.
[86] for discussions about the mass mixing induced by the electroweak (EW)
symmetry breaking for $m_{E}\sim m_{L}$. . Under this assumption, the mass of
the neutral VL lepton $L^{0}$ is the same as that of the charged counterpart
$L^{-}$ approximately. We further assume that the VLLs exclusively couple to
the second generation of the SM leptons to avoid lepton flavor violations
strongly constrained by the experiments. In Fig. 1, we show the next-to-
leading order (NLO) cross sections of the VLL pair production, $pp\to
E\overline{E}$ or $L\overline{L}$, at $\sqrt{s}=13$ TeV. The cross sections at
the fixed-order NLO are calculated by using MadGraph5 [87] based on the UFO
model file [88] developed in Ref. [89].

The lightest charged VLL decays to a SM boson and a lepton in the second
generation $E^{-}/L^{-}\to W^{-}\nu_{\mu},Z\mu^{-},h\mu^{-}$ if there is no
new particle lighter than the VLL. The neutral one $L^{0}$ decays as $L^{0}\to
W^{+}\mu^{-}$ 555 The decays $L^{0}\to Z\nu_{\mu}$, $h\nu_{\mu}$ will be open
if we introduce additional SM singlet VLL, but we do not consider these modes
for simplicity. . If the mass of the lightest charged VLL is large enough
compared to the SM boson masses, the branching ratios are given by

$\displaystyle\mathrm{BR}\left({E^{-}}\to{W^{-}\nu_{\mu}}\right):\mathrm{BR}\left({E^{-}}\to{Z\mu^{-}}\right):\mathrm{BR}\left({E^{-}}\to{h\mu^{-}}\right)=$
$\displaystyle\ 2:1:1\,,$ (2)
$\displaystyle\mathrm{BR}\left({L^{-}}\to{W^{-}\nu_{\mu}}\right):\mathrm{BR}\left({L^{-}}\to{Z\mu^{-}}\right):\mathrm{BR}\left({L^{-}}\to{h\mu^{-}}\right)=$
$\displaystyle\ 0:1:1\,.$ (3)

as expected from the Goldstone boson equivalence theorem. Nevertheless, we
shall allow other values of the branching ratios since the above patterns can
be broken when some of the new Yukawa couplings have a certain hierarchical
structure [75, 32].

##  3 Analysis

In this section, we summarize our strategy to recast the ATLAS analyses
searching for the triplet leptons in the type-III seesaw model [84, 85]. In
Ref. [84], the signals are composed of two light leptons, i.e. $e$ or $\mu$,
with at least two jets and large missing energy, where the two jets are from a
$Z/W$ boson decay. Both opposite sign (OS) and same sign (SS) leptons are
searched in the analysis, but in our case the signal will not contribute to
the signal regions (SRs) with the SS leptons because our neutral VLL, $L^{0}$,
is a Dirac but not a Majorana fermion. We also note that the SRs containing
muons will be the most important to identify the muon-philic VLL. In Ref.
[85], the signals with 3 or 4 leptons are studied. Overall, our final states
should be either di-leptons with a hadronic $Z/W$ boson or 3/4 charged leptons
involving muons.

We simulate the signal events by using MadGraph5 [87] based on the UFO model
file generated by FeynRules [90, 91]. The decays of the VLLs are handled by
MadSpin [92], and then the parton-level events are showered and hadronized by
PYTHIA8 [93]. The pair productions of the VLLs are simulated up to one
additional parton using the MLM matching [94] with
$\texttt{xqcut}=m_{\mathrm{VLL}}/10$. The fast detector simulation is
performed by Delphes3.4.2 [95] with the default ATLAS card. Jets are
reconstructed using the anti-$k_{T}$ algorithm [96, 97] with $\Delta R=0.4$.
The same cut conditions are applied for the generated events as the ATLAS
analyses [84, 85], but the object-based missing transverse momentum
significance $S(E_{T}^{\mathrm{miss}})$ is calculated by the approximated one
$S(E_{T}^{\mathrm{miss}}):=E_{T}^{\mathrm{miss}}/\sqrt{H_{T}}$, where
$E_{T}^{\mathrm{miss}}$ is the missing transverse momentum and $H_{T}$ is the
sum of transverse momenta of visible particles [98]. Note that the Higgs
decays are irrelevant for the SRs where the b-tagged jets are vetoed. 666 The
b-veto is not required in the $Q2$ SR where the total charge of the leptons
are $\pm 2$, but the rate is tiny in the SR because the VLLs are Dirac
fermions.

We shall study the pair productions of the VLLs:

$\displaystyle pp\to E^{+}E^{-},L^{+}L^{-},L^{\pm}L^{0},L^{0}L^{0},$ (4)

followed by the decays of the VLLs to a lepton in the second generation and a
$Z$ or $W$ boson. We label the processes as $XY$-$VV^{\prime}$, where $X,Y$
are pair-produced VLLs. In this labeling, $X,Y=F^{\pm},L^{0}$, where
$F^{\pm}=E^{\pm},L^{\pm}$. The kinetic distributions will be common for
$E^{\pm}$ and $L^{\pm}$, while the production cross sections are different.
$V~{}(V^{\prime})$ is a SM EW gauge boson from the decay of $X$ ($Y$). The
processes labeled by $F^{+}F^{-}$-$WW$ will not contribute to the SRs because
the leptons produced from the VLL decays are neutrinos, and hence the numbers
of leptons are not enough to pass the cuts. Altogether, the relevant processes
are $F^{+}F^{-}$-$ZZ$, $F^{+}F^{-}$-$WZ$, $L^{0}F^{-}$-$WW$, $L^{0}F^{-}$-$WZ$
and $L^{0}L^{0}$-$WW$, so we simulate these 5 processes.

We use the statistic variable $q_{\mu}$ defined as [99],

$\displaystyle
q_{\mu}=-2\log\frac{L(\mu,\hat{\hat{b}})}{L(\hat{\mu},\hat{b})},$ (5)

where the likelihood function is given by

$\displaystyle L(\mu,b):=\prod_{i}\frac{\left(\mu
s_{i}+b_{i}\right)^{n_{i}}}{n_{i}!}e^{-(\mu
s_{i}+b_{i})}\times\frac{1}{\sqrt{2\pi\sigma_{i}^{2}}}\mathrm{exp}\left[-\frac{(b_{i}-b_{i}^{0})^{2}}{2\sigma_{i}^{2}}\right].$
(6)

Here, $n_{i}$, $s_{i}$, $b_{i}^{0}$ and $\sigma_{i}$ are respectively the
number of observed, signal, background events and the error of the number of
backgrounds in the bin labeled by $i$. Note that the SRs in the analysis are
exclusive with each other because these have the different number of leptons.
In Eq. (5), $\hat{\hat{b}}$ is the number of backgrounds which maximizes $L$
for a given $\mu$, and $(\hat{\mu},\hat{b})$ are the values of $(\mu,b)$ which
maximizes $L$. We only consider the uncertainties from the backgrounds which
are typically the dominant sources for the errors. The 95% C.L. limit is
defined as [100]

$\displaystyle
CL_{s}:=\left.\frac{1-\Phi\left(\sqrt{q_{\mu}}\right)}{\Phi\left(\sqrt{q_{\mu}^{A}}-\sqrt{q_{\mu}}\right)}\right|_{\mu=1}<0.05,$
(7)

where $\Phi$ is the cumulative distribution function of the normal
distribution. Here, $q^{A}_{\mu}$ is the test statistics in Eq. (5) by
replacing $n_{i}\to b^{0}_{i}$. For the future limits, the significance is
given by [99]

$\displaystyle Z=\sqrt{q_{\mu}|_{\mu=1,n=b}}.$ (8)

The number of signal and background events (its error) are rescaled as
$R_{\mathcal{L}}$ ($\sqrt{R_{\mathcal{L}}}$), where $R_{\mathcal{L}}$ is the
ratio of the integrated luminosity $\mathcal{L}$. The number of observed
events and the estimated backgrounds are shown in the ATLAS papers [84, 85].
The number of signal events is calculated as

$\displaystyle
s_{i}=\mathcal{L}\sum_{P}\sum_{D}\sigma_{P}\times\mathrm{Br}_{D}\times\epsilon^{P,D}_{i},$
(9)

where $P=F^{-}F^{+},L^{0}F^{\pm},L^{0}L^{0}$ labels the production of the VLL
pairs whose cross section is $\sigma_{P}$, and $D=ZZ,WW,WZ$ labels the decays
whose branching fraction is $\mathrm{Br}_{D}$ given by
$\mathrm{BR}\left({X}\to{V\ell_{2}}\right)$ $\times$
$\mathrm{BR}\left({Y}\to{V^{\prime}\ell_{2}}\right)$, where
$\ell_{2}=\mu,\nu_{\mu}$. $\epsilon^{P,D}_{i}$ is the number of events passing
the cut in the bin $i$, per the number of events generated. We obtain the
values of $\epsilon^{P,D}_{i}$ by the simulation.

##  4 Results

Figure 2:  Efficiencies from the processes studied in this work. Here,
$F^{\pm}=E^{\pm},L^{\pm}$ is one of the charged VLLs. In addition to the
efficiencies, we show the $95\%$ C.L. upper bounds on the production cross
sections based on the individual SRs on the bottom-right panel. The gray lines
show the production cross sections of the VLL pairs.

In this section, we explain our analysis results on the signal efficiencies,
the current limits, and the future limits after the full running of the HL-
LHC. Figure 2, except the bottom-right panel, shows the efficiencies of the
processes in the relevant SRs: $e\mu$-OS and $\mu\mu$-OS for the two lepton
final states defined in Table 2 of Ref. [84], and $3\ell$-ZL, $3\ell$-ZLveto,
$3\ell$-JNLow, $4\ell$-Q0 and $4\ell$-Q2, defined in Tables 2 and 3 of Ref.
[85]. Here, the label ZL represents leptonically decaying $Z$ boson, and it is
required to exist in $3\ell$-ZL and $3\ell$-JNLow, while it is vetoed in
$3\ell$-ZLveto. The jet multiplicity is required to be less than one in
$3\ell$-JNLow. For the $4\ell$ SRs, $Q$ is the absolute value of the sum of
the lepton charges in final states and it is required to be zero or two. Note
that the $4\ell$-Q2 SR is possible only when some of the leptons are missing
and hence the corresponding efficiencies are lower than the others. The
efficiencies are smaller than $10^{-3}$ in the other SRs, namely $ee$-OS and
$\ell\ell^{\prime}$-SS ($\ell,\ell^{\prime}=e,\mu$), because the VLLs decay to
the second generation leptons and they preserve the lepton numbers. Note that
the efficiencies in the $2\ell$ SRs decrease for the heavier VLL mass since
the event selection requirement $60<m_{jj}<100$ GeV in Ref. [84] is less
likely to be satisfied. Whereas, those in the SRs with $3$ or $4$ leptons
increase for heavier VLL due to larger lepton transverse momentum. All of the
efficiencies drop for light VLL where $m_{\mathrm{VLL}}\lesssim 200$ GeV, due
to the requirement of two leptons with $p_{T}>40\mathrm{GeV}$. We point out
that the $4\ell$-Q0 SR and $\mu\mu$-OS are the dominant event selection
category for the $F^{+}F^{-}$-$ZZ$ channel where the even numbers of leptons
are predicted. For the channels involving a $W$ boson, the $3\ell$-JNLow SR is
the most important, particularly for heavier VLL with $m_{\rm VLL}\gtrsim 700$
GeV. There are three leptons if a $W$ boson decays leptonically with about 20%
BR, and the cuts for the invariant masses of 2-4 leptons systems and
$H_{T}+E_{T}^{\mathrm{miss}}$ are not applied in the $3\ell$-JNLow. The
tighter cut for $m_{T}(\ell_{1,2})$ than the other $3\ell$-SRs are easily
satisfied for heavier VLLs, and thus the efficiency increases for heavier VLLs
in this SR.

In addition to the efficiencies, we show the $95\%$ C.L. upper bounds on the
production cross sections based on the individual SRs on the bottom-right
panel of Fig. 2. The total cross sections of the charged VLL pair-productions
are shown for comparison. The limits can be obtained by multiplying the factor
from the branching fractions to the SM gauge bosons, i.e. $\mathrm{Br}_{D}$ in
Eq. (9). Note that the limits are from the individual processes, and the
actual limits on the VLLs shown in the following will be obtained by combining
these processes.

Figure 3:  Current $95\%$ C.L. limits on the VLL masses for the singlet VLL
(left), doublet without $L^{0}$ (middle) and doublet with $L^{0}$ (right).
There is no limits on the mass in the white region where the corresponding BRs
are too small. The orange points indicate the predictions in the minimal
scenarios.

We obtain the current limits on the VLLs with the signal efficiencies. Figure
3 shows the 95% C.L. lower limits on the VLL mass in the two dimensional plane
of the branching ratios into $W$ boson and $Z$ boson. The VLL is a $SU(2)_{L}$
singlet (doublet) in the left (middle and right) panel. In the middle panel,
only the charged VLL production is taken into account, which can be applied to
the scenarios with large mass splitting between the neutral and charged VLLs,
$m_{L_{0}}\gg m_{L^{\pm}}$, in a theoretical set-up such as mixing with
additional neutral particles [29]. If the masses of charged and neutral
leptons are sufficiently degenerate, the results on the right panel are
applied.

In the left panel of Fig. 3, we see that the current limit on the singlet VLL
mass reaches about 400 GeV for large BR($E\to Z\mu)\sim 1$ due to the leptonic
decay modes of $Z$ boson. However, there is no limit for BR($E\to
Z\mu)\lesssim 0.5$ due to the small production cross sections of $pp\to EE$.
In particular, no limit beyond about 100 GeV exists in the case of the minimal
singlet VLL model shown in Eq. (2), pointed by the orange marker. On the other
hand, the limit is stronger for the doublet case because of the larger cross
sections. Without the production of the neutral VLL the limit is about $m_{\rm
VLL}\gtrsim 500$ GeV for $\mathrm{BR}\left({L}\to{Z\mu}\right)\sim 1$ and
$m_{\rm VLL}\gtrsim 150$ GeV up to
$\mathrm{BR}\left({L}\to{Z\mu}\right)\gtrsim 0.3$. With the production of the
neutral VLL, we obtain the lower limit of VLL mass at least $620$ GeV,
dominantly determined by $L^{0}L^{0}$-$WW$. The limit reaches to 800 GeV for
$\mathrm{BR}\left({L}\to{Z\mu}\right)\sim 1$, and the limit for the minimal
doublet VLL scenario in Eq. (3) is about 750 GeV. Note that the sensitivities
are mostly determined by the final states belonging to $3\ell$-ZNLow or
$4\ell$-Q0 SRs as seen in Fig. 2, requiring high lepton multiplicity.
Therefore, the contours in Fig. 3 depend more on BR($E/L\to Z\mu$).

Figure 4:  The same figures as Fig. 4, but showing the future exclusion limits
at the HL-LHC.

Figure 4 shows the same figure as Fig. 3 but with the $3~{}\mathrm{ab}^{-1}$
data at the future HL-LHC. Interestingly, we expect the full running of the
HL-LHC can probe a much wider range of the singlet VLL, as shown in the left
panel of the figure. The expected sensitivity on the lower limit of $m_{VLL}$
is about 650 GeV for $\mathrm{BR}\left({E}\to{Z\mu}\right)\sim 1$, and it
reduces to $100~{}\mathrm{GeV}$ for $\mathrm{BR}\left({E}\to{Z\mu}\right)\sim
0.1$. It is remarkable that the minimal scenario in Eq. (3) can be constrained
(discovered) up to about 420 (250) GeV. From the middle panel, without
including the neutral VLL productions, the limit is at most 750 GeV for
$\mathrm{BR}\left({L}\to{Z\mu}\right)\sim 1$ and the limit decreases to 100
GeV as $\mathrm{BR}\left({L}\to{Z\mu}\right)$ decreases up to about
$0.1$-$0.2$ depending on $\mathrm{BR}\left({L}\to{W\nu}\right)$. With the
neutral component productions, from the right panel, the expected limit is
about 1050 GeV for $\mathrm{BR}\left({L}\to{Z\mu}\right)\sim 1$, and decreases
to 880 GeV for $\mathrm{BR}\left({L}\to{Z\mu}\right)\sim 0$.

##  5 Summary

In this work, we studied the LHC searches for the pair productions of the
fourth generation VLLs which decay to a SM EW boson and a lepton in the second
generation. We recast the ATLAS analyses searching for the triplet lepton in
the type-III seesaw [84, 85]. In general, the limits for the VLLs are weaker
than those for the triplets because the production cross sections are smaller
and our VLLs are Dirac fermions not producing the lepton number violating
signals. In the minimal scenario, shown in Eqs. (2) and (3), there is no limit
beyond about $100~{}\mathrm{GeV}$ on the singlet VLL even with the Run-2 data,
while the limit is at about 750 GeV for the doublet case. To obtain the limits
for the singlet case using the Run-2 data, we should consider more dedicated
analysis to overcome the small production cross sections, e.g. by lowering the
threshold for lepton $p_{T}$’s, to be sensitive for the VLL masses within
$100$-$200$ GeV, and searching for a $3\ell$ resonance decayed from a charged
VLL. For the singlet case, it would also be efficient to search for the decays
via the SM Higgs boson which dominantly decays to a pair of bottom quarks.

We also obtained the expected sensitivities at the HL-LHC with
$3~{}\mathrm{ab}^{-1}$ data and found the future limits would reach to 420
(950) GeV for the singlet (doublet) case. It is remarkable that the future
experiment will put the limit on the minimal scenario of the singlet VLL. Note
that our analysis covers the general cases where the branching fractions can
be different from those in Eqs. (2) and (3). The most significant decay mode
is $E^{\pm},L^{\pm}\to Z\mu^{\pm}$, and hence the limits become the most
stringent for $\mathrm{BR}\left({E,L}\to{Z\mu}\right)\sim 1$, and are weaker
for smaller $\mathrm{BR}\left({E,L}\to{Z\mu}\right)$.

The muon-philic VLLs are known to be a good candidate to explain the muon
$g-2$ anomaly. In such a case, the chiral enhancement of the muon $g-2$ by the
muon-philic VLLs is correlated to the rate for $h\to\mu^{+}\mu^{-}$, and the
VLLs lighter than about 500 GeV may be excluded [86, 60, 61]. Our searches can
constrain the scenarios of explaining the muon $g-2$ by the chiral enhancement
of light VLL correlated with the rate for $h\to\mu^{+}\mu^{-}$ as well. The
comprehensive study for the VLL explanation of the muon $g-2$ is our future
work.

## Acknowledgments

The authors thank Radovan Dermisek for useful discussion. This work was
supported by IBS under the project code, IBS-R018-D1. S.S. acknowledges
support from the National Research Foundation of Korea (NRF-2020R1I1A3072747
and NRF-2022R1A4A5030362).

## References

  * [1] A. Djouadi and A. Lenz, Sealing the fate of a fourth generation of fermions, Phys. Lett. B 715 (2012) 310–314, [arXiv:1204.1252].
  * [2] O. Eberhardt, G. Herbert, H. Lacker, A. Lenz, A. Menzel, U. Nierste, and M. Wiebusch, Impact of a Higgs boson at a mass of 126 GeV on the standard model with three and four fermion generations, Phys. Rev. Lett. 109 (2012) 241802, [arXiv:1209.1101].
  * [3] S. P. Martin, Extra vector-like matter and the lightest Higgs scalar boson mass in low-energy supersymmetry, Phys. Rev. D 81 (2010) 035004, [arXiv:0910.2732].
  * [4] P. W. Graham, A. Ismail, S. Rajendran, and P. Saraswat, A Little Solution to the Little Hierarchy Problem: A Vector-like Generation, Phys. Rev. D 81 (2010) 055016, [arXiv:0910.3020].
  * [5] T. Moroi, R. Sato, and T. T. Yanagida, Extra Matters Decree the Relatively Heavy Higgs of Mass about 125 GeV in the Supersymmetric Model, Phys. Lett. B 709 (2012) 218–221, [arXiv:1112.3142].
  * [6] M. Endo, K. Hamaguchi, S. Iwamoto, and N. Yokozaki, Higgs mass, muon g-2, and LHC prospects in gauge mediation models with vector-like matters, Phys. Rev. D 85 (2012) 095012, [arXiv:1112.5653].
  * [7] M. Endo, K. Hamaguchi, S. Iwamoto, and N. Yokozaki, Higgs Mass and Muon Anomalous Magnetic Moment in Supersymmetric Models with Vector-Like Matters, Phys. Rev. D 84 (2011) 075017, [arXiv:1108.3071].
  * [8] R. Dermisek and N. McGinnis, Mass scale of vectorlike matter and superpartners from IR fixed point predictions of gauge and top Yukawa couplings, Phys. Rev. D 97 (2018), no. 5 055009, [arXiv:1712.03527].
  * [9] R. Dermíšek and N. McGinnis, Top-bottom-tau Yukawa coupling unification in the MSSM plus one vectorlike family and fermion masses as IR fixed points, Phys. Rev. D 99 (2019), no. 3 035033, [arXiv:1810.12474].
  * [10] J. Y. Araz, S. Banerjee, M. Frank, B. Fuks, and A. Goudelis, Dark matter and collider signals in an MSSM extension with vector-like multiplets, Phys. Rev. D 98 (2018), no. 11 115009, [arXiv:1810.07224].
  * [11] M. Dine and W. Fischler, A Phenomenological Model of Particle Physics Based on Supersymmetry, Phys. Lett. B 110 (1982) 227–231.
  * [12] L. Alvarez-Gaume, M. Claudson, and M. B. Wise, Low-Energy Supersymmetry, Nucl. Phys. B 207 (1982) 96.
  * [13] C. R. Nappi and B. A. Ovrut, Supersymmetric Extension of the SU(3) x SU(2) x U(1) Model, Phys. Lett. B 113 (1982) 175–179.
  * [14] M. Dine and A. E. Nelson, Dynamical supersymmetry breaking at low-energies, Phys. Rev. D 48 (1993) 1277–1287, [hep-ph/9303230].
  * [15] M. Dine, A. E. Nelson, and Y. Shirman, Low-energy dynamical supersymmetry breaking simplified, Phys. Rev. D 51 (1995) 1362–1370, [hep-ph/9408384].
  * [16] M. Dine, A. E. Nelson, Y. Nir, and Y. Shirman, New tools for low-energy dynamical supersymmetry breaking, Phys. Rev. D 53 (1996) 2658–2669, [hep-ph/9507378].
  * [17] G. Panico and A. Wulzer, The Composite Nambu-Goldstone Higgs, vol. 913.  Springer, 2016.
  * [18] G. Cacciapaglia, C. Pica, and F. Sannino, Fundamental Composite Dynamics: A Review, Phys. Rept. 877 (2020) 1–70, [arXiv:2002.04914].
  * [19] G. Cacciapaglia, A. Deandrea, and K. Sridhar, Review of fundamental composite dynamics, Eur. Phys. J. ST 231 (2022), no. 7 1221–1222.
  * [20] J. E. Kim, Weak Interaction Singlet and Strong CP Invariance, Phys. Rev. Lett. 43 (1979) 103.
  * [21] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Can Confinement Ensure Natural CP Invariance of Strong Interactions?, Nucl. Phys. B 166 (1980) 493–506.
  * [22] J. E. Kim and G. Carosi, Axions and the Strong CP Problem, Rev. Mod. Phys. 82 (2010) 557–602, [arXiv:0807.3125]. [Erratum: Rev.Mod.Phys. 91, 049902 (2019)].
  * [23] D. J. E. Marsh, Axion Cosmology, Phys. Rept. 643 (2016) 1–79, [arXiv:1510.07633].
  * [24] L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli, The landscape of QCD axion models, Phys. Rept. 870 (2020) 1–117, [arXiv:2003.01100].
  * [25] K. Choi, S. H. Im, and C. Sub Shin, Recent Progress in the Physics of Axions and Axion-Like Particles, Ann. Rev. Nucl. Part. Sci. 71 (2021) 225–252, [arXiv:2012.05029].
  * [26] R. N. Mohapatra and G. Senjanovic, Natural Suppression of Strong p and t Noninvariance, Phys. Lett. B 79 (1978) 283–286.
  * [27] K. S. Babu and R. N. Mohapatra, A Solution to the Strong CP Problem Without an Axion, Phys. Rev. D 41 (1990) 1286.
  * [28] R. Dermisek, E. Lunghi, and S. Shin, New decay modes of heavy Higgs bosons in a two Higgs doublet model with vectorlike leptons, JHEP 05 (2016) 148, [arXiv:1512.07837].
  * [29] R. Dermisek, E. Lunghi, and S. Shin, Two Higgs doublet model with vectorlike leptons and contributions to $pp\to WW$ and $H\to WW$, JHEP 02 (2016) 119, [arXiv:1509.04292].
  * [30] R. Dermisek, E. Lunghi, and S. Shin, Contributions of flavor violating couplings of a Higgs boson to $pp\to WW$, JHEP 08 (2015) 126, [arXiv:1503.08829].
  * [31] R. Dermisek, E. Lunghi, and S. Shin, New constraints and discovery potential for Higgs to Higgs cascade decays through vectorlike leptons, JHEP 10 (2016) 081, [arXiv:1608.00662].
  * [32] R. Dermíšek, E. Lunghi, and S. Shin, Hunting for Vectorlike Quarks, JHEP 04 (2019) 019, [arXiv:1901.03709]. [Erratum: JHEP 10, 058 (2020)].
  * [33] R. Dermisek, E. Lunghi, and S. Shin, Cascade decays of heavy Higgs bosons through vectorlike quarks in two Higgs doublet models, JHEP 03 (2020) 029, [arXiv:1907.07188].
  * [34] R. Dermisek, E. Lunghi, N. McGinnis, and S. Shin, Signals with six bottom quarks for charged and neutral Higgs bosons, JHEP 07 (2020) 241, [arXiv:2005.07222].
  * [35] R. Dermisek, J. Kawamura, E. Lunghi, N. McGinnis, and S. Shin, Combined signatures of heavy Higgses and vectorlike fermions at the HL-LHC, in Snowmass 2021, 3, 2022.  arXiv:2203.03852.
  * [36] R. Dermisek, J. Kawamura, E. Lunghi, N. McGinnis, and S. Shin, Leptonic cascade decays of a heavy Higgs boson through vectorlike leptons at the LHC, JHEP 10 (2022) 138, [arXiv:2204.13272].
  * [37] R. Dermisek, K. Hermanek, N. McGinnis, and S. Yoon, Ellipse of Muon Dipole Moments, Phys. Rev. Lett. 129 (2022), no. 22 221801, [arXiv:2205.14243].
  * [38] R. Dermisek, K. Hermanek, N. McGinnis, and S. Yoon, Predictions for Muon Electric and Magnetic Dipole Moments from $h\rightarrow\mu^{+}\mu^{-}$ in Two-Higgs-Doublet Models with New Leptons, arXiv:2306.13212.
  * [39] R. Dermisek, Insensitive Unification of Gauge Couplings, Phys. Lett. B 713 (2012) 469–472, [arXiv:1204.6533].
  * [40] R. Dermisek, Unification of gauge couplings in the standard model with extra vectorlike families, Phys. Rev. D 87 (2013), no. 5 055008, [arXiv:1212.3035].
  * [41] Y. Bai and J. Berger, Lepton Portal Dark Matter, JHEP 08 (2014) 153, [arXiv:1402.6696].
  * [42] S. Chang, R. Edezhath, J. Hutchinson, and M. Luty, Leptophilic Effective WIMPs, Phys. Rev. D 90 (2014), no. 1 015011, [arXiv:1402.7358].
  * [43] J. Kawamura, S. Okawa, and Y. Omura, Interplay between the b$\to s\ell\ell$ anomalies and dark matter physics, Phys. Rev. D 96 (2017), no. 7 075041, [arXiv:1706.04344].
  * [44] L. Calibbi, R. Ziegler, and J. Zupan, Minimal models for dark matter and the muon g$-$2 anomaly, JHEP 07 (2018) 046, [arXiv:1804.00009].
  * [45] J. Kawamura, S. Okawa, and Y. Omura, Current status and muon $g-2$ in explanation of lepton portal dark matter, JHEP 08 (2020) 042, [arXiv:2002.12534].
  * [46] J. Kawamura, S. Okawa, and Y. Omura, W boson mass and muon g-2 in a lepton portal dark matter model, Phys. Rev. D 106 (2022), no. 1 015005, [arXiv:2204.07022].
  * [47] J. Kawamura, S. Okawa, Y. Omura, and Y. Tang, WIMP dark matter in the parity solution to the strong CP problem, JHEP 04 (2019) 162, [arXiv:1812.07004].
  * [48] D. Dunsky, L. J. Hall, and K. Harigaya, Higgs Parity, Strong CP, and Dark Matter, JHEP 07 (2019) 016, [arXiv:1902.07726].
  * [49] T. H. Jung, Spontaneous Twin Symmetry Breaking, Phys. Rev. D 100 (2019), no. 11 115012, [arXiv:1902.10978].
  * [50] K. Kannike, M. Raidal, D. M. Straub, and A. Strumia, Anthropic solution to the magnetic muon anomaly: the charged see-saw, JHEP 02 (2012) 106, [arXiv:1111.2551]. [Erratum: JHEP 10, 136 (2012)].
  * [51] R. Dermisek, A. Raval, and S. Shin, Effects of vectorlike leptons on $h\to 4\ell$ and the connection to the muon g-2 anomaly, Phys. Rev. D 90 (2014), no. 3 034023, [arXiv:1406.7018].
  * [52] Z. Poh and S. Raby, Vectorlike leptons: Muon g-2 anomaly, lepton flavor violation, Higgs boson decays, and lepton nonuniversality, Phys. Rev. D 96 (2017), no. 1 015032, [arXiv:1705.07007].
  * [53] S. Raby and A. Trautner, Vectorlike chiral fourth family to explain muon anomalies, Phys. Rev. D 97 (2018), no. 9 095006, [arXiv:1712.09360].
  * [54] P. Arnan, A. Crivellin, M. Fedele, and F. Mescia, Generic Loop Effects of New Scalars and Fermions in $b\to s\ell^{+}\ell^{-}$, $(g-2)_{\mu}$ and a Vector-like $4^{\rm th}$ Generation, JHEP 06 (2019) 118, [arXiv:1904.05890].
  * [55] J. Kawamura, S. Raby, and A. Trautner, Complete vectorlike fourth family and new U(1)’ for muon anomalies, Phys. Rev. D 100 (2019), no. 5 055030, [arXiv:1906.11297].
  * [56] J. Kawamura, S. Raby, and A. Trautner, Complete vectorlike fourth family with U(1)’ : A global analysis, Phys. Rev. D 101 (2020), no. 3 035026, [arXiv:1911.11075].
  * [57] E. J. Chun and T. Mondal, Explaining $g-2$ anomalies in two Higgs doublet model with vector-like leptons, JHEP 11 (2020) 077, [arXiv:2009.08314].
  * [58] M. Frank and I. Saha, Muon anomalous magnetic moment in two-Higgs-doublet models with vectorlike leptons, Phys. Rev. D 102 (2020), no. 11 115034, [arXiv:2008.11909].
  * [59] A. Crivellin, F. Kirk, C. A. Manzari, and M. Montull, Global Electroweak Fit and Vector-Like Leptons in Light of the Cabibbo Angle Anomaly, JHEP 12 (2020) 166, [arXiv:2008.01113].
  * [60] R. Dermisek, K. Hermanek, and N. McGinnis, Highly Enhanced Contributions of Heavy Higgs Bosons and New Leptons to Muon $g$$-$2 and Prospects at Future Colliders, Phys. Rev. Lett. 126 (2021), no. 19 191801, [arXiv:2011.11812].
  * [61] R. Dermisek, K. Hermanek, and N. McGinnis, Muon g-2 in two-Higgs-doublet models with vectorlike leptons, Phys. Rev. D 104 (2021), no. 5 055033, [arXiv:2103.05645].
  * [62] H. M. Lee, J. Song, and K. Yamashita, Seesaw lepton masses and muon $g-2$ from heavy vector-like leptons, J. Korean Phys. Soc. 79 (2021), no. 12 1121–1134, [arXiv:2110.09942].
  * [63] H. M. Lee and K. Yamashita, A model of vector-like leptons for the muon $g-2$ and the W boson mass, Eur. Phys. J. C 82 (2022), no. 8 661, [arXiv:2204.05024].
  * [64] J. Kawamura and S. Raby, W mass in a model with vectorlike leptons and U(1)’, Phys. Rev. D 106 (2022), no. 3 035009, [arXiv:2205.10480].
  * [65] W. Abdallah, M. Ashry, J. Kawamura, and A. Moursy, Semi-visible dark photon in a model with vector-like leptons for the $(g-2)_{e,\mu}$ and $W$-boson mass anomalies, arXiv:2308.05691.
  * [66] Muon g-2 Collaboration, D. P. Aguillard et al., Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm, arXiv:2308.06230.
  * [67] Muon g-2 Collaboration, G. W. Bennett et al., Final Report of the Muon E821 Anomalous Magnetic Moment Measurement at BNL, Phys. Rev. D73 (2006) 072003, [hep-ex/0602035].
  * [68] Muon g-2 Collaboration, B. Abi et al., Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm, Phys. Rev. Lett. 126 (2021) 141801, [arXiv:2104.03281].
  * [69] T. Aoyama et al., The anomalous magnetic moment of the muon in the Standard Model, Phys. Rept. 887 (2020) 1–166, [arXiv:2006.04822].
  * [70] CMS Collaboration, A. M. Sirunyan et al., Search for vector-like leptons in multilepton final states in proton-proton collisions at $\sqrt{s}$ = 13 TeV, Phys. Rev. D 100 (2019), no. 5 052003, [arXiv:1905.10853].
  * [71] CMS Collaboration, A. Tumasyan et al., Inclusive nonresonant multilepton probes of new phenomena at $\sqrt{s}$=13 TeV, Phys. Rev. D 105 (2022), no. 11 112007, [arXiv:2202.08676].
  * [72] ATLAS Collaboration, Search for third-generation vector-like leptons in $pp$ collisions at $\sqrt{s}=13\,\text{TeV}$ with the ATLAS detector, arXiv:2303.05441.
  * [73] P. N. Bhattiprolu and S. P. Martin, Prospects for vectorlike leptons at future proton-proton colliders, Phys. Rev. D 100 (2019), no. 1 015033, [arXiv:1905.00498].
  * [74] ATLAS Collaboration, G. Aad et al., Search for heavy lepton resonances decaying to a $Z$ boson and a lepton in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector, JHEP 09 (2015) 108, [arXiv:1506.01291].
  * [75] R. Dermisek, J. P. Hall, E. Lunghi, and S. Shin, Limits on Vectorlike Leptons from Searches for Anomalous Production of Multi-Lepton Events, JHEP 12 (2014) 013, [arXiv:1408.3123].
  * [76] N. Kumar and S. P. Martin, Vectorlike Leptons at the Large Hadron Collider, Phys. Rev. D 92 (2015), no. 11 115018, [arXiv:1510.03456].
  * [77] G. Guedes and J. Santiago, New leptons with exotic decays: collider limits and dark matter complementarity, JHEP 01 (2022) 111, [arXiv:2107.03429].
  * [78] J. Kawamura and S. Raby, Signal of four muons or more from a vector-like lepton decaying to a muon-philic Z’ boson at the LHC, Phys. Rev. D 104 (2021), no. 3 035007, [arXiv:2104.04461].
  * [79] Muon g-2 Collaboration, G. W. Bennett et al., Measurement of the negative muon anomalous magnetic moment to 0.7 ppm, Phys. Rev. Lett. 92 (2004) 161802, [hep-ex/0401008].
  * [80] Muon g-2 Collaboration, B. Abi et al., Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm, Phys. Rev. Lett. 126 (2021), no. 14 141801, [arXiv:2104.03281].
  * [81] LHCb Collaboration, R. Aaij et al., Test of lepton universality with $B^{0}\rightarrow K^{*0}\ell^{+}\ell^{-}$ decays, JHEP 08 (2017) 055, [arXiv:1705.05802].
  * [82] LHCb Collaboration, Measurement of lepton universality parameters in $B^{+}\to K^{+}\ell^{+}\ell^{-}$ and $B^{0}\to K^{*0}\ell^{+}\ell^{-}$ decays, arXiv:2212.09153.
  * [83] Belle-II Collaboration, F. Abudinén et al., Measurement of the branching fraction for the decay $B\to K^{\ast}(892)\ell^{+}\ell^{-}$ at Belle II, arXiv:2206.05946.
  * [84] ATLAS Collaboration, G. Aad et al., Search for type-III seesaw heavy leptons in dilepton final states in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector, Eur. Phys. J. C 81 (2021), no. 3 218, [arXiv:2008.07949].
  * [85] ATLAS Collaboration, G. Aad et al., Search for type-III seesaw heavy leptons in leptonic final states in pp collisions at $\sqrt{s}=13~{}\text{TeV}$ with the ATLAS detector, Eur. Phys. J. C 82 (2022), no. 11 988, [arXiv:2202.02039].
  * [86] R. Dermisek and A. Raval, Explanation of the Muon g-2 Anomaly with Vectorlike Leptons and its Implications for Higgs Decays, Phys. Rev. D 88 (2013) 013017, [arXiv:1305.3522].
  * [87] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro, The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [arXiv:1405.0301].
  * [88] C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, and T. Reiter, UFO - The Universal FeynRules Output, Comput. Phys. Commun. 183 (2012) 1201–1214, [arXiv:1108.2040].
  * [89] A. H. Ajjath, B. Fuks, H. S. Shao, and Y. Simon, Precision predictions for exotic lepton production at the Large Hadron Collider, arXiv:2301.03640.
  * [90] A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, FeynRules 2.0 - A complete toolbox for tree-level phenomenology, Comput. Phys. Commun. 185 (2014) 2250–2300, [arXiv:1310.1921].
  * [91] N. D. Christensen and C. Duhr, FeynRules - Feynman rules made easy, Comput. Phys. Commun. 180 (2009) 1614–1641, [arXiv:0806.4194].
  * [92] P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk, Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations, JHEP 03 (2013) 015, [arXiv:1212.3460].
  * [93] T. Sjostrand, S. Mrenna, and P. Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput. Phys. Commun. 178 (2008) 852–867, [arXiv:0710.3820].
  * [94] F. Caravaglios, M. L. Mangano, M. Moretti, and R. Pittau, A New approach to multijet calculations in hadron collisions, Nucl. Phys. B 539 (1999) 215–232, [hep-ph/9807570].
  * [95] DELPHES 3 Collaboration, J. de Favereau, C. Delaere, P. Demin, A. Giammanco, V. Lemaître, A. Mertens, and M. Selvaggi, DELPHES 3, A modular framework for fast simulation of a generic collider experiment, JHEP 02 (2014) 057, [arXiv:1307.6346].
  * [96] M. Cacciari, G. P. Salam, and G. Soyez, The anti-$k_{t}$ jet clustering algorithm, JHEP 04 (2008) 063, [arXiv:0802.1189].
  * [97] M. Cacciari, G. P. Salam, and G. Soyez, FastJet User Manual, Eur. Phys. J. C 72 (2012) 1896, [arXiv:1111.6097].
  * [98] ATLAS Collaboration, Object-based missing transverse momentum significance in the ATLAS detector, tech. rep., CERN, Geneva, 2018.  All figures including auxiliary figures are available at https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2018-038.
  * [99] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymptotic formulae for likelihood-based tests of new physics, Eur. Phys. J. C 71 (2011) 1554, [arXiv:1007.1727]. [Erratum: Eur.Phys.J.C 73, 2501 (2013)].
  * [100] ATLAS, CMS, LHC Higgs Combination Group Collaboration, Procedure for the LHC Higgs boson search combination in Summer 2011, .