Lepton flavor violation induced by neutral and doubly-charged scalars at future lepton colliders

# Lepton flavor violation induced by neutral and doubly-charged scalars at future lepton colliders

P. S. Bhupal Dev Department of Physics and McDonnell Center for the Space Sciences, Washington University, St. Louis, MO 63130, USA    Rabindra N. Mohapatra Maryland Center for Fundamental Physics, Department of Physics, University of Maryland, College Park, MD 20742, USA    Yongchao Zhang111Talk presented at the International Workshop on Future Linear Colliders (LCWS2018), Arlington, Texas, 22-26 October 2018. C18-10-22. Department of Physics and McDonnell Center for the Space Sciences, Washington University, St. Louis, MO 63130, USA Center for High Energy Physics, Peking University, Beijing 100871, China
July 3, 2019
###### Abstract

New physics scenarios beyond the Standard Model (SM) for neutrino mass mechanism often necessitate the existence of a neutral scalar and/or doubly-charged scalar , which couple to the SM charged leptons in a flavor violating way, while evading all existing constraints. Such scalars could be effectively produced at future lepton colliders like CEPC, ILC, FCC-ee and CLIC, either on-shell or off-shell, and induce striking charged lepton flavor violating (LFV) signals. We find that a large parameter space of the scalar masses and the LFV couplings can be probed at lepton colliders, well beyond the current low-energy constraints in the lepton sector. The neutral scalar explanation of the muon anomaly could also be directly tested.

## I Introduction

There are various theoretical models of new physics which lead to charged lepton flavor violation (cLFV) effects at an observable level, which is strongly implied by the observation of flavor violation in the neutral lepton sector, i.e. the solar, atmospheric and reactor neutrino oscillations PDG (). In the minimal extension of Standard Model (SM) with Dirac neutrinos, cLFV rates are highly suppressed due to the tiny neutrino masses. This makes the experimental searches for cLFV all the more interesting, because any observable effect must come from physics beyond the minimally extended SM related to the origin of neutrino mass. LFV often arises from the extension of the Higgs sector, which allows flavor-violating Yukawa couplings of new neutral and/or doubly-charged scalars beyond the SM. In this report, we summarize the main results of model-independent LFV prospects at future lepton colliders, such as the Circular Electron-Positron Collider (CEPC) CEPC-SPPCStudyGroup:2015csa (), International Linear Collider (ILC) Baer:2013cma (), Future Circular Collider (FCC-ee) Gomez-Ceballos:2013zzn () and Compact Linear Collider (CLIC) Battaglia:2004mw (), based on our recent papers Dev:2017ftk (); Dev:2018upe (). This complements previous studies of LFV at lepton colliders that are performed in the framework of effective four-fermion couplings Kabachenko:1997aw (); Cho:2016zqo (); Ferreira:2006dg (); Aranda:2009kz (); Murakami:2014tna () or in the context of flavor-violating SM Higgs decays Banerjee:2016foh (); Chakraborty:2016gff (); Chakraborty:2017tyb (); Qin:2017aju () and tau decays Hays:2017ekz () or with doubly-charged scalars Rodejohann:2010bv (); Nomura:2017abh (). Compared to the hadron colliders, the lepton colliders are generally very “clean” and the SM processes therein are well understood, which render them primary facilities to search for new physics via the LFV signals.

If any of the new neutral scalars is hadrophobic (i.e. couples dominantly to leptons), it could remain sufficiently light and contribute sizably to cLFV with the Yukawa couplings ( being the lepton flavor indices), while easily evading the direct searches at hadron colliders, as well as the low-energy quark flavor constraints, such as the rare flavor-changing decays and oscillations of and mesons. Some well-motivated examples include supersymmetric models with leptonic -parity violation susy (); susy2 (); susy3 (); susy4 (), left-right symmetric models Dev:2016dja (); Dev:2016vle (); Dev:2017dui (); Maiezza:2016ybz (), two-loop models for neutrino masses babu (), mirror models mirror (); mirror2 (); mirror3 (); mirror4 (); mirror5 (), and two-Higgs doublet models 2HDM (); 2HDM2 (), where the cLFV coupling might arise at tree and/or loop level.

