Multi-component Fermionic Dark Matter and IceCube PeV scale Neutrinos in Left-Right Model with Gauge Unification

We consider a simple extension of the minimal left-right symmetric model (LRSM) in order to explain the PeV neutrino events seen at the IceCube experiment from a heavy decaying dark matter. The dark matter sector is composed of two fermions: one at PeV scale and the other at TeV scale such that the heavier one can decay into the lighter one and two neutrinos. The gauge annihilation cross sections of PeV dark matter are not large enough to generate its relic abundance within the observed limit. We include a pair of real scalar triplets $\Omega_{L,R}$ which can bring the thermally overproduced PeV dark matter abundance into the observed range through late time decay and consequent entropy release thereby providing a consistent way to obtain the correct relic abundance without violating the unitarity bound on dark matter mass. Another scalar field, a bitriplet under left-right gauge group is added to assist the heavier dark matter decay. The presence of an approximate global $U(1)_X$ symmetry can naturally explain the origin of tiny couplings required for long-lived nature of these decaying particles. We also show, how such an extended LRSM can be incorporated within a non-supersymmetric $SO(10)$ model where the gauge coupling unification at a very high scale naturally accommodate a PeV scale intermediate symmetry, required to explain the PeV events at IceCube.


Introduction
There have been several irrefutable experimental observations suggesting the presence of dark matter (DM) in the Universe, starting from the galaxy cluster observations by Fritz Zwicky [1] back in 1933, observations of galaxy rotation curves in 1970's [2], the more recent observation of the bullet cluster [3] to the latest cosmology data provided by the Planck satellite [4]. The latest Planck data suggest that around 26% of the present Universe's energy density consist of non-baryonic or dark matter. In terms of density parameter and h = (Hubble Parameter)/100, the present dark matter abundance is conventionally reported as [4] Ω DM h 2 = 0.1187 ± 0.0017. (1.1) Even though there are enough evidences from astrophysics and cosmology confirming the presence of dark matter in the Universe, no experiment has so far been able to probe it directly. The most recent dark matter direct detection experiments like LUX, PandaX-II have also reported their null results [5][6][7]. There have been many efforts to look for indirect dark matter signatures at different experiments with the hope that even though dark matter may not scatter off nuclei significantly in order to be consistent with null results at direct detection experiments, but they may decay or annihilate into the standard model (SM) particles on cosmological scales and leave some indirect signatures. The IceCube experiment at the south pole is one such place where such indirect dark matter signatures can be looked for.
The IceCube collaboration has in fact reported 54 ultra-high energy (UHE) neutrino events corresponding to deposited energy in the range 20 TeV to 2 PeV in their 4 year dataset [8][9][10], with a 6.4σ excess over the atmospheric background. Recently, IceCube has observed a track-like event which extends its deposited energy to ∼ 3 PeV [11]. With no significant evidence for astrophysical sources [12,13], this observation has led to several beyond standard model (BSM) proposals as possible interpretation of these high energy neutrino events. Among them, probably the simplest BSM interpretation of these high energy events is in terms of a PeV scale long-lived dark matter candidate that can decay into neutrinos . Among these, we are particularly interested in the framework adopted by [27,37]. The authors in both of these works considered a type of Left-Right Symmetric Model with relevant field content to generate the high energy IceCube events from the decay of a PeV dark matter candidate. The LRSM is one of most highly motivated BSM frameworks which in its generic form [39][40][41][42][43], not only explains the origin of parity violation in weak interactions but also explains the origin of tiny neutrino masses naturally. The gauge symmetry group and the field content of the generic LRSM can also be embedded within grand unified theory (GUT) symmetry groups like SO(10) providing a non-supersymmetric route to gauge coupling unification. The right handed fermions of the SM forms doublet under a new SU (2) R group in LRSM such that the theory remains parity symmetric at high energy. This necessitates the inclusion of the right handed neutrino as a part of the right handed lepton doublet. Both the works [27,37] considered the lightest right handed neutrino as the PeV dark matter candidate. However, in generic LRSM, such a PeV right handed neutrino will be extremely short lived compared to the age of the Universe, due to its decay into the SM fermions mediated by the SU (2) R vector bosons. Therefore the authors of [27,37] considered two different types of LRSM which can explain the IceCube result well but not as motivating as the generic LRSM from grand unification point of view. For example, the work [27] considered the charged isospin partner of the right handed neutrino to be a super-heavy fermion instead of the usual right handed charged lepton, in order to prevent the fast decay of dark matter into the SM particles. On the other hand, the authors of [37] considered a hadrophobic SU (2) R to prevent the W R mediated fast decay of the dark matter candidate.
Instead of pursuing this non-conventional LRSM route taken in [27,37], here we consider a simple extension of the minimal LRSM to incorporate the IceCube observation while retaining other generic features of LRSM including the possibility of embedding it within SO(10) GUT. The possibility of dark matter within GUT has been explored quite extensively within supersymmetric frameworks. Within non-supersymmetric SO(10) GUT also, similar studies have appeared in some recent works including [44,45]. The interesting feature of these scenarios is the natural origin of the symmetry stabilising the dark matter candidate within a UV complete theory that can achieve gauge coupling unification and can also predict the proton decay lifetime. While we do not follow a general top-down approach here similar to [44,45], we start with a left-right symmetric model having necessary particle content to produce the correct dark matter relic abundance along with the IceCube high energy neutrino events and then consider the possible embedding of these particles within SO(10) multiplets. This approach is similar to the bottom-up approach followed in LRSM dark matter studies with gauge coupling unification discussed in the recent works [46][47][48]. In our model, the neutral components of two fermion triplets are dark matter candidates with the heavier of them having PeV scale mass and the lighter one having mass at TeV scale. Although the abundance of TeV scale DM can be kept within the observed DM abundance in the usual thermal freeze-out mechanism, the PeV scale DM becomes over-abundant due to insufficient annihilations or the unitarity bound [49,50]. We introduce two different types of scalar multiplets: a pair of triplet and one bi-triplet which serve two purposes. One of them can be long lived and decay after the freeze-out of PeV DM so that its abundance can be diluted by the release of entropy [51] while the other can assist the PeV DM decay into the lighter one and a pair of high energy neutrinos. The long-lived nature of such scalars as well as the heavier DM requires arbitrary fine-tuning of different parameters of the model. However, the presence of an approximate global U (1) X symmetry can naturally explain the origin of such tiny couplings required for long-lived nature of these decaying particles.
We show that the additional fields included to the minimal LRSM in order to accommodate a heavy long lived dark matter decaying into the high energy neutrinos observed at the IceCube experiment also assist in achieving successful gauge coupling unification at a high scale. Apart from the usual LRSM particles, the additional fields included in this work also can be accommodated within different SO(10) multiplets. The gauge coupling unification at a high energy scale is consistent with an intermediate scale around PeV scale. Although this could well be a coincidence, but it could also be an indication suggesting a more fundamental origin of a dark matter candidate having mass around a few PeV. Interestingly, the intermediate Pati-Salam symmetry in generic SO(10) GUT models has a lower bound from phenomenological considerations of rare meson decays, which lie near the PeV scale [52]. Also, such PeV scale intermediate left-right symmetry has other cosmological motivations [53]. Therefore, the PeV dark matter in our model is not just motivated from phenomenological considerations, but it could also have a deeper theoretical origin as this PeV intermediate scale arises naturally from the demand of successful gauge coupling unification and consistent cosmology. This paper is organised as follows. In section 2, we discuss our model in details. We discuss the gauge coupling unification along with SO(10) embedding in section 3. In section 4, we discuss the details of dark matter relic abundance calculation and constraints on the parameter space from the requirement of satisfying the relic density bound. In section 5, we discuss the way the IceCube high energy neutrino events can be explained from PeV dark matter decay in our model. We finally conclude in section 6.

