The signatures of the new particles $h_2$ and $Z_{\mu\tau}$ at e-p colliders in the $U(1)_{L_\mu-L_\tau}$ model

TL;DR

This paper explores the potential to detect new particles predicted by the ${U(1)}_{L_-L_}$ model at future electron-proton colliders, focusing on their production channels and observability under experimental constraints.

Contribution

It provides a detailed analysis of the production and detection prospects of the scalar $h_2$ and gauge boson $Z_{}$ at e-p colliders, considering beam polarization and parameter constraints.

Findings

Signal significance above 5$$ can be achieved for $h_2$ and $Z_{}$ detection.

$Z_{}$ detection sensitivity reaches 5$$ at certain parameter values.

Detection of $h_2$ decays into SM particles is challenging.

Abstract

Considering the superior performances of the future e-p colliders, LHeC and FCC-eh, we discuss the feasibility of detecting the extra neutral scalar $h_{2}$ and the light gauge boson $Z^{}_{\mu\tau}$, which are predicted by the ${U(1)}_{L^{}_{\mu} - L^{}_{\tau}}$ model. Taking into account the experimental constraints on the relevant free parameters, we consider all possible production channels of $h_{2}$ and $Z^{}_{\mu\tau}$ at e-p colliders and further investigate their observability through the optimal channels in the case of the beam polarization P($e^{-}$)= -0.8. We find that the signal significance above 5$\sigma$ of $h_{2}$ as well as $Z^{}_{\mu\tau}$ detecting can be achieved via $e^{-}p\to{e^{-}jh_{2}(\to{Z_{\mu\tau}Z_{\mu\tau}})}\to~e^{-}j+/\!\!\!\!{E}^{}_{T}$ process and a 5$\sigma$ sensitivity of $Z^{}_{\mu\tau}$ detecting can be gained via $e^{-}p\to{e^{-}jh_{1}(\to{Z^{}_{\mu\tau}Z^{}_{\mu\tau}})\to}~e^{-}j+/\!\!\!\!{E}^{}_{T}$ process at e-p colliders with appropriate parameter values and a designed integrated luminosity. However, the signals of $h_{2}$ decays into pair of SM particles are difficult to be detected.

Figures

Tables

Table 1. Table 1: Particle contents and corresponding charge assignments under various symmetry groups in the U ​ ( 1 ) L μ − L τ 𝑈 subscript 1 subscript 𝐿 𝜇 subscript 𝐿 𝜏 {U(1)}_{L_{\mu}-L_{\tau}} model.

Gauge Group	Scalar Fields			Lepton Fields
	${\phi }_h$	${\phi }_H$	${\phi }_{DM}$	$L_e$	$L_{\mu }$	$L_{\tau }$	$e_R$	${\mu }_R$	${\tau }_R$	$N_R^e$	$N_R^{\mu }$	$N_R^{\tau }$
$SU{\left(2\right)}_L$	2	1	1	2	2	2	1	1	1	1	1	1
$U{\left(1\right)}_Y$	1/2	0	0	-1/2	-1/2	-1/2	-1	-1	-1	0	0	0
$U{\left(1\right)}_{L_{\mu }-L_{\tau }}$	0	1	$n_{\mu \tau }$	0	1	-1	0	1	-1	0	1	-1

Gauge Group	Baryon Fields
	$u_L$	$d_L$	$c_L$	$s_L$	$t_L$	$b_L$	$u_R$	$d_R$	$c_R$	$s_R$	$t_R$	$b_R$
$SU{\left(2\right)}_L$	2	2	2	2	2	2	1	1	1	1	1	1
$U{\left(1\right)}_Y$	1/6	1/6	1/6	1/6	1/6	1/6	2/3	-1/3	2/3	-1/3	2/3	-1/3
$U{\left(1\right)}_{L_{\mu }-L_{\tau }}$	0	0	0	0	0	0	0	0	0	0	0	0

Table 3. Table 2: Effect of individual kinematical cuts on the signal-1 for M h 2 = 300 subscript 𝑀 subscript ℎ 2 300 {M}_{h_{2}}=300 GeV, M Z μ ​ τ = 0.1 subscript 𝑀 subscript 𝑍 𝜇 𝜏 0.1 {M}_{Z_{\mu\tau}}=0.1 GeV, sin α = 0.2 𝛼 0.2 \alpha=0.2 and g μ ​ τ = 1 × 10 − 3 subscript 𝑔 𝜇 𝜏 1 superscript 10 3 g_{\mu\tau}=1\times 10^{-3} and the SM backgrounds at the LHeC with E e − = 140 ​ ( 60 ) subscript 𝐸 superscript 𝑒 140 60 E_{e^{-}}=140~{}(60)  GeV. The statistical significance is computed for the integrated luminosity as 1 ab -1 .

LHeC,  $E_{e^{-}}=140\left(60\right)$ GeV,  $E_p=7$ TeV
cuts	signal (S)	total background (B)	$S/\sqrt{S+B}$
initial (no cut)	417.0 (104.0)	$1.01\times {10}^6$ ( $5.09\times {10}^5$)	0.41 (0.15)
basic cuts	392.1 (97.6)	$8.53\times {10}^5$ ( $4.11\times {10}^5$)	0.42 (0.16)
() GeV	238.6 (70.7)	$5.17\times {10}^4$ ( $4.89\times {10}^4$)	1.05 (0.32)
$\mathrm{TET}>300$ ($\mathrm{TET}>260$) GeV	235.9 (64.6)	$1.65\times {10}^4$ ( $4.63\times {10}^3$)	1.83 (0.94)
$\mathrm{THT}>200$ ($\mathrm{THT}>200$) GeV	235.7 (63.8)	$9.91\times {10}^3$ ( $1.42\times {10}^3$)	2.30 (1.65)
$E_T\left(jj\right)>100$ ($E_T\left(jj\right)>90$) GeV	232.7 (60.4)	$7.83\times {10}^3$ ( $1.08\times {10}^3$)	2.59 (1.79)
$E_T\left(e^{-}j\right)>150$ ($E_T\left(e^{-}j\right)>120$) GeV	232.3 (59.5)	$7.44\times {10}^3$ ( $9.93\times {10}^2$)	2.65 (1.84)

Table 4. Table 3: Effect of individual kinematical cuts on the signal-2 for M Z μ ​ τ = 0.1 subscript 𝑀 subscript 𝑍 𝜇 𝜏 0.1 {M}_{Z_{\mu\tau}}=0.1 GeV and χ = 9 × 10 − 5 ​ GeV − 1 𝜒 9 superscript 10 5 superscript GeV 1 \chi=9\times 10^{-5}~{}\rm{GeV^{-1}} and the SM backgrounds at the LHeC with E e − = 140 ​ ( 60 ) subscript 𝐸 superscript 𝑒 140 60 E_{e^{-}}=140~{}(60)  GeV. The statistical significance is computed for the integrated luminosity as 1 ab -1 .

LHeC,  $E_{e^{-}}=140\left(60\right)$ GeV,  $E_p=7$ TeV
cuts	signal (S)	total background (B)	$S/\sqrt{S+B}$
initial (no cut)	1288.0 (544.0)	$1.01\times {10}^6$ ( $5.09\times {10}^5$)	1.28 (0.76)
basic cuts	1205.9 (508.1)	$8.53\times {10}^5$ ( $4.11\times {10}^5$)	1.30 (0.79)
() GeV	980.9 (386.0)	$8.79\times {10}^4$ ( $4.05\times {10}^4$)	3.29 (1.91)
$\mathrm{TET}>200$ ($\mathrm{TET}>160$) GeV	827.7 (361.9)	$2.26\times {10}^4$ ( $1.27\times {10}^4$)	5.41 (3.16)
$\mathrm{THT}>140$ ($\mathrm{THT}>120$) GeV	808.1 (354.3)	$1.20\times {10}^4$ ( $6.09\times {10}^3$)	7.14 (4.41)
$E_T\left(jj\right)>60$ ($E_T\left(jj\right)>60$) GeV	803.5 (374.8)	$9.65\times {10}^3$ ( $4.46\times {10}^3$)	7.86 (5.02)
$E_T\left(e^{-}j\right)>100$ ($E_T\left(e^{-}j\right)>80$) GeV	796.1 (343.7)	$8.96\times {10}^3$ ( $4.19\times {10}^3$)	8.06 (5.10)

