750 GeV Diphoton excess from Gauged $B-L$ Symmetry

We show that the recently observed 750 GeV diphoton excess at LHC can be due to the decay of a $SU(2)_L$ singlet scalar particle having 3 units of charge under gauged $B - L$ symmetry. Such a particle arises as an essential ingredient of recently studied gauged $B -L$ extension of the Standard Model with unconventional charge assignment for right handed neutrinos. Apart from being one of the simplest extensions of the Standard Model, the model also contains a dark matter candidate and Dirac neutrinos with naturally small masses.


I. INTRODUCTION
The ATLAS and CMS collaborations at the LHC have recently reported an excess of events in the invariant mass distribution containing two photons at √ s = 13 TeV [1][2][3]. The ATLAS Collaboration [2], with 3.2 fb −1 data, has reported an excess of 3.9σ at diphoton invariant mass around 750 GeV. The significance reduces to 2.3σ once the Look Elsewhere Effect is included. This corresponds to an excess in signal σ(pp → γγ) of about (10 ± 3 fb) with a best fit width ∼ 45 GeV. The value of the experimental acceptance in ATLAS is about 0.4.
The CMS collaboration has also found an excess in diphoton events with local significance 2.6σ [3] at √ s = 13 TeV with 2.6 fb −1 data at mass around 750 GeV. This significance reduces to 2.0σ if large width (∼ 45 GeV) is assumed. This translates to an excess in signal cross section σ(pp → γγ) of about (6 ± 3 fb).
These excess events do not have any significant missing energy, leptons or jets associated with them. No excess of events have been found in ZZ, dilepton, dijet channels in the same invariant mass region for new data. Although this excess could well be a statistical fluctuation, it has drawn significant attention as it can also arise from decay of new particle with mass around 750 GeV [4].
If the observed diphoton excess indeed corresponds to decay of a hitherto unknown particle then this will be the first confirmation of new physics beyond the Standard Model (SM). If the observed excess is due to a resonance it has to be a boson and it cannot be a spin-1 particle [5,6]. This leaves the possibility of it being either a spin-0 or spin-2 particle. If it is indeed a new particle then one must wonder what kind of new physics incorporates it. This issue can only be settled by looking at various possible new physics scenarios that can potentially lead to such a particle. It is generally expected that this new physics will also be related to other open problems in high energy physics which do not have a satisfactory explanation within SM. Chief among them is the problem of neutrino masses (and their relative smallness) and the nature of dark matter. It will be quite satisfying if the observed new particle has a natural connection with the models addressing at least one of these issues. In this work we assume this resonance to be a spin-0 particle and look at one such promising model based on gauge B − L symmetry. This model was recently proposed to explain the smallness of neutrino mass if neutrinos are Dirac particles and also has an interesting dark matter candidate [7][8][9].
The gauged B − L symmetry is one of the simplest and most well studied extension of SM [10,11]. In SM, Baryon number B and Lepton number L are accidentally conserved classical symmetries. However, both B and L currents are anomalous and only the combination B − L is anomaly free. In the conventional gauged B − L model, the B − L symmetry is promoted to an anomaly free gauge symmetry by addition of three right handed neutrinos ν i R each transforming as −1 under the U (1) B−L [10,11]. It was shown that if this B − L symmetry is spontaneously broken by a SU (2) L singlet scalar χ 2 having two units of B − L charge, then the right handed neutrinos can acquire a Majorana mass term M R proportional to the vacuum expectation value (vev) u 2 of the singlet scalar. Moreover if the B − L breaking scale is far greater than the electroweak scale then the right handed neutrinos acquire a large mass leading to a natural implementation of Type I seesaw mechanism. However, in this scenario the B −L breaking scale is expected to be very high (same as seesaw scale) and it is very difficult to test this model at LHC.
Recently, another simple choice of B − L charges for right handed neutrinos which leads to anomaly free U (1) B−L gauge symmetry has been proposed. Unlike the previous case, here the three right handed neutrinos transform as ν i R = (+5, −4, −4) under B − L symmetry [7,12]. It was shown that such a charge assignment can lead to Dirac neutrinos with naturally small masses if the B − L symmetry is spontaneously broken by a SU (2) L singlet scalar χ 3 transforming as ∼ 3 under U (1) B−L symmetry. This new B − L model can also have a candidate for long lived scalar dark matter if two new singlet scalars transforming as χ 2 ∼ 2 and χ 6 ∼ −6 under U (1) B−L symmetry are added to it [9]. Unlike the conventional B − L model, here the B − L breaking scale need not be high and can be well within the reach of LHC. This opens up the possibility of testing various features of this model in the present run of LHC. Moreover, the dark matter in this model has a significant interaction with the nuclei and can be detected in present or near future dark-matter direct-search experiments. We refer the interested readers to [7][8][9] for further details. Thus, apart from providing a explanation for the nature and small masses for neutrinos as well as a candidate for dark matter, various aspects of this new B − L model are quite testable both in colliders and dark matter direct detection experiments. Furthermore, as we will discuss in the subsequent sections, owing to presence of singlet scalars, this model is ideally suited to explain the recently observed 750 GeV diphoton excess and the aim of this paper is to look at this possibility in details.
The plan of the paper is as follows. In Section II we look at the details of the gauged B − L symmetry model which is a slightly extended version of the previously discussed model. We show that the modified model is also free of anomalies. In Section III we look at the details of the scalar and Yukawa sectors and identify a viable candidate which can explain the observed 750 GeV diphoton excess. In Section IV we discuss the production and decay of the 750 GeV particle and compare our computation with the experimental results. We finally conclude in V.