The Model
The left-right symmetric model [39][40][41][42][43] is one of the most widely studied BSM framework that can simultaneously explain the origin of tiny neutrino masses and parity violation at weak interactions. The gauge symmetry of the standard model namely, such that the right handed fermions transform as doublets under SU (2) R , making the theory left-right symmetric. The model also has an in-built discrete Z 2 symmetry or D-parity which ensures the equality of couplings in SU (2) L,R sectors. The effective parity violating electroweak physics at low energy arises as a result of spontaneous breaking of the SU (2) In the manifest left-right symmetric model SU (3) c × SU (2) L × SU (2) R × U (1) B−L , the scale of SU (2) R gauge symmetry breaking and parity breaking are identical which is not necessary. The original idea of D-parity invariance in left-right symmetric models is proposed in [54][55][56] where the discrete parity symmetry gets broken much before the SU (2) R gauge symmetry breaking. The key difference between Lorentz parity and D-parity is that Lorentz parity acts on the Lorentz group and interchanges left-handed fermions with the right-handed ones but the bosonic fields remain the same whereas D-parity acts on the gauge groups SU (2) L × SU (2) R interchanging the SU (2) L scalar fields with the SU (2) R scalar fields in addition to the interchange of fermions. The effect of the spontaneous breaking of D-parity results in an asymmetry between left and right-handed scalar fields making the coupling constants of SU (2) R and SU (2) L evolve separately under the renormalization group leading to unequal SU (2) L and SU (2) R gauge couplings. To illustrate the idea, consider the scalar sector of the left-right model with spontaneous D-parity breaking mechanism consists of a SU (2) singlet scalar field σ which is odd under discrete D-parity, two scalar triplets ∆ L , ∆ R and a bidoublet Φ. After assigning a non-zero vacuum expectation value (vev), σ to D-parity odd scalar singlet σ, the left-right symmetry with D-parity where µ ∆ is the mass term for triplets i.e, µ 2 ∆ Tr ∆ † L ∆ L + ∆ † R ∆ R and λ ∆ is the trilinear coupling in the term M σTr ∆ † L ∆ L − ∆ † R ∆ R . With a large parity breaking vev to the D-parity odd scalar singlet; σ ∼ M Pl , the model yields all the left-handed (LH) scalars to have heavy masses i.e., O(M Pl ) while those of the right-handed (RH) scalars can have much lighter masses near the TeV scale with M 2 In fact M ∆ R can have any value below M Pl depending on the degree of fine-tuning in the scalar coupling λ ∆ . Thus, the spontaneous D-parity breaking is not only gives mass splitting between left and right handed scalar fields, but the asymmetry in the scalar sector at the energy scales below D-parity breaking scale M Pl creates an asymmetry in the gauge couplings, g L = g R for the surviving left-right gauge group. We implement this nice idea in the present model to provide required mass splitting between Dark matter fermions as well as scalars.
The usual fermion content of the LRSM along with the dark matter candidates Σ L,R added in the present model are shown in table 1. The scalar content of the minimal LRSM consist of a bidoublet Φ responsible for generating Dirac mass terms of all fermions and also for breaking the electroweak gauge symmetry spontaneously. The minimal model also has a pair of complex triplet scalars ∆ L,R in order to break the LRSM gauge symmetry spontaneously to that of the SM and also to generate Majorana mass terms of light and heavy neutrinos. The scalar content of the present model with the addition of a pair of real triplets, one bitriplet and one D-parity odd singlet, is shown in table 2. The parity (1, 1, 1, 0) Table 2: Scalar content of the model odd scalar not only helps in gauge coupling unification as we will discuss later, but also splits the masses of left and right handed dark matter candidates. At a very high energy scale, the parity odd singlet σ can acquire a vev to break Dparity spontaneously while the neutral component of ∆ R acquires a non-zero vev at a later stage to break the gauge symmetry of the LRSM into that of the SM which then finally gets broken down to the U (1) em of electromagnetism by the vev of the neutral component of Higgs bidoublet Φ. Thus, the symmetry breaking chain is Denoting the vev of the neutral components of the bidoublet as k 1,2 and that of triplet ∆ R as v R , the gauge boson masses after spontaneous symmetry breaking can be written as where θ w is the Weinberg angle. The neutral components of the other scalar fields ψ, Ω L,R do not acquire any vev. However, the neutral component of the scalar triplet ∆ L can acquire a tiny but non-zero induced vev after the electroweak symmetry breaking as with M W L ∼ 80.4 GeV being the weak boson mass and γ is a function of various couplings in the scalar potential. The bidoublet also gives rise to non-zero W L − W R mixing parameterised by ξ as which is constrained to be ξ ≤ 7.7 × 10 −4 [57,58]. Such tiny left-right mixing simplifies the calculation of dark matter relic abundance as we discuss in detail later. The relevant Yukawa couplings for the standard model fermion masses can be written as C is the charge conjugation operator and the indices i, j = 1, 2, 3 correspond to the three generations of fermions. The light neutrino mass after the spontaneous symmetry breaking can be written as where (M D ) ij = y ij k 1 +y ij k 2 are the elements of Dirac neutrino mass matrix and (M R ) ij = (f R ) ij v R are the right handed neutrino mass matrix elements. The fermion triplets Σ L,R shown in table 1 can be written in the component form as the following matrix representation where the neutral component of the each fermion triplet can be a stable dark matter candidate whose stability is ensured by either high SU (2) dimensions forbidding any decay at tree level or due to remnant discrete symmetry arising after spontaneous symmetry breaking of LRSM [59]. There exists a parity odd singlet scalar σ that can introduce a mass splitting between Σ L and Σ R . The terms involving Σ L,R and σ can be written as A non-zero vev of the parity odd singlet scalar σ lead to the masses of fermion triplets as This generates a mass splitting of 2λ σ σ between the two dark matter particles. Since we want our heavier dark matter candidate Σ R to decay into the lighter one and other standard model particles, we introduce a scalar bitriplet into the model which can mediate such decays. Such a scalar bitriplet field ψ ∼ (1, 3, 3, 0) can be written in matrix form as (2.8) The relevant interaction Lagrangian is given by where the contribution of the parity odd singlet in the desired mass splitting M R = M L is already taken into account. A pair of scalar triplets Ω L,R is introduced in order to achieve the desired relic abundance of the heavier dark matter candidate as we discuss below. These scalar triplets couple to the fermion dark matter triplets as The complete Yukawa Lagrangian of the model therefore, can be written as , (2.9) and (2.10) respectively. For the sake of completeness we also write down the scalar potential of the model by dividing it into two parts: one for minimal LRSM and the other for the relevant part of the additional scalar content (Ω L,R , ψ, σ). The scalar potential for the minimal LRSM is where the bilinear terms in Higgs fields are The self-interaction terms of Φ are: and the ∆ L,R self-and cross-couplings are as follows: In addition, there are also Φ − ∆ L and Φ − ∆ R interactions present in the model, The scalar potential involving the newly introduced scalar fields beyond the minimal LRSM is The details of different terms on the right-hand side of the above equation can be written as follows, It should be noted in the scalar potential terms involving the scalar bitriplet, both the two and three dimensional representations of SU (2) generators have to be used at appropriate places in order to construct an invariant term. The details of this is skipped as these terms are shown for the sake of completeness here and are not relevant for subsequent discussions. The complete scalar potential of the model is, therefore On the other hand, the kinetic Lagrangian of the model can be written as In L Scalar kin the covariant derivatives D µ for the scalars can be written in a way similar to that of the fermions.