Table 5. Table 4: Effect of individual kinematical cuts on the signal-1 for M h 2 = 300 subscript 𝑀 subscript ℎ 2 300 {M}_{h_{2}}=300 GeV, M Z μ ​ τ = 0.1 subscript 𝑀 subscript 𝑍 𝜇 𝜏 0.1 {M}_{Z_{\mu\tau}}=0.1 GeV, sin α = 0.2 𝛼 0.2 \alpha=0.2 and g μ ​ τ = 1 × 10 − 3 subscript 𝑔 𝜇 𝜏 1 superscript 10 3 g_{\mu\tau}=1\times 10^{-3} and backgrounds at the FCC-eh. The statistical significance S ​ S 𝑆 𝑆 SS is computed for an integrated luminosity of 1 ab -1 .

FCC-eh,  $E_{e^{-}}=60$ GeV,  $E_p=50$ TeV
cuts	signal (S)	total background (B)	$S/\sqrt{S+B}$
initial (no cut)	925.0	$1.81\times {10}^6$	0.69
basic cuts	863.3	$1.30\times {10}^6$	0.76
GeV	303.2	$1.29\times {10}^5$	0.84
$\mathrm{TET}>280$ GeV	279.1	$1.35\times {10}^4$	2.38
$\mathrm{THT}>200$ GeV	277.7	$7.34\times {10}^3$	3.18
$E_T\left(jj\right)>100$ GeV	263.7	$5.84\times {10}^3$	3.38
$E_T\left(e^{-}j\right)>120$ GeV	261.9	$5.61\times {10}^3$	3.42

Table 6. Table 5: Effect of individual kinematical cuts on the signal-2 for M Z μ ​ τ = 0.1 subscript 𝑀 subscript 𝑍 𝜇 𝜏 0.1 {M}_{Z_{\mu\tau}}=0.1 GeV and χ = 9 × 10 − 5 ​ GeV − 1 𝜒 9 superscript 10 5 superscript GeV 1 \chi=9\times 10^{-5}~{}\rm{GeV^{-1}} and backgrounds at the FCC-eh. The statistical significance S ​ S 𝑆 𝑆 SS is computed for an integrated luminosity of 1 ab -1 .

FCC-eh,  $E_{e^{-}}=60$ GeV,  $E_p=50$ TeV
cuts	signal (S)	total background (B)	$S/\sqrt{S+B}$
initial (no cut)	2412.0	$1.82\times {10}^6$	1.79
basic cuts	2246.7	$1.30\times {10}^6$	1.96
GeV	1086.9	$1.29\times {10}^5$	3.01
$\mathrm{TET}>180$ GeV	981.0	$3.76\times {10}^4$	5.00
$\mathrm{THT}>120$ GeV	972.0	$2.46\times {10}^4$	6.08
$E_T\left(jj\right)>60$ GeV	959.3	$1.92\times {10}^4$	6.76
$E_T\left(e^{-}j\right)>80$ GeV	952.9	$1.79\times {10}^4$	6.94

Equations

$$\mathcal{L}=$$

$${\mathcal{L}}_{DM}=(D^{\nu }{\phi }_{DM}{)}^{†}(D_{\nu }{\phi }_{DM})-{\mu }_{DM}^2{\phi }_{DM}^{†}{\phi }_{DM}-{\lambda }_{DM}({\phi }_{DM}^{†}{\phi }_{DM}{)}^2$$

$$-{\lambda }_{Dh}({\phi }_{DM}^{†}{\phi }_{DM})({\phi }_h^{†}{\phi }_h)-{\lambda }_{DH}({\phi }_{DM}^{†}{\phi }_{DM})({\phi }_H^{†}{\phi }_H),$$

$$V={\mu }_H^2{\phi }_H^{†}{\phi }_H+{\lambda }_H({\phi }_H^{†}{\phi }_H{)}^2+{\lambda }_{hH}({\phi }_h^{†}{\phi }_h)({\phi }_H^{†}{\phi }_H).$$

$$\displaystyle\!\!\!\!\!\varphi_{h}=\left(\begin{array}[]{c}\tilde{H}\cr\frac{v+H+iA}{\sqrt{2}}\cr\end{array}\right)\;,$$

$$\displaystyle\quad\mathcal{M}^{2}_{scalar}=\left(\begin{array}[]{cc}2{\lambda}_{h}v^{2}&{\lambda}_{hH}v_{\mu\tau}v\cr{\lambda}_{hH}v_{\mu\tau}v&2{\lambda}_{H}v^{2}_{\mu\tau}\cr\end{array}\right)\;.$$

$$h_1=H\cos \alpha +H_{\mu \tau }\sin \alpha ,h_2=-H\sin \alpha +H_{\mu \tau }\cos \alpha ,$$

$$\tan 2\alpha =\frac{{\lambda }_{hH}{v}_{\mu \tau }v}{{\lambda }_h{v}^2-{\lambda }_H{v}_{\mu \tau }^2}.$$

$$M_{h_1}^2$$

$$M_{h_2}^2$$

$$g_{Z_{\mu \tau }{Z}_{\mu \tau }{h}_1}=\frac{2M_{Z_{\mu \tau }}^2}{v_{\mu \tau }}\sin \alpha ,g_{f\bar{f}{h}_1}=-\frac{M_f}{v}\cos \alpha ,g_{VV{h}_1}=\frac{2M_V^2}{v}\cos \alpha ,$$

$$g_{{\phi }_{DM}^{†}{\phi }_{DM}{h}_1}=-(v{\lambda }_{Dh}\cos \alpha +v_{\mu \tau }{\lambda }_{DH}\sin \alpha ),$$

$$g_{Z_{\mu \tau }{Z}_{\mu \tau }{h}_2}=\frac{2M_{Z_{\mu \tau }}^2}{v_{\mu \tau }}\cos \alpha ,g_{f\bar{f}{h}_2}=\frac{M_f}{v}\sin \alpha ,g_{VV{h}_2}=-\frac{2M_V^2}{v}\sin \alpha ,$$

$$g_{{\phi }_{DM}^{†}{\varphi }_{DM}{h}_2}=(v{\lambda }_{Dh}\sin \alpha -v_{\mu \tau }{\lambda }_{DH}\cos \alpha ),$$

$$g_{h_1{h}_1{h}_2}=6v{\lambda }_h{\cos }^2\alpha \sin \alpha -6v_{\mu \tau }{\lambda }_H{\sin }^2\alpha \cos \alpha -2v{\lambda }_{hH}\sin \alpha$$

$$+6v{\lambda }_{hH}{\sin }^3\alpha -v_{\mu \tau }{\lambda }_{hH}\cos \alpha +3v_{\mu \tau }{\lambda }_{hH}{\sin }^2\alpha \cos \alpha .$$

$$g_{Z_{\mu \tau }{v}_{\mu }{v}_{\mu }}=\frac{M_{Z_{\mu \tau }}}{v_{\mu \tau }},g_{Z_{\mu \tau }{v}_{\tau }{v}_{\tau }}=-\frac{M_{Z_{\mu \tau }}}{v_{\mu \tau }}.$$

$${\Gamma }_{Z_{\mu \tau }}=\frac{g_{\mu \tau }^2{M}_{Z_{\mu \tau }}}{12\pi },$$

$$g_{\mu \tau }\simeq [2\times 10^{-4},2\times 10^{-3}],M_{Z_{\mu \tau }}\simeq [5,210]\mathrm{MeV}.$$

$$v_{\mu \tau }=\frac{M_{Z_{\mu \tau }}}{g_{\mu \tau }}\simeq [10,1000]\mathrm{GeV},$$

$$\Gamma (h_1\to Z_{\mu \tau }{Z}_{\mu \tau })=\frac{g_{Z_{\mu \tau }{Z}_{\mu \tau }{h}_1}^2(M_{h_1}^4-{4M}_{Z_{\mu \tau }}^2{M}_{h_1}^2+12M_{Z_{\mu \tau }}^4)\sqrt{M_{h_1}^2-{4M}_{Z_{\mu \tau }}^2}}{128\pi M_{h_1}^2{M}_{Z_{\mu \tau }}^4},$$

