Non-Standard Neutrino Interactions in Modified \nu2HDM

# Non-Standard Neutrino Interactions in Modified ν2hdm

Ujjal Kumar Dey, b,c    Newton Nath, d    Soumya Sadhukhan Centre for Theoretical Studies, Indian Institute of Technology Kharagpur,
Kharagpur 721302, IndiaInstitute of High Energy Physics, Chinese Academy of Sciences,
Beijing 100049, ChinaSchool of Physical Sciences, University of Chinese Academy of Sciences,
Beijing 100049, ChinaPhysical Research Laboratory,
###### Abstract

In the traditional neutrinophilic two-Higgs doublet model (2HDM), there is no non-standard neutrino interaction (NSI) as the interactions between the Standard Model fermions with neutrinos are negligibly small due to the tiny mixing of the two scalar doublets. In this work, we propose that if 2HDM is modified by considering the right-handed electron, is negatively charged under a global -symmetry then one can generate significant amount of NSI along with the tiny Dirac neutrino mass. Depending on different constraints from the LEP experiment, tree level lepton flavor violating processes, big-bang neucleosynthesis etc., we observe that this model significantly restricts the range of permissible NSI parameters, putting a strict upper bound on different NSIs. Furthermore, considering these model-dependent NSIs, we study their impact on the next-generation superbeam experiment, DUNE. We present a detailed discussion on the mass hierarchy sensitivity and the CP-violation discovery study considering the impact of both diagonal as well as off-diagonal NSIs.

## 1 Introduction

The Standard Model (SM) of particle physics epitomizes our knowledge of fundamental interactions of subatomic world with all its grandeur. Barring from a few minor disagreements the SM is hitherto untarnished by the direct observations from high energy colliders like LHC. The most elusive of the particles that constitute the SM is undoubtedly the neutrinos which can very well usher us towards the physics beyond the SM owing to its perplexing properties that they are massive (though tiny) and their different flavors are substantially mixed OHLSSON20161 . These distinct attributes are exotic to the basic tenets of SM but are experimentally observed. The discovery of neutrino oscillation indubitably established the existence of mass for the neutrinos, whereas in the SM the neutrinos are supposed to be massless. The Super-Kamiokande experiment in Japan published their result to establish the phenomenon of neutrino oscillation in 1998 Fukuda:1998mi . Later, various other experiments such as solar, reactor and most recently long baseline experiments like, T2K Abe:2017uxa and NOA Adamson:2017gxd confirm the phenomenon of neutrino oscillation. In addition to the confirmation of the presence of neutrino mass, neutrino oscillation results bring forth some pertinent questions, e.g., (i) the neutrino mass hierarchy, i.e. whether neutrinos obey normal hierarchy (NH, ) or inverted hierarchy (IH, ), (ii) the octant of , i.e. lies in the lower octant (LO, ) or in the higher octant (HO, ) and (iii) the determination of Dirac CP phase, . As the current experiments are striving to resolve the degeneracies in hierarchy and octants and to assess the sensitivity in CP violation, new theoretical inputs which can possibly help to address these issues are worth exploring. It is thus quite interesting to reconnoiter models which can provide natural explanation to the tiny neutrino mass along with other issues mentioned above.

The enigma of neutrino mass lead the particle physicists to explore various theoretical models which can explain the neutrino mass as well as the observed neutrino mixing angles. Over the period the literature grew (and still growing) on neutrino mass and today there are a number of phenomenological models addressing this issue (for a review see King:2015aea ; Cai:2017jrq ). Some of the new physics scenario could lead to corrections to the effective neutrino interactions through higher dimensional operator. Non-standard interactions (NSI) of neutrinos can be induced by the new physics beyond the SM (BSM). In the BSM scenarios NSI can arise when the heavier messenger fields are integrated out which can generate the dimension-6 Buchmuller:1985jz ; Bergmann:1998ft ; Bergmann:1999pk and dimension-8 Berezhiani:2001rs ; Davidson:2003ha effective operators. For a detailed review and phenomenological consequences see Refs. Ohlsson:2012kf ; Miranda:2015dra and the references therein. These were originally discussed even before the establishment of neutrino oscillation phenomena Wolfenstein:1977ue ; Valle:1987gv ; Roulet:1991sm ; Guzzo:1991hi . Over the year many BSM scenarios have been studied where NSIs can be realized. Some of the popular models where this can be present are the extended models with new particle as the messenger, models with single or multiply charged heavy scalars, leptoquark, -parity violating supersymmetry etc. Farzan:2017xzy ; Bilenky:1993bt ; Barranco:2007tz ; Antusch:2008tz ; Malinsky:2008qn ; Ohlsson:2009vk ; Forero:2016ghr . Moreover, there are studies which considers the effect of NSI in collider experiments Anchordoqui:2012qu ; Franzosi:2015wha ; Davidson:2011kr ; Davidson:2011xz ; Friedland:2011za ; Choudhury:2018xsm . Normally in the extensions of SM which gives rise to NSIs can get stringent constraints from the charged lepton flavor violating (LFV) processes. Therefore, it is imperative to maintain these LFV constraints while building any model from the NSI perspective.

We construct a variant of neutrinophilic two-Higgs doublet model (2HDM) where we can have a significant NSI while maintaining various phenomenological constraints. In the standard 2HDM there is no NSI due to the absence of interaction between left-handed neutrinos with the other leptons and quarks via the extra scalars. We modify the standard 2HDM such that the second scalar doublet couples only to the electron and neutrinos. This is achieved by assigning a negative charge of the under a global symmetry. The charged Higgs present in the spectrum couples the left-handed electron and the charged leptons and thus can give rise to NSI by playing the role of the messenger. In this set-up we can have different NSI parameters. Apart from the presence of NSI the hierarchy in the fermion sector is somewhat less drastic in this model compared to the SM. This is because of the fact that the second Higgs doublet couples only to the electron and neutrinos and thus it is responsible for the mass of neutrinos and electron. The first Higgs doublet takes care of all the other SM fermions. Clearly, the hierarchy in the Yukawa couplings is minimized owing to the two vacuum expectation values (vev) of the two doublets. This betterment in the fermion hierarchy is one of the novel features of this model. Phenomenological constraints characteristic to any 2HDM will also be applicable in this scenario too. We consider the constraints on the model parameter space from lepton flavor violating processes (LFV), oblique parameters, , big-bang neucleosynthesis (BBN) etc.

Furthermore, we study the impact of model-dependent NSIs in case of the long baseline (LBL) experiment like, DUNE. For DUNE, the mass hierarchy sensitivity considering model-independent bounds on NSI has been studied in Liao:2016hsa ; Coloma:2016gei ; Deepthi:2016erc ; Deepthi:2017gxg . These studies show that the mass hierarchy sensitivity of LBL experiments get seriously compromised due to the presence of intrinsic , degeneracy. This degeneracy remains true irrespective of baseline and energy. On the other hand, if LBL experiment, like DUNE observes sensitivity then it will be able to rule out certain parameter space of model-independent NSIs. In the modified 2HDM model considered here, parameter space of different NSIs are constrained which helps to avert mass hierarchy degeneracy. We further examine the mass hierarchy sensitivity of DUNE considering model-based NSI parameters. Later, we also address the issue of CP-violation (CPV) sensitivity for DUNE. We illustrate the role of new CP-phase on the measurement of considering both CP-conserving as well as CP violating values for the new phases. In our analysis, we assume NSI both in data as well as in theory. Some of the recent phenomenological studies, considering model-independent bounds of NSI, in the context of DUNE can be found in Masud:2015xva ; deGouvea:2015ndi ; Coloma:2015kiu ; Liao:2016hsa ; Soumya:2016enw ; Blennow:2016etl ; Forero:2016cmb ; Huitu:2016bmb ; Masud:2016bvp ; Coloma:2016gei ; Masud:2016gcl ; Agarwalla:2016fkh ; Liao:2016bgf ; Fukasawa:2016gvm ; Blennow:2016jkn ; Liao:2016orc ; Deepthi:2016erc ; Ghosh:2017ged ; Ghosh:2017lim ; Deepthi:2017gxg and the references therein.

We organize our paper as follows. In section 2, we introduce the concept of non-standard neutrino interaction and their model-independent bounds. Section 3 is devoted to a detailed description of the model that has been considered in this work. Traditional 2HDM model is discussed in the first part whereas the modification of the model which leads to NSI is discussed in the second part. Various phenomenological constraints coming from LEP data, LFV processes, BBN etc. are discussed in section 4. We illustrate the effect of diagonal and off-diagonal NSI for DUNE in section 5. First, we present our result considering appearance channel probability (). Later, we discuss the impact of NSI on the determination of mass hierarchy sensitivity and the CPV sensitivity in case of DUNE. Finally, we summarize and conclude our noteworthy findings in section 6.

## 2 General Description of NSI

In this study, we consider the effect of neutral-current NSI in presence of matter which is describe by the dimension-6 four- fermion operators of the form Wolfenstein:1977ue ,

 LNCNSI=−2√2GF(¯¯¯ναγρPLνβ)(¯fγρPCf)ϵfCαβ+h.c. (1)

where are NSI parameters, , , , and is the Fermi constant111Note that here we neglect NSI due to charges-current interactions which mainly affect neutrino production and detection Khan:2013hva ; Ohlsson:2013nna ; Girardi:2014kca ; DiIura:2014csa ; Agarwalla:2014bsa ; Blennow:2015nxa .. Note that in general, the elements of are complex for and real for due to the hermitian nature of the interactions. For the matter NSI, are defined as,

 ϵαβ=∑f,CϵfCαβNfNe , (2)

where is the number density of fermion and . For the Earth matter, we can assume that the number densities of electrons, protons, and neutrons are equal (i.e. ), in such a case and one can write,

 ϵαβ=√∑C((ϵeCαβ)2+(3ϵuCαβ)2+(3ϵdCαβ)2). (3)

The modified Hamiltonian in presence of propagation NSI, in the flavor basis, can be written as,

 H=12E[Udiag(0,Δm221,Δm231)U†+diag(A,0,0)+Aϵαβ], (4)

where is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix Patrignani:2016xqp , (), represents the potential arising from standard interactions (SI) of neutrinos in matter, and can be written as

 ϵαβ =⎛⎜⎝ϵeeϵeμϵeτϵ∗eμϵμμϵμτϵ∗eτϵ∗μτϵττ⎞⎟⎠, (5)

where for .

The model-independent bounds Blennow:2014sja ; Ohlsson:2012kf on these parameters are

 |ϵee|<4.2,  |ϵeμ|<0.33,  |ϵeτ|<3.0, |ϵμμ|<0.07,  |ϵμτ|<0.33,  ϵττ|<21. (6)

Having introduced general descriptions of NSI and its model-independent bounds, in next section we describe our model and calculate the model-dependent bounds of different NSI parameters persistence with different experimental bounds.

## 3 Model Description

For completeness we first lay down the details of the traditional neutrinophilic two Higgs doublet model (2HDM), for more details see e.g., Davidson:2009ha ; Machado:2015sha and references therein. In passing we note the reasons for the absence of NSI in the standard scenario. We then modify the vanilla 2HDM which gives rise to NSI effects.

### 3.1 Standard ν2hdm

We start with an extension of the SM by adding an extra scalar doublet having similar gauge quantum numbers as the SM scalar doublet . Three right handed neutrinos (RHN), are introduced and they will constitute the Dirac masses with the left handed neutrinos of the SM. There are mainly two ways to obtain the masses for the neutrinos with this particle spectrum. In one set-up Ma:2000cc ; Gabriel:2006ns a symmetry is imposed under which the fields and are odd and all the SM fields are even. In this set-up the Majorana mass terms for are not forbidden. In another set-up Davidson:2009ha , a global symmetry is imposed and the fields and have +1 charge under this and all the SM fields are neutral. In this scenario the neutrinos exclusively have Dirac masses. For the reasons that will be discussed in the later part we will follow the case closely. In both of these cases the Yukawa interaction of the neutrinos is of the form where is the left-handed lepton doublet and . The Yukawa structure of the other SM fermions is of the usual form involving . Note that when the is unbroken, has no vev and the neutrinos remain massless.

The most general scalar potential for the exact symmetric222If instead of the global one imposes then another quartic term of the form has to be added to . case is given by

 V(Φ1,Φ2) =m211Φ†1Φ1+m222Φ†2Φ2 +λ12(Φ†1Φ1)2+λ22(Φ†2Φ2)2+λ3Φ†1Φ1Φ†2Φ2+λ4Φ†1Φ2Φ†2Φ1. (7)

The complex scalar doublets and carry hypercharge . Since in the exact symmetric case neutrinos are massless, to give masses to neutrinos this symmetry is to be broken. This is done by introducing a soft-breaking term of the form . Clearly, the smallness of makes it technically natural in the ’t-Hooft sense. The two scalar doublets can be presented as

 Φa=(ϕ+a(va+ha+iηa)/√2),a=1,2. (8)

The vacuum expectation values (vev) of the two scalars can be denoted as , , and responsible for neutrino masses. Therefore, while . In general the physical scalar fields are given by

 H+=ϕ+1sinβ−ϕ+2cosβ≃−ϕ+2,A=η1sinβ−η2cosβ≃−η2, (9a) h=−h1cosα−h2sinα≃−h1,H=h1sinα−h2cosα≃−h2, (9b)

where is the rotation angle for CP-even states and is of the CP-odd and charged scalar states. These rotations diagonalize the mass matrices with,

 tan2α =2(−m212+λ34v1v2)m212(v2/v1−v1/v2)+λ1v21−λ2v22, (10) tanβ =v2v1, (11)

where . From the requirement of tiny neutrino mass we have . Thus for small , both and will be very small. From Eqs. (9) it is evident that the SM-like 125 GeV Higgs comes from the doublet whereas the BSM scalars are dominantly comprised of components of the doublet . The BSM scalars develop neutrinophilic interactions in the Yukawa sector. The Yukawa couplings of the new scalars in the limit are described as,

 LY⊃mνiv2H¯νiνi−imνiv2A¯νiγ5νi−√2mνiv2[U∗ℓiH+¯νiPLℓ+h.c.], (12)

where are neutrino masses and is the PMNS matrix.

In this set-up due to the above-mentioned charge assignments of the respective fields, left-handed neutrinos can not couple to right-handed charged leptons through the scalar doublet . The possibility of such coupling via the mixing of the other doublet is negligible since that coupling is proportional to which is negligibly small. Since this type of coupling is needed to generate the NSI effect, in this set-up the NSI effect will be negligible.

### 3.2 Modified ν2hdm

The standard 2HDM scenario discussed above can not give rise to non-standard interactions in the neutrino sector. Since two different Higgs doublets are responsible for providing mass to the neutrinos and other massive SM particles, there is no interaction of left-handed neutrinos with the SM massive leptons and quarks. Here we propose to give mass to the electron along with the neutrinos through the second scalar doublet . For that purpose, among the SM fields only the right-handed electron is endowed with the charge () under the global . With this quantum number assignment the relevant terms for the neutrino and electron sector become

 Lmν2HDM⊃ye¯LeΦ2eR+yν¯Le~Φ2νR+h.c., (13)

where is the SM electron doublet . Expanding the first term of the Lagrangian gives us terms involving charged Higgs () as,

 LYukH±⊃ye¯νeLH+eR+h.c. (14)

Presence of this interaction results in the -channel process of Fig. 1 which can give rise to non-standard interaction terms between SM neutrinos and electrons. The effective Lagrangian, after integrating out the heavy charged Higgs, looks like

 Leff⊃y2e4m2H±(¯νeLγρνeL) (¯eRγρeR)+h.c. (15)

From this effective Lagrangian, comparing it with the defined form of the NSI Lagrangian, Eq. (1), the NSI parameter can be written as

 ϵee=12√2GFy2e4m2H±. (16)

Note that, in this model, only the right handed electron contributes in Eq. (3) to provide . With the same quantum number assignment, the extra relevant terms for the lepton sector that can be added are:

 Lmν2HDM⊃y1¯LμΦ2eR+y2¯LτΦ2eR+h.c., (17)

where are the SM lepton doublets . These terms will provide the following NSI parameters,

 ϵeμ=12√2GFyey14m2H±,ϵeτ=12√2GFyey24m2H±, (18a) ϵμτ=12√2GFy1y24m2H±,ϵμμ=12√2GFy214m2H±. (18b)

Note that here the Yukawa couplings can be complex. The complex nature of these Yukawa results in generating phases for the NSI parameters.

Before delving into the discussion on phenomenological constraints on this model and the NSI studies we mention a few features of this model. The first term of Eq. (13) appears due to modification in the 2HDM. The vacuum expectation value of second doublet gives mass to the electron here. Thus with an order one Yukawa coupling to get the electron mass  MeV we need MeV. With such a large  MeV the neutrino Yukawas will be of the order of . This re-introduces hierarchy in the fermion Yukawa couplings i.e. the hierarchy of the Yukawas to accommodate mass of different fermions ranging from to , with the same vev. In the SM if the neutrinos are to be given mass through the SM Higgs mechanism, a hierarchy in the Yukawas is required to accommodate the top mass () along with the Dirac neutrino () mass, using a vev of  GeV. In the general 2HDM since the neutrinos are given mass through the second doublet, a hierarchy is only required to arrange for masses from top quark to electron. In our modified 2HDM, we have a hierarchy in Yukawa couplings to account for the mass of the electron and neutrinos using the same  MeV. There remains another hierarchy to arrange for the masses of the rest of the fermions, i.e., from top quark to the muon, using the vev  GeV. Distributing the Yukawa hierarchies into two sectors, i.e. in sector and in sector, reduces the overall hierarchy of Yukawas in the model in comparison to the SM (where it is ). However, the modified 2HDM has two different hierarchies, which is somewhat less likable feature than the requirement of only one hierarchy of order in the original 2HDM.

In this modified 2HDM, leptophilic BSM scalar couplings arise where () couples through the Lagrangian terms

 LLP⊃ye√2HeLeR+iye√2AeLeR+h.c. (19)

This coupling has stringent bounds from LEP data which we discuss in the next section.

Before ending this section we make a few remarks about the choice of -symmetric 2HDM over its symmetric counterpart. Firstly, in the symmetric 2HDM CP-even new scalar is very light. For such a light scalar of mass up to MeV scale the constraint on the leptophilic coupling is very tight. If 2HDM with is modified to have an electrophilic scalar coupling, from Knapen:2017xzo , the upper limit on the Yukawa will be of the order of . With a tiny Yukawa coupling which also determines the magnitude of the non-standard neutrino interaction, NSI effects of the modified model will be negligible. Secondly, the symmetric 2HDM has been severely constrained from the oblique parameters. The modification of the oblique parameter due to the presence of a very light neutral scalar, , pushes it to a large negative value to rule out the model at confidence level333It has been shown that the introduction of vector-like leptons can relax this constraint Sadhukhan:2018nsk . Machado:2015sha . In the case, however, sufficiently heavy CP-even scalar is possible and also the mass degeneracy of CP-even and CP-odd scalars helps in satisfying parameter bounds easily.

## 4 Phenomenological Constraints

Typical to any 2HDM scenario the parameter space of the present set-up will also be subject to constraints from the data pertaining to heavy scalars, flavor violating issues etc. Note that in our model the second doublet couples only to the leptons and thus the constrains from the LEP and LFV is of more relevance. In this section we discuss constraints on the heavy scalars from the LEP searches, LFV decays, as well as BBN.

### 4.1 Constraints from LEP

#### 4.1.1 Charged scalar mass

The boson decay width measurement at LEP points out that there is a tiny room for decays to BSM particles. This suggests that the new scalars of the modified 2HDM should be heavier than half of boson mass to kinematically forbid the decays. The search of a charged Higgs at LEP through the process with further decaying to puts a lower bound on the charged Higgs mass as  GeV Abbiendi:2013hk .

#### 4.1.2 Constraint from e+e−→l+l−

Measurement of cross section at the LEP experiment can be expressed in terms of a limit on the scale of an effective four-electron interaction as  TeV LEP:2003aa . In modified 2HDM, process will take place through both the CP-even () and CP-odd () scalar mediators and the effective four lepton operator will look like

 Leff⊃ye28m2H(¯eLγρeL)(¯eRγρeR)+ye28m2A(¯eLγρeL)(¯eRγρeR). (20)

As in the case of a global -symmetric 2HDM, the degeneracy of CP-even and CP-odd scalar masses reduce the effective operator as,

 Leff⊃ye24m2H(¯eLγρeL)(¯eRγρeR),

which being compared to the effective operator form with a scale provide the bound on as In our modified 2HDM with global symmetry, the constraints coming from oblique parameters () allow mass splitting between the charged scalars and the neutral ones to be largest when is small ( GeV) Machado:2015sha . For a light charged Higgs the neutral scalar masses can go up to several hundred GeV’s. We take  GeV. If higher values of charged Higgs masses are taken then mass splitting between charged and neutral scalars decreases, which along with large constrains the Yukawa coupling tightly enough to allow very small values of NSI parameters. With  GeV and  TeV the constraint translates to which further translates to a limit on the vev of as

 v2≥2.5 MeV. (21)

Using this bound we fix the value at MeV to get tightest limits on other Yukawas from LEP measurement in other processes like, .

#### 4.1.3 Mono-photon constraint

Another LEP constraint arises from the dark matter search in the mono-photon signal , where is either initial state radiation or it appears due to internal bremsstrahlung. In modified 2HDM similar mono-photon processes can occur, through the charged Higgs exchange which couples right-handed electron with left-handed neutrinos and vice versa. The LEP DM search limit is thus rewritten as

 y4e+2y2ey22+y42≤16m4H±Λ4DM, (22)

with DM scale  GeV Fox:2011fx for light DM particles.

### 4.2 Lepton Flavor Violation Constraints

#### 4.2.1 ℓα→ℓβγ

For a generic process , the scalar mediated branching ratio reads Fukuyama:2008sz

 BR(ℓα→ℓβγ)=BR(ℓα→e¯νν)αEM192π|⟨m2αβ⟩|2ρ2. (23)

The strongest bound on this type of lepton flavor violating decay comes from the MEG experiment which gives the upper limit as  Adam:2013mnn , while bounds on the other LFV decay channels (), obtained from the BaBar Collaboration, are weaker. Though the lepton flavor violating processes like happens through the loop driven processes in this model, experimental constraints are strong enough to severely constrain that. The decay branching ratio for this process which occurs through the charged scalar mediator can be written as Bertuzzo:2015ada

 BR(μ→eγ)=BR(μ→e¯νν)αEM192π|⟨m2μe⟩|2ρ2, (24)

where . In terms of defined here, the confidence level bound read:

 ρ≲1.2eV−2     [μ→eγ], (25)

The limit on translates into a limit where for eV one must have GeV. This is the tightest limit upto now on the and parameter space. Limits from similar other LFV decay modes are less constraining. With the sensitivity of the MEG expected to be improved further, the bound on is expected to be improved by about one order of magnitude. The limits imposed by the MEG bound on the plane are shown in Fig. 1 of Bertuzzo:2015ada .

#### 4.2.2 τ(μ)→3e

The lepton flavor violating decays like, and is also possible in the modified 2HDM through the neutral scalar () mediation at tree level and thus putting stringent constraints on the Yukawa couplings and . Belle results Hayasaka:2010np of decay normalized to decay as

 Γ(τ→3e)Γ(τ→μνμντ)≤1.58×10−7,

which implies the constraints on 2HDM parameter space as

 yey2≤(0.16 mH)2(1TeV)2. (26)

This is a tighter bound on - plane compared to the LEP and LEP mono-photon constraints. Similarly for decay there is even stronger experimental limit as will put bound on - plane as

 yey1≤(8.12×10−3 mH)2(1TeV)2. (27)

For moderate and allowed values this bound reduces the allowed Yukawa values to a tiny value () which makes any NSI parameter involving like insignificant. We have not considered this Yukawa coupling for our further analysis.

In Fig. 2 we show the parameter space allowed by various constraints from LEP and rare lepton decays. Here the neutral scalar masses are kept at  GeV where the charged Higgs mass is around  GeV. As discussed earlier, maximal mass difference between neutral and charged scalars is required to enhance NSI parameter values. This is the amount of maximal mass splitting that is allowed from the parameters.

### 4.3 g−2 Constraints

The charged scalars can contribute to muon and electron  Dedes:2001nx at one loop, but the contribution is negligible due to a suppression factor in the amplitude. Unlike a general 2HDM where the two loop contributions are dominant, charged lepton couplings to are suppressed here by a factor in 2HDM, making constraints insignificant in this scenario, as shown in Grifols:1979yk .

### 4.4 BBN Constraints

The new particles that are introduced in this model are the right-handed neutrinos which are very light, with mass of the scale of eV. These new degrees of freedom, ultra light neutrinos can be created in the early Universe to populate it through the process that occur through the exchange of the charged scalar . From the big bang neucleosynthesis (BBN) there is a limit on new relativistic degrees of freedom which stands as at 95 confidence level with He + Planck TT + lowP + BAO data Ade:2015xua . From this one can obtain the decoupling temperatures of neutrinos which for right-handed neutrinos is  MeV and for left-handed ones  MeV Olive:1999ij ; Anchordoqui:2012qu ; Zhang:2015wua . These considerations put an experimental upper limit on the charged Higgs mediated process as  Davidson:2009ha which can be translated to a bound on the neutrino Yukawa coupling as

 yνi≤0.05×[mH±100 % GeV][1/√2|Uei|]. (28)

In the usual 2HDM this constraint is used to put tight constraints on Yukawas. In the modified 2HDM this constraint is trivially satisfied due to the larger values of  MeV the preferred values are much smaller.

## 5 Study of NSI

In this section, we first present a discussion on the constraint range of NSI parameters in this model. Furthermore, we proceed with a detailed study of these bounds considering DUNE where our initial discussion is on the impact of model-dependent NSI considering appearance channel probability for both neutrinos and antineutrinos. Moreover, we also perform substantially detailed study at the level where we illustrate the impact of NSI on mass hierarchy sensitivity and CP violation discovery study.

### 5.1 NSI in the modified ν2hdm

Here we fix different benchmark values of and vary in the range of  GeV. The LEP bound sets the lower limit of charged Higgs mass at 80 GeV. To satisfy the LEP constraint we need to have heavier electrically neutral BSM scalars as that can keep the cross section low. On the other hand, the heavier charged Higgs will suppress the values of non-standard neutrino interaction parameters, as can be seen from Eq. (16). Therefore, it is optimal to have a maximal splitting between charged and neutral BSM scalars and the splitting is controlled by oblique parameter constraints. Maximal splitting is possible for smaller mass of charged scalar upto  GeV, when the neutral scalar masses are at  GeV. So, the upper limit of is chosen at  GeV. From the LEP upper limit on process, a lower bound on value is given as  MeV. For three different values of  MeV we choose three benchmark points and calculate for those values. Allowed values of parameter for different benchmark points are given in Table 1.

Constraints from lepton flavor violating decay forces the Yukawa coupling to be small, i.e. of the order of so that all the values involving i.e. will be negligible to have any significant effect on the observables. So we do not consider those NSI parameters for our analysis. Only the NSI parameters that can be explored are the ones like . In general, the model independent constraint on NSI parameter is too loose to put any significant constraint on model parameters. For the case of , the tightest bound comes from LFV decay which restricts to be below for  MeV and this upper limit increases with . We have chosen different values of Yukawa coupling within this limit. For these different values, numbers at three benchmark points are given in Table 2.

### 5.2 Probability level descriptions

We now study the appearance channel () probability in case of DUNE considering model-dependent NSI parameters, and . The relevant analytical expression for the appearance channel probability in case of normal hierarchy (NH) considering , and as a small parameters, except , can be written as Liao:2016hsa

 Pμe = x2f2+2xyfgcos(Δ+δCP)+y2g2 (29) + 4^Aϵeτs23c23{xf[fcos(ϕeτ+δ)−gcos(Δ+δ+ϕeτ)]−yg[gcosϕeτ−fcos(Δ−ϕeτ)]} + 4^A2(g2+f2)c223s223|ϵeτ|2−8^A2fgs23c23c23ϵ2eτcosΔ + O(s213ϵ,s13ϵ2,ϵ3),

where

 x=2s13s23, y=2rs12c12c23, (sij=sinθij,cij=cosθij,ij=12,13,23), Δ = Δm231L4E, ^A=AΔm231, f,¯f=sin[Δ(1∓^A(1+ϵee))](1∓^A(1+ϵee)), g=sin[^A(1+ϵee)Δ]^A(1+ϵee) (30)

Similar expression for the inverted hierarchy (IH) can be calculated by replacing (i.e. , (i.e. and ), ). Also the expressions for antineutrino probability () can be obtained by replacing ( ), .

DUNE is proposed as a next generation superbeam experiment at Fermilab, USA Adams:2013qkq ; Acciarri:2016crz ; Strait:2016mof ; Acciarri:2016ooe ; Acciarri:2015uup . It will use existing NuMI (Neutrinos at the Main Injector) beamline design at Fermilab as a source of neutrinos. Its far detector will be placed at Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a distance of 1300 km from the source. DUNE collaboration has planned to use LArTPC (liquid argon time-projection chamber) detector with the volume of 10 kton and 40 kton corresponding to first phase and final phase respectively. In our simulation, we consider 10 kton detector volume. We also consider the flux corresponding to 1.2 MW beam power with 120 GeV protons which gives protons on target (POT) per year. In the analysis, we use method of pulls GonzalezGarcia:2004wg ; Fogli:2002pt as outlined in Gandhi:2007td to account for systematic errors. For the numerical simulation, we use the GLoBES package Huber:2004ka ; Huber:2007ji along with the required auxiliary files messier_xsec ; Paschos:2001np . The remaining details of the DUNE and numerical specifications considered here are same as Nath:2015kjg . We have also added prior on sin.

In Fig. 4, we discuss appearance channel probability () considering diagonal NSI parameter in top row whereas in bottom row, we show our results considering off-diagonal NSI parameter considering 2 GeV neutrino beam energy. We present our results considering neutrinos (antineutrinos) in left (right) panel. Depending on true hierarchy and true octant, we have four combinations of hierarchy - octant namely, NH-HO, NH-LO, IH-HO and IH-LO. We discuss our results for all the four possibilities. The tiny bands in top row are over the range of model-dependent NSI parameter as given by first column of Table 1 444Note that for the illustrative purpose, we describe our results considering one set of solution.. In bottom row, we describe the role of off-diagonal NSI parameter where different bands are over the range of new CP-phase (see figure legends for details). Note that in bottom row, to have better understanding of new CP phase, we fix the value of NSI parameters and the representative value of is taken from Table 2. A general probability level discussion, considering similar set-up, in absence of NSI on various degenerate solutions are noticed in Ref. Nath:2015kjg whereas some noteworthy issues considering model-independent NSI parameters are presented in Ref. Deepthi:2016erc considering DUNE.

We concentrate here on the impact of model-dependent NSI parameters. In literature, (see Liao:2016hsa ; Coloma:2016gei ; Deepthi:2016erc ; Deepthi:2017gxg ) it was pointed out that appearance channel probability suffers with degeneracy due to the presence of intrinsic ), and degeneracy in presence of model-independent NSI parameters. Here we notice considering our model-dependent constrained parameter space of that DUNE has no hierarchy degeneracy as shown in the first plot of top row in Fig. 4. This can be understood by comparing NH bands (blue, brown) with IH bands (yellow, red), as the former have no intersection with latter. From the first plot of bottom row, we also observe that has no hierarchy degeneracy even in the presence of off-diagonal NSI parameter, . Note that with the inclusion of , width of different bands become boarder, this is due to unconstrained range of new CP-phase which effects the measurement of Dirac CP-phase. Similar results are also observed for antineutrinos as shown by right panel and conclusion made for neutrinos remain same for antineutrinos. Moreover, to understand the impact of NSI on the determination of hierarchy and CP violation sensitivity, we extend our study considering our model-dependent NSI parameters at the level. In next section, we illustrate our details.

### 5.3 Sensitivity study

We now discuss the mass hierarchy sensitivity as well as the CP-violation sensitivity of DUNE with model-dependent NSI parameters. For the comparison, we also describe our results considering the standard interactions (SI). Throughout the study, we consider [5+5] years run time for DUNE 555Note that [5+5] means we divide total exposures with 5 years of neutrino run and another 5 years of antineutrino run.. In Fig. 5, we describe our results for the mass hierarchy sensitivity considering true as for lower octant (LO) and for higher octant (HO) in the first and second panel respectively. Whereas in top (bottom) row, we present out results considering true hierarchy as NH (IH) and marginalized over test hierarchy. We take true values of other neutrino parameters as, , , eV and eV which are compatible with the current global-fit data deSalas:2017kay ; Capozzi:2016rtj ; Gonzalez-Garcia:2015qrr . In test, we marginalize over their ranges 666Note that we do not marginalize over .. For the NSI parameters, we marginalize over test considering the range as given in BP-I of Table 1 whereas we keep fixed value for the off-diagonal NSI parameter as given in Table 2 (see figure legends for details). This benchmark point is chosen as value which satisfies all other constraints and is well within allowed parameter region. Note that the four different curves correspond to SI (gray solid), diagonal NSI parameter (black long-dotted), off-diagonal NSI parameter with (blue dash-dotted) and (yellow dotted) respectively. We quantify our mass hierarchy sensitivity as below,

 χ2NH−IH=min∑i[Ni(NHtr,ϵtr,ϕtr)−Ni(IHte,ϵte,ϕtr)]2σ[Ni(NHtr,ϵtr,ϕtr)]2, (31)

where represents the number of events for the oscillation parameters. Also, are marginalized for and kept fixed in both true and test for whereas are considered fixed in both true and test to calculate N. We discuss our results considering a benchmark of confidence level (C.L.) as shown by the horizontal line.

From the gray curve, we observe that DUNE can reach hierarchy sensitivity with the 10 kton detector mass, independent of the true values of and irrespective of the nature of true hierarchy NH or IH, in case of HO (right panel) for SI. Whereas for LO (left panel), we notice it achieves hierarchy sensitivity for all the true values of except the regions around true for neutrinos. In the presence of (see black long-dotted curve), we find that DUNE achieves greater than or almost equivalent to sensitivity for all the mentioned cases. Considering , we describe our results for both CP conserving (i.e., when ) as well as CP violating (i.e., when ) values of new CP-phase. In the case of HO, we notice that even in the presence of off-diagonal NSI parameters DUNE can achieve hierarchy sensitivity for all the true values of . Whereas in case of LO, we find that DUNE attains hierarchy sensitivity irrespective of the true values of for all the cases except the regions around true . Similar results hold for the case of true IH as shown by bottom row, except that for IH over all minima comes near true .

In Fig. 6, we discuss CP-violation (CPV) discovery for DUNE considering both SI and NSI. The left (right) panel describes our result for LO (HO) whereas top (bottom) row discusses our results for true hierarchy as NH (IH). The true values that we considered corresponding to , are based on our model as given in Tables (1, and 2); whereas for the new CP-phase due to NSI we take two cases, namely CP conserving () and CP violating () values. We describe CPV discovery as,

 χ2CPV=min∑i[Ni(δtrCP,ϵtr,ϕtr)−Ni(δteCP(0∘,±180∘),ϵte,ϕtr)]2σ[Ni(δtrCP,ϵtr,ϕtr)]2, (32)

where represents the number of events for the oscillation parameters. Also are defined as Eq. (31).

We draw a line as a benchmark for the discussion of our results. Comparing all curves, we notice that DUNE achieves maximum CPV discovery sensitivity for the SI compare to NSI parameters. Whereas, considering diagonal NSI parameter, we find that it gives almost equal CP sensitivity as SI at 3 C.L for both the hierarchies. Further, when we add off-diagonal NSI parameters with the diagonal NSI parameter, as shown by blue dotted-dash and yellow dotted curves, we observe that CPV sensitivity reasonably decreases for NH. Whereas for IH, we find that CPV sensitivity gets enhanced for CP violating value compare to CP conserving value. Thus, an extra CP phase confuses the measurement of Dirac-CP phase () and hence effects the measurement of overall CP sensitivity even in the case of constraint parameter space of NSIs.

## 6 Summary and Conclusion

Neutrino oscillation experiments opened a new vista to probe the fundamental properties of neutrinos. New physics models, incorporating neutrino masses, are testable in these experiments via the oscillation data. Many of these BSM scenarios give rise to NSI that can be tested in the oscillation experiments. But such models are constrained from lepton flavor violation issues.

The traditional 2HDM is one of these popular models which tries to explain neutrino masses by extending SM with two scalar doublets and right-handed neutrinos. However, this model produces negligible amount of NSI due to almost non-existent interaction of the SM charged leptons and quarks with the neutrinos, determined by the tiny mixing () of the two scalar doublets. In this study, we propose a modified 2HDM which is an improvement over the usual 2HDM. We can have sizable NSI parameter while maintaining LFV constraints. We achieve this by assigning a charge to under a global symmetry. This can lead to an observationally significant NSI along with the presence of a tiny Dirac neutrino mass, keeping the original motivation of 2HDM intact. This modification of 2HDM reintroduces hierarchy in the Yukawas in the sector, but simultaneously eases the hierarchy in the side, compared to those in the the traditional 2HDM. Softly broken global 2HDM allows the presence of heavy neutral BSM scalar, that helps to address the stringent bound on electrophilic Yukawa of an ultra-light neutral scalar as well as the tight constraints from the oblique parameter () measurements. Combined effects of the LEP constraints like , mono-photon search along with the bound on mass cuts short a significant amount of allowed parameter space of the modified 2HDM case, therefore putting stringent bounds on the NSI parameters. Presence of lepton flavor violating (LFV) decays like , put stringent upper bound on the Yukawa couplings and that results in any NSI parameter involving negligible apart from significant reduction of upper bound of other NSI parameters.

Depending on these constraints, this model predicts the range of permissible NSI parameter . We also find that the only possible off-diagonal NSI parameter in this model is whereas remaining NSIs become insignificant under model constraint. Thus, the effects of NSI parameters involving such as etc. is not studied here. Later, we study the impact of these NSIs considering LBL experiment like, DUNE. At the probability level, considering model-dependent NSIs, we observe no wrong hierarchy degeneracy even in the presence of off-diagonal NSI parameter. Furthermore, at the level, we find that DUNE shows around 5 hierarchy sensitivity when one adopt NH as well as IH both as a true hierarchy considering one at a time. These results remain valid irrespective of the value of true Dirac CP phase, . From our study of CP discovery, we notice that CP violation in leptonic sector gets affected even in the presence of model-dependent diagonal NSI parameter. Further, we observe that extra CP-phase, due to off-diagonal NSI parameter, confuses the measurement of Dirac-CP phase and hence affects the measurement of overall CP sensitivity.

###### Acknowledgements.
UKD acknowledges the support from Department of Science and Technology, Government of India under the fellowship reference number PDF/2016/001087 (SERB National Post-Doctoral Fellowship). UKD thanks Dr. Tirtha Sankar Ray for useful discussion. The research work of NN was supported in part by the National Natural Science Foundation of China under grant No.11775231. The authors would like to thank Prof. Srubabati Goswami for her insightful comments and careful reading of the manuscript.