II. THE GAUGE B − L SYMMETRY MODEL
The anomaly free gauged B − L model with unconventional B − L charges for right handed neutrinos was originally constructed to obtain Majorana neutrinos [12] or Dirac neutrinos [7,8] with naturally small masses. The model for Dirac neutrinos was further extended in [9] to accommodate a long lived dark matter particle. The SU (2) L singlet scalars of the model are required to break the gauge B − L symmetry and we show in this work, that a linear combination of these scalar can be a viable candidate for 750 GeV resonance, whose decay can provide a possible explanation for the recently observed diphoton excess [2,3]. In this work we study a slightly extended version of the model discussed in [9], where we have also included two SU (2) L singlet vector "quarks" X L,R , Y L,R . Although they are SU (2) L singlets, these exotic quarks do carry SU (3) c colour charge as well as U (1) Y , U (1) B−L charges. The SU (3) c × SU (2) L × U (1) Y and U (1) B−L charge assignment for the fermions in the model are as follows: In Table I apart from the SM particles we have also included three right handed neutrinos ν i R , three SU (2) L singlet heavy fermions N i L,R (as in the previous model [9]) and two pair of exotic "quarks" X L,R , Y L,R which carry color and electromagnetic charges but are singlet under SU (2) L .
The charge assignment for the scalars in this model (which are same as in [9]) are as follows:  In Table II, Φ = (φ + , φ 0 ) T is the usual SU (2) L doublet scalar and χ i are SU (2) L singlet scalars. The new fermions introduced in the model can potentially lead to anomalies. Thus, it is important to ensure that the model is anomaly free. The new particles can induce following triangular anomalies: It has already been shown in [7] that for the case when exotic quarks X, Y are not present, the model is completely anomaly free. It can also be easily seen that the addition of the X, Y quarks does not spoil the anomaly cancellation and hence the model remains anomaly free.