The doubly-charged scalar and its LFV Yukawa couplings might be closely related to neutrino mass generation, e.g. in the type-II seesaw type2a (); type2b (); type2c (); type2d (); type2e () and its left-right extension LR (); Mohapatra:1974gc (); Senjanovic:1975rk (). The magnitudes of might be sizable and hence accessible at lepton colliders. There are of course constraints on some of these couplings from rare lepton decays like and , but they leave enough room for some of them being of order one. Furthermore, rare lepton decays generally probe products of two different couplings whereas the lepton collider probes them individually Dev:2017ftk (); Dev:2018upe (). We show the interesting ranges of effective Yukawa couplings and that can be measured in the planned lepton colliders such as CEPC and ILC and can provide new ways to test the underlying extended Higgs models. They will in any case provide complementary information to rare lepton decay constraints on the and couplings, which makes such studies interesting from the synergistic viewpoint of energy and intensity frontiers.

This paper is organized as follows: Section II is devoted to the beyond SM neutral scalar , where we consider the LFV from both the on-shell and off-shell production of at lepton colliders. As two benchmark setups for the future lepton colliders, we show the corresponding prospects at CEPC 240 GeV and ILC 1 TeV, with total integrated luminosity of 5 ab and 1 ab, respectively. In Section III, we focus on the single, pair and off-shell production of doubly-charged scalar at lepton colliders through the Yukawa couplings and the resultant LFV prospects at ILC and CEPC. We conclude in Section IV. More details can be found in Refs. Dev:2017ftk (); Dev:2018upe ().

## Ii LFV induced by new neutral scalar

Without loss of generality, we can write the effective Yukawa couplings of the neutral scalar to the charged leptons as

 LY = hαβ¯ℓα,LHℓβ,R + H.c.. (1)

Here for simplicity we assume the couplings are all real and chirality-independent and thus symmetric. The scalar may or may not be responsible for symmetry breaking and/or mass generation of other particles in realistic models, where it could be part of a singlet, doublet or triplet scalar field. We assume further that is CP-even and its mixing with and coupling to the SM Higgs is small. If the scalar is CP-odd, the limits and prospects derived here would not change significantly. Even though there are all varieties of stringent low-energy cLFV constraints, such as , only a few of them are directly relevant to the LFV prospects discussed below. With an ab level of integrated luminosity at future lepton colliders, a large parameter space of the scalar mass and could be probed, well beyond the current cLFV constraints, and complementary to the projected low-energy constraints from future experiments at the intensity frontier Calibbi:2017uvl (). In addition, the Lagrangian in Eq. (1) also gives rise to a one-loop contribution to the lepton anomalous magnetic moment. In particular, the scalar field coupling whose induced loop graph helps resolve the longstanding muon discrepancy can be directly tested at lepton colliders.

### ii.1 On-shell LFV

If kinematically allowed, the neutral scalar can be directly produced at lepton colliders, in association with a pair of flavor-changing leptons through the couplings in Eq. (1), i.e.

 e+e−→ℓ±αℓ∓β+H (2)

mediated by the SM photon or boson,222We have also the production channel  Dev:2018upe (). However, this process involves only the flavor-conserving coupling , and we will not consider it here. and

 e+e−→να¯νβ+H(with α=e or β=e) (3)

mediated by the SM boson. In future lepton colliders, high luminosity photon beams can also be obtained by Compton back-scattering of low-energy, high-intensity laser beam off the high-energy electron beam Ginzburg:1981vm (); Ginzburg:1982yr (); Telnov:1989sd (), and the neutral scalar can also be produced through the processes

 e±γ→ℓ±α+Handγγ→ℓ±αℓ∓β+H. (4)