$$\Gamma (h_1\to {\phi }_{DM}^{†}{\phi }_{DM})=\frac{g_{{\phi }_{DM}^{†}{\phi }_{DM}{h}_1}^2}{32\pi M_{h_1}}\sqrt{1-\frac{{4M}_{DM}^2}{M_{h_1}^2}},$$

$$\Gamma (h_1)={\cos }^2\alpha {\Gamma }_{SM}+\Gamma (h_1\to Z_{\mu \tau }{Z}_{\mu \tau })+\Gamma (h_1\to {\phi }_{DM}^{†}{\phi }_{DM}),$$

$$BR(h_1\to \mathrm{invisibles})=\frac{\Gamma (h_1\to Z_{\mu \tau }{Z}_{\mu \tau })+\Gamma (h_1\to {\phi }_{DM}^{†}{\phi }_{DM})}{\Gamma (h_1)}.$$

$$\chi =\frac{\alpha }{v_{\mu \tau }}\le 2.2\times 10^{-4}\mathrm{Ge}{\mathrm{V}}^{-1}.$$

$$\Gamma (h_2\to Z_{\mu \tau }{Z}_{\mu \tau })=\frac{g_{Z_{\mu \tau }{Z}_{\mu \tau }{h}_2}^2(M_{h_2}^4-{4M}_{Z_{\mu \tau }}^2{M}_{h_2}^2+12M_{Z_{\mu \tau }}^4)\sqrt{M_{h_2}^2-{4M}_{Z_{\mu \tau }}^2}}{128\pi M_{h_2}^2{M}_{Z_{\mu \tau }}^4}.$$

$$\Gamma (h_2\to Z_{\mu \tau }{Z}_{\mu \tau })=\frac{M_{h_2}^3{\cos }^2\alpha }{32\pi v_{\mu \tau }^2}.$$

$$\Gamma (h_2\to {\phi }_{DM}^{†}{\phi }_{DM})=\frac{g_{{\phi }_{DM}^{†}{\phi }_{DM}{h}_2}^2}{16\pi M_{h_2}}\sqrt{1-\frac{{4M}_{DM}^2}{M_{h_2}^2}}.$$

$$\Gamma (h_2\to VV)=\frac{g_{VV{h}_2}^2(M_{h_2}^4-{4M}_V^2{M}_{h_2}^2+12M_V^4)\sqrt{M_{h_2}^2-{4M}_V^2}}{64\pi S_V{M}_{h_2}^2{M}_V^4},$$

$$\Gamma (h_2\to h_1{h}_1)=\frac{g_{h_1{h}_1{h}_2}^2\sqrt{M_{h_2}^2-{4M}_{h_1}^2}}{32\pi M_{h_2}^2}.$$

$${\lambda }_H$$

$${\lambda }_h$$

$${\lambda }_{hH}$$

$$\Gamma (h_2\to f\bar{f})=\frac{n_c{g}_{f\bar{f}{h}_2}(M_{h_2}^2-{4M}_f^2{)}^{\frac{3}{2}}}{8\pi M_{h_2}^2},$$

Full text

I. Introduction

As is known, the standard model (SM) of particle physics is one of the most successful theories over the past decades which describes a variety of experimental results over the wide range of energy scale from eV to TeV. Discovery of the 125 GeV Higgs boson at the Large Hadron Collider (LHC) in 2012 [1, 2] proves the success of the SM once again. However, so far the SM still has certain limitations. Some experimental facts have been plaguing people, and there is an urgent need to extend the SM. For instance, the sub-eV masses and peculiar mixing pattern of neutrinos  [3, 4], the muon $(g-2)$ anomalous magnetic moment [5], the exploration of dark matter (DM) [6] and dark energy [7, 8], the baryon asymmetry of the Universe and so on. Furthermore, discovery of Higgs boson provides an outstanding portal to new physics (NP) beyond the SM. Precision measurements of the Higgs boson properties are also one of the most important tasks of high-energy particle physics due to its possible role as portal to beyond the SM (BSM) sectors [9, 10, 11, 12, 13].

So far, there are many well motivated extensions of the SM, such as SUSY [14, 15, 16, 17], two Higgs doublet model [18, 19, 20, 21, 22, 23], and extension of the SM with an extra $U(1)$ gauge group [24, 25, 26, 27, 28, 29, 30]. In this work, we will consider the gauged ${U(1)}_{L_{\mu}-L_{\tau}}$ extension of the SM due to its relatively simple theoretical structure, which has a complete gauge group $SU(3)_{C}\times SU(2)_{L}\times U(1)_{Y}\times{U(1)}_{L_{\mu}-L_{\tau}}$ and is called the ${U(1)}_{L_{\mu}-L_{\tau}}$ model [31, 32, 33, 34, 35, 36, 37, 38, 39]. One of the advantages of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model is that the anomaly cancellation does not require any extra chiral fermionic degrees of freedom. In this model, the breaking of ${L_{\mu}-L_{\tau}}$ symmetry conduces to additional terms in the neutrino mass matrix, which offers an explanation for the neutrino masses and mixing simultaneously [37, 38]. Besides, the scalar sector has been expanded by two additional complex scalar singlets ($\varphi_{H}$ and $\varphi_{DM}$) with nonzero ${L_{\mu}-L_{\tau}}$ charge. The scalar $\varphi_{DM}$ can act as a viable DM candidate [40]. The other scalar $\varphi_{H}$ acquires a vacuum expectation value (VEV) $v_{\mu\tau}$ and thereby making an additional neutral scalar $h_{2}$ after spontaneous breaking of ${U(1)}_{L_{\mu}-L_{\tau}}$, which indicates that $h_{2}$ has a mass of the same order with $v_{\mu\tau}$ about 10 GeV - 1000 GeV [36, 41]. In principle, the additional neutral scalar can be produced and decay via their mixing with the SM-like Higgs boson $h_{1}$ [51].

On the other hand, an extra neutral gauge boson $Z_{\mu\tau}$ is also introduced and obtains a mass after spontaneous symmetry breaking of ${U(1)}_{L_{\mu}-L_{\tau}}$. $Z_{\mu\tau}$ does not couple to the SM quarks and the first generation leptons, which makes it avoid restrictions coming from lepton and hadron colliders such as LEP and LHC. Therefore, the mass of $Z_{\mu\tau}$ can be as light as 100 MeV for a low value of gauge coupling $g_{\mu\tau}\leq 10^{-3}$, which is required to meet the limits arising from neutrino trident production. The $Z_{\mu\tau}$ with an MeV-scale mass can resolve the muon $(g-2)$ anomaly, explain the deficit of cosmic neutrino flux [37, 43, 44] and resolve the problem of relic abundance of DM in the scenario with a light weakly interacting massive particle [45, 46, 47] simultaneously. Therefore, searching for its possible collider evidences plays a vital role in exploring NP. Many attempts to discover this kind of new particles have been made in the meson decay experiment [48], beam dump experiment [49], electron-positron collider experiments [50] and so on.

Searches for the new particles predicted by the ${U(1)}_{L_{\mu}-L_{\tau}}$ model are presently being conducted at the LHC and ILC [36]. While, another Higgs factory besides the LHC and ILC, such as the LHeC (Large Hadron electron Collider) and FCC-eh (Future Circular Collider in hadron-electron mode) [51, 52, 54, 53], could precisely determine their specific properties. In this paper, we mainly devote to study of the $h_{2}$ and $Z_{\mu\tau}$ productions and further explore the possibility of detecting their signatures at e-p colliders. We present a full simulation study of the production cross sections of $h_{2}$ and $Z_{\mu\tau}$ with the beam polarization P($e^{-}$)= -0.8. Then, we investigate their observability through the processes $e^{-}p\to{e^{-}jh_{2}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$, $e^{-}p\to{e^{-}jh_{1}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ and $e^{-}p\to{\nu j}{h_{2}(\to{ZZ})\to}~{}2l^{+}2l^{-}j+\cancel{E}_{T}$, respectively. We further analyze the signal significance of $h_{2}$ and $Z_{\mu\tau}$ detecting which depends on the free parameters. Our numerical results show that the signals of $h_{2}\to Z_{\mu\tau}Z_{\mu\tau}$ and $h_{1}\to Z_{\mu\tau}Z_{\mu\tau}$ are promising to be detected at e-p colliders with appropriate parameter values and high integrated luminosity. But, due to the interference of substantial backgrounds and the low number of events, searching for the signal of the decay channel $h_{2}\to ZZ$ are harder to achieve at e-p colliders.

Rest of the paper has been arranged in the following manner. In Sec. II, we briefly review the basic features of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model and show the allowed parameter space of this model. In Sec. III, we not only give the partial widths of the main decay channels of the scalar $h_{2}$, but also calculate its production cross sections via the $W^{+}W^{-}$ and $ZZ$ fusion processes. The production cross sections of the new gauge boson $Z_{\mu\tau}$ through $h_{2}$ and $h_{1}$ decays are calculated in Sec. IV. We estimate the numbers of the signal and background events, and investigate the signal observability and discovery potentiality of $h_{2}$ and $Z_{\mu\tau}$ through their respective promising production channels in Sec. V. Finally, our conclusions are given in Sec. VI.

II. The Basic Features of the ${U(1)}_{L_{\mu}-L_{\tau}}$ Model

The gauged ${U(1)}_{L_{\mu}-L_{\tau}}$ extension of the SM is one of the most extensively studied NP models, which can successfully solve the origin of tiny neutrino masses, the DM relic abundance and the muon $(g-2)$ anomalous magnetic moment. Refs.[37, 38] have made a detailed analysis about solving these puzzles in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model. In this model, the gauge sector of the SM is enhanced by imposing a local ${U(1)}_{L_{\mu}-L_{\tau}}$ symmetry to the SM Lagrangian, where $L_{\mu}$ and $L_{\tau}$ are the muon and tau lepton numbers, respectively. Therefore, the complete gauged group is $SU(3)_{C}\times SU(2)_{L}\times U(1)_{Y}\times{U(1)}_{L_{\mu}-L_{\tau}}$. The SM particle content has been extended by including three extra right-handed (RH) neutrinos and two SM gauge singlet scalars. All particles included in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model and their charge assignments under various symmetry groups are listed in Table 1.

The Lagrangian of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model is as follows

[TABLE]

In above equation we have ignored the kinetic-mixing term between the groups $U(1)_{Y}$ and ${U(1)}_{L_{\mu}-L_{\tau}}$ in the case of assuming that the mixing is very small. The terms $\mathcal{L}_{SM}$, $\mathcal{L}_{N}$ and $\mathcal{L}_{DM}$ represent the SM, right hand (RH) neutrino and DM sectors, respectively. Since the processes we are studying do not involve RH neutrinos, its specific form is not given here. $\mathcal{L}_{DM}$ represents the dark sector Lagrangian including the kinetic term of the DM candidate ${\varphi}_{DM}$ and the interaction terms of ${\varphi}_{DM}$ with the scalars fields ${\varphi}_{h}$ and ${\varphi}_{H}$. The expression of $\mathcal{L}_{DM}$ is given by

[TABLE]

where the parameters $\lambda_{DM}$, $\lambda_{Dh}$ and $\lambda_{DH}$ are quartic couplings of the scalar fields. As these couplings are feeble ($\sim 10^{-12}$) [37], the DM can not attain thermal equilibrium with the thermal soup, which is called the Feebly Interacting Massive Particle (FIMP). In Eq. (1), the covariant derivatives involving in the kinetic energy term $\left|{D}_{\nu}{\varphi}_{H}\right|^{2}$ of the extra Higgs singlet ${\varphi}_{H}$ can be expressed in a generic form $D_{\nu}\phi=(\partial_{\nu}+ig_{\mu\tau}Q_{\mu\tau}(\phi)Z_{\mu\tau\nu})\phi$, where $\phi$ is any SM single field which has ${U(1)}_{L_{\mu}-L_{\tau}}$ charge $Q_{\mu\tau}(\phi)$ (listed in Table 1) and $g_{\mu\tau}$ represents ${U(1)}_{L_{\mu}-L_{\tau}}$ group’s gauge coupling constant. The scalar potential $V$ contains all the self interactions of $\varphi_{H}$ and its interactions with SM Higgs doublet. Its expression form is given by

[TABLE]

The last term in Eq. (1) represents the kinetic term for the additional gauge boson $Z_{\mu\tau}$ in terms with field strength tensor $F^{\rho\sigma}_{\mu\tau}=\partial^{\rho}Z^{\sigma}_{\mu\tau}-\partial^{\sigma}Z^{\rho}_{\mu\tau}$ of the ${U(1)}_{L_{\mu}-L_{\tau}}$ gauge group. When the scalar field $\varphi_{H}$ has a non-zero of VEV, the ${U(1)}_{L_{\mu}-L_{\tau}}$ symmetry breaks spontaneously and consequently the corresponding new gauge boson $Z_{\mu\tau}$ obtains the mass ${M}_{Z_{\mu\tau}}=g_{\mu\tau}v_{\mu\tau}$. The SM Higgs doublet ${\varphi}_{h}$ and the new scalar ${\varphi}_{H}$ take the following form

[TABLE]

where $\tilde{H}$, $A$ and $a$ are the massless Nambu-Goldstone Bosons (NGBs) absorbed by the gauge bosons $W^{\pm}$, $Z$ and $Z_{\mu\tau}$, while $v$ and $v_{\mu\tau}$ are the VEVs of the scalars $\varphi_{h}$ and $\varphi_{H}$, respectively. Furthermore, $H$ and $H_{\mu\tau}$ represent the physical CP-even scalar bosons. When both $\varphi_{h}$ and $\varphi_{H}$ obtain their respective VEVs, there will be a mass mixing between the states $H$ and $H_{\mu\tau}$. The square of scalar mass matrix with off-diagonal elements proportional to $\lambda_{hH}$ is given by

[TABLE]

Rotating the basis states $H$ and $H_{\mu\tau}$ by a suitable angle $\alpha$, we can make the above mass matrix diagonal. The new basis states ($h_{1}$ and $h_{2}$), now representing two physical states, are the linear combinations of $H$ and $H_{\mu\tau}$ with the mixing angle $\alpha$ between $H$ and $H_{\mu\tau}$, which can be expressed as

[TABLE]

[TABLE]

When $\alpha\ll 1$, $h_{1}$ can be identified as the SM-like Higgs boson which has already been discovered by the CMS [1] and ATLAS [2] collaborations in 2012. $h_{2}$ is a new scalar particle. The masses of these two physical scalars $h_{1}$ and $h_{2}$ are given by

[TABLE]

In this paper, we will assume that the values of $M_{h_{1}}$ and $v$ are fixed at 125 GeV and 246 GeV, respectively.

On one hand, compared with the SM Higgs boson, the couplings of $h_{1}$ with the SM particles are suppressed by a factor $\cos\alpha$ in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model. Some relevant couplings of the SM-like Higgs boson $h_{1}$ with the SM particles, the new gauge boson $Z_{\mu\tau}$ and the DM candidate ${\varphi}_{DM}$ are given by

[TABLE]

where $f$ represents all of the SM fermions and $V$ represents the electroweak gauge bosons $W^{\pm}$ or $Z$. On the other hand, similar with $h_{1}$, the scalar $h_{2}$ can couple to all the SM particles and other new particles, such as new gauge boson $Z_{\mu\tau}$ and the DM particle ${\varphi}_{DM}$. Here, we also list the $h_{2}$ couplings, which are related our calculation

[TABLE]

In the ${U(1)}_{L_{\mu}-L_{\tau}}$ model, $Z_{\mu\tau}$ has a light mass and no couplings to the SM quarks and the first generation leptons, so it can only decay to neutrinos. The couplings of $Z_{\mu\tau}$ with neutrinos are expressed as

[TABLE]

Taking no account of the neutrino masses, the expression form of the total decay width of $Z_{\mu\tau}$ is given by

[TABLE]

where we have also ignored the neutrino mixing.

To produce the appropriate neutrino mass and explain the muon $(g-2)$ anomaly, Ref.[36] has show the favored regions of the gauge coupling $g_{\mu\tau}$ and the $Z_{\mu\tau}$ mass, which are summarized as

[TABLE]

According to Eq. (18), the range of $v_{\mu\tau}$ is given by

[TABLE]

which indicates that, after the spontaneous symmetry breaking, the new scalar $h_{2}$ obtains mass as the same order with $v_{\mu\tau}$.

From above discussions we can see that, besides decaying to SM particles, the SM-like Higgs boson $h_{1}$ has extra decay modes $Z_{\mu\tau}Z_{\mu\tau}$ and $\varphi^{\dagger}_{DM}\varphi_{DM}$ for $M_{DM}<M_{h_{1}}$. The expressions of these decay channels are

[TABLE]

[TABLE]

where the $M_{DM}$ is Dark Matter mass. The total decay width of $h_{1}$ can be written as

[TABLE]

where $\Gamma_{SM}$ is the total width of the Higgs boson in the SM. In the ${U(1)}_{L_{\mu}-L_{\tau}}$ model, the decays of $\varphi_{DM}$ and $Z_{\mu\tau}$ are invisible, the branching ratio of the invisible decays is given by

[TABLE]

As we can see from Eq. (14) and Eq. (15), the couplings of ${\varphi}_{DM}$ with scalar bosons $h_{1}$ and $h_{2}$ depend on the parameters $\lambda_{Dh}$ and $\lambda_{DH}$. In this work, we take $\lambda_{Dh}=9.8\times 10^{-13}$ and $\lambda_{DH}=1.3\times 10^{-11}$ [37], which makes the $\Gamma(h_{1}\to\varphi^{\dagger}_{DM}\varphi_{DM})$ and $\Gamma(h_{2}\to\varphi^{\dagger}_{DM}\varphi_{DM})$ so feeble that they can be ignored. Using the constraint on the branching ratio of the Higgs invisible decay, $\rm BR_{invis}\leq 0.24$ at $95\%$ C. L. from the LHC data [55], the sine of scalar mixing angle $\sin\alpha$ must be satisfied $\sin\alpha\leq 0.3$. Then, for the factor $\chi$, there is

[TABLE]

To summarize, in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model, three new free parameters are introduced, which are the new gauge coupling constant $g_{\mu\tau}$, the $Z_{\mu\tau}$ mass $M_{Z_{\mu\tau}}$ and the scalar mixing angle $\alpha$, respectively. In the following, we will focus our attention on the phenomenology of the new particles $h_{2}$ and $Z_{\mu\tau}$ at e-p colliders in the above allowed parameter space.

III. Decays and Productions of the Scalar $h_{2}$

**3.1. Decays of the scalar $h_{2}$ **

In the ${U(1)}_{L_{\mu}-L_{\tau}}$ model, the scalar $h_{2}$ can not only decay to the SM particles but also decay to the new particles $Z_{\mu\tau}$ and $\varphi_{DM}$. Here, we give the decay width expressions of its several major decay modes. The expression form of the decay width for the decay channel $h_{2}\to Z_{\mu\tau}Z_{\mu\tau}$ is given by

[TABLE]

Under the assumption $\frac{{M}_{Z_{\mu\tau}}}{{M}_{h_{1}}}\to 0$, we can obtain the following form

[TABLE]

The decay width of the channel $h_{2}\to\varphi^{\dagger}_{DM}\varphi_{DM}$ is

[TABLE]

As we already mentioned, the value of $\Gamma(h_{2}\to\varphi^{\dagger}_{DM}\varphi_{DM})$ is small enough and we will neglect it.

The width of $h_{2}$ decaying to vector bosons is given as

[TABLE]

where $S_{V}$ represents the statistical factor. Its value equals to 1 for $W^{\pm}$ boson and 2 for $Z$ boson.

The width for the decay process $h_{2}\to h_{1}h_{1}$ can be written as

[TABLE]

From Eq. (15) one can see that the coupling constant $g_{h_{1}h_{1}h_{2}}$ depends on the couplings $\lambda_{h}$, $\lambda_{H}$ and $\lambda_{hH}$, which are expressed as

[TABLE]

When we take $\sin\alpha=0.2$ and $v_{\mu\tau}=\frac{M_{Z_{\mu\tau}}}{g_{\mu\tau}}=100$ GeV, the decay width $\Gamma(h_{2}\to h_{1}h_{1})$ only depends on $M_{h_{2}}$.

The width of $h_{2}$ decaying to the SM fermion pair is given as

[TABLE]

where the color charge $n_{c}=1$ for leptons and 3 for quarks. Fig. 1 shows the branching ratios for the main decay modes of the scalar $h_{2}$ as functions of the mass parameter $M_{h_{2}}$ for the fixed values sin$\alpha$ = 0.2, $M_{Z_{\mu\tau}}=0.1$ GeV and $g_{\mu\tau}=1\times 10^{-3}$, where the curves from high to low correspond the $Z_{\mu\tau}Z_{\mu\tau}$ decay modes, the $W^{+}W^{-}/ZZ$ decay modes, di-higgs decay mode and di-top decay mode, respectively. One can see from this figure that the value of the branching ratio $BR(h_{2}\to Z_{\mu\tau}Z_{\mu\tau})$ is about $98\%$ and only is $2\%$ for the rest decay channels. Certainly, the values of these branching ratios would vary as the values of the parameters sin$\alpha$ and $g_{\mu\tau}$ changing. However, in the allowed parameter space of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model, the decay process $h_{2}\to Z_{\mu\tau}Z_{\mu\tau}$ is the main decay channel of the scalar $h_{2}$.

**3.2. Productions of the scalar $h_{2}$

**

Like the SM Higgs boson, the additional scalar $h_{2}$ in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model is produced via two major channels: the charged current (CC) production channel via $W^{+}W^{-}$ fusion and the neutral current (NC) production channel via $ZZ$ fusion [56, 57] at e-p colliders. Fig. 2 gives the corresponding Feynman diagrams for the $h_{2}$ production via CC production channel and NC production channel at e-p colliders, respectively.

Then, employing Madgraph5/aMC@NLO [58], we calculate the production cross sections of the processes $e^{-}p\to{\nu_{e}jh_{2}}$ and $e^{-}p\to{e^{-}jh_{2}}$ as functions of $M_{Z_{\mu\tau}}$ at the LHeC. It is well known that polarization of the initial state electron can affect the production cross sections. Our numerical results show that the beam polarization P($e^{-}$)= -0.8 can maximize the cross sections. Therefore, we will take $P(e^{-})=-0.8$ in following numerical calculation. Since $g_{\mu\tau}$ does not affect the cross sections of $h_{2}$ production via the CC and NC processes, we do not consider it here. In Figs. 3 (a) and  3 (b), the curves show the cross sections of the $e^{-}p\to{e^{-}jh_{2}}$ and $e^{-}p\to{\nu_{e}jh_{2}}$ processes with $E_{e^{-}}=140$ GeV and different values of mixing angle $\sin\alpha=0.2$ (solid), $0.05$ (dashed) and $0.01$ (dotted). One can see from these figures that the values of the production cross section $\sigma$ decrease as the $h_{2}$ mass increases. For the $e^{-}p\to{e^{-}jh_{2}}$ process and 10 GeV $\leq{M}_{h_{2}}\leq$ 1000 GeV, its values are in the ranges of $1.73~{}\rm{pb}\times 10^{-6}\leq\sigma\leq 4.33\times 10^{-3}$ pb (solid), $1.08\times 10^{-7}~{}\rm{pb}\leq\sigma\leq 2.71\times 10^{-4}$ pb  (dashed) and $4.33\times 10^{-9}~{}\rm{pb}\leq\sigma\leq 1.08\times 10^{-5}$ pb (dotted), respectively. For the $e^{-}p\to{\nu_{e}jh_{2}}$ process and 10 GeV $\leq{M}_{h_{2}}\leq$ 1000 GeV, its values are in the ranges of $1.77\times 10^{-5}~{}\rm{pb}\leq\sigma\leq 3.73\times 10^{-2}$ pb (solid), $1.11\times 10^{-6}~{}\rm{pb}\leq\sigma\leq 2.33\times 10^{-3}$ pb  (dashed) and $4.43\times 10^{-8}~{}\rm{pb}\leq\sigma\leq 9.32\times 10^{-5}$ pb (dotted) respectively. It is worth mentioning that the cross section of the $e^{-}p\to{\nu_{e}jh_{2}}$ process is larger than that of the $e^{-}p\to{e^{-}jh_{2}}$ process by about one order of magnitude.

IV. Productions of the New Gauge Boson $Z_{\mu\tau}$

Now, we turn our attention to the new gauge boson $Z_{\mu\tau}$. As mentioned in the previous section, $Z_{\mu\tau}$ can not establish couplings with all the SM quarks and the first generation leptons, making it very difficult to be produced directly. So it is a attractive scheme to obtain $Z_{\mu\tau}$ by considering its indirect production. Similar with the new scalar $h_{2}$, besides decaying to the SM particles, the SM-like Higgs boson $h_{1}$ can also decay to a pair of $Z_{\mu\tau}$. Eq. (15) has given the expression form of the decay width $\Gamma(h_{1}\to Z_{\mu\tau}Z_{\mu\tau})$, which can be simplified to

[TABLE]

From above equation we can see that the production rate of the $Z_{\mu\tau}$ pair from $h_{1}$ decaying is actually determined by the factor $\chi^{2}\simeq\sin^{2}\alpha/v^{2}_{\mu\tau}$. So, in this work, all the results for the $Z_{\mu\tau}$ production via $h_{1}$ decaying can be expressed as functions of the factor $\chi$. Next, we will consider its indirect productions via the decays of $h_{2}$ and $h_{1}$, respectively.

As can be seen from Fig. 3, the CC production of scalar $h_{2}$ has larger cross section than that for its NC production. However, as mentioned earlier, $Z_{\mu\tau}$ can only decay to neutrinos in the ${U(1)}_{L_{\mu}-L_{\tau}}$ model. So the final states of the $e^{-}p\to{\nu_{e}jh_{2}}(h_{2}\to{Z_{\mu\tau}Z_{\mu\tau}})$ and $e^{-}p\to{\nu_{e}jh_{1}}(h_{1}\to{Z_{\mu\tau}Z_{\mu\tau}})$ processes would be jets and missing energy, which are difficult to be distinguished from the deeply inelastic scattering (DIS) backgrounds. Moreover, lack of kinematic handles in the final state makes it extremely difficult to filter signal from many backgrounds. Therefore, in this work we will focus on the NC production channels $e^{-}p\to{e^{-}jh_{1}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ and $e^{-}p\to{e^{-}jh_{2}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ to study the feasibility of detecting $h_{2}$ and $Z_{\mu\tau}$ . In Fig. 4, we show the leading order Feynman diagrams of the $Z_{\mu\tau}$ productions by the decays of $h_{1}$ and $h_{2}$ at e-p colliders.

Employing the Madgraph5/aMC@NLO [58], we calculate the cross sections of the $Z_{\mu\tau}$ production processes $e^{-}p\to{e^{-}jh_{2}\to}~{}e^{-}jZ_{\mu\tau}Z_{\mu\tau}$ and $e^{-}p\to{e^{-}jh_{1}\to}~{}e^{-}jZ_{\mu\tau}Z_{\mu\tau}$. Considering the favored region of the parameter space to resolve $(g-2)_{\mu}$ discrepancy, $g_{\mu\tau}$ is fixed to $g_{\mu\tau}=1\times 10^{-3}$ for reference. Due to $BR(h_{2}\to Z_{\mu\tau}Z_{\mu\tau})\sim 1$, the cross section of the $Z_{\mu\tau}$ production via $h_{2}$ decaying as a function of the mass $M_{h_{2}}$ for $\sin\alpha=0.2$, $0.05$ and $0.01$ at LHeC with $E_{e^{-}}=140$ GeV is essentially the same as Fig. 3 (a). So we don’t show it again. Fig. 5 shows the cross sections of the $Z_{\mu\tau}$ production by $h_{1}$ decaying as functions of the factor $\chi$ at e-p colliders, where the different curves show the $Z_{\mu\tau}$ production cross sections at different colliders: FCC-eh (solid), LHeC with $E_{e^{-}}=140~{}\rm{GeV}$ (dashed) and  LHeC with $E_{e^{-}}=60~{}\rm{GeV}$ (dotted). From Fig. 5, we can see that, for $1\times 10^{-5}$ GeV $\leq\chi\leq$ $2\times 10^{-4}$ GeV, the values of the $Z_{\mu\tau}$ production cross sections are in the ranges $3.61\times 10^{-5}~{}\rm{pb}\leq\sigma\leq 1.75\times 10^{-2}$ pb (solid), $1.94\times 10^{-5}~{}\rm{pb}\leq\sigma\leq 9.37\times 10^{-3}$ pb  (dashed) and $8.16\times 10^{-6}~{}\rm{pb}\leq\sigma\leq 3.97\times 10^{-3}$ pb (dotted).

V. Signatures of the New Particles $h_{2}$ and $Z_{\mu\tau}$ at e-p Colliders

In this section, we analyze the observation potential by performing a Monte Carlo simulation of the signal and background events and explore the observability of the additional scalar $h_{2}$ and the new gauge boson $Z_{\mu\tau}$ at e-p colliders with the integrated luminosity of 1 $\rm{ab^{-1}}$. On one hand, we will explore the observability of $h_{2}$ as well as $Z_{\mu\tau}$ via the processes $e^{-}p\to{e^{-}jh_{2}\to}~{}e^{-}jZ_{\mu\tau}Z_{\mu\tau}$ and $e^{-}p\to{e^{-}jh_{1}\to}~{}e^{-}jZ_{\mu\tau}Z_{\mu\tau}$. For the second process, we mainly use it to explore the signature of $Z_{\mu\tau}$. On the other hand, for more comprehensive, we will also analysis the possibility of detecting $h_{2}$ through CC production channel followed by $h_{2}$ decaying into vector boson pair ($ZZ$).

We use Madgraph5/aMC@NLO [58] to calculate the relevant production cross sections and generate the signal and background events, where the UFO format of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model has been obtained by using FeynRules [59]. Moreover, the parton distribution function (PDF), NNPDF2.3 [60], is used at leading order and Pythia-pgs [61] is employed for parton showering, hadronization and fast detector simulation. Finally, MadAnalysis5 [62] is applied for data analysis and plotting. All of the SM input parameters are taken from Particle Data Group (PDG) [63].

**A. $h_{2}\to Z_{\mu\tau}Z_{\mu\tau}$ and $h_{1}\to Z_{\mu\tau}Z_{\mu\tau}$ Channels

**

In this subsection, we take both $h_{2}$ and $h_{1}$ productions at e-p colliders through NC production channels followed by $h_{2}\to{Z_{\mu\tau}Z_{\mu\tau}}$ and $h_{1}\to{Z_{\mu\tau}Z_{\mu\tau}}$ as our signals, signal-1 and signal-2, respectively. Since $Z_{\mu\tau}$ has invisible final state in the detector, these two processes provide the same final state that includes one electron, one jet and a large missing transverse energy $\cancel{E}_{T}$

[TABLE]

[TABLE]

in which $\cancel{E}_{T}$ comes from $Z_{\mu\tau}\to\nu\bar{\nu}$. For $\sin\alpha=0.2$ and $g_{\mu\tau}=1\times 10^{-3}$, the values of the cross section for the signal-1 are $4.179\times 10^{-4}$ pb ($1.046\times 10^{-4}$ pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and $9.251\times 10^{-4}$ pb at the FCC-eh . While the values of the cross section for the signal-2 are $1.898\times 10^{-3}$ pb ($5.441\times 10^{-4}$ pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and $2.412\times 10^{-3}$ pb at the FCC-eh for $\chi=9\times 10^{-5}~{}\rm{GeV^{-1}}$.

For the signal $e^{-}j\cancel{E}_{T}$, the leading irreducible SM backgrounds can be classified into two general categories. The first category has a final state $e^{-}j\nu_{e}\bar{\nu}_{e}$ which comes from the following two processes

[TABLE]

[TABLE]

The total cross section for this kind of irreducible backgrounds is 0.4334 pb (0.205 pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and 0.8116 pb at the FCC-eh, which will severely pollute the physical signal.

The second category has a final state $e^{-}j\nu_{\mu,\tau}\bar{\nu}_{\mu,\tau}$

[TABLE]

Its production cross section is 0.05685 pb (0.03422 pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and 0.1052 pb at the FCC-eh. Besides, more remarkably, the photoproduction of the state $W+j$, which has a larger cross section, is also an irreducible SM background if the $W$ boson decays to an electron and neutrino. But it can be negligible after all selection cuts, because of its unique kinematic features.

There are also some reducible backgrounds which come from various sources. The most threatening reducible backgrounds result from the production of $\tau$ in the final state. One is

[TABLE]

Its cross section is 0.264 pb (0.1331 pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and 0.3957 pb at the FCC-eh. The other one is

[TABLE]

Its cross section is 0.2816 pb (0.1362 pb) for $E_{e^{-}}=140~{}(60)$ GeV at the LHeC and 0.5015 pb at the FCC-eh. The main reasons why the above two processes (Eq. (38) and Eq. (39)) can be viewed as reducible backgrounds are: (I) The $\tau$-jets may be misidentified as hadronic jets. (II) The detection of hadronic decay products of $\tau$ cannot be expected to be fully efficient due to the products being too soft, which will lead to generation of the missing energy $\cancel{E}_{T}$. Furthermore, we can even consider the case (II) as a source of partial irreducible background. The $e^{-}p\to\nu_{e}j\tau^{+}\nu_{\tau}$ and $e^{-}p\to\nu_{e}j\tau^{-}\bar{\nu}_{\tau}$ processes are reducible backgrounds in which the $\tau$ decays to an electron. Fortunately, we could suppress them to an insignificant order because of the totally different kinematic distribution of the final electron. Some other reducible backgrounds are $e$+multijet productions in which the $\cancel{E}_{T}$ comes from jet’s mismeasurement and $jj\nu$ production in which one jet is misidentified as an electron. In this work, we do not simulate both of them because their contributions can be negligible after all selection cuts. The signal and background events are generated with following basic cuts [56] in Madgraph5/aMC@NLO [58]

• lepton transverse momentum   $p_{T}(l^{\pm})>5$ GeV,

• jet transverse momentum   $p_{T}(j)>20$ GeV,

•   lepton pseudorapidity in the range $|\eta(l^{\pm})|<5$ ,

•   jet pseudorapidity in the range $|\eta(j)|<5$ ,

•   angular separation between jet and lepton $\Delta R(jl^{\pm})>0.4$ ,

where $\eta=1/2\ln(\tan\theta)$ is the pseudorapidity, where $\theta$ indicates the scattering angle in the laboratory frame. $\Delta R=\sqrt{(\Delta\phi)^{2}+(\Delta\eta)^{2}}$ is the particle separation, where $\Delta\phi$ and $\Delta\eta$ represent the rapidity gap and the azimuthal angle gap between the particle pair, respectively.

After the basic cuts, we further employ optimized kinematical cuts on separating the signals from the SM backgrounds. In our theoretical framework, although the SM backgrounds have a huge effect on the signals, there are many kinematical differences between them that can be exploited. In Fig. 6, we show the normalized distributions of the total missing transverse energy $\cancel{E}_{T}$, the visible transverse energy TET, the missing transverse hadronic energy THT, jet pair transverse energy $E_{T}(jj)$ and electron jet transverse energy $E_{T}(e^{-}j)$ for the signals and backgrounds at the LHeC with $E_{e^{-}}=140$ GeV and an integrated luminosity of 1 ab*-1*. From these figures, we can see that the distributions of signals have good distinctions from the distributions of the relevant backgrounds (peaks locate in different locations). In principle, there are other variables which we can use to discriminate the signals from backgrounds. But, these variables are remarkably similar and can not work significantly better than above kinematic variables. After all these kinematical cuts are applied, the event numbers of signal-1, signal-2 and corresponding backgrounds are summarized in Table 2 and Table 3 for the LHeC with $E_{e^{-}}=140~{}(60)$ GeV, respectively. The values of the statistical significance $SS$ are also shown in these tables, which is defined as $SS=S/\sqrt{S+B}$ with $S$ and $B$ being the number of signal and background events, respectively.

On the other hand, it is well known that the FCC-eh collides electrons to protons with $E_{e^{-}}=60$ GeV and $E_{p}=50$ TeV, which is a typical deep inelastic facility with $\sqrt{s}\approx 3.5$ TeV. Therefore, we need to modify the above veto criteria and kinematic cuts to adjust the progressive detector simulation, because the FCC-eh has a higher proton beam energy than LHeC. The modified values for the kinematic cuts and the event numbers of signal-1, signal-2 and backgrounds are presented in Table 4 and Table 5, respectively.

From these tables, we can see that, for sin$\alpha=0.2$, ${M}_{Z_{\mu\tau}}=0.1$ GeV, $M_{h_{2}}=300$ GeV, $g_{\mu\tau}=1\times 10^{-3}$ and the integrated luminosity being 1 ab*-1*, the values of $SS$ for signal-1 can reach 2.65 (1.84) at the LHeC with $E_{e^{-}}=140~{}(60)$ GeV and 3.42 at the FCC-eh with $E_{e^{-}}=60~{}\rm{GeV}$, $E_{p}=50~{}\rm{TeV}$. The values of $SS$ for signal-2 can reach 8.06 (5.10) at the LHeC with $E_{e^{-}}=140~{}(60)$ GeV and 6.94 at the FCC-eh with $E_{e^{-}}=60~{}\rm{GeV}$, $E_{p}=50~{}\rm{TeV}$ when we take ${M}_{Z_{\mu\tau}}=0.1$ GeV, $\chi=9\times 10^{-5}~{}\rm{GeV^{-1}}$ and an integrated luminosity of 1 ab*-1*.

In Figs. 7 (a), 7 (b) and 7 (c), performing the scan over the parameter spaces of ${M}_{h_{2}}$ and sin$\alpha$ for $M_{Z_{\mu\tau}}=0.1$ GeV, we show the experimental evidence region ($3\leq\rm{SS}<5$) and experimental discovery region ($5\leq\rm{SS}$) of signal-1 at different e-p colliders with integrated luminosity being 1 ab*-1*. For $\sin\alpha\leq 0.3$, from Fig. 7 (a), we obtain the $h_{2}$ mass region of above $3\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 320~{}\rm{GeV}$ and above $5\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 270~{}\rm{GeV}$ at the LHeC with $E_{e^{-}}=60$ GeV. From Fig. 7 (b), we obtain the $h_{2}$ mass region of above $3\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 400~{}\rm{GeV}$ and above $5\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 310~{}\rm{GeV}$ at the LHeC with $E_{e^{-}}=140$ GeV. From Fig. 7 (c), we obtain the $h_{2}$ mass region of above $3\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 480~{}\rm{GeV}$ and above $5\sigma$ confidence level as $10{~{}\rm{GeV}}\leq{M}_{h_{2}}\leq 360~{}\rm{GeV}$ at the FCC-eh. Based on these numerical results, we can say that the possible signatures of $h_{2}$ and $Z_{\mu\tau}$ from signal-1 is limited in the lower ${M}_{h_{2}}$ range and could be detected at e-p colliders with an integrated luminosity of 1 ab*-1*. On the other side, the FCC-eh could offer a better detection capabilities than LHeC under the same integrated luminosity.

The required integrated luminosities for observing the new gauge boson $Z_{\mu\tau}$ from signal-2 at the 3$\sigma$ and 5$\sigma$ confidence levels at different e-p colliders are plotted as functions of $\chi$ in Fig. 8. We can see that one can obtain larger statistical significance for larger $\chi$ value within allowed parameter space. We can easily obtain $5\sigma$ statistical significance for taking $\chi\geq 5\times 10^{-5}$ within the designed luminosity region. Furthermore, the LHeC with $E_{e^{-}}=140$ GeV has the best sensitivity to the signal-2. However, it is likely that the properties of the SM-like Higgs boson will be narrowed before 1 ab*-1* of data are collected at an e-p collider. So we present Fig. 9 with no thought of $BR(h_{1}\to\rm{inv})$. In order to illustrate excluded regions of the free parameters ${M}_{Z_{\mu\tau}}$ and sin$\alpha$ for reaching a given statistical significance, in Figs. 9 (a–b), we plot the $3\sigma$ and $5\sigma$ discovery of signal-2 for $g_{\mu\tau}=1\times 10^{-3}$ in the sin$\alpha$-${M}_{Z_{\mu\tau}}$ plane at different e-p colliders with an integrated luminosity 1 ab*-1*. From these figures, we can see that one can easily obtain $5\sigma$ statistical significance for taking lower values of sin$\alpha$ in the parameter space of ${M}_{Z_{\mu\tau}}$. Thus, from a phenomenological point of view, the signal-2 is more likely to result in a detection of the new gauge boson $Z_{\mu\tau}$ at a lower integrated luminosity and more achievable experimental conditions.

**B. $h_{2}\to ZZ\to 4l$ Channel

**

In this subsection we proceed to investigate the prospects of e-p colliders in searching for the scalar boson $h_{2}$ by focusing on its leading decay mode $ZZ$

[TABLE]

In above equation we have required the two $Z$ bosons to decay leptonically, then the signal consists of two lepton pairs, one jet and a large missing transverse energy $\cancel{E}_{T}$.

For the signal $2l^{+}+2l^{-}+j+\cancel{E}_{T}$, the main sources of the irreducible backgrounds come from the processes $e^{-}+p\to\nu jZZ$, $e^{-}jZZ$, $e^{-}jW^{+}W^{-}$, $e^{-}jZW^{+}$, $\nu jW^{+}W^{-}$ and $\nu jZW^{-}$. Similar as above we can calculate the statistical significance easily for the luminosity of 1 ab*-1* at e-p colliders. Our numerical results reveal that it is challenging to discover the $h_{2}$ signature through the $h_{2}\to ZZ\to 4l$ channel at the e-p colliders for $M_{h_{2}}$ = 200 (600) GeV, with $\sigma 3$ ($\sigma 5$) level statistical significance. We find that detecting this kind of signal at $3\sigma$ level requires the integrated luminosity be larger than $10^{5}$ ab*-1* for the benchmark point defined by $M_{h_{2}}=200$ GeV, $\sin\alpha=0.2$ and $g_{\mu\tau}=1\times 10^{-3}$, which extremely outreaches the designed luminosity. Thus, we do not show the relevant numerical results.

Certainly, the scalar $h_{2}$ can also decay to the modes $W^{+}W^{-}$, $t\bar{t}$ and $h_{1}h_{1}$ as long as its mass is large enough. However, the $h_{2}\to WW\to\nu l2j$ or $h_{2}\to WW\to 2\nu 2l$ channel is much more difficult to reconstruct compared to the $h_{2}\to ZZ\to 4l$ channel due to the final state neutrino which escapes from the detector and makes it impossible to fully reconstruct the $h_{2}$ system. So we’re not going to consider the channel here even though it has a larger branching ratio than that of $h_{2}\to ZZ$. We have to say that the detection of the new scalar boson $h_{2}$ via it decaying to pair of SM particles at e-p colliders is difficult to achieve at present. Of course, we also don’t rule out the existence of some unique kinematic cuts that we do not take into account to optimize background suppression and improve signal observability.

VI. Conclusions

The ${U(1)}_{L_{\mu}-L_{\tau}}$ model, which can explain the muon $(g-2)$ anomaly, small neutrino masses and provide a candidate of DM, is phenomenologically rich and predictive. In this model, the additional scalar $h_{2}$ and gauge boson $Z_{\mu\tau}$ are obtained after spontaneous breaking of $L_{\mu}-L_{\tau}$ symmetry. New scalar $h_{2}$ mixing with the SM-like Higgs boson is helpful to improve the precision of Higgs boson measurements. Furthermore, the gauge boson $Z_{\mu\tau}$ possessing a mass around the MeV scale can explain the deficit of cosmic neutrino flux and resolve the problem of muon $(g-2)$ anomaly and relic abundance of DM simultaneously. So, studying these two new particles is of great significance for exploring this kind of new physics models. In this paper, we have studied the possibility of searching for the new particles $h_{2}$ and $Z_{\mu\tau}$ at e-p colliders. Since $Z_{\mu\tau}$ can not couple with the SM quarks and the first generation leptons, it is very difficult to be produced directly at colliders. So we consider its productions via decays of $h_{1}$ and $h_{2}$. Although the CC production of $h_{1}$ and $h_{2}$ have larger cross sections, their final states will generate mono-jet plus missing energy, which accidentally coincides with the DIS backgrounds. Therefore we focus on NC production channels $e^{-}p\to{e^{-}jh_{1}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ and $e^{-}p\to{e^{-}jh_{2}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$, which provide good kinematic handles to distinguish the signals from the SM backgrounds. In addition to this, we also study the CC production of $h_{2}$ and further consider its non-invisible decays(e.g. $Z$ gauge boson) as complementary.

After giving the decay width expressions of several main decay channels of new scalar $h_{2}$, we calculate the production cross sections of the processes $e^{-}p\rightarrow e^{-}jh_{2}$ and $e^{-}p\rightarrow\nu_{e}jh_{2}$ with the beam polarization P($e^{-}$)= -0.8 in the context of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model. The production cross sections of $Z_{\mu\tau}$ are further calculated. Then, we investigate the observability of $h_{2}$ and $Z_{\mu\tau}$ through the signal-1 from the process $e^{-}p\to{e^{-}jh_{2}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ and the signal-2 from the process $e^{-}p\to{e^{-}jh_{1}(\to{Z_{\mu\tau}Z_{\mu\tau}})\to}~{}e^{-}j+\cancel{E}_{T}$ at e-p colliders with 1 ab*-1* integrated luminosity. After simulating the signals as well as the relevant backgrounds, and applying suitable kinematic cuts on the variables $\cancel{E}_{T}$, $\rm{TET}$, $\rm{THT}$, $E_{T}(jj)$ and $E_{T}(e^{-}j)$, the values of the statistical significance $SS$ for signal-1 can reach 2.65 (1.84) at the LHeC with $E_{e^{-}}=140~{}(60)$ GeV and 3.42 at the FCC-eh with $E_{e^{-}}=60~{}\rm{GeV}$, $E_{p}=50~{}\rm{TeV}$ when we take sin$\alpha=0.2$, $g_{\mu\tau}=1\times 10^{-3}$ GeV, ${M}_{Z_{\mu\tau}}=0.1~{}\rm{GeV}$ and ${M}_{h_{2}}=300$ GeV. While for signal-2, its values can reach 8.06 (5.10) at the LHeC with $E_{e^{-}}=140~{}(60)$ GeV and 6.94 at the FCC-eh when we take $\chi=9\times 10^{-5}~{}\rm{GeV^{-1}}$ and ${M}_{Z_{\mu\tau}}=0.1~{}\rm{GeV}$. Performing the scan over all parameter space, we find that the signals of $h_{2}$ and $Z_{\mu\tau}$ from signal-1 is limited in the lower ${M}_{h_{2}}$ range and could be detected at e-p colliders with an integrated luminosity of 1 ab*-1*. The signal of $Z_{\mu\tau}$ might be easily detected via signal-2 at e-p colliders, while the LHeC with $E_{e^{-}}=140$ GeV has the best sensitivity to signal-2. In the end, we analysis the signals of $h_{2}$ through the CC production channel via its decaying into a pair of gauge bosons. However, due to the interference of many backgrounds and the low number of events, searching for these kind of signals are harder to achieve at e-p colliders. Thus, we expect that the possible signals of the ${U(1)}_{L_{\mu}-L_{\tau}}$ model might be detected at future e-p colliders via $h_{2}\to Z_{\mu\tau}Z_{\mu\tau}$ and $h_{1}\to Z_{\mu\tau}Z_{\mu\tau}$ channels.

ACKNOWLEDGEMENT

This work was supported in part by the National Natural Science Foundation of China under Grant Nos.11875157, 11847303 and 11605081.