III. THE SCALAR AND YUKAWA SECTORS
In this section we look at the details of the scalar and Yukawa sector of our model and identify the candidate for 750 GeV resonance. The scalar potential of our model is given by The Minimum of V is given by: where φ 0 = v, χ 3 = u 3 , χ 6 = u 6 , are the vev of the scalar fields. Moreover, just like in [9], here also the singlet scalar χ 2 does not acquire any vev i.e. χ 2 = 0. The minimum of V is determined by Since χ 2 = 0, there is one dark-matter scalar boson χ 2 with mass given by There is also one physical pseudoscalar boson with mass given by There are three physical scalar bosons spanning the basis [h, √ 2Re(χ 3 ), √ 2Re(χ 6 )], with 3 × 3 mass-squared matrix given by A. Simplifying Scenario The mass matrix in Eq. 15 can be diagonalized to give three CP even scalars which will be linear combinations of Φ, χ 3 , χ 6 scalars. However for sake of illustration, we look at a special case of the generic mass matrix in Eq. 15 which takes a simple form if we assume where a and b are two independent parameters. With these simplifying assumptions, the mass matrix of Eq. 15 becomes The eigenvalues of the mass matrix in Eq. 17 are given by The masses of the scalars are then given by with one scalar having the mass same as the 125 GeV resonance and another having mass 750 GeV. For sake of definiteness we will identify the first eigenstate with the 125 GeV scalar (henceforth called "Higgs") and the second eigenstate as 750 GeV scalar i.e. we demand m 1 = 125 GeV and m 2 = 750 GeV. The mass of the third scalar then depends on the value of a and b. Solving for a and b we find that a = 108.5 GeV and b = 750 GeV leads to desired masses for m 1 and m 2 scalars. The mass of the third scalar m 3 then becomes m 3 = 1.30 TeV.
The masses of the pseudoscalar, dark matter and Z are dependent on the value of other free parameters e.g. the value of vevs u 3 , u 6 , the U (1) B−L coupling g X as well as on the quartic coupling of scalars λ ij . We also like to note that in the simplified mass matrix of Eq. 17, not all of λ ij are independent parameters owing to Eq. 16. The mass of dark matter χ 2 is also dependent on the additional parameter m 2 2 and the quartic couplings λ i2 ; i = 0, 3, 6. Thus for a large range of parameter space, we can also have heavy Z as required by constraints from LUX dark matter direct detection experiment [9,13]. Since the mass of the dark matter m χ2 depends on additional free parameters therefore it can either be greater than or less than 750/2 = 375 GeV. This leads to two distinct cases; if m χ2 ≤ 375 GeV then the 750 GeV resonance can decay into dark matter and it can lead to significant invisible decay width. If m χ2 > 375 GeV then this decay is kinematically forbidden. In later sections we will study both these cases in details.

B. Higgs Coupling to SM Gauge Bosons
In Section III A, we showed that in a simplified mass matrix, one of the scalar combinations can be identified with the 125 GeV Higgs recently discovered at LHC. The recent data from both ATLAS and CMS experiments suggest that this scalar has couplings very similar to SM Higgs couplings with the SM gauge bosons. In this section we show that in our model also, there exists a decoupling limit where one of the scalars will have almost SM like couplings with the SM gauge bosons.
For sake of simplicity we will work with the simplified mass matrix of Eq. 17. The mass matrix mixes φ R 0 , χ R 3 , χ R 6 states with each other. This matrix can be diagonalized by an orthogonal matrix with the diagonal matrix corresponding to the masses of the three physical scalars as shown in Eq. 18. The physical scalars are then given by where tan 2θ = 2 √ 2ab 3b 2 − 2a 2 . It is clear from Eq. 20 that as sin θ → 0, h couplings to SM gauge bosons become SM like. For the case of a = 108.5 GeV and b = 750 GeV we find that cos θ = 0.997 sin θ = 0.069. This implies that in our model the couplings of the scalar h of mass 125 GeV with W, Z gauge bosons are almost SM like. The small deviations from the SM like couplings are well within the experimental limits [14]. It should be noted that in this limit the other scalars as well as the Z boson can be made heavy (assuming all couplings to be O(1) ) in congruence with the experimental bounds for these particles [9].

C. The Yukawa Sector
Apart from its coupling to the scalars, χ 3 has following Yukawa couplings As evident from Eq. 21 both quarks X, Y acquire mass after spontaneous breaking of B − L symmetry and the masses are proportional to the vev u 3 of χ 3 . The Yukawa coupling of χ 3 translates into Yukawa coupling of the scalar H 1 ≡ χ 3 − χ 6 . The Yukawa couplings of other scalars with fermions are same as in [9] and we refer the interested reader to [9] for further details. Owing to the coupling of χ 3 with quarks X, Y ; the 750 GeV scalar H 1 can be efficiently produced through gluon-gluon fusion at LHC. The production and decay of this scalar are discussed in details in the next section.