Let us first consider the coupling , switching off all other (LFV) couplings. It should be emphasized that all the production amplitudes in Eqs. (2) to (4) depend only on the LFV couplings (here ), and thus could be easily made to satisfy the rare lepton decay constraints, such as those from and , which depend on the product . With vanishing or suppressed couplings to the quark sector, the conversion limits are irrelevant. Furthermore, for real Yukawa couplings, we do not either have any limits from electric dipole moment. Only the following constraints are relevant:

1. Muonium-antimuonium oscillation: This could occur in both and -channels, with the oscillation probability , and is constrained by the MACS experiment Willmann:1998gd ().

2. : The anomalous magnetic moment of electron receives a contribution from the loop, which is strongly constrained from one-electron quantum cyclotron experiments Hanneke:2008tm ().

3. : A -channel could mediate the scattering , which interferes with the SM diagrams in the -channel and is constrained from LEP data Abdallah:2005ph ().

All these constraints on and are shown as the shaded regions in the upper two panels of Fig. 1 (which correspond to two benchmark configurations of future lepton colliders). To explain the longstanding theoretical and experimental discrepancy of the muon , i.e.  PDG (), the LFV coupling is required to be larger, as shown by the green band in the upper panels of Fig. 1, which is already excluded by the electron and muonium oscillation constraints.

As for the future lepton collider sensitivity, we specifically consider two benchmark configurations, namely, the CEPC 240 GeV with an integrated luminosity of 5 ab and ILC 1 TeV with a luminosity of 1 ab. The SM background is dominated by particle mis-identification, e.g. electron mis-identified as muon or vice versa, and is very small (see e.g. Refs. bkg1 (); bkg2 (); bkg3 ()). After being produced, could decay back into the charged lepton pairs or other SM particles. To be concrete, we consider the most optimistic case, i.e. the neutral scalar decays predominantly into a pair of leptons, i.e. (here decays only into pair). We apply nominal cuts GeV, and using CalcHEP Belyaev:2012qa (). The corresponding prospects of and in the channels of Eqs. (2) to (4) at CEPC and ILC are shown respectively in the upper left and upper right panels of Fig. 1, where we have assumed a minimum of 10 signal events at both colliders. It is clear from Fig. 1 that a large region of and can be probed in future lepton colliders, which significantly extends the current limits. Reconstructing the peak from the decay products could further improve the significance of the LFV signals, which however will be rather model-dependent.

Turning now to the coupling , the most stringent limits come from the electron  Hanneke:2008tm (), and the LEP data Abdallah:2005ph (). The reconstruction of lepton is more challenging than ,333Following Ref. Baer:2013cma (), we have adopted an efficiency factor of for the tau lepton reconstruction. but there is still ample parameter space to probe at both CEPC and ILC, as shown in the middle panels of Fig. 1, where we have assumed is the predominant decay mode.

Regarding the coupling , there are currently no experimental limits, except for the muon discrepancy PDG (). This could be explained in presence of when it couples to muon and tau, as shown by the brown line (for the central value) and the green and yellow bands for the and confidence level (CL), respectively, in the lower panels of Fig. 1, while the shaded region is excluded by the current muon data at the 5 CL. It is clear in Fig. 1 that the anomaly can be directly tested at CEPC (ILC) up to a scalar mass of GeV, in particular in the channel. With a larger luminosity being planned Gomez-Ceballos:2013zzn () at FCC-ee and a higher energy at CLIC, these prospects could be even better.

### ii.2 Off-shell LFV

The LFV signals could also be produced from an off-shell , i.e. . This could occur in both and channels, whereas in the -channel, is on-shell if the colliding energy , which corresponds to the resonant production with the subsequent decay of . Different from the on-shell case, the off-shell production amplitudes have a quadratic dependence on the Yukawa couplings (some of them might be flavor conserving), and thus largely complementary to the on-shell LFV searches.

The amplitude is proportional to . This is tightly constrained by the data PDG (), leaving no hope to see any collider LFV signal in this channel. In the tau sector, the LFV decay constraints are comparatively weaker. The limits on from and are almost constants, as in effective field theories with superheavy mediators. As for the on-shell case above, the coupling are constrained by electron and the LEP data Abdallah:2005ph (). All the constraints are presented in the top panel of Fig. 2, with the shaded regions excluded. Note that the muon can not be used to set unambiguous limits on the combinations and , although the couplings and could contribute to by themselves.

