A comparative study of the S 1 and U 1 leptoquark effects in the light quark regime

We study the phenomenology of two leptoquarks, the Standard Model SU ( 2 ) singlets S 1 and U 1 , with regard to the latest experimental data from the low-energy ﬂavor physics measurements, LHC, and the IceCube neutrino experiment. We consider a scenario when scalar (vector) leptoquark S 1 ( U 1 ) couples exclusively to the down quark and the neutrinos (charged leptons) of all ﬂavors, where the leptoquark in question couples to the SM lepton doublets. The couplings of S 1 ( U 1 ) to the up-type quarks and the charged leptons(neutrinos)areinturnuniquelydeterminedvia SU ( 2 ) symmetry. We ﬁnd that the most important constraints on the leptoquark parameter space originate from ﬂavor physics measurements, followed by the LHC search limits that take over the ﬂavor physics ones in the large LQ mass regime. We furthermore show that S 1 ( U 1 ) marginally improves (spoils) the ﬁt of the current IceCube data with respect to the SM case within the region of parameter space that is otherwise consistent with various low-energy ﬂavor physics measurements and the latest LHC input. Our study offers an up-to-date analysis for these two leptoquarks in view of the latest experimental data.

We are interested, in this manuscript, in the scenarios when LQs primarily couple to the quarks of the first generation and neutrinos of all flavors and investigate the viability of the associated parameter space spanned by the LQ masses and coupling strengths in view of the latest experimental data from flavor physics, LHC, and the South Pole situated IceCube detector. We accordingly study the implications of the most recent and the most relevant experimental results on the parameter space for two representative LQ scenarios. One scenario features scalar LQ S 1 and the other uses vector LQ U 1 , where both fields are singlets under the SM SU (2) group, allowing them to couple to both the left-and the right-handed quarks and leptons. We perform, in particular, a thorough analysis of the viability of the S 1 (U 1 ) scenario assuming non-zero couplings between S 1 (U 1 ), down quark, and neutrinos (charged leptons) of all three generations, where the LQ in question couples to the SM lepton doublets. Consequentially, S 1 (U 1 ) couples up-type quarks to charged leptons (neutrinos). We also entertain the possibility that S 1 (U 1 ) couples down quark (up quark) to the righthanded neutrinos to investigate the sensitivity of the latest IceCube data to constrain the associated parameter space.
The outline of the paper is as follows. We describe the two LQ scenarios and the flavor ansatz considered in our work in Sect. 2. The constraints from the low-energy flavor physics experiments for these two LQ scenarios are presented in Sect. 3. The LHC constraints from the single LQ and the LQ pair productions are discussed in Sect. 4. We then perform the data analysis of the IceCube PeV events within these two frameworks in Sect. 5. The combined analysis using the lowenergy flavor observables, along with the LHC results and the latest IceCube data, for both S 1 and U 1 , is presented in Sect. 6. Finally we conclude in Sect. 7.

Leptoquark scenarios
We briefly review in this section the LQ scenarios we consider in our work. The two representative scenarios that are addressed in our analyses are the scalar LQ S 1 and the vector LQ U 1 .

Scalar leptoquark S 1 = (3, 1, 1/3)
We study the signatures of S 1 , whose SU (3)× SU (2)×U (1) quantum numbers are (3, 1, 1/3), on the flavor, LHC, and IceCube observables. In our normalisation the electric charge of S 1 is 1/3 in the absolute units of the electron charge. The relevant Lagrangian terms, in the mass eigenstate basis, are of the form where the subscripts i, j (= 1, 2, 3) denote the flavor of the quarks and leptons, V is the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix and U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix. We work under the assumption that the only non-zero S 1 couplings are y L 11 ≡ y L dν e , y L 12 ≡ y L dν μ , and y L 13 ≡ y L dν τ . We also entertain the possibility that S 1 couples to the right-handed neutrinos and we set the couplings y R 11 ≡ y R dν e , y R 12 ≡ y R dν μ , and y R 13 ≡ y R dν τ to be equal to each other, if and when switched on. All the other LQ Yukawa couplings are set to zero. Note that the S 1 couplings with the up-type quarks and charged leptons are fixed by the CKM mixing matrix.