IV. PRODUCTION AND DECAY OF H1
In this section we look at the details of the production and decay channels for the 750 GeV scalar at the LHC. In particular, we show that the decay of this scalar to two photons can lead to the observed diphoton excess. Moreover, as we will show, the decay of H 1 to SM fermions and the Higgs are suppressed thus explaining the non-observation of any excess in other channels.
Since, χ 3 couples to quarks X, Y through its Yukawa couplings Eq. 21, therefore, at LHC H 1 = (χ R 6 − χ R 3 ) can be efficiently produced by gluon-gluon fusion through triangular loop involving X, Y . If both X, Y are heavier than m H 1 2 then the tree level decay of H 1 to both X, Y is kinematically forbidden. In such a case its decay to two photons through triangular loop involving X, Y as shown in fig 1, can be significant leading to the observed anomaly in diphoton channel. then it can also decay into a pair of dark matter particles as shown in 4. This can lead to appreciable invisible decay width for H 1 .
As the SM like Higgs (h) couples to both χ 3 and χ 6 with same coupling strength, this will result in cancellation of its interaction strength with H 1 ≡ χ 3 − χ 6 , resulting in very small H 1 hh coupling. Hence the H 1 → hh decay will be vanishingly small. Furthermore, since H 1 is primarily a mixed state of SU (2) L singlets χ 3 and χ 6 as shown in Eq. 20 therefore it has negligible tree level coupling to SM fermions as well as W and Z gauge bosons. Therefore, its decay in other channels like dilepton, dijet and diboson is extremely suppressed. This observation is also in congruence with the experimental results which show lack of any statistically significant excess in these channels.
Since, neither ATLAS nor CMS has seen any hint or anomalous excess in any channel for masses below 750 GeV, we further require that all the other new particles (except dark matter) should be sufficiently massive. Thus we require that the mass of the pseudoscalar m A > 1 TeV. The mass of other CP even scalar H 2 is also greater than 1 TeV as obtained in Section III B. The mass of the dark matter m χ2 is left as a free parameter and depending on its mass, the H 1 → χ 2 χ * 2 decay may or may not be kinematically forbidden 1 . We will consider both cases in next section. Also, we limit the mass of the Z boson m Z > 12 TeV which is well above the dark matter direct detection constraints from the LUX experiment [9,13].
Thus the only prominent decay modes of interest are H 1 → γγ, H 1 → gg and if m χ2 ≤ m H1 /2 then H 1 → χ 2 χ 2 also. The partial decay widths of H 1 in these modes are given as, 1 Here χ * 2 denotes the anti-particle and should not be confused with an off-shell particle.
with m i , Q i being corresponding fermion (X, Y ) masses and electromagnetic charges respectively.
The f i s here denote fermion Yukawa couplings with the scalar H 1 whereas α s , α denote strong and electromagnetic interaction gauge coupling strengths. N c is the color factor which is 3 for the quarks. Also without loss of generality we have normalized the dimension full coupling between χ 2 and H 1 by the vev u 3 with κ χ2 being a dimensionless parameter. The parameter κ χ2 is a function of the vevs u 3 , u 6 as well as the quartic couplings between χ 2 with χ 3 , χ 6 fields as given in Eq. 7. The f (τ i ) for our case, where m X,Y > m H1 /2 is given as In addition to these decay modes, H 1 can also decay to Zγ and ZZ through triangular loops involving the X, Y quarks. For our case of m 2 Z << m 2 H1 the decay width to Zγ and ZZ are given by where s W = sin θ W , c W = cos θ W and θ W is the electroweak angle. Since the vector quarks X, Y are both SU (2) L singlets and all the three decays namely H 1 → γγ, Zγ, ZZ proceed through the same triangle loop the ratio of partial decay width in these three channel are given by As clear from Eq. 28, the loop decays of H 1 to Zγ and ZZ are suppressed compared to the γγ decays by a factor proportional to the electroweak angle. Thus, H 1 is an ideal candidate to explain the observed diphoton excess and lack of significant excess in other decay channels.