Gauge Coupling Unification
In the previous discussion, we demonstrated that PeV scale decaying dark matter has a potential to explain the recently observed IceCube data. We intend here to examine whether such a framework under consideration can be embedded in a non-SUSY SO(10) GUT theory leading to successful gauge coupling unification. The symmetry breaking pattern of SO(10) such that it has the Pati-Salam gauge group as its intermediate symmetry breaking step as follows The important stage of symmetry breaking G 224 → G SM is happened when righthanded Higgs field Ω R (1, 3, 15) gets its non-zero vev. This symmetry breaking scale is fixed at few PeV scale such that decaying dark matter Σ R , right-handed charged gauge boson W R and relevant scalars playing an important role in dark matter phenomenology and IceCube explanation -all lie around that scale. The last stage of symmetry breaking G SM → U (1) em × SU (3) C is done via SM Higgs doublet reproducing fermion masses and mixing.
Here both manifest left-right theory with gauge group SU  The one-loop renormalisation group equations (RGEs) for gauge couplings g i for dif-  ferent stages of symmetry breaking scales is given by where the classic formula for one-loop beta-coefficients 1 b i are given by The parameters used in the above formula for one-loop beta coefficients have their usual meaning as follows representation, • T (R f ) and T (R s ): traces of the irreducible representation for a fermion and scalar R f,s , respectively, (3.5) • d(R f,s ): dimension of a given fermion or scalar representation R f,s under non-abelian gauge group SU (N ) except the considered i-th gauge group.
• Extra factor of 1/2 should be multiplied if the Higgs belongs to a real representation. The