Vector leptoquark
The relevant Lagrangian terms for the U 1 LQ, in the mass eigenstate basis, are We consider the scenario where U 1 only couples to the down quark and charged leptons of all three generations with We also analyse the possibility when the couplings χ R 11 ≡ χ R uν e , χ R 12 ≡ χ R uν μ , and χ R 13 ≡ χ R uν τ of U 1 with the up quark and the right-handed neutrinos are switched on and equal to each other. All other U 1 couplings are set to zero.

Low-energy constraints
The LQ interaction ansatz defined in the previous section can lead to leptonic decays of pseudoscalar mesons or flavor changing processes at both the tree and the one-loop levels. The LQ couplings to the first generation quarks and electron are strongly constrained by the atomic parity violation (APV) experimental results. The experimental upper bounds on the → γ decay branching ratios, with the LQ contribution coming in the loop, will also constrain the couplings of S 1 with the quarks and leptons. These branching ratios receive contribution from both the left-handed as well as the right-handed couplings of the quarks to the leptons. The upper limits on the lepton flavor violating decays of μ and τ leptons are obtained from various experiments with BR(μ → eγ ) < 4.2 × 10 −13 [27], BR(τ → eγ ) < 3.3 × 10 −8 , and BR(τ → μγ ) < 4.4 × 10 −8 [28] @ 90% C.L.. S 1 also contributes at the tree level to the rare flavor process D 0 → μ + μ − . The most recent measurement of this branching ratio comes from LHCb [29] and reads BR(D 0 → μ + μ − ) < 7.6 × 10 −9 . The same couplings contribute to the D 0 −D 0 oscillations. Following the study of the LQ effects in the D 0 −D 0 oscillations explained in detail in Refs. [1,30], we require that the LQ contributions are smaller than the current bounds on the D 0 −D 0 mixing amplitude.
Since the LQs, in our case, yield new contributions to → γ , the APV measurements, the rare meson decays, and the ratio of the leptonic decays of the pseudoscalar meson, we take into account all these constraints.
Lepton flavor violation in the pion sector The contribution of weak singlets S 1 and U 1 to the pion muonic decays is different from the pion electron decays due to the different values of e and μ couplings with the first generation quarks as well as the dependence on m e and m μ . The effects of this type can be exposed by the lepton flavor universality ratios R π e/μ and R π τ/μ , where  [31] (see Eq. (A.2). The measured ratio is R π τ/μ | exp = 0.1082 ± 0.0005 [32], while the SM value is found to be R π τ/μ | SM = 0.1088 ± 0.0002, using Eq. (A.2). We list the formulas for the branching ratios of the pion, D 0 meson, and the τ lepton for the S 1 and U 1 cases in Appendix A and specify numerical values of input parameters relevant for our analysis in Table 1. We furthermore summarize in Fig. 1 results of a randomized scan within the parameter space (m S 1 ∈ (300 GeV, 1.5 TeV), y L dν e , y L dν μ , y L dν τ ∈ (0.0, 0.8)) and (m U 1 ∈ (500 GeV, 2.5 TeV), χ L de , χ L dμ , χ L dτ ∈ (0.0, 0.8)) that takes into account the constraints from the pion sector, APV, the rare radiative decays → γ , and D 0 → μ + μ − decays. The plots in Fig. 1 show currently allowed parameter spaces of the lefthanded couplings of the down quark and the first generation leptons as a function of the LQ mass as well as the correlations between different left-handed couplings. We find that y L dν e and y L dν μ (χ L de and χ L dμ ) cannot be simultaneously large due to conflict with the current results from the low-energy sector in the S 1 (U 1 ) case. This can be clearly seen in the panels of the second column of Fig. 1. We therefore mainly work, in what follows, in the presence of the left-handed coupling of the down quark and the first generation leptons, with the other couplings being set to zero.
The flavor experiments constrain the parameter space of the vector LQ more tightly than that of the scalar one. The S 1 LQ also contributes at the loop level to the Z → ¯ decay amplitude, with S 1 and the up-type quarks running in the loop. The Z branching ratio to a pair of leptons has been precisely measured at LEP [32], thereby imposing constraints on the S 1 parameter space. We have used formula for the one-loop contribution of S 1 computed in Ref. [33] and found that the bounds on the LQ couplings from the Z leptonic branching ratio are not up to par with the other experimental constraints considered before. The loop level contribution to the Z → ¯ decay amplitude in case of U 1 is also negligible for the parameter space that survives the other low-energy flavor physics experiments.
We next discuss the LHC limits on the LQ masses and their couplings.

Constraints from the Large Hadron Collider
The LQ couplings to the quark-lepton pairs have been constrained through both direct and indirect searches in a number of collider experiments. Prior to the LHC era, the LQs were searched for at LEP [34], HERA [35,36] and Tevatron [37,38]. The LQs have been hunted for at the LHC mainly through pair production [39][40][41] but there are also several searches/recasts that rely on the single LQ production [42] as well as dilepton [43] and monolepton [26] Drell-Yan processes to generate constraints. For a summary of currently available bounds on the LQ masses and associated couplings from the LHC searches for various flavor final states see, for example, Refs. [44,45].
The LQs are pair-produced through either gg or qq fusion that is primarily dominated by the QCD interactions. There is also a Yukawa coupling contribution to the LQ pair production, corresponding to a t-channel process, with its amplitude being proportional to the product of the two relevant Yukawa couplings. This t-channel process is highly suppressed compared to the QCD driven one, unless the Yukawa couplings are rather large. The production cross section at the LHC for a pair of vector LQs, when compared to the scalar ones, additionally depends on the underlying theory for the origin of vector LQs. The relevant trilinear and quartic couplings of vector LQs to a gluon or a pair of gluons is then completely fixed by the extended gauge invariance of the model wherein the vector LQs appear as fundamental objects. We here work in the limit, where the vector LQ is some low energy manifestation of a more fundamental theory at high energy scale. The vector LQ-gluon interaction terms can then be obtained from the most general SU (3) invariant effective Lagrangian given by where G μν is the gluon field strength tensor and κ is a dimensionless parameter which we consider to be 1 for our calculations. Dedicated studies have been performed at the LHC, assuming LQ pair production and a 100% branching ratio (BR) of LQ decaying into a charged lepton and jet ( j j + − , where = e, μ) or to jet and missing energy ( j jνν). Considering the model Lagrangian in Eq. (2.1) (Eq. (2.2)) in the S 1 (U 1 ) case, the final states relevant for our analysis are j jee, j jμμ, and j jνν, with jν having the dominant branching ratio. There being no distinction made in the LHC LQ searches in case of the light quark jets (u, d, c, s), we can consider the LHC limit directly. The upper limits on the LQ production cross section times BR 2 for these final states are provided by the LHC collaborations [39,46,47]. The LQs coupling to first-generation quarks and electrons or muons are also sought in single production processes, i.e., pp → + − j [42]. This process occurs via s-and t-channel quark-gluon fusion and is directly proportional to the Yukawa coupling of the u(d) quarks to the leptons. The single production of LQs at 8 TeV LHC [42] is also considered in our analysis.
We show in Fig. 2 the allowed parameter space for two different choices of values of relevant couplings of S 1 (U 1 ) as a function of m S 1 (m U 1 ). In an effort to confront the LHC constraints with the flavor physics measurements, we consider two particular scenarios: The region below the purple line in Fig. 2 is currently allowed by the low-energy experiments at the 2 σ level, as discussed in Sect. 3. The parameter space allowed by the LHC, at 95% C.L., is the region below the brown line in case of LQ pair production and the subsequent decay to light quarks and leptons. The allowed region from the single LQ production is the area below the green (for qe + e − ) and blue (for qμ + μ − ) lines. In the single LQ production, the process pp → qe + e − is proportional to the Yukawa coupling y L dν e (χ L de ), whereas pp → qμ + μ − is proportional to y L dν μ (χ L dμ ) for the S 1 (U 1 ) case. Therefore no limit is obtained from the process pp → qμ + μ − , in case of . The low-energy experiments, as discussed in the previous section, do not allow simultaneous presence of large values for the couplings of the first generation quarks to electrons and muons. Therefore the choice y L dν e (χ L de ) = y L dν μ (χ L dμ ) is strongly constrained by the flavor physics measurements. The single LQ production has more stringent limits than the pair production for ee final state at large LQ masses. This is due to the final state phase space. The blue dotted line in Fig. 2c is incomplete as the experimental result from pp → qμ + μ − is provided up to the LQ mass of 1.8 TeV. The most stringent constraint on the available parameter space currently comes from the low-energy flavor experiments denoted by the purple line. The LHC direct searches from the LQ pair production currently exclude m S 1 < 700 GeV, and m U 1 < 1700 GeV, irrespective of the choice of couplings as long as LQ decays promptly. A study done in Ref. [26] has shown that stringent limits can also be obtained on the strength of the LQ coupling to the first generation quarks and left-handed electrons and muons through the monolepton searches.
We study next in details the effects of S 1 and U 1 on the observed IceCube PeV events.

PeV events in IceCube
The observation of the High Energy Starting Events above 100 TeV at the South Pole situated IceCube detector [48][49][50][51], consistent with a flux of high-energy astrophysical neutrinos from outside the galaxy, has motivated a large number of studies that explore the IceCube potential to test various NP models. Since some of the most studied NP sources are various LQ scenarios our intention is to investigate whether the latest IceCube data [51] might offer an independent way to constrain the LQ mass m LQ and the strength of its couplings to the quark-lepton pairs that would be on par with the flavor and collider physics constraints.
The high-energy neutrinos coming from outside the atmosphere are detected in the IceCube detector by observing the Cherenkov light emitted by the secondary charged particles produced in the interaction of the neutrinos with the nucleus present in the ice. The charged current (CC) and the neutral current (NC) interactions have distinctive topologies depend- where T is the exposure time in seconds comprising 2635 days of data taking [51], N A is the Avogadro's number 6.022×10 23 , k is showers and tracks for each channel ν = {e, μ, τ } induced by the charged and neutral current interactions (ch) for an incoming neutrino flux of type f (astrophysical (a), conventional atmospheric (ν) or prompt atmospheric flux ( p)). The effective mass of the detector, M eff (E k,ch true ), is a function of the true electromagnetic equivalent energy and is defined as the mass of the target material times the efficiency of converting the true deposited energy of the event into an observed signal. The energy resolution function is given by R(E k,ch true , E dep , σ (E k,ch true )) and is represented by a Gaussian distribution [52]. The effect of the earth's attenuation, in case of neutrino's energy above a few TeV, where the mean free path inside the earth becomes comparable to the distance travelled by the Fig. 3 The relevant Feynman diagrams for the neutrino-quark interactions mediated by the scalar LQ S 1 Fig. 4 The relevant Feynman diagrams for the neutrino-quark interactions mediated by the vector LQ U 1 neutrino, is denoted by Att ν (E ν ). The incoming neutrino flux is given by dφ f ν /d E ν , where the incoming astrophysical neutrino flux follows the isotropic single unbroken power-law spectrum. This spectrum is given by [51] dφ astro where f is the fraction of neutrinos of each flavor . The fit is performed assuming a (1/3 : 1/3 : 1/3) ⊕ flavor ratio, which yields the best fit value for the spectral index γ = 2.89 +0.20 −0.19 , with a normalization 0 = 6.45 +1.46 −0.46 × 10 −18 GeV cm −2 s −1 sr −1 at 1σ significance.
The neutrino-nucleon differential cross section for different channels in case of the CC and the NC interactions is given by dσ ch ν (E ν , y)/dy. The SM differential cross section is given by Eq. (B.1) in Appendix B. At the IceCube detector the neutrinos interact with the nucleons present in the ice. We assume that the natural ice nucleus can be treated as an isoscalar with 10 protons and 8 neutrons. We calculate the event spectra of showers and tracks for each flavor in case of SM assuming an isotropic power-law spectrum. Since we find that the largest contribution to the event spectra comes from the ν e showers, the electron neutrino should be sensitive to the NP effects if one is to have an enhanced effect compared to the SM. We therefore study next the effect of S 1 and U 1 on the IceCube spectrum when these LQs couple the first generation quark to the electron.
The scalar LQ S 1 mediates the NC interactions ν d → ν d and ν d → ν d and the CC interactions ν ū → d , ν ū → d , and dν → u , dν → u , where = . The Feynman diagrams for the relevant processes are shown in Fig. 3. The charm contribution towards the tchannel CC process depicted in Fig. 3(iii), due to small PDFs, is maximally around 0.001% for the choice of the mass and the couplings considered here and is therefore neglected. The differential ν j N cross sections, in the presence of the S 1 interactions, are given by Eq. (B.8) in Appendix B.1.
The vector LQ U 1 contributes to both the NC and CC interactions. The relevant Feynman diagrams, in the presence of U 1 , are shown in Fig. 4. The u i in Fig. 4 represents the contributions from all three generations of up-type quarks. Note, however, that the charm contributions towards the NC processes, due to small PDFs, are of the order of 0.001% and can be safely neglected. The differential ν j N cross sections in the presence of the U 1 interaction are given by Eq. (B.10) in Appendix (B.2). The U 1 LQ compared to the S 1 case interferes with the SM leading to interesting features. Since the U 1 LQ interferes with the SM contribution, we show in Fig. 5 the ratio of the ν e N total cross sections for the SM + U 1 and the SM for different values of masses and couplings. The interference effect is clearly visible for low values of mass and large values of χ L de . The seven years of IceCube data have fewer events when compared to the SM in the 200-300 TeV energy range whereas for energies above 1000 TeV there are more events when compared to the SM. It can be seen from Fig. 5 that there is a crossover in the relevant energy range making it an interesting feature for a more detailed study. We would like to point out that the inclusion of χ L dμ and/or χ L dτ will push the crossover away from the interesting energy range. The three couplings then have to be adjusted so as to get the required effect.
We study whether the SM + LQ scenarios result in a better or a worse fit of the IceCube data compared to the SM case by calculating parameter δ that corresponds to the percent change in χ 2 [53]. We accordingly define Fig. 5 The ratio of ν e N total cross sections between the SM + U 1 and the SM for different masses and couplings of U 1 Fig. 6 The contour lines of constant δ in the m S1(U1) -y L dνe (χ L de ) plane showing the percent change of χ 2 compared to the SM case. The region to the right of the dotted curve in case of S 1 and above the dotted curve in case of U 1 results in a fit comparable or better than the SM one. The region above the blue line is excluded at the 2 σ level by the limits from the APV experiments where the observed number of events N data i in each bin i is compared with the LQ scenario prediction and, in our case, LQ = S 1 , U 1 . We consider the events in the neutrino deposited energy range [60 TeV, 10 PeV] that is divided in 20 logarithmic energy bins. We initially use only the data for the bins with the non-zero number of events. The SM value with the best fit value of γ and C 0 from the IceCube data results in a χ 2 value of 0.15. This shows that the current IceCube data is quite compatible with the SM. The NP contribution to the number of events in each bin depends on the values of y L dν e ,dν μ ,dν τ , y R dν , and m S 1 (χ L de,dμ,dτ , χ R uν , and m U 1 ) in the S 1 (U 1 ) case. We present in Fig. 6a contours of constant δ for the SM + S 1 scenario in the m S 1 -y L dν e plane. Since the S 1 contribution simply adds to the SM one, a small mass and a large value for the LQ-neutrino-quark coupling will lead to an enhanced number of events in each bin. This is beneficial for (detrimental to) the bins where there is an observed excess (lack) of events compared to the SM case.
The hatched region above the blue line in Fig. 6a is currently excluded at the 2 σ level by the APV results. The region to the right of the black dotted curve in Fig. 6a results in a fit to the IceCube data that is marginally better than the SM one.
The contours of constant value of δ in the SM + U 1 scenario are shown in Fig. 6b for the m U 1 -χ L de parameter space. The region above the black dotted curve results in a fit of the IceCube data that is better than the SM one. It can be seen, through comparison of Fig. 6a, b, that the U 1 scenario has much larger region of parameter space that results in a better fit compared to the SM case than the S 1 scenario due to the fact that U 1 signatures interfere with the SM. The 2 σ limit on χ L de from the APV experiment, i.e., χ L de ≤ 0.34 × m U 1 /(1 TeV), as a function of m U 1 is shown by the blue line, with the region above the blue line being excluded.
The couplings of S 1 and U 1 with the right-handed neutrinos are first fixed to zero, for simplicity. We find that the best fit to the recent IceCube data in the U 1 case is obtained for m U 1 = 710 GeV and χ L de = 1.25 and results in a 9.5% improvement over the SM fit. This is in contrast to the S 1 case, which for most of the parameter space considered in The m S1 -y L dνe (left panel) and m U1 -χ L de (right panel) parameter space where the region above the blue line is disallowed at 2 σ level from 2635 days of IceCube data. We also show the low energy disallowed region (space above purple line) at the 2 σ level. The region to the left of the brown line is currently disallowed at 95% C.L. by the 13 TeV LHC data on the LQ direct searches via the j jνν, j jee, and j jμμ final states. The single LQ production is also included in the S 1 case our work either results in a fit worse or comparable to the SM. Therefore we show in the left (right) plot of Fig. 7 the contribution of U 1 (S 1 ) for mass of 710 GeV and χ L de = 1.25 (y L dν e = 1.25), which gives δ = 9.5% (δ = −2.9%). All other couplings are set to zero.
We see from the above analysis that no significant deviation from the SM prediction is seen in the current IceCube data. We use this information to put an upper bound on the y L dν e (χ L de ) coupling as a function of m S 1 (m U 1 ) through a binned likelihood analysis with the Poisson likelihood function [54]. This constraint obtained on the S 1 and U 1 parameter space is then compared with the results from the low-energy flavor experiments and the LHC in the next section.

Combined analysis of the low-energy flavor physics, LHC, and IceCube constraints
Our goal is to combine the low-energy flavor physics, LHC, and IceCube constraints on the parameter spaces associated with the S 1 and U 1 scenarios. The summary of our analysis of these constraints on the m S 1 (m U 1 )-y L dν e (χ L de ) parameter space is shown in Fig. 8. The LHC constraints from the LQ pair production with the dijet + MET, j jee, and j jμμ final states are considered and the currently allowed space, at 95% C.L., is shown by the area below (right of) the brown line in case of S 1 (U 1 ). The region below 850 GeV (1.6 TeV) in the S 1 (U 1 ) case is completely excluded by the LHC data. The region below the purple line in Fig. 8 is allowed by  (a) (b) Fig. 9 The allowed parameter space at 95% C.L. from the 13 TeV LHC direct LQ searches is shown in brown in the y L dνe -y R dν plane, with y L dνμ = y L dντ = 0 for m S1 = 800 GeV (left panel) and m S1 = 1 TeV the flavor observables and the radiative decays of leptons. The region above the blue line is currently disallowed by IceCube data at 2 σ level. For our statistical analysis of the IceCube constraints, the LQ mass and the couplings are kept as free parameters, with γ and C 0 fixed to the IceCube best fit data. We have used, for our numerical calculation, the (1/3 : 1/3 : 1/3) flavor ratio for the incoming flux. Even though the IceCube data alone favors U 1 over S 1 the actual parameter space allowed by the low-energy flavor experiments and LHC corresponds to the region where S 1 (U 1 ) marginally improves (spoils) the IceCube data fit when compared to the SM case.
Overall we find that the limits obtained from the most recent IceCube data are considerably weaker when compared to the constraints from the low-energy observables and direct LQ searches at the LHC. This is mostly due to the current lack of statistics in the high-energy bins of the IceCube spectrum. We next discuss the effects of inclusion of y R dν on our analysis. The leading processes for the LQ pair production at the LHC will be via the initial state of gg, uū, and dd. The cross section in the S 1 case will be particularly enhanced through the dd initial state for large values of the right-handed y R dν couplings. That regime will generate large branching ratio of S 1 to jν and will, therefore, be strongly constrained by the j jνν final state searches at the LHC. We show in Fig. 9 the allowed parameter space from the LHC in brown in the y L dν e -y R dν plane with y L dν μ = y L dν τ = 0 for m S 1 = 800 GeV (left panel) and 1 TeV (right panel) and contours of constant values of δ in green. Since the IceCube data are more sensitive to the ν e coupling we consider this particular choice to check the effect of y R dν on y L dν e . Clearly, the inclusion of y R dν slightly spoils the fit of the IceCube data in the region of interest. This is expected since the coupling y R dν only appears quadratically in the final state of the ν N cross section. Note that y R dν is also tightly constrained by the flavor observables and we present in Fig. 9 the allowed region by purple points.
The inclusion of the U 1 couplings to the right-handed neutrinos follows the same pattern as in the S 1 case. Provided that the right-handed neutrinos do not contribute to the initial state, the current IceCube data are not sensitive to these couplings. The consideration of the righthanded neutrinos in the initial state will lead to a change in the initial flux at the source. This flux will depend on the mass of ν R and also on the possible decay channels in case of heavy ν R . This is beyond the scope of present analysis. The resulting final high-energy cosmic neutrino flux ratios on earth by the possible mixing between the three active neutrinos and the fourth sterile neutrino have been studied in Ref. [55]. The explanation of the PeV neutrinos at IceCube, with the consideration of heavy righthanded neutrino, acting as a dark matter has been studied in Ref. [56].

Conclusions
We analyse the impact of the latest low-energy flavor physics measurements, LHC search limits, and IceCube data on the parameter space of the electroweak SU (2) singlet scalar (vector) LQ S 1 (U 1 ). We perform, in particular, a thorough analysis of the viability of the S 1 (U 1 ) scenario assuming non-zero couplings between S 1 (U 1 ), down quark, and neu-trinos (charged leptons) of all three generations, where the LQ in question couples to the SM lepton doublets. Consequentially, the SU (2) symmetry of the SM requires that S 1 (U 1 ) couples up-type quarks to charged leptons (neutrinos).
We find that the limits obtained from the most recent Ice-Cube data are considerably weaker when compared to the constraints from the low-energy observables and direct LQ searches at the LHC and we quantify this inferiority. We attribute this disparity in constraining power to the current lack of statistics in the high-energy bins of the IceCube spectrum. Even though the IceCube data alone favors U 1 over S 1 , the viable parameter space allowed by the low-energy flavor physics and the LHC data analyses singles out the region where S 1 (U 1 ) marginally improves (spoils) the IceCube data fit when compared to the SM case. In this region, in the S 1 (U 1 ) case, m S 1 ≥ 900 GeV (m U 1 ≥ 1.6 TeV) and the relevant S 1 -neutrino-down quark (U 1 -neutrino-up quark) coupling is small. We have also verified that the couplings of both S 1 and U 1 to the right-handed neutrinos are not being sensitive to the current IceCube data provided that the right-handed neutrinos only contribute to the final state. Since we investigate scenarios where LQs primarily couple to the first generation quarks, the most important constraints originate from flavor physics measurements, followed by the LHC search limits that take over flavor physics limits in the large LQ mass regime.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The experimental data used in the present study has been published by the IceCube collaboration [51]. The LHC results used throughout the manuscript are published by the ATLAS [39] and CMS [46,47] collaborations.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

A Formulas
We list here the different branching ratios used for our analysis in Sect. 3.
The π → ν and the τ − → π − ν branching ratios at the leading order in SM are given by: The electroweak corrections to BR(π → ν) were calculated in Ref. [31] and for BR(τ − → π − ν) in Ref. [57]. The relevant branching ratios in the LQ models is given by,

B Neutrino-nucleon differential cross sections
The neutrino-nucleon scattering in the SM gives rise to the charged current (ν N → X ) and the neutral current (ν N → ν X ) interactions mediated by W and Z bosons, respectively. The target nucleon N is an isoscalar nucleon with N = (n + p)/2, X is the hadronic final state, and = e, μ, τ . The SM differential cross sections in terms of the scaling variables are given as E denotes the energy carried away by the outgoing lepton or the neutrino in the laboratory frame while x is the fraction of the initial nucleon momentum taken by the struck quark.
Here, q(x, Q 2 ) andq(x, Q 2 ) (q 0 (x, Q 2 ) andq 0 (x, Q 2 )) are the quark and anti-quark density distributions in a proton, respectively, summed over valence and sea quarks of all flavors relevant for CC (NC) interactions: with the chiral couplings given by The parton distribution functions (PDFs) of the quarks are evaluated at energy Q 2 , and the Mathematica package MSTW [58] is used throughout this work. There are also neutrino-electron interactions, but they can be generally neglected with respect to the neutrino-nucleon cross section because of the smallness of electron's mass, except for the resonant formation of the intermediate W − boson in theν e e interactions at around E ν = m 2 W /(2m e ) = 6.3 × 10 6 GeV, known as the Glashow resonance. The differential cross sections for all the neutrino electron reactions are listed in Ref. [59].
B.1 ν j N cross sections in the presence of the S 1 The differential ν j N cross sections in the presence of the S 1 interaction are given by where j = 1,2,3 with ch = CC, N C. The coefficients are given by N CC = 3 k=1 |(V * y L ) 1k | 2 , N NC = 3 k=1 |(y L U ) 1k | 2 + 3 k=1 |y R 1k | 2 . The decay width S 1 of S 1 given by Note that the effect of the right-handed couplings to the neutrinos is only visible in NC interactions.
B.2 ν j N cross sections in the presence of the U 1 The modified q 0 (x, Q 2 ) andq 0 (x, Q 2 ) listed in Eqs. (B.6) and (B.7), in the presence of U 1 , are given below.
where j = 1, 2, 3. Eqs. (B.10) and (B.11) correspond to the case when the same flavor neutrino is in the initial and the final states and there is an interference with the SM contribution. There will be additional contributions from the cases where the final state will consist of the right-handed neutri-nos or will have a neutrino of different flavor from the initial one and are given by In case of CC interactions, the coefficient m 2 taking into account the interference terms only when the final state is similar to the SM. The decay width of U 1 is given by