The dominant SM background is from the process , which is expected to be small, if we require the two charged leptons to be back-to-back and their reconstructed energy . Assuming 10 signal events as above, the prospects of the coupling are shown in the left panel of Fig. 2. At the resonance , the production cross section can be greatly enhanced by . To be specific, we have set the width (30) GeV at GeV (1 TeV). For , the production rate diminishes rapidly as becomes heavier. An off-shell could however be probed up to a few-TeV range, as shown in Fig. 2, and ILC is expected to be more promising than CEPC in this mass range, as a result of higher .

The process could proceed via both and -channels, which depend on different couplings, namely and , and are constrained respectively by the rare decays and . Analogous to the case above, a broad range of and could be probed in the -channel, especially in vicinity of the resonance, as shown by the bottom left panel of Fig. 2. In the -channel, the cross sections are comparatively smaller, and the detectable regions are much narrower, as shown by the bottom right panel of Fig. 2. One can see that orders of magnitude of the couplings can be probed in the off-shell channels at future lepton colliders, i.e. with couplings from up to for a mass range from 100 GeV to few TeV.

## Iii LFV induced by doubly-charged scalar

In addition to the neutral scalar case, LFV could also be induced by the presence of a doubly-charged scalar, with the model-independent Yukawa couplings

 LY = fαβH++¯ℓcαℓβ+H.c., (5)

with , the flavor indices and superscript denoting charge conjugate. The doubly-charged scalar can be either left-handed or right-handed, i.e. coupling either to the left-handed or right-handed charged leptons. We take the doubly-charged scalar in Eq. (5) to be right-handed, unless otherwise specified. The main results in this section hold also for left-handed .

At high-energy colliders, the “smoking-gun” signal of are same-sign charged lepton pairs from the decay , mediated by the interaction in Eq. (5). The current dilepton limits at the Large Hadron Collider (LHC) exclude the right-handed (left-handed) with mass up to roughly 650 GeV (800 GeV) ATLAS:2017iqw (); CMS:2017pet (), depending on the specific lepton flavors. The final states involving only the and flavors are the most stringent ATLAS:2017iqw () and those with the flavor are much weaker CMS:2017pet (), as a result of the poor reconstruction efficiency of the lepton at hadron colliders.

There are also direct searches of doubly-charged scalars at LEP in the single Abbiendi:2003pr () and pair Achard:2003mv () production modes. However, limited by the center-of-mass energy, these constraints are very weak. An off-shell in the -channel could mediate the Bhabha scattering and interfere with the SM diagrams. This alters both the total cross section and the differential distributions Abbiendi:2003pr (); Achard:2003mv (). By Fierz transformation, the doubly-charged scalar contributes to the effective four-fermion contact interaction

 1Λ2eff(¯eRγμeR)(¯fRγμfR), (6)

where corresponds to the effective cutoff scale related to the doubly-charged scalar mass and the Yukawa couplings, and is constrained by the (with ) data Abdallah:2005ph (). It turns out that the cutoff scale has been excluded up to 7.6 TeV by the LEP data.

For , the effective interaction in Eq. (6) would also induce an additional contribution to the Møller scattering and can be constrained by the upcoming MOLLER experiment Benesch:2014bas (), which could probe the effective scale TeV, slightly stronger than the current limit from LEP data Dev:2018sel (). Other low-energy LFV constraints on the couplings include the rare decays of charged leptons like and , the anomalous magnetic moments of electron and muon, the LEP data (with Abdallah:2005ph (), and muonium-anti-muonium oscillation, which are all highly suppressed in the SM. The calculation details can be found in Ref. Dev:2018upe () and references therein.

The doubly-charged scalar can be pair-produced at lepton colliders through the Drell-Yan process via the gauge interactions with the SM photon and boson. Drell-Yan production is, however, not a useful probe of the Yukawa interactions unless is long-lived, because for Drell-Yan pair production of , both larger and smaller Yukawa couplings give rise to the same line-shape over the background, and thus, we can not determine the absolute values of from the flavor-dependent branching fractions , as long as the leptonic decays dominate. Only when the Yukawa couplings enter the production of , could we measure these couplings in the (LFV) processes. All the relevant production processes can be categorized into two groups: i) the on-shell production including both the pair and single production modes and ii) the off-shell production of .