Dark Matter
In this model, we have two dark matter candidates namely the neutral components of fermion triplets Σ L,R . We consider the right handed DM to be at PeV scale while the left handed one can be as light as a TeV. Usually, the relic density of Cold Dark Matter (CDM) in the Weakly Interacting Massive Particle (WIMP) scenario is calculated using the formalism discussed in ref. [69] which can be generalised to our present case as already discussed in [46,59,70] as However, for generic thermal relic dark matter there exists an upper bound on the dark matter mass (M DM ≤ O(100) TeV) coming from the unitarity limit on the annihilation cross section [49]. However, such a limit on dark matter mass can be relaxed if dark matter annihilation receives a Breit-Wigner enhancement [49,71]. As we will see below, the usual gauge interactions of the right handed fermion dark matter are not sufficient to give such an enhancement. Even if we include a separate scalar mediated annihilation channel, it is not sufficient the avoid the unitarity bound as we discuss below. The total abundance is constrained by the observed number mentioned in Eq. (1.1) given by the observations from the Planck experiment [4]. We can calculate the relic abundance of both the DM candidates independently in the limit of negligible left-right mixing of gauge bosons and negligible coannihilations between the two DM candidates mediated by the bitriplet scalar. The assumption of negligible coannihilation is justified Figure 2: The annihilation and coannihilation diagrams for fermion triplet dark matter relevant for relic density computation within the framework of minimal left-right dark matter [59].
in the limit of large mass splitting between the two DM candidates as we discuss below. In this decoupled limit, the individual relic abundance can be calculated by solving the corresponding Boltzmann equation where n χ is the dark matter number density and n eqb χ is the corresponding equilibrium number density. H is the Hubble expansion rate of the Universe and σv is the thermally averaged annihilation cross section of the dark matter particle χ. The approximate analytical solution of the above Boltzmann equation leads to the following expression for density parameter [72,73] where x F = m χ /T F , m χ is the mass of dark matter particle, T F is the freeze-out temperature, g * is the number of relativistic degrees of freedom at the time of freeze-out and M Pl ≈ 10 19 GeV is the Planck mass. Here a, b are the coefficients of partial wave expansion σv = a + bv 2 . The freeze-out temperature or x F can be calculated iteratively from the following relation which follows from equating the interaction rate with the Hubble expansion rate of the Universe. The thermal averaged annihilation cross section σv is given by [74] σv = 1 8m 4 where K i 's are modified Bessel functions of order i. Although we are assuming negligible coannihilations between the two DM candidates justified by their large mass splitting, there can be sizeable coannihilations within each DM multiplet, between the neutral and the charged components. At the classical level, the mass splitting between the charged and the neutral components of each DM multiplet is zero. However, a non-zero mass splitting can arise at one loop level given by [59,70], where s M = sin θ M ≡ tan θ W , r X = M X /M and Here the one loop self-energy corrections through mediations of gauge bosons are presented within the square bracket of the second expression. For example, the mass splitting with the approximation M Σ M W R goes as α 2 M W R − c 2 M M Z R /2. Due to such tiny one loop mass splitting, the next to lightest component of each DM multiplet can be thermally accessible during the dark matter freeze-out. In such a situation, the dark matter can coannihilate with the heavier components which then affects the relic abundance of dark matter. In the presence of coannihilations, the effective cross section is given by [75] The thermally averaged cross section can be written as Since the effective degrees of freedom g eff decreases exponentially with increasing mass splitting between DM and heavier components, one can ignore such coannihilation effects for scenarios with large mass splitting. The effective annihilation cross-section for the fermion triplet can be written as is the mass splitting ratio and x = M Σ 0 A /T . Here A = L, R denotes the dark matter candidates.
Considering the heavier dark matter Σ 0 R to be in the PeV regime we calculate its relic abundance by considering only gauge interactions. The gauge boson mediated annihilation and coannihilation channels are shown in figure 2. However, it was found for PeV right handed fermion dark matter that the gauge boson mediated channels were not sufficient to provide the resonance enhancement necessary to produce correct relic abundance of such superheavy dark matter candidates, as clearly seen from figure 3. We then incorporate the presence of scalar triplet Ω R through which the dark matter can (co)annihilate into a pair of light scalars from the bidoublet Φ. One such possible process is the coannihilation of Σ 0 R , Σ ± R into Φ 0 , Φ ± with s-channel mediation of Ω ± R . Taking degenerate masses of Σ 0 R , Σ ± R as well as Φ 0 , Φ ± for illustrative purposes , the corresponding annihilation cross section is given by where λ, f S are dimensionless couplings of Ω R with Σ R and scalar bidoublet Φ respectively. As seen from the scalar potential in Eq. (2.22), f S can be identified with f 15 . v S is the vev of the neutral component of Ω R which is taken to be similar to v R . m Ω ± R , M h are the masses of Ω ± R , Φ respectively. Since Ω R does not couple directly to the SM fermions, it may be possible to tune its decay width Γ Ω ± R in such a way to give the required resonance enhancement to dark matter annihilations. This possibility was utilised in the previous work [27] to bring the dark matter abundance into the observed range. However, as pointed out by the same authors later [76,77], one can not tune the decay width arbitrarily to enhance the cross section and hence the unitarity bound on the dark matter mass still prevails, in agreement with an earlier work [50]. This can be understood by approximating the above annihilation cross section for the resonance just above the threshold as , BR Φ 0 ,Φ ± being the corresponding branching ratios. Therefore, even if we adjust the resonance, the annihilation cross section decreases with the mass of dark matter and hence it again leads to a unitarity bound M Σ 0 R ≤ 100 TeV [49]. The authors of [50] found a slightly stronger bound M Σ 0 R ≤ 20 TeV from the unitarity requirement. Thus for PeV scale heavy dark matter, we will overproduce it at least by a factor of around 10 6 , if the unitarity bound on the annihilation cross section is respected. This problem of producing the correct thermal relic abundance of PeV DM is also summarised in the plot shown in figure 4. It shows that even near the resonance, the annihilation cross section remains much below the required one for producing the correct abundance. On the other hand, the required annihilation for correct relic abundance also lies above the maximum allowed cross section of DM from unitarity arguments.
The thermally overproduced PeV dark matter in our model can be reconciled with the observed dark matter relic abundance if there exists a long lived particle that decays into the standard model particles after the PeV dark matter freezes out. Such a decay should however occur before the big bang nucleosynthesis (BBN) temperature T BBN ∼ O(MeV) in order to be consistent with successful BBN predictions. Such late decay of long lived particles can release extra entropy and dilute the abundance of PeV dark matter to bring it into the observed limit [51]. The dilution factor due to the decay of such a heavy long lived particle N is given by [51] where Γ N is the decay width of the heavy particle with mass M N and Y N = n s is the initial abundance of the particle N before it started to decay. Also, g * (T r ) is the relativistic degrees of freedom at a temperature T r just after the decay of N . This temperature to which the Universe cools down to following the release of entropy due to the decay of N can be approximated as Although there can be several different possibilities of long lived particle in our model, we consider the late decay of Ω R into the standard model particles.
We can now find the decay diagrams through which the lightest component of Ω R can decay. If the neutral component of Ω R does not acquire any vev, and M Ω R < M Σ 0 R , then the lightest component of Ω R can decay only at one-loop into the standard model particles. We consider the lightest component of Ω R to be the neutral particle which can decay into a pair of photons through Σ ± R in loop. The Decay of the neutral component of the triplet Ω R to two photons are given as 14) (4.15) where t = M 2 Σ /m 2 Ω and λ is the dimensionless couplings of Ω R with Σ R appearing in Eq. (2.10). We consider the initial abundance of the decaying particle Here T f,Ω R is the freeze-out temperature of Ω R which is assumed to take place much before it starts decaying so that between freeze-out 10. 15 10. 14 10. 13 10. 12 10. 11 10. 10  and decay, the Ω R particles are neither being created nor destroyed. If we assume the freeze-out to be above the electroweak scale, then g * (T f,Ω R ) = 106.75 corresponding to all the SM degrees of freedom. It should be noted that here we are assuming equality between g * s appearing in entropy density expression and g * here which is valid at high temperatures where all species are in equilibrium with each other. We also take g * (T r ) = 10.75, assuming that the decay is taking place very late, towards the end of QCD phase transition. This also automatically satisfies the criteria that the decay should take place after PeV dark matter freezes out. We now demand the dilution factor to be around 10 −7 − 10 −5 , the decay of the heavy particle to take place before BBN that corresponds to approximately 1 second after the big bang. These criteria constrain the parameter space in terms of Ω R mass and its coupling to dark matter, i.e. λ, as seen from the plot shown in figure 5. The upper limit of Ω R mass comes from the requirement of forbidding tree level decay into Σ ± R pairs. It can be seen from the plot in figure 5 that the requirement of producing such a large entropy at late epochs, requires very long lifetime or equivalently, small decay width of the decaying particle forcing the coupling λ to be very small in spite of one-loop suppression. Such fine-tunings can be relaxed to some extent if we consider Ω ± R as the lightest component of Ω R which can decay into a photon and an off-shell W ± R .
Naturalness of parameters: In the above calculation, we have seen that the coupling between fermion dark matter candidates Σ L,R and the scalar triplets Ω L,R has to be very Table 4: Fields having non-trivial transformation under U (1) X global symmetry small in order to have the correct entropy dilution mechanism at work or equivalently to make sure that Ω R is long lived enough to release the required entropy during late epochs. Apart from this, one also requires additional fine tunings in order to make the heavier dark matter long lived and its decay into leptons plus the lighter dark matter more dominant over other decay channels. The additional scalar multiplets Ω L,R , ψ do not couple to the SM fermions, but they can still decay fast into the SM Higgs by virtue of scalar couplings shown in Eqs.  table 4 transform non-trivially. Such a symmetry prevents the fast decay of Ω 0 R , Σ 0 R , both of which we have considered to be long-lived in this work. Clearly, such an approximate global symmetry can naturally suppress the Yukawa couplings like f ψ , λ appearing in equations (2.9), (2.10) respectively as these trilinear terms violate the U (1) X global symmetry explicitly. Similarly, the trilinear mass term µ 16 appearing in equation (2.24) can also be naturally small, preventing the fast decay of Σ 0 R → Σ 0 L ΦΦ. One can also have decays like Σ 0 the neutral components of the scalar bitriplet acquire non-zero vev's. However, one needs to keep the vev of the bitriplet scalar small as the constraints on the electroweak ρ parameter restricts it to v ψ ≤ 2 GeV [78]. In the standard model, the ρ parameter is unity at tree level, given by where θ W is the Weinberg angle. Any contribution to the masses of electroweak gauge bosons from scalar multiplets apart from the scalar doublet will give rise to deviation from the above tree level expression and hence will be constrained from electroweak precision measurements. The tiny vev of the neutral components of ψ can be guaranteed if ψ also transforms non-trivially under the approximate global symmetry U (1) X , as shown in the table 4. Such symmetry will naturally explain the vanishingly small trilinear mass term µ 16 in Eq. (2.24) which can automatically guarantee a tiny vev to the neutral component of bitriplet scalar ψ. From the minimisation of the scalar potential with respect to ψ, it is straightforward to find the induced vev to be where µ 2 ψ > 0 is the bare mass squared term of the bitriplet scalar. Therefore, a tiny mass term µ 16 will naturally ensure the smallness of the vev. Such an approximate symmetry will also give rise to a small trilinear term µ 236 Eq. (2.24) which will appear in the decay of heavier DM to lighter DM as we discuss in the next section. Therefore, an approximate global U (1) X symmetry can suppress the fast decay of heavier DM into the lighter one and a pair of Φ or W ± L or Z L bosons. Similarly, one can also have fast decay of Ω 0 R into the standard model fields if it acquires a non-zero vev. In the discussion on thermal abundance of heavier dark matter, we considered non-zero vev of Ω 0 R in order to allow the s-channel coannihilation into light scalars. However, since that was not sufficient to produce the correct thermal abundance, we abandon that setup and have moved on to discussing the decay of Ω 0 R as mentioned above. In fact, a non-zero vev of Ω 0 R will allow it to decay very fast into the SM fermions mediated by W R bosons. Even if we assume a zero vev of Ω 0 R by assuming its bare mass squared term to be positive definite, it can still acquire a non-zero induced vev by virtue of its coupling with the scalar bidoublet shown in Eq. (2.22). This induced vev can be found from the minimisation of the scalar potential as where µ 2 Ω > 0 is the bare mass squared term of the triplet scalar Ω R . The same trilinear interaction involving µ 15 also induces the decay of Ω 0 R into a pair of SM Higgs bosons. In order to prevent Ω 0 R from acquiring a non-zero vev and from decaying into a pair of SM Higgs bosons we need to prevent its trilinear interactions with the bidoublet or make those interactions very small. This is possible by virtue of the non-trivial charge of Ω L,R under the approximate global symmetry, as shown in table 4. In particular, the trilinear mass terms µ 14 , µ 15 in Eq. (2.22) can be vanishingly small as they explicitly break the U (1) X global symmetry, in accordance with the naturalness criteria. Therefore, the presence of an approximate global symmetry U (1) X can naturally give rise to long lived heavier DM, long lived Ω 0 R required to release entropy at late epochs. Although such an approximate symmetry can naturally explain the long life of these two particles, it however can not explain the dominance of leptonic final state decay mode of heavier dark matter over other decay modes with non-leptonic final states like a pair of Φ or W ± L or Z L bosons in the final state. Since such approximate global symmetries are likely to originate from an ultraviolet complete theory, we leave this question to the details of such a theory at high energy scale.
We also calculate the relic abundance of left fermion triplet dark matter. Since the mass splitting between different components of the left fermion triplet is a function of electroweak gauge boson mass, the only free parameter in this case is the mass of dark matter. We show the variation of its relic abundance with mass in figure 6. It can be seen  Figure 7: Feynman diagrams for next to lightest dark matter particle decaying into lightest dark matter particle and two light neutrinos.
from the figure that correct relic abundance can be generated for left fermion triplet mass around 2.5 TeV. For heavier masses, the relic is overproduced whereas it is under-abundant for lower masses. Since the left fermion triplet is the lighter dark matter candidate in our model and hence cosmologically stable, it is important to make sure that it does not get overproduced. In our subsequent analysis, we will assume the heavier dark matter to be the most dominant component which restricts the left triplet dark matter masses to be smaller than 2.5 TeV. Such TeV-scale mass for left fermion triplet is also taken into account in the RGE analysis discussed in the previous section.
Apart from the total relic abundance criteria to be satisfied by the two dark matter candidates together, the heavier dark matter decay should be consistent with the observed neutrino flux at IceCube experiment. The amplitude of the heavier dark matter decay diagram 7 is given as, where with f ψ , f L being Yukawa couplings appearing in Eqs. (2.9), (2.4) respectively, θ being the scalar mixing parameter, and M S being the mass of the neutral physical scalar as mediator. The mixing between the scalars ψ, ∆ L can be given by the angle θ as where µ 236 is the trilinear mass term appearing in Eq. (2.24). After doing further calculation, we get the following differential decay width, where Integrating the differential decay width over the neutrino energy gives the following decay width of dark matter, Assuming M S M Σ 0 R , we get the neutrino energy distribution (for massless final states), (4.20) The normalised energy distribution of the primary neutrinos is shown in figure 8 for three different values of heavier dark matter masses. From figure 8 it is clear that PeV events at IceCube can be described by DM decay with its mass of about 5 PeV.

IceCube PeV Events
The 1347-day IceCube high energy starting events (HESE) neutrino data shows 54 events with deposited energies ranging between 30 TeV to 2.1 PeV [8,10,79]. That the origin of the events above a few tens of TeV is non-atmospheric and extra-terrestrial in nature has been established with a high level of significance [8]. The three highest energy events of deposited energies 1.04 PeV, 1.14 PeV and 2.0 PeV are all shower events [9]. In addition, recently a track-like event with deposited energy ∼ 3 PeV is reported by IceCube [11]. Its explanation as a signal from a new astrophysial neutrino flux is explored in [80], suggesting a window into astrophysical neutrinos at E ν ∼ 100 PeV. Moreover, the highly sought after Glashow resonance, with the energy around 6.3 PeV, has not yet been observed 2 . A single power-law fit to the IceCube flux in the expected E −2 spectrum, coming from Fermi shock acceleration considerations, predicts more number of multi-PeV events, it is thus disfavoured.
It has been shown that the scenarios with superheavy decaying dark matter can ameliorate some of the previously mentioned tensions [14-25, 27, 28]. In the present model that we are considering we have a two-components (Σ 0 L and Σ 0 R ) DM scenario. The mass of the heavier one among these two species, the Σ 0 R is in PeV range. We want to test its feasibility as a possible candidate which can explain the IceCube PeV events through its three-body decay as Σ 0 R → Σ 0 L ν ν , where ≡ e, µ, τ . Clearly, the primary production of neutrinos will contribute in the incoming neutrino flux. The other decay channels of Σ 0 R includes Σ 0 R → Σ 0 L , where ≡ e, µ, τ . It is clear that in this model, the secondary production of neutrinos is possible through the decays of primary leptons (µ and τ ) which contributes in the lower energy range of the event distribution at IceCube. It is important to mention that the astrophysical contribution provides a good fit to the low-energy data and any additional contribution is tightly constrained from IceCube data [10]. In our considered model, by suppressing the decay of Σ 0 R into quarks and bosons, it is possible to avoid these constraints.
The neutrino flux coming from the heavy DM decay has two components, namely the galactic and extra-galactic components. In the subsequent, we briefly elaborate the standard method to calculate this flux following the procedure prescribed in [15,19,84].
The galactic component comes from the decay of DM in the Milky Way halo and the differential flux corresponding to this is given by, The absence of the Glashow resonance is discussed in details in refs. [81][82][83] where m Σ 0 R and τ Σ 0 R denote the mass and lifetime of the DM particle Σ 0 R and dN (E ν )/dE ν is the energy spectrum of neutrinos produced in the decay of the DM. We compute dN (E ν )/dE ν for primary as well as secondary neutrinos using the methods outlined in ref. [85,86]. For the considered spectrum we implicitly assumed the sum over all neutrino and anti-neutrino flavors. The quantity ρ Σ 0 R (r) is the density profile of DM particles in our galaxy as a function of distance from the galactic centre. In our calculations we take the density profile to be of the NFW form [87] ρ Σ 0 R (r) = ρ 0 (r 0 /r)/(1 + r/r 0 ) 2 with ρ 0 = 0.33 GeV cm −3 and r 0 = 20 kpc. The integral over the parameter s is basically the line-of-sight integral to obtain the flux at the Earth and s is related to r via the relation r(s, l, b) = s 2 + R 2 − 2sR cos b cos l, where R 8.5 kpc is the distance of the Sun from the galactic centre; l and b represent the longitude and latitude in galactic coordinates.
The second component, i.e., the extra-galactic component comes from the DM decays at cosmological scales and this produces an isotropic diffuse flux which in its differential from is given by, where ρ c = 5.5 × 10 −6 GeV cm −3 is the critical density of the universe and the Hubble expansion rate as a function of the redshift z is given by H(z) = H 0 Ω Λ + Ω m (1 + z) 3 . For the considered parameters we use the results from Planck collaboration [4] which, following the standard ΛCDM cosmology gives Ω DM = 0.27, Ω Λ = 0.68, Ω m = 0.32, and h ≡ H 0 /100 km s −1 Mpc −1 = 0.67. The total neutrino flux coming from DM decay is, can be written as the sum of galactic and extra-galactic contributions and given as, For completeness we also take into account the astrophysical neutrino flux, the contributions of which comes from extragalactic supernova remnants (SNR), hypernova remnants (HNR) [88], AGN [89,90], and GRB [91]. In order to parametrize the astrophysical flux, we assume a single unbroken power-law astrophysical flux: where Φ 0 = 2.2 × 10 −8 GeVcm −2 sec −1 sr −1 and γ = 0.58 which correspond to the best-fit value at Icecube [10] assuming equal composition of neutrino flavor on Earth 3 . The total neutrino flux will thus be the sum of the astrophysical flux and the flux coming from DM  decay and this can be written as, After getting the total neutrino flux, it is now straightforward to compute the number of events and compare with the IceCube data. To this end we define the quantity N i which gives the number of events in a given deposited energy bin [E i , E i+1 ] as follows, where t is the exposure time which is 1347 days, A(E ν ) is the neutrino effective area [8].
To compare the prediction of (5.6) with the IceCube observations, one needs to transmute the number of events in the terms of deposited energy. But as discussed in [24] due to low statistics collected till now, it is acceptable to assume that the two energies coincide. It is important to mention that in the calculation of N i we have summed over all neutrino flavors. For a characteristic benchmark point with m Σ 0 R = 5 PeV and τ Σ 0 R = 1.2 × 10 29 sec, we present the event distribution at the IceCube in figure 9. From this figure it can be seen that the contribution in the neutrino flux coming from the decay of the DM, Σ 0 R can explain the excess in the observed IceCube PeV neutrino events and the subsequent cut-off in the distribution. The lower energy region of the distribution is, however, in agreement with the astrophysical and atmospheric neutrino flux which constitutes the main background events and the contribution in the flux from the DM decay is negligible in this region.

Summary and Conclusion
Motivated by the possibility of dark matter interpretation of the IceCube high energy neutrino events we study a non-supersymmetric left-right model that not only explains the IceCube events from a heavy PeV dark matter decay, but also achieves gauge coupling unification at high energy scale with SO(10) embedding. Although PeV scale heavy dark matter does not fall in the typical weakly interacting massive particle ballpark due to the unitarity bound on dark matter mass, our model can evade such bounds through resonance enhancement of dark matter annihilations. Also, the particle content of the model is chosen in such a way that predicts a PeV scale intermediate SU (2) R × U (1) B−L symmetry, suggesting a fundamental reason behind the origin of PeV scale dark matter, apart from its motivation from IceCube point of view.
The model has a pair of left and right fermion triplet dark matter candidates. Their equal masses in the limit of exact left-right symmetry acquire a large splitting due to the breaking of discrete left-right symmetry at a very high energy scale, resulting in a PeV scale heavy right handed fermion dark matter and a light TeV scale left handed fermion dark matter. The PeV scale dark matter can decay into the lighter one along with a pair of neutrinos on cosmological scale so that the high energy IceCube events can be explained. The large mass splitting between left and right fermion dark matter keeps their coannihilations sub-dominant. Also, in the limit of negligible left-right mixing, the annihilation of heavier dark matter into lighter ones can be ignored. This simplifies the dark matter relic abundance calculations as it can be calculated separately for individual dark matter candidates. The individual dark matter multiplets can however, have coannihilations within its components due to tiny mass splittings that arise at one loop level. We first calculate the heavier dark matter abundance by incorporating these coannihilations through right handed gauge bosons and show that these gauge interactions are not sufficient to give the required resonance enhancement to dark matter annihilations in order to bring the relic abundance down to the observed one. We then show that even if we introduce another scalar to give further s-channel enhancement, it is not possible to tune its decay width arbitrarily in order to generate the required annihilation cross section. The unitarity limit on the cross section forces us to restrict the dark matter mass to be at most 20 TeV. This overproduces the heavy dark matter by a factor of 10 5 − 10 6 compared to the observed abundance. We then show one possible way of diluting this abundance by late decay of a heavy scalar which releases entropy through its decay before BBN and dilutes the thermally overproduced PeV dark matter by appropriate amount. Here we assume the heavier dark matter to be the dominant dark matter component to avoid any suppression in its decay flux into the high energy neutrinos. The lighter dark matter relic abundance can be easily kept at a suppressed value by choosing its mass suitably in the TeV regime. Unlike the PeV right handed dark matter, the TeV left handed dark matter annihilation can be sufficiently enhanced to keep the relic abundance suppressed even with electroweak scale gauge bosons as mediators.
The decay of PeV right handed dark matter to the TeV left handed dark matter with a pair of neutrinos can explain the high energy IceCube events. Using 4-yr IceCube data set of 54 events we show that together with astrophysical neutrino flux for lower energy events and the contribution in neutrino flux coming from the decay of the PeV right handed dark matter can fit the IceCube data from TeV-PeV energy range. We illustrate the fitting with a benchmark mass of the 5 PeV and lifetime of ∼ 10 29 s. Clearly, future data with more statistics and precision will put more restrictions on this model. In passing it is worth mentioning that our benchmark evades the recent strong constraints that normally come from the gamma ray observations. In the present scenario, PeV right handed dark matter has a three body decay into ν, µ and τ leptons which give rise to a soft spectrum for photon in comparison to the usual two body decay. As a result, bounds from this gamma ray observations [95] do not have one to one correspondence with the present scenario for τ ∼ 10 29 s and they can be exempted safely.
It is worth mentioning that, in order to have a long lived PeV scale dark matter candidate along with another long lived particle in order to generate the correct abundance of heavier DM by the mechanism of entropy dilution, we need to fine-tune several parameters of the model. In this minimal setup we have incorporated an approximate U (1) X global symmetry for this purpose. Another possible way is to go beyond this minimal setup in such a way that the decay of these long lived particles occur only at radiative level, naturally suppressing the decay widths. We leave such a study to future works.