Exploring collider aspects of a neutrinophilic Higgs doublet model in multilepton channels

# Exploring collider aspects of a neutrinophilic Higgs doublet model in multilepton channels

Katri Huitu Department of Physics, and Helsinki Institute of Physics, P. O. Box 64, FI-00014 University of Helsinki, Finland    Timo J. Kärkkäinen Department of Physics, and Helsinki Institute of Physics, P. O. Box 64, FI-00014 University of Helsinki, Finland    Subhadeep Mondal Department of Physics, and Helsinki Institute of Physics, P. O. Box 64, FI-00014 University of Helsinki, Finland    Santosh Kumar Rai Regional Centre for Accelerator-based Particle Physics,
Harish-Chandra Research Institute, HBNI, Jhusi, Allahabad - 211019, India
###### Abstract

We consider a neutrinophilic Higgs scenario where the Standard Model is extended by one additional Higgs doublet and three generations of singlet right-handed Majorana neutrinos. Light neutrino masses are generated through mixing with the heavy neutrinos via Type-I seesaw mechanism when the neutrinophilic Higgs gets a vacuum expectation value (VEV). The Dirac neutrino Yukawa coupling in this scenario can be sizable compared to those in the canonical Type-I seesaw mechanism owing to the small neutrinophilic Higgs VEV giving rise to interesting phenomenological consequences. We have explored various signal regions likely to provide a hint of such a scenario at the LHC as well as at future colliders. We have also highlighted the consequences of light neutrino mass hierarchies in collider phenomenology that can complement the findings of neutrino oscillation experiments.

preprint: HIP-2017-35/TH, HRI-RECAPP-2017-015

## I Introduction

The discovery of the 125 GeV Higgs boson Aad et al. (2012); Chatrchyan et al. (2012a) has been a remarkable achievement of the Large Hadron Collider (LHC). This has provided us a closure regarding the predictions of the Standard Model (SM). While our quest towards understanding the physics beyond the Standard Model (BSM) continues, the 13 TeV run of the LHC is expected to make a big impact both in terms of higher energy reach and better precision by accumulating huge amount of data at large luminosity. The enigma of non-zero neutrino mass has pushed the theorists as well as experimentalists to develop new theories and experimental techniques in order to establish the right theoretical pathway towards unveiling the true nature of neutrino mass generation. The neutrino oscillation experiments have established the fact that at least two of the three light neutrinos are massive, and that they have sizable mixing among themselves (for a review, see Gonzalez-Garcia and Maltoni (2008)). The SM, lacking any right-handed neutrinos, is unable to account for these phenomena. This has led to a plethora of scenarios leading to neutrino mass generation Minkowski (1977); Mohapatra and Senjanovic (1980); Gell-Mann et al. (); Yanagida (); Glashow (1980); Schechter and Valle (1982, 1980); Weinberg (1979, 1980). As the resulting neutrino mass eigenstates may be either Dirac or Majorana type, both scenarios have potentially unique signatures Keung and Senjanovic (1983); Datta et al. (1994); Almeida et al. (2000); Panella et al. (2002); Han and Zhang (2006); del Aguila et al. (2007); Huitu et al. (2008); Atre et al. (2009); Keung and Senjanovic (1983); Datta et al. (1994); Almeida et al. (2000); Panella et al. (2002); Han and Zhang (2006); Chen and Dev (2012); Dev et al. (2014); Das et al. (2017) in the collider experiments. The LHC collaborations have put significant effort to extract any possible information about such scenarios from the accumulated data and the null results so far have only been able to constrain the parameter space of various neutrino mass models Abulencia et al. (2007); Chatrchyan et al. (2012b); ATL (2012); Khachatryan et al. (2015); Bhupal Dev et al. (2012); Deppisch et al. (2015).

In the post Higgs discovery LHC era, the true nature of the scalar sector remains another vital area of interest. The natural question that arises is whether the 125 GeV Higgs is the only scalar as predicted by the SM or other exotic scalars exist alongside, as predicted by various BSM theories including some of the neutrino mass models Schechter and Valle (1980); Magg and Wetterich (1980). The measurements of couplings of the 125 GeV Higgs with known SM particles have so far been consistent with the SM predictions CMS (2015). Thus, even if this Higgs boson were indeed part of a larger scalar sector, its mixing with the other states would be small. There are still enough uncertainties in these measurements to allow new exotic scalar multiplets. Unless the LHC observes some hint of a new scalar, our only hope lies in the precision measurements of the Higgs couplings in order to constrain the BSM physics scenarios. Meanwhile, there has been a long term interest in the simplest two-Higgs doublet models (2HDM) (for a review, see Branco et al. (2012)) which are also strongly motivated by supersymmetric scenarios. A two-Higgs doublet model predicts the presence of two CP-even, one CP-odd and two charged Higgses, one of the CP-even Higgs states being the 125 GeV Higgs boson. Despite the presence of these additional scalar states, the mixing between the two doublets can be arranged so that the other scalars are practically decoupled from the SM Higgs. In such cases, the interaction of the SM-like Higgs with the exotic scalars may be so suppressed that any hint of such interactions can be very hard to pick up even with the precision measurements at the LHC. The hope of finding these scalars, therefore, lies in their direct search. While the increasing center-of-mass energy at the LHC can probe heavier exotic particles, extracting any new physics information from the tremendous amount of collected data also faces the increasing challenge of tackling the QCD background. Hence looking for lepton-enriched final states is understandably efficient in suppressing the SM background contributions and probing new physics scenarios which can potentially give rise to lepton-rich final states.