### iii.1 On-shell LFV

Besides the Drell-Yan process,444The doubly-charged scalar can also be pair-produced from the fusion of laser photons , induced by both the trilinear and quartic gauge-scalar interactions Chakrabarti:1998qy (). the doubly-charged scalar can be pair-produced at lepton colliders through the Yukawa interactions to the SM charged fermions Kuze:2002vb (),

 e+e−→H++H−−(through Yukawa couplings feℓ) (7)

with covering all three flavors , and . Given the Yukawa couplings to the SM charged fermions, the doubly-charged scalar can be singly produced from Barenboim:1996pt (); Kuze:2002vb (); Lusignoli:1989tr (); Yue:2007kv (); Dev:2018upe ()

 e+e−,γγ→H±±ℓ∓αℓ∓β (8)
 e±γ→H±±ℓ∓α, (9)

with the laser photon beam as in the neutral scalar case. With the same cuts as in Section II, we have estimated the prospects of all the LFV couplings at ILC 1 TeV with the integrated luminosity of 1 ab, as functions the doubly-charged mass , in all the pair and single production modes in Eqs. (7) to (9), which are collected in Fig. 3. For simplicity, we have assumed that decays predominately back to the same-sign charged leptons, and the electron and positron beams are both unpolarized. There is a lower cut of 500 GeV on the single production mode in the and single production processes, as below GeV, the processes are dominated by the pair production with from the decay of an on-shell . The limits on the LFV couplings are mainly from the LEP data, as a -channel could mediate the processes with the LFV couplings , which exclude the shaded pink regions in Fig. 3. The limits from electron and muon are highly suppressed by the charged lepton masses and are not shown in these plots. One should note that the contributions from doubly-charged scalar loops to the magnetic dipole moment are always negative Lindner:2016bgg (), so the muon anomaly can not be explained by the doubly-charged scalar . The vertical dashed gray lines indicate the current same-sign dilepton limits on the doubly-charged scalar mass from LHC ATLAS:2017iqw (); CMS:2017pet (), assuming the dilepton branching of . More details can be found in Ref. Dev:2018upe (). For a smaller dilepton branching ratio, the LHC limits on doubly-charged scalar mass will be weaker Dev:2018kpa ().

### iii.2 Off-shell LFV

With a -channel doubly-charged scalar , we can have the off-shell LFV processes Godfrey:2001xb (); Godfrey:2002wp (); Rizzo:1981xx (); Dev:2018upe ()

 e+e−→ℓ±αℓ∓βand% e±γ→ℓ∓αℓ±βℓ±γ (10)

which depend quadratically on the couplings . As in the neutral scalar case, the limit from is so stringent that it has precluded the doubly-charged scalar-mediated signal at future lepton colliders. The processes () and () involve the lepton and depend on the coupling (), and their prospects at CEPC 240 GeV and ILC 1 TeV are shown in the left (right) panel of Fig. 4. Suppressed by the three-body phase space, the sensitivities in the processes are comparatively weaker than those from processes. The most important constraints on the couplings and are from the LFV decays , , and the LEP data Abdallah:2005ph (), and these are depicted as the shaded regions in Fig. 4 Dev:2018upe (). It is also possible to probe the LFV couplings

 |f†eefμτ|,|f†eμfμμ|,|f†eμfμτ|,|f†eμfττ|,|f†eτfμμ|,|f†eτfμτ|,|f†eτfττ| (11)

in the process, however, as a result of the small production cross section, the probable parameter space is very limited for these couplings at both CEPC and ILC.