A. Numerical Results
In this section we present the numerical results and the allowed parameter range for the masses of the quarks X, Y and the Yukawa coupling which can explain the excess observed at the LHC in the diphoton channel. For this part we have used MadGraph5aMC@NLO [15] with NN23LO1 PDF set [16] to obtain the numerical estimates taking K factor of 1.5 into account for NLO correction [17].
For H 1 to be a viable candidate to explain the observed diphoton excess it not only has to explain the LHC data for 13 TeV run but should also satisfy the nonobservance of any statistically significant excess in the previous 8 TeV run in various channels. In our model the only significant decay channels for H 1 are the loop induced gg, γγ, Zγ and ZZ decays. Moreover, if m χ2 ≤ m H1 /2 then it can also decay to two dark matter particles through H 1 χ 2 χ * 2 tree level couplings. The 8 TeV constraints on production × branching fraction σ × Br(h 1 → f i f j ); f i,j ≡ g, γ, Z, χ 2 on these channels are [19,20]: As mentioned before, since in our model the dark matter mass m χ2 is not fixed so there arise two distinct possibilities; either m χ2 > m H1 /2 or m χ2 ≤ m H1 /2. For first case, H 1 decay to two dark matter particles is kinematically forbidden and the only prominent channels are its loop decays to gg, γγ as well as tan θ W suppressed loop decays to Zγ and ZZ. In the second case, H 1 can also decay to two dark matter particles. Here we analyze both these possibilities in Subsection IV A 1 and Subsection IV A 2 respectively.
1. Case-I : mχ 2 > mH 1 2 In this case the only important decay modes for H 1 are H 1 → γγ and H 1 → gg along with H 1 → Zγ and H 1 → ZZ both of which are θ W suppressed. All of these decay modes are loop level, going through triangle loops involving X, Y quarks. Since, g, γ and Z all couple to the quarks through gauge interactions so their interaction strengths are fixed, and are proportional to α s , α, the strong and electromagnetic coupling constants respectively. Hence, for this case of our model the production and decay rate of H 1 depends on only two free parameters, the masses of X, Y quarks and the Yukawa coupling between H 1 and quarks. Moreover, since the quarks X, Y acquire mass through the vev of χ 3 so the Yukawa coupling can be equivalently replaced by the vev u 3 as a free parameter. The allowed mX -u3 range corresponding to CMS (green), ATLAS (deep blue) and the overlap (light blue) ranges with 95% confidence level. Also, shown is the 95% confidence level H1 → γγ exclusion line (red dashed) from 8 TeV run with the regions on the left of the line being incompatible with it. The black shaded region is also excluded by the perturbativity constraints.
In Fig. 5 we show the allowed ranges of the exotic quark masses and the value of the vev, χ 3 = u 3 , that can explain the observed 750 GeV diphoton excess for both the CMS and ATLAS experiments within 95% confidence level. In obtaining the numerical results, for simplicity we assume that the masses of the exotic quarks X, Y are degenerate i.e. m X = m Y and treat them as a single parameter m X . In addition we require that all the couplings in our model remain perturbative. The region of the parameter space excluded due to non-perturbativity of the couplings is explicitly shown in Fig. 5.
Furthermore, in plotting Fig. 5 we have imposed the 8 TeV exclusion limits for the heavy scalar of mass 750 GeV in all other channels. The strongest constraint from 8 TeV exclusion limits actually comes from non-observance of any statistically significant excess in the γγ decay channel. In Fig. 5 the thick red line corresponds to γγ exclusion limit of Eq. 29. The parameter space on the left of the red line is incompatible with the 8 TeV data.
As mentioned before, the scalar H 1 that we are considering here, does not couple to SM fermions at tree level. Therefore the limits given in [18] can be easily satisfied. The coupling H 1 hh is also negligibly small and σ(pp → H 1 → hh) is well under the experimental limit [21]. Moreover, the scalar is a neutral SU(2) singlet. Therefore, it does not have any tree level coupling to either W or Z bosons. As the exotic fermions are SU(2) singlet even the H 1 → W W decay through the triangle loop is not possible. However, it can couple to ZZ, Zγ at loop level through triangle loop of the exotic fermions and it has to be taken into account.
Compared to the H 1 → γγ the ZZ, Zγ decays are suppressed. Furthermore, as shown in Eq. 29 the exclusion limits on these decay channels are also relatively weaker. Thus the constraints from these decay channels are quite weak and do not impose any additional constraints on the allowed parameter range shown in the plot. Finally the H 1 decay to gluons is also well below the experimental limit and does not impose any addition constraints on the allowed parameter range. As an example the values for these decay channels for a benchmark point (m X = 1 TeV and u 3 = 205 GeV) on the γγ exclusion line of Fig. 5 are given as As clear from Eq. 30, for this case of our model, apart from the γγ decay channel, the constraints from all other decay channels are easily satisfied. Even for the γγ channel, our model has enough parameter space compatible with both the observed 13 TeV excess and the 8 TeV constraints.
In this case, in addition to the decay channels discussed in previous section, H 1 decay to dark matter is also kinematically allowed and it can have appreciable invisible decay width. In Fig. 6 we show the allowed parameter range for the exotic quark masses and u 3 , that can explain the observed 750 GeV diphoton excess for both the CMS and ATLAS experiments within 95% confidence level. In plotting Fig. 6 we have taken κ χ2 = 0.5 and have also imposed all the constraints from 8 TeV run listed in Eq. 29.
As in the previous case, again the only significant constraint from 8 TeV run comes from the γγ decay channel which is also plotted in Fig. 6. The parameter space on the left of the red line is incompatible with the 8 TeV data. As before, the constraints from other channels including the invisible decay to dark matter are rather weak and do not give any additional constraint. As an example the values for these decay channels for a benchmark point (m X = 1 TeV and u 3 = 180 GeV) on the γγ exclusion line of Fig. 6 are given as As clear from Eq. 31, like the previous case here also only the constraints from γγ channel for 8 TeV run are important. The constraints from all other channels are comfortably satisfied. Furthermore, just like the previous case, in this case also our model has enough parameter space compatible with both the observed 13 TeV excess and the 8 TeV constraints. Thus the 750 GeV excess can be understood in our model as the decay of H 1 to a pair of photons.
Finally, before ending the section we like to discuss briefly about the total decay width of H 1 for the two cases. The first thing to note is that given the current low statistics, the estimates of decay width are very poor. This aspect is highlighted by the fact that while CMS data prefers narrow decay width of around a few GeV for the resonance, the ATLAS prefers a relatively broader resonance with decay width ∼ 45 GeV. Thus the current estimates of decay width are highly uncertain and are likely to change significantly in the future runs.
In our model, if the H 1 decay to dark matter is kinematically forbidden then the dominant decay channels will all be loop induced with H 1 → gg being the most significant. In such a scenario H 1 will be a narrow resonance with total decay width up to a few GeVs. However, if H 1 decay to dark matter is kinematically allowed then it can have significant invisible decay width owing to the fact that such a decay is not loop suppressed. In this case the H 1 can be a broad resonance.
If in future runs the ATLAS experiments estimates of a broad resonance persists then for our model it will imply a significant invisible decay width. Depending on the value of κ χ2 , H 1 can have decay width upto 40 GeV, albeit for a small parameter range. In such a case, a better solution can be obtained by adding a pair of SU (2) L singlet charged leptons to our model. However, at this stage we feel that such an extension of our model is premature and not necessary.