In this work, we consider a 2HDM where the additional Higgs doublet has an odd symmetry charge opposite to all the SM particles, preventing it from interacting directly with the leptons and quarks. One can additionally incorporate right-handed neutrinos in the model with similar transformation property under symmetry as the new Higgs doublet. One can thus generate Dirac neutrino mass terms when the breaks spontaneously and the new Higgs doublet gets a vacuum expectation value (VEV). This class of models, known as neutrinophilic Higgs doublet models (HDM) have been proposed long ago Ma (2001); Gabriel and Nandi (2007); Davidson and Logan (2009) and the relevant phenomenology has been studied quite extensively Gabriel et al. (2008); Haba and Tsumura (2011); Haba and Horita (2011); Chao and Ramsey-Musolf (2014); Maitra et al. (2014); Chakdar et al. (2014); Seto (2015); Bertuzzo et al. (2016); Guo et al. (2017). In principle, one can also generate Majorana neutrino mass terms in such a scenario, since a Majorana mass term for the additional right-handed neutrinos does not break the symmetry but breaks the accidental lepton number symmetry by two units (). Such neutrino mass generation mechanism looks very similar to the Type-I seesaw Minkowski (1977); Mohapatra and Senjanovic (1980); Yanagida (); Gell-Mann et al. () case, save for the fact that one uses the neutrinophilic Higgs VEV instead of electroweak VEV in order to generate the light-heavy neutrino mixing. The advantage of having the additional Higgs doublet to generate non-zero neutrino masses is that the additional VEV can be very small111This is also preferred from naturalness argument ’t Hooft (1980). in order to counter the smallness of the light neutrino masses which would otherwise be fit with a very small Dirac neutrino Yukawa coupling that has no significant collider phenomenological aspects.

Depending on whether the non-zero neutrinos are Dirac or Majorana type, the collider signals of a HDM scenario can be very different. When Majorana neutrinos exist, smoking gun signal would be lepton number violating final states. In this work, instead of looking for direct heavy neutrino production, we have considered the production of the neutrinophilic charged Higgs () and explored its various possible decay modes. There are some earlier studies on the charged Higgs in similar scenarios emphasising its decay into a charged lepton and a heavy neutrino in the process Guo et al. (2017). We show that even cleaner signals can be obtained using this decay mode with higher lepton multiplicity where the SM background is practically non-existent. We also show that sizable signal event rates can be obtained with other possible decay modes of the , which can serve as complementary channels in probing a HDM-like scenario. We perform our analysis using the 13 TeV LHC as well as an collider with 1 TeV center-of-mass energy. In the process, one can extract information on the neutrino sector parameters also. We show that a very clean indication of the neutrino mass hierarchy can be obtained from the multiplicity of the charged leptons in the final state even after a rigorous collider simulation. Such information can be very useful in complementing the neutrino oscillation experiments.

## Ii Model

In the HDM model, the particle content of the SM is extended by one additional Higgs doublet () and three generations of SM gauge singlet right-handed neutrinos (). A discrete symmetry is introduced, under which both and , , are odd while all the SM fields are even. The most general scalar potential involving the two Higgs doublets is given by

 Vsc=−m21ϕ†ϕ+m22ϕ†νϕν−m23(ϕ†ϕν+ϕ†νϕ)+λ12(ϕ†ϕ)2+λ22(ϕ†νϕν)2 +λ3(ϕ†ϕ)(ϕ†νϕν)+λ4(ϕ†ϕν)(ϕ†νϕ)+λ52[(ϕ†ϕν)2+(ϕ†νϕ)2], (1)

where a non-zero explicitly breaks the symmetry in the model. In the absence of this term, the symmetry can be broken spontaneously by a vacuum expectation value (VEV) of the field , while the standard electroweak symmetry is broken when acquires a VEV, .

Let us first discuss a framework, where , i.e. symmetry is broken only spontaneously in order to generate light neutrino masses and mixing. The model is constrained by sterile neutrino searches, effective number of neutrinos and amount of He required in big bang nucleosynthesis (BBN), observed temperature anisotropies of cosmic microwave background (CMB) and astrophysical limits.

Due to an instability of right-handed neutrinos induced by their mixing to left-handed neutrinos, the mixing strength between , , and , that is, , can be probed by sterile neutrino searches. In semileptonic meson decays, are produced, and can subsequently decay to charged leptons and mesons. Present constraints on and allow a region where their magnitude is of order to , assuming GeV Deppisch et al. (2015). For tau-sterile mixing, , assuming GeV.

In HDM, however, we found the model favoring even lower values of active-sterile mixing, of order to , at  GeV, and even lower for higher Majorana neutrino masses (see Fig. 1).

The largest and smallest active-sterile mixings are driven by and elements. Therefore all the active-sterile mixing elements fall between them: . The matrix elements are proportional to , therefore with some constants and . They are deduced from Fig. 1, having values MeV and MeV. The matrix elements then belong to the following interval:

 1.4⋅10−10√GeVMNi≲|Uℓi|≲7.4⋅10−8√GeVMNi (2)

In addition, the constraints for were derived from assumption that the branching ratios for decay are dominant. This is not applicable for HDM, since then the decay modes of right-handed neutrinos are dominated by decays to invisible particles.

As the model is unconstrained by semileptonic and leptonic decay modes, the lower bound for arises from BBN. In the early universe the right-handed neutrinos must be heavy enough to fall off from the thermal equilibrium before BBN. This is due to the latest results for effective number of neutrinos () by PLANCK Planck Collaboration (2016), which forbids large interference from right-handed neutrinos. This leads to a constraint MeV.

In addition neutrinophilic VEV is constrained from both above and below. Ultralight VEV is forbidden by astrophysical constraints: (eV) Zhou (2011); Sher and Triola (2011). On the other hand, the surface energy density associated with the domain wall arising from discrete symmetry breaking is Zeldovich et al. (1974). The effect of these domain walls to the temperature anisotropies of CMB is

 ΔTT≈GηH0 (3)

where is Newton’s gravitational constant, is Hubble constant and we have assumed et al (2006). Since the observed temperature anisotropies by PLANCK are , the birth of a domain wall will not contradict cosmological data if the VEV is small. If we require the contribution to CMB temperature anisotropies not to exceed the experimental limit, together with the astrophysical constraints, we get

 O(eV)≲vν≲O(MeV) (4)

In order to apply perturbative theory to HDM, the absolute values of the elements of the light neutrino Yukawa coupling matrices must be at most. We performed a global fit to available neutrino oscillation data to calculate the matrix elements, assuming normal neutrino mass ordering, higher octant and no CP violation. We found the dependence of the largest Yukawa coupling of and to be

 max|Y(vν,MN)|≈0.629×100keVvν×√MN100GeV (5)

The dependence is illustrated in Fig. 2.

The breaking of symmetry is necessary in order to generate light neutrino masses within the framework of this model by means of their mixing with heavy right-handed neutrinos. One can add the following Yukawa interaction and Majorana neutrino mass terms to the Lagrangian while keeping the parity unbroken:

where, represents the Majorana mass terms corresponding to the right-handed neutrinos. Once acquires a VEV, the Yukawa term gives rise to Dirac neutrino mass terms, .

The physical Higgs sector now consists of two neutral CP-even (, ), one neutral CP-odd () and the charged Higgs ()222We have assumed the scalar potential to be CP invariant.. In the case when , the physical mass eigenvalues at tree level are given by:

 mh=√λ1v2ϕ, mHν=√λ2v2ν, mAν=√−λ5v2, mH±ν=√−v22(λ4+λ5) (7)

where, . being small, terms proportional to have been neglected. Note that the mixing angle between the SM and neutrinophilic Higgs states are proportional to the ratio and can be safely neglected since we assume . Under this circumstance, the CP-even neutrinophilic Higgs () is always light and the heavy neutrino almost always decay into and a light neutrino resulting in an opposite-sign dilepton signal for a charged Higgs pair production channel Maitra et al. (2014). However, if the explicit symmetry breaking term is present in the Lagrangian, i,e, , the mass eigenvalues are given by:

 mh=√λ1v2ϕ, mHν=√m23vϕvν+λ2v2ν, mAν=√m23vϕvν−λ5v2, mH±ν=√m23vϕvν−v22(λ4+λ5) (8)

Now the neutrinophilic CP-even Higgs can be heavy depending on our choice of , . A heavy and (or) opens up the possibility of a cascade decay via heavy neutrinos resulting in multi-lepton signals of such a scenario that we intend to explore. In the limit , the symmetry of the theory is enhanced. Thus, can be assumed to be naturally small. Besides, a large can also give rise to significant mixing between the two Higgs doublets, which is strictly constrained from the present Higgs data.

### ii.1 Neutrino Mass Generation

The neutrino oscillation data Gonzalez-Garcia and Maltoni (2008); Gonzalez-Garcia et al. (2016); de Salas et al. (2017) indicates that at least two of the three light neutrinos have non-zero mass. One of the most natural ways to generate tiny neutrino mass is via seesaw mechanism Minkowski (1977); Mohapatra and Senjanovic (1980); Yanagida (); Gell-Mann et al. (); Glashow (1980); Schechter and Valle (1982, 1980); Weinberg (1979, 1980). In HDM the mechanism is very similar to that of Type-I seesaw Minkowski (1977); Mohapatra and Senjanovic (1980); Yanagida (); Gell-Mann et al. (). The mixing between light and heavy neutrinos is introduced via the term in the aftermath of symmetry breaking, when gets a VEV. In the basis the neutrino mass matrix looks like

 M6×6=(03×3mD3×3mTD3×3mR3×3), (9)

where . The light effective neutrino mass matrix in the approximation is given by

 Mν=mDm−1RmTD. (10)

The above equation looks exactly similar to what we obtain in canonical Type-I seesaw scenario. The only difference is that in the present framework can be quite small and as a result one can have larger compared to the canonical Type-I seesaw scenario, thus making this model phenomenologically more interesting. In order to fit the oscillation data, one also needs to account for the mixing among the three light neutrino states constrained by the PMNS matrix. One can rewrite in Eq. (10) as

 Mν=UTmdiagνU, (11)

where is the diagonal light neutrino mass matrix and is the PMNS mixing matrix. In order to produce proper mixing satisfying the experimental bounds on the PMNS matrix elements, one of the matrices, or , has to be off-diagonal. Here we choose to keep diagonal and fit the PMNS matrix via an off-diagonal . Thus is obtained using Casas-Ibarra parameterization Casas and Ibarra (2001)

 yν=1vν√mdiagRR√mdiagνUT, (12)

where can be any orthogonal matrix and complex provided . For simplicity, we have chosen to be an identity matrix.

Thus with correct choices of the parameters and , Eq. (9) is capable of explaining the neutrino oscillation data at the tree level itself. There is a potential source of large correction Pilaftsis (1992); Grimus and Lavoura (2002) to the neutrino states at one loop arising from the loops. These mass corrections can be sizeable enough to violate the experimental limits. However, the loop contributions to the neutrino masses corresponding to and have a mutual sign difference and can exactly cancel each other if they are mass degenerate Aristizabal Sierra and Yaguna (2011); Aristizabal Sierra (2011, 2012); Haba and Tsumura (2011). As can be seen both from Eq. (7) and (8), the mass splitting between these two states is driven by the parameter which is therefore set equal to zero throughout this work.

## Iii Constraints and Benchmark Points

Constraints on the charged Higgs mass and its couplings may arise from direct collider search results, neutrino oscillation data and lepton flavor violating decay branching ratios. The LHC collaborations have looked for signatures of exotic scalars in various channels and put bounds on the charged Higgs mass in the range 300 - 1000 GeV provided it can decay only into a top and a bottom quark ATL (2016); McCarn (2016); Ohman (2016); Akeroyd et al. (2017). However, in our present scenario, the charged Higgs, being a neutrinophilic one, does not couple to the quarks. In such scenarios, there are no direct search constraints on . In principle, the constraints derived from slepton searches at the LHC can be reinterpreted to put bounds on the neutrinophilic charged Higgs masses although only in the massless limit of the lightest neutralino. Two body decay of the sleptons into a charged lepton and lightest neutralino gives rise to a dilepton signal which can be relevant for the the present scenario. Existing data excludes slepton masses upto 450 GeV in presence of massless neutralino Aad et al. (2014); Collaboration (2017). However, one always obtains same-flavor-opposite-sign(SFOS) lepton pairs from such slepton pair production processes. The signal requirement also demands a jet veto in the central region alongside the SFOS lepton pair for such analyses. In the present scenario, largest event rate in such a signal region can be obtained when decays into a charged lepton and a heavy neutrino. Heavy neutrino further decays into a light neutrino and Z-boson which further decays invisibly. Clearly, the resulting signal cross-section is rendered small due to branching suppressions. Demand of SFOS lepton pairs makes this cross-section even smaller333The obtained signal cross-section for our lightest benchmark point even before the detector simulation is less than the observed number as quoted in Aad et al. (2014); Collaboration (2017).. Thus, the existing slepton mass limit when reinterpreted for proves to be much weaker. Its couplings with the heavy neutrinos on the other hand, can be constrained from neutrino oscillation data and lepton flavor violating decay branching ratios Guo et al. (2017). As mentioned in section II.1, we have used off-diagonal while fitting the PMNS matrix. These off-diagonal entries are severely constrained from LFV decay branching ratio constraints Adam et al. (2013); Baldini et al. (2013); Aubert et al. (2010); Aushev et al. (2010); Hayasaka et al. (2010); Bertl et al. (2006); Bartoszek et al. (2014). These constraints are also reflected upon our choice of the neutrinophilic Higgs VEV, . It has been observed and also verified by us that can be GeV Guo et al. (2017) at the smallest, if the neutrino oscillation data and the LFV constraints are to be satisfied simultaneously, the most stringent constraint arising from the non-observation of BR() Adam et al. (2013); Baldini et al. (2013). This constraint puts the spontaneously breaking scenario in jeopardy. As evident from Fig. 2, such a choice of is clearly ruled out from restrictions on CMB temperature anisotropies induced by domain walls. However, if the symmetry is broken explicitly, this domain wall problem can be averted. Hence for this work, we choose to work with the scenario only.

### iii.1 Charged Higgs branching ratios and pair production cross-section

The possible decay modes of the neutrinophilic charged Higgs () in our present scenario are , , and . The relevant interaction vertices are given in the Appendix. Depending on the mass hierarchy of , () and , and the choice of neutrino mass hierarchy one (or two) of these decay modes determine the event rates of the different possible final states at the collider. Note that the branching ratios of the decays into the neutral CP-even and CP-odd Higgs states are always the same since they are mass degenerate by our choice of the parameters. These two decay modes dominate over the heavy neutrino decay modes always, if the mass difference, is larger than that of the W-boson mass, . This is an artefact of the small Dirac neutrino Yukawa parameters, which are otherwise constrained by neutrino oscillation data and the non-observation of LFV decays. The being smaller by orders of magnitude from the competitive gauge coupling, a large branching ratio into the or decay modes are not ensured even if . In spite of the additional phase space suppression, three-body decays of via off-shell -decay, dominate over these two-body modes unless .

This behavior is depicted in Fig. 3 where the competitive nature of BR() and BR(), where , is clearly visible through the distributions of the starred and circular points respectively. BR() overtakes the three-body decay branching ratio only if GeV. For , BR() takes over and remains the only dominant decay mode.

In Fig. 4, we have shown variation of the production cross-section at the LHC and an collider. The figure on the left shows the variation of the cross-sections as a function of at the 13 TeV LHC and different center-of-mass energies (500 GeV, 1 TeV and 3 TeV) at an collider. Note that, at the LHC, the production channels include , and while for the collider, pair production is the only viable option. Since we have assumed for our study, the cross-sections of the above mentioned second and third production channels exactly equal. Hence we have shown their combined cross-section in the figure and evidently, it dominates over the pair production cross-section throughout the entire charged Higgs mass range. However, both these cross-sections fall rapidly with increasing mass. On the other hand, at an collider the cross-section falls far less rapidly implying the fact that such a collider will be more effective than the LHC in order to probe heavier charged Higgs masses. The figure on the right shows the variation of the pair production cross-section at an collider with varying center-of-mass energies for our chosen benchmark points. Moreover, a lepton collider is likely to be much cleaner in terms of the SM background contributions. In this work, we have taken into account all the aforementioned production channels for LHC and just the pair production for the collider analysis.

### iii.2 Choice of benchmark points

We now proceed to choose some benchmark points representing the different interesting features of the present scenario for further collider studies. As discussed earlier, one can obtain different possible final states depending upon the mass hierarchies of , () and . Since we also aim to correlate the light neutrino mass hierarchy with the multiplicity of different lepton flavor final states, we will study cases in which at least one of the heavy neutrinos is lighter than the neutrinophilic Higgs states so that it can appear in the cascade. In Table 1 below we present the input parameters, relevant masses and the resulting for the four benchmark points of our choice. We have incorporated the complete model in SARAH Staub (2008, 2010, 2011, 2014, 2015), and subsequently imported in SPheno Porod (2003); Porod and Staub (2012) in order to perform the analytical and numerical computation of the masses and mixings of the particles, their branching ratios and other relevant constraints. See Appendix for LFV constraints for our benchmarks.

The four benchmark points are chosen such that all the dominant decay modes of the neutrinophilic Higgs and the heavy neutrinos are highlighted by different mass hierarchies. The relevant branching ratios are shown in Table 2. The two most dominant decay modes of are , where , and . The first decay mode is driven by the Dirac neutrino Yukawa couplings, , whereas the second one is driven by gauge couplings. As discussed above, the elements of are already constrained from the neutrino oscillation data as well as from the LFV constraints, and thus are in general weaker than the competitive gauge coupling. Hence, if the mass splittings among the neutral and charged neutrinophilic Higgs and the heavy neutrino states are such that both and decay modes are kinematically accessible for , the gauge boson associated one becomes its only relevant decay mode. However, if at least one of the heavy neutrinos is lighter than the and the states are almost degenerate to it, then the decay via heavy neutrinos becomes important. The latter scenario is highlighted in BP1 and BP3 while BP2 represents the former scenario. BP4, on the other hand, highlights the situation where the two-body mode competes with the three-body decay into (or ) alongside an off-shell -boson. However, being on the larger side, the three-body decay dominates as discussed earlier in Section III.1. The heavy neutrinos () in this scenario can decay either via the SM gauge bosons () or the different Higgs states. Note that decays of into , and can only occour through their mixing with the light neutrinos which are suppressed in the present scenario. Hence, these decay modes become relevant for only if the neutrinophilic Higgs states are kinematically inaccessible to it. The choice of neutrino mass hierarchy clearly reflects in the branching ratios of both and and is also expected to be reflected in the final event rates of the multi-lepton signals we intend to explore.

In addition, we checked the effect of neutrinophilic Higgses on the oblique parameters . We have ensured the corrections induced by our benchmark points do not exceed the uncertainties given in He et al. (2001). See Fig. 5 for the allowed values.

## Iv Collider Analysis

Charged Higgs in our benchmark scenarios can give rise to novel signatures in lepton enriched final states. Majorana neutrinos, if produced via cascade from the charged Higgs can further decay resulting in same-sign leptonic final states, which are characteristic to seesaw models and also have much less SM background. The gauge bosons resulting from the decays of the neutrinophilic Higgs and heavy neutrinos may also decay leptonically and thus one can easily obtain a multi-lepton final state associated with missing energy. Leptonic branching ratios of the gauge bosons being small, one would expect smaller event rates in the final state with increasing lepton multiplicity. However, it also means less SM background to deal with resulting in cleaner signals. In this section we explore the different possible multi-lepton final states with or without the presence of additional jets along with detailed signal to background simulation in order to ascertain the discovery potential of the charged Higgs for our chosen benchmark points in the context of 13 TeV LHC as well as future lepton colliders.

### iv.1 Identifying signal regions

In the context of LHC, we aim to study cleaner multi-lepton channels with no tagged jets in the final state. The possible final states that we probe in the present context are (SR1), (SR2), and same-sign trilepton () + (SR3), where represents everything else (jets, photons or leptons) 444For SR3, consists of no leptons since in this case, we demand exactly three leptons with same sign in the final state. in the final states. As mentioned earlier, the various branching ratios of and hence the final signal event rates depend on the mass difference factor . Thus, it is interesting to study how the signal rates vary depending on which in turn can also provide an indirect hint about the masses of and (or) .

In Fig. 6 we have shown the variation of the cross-sections corresponding to the three signal regions mentioned above as a function of with color-coded . Note that these cross-sections are theoretical estimates obtained after combining contributions from all the three relevant production modes of at the LHC prior to detector simulation and do not include the cut efficiencies. The two rows of the figures correspond to normal and inverted hierarchy of neutrino masses respectively. While SR1 only receives contribution from pair-production, both SR2 and SR3 are enriched with contributions from pair production as well as associated production of the . Most of the signal events corresponding to SR1 and SR2 are expected to arise from decay into a charged lepton and a heavy neutrino followed by the heavy neutrino decay into a charged lepton and . Depending on the leptonic or hadronic decays of the -bosons, one can obtain various lepton multiplicities as represented by these signal regions. The signal cross-sections are largest when and () are mass degenerate for any given and they drop with increasing and . SR3, on the other hand, receives more contribution when decays into (or ) along with an on-shell or off-shell -boson. Moreover, in the case of pair production, same-sign leptons can not be obtained if both the decays via resulting the cross-section of SR3 being smaller with smaller . However, the associated production channels contribute dominantly to this signal region throughout the whole range of . As evident from Fig. 6, SR3 is the most favorable channel to look for such scenarios. In general, the inverted hierarchy of the light neutrino masses is expected to generate more multi-leptonic events owing to the larger branching ratio, BR() as reflected by the plots on the bottom row.

In a similar way, we now proceed to choose some signal regions for our analyses in the context of an collider. The possible final states we probe in this context are (SR4), (SR5), and (SR6) 555Just like SR3, in SR6 also does not contain any leptons..

The corresponding signal rates are showcased as a function of in Fig. 7. Trends of the distributions are similar to what we obtained for the LHC case. However, the difference in the production cross-section is manifested by the signal cross-sections indicating a larger event rate at the LHC for similar final states at low region. The rapid fall in production cross-section with increasing at the LHC makes it less relevant for heavier charged Higgs masses. An collider can be more effective provided the center-of-mass energy is large enough for the production. Here, the signal rates drop alarmingly close to due to the choice of center-of-mass energy as 1 TeV.

### iv.2 Analysis

In order to carry out the simulation, events were generated at the parton level using MadGraph5 Alwall et al. (2011, 2014) with nn23lo1 parton distribution function Ball et al. (2013, 2015) and the default dynamic factorisation and renormalisation scales https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/FAQ-General-13 (). We have used PYTHIA Sjostrand et al. (2006) for the subsequent decay of the particles, showering and hadronization. After that the events are passed through Delphes de Favereau et al. (2014); Selvaggi (2014); Mertens (2015) for detector simulation. Jets have been reconstructed using anti-kT algorithm via FastJet Cacciari et al. (2012, 2008). The -jet and -jet tagging efficiencies as well as the mistagging efficiencies of the light jets as or -jet have been incorporated according to the latest ATLAS studies in this regard ATL (2015).

#### Primary selection criteria

We have applied the following cuts (C0) on the jets, leptons and photons in order to identify them as final state particles:

• All the charged leptons are selected with a transverse momentum threshold GeV and in the pseudo-rapidity window .

• All the jets including -jets and -jets must have GeV and .

• We demand between all possible pairs of the final state particles to make sure they are well separated.

As discussed in Section III, the choice of neutrino mass hierarchy affects the branching ratios of the neutrinophilic Higgs as well as the heavy neutrinos in certain flavor specific decay modes. Thus, the hierarchical effect is reflected by the abundance of certain flavor of leptons in the signal events. As we have seen, one would expect less abundance of electrons in the final states for normal hierarchy scenario compared to that for inverted hierarchy.

This feature is evident in Fig. 8 which shows the electron multiplicity in the final state with at least four leptons for BP1 in normal as well as inverted hierarchy scenarios. Such lepton multiplicity distributions can thus provide indirect probe of existing neutrino mass hierarchy.

### iv.3 Results@LHC13

In the context of LHC, we have studied the final states corresponding to SR1, SR2 and SR3 as defined in Section IV.1. Although the choice of our signal regions ensure small or no SM background, we have checked the relevant production channels, , , , , , and + jets nevertheless in this regard. In Table 3 we show the expected number of different signal events at 13 TeV run of the LHC with an integrated luminosity () of 1000 after imposing a transverse missing energy cut, GeV and -jet veto (C1) 666These cuts help reduce some of the surviving SM background contributions., in addition to the primary selection criteria, C0. The choice of our signal regions combined with the cuts C1 render the SM backgrounds to negligible event numbers. We have observed that our SR1 is nearly backgroundless whereas SR2 and SR3 are left with 2 and 1 SM-background events, respectively, at 1000 integrated luminosity. As for the obtained signal event numbers, one can easily get an estimate of the expected rate from Fig. 6 for the different final states. However, note that in these figures the heavy neutrino mass is kept fixed at 100 GeV and if this mass is changed, so does the heavy neutrino branching ratios and hence the signal cross-sections. However, the cross-sections shown in these figures are good enough for order of magnitude estimation for a given .

As expected SR1 has the smallest event rate owing to its large lepton multiplicity, but with negligible SM background. Thus it can be a very clean signal but only if the charged Higgs mass is on the lighter side, as in BP1, and at least one of the heavy neutrinos is lighter than the charged Higgs. The situation, however, worsens considerably with increasing charged Higgs mass, as indicated by BP3. SR2 has a much better event rate and can probe BP1 and BP3 at much lower luminosity than SR1. As the numbers in Table 3 indicate, the inverted hierarchy scenario for BP1 can be probed with a 3 statistical significance at an integrated luminosity of 30 in both SR2 and SR3, i.e, if these signal regions are studied, in this mass range can be probed and possibly be excluded with the LHC data already accumulated. For the corresponding normal hierarchy case however, for the same benchmark point, one needs for similar discovery significance in SR3. BP3 requires an integrated luminosity of (for inverted hierarchy) or more. For the benchmark points like BP4 and BP2, the decay is either suppressed or absent altogether. Thus for such points SR1 ceases to be a viable signal region while SR2 is relevant only at large luminosities. In this case SR3 turns out to be the most viable signal region. In this signal region, to achieve 3 statistical significance in the inverted hierarchy case of BP4 and BP2 one requires and respectively. For all the benchmark points, the choice of light neutrino mass hierarchy is clearly manifested through the different signal event rates. Evidently, with multi-leptonic final states, an inverted hierarchy scenario is more likely to be probed at lower luminosities at the LHC.

### iv.4 Results@e+e− collider

In Table 4 below we have presented the expected number of different signal events at 1 TeV run of an collider with an integrated luminosity of 100 after imposing a missing energy cut, GeV and b-jet veto (D1), in addition to the primary selection criteria, C0. The event rates are quite good and devoid of any direct SM background, which makes it an ideal platform to look for a neutrinophilic charged Higgs. Although the number of events shown in Table 4 correspond to 100 , the inverted hierarchy scenarios in BP1 and BP3 can be probed with a statistical significance of 3 at a much lower luminosity ( 10 ). Note the improved event rates in signal regions SR5 and SR6 despite of the smaller pair production cross-section in BP3 over those of BP1. This is a consequence of increased hadronic branching ratio of and improved cut efficiency due to the larger mass gap between and in BP3. Even BP2 which can be probed at the LHC only at very high luminosity can be probed here at around with similar statistical significance via SR5 which turns out to be the most favored signal in general for all the benchmark points. The overall signal rate is relatively weaker in SR6 due to better lepton tagging efficiency at a lepton collider, which results in smaller number of events with exactly three same-sign leptons as demanded.

## V Summary and Conclusions

We have considered a simple extension of the SM with one additional scalar doublet and three generations of singlet right-handed Majorana neutrinos, where the additional Higgs states interact with the SM sector only via the right-handed neutrinos. The model, known as the neutrinophilic Higgs doublet model, is a well-motivated framework from the viewpoint of neutrino mass generation. The light neutrinos gain tiny non-zero masses via Type-I seesaw mechanism when the neutrinophilic Higgs obtains a VEV to break the symmetry. We have discussed in brief why the spontaneous breaking of the symmetry is disfavored, if one imposes the constraints derived from the CMB temperature anisotropies induced by domain walls as well as LFV decay branching ratios. We have, therefore, considered a scenario where the parity is broken explicitly and thus is devoid of the domain wall problem. In such a scenario, the charged Higgs can have interesting collider phenomenology, explored in this work. Depending on the different decay modes of the neutrinophilic charged Higgs, we have identified some particularly clean signal regions likely to provide a hint of HDM scenarios at the collider experiments.

We have also highlighted the interesting roleplay of the light neutrino mass hierarchy. Whether the neutrinos follow normal or inverted hierarchy, is likely to be manifested via multiplicity of different flavored leptons in the final state. Thus such a finding at the collider experiments can complement the neutrino oscillation experiments which are yet to ascertain the correct mass hierarchy of the three light neutrinos.

The fact that the charged Higgs pair production cross-section falls quite rapidly at the LHC with increasing mass, led us to perform a comparative study between the LHC and a future machine in order to probe such scenarios. We observed that although LHC is quite efficient to probe light charged Higgs masses, an collider will be able to probe a much larger parameter space with heavier states.

## Vi Acknowledgements

KH and SM acknowledge the H2020-MSCA-RICE-2014 grant no. 645722 (NonMinimalHiggs). TK expresses his gratitude to Magnus Ehrnrooth foundation for financial support. The work of SKR was partially supported by funding available from the Department of Atomic Energy, Government of India, for the Regional Centre for Accelerator-based Particle Physics (RECAPP), Harish-Chandra Research Institute.

## Vii Appendix

### vii.1 Relevant interaction vertices

 ℓ∓iNjH±kν i3∑a=13∑b=1yabνUiaLUj3+bN HiνH±jνW∓ i2gZi2(pHiν−pH±jν)μ νiNjHkν − i3∑a=13∑b=1UiaNUj3+bNyabνZk2

where , and are the charged lepton, neutrino and CP-even neutral Higgs mixing matrices, the bases of the mass matrices being , and respectively. Note that, is a diagonal matrix.

### vii.2 Lepton Flavor Violating Branching Ratios

In Table 5 we have shown the obtained branching ratios for various lepton flavor violating processes corresponding to the four benchmark points. BR() is projected to be probed experimentally up to in near future Baldini et al. (2013). As indicated by the numbers, the obtained branching ratios for this process are at least one order of magnitude smaller for our benchmark points. The rest of these obtained LFV branching ratios are several orders of magnitude below the present experimental sensitivity in the respective channels.

## References

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