One can see in Figs. 3 and 4 that in both on-shell and off-shell production modes of the doubly-charged scalar , a large parameter space of the mass and the LFV couplings can be probed at ILC 1 TeV (and CEPC 240 GeV). With a luminosity of order of ab, the couplings can be probed up to the order of in the on-shell mode, and the effective cutoff scale could go up to few 10 TeV in the off-shell processes.

## Iv Conclusion

We have shown that a neutral scalar and doubly-charged scalar , which are both well-motivated in a large class of new physics scenarios, can be probed at future lepton colliders via their LFV couplings to the SM charged lepton sector. We present a model-independent analysis of how far the LFV coupling strengths and the scalar masses can be probed beyond the existing limits from the low-energy sector. In particular, we find that the full mass and coupling range of the neutral scalar, that can explain the muon anomaly, can be tested in the future lepton colliders. This is largely complementary to the searches of LFV in the low-energy experiments and at hadron colliders.

## Acknowledgments

Y.Z. would like to thank the Center for High Energy Physics, Peking University for generous hospitality where the paper was finished. The work of B.D. and Y.Z. is supported by the US Department of Energy under Grant No. DE-SC0017987. The work of R.N.M was supported by the US National Science Foundation Grant No. PHY1620074.

## References

• (1) M. Tanabashi et al. [Particle Data Group], Phys. Rev. D 98, no. 3, 030001 (2018).
• (2) CEPC-SPPC Study Group, IHEP-CEPC-DR-2015-01, IHEP-TH-2015-01, IHEP-EP-2015-01.
• (3) H. Baer et al., arXiv:1306.6352 [hep-ph].
• (4) M. Bicer et al. [TLEP Design Study Working Group], JHEP 1401, 164 (2014) [arXiv:1308.6176 [hep-ex]].
• (5) E. Accomando et al. [CLIC Physics Working Group], hep-ph/0412251.
• (6) P. S. B. Dev, R. N. Mohapatra and Y. Zhang, Phys. Rev. Lett. 120, no. 22, 221804 (2018) [arXiv:1711.08430 [hep-ph]].
• (7) P. S. B. Dev, R. N. Mohapatra and Y. Zhang, Phys. Rev. D 98, no. 7, 075028 (2018) [arXiv:1803.11167 [hep-ph]].
• (8) V. V. Kabachenko and Y. F. Pirogov, Eur. Phys. J. C 4, 525 (1998) [hep-ph/9709414].
• (9) G. C. Cho and H. Shimo, Mod. Phys. A 32, no. 24, 1750127 (2017) [arXiv:1612.07476 [hep-ph]].
• (10) P. M. Ferreira, R. B. Guedes and R. Santos, Phys. Rev. D 75, 055015 (2007) [hep-ph/0611222].
• (11) J. I. Aranda, A. Flores-Tlalpa, F. Ramirez-Zavaleta, F. J. Tlachino, J. J. Toscano and E. S. Tututi, Phys. Rev. D 79, 093009 (2009) [arXiv:0905.4767 [hep-ph]].
• (12) B. Murakami and T. M. P. Tait, Phys. Rev. D 91, 015002 (2015) [arXiv:1410.1485 [hep-ph]].
• (13) S. Banerjee, B. Bhattacherjee, M. Mitra and M. Spannowsky, JHEP 1607, 059 (2016) [arXiv:1603.05952 [hep-ph]].
• (14) I. Chakraborty, A. Datta and A. Kundu, J. Phys. G 43, no. 12, 125001 (2016) [arXiv:1603.06681 [hep-ph]].
• (15) I. Chakraborty, S. Mondal and B. Mukhopadhyaya, Phys. Rev. D 96, no. 11, 115020 (2017) [arXiv:1709.08112 [hep-ph]].
• (16) Q. Qin, Q. Li, C. D. Lü, F. S. Yu and S. H. Zhou, Eur. Phys. J. C 78, no. 10, 835 (2018) [arXiv:1711.07243 [hep-ph]].
• (17) C. Hays, M. Mitra, M. Spannowsky and P. Waite, JHEP 1705, 014 (2017) [arXiv:1701.00870 [hep-ph]].
• (18) W. Rodejohann and H. Zhang, Phys. Rev. D 83, 073005 (2011) [arXiv:1011.3606 [hep-ph]].
• (19) T. Nomura, H. Okada and H. Yokoya, Nucl. Phys. B 929, 193 (2018) [arXiv:1702.03396 [hep-ph]].
• (20) C. S. Aulakh and R. N. Mohapatra, Phys. Lett. 119B, 136 (1982).
• (21) L. J. Hall and M. Suzuki, Nucl. Phys. B 231, 419 (1984).
• (22) G. G. Ross and J. W. F. Valle, Phys. Lett. 151B, 375 (1985).
• (23) R. Barbier et al., Phys. Rept. 420, 1 (2005) [hep-ph/0406039].
• (24) P. S. B. Dev, R. N. Mohapatra and Y. Zhang, JHEP 1605, 174 (2016) [arXiv:1602.05947 [hep-ph]].
• (25) P. S. B Dev, R. N. Mohapatra and Y. Zhang, Phys. Rev. D 95, no. 11, 115001 (2017) [arXiv:1612.09587 [hep-ph]].
• (26) P. S. B. Dev, R. N. Mohapatra and Y. Zhang, Nucl. Phys. B 923, 179 (2017) [arXiv:1703.02471 [hep-ph]].
• (27) A. Maiezza, G. Senjanović and J. C. Vasquez, Phys. Rev. D 95, no. 9, 095004 (2017) [arXiv:1612.09146 [hep-ph]].
• (28) K. S. Babu, Phys. Lett. B 203, 132 (1988).
• (29) P. Q. Hung, Phys. Lett. B 649, 275 (2007) [hep-ph/0612004].
• (30) P. Q. Hung, Phys. Lett. B 659, 585 (2008) [arXiv:0711.0733 [hep-ph]].
• (31) J. P. Bu, Y. Liao and J. Y. Liu, Phys. Lett. B 665, 39 (2008) [arXiv:0802.3241 [hep-ph]].
• (32) C. F. Chang, C. H. V. Chang, C. S. Nugroho and T. C. Yuan, Nucl. Phys. B 910, 293 (2016) [arXiv:1602.00680 [hep-ph]].
• (33) P. Q. Hung, T. Le, V. Q. Tran and T. C. Yuan, Nucl. Phys. B 932, 471 (2018) [arXiv:1701.01761 [hep-ph]].
• (34) G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, Phys. Rept. 516, 1 (2012) [arXiv:1106.0034 [hep-ph]].
• (35) A. Crivellin, J. Heeck and P. Stoffer, Phys. Rev. Lett. 116, no. 8, 081801 (2016) [arXiv:1507.07567 [hep-ph]].
• (36) M. Magg and C. Wetterich, Phys. Lett. B 94, 61 (1980).
• (37) J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980).
• (38) T. P. Cheng and L. -F. Li, Phys. Rev. D 22, 2860 (1980).
• (39) G. Lazarides, Q. Shafi and C. Wetterich, Nucl. Phys. B 181, 287 (1981).
• (40) R. N. Mohapatra and G. Senjanovic, Phys. Rev. D 23, 165 (1981).
• (41) J. C. Pati and A. Salam, Phys. Rev. D 10, 275 (1974) Erratum: [Phys. Rev. D 11, 703 (1975)].
• (42) R. N. Mohapatra and J. C. Pati, Phys. Rev. D 11, 2558 (1975).
• (43) G. Senjanovic and R. N. Mohapatra, Phys. Rev. D 12, 1502 (1975).
• (44) L. Calibbi and G. Signorelli, Riv. Nuovo Cim. 41, no. 2, 1 (2018) [arXiv:1709.00294 [hep-ph]].
• (45) I. F. Ginzburg, G. L. Kotkin, V. G. Serbo and V. I. Telnov, Nucl. Instrum. Meth. 205, 47 (1983).
• (46) I. F. Ginzburg, G. L. Kotkin, S. L. Panfil, V. G. Serbo and V. I. Telnov, Nucl. Instrum. Meth. A 219, 5 (1984).
• (47) V. I. Telnov, Nucl. Instrum. Meth. A 294, 72 (1990).
• (48) L. Willmann et al., Phys. Rev. Lett. 82, 49 (1999) [hep-ex/9807011].
• (49) D. Hanneke, S. Fogwell and G. Gabrielse, Phys. Rev. Lett. 100, 120801 (2008) [arXiv:0801.1134 [physics.atom-ph]].
• (50) J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 45, 589 (2006) [hep-ex/0512012].
• (51) C. Milstene, G. Fisk and A. Para, JINST 1, P10003 (2006) [physics/0609018].
• (52) D. Yu, M. Ruan, V. Boudry and H. Videau, Eur. Phys. J. C 77, no. 9, 591 (2017) [arXiv:1701.07542 [physics.ins-det]].
• (53) A. Hammad, S. Khalil and C. S. Un, Phys. Rev. D 95, no. 5, 055028 (2017) [arXiv:1605.07567 [hep-ph]].
• (54) A. Belyaev, N. D. Christensen and A. Pukhov, Comput. Phys. Commun. 184, 1729 (2013) [arXiv:1207.6082 [hep-ph]].
• (55) The ATLAS collaboration [ATLAS Collaboration], ATLAS-CONF-2017-053.
• (56) CMS Collaboration [CMS Collaboration], CMS-PAS-HIG-16-036.
• (57) G. Abbiendi et al. [OPAL Collaboration], Phys. Lett. B 577, 93 (2003) [hep-ex/0308052].
• (58) P. Achard et al. [L3 Collaboration], Phys. Lett. B 576, 18 (2003) [hep-ex/0309076].
• (59) J. Benesch et al. [MOLLER Collaboration], arXiv:1411.4088 [nucl-ex].
• (60) P. S. B. Dev, M. J. Ramsey-Musolf and Y. Zhang, Phys. Rev. D 98, no. 5, 055013 (2018) [arXiv:1806.08499 [hep-ph]].
• (61) S. Chakrabarti, D. Choudhury, R. M. Godbole and B. Mukhopadhyaya, Phys. Lett. B 434, 347 (1998) [hep-ph/9804297].
• (62) M. Kuze and Y. Sirois, Prog. Part. Nucl. Phys. 50, 1 (2003) Erratum: [Prog. Part. Nucl. Phys. 53, 679 (2004)] [hep-ex/0211048].
• (63) G. Barenboim, K. Huitu, J. Maalampi and M. Raidal, Phys. Lett. B 394, 132 (1997) [hep-ph/9611362].
• (64) M. Lusignoli and S. Petrarca, Phys. Lett. B 226, 397 (1989).
• (65) C. X. Yue and S. Zhao, Eur. Phys. J. C 50, 897 (2007) [hep-ph/0701017].
• (66) S. Godfrey, P. Kalyniak and N. Romanenko, Phys. Rev. D 65, 033009 (2002) [hep-ph/0108258].
• (67) S. Godfrey, P. Kalyniak and N. Romanenko, Phys. Lett. B 545, 361 (2002) [hep-ph/0207240].
• (68) T. G. Rizzo, Phys. Rev. D 25, 1355 (1982) Addendum: [Phys. Rev. D 27, 657 (1983)].
• (69) C. X. Yue, S. Zhao and W. Ma, Nucl. Phys. B 784, 36 (2007) [arXiv:0706.0232 [hep-ph]].
• (70) M. Lindner, M. Platscher and F. S. Queiroz, Phys. Rept. 731, 1 (2018) [arXiv:1610.06587 [hep-ph]].
• (71) P. S. B. Dev and Y. Zhang, JHEP 1810, 199 (2018) [arXiv:1808.00943 [hep-ph]].
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters