
2
concrete example is the lepton-specific two Higgs doublet
model (2HDM) [19–21].
Without loss of generality and referring to any specific
model, we write the effective LFV Yukawa couplings of
the (B)SM Higgs boson to the charged leptons as
LY⊃ −Yij ¯
`LieRj h(H)+H.c.(4)
with Yij (i6=j) as free parameters. We set the diago-
nal couplings Yii to their respective SM values, i.e. Yii =
mi/v. The new interactions in Eq. (4) can lead to new
LFV Higgs decay modes that may be directly observable
in current and future collider experiments. The ATLAS
and CMS collaborations have performed several searches
to study the LFV decays of the SM Higgs boson in the
h→eµ [22], eτ [22–26] and µτ [23–28] channels at the
LHC; however, any significant excess over SM expecta-
tions is yet to be observed. LFV decays of neutral heavy
resonances and heavy Higgs bosons at the LHC have also
been investigated in Refs. [29–33] and [34,35], respec-
tively. The prospects of probing LFV signals induced by
Higgs at the future lepton [36–40] and hadron [41–50]
colliders have also been widely explored.
In this work, we focus on the LFV Higgs signal in the
µτ channel that can be effectively probed in a model-
independent way at the HL-LHC using the vector bo-
son fusion (VBF) channel. We first study the pro-
jected sensitivity for LFV decays of the SM-like Higgs
boson, h→µτ, in the VBF Higgs production chan-
nel, pp →(h→µτ)jj, at the HL-LHC (√s= 14 TeV,
L= 3 ab−1). We then perform a detailed collider analy-
sis to explore LFV decays of a BSM Higgs boson Hin the
VBF production channel, pp →(H→µτ )jj for several
BSM Higgs masses mH, and derive ‘model-agnostic’ pro-
jected upper limits on the production cross-section times
Br(H→µτ) for mH∈[20,800] GeV.
Typically, at the hadron colliders, the major impetus
has been on gluon gluon fusion (ggF) Higgs production
mode, while LFV Higgs decays in the VBF channel have
been much less explored. This bias is understandable
since the ggF production rates are much larger than VBF
in a typical SM Higgs-like scenario with mh∼125 GeV.
As expected, the leading sensitivity in searches for LFV
decays of harises from the non-VBF Higgs production
category, largely constituted by ggF production. This is
also highlighted in the recent searches by CMS [26] and
ATLAS [25], where the limits from VBF signal regions
alone are weaker by a factor of few than their non-VBF
counterparts. However, the VBF production channel be-
comes extremely relevant in new physics scenarios with
extended Higgs sectors like the singlet and 2HDM ex-
tensions. The ggF and VBF production rates become
comparable for heavier BSM Higgs states Hat masses
closer to O(1) TeV [51,52]. In addition, the distinct phe-
nomenological features of the VBF topology offer better
control for signal-to-background discrimination than the
ggF signal. Overall, the VBF channel can play a comple-
mentary role, if not leading, in the search for LFV decays
of BSM Higgs extensions and may lead to exciting the-
oretical implications for new physics, as we show in this
work.
It is worth noting that LFV decays of the (B)SM Higgs
boson can also be realized in other Higgs production
channels, such as Higgstrahlung process pp →Zh(H).
For a SM-like Higgs scenario with mh∼125 GeV, the
VBF production cross-section is 4.4 times larger than
the Zh production rate at the √s= 14 TeV LHC [51].
The disparity grows wider at higher Higgs masses; for
instance, at mH∼1 TeV, the VBF to ZH cross-section
ratio is ∼292 [51]. Because of considerably smaller cross-
sections, especially at heavier Higgs masses, the sensitiv-
ity for LFV decays of BSM Higgs bosons in the Hig-
gstrahlung channel is expected to be sub-leading than
in ggF and VBF modes in a typical 2HDM extension.
Therefore, the Zh mode is not considered in the present
analysis.
As for the LFV signal itself, ideally, all three LFV de-
cay channels, h/H →eµ, eτ and µτ , should be considered
in the search for LFV decays of the (B)SM Higgs bosons.
However, the partial decay width of an SM Higgs-like
boson into the LFV final states is typically proportional
to the mass of the heavier lepton, resulting in a usu-
ally suppressed signal rate for the eµ channel than eτ
and µτ. Additionally, the rare µdecay processes typi-
cally impose stringent upper limits on Yeµ, thus further
restricting the search potential in the eµ channel. The
eτ and µτ decay channels result in roughly similar sen-
sitivity at the LHC [26]; however, the background sim-
ulation for the eτ channel is more challenging due to
relatively more significant contamination from the non-
prompt-lepton backgrounds. Moreover, we would also
like to connect the LFV signal at HL-LHC with the preci-
sion low-energy observable of muon anomalous magnetic
moment [53], for which the loop contribution involving a
µτ flavor-violating Higgs is typically enhanced by a factor
of mτ/mµ[12]. In addition, we find that the low-energy
constraint from muon electric dipole moment (EDM) on
Yµτ is 10 orders of magnitude less stringent than that
on Yeτ from electron EDM [54], thus making the collider
study of the µτ channel more relevant. Due to the above
reasons, we only focus on LFV decays of the Higgs boson
in the µτ channel.
The rest of the paper is organized as follows. In Sec. II
we present various low energy LFV constraints on the
Yukawa couplings Yµτ and Yτµ. Sec. III discusses the pro-
jected sensitivity of the LFV couplings of the SM Higgs
boson at the HL-LHC in the VBF channel. Sec. IV dis-
cusses the HL-LHC reach for a generic BSM Higgs LFV
decay, as well as a specific example in 2HDM. We con-
clude in Sec. V.
II. LOW-ENERGY CONSTRAINTS
The LFV couplings in Eq. (4) are subject to various
low energy constraints discussed below.