V. CONCLUSIONS
The recently observed 750 GeV diphoton excess has drawn significant attention as it could be the first signs of new physics at LHC. Although it is too early to conclude that this is a definite sign of new physics but nonetheless it raises an intriguing possibility that it might originate from decay of a hitherto unknown particle. In this work we have looked at the extended gauged B − L symmetry model with unconventional charges for the right handed neutrinos as a possible candidate new physics model to explain the 750 GeV excess. The model was originally constructed to obtain Dirac neutrinos with naturally small masses and also has a long lived dark matter particle. Unlike the conventional gauged B − L symmetry model where the B − L scale is expected to be quite high, being related with the seesaw scale, in our model the B − L scale can be well within the LHC range thus opening up the possibility of testing its various aspects at LHC.
We have looked at the possibility that the observed diphoton excess can arise due to decay of the scalar particle H 1 into two photons. This scalar in our model is predominantly composed of the singlet scalars χ 3 and χ 6 which are essential ingredients of the model. They are required in order to spontaneously break the gauged B − L symmetry as well as to obtain Dirac neutrinos with small masses. We have further shown that the model not only explains the diphoton excess but also satisfies all the other experimental constraints like non-observation of any excess in dilepton, dijet, diboson and invisible channels. It also has a 125 GeV particle h which has almost SM Higgs like couplings to the other SM particles and satisfies all the other experimental constraints for the 125 GeV scalar. Moreover, since in our model h is predominantly composed of the SU (2) L doublet scalar with very small admixture from SU (2) L singlet scalars, it naturally explains why the 125 GeV particle has almost SM like couplings.
Thus to conclude, the gauged B − L model considered here appears to be a promising candidate for new physics. It has all the right ingredients to explain not only the 750 GeV diphoton excess but all the other experimental results both for the 750 GeV resonance as well as the 125 GeV resonance. Moreover, the model also connects the observed new physics with the already well know and long standing problems of neutrino masses and dark matter and attempts to provide a unifying solution to all of them. Also, it has several testable predictions like existence of heavier particles in ∼ 1 TeV range, Dirac nature of neutrinos and candidate for dark matter. These aspects can be tested in future run of LHC as well as in dark matter direct detection experiments and in various neutrino physics experiments.