Two unusual neutrino events in the Antarctic Impulse Transient Antenna (ANITA) appear to have been generated by air showers from a particle emerging from the Earth at angles 25-35 degrees above the horizon. We evaluate the effective aperture for ANITA with a simplified detection model to illustrate the features of the angular dependence of expected events for incident standard model tau neutrinos and for sterile neutrinos that mix with tau neutrinos. We apply our sterile neutrino aperture results to a dark matter scenario with long-lived supermassive dark matter that decay to sterile neutrino-like particles. We find that for up-going air showers from tau decays, from isotropic fluxes of standard model, sterile neutrinos or other particles that couple to the tau through suppressed weak interaction cross sections cannot be responsible for the unusual events.

In view of the IceCube's 6-year high-energy starting events (HESE) sample, we revisit the possibility that the updated data may be better explained by a combination of neutrino fluxes from dark matter decay and an isotropic astrophysical power-law than purely by the latter. We find that the combined two-component flux qualitatively improves the fit to the observed data over a purely astrophysical one, and discuss how these updated fits compare against a similar analysis done with the 4-year HESE data. We also update fits involving dark matter decay via multiple channels, without any contribution from the astrophysical flux. We find that a DM-only explanation is not excluded by neutrino data alone. Finally, we also consider the possibility of a signal from dark matter annihilations and perform analogous analyses to the case of decays, commenting on its implications.

Sterile neutrinos with mass in the eV-scale and large mixings of order $\theta_0\simeq 0.1$ could explain some anomalies found in short-baseline neutrino oscillation data. Here, we revisit a neutrino portal scenario in which eV-scale sterile neutrinos have self-interactions via a new gauge vector boson $\phi$. Their production in the early Universe via mixing with active neutrinos can be suppressed by the induced effective potential in the sterile sector. We study how different cosmological observations can constrain this model, in terms of the mass of the new gauge boson, $M_\phi$, and its coupling to sterile neutrinos, $g_s$. Then, we explore how to probe part of the allowed parameter space of this particular model with future observations of the diffuse supernova neutrino background by the Hyper-Kamiokande and DUNE detectors. For $M_\phi \sim 5-10$~keV and $g_s \sim 10^{-4}-10^{-2}$, as allowed by cosmological constraints, we find that interactions of diffuse supernova neutrinos with relic sterile neutrinos on their way to the Earth would result in significant dips in the neutrino spectrum which would produce unique features in the event spectra observed in these detectors.

After the first four years of data taking, the IceCube neutrino telescope has observed 54 high-energy starting events (HESE) with deposited energies between 20 TeV and 2 PeV. The background from atmospheric muons and neutrinos is expected to be of about 20 events, all below 100 TeV, thus pointing towards the astrophysical origin of about 8 events per year in that data set. However, their precise origin remains unknown. Here, we perform a detailed analysis of this event sample (considering simultaneously the energy, hemisphere and topology of the events) by assuming two contributions for the signal events: an isotropic power-law flux and a flux from decaying heavy dark matter. We fit the mass and lifetime of the dark matter and the normalization and spectral index of an isotropic power-law flux, for various decay channels of dark matter. We find that a significant contribution from dark matter decay is always slightly favored, either to explain the excess below 100 TeV, as in the case of decays to quarks or, as in the case of neutrino channels, to explain the three multi-PeV events. Also, we consider the possibility to interpret all the data by dark matter decays only, considering various combinations of two decay channels. We show that the decaying dark matter scenario provides a better fit to HESE data than the isotropic power-law flux.

We consider propagation of high energy earth-skimming taus produced in interactions of astrophysical tau neutrinos. For astrophysical tau neutrinos we take generic power-law flux, $E^{-2}$ and the cosmogenic flux initiated by the protons. We calculate tau energy loss in several approaches, such as dipole models and the phenomenological approach in which parameterization of the $F_2$ is used. We evaluate the tau neutrino charged-current cross section using the same approaches for consistency. We find that uncertainty in the neutrino cross section and in the tau energy loss partially compensate giving very small theoretical uncertainty in the emerging tau flux for distances ranging from $2$ km to $100$ km and for the energy range between $10^6$ GeV and $10^{11}$ GeV, focusing on energies above $10^8$ GeV. When we consider uncertainties in the neutrino cross section, inelasticity in neutrino interactions and the tau energy loss, which are not correlated, i.e. they are not all calculated in the same approach, theoretical uncertainty ranges from about $30\%$ and $60 \%$ at $10^8$ GeV to about factors of 3.3 and 3.8 at $10^{11}$ GeV for the $E^{-2}$ flux and the cosmogenic flux, respectively, for the distance of 10 km rock. The spread in predictions significantly increases for much larger distances, e.g., $\sim 1,000$ km. Most of the uncertainty comes from the treatment of photonuclear interactions of the tau in transit through large distances. We also consider Monte Carlo calculation of the tau propagation and we find that the result for the emerging tau flux is in agreement with the result obtained using analytic approach. Our results are relevant to several experiments that are looking for skimming astrophysical taus, such as the Pierre Auger Observatory, HAWC and Ashra. We evaluate the aperture for the Auger and discuss briefly application to the the other two experiments.

We evaluate the prompt atmospheric neutrino flux using the different QCD models for heavy quark production including the $b$ quark contribution. We include the nuclear correction and find it reduces the fluxes by $10 \% - 50\%$ according to the models. Our heavy quark results are compared with experimental data from RHIC, LHC and LHCb.

We evaluate the prompt atmospheric neutrino flux including nuclear correction and $B$ hadron contribution in the different frameworks: NLO perturbative QCD and dipole models. The nuclear effect is larger in the prompt neutrino flux than in the total charm production cross section, and it reduces the fluxes by $10\% - 30\%$ depending on the model. We also investigate the uncertainty using the QCD scales allowed by the charm cross section data from RHIC and LHC experiments.

We evaluate the prompt atmospheric neutrino flux at high energies using three different frameworks for calculating the heavy quark production cross section in QCD: NLO perturbative QCD, $k_T$ factorization including low-$x$ resummation, and the dipole model including parton saturation. We use QCD parameters, the value for the charm quark mass and the range for the factorization and renormalization scales that provide the best description of the total charm cross section measured at fixed target experiments, at RHIC and at LHC. Using these parameters we calculate differential cross sections for charm and bottom production and compare with the latest data on forward charm meson production from LHCb at $7$ TeV and at $13$ TeV, finding good agreement with the data. In addition, we investigate the role of nuclear shadowing by including nuclear parton distribution functions (PDF) for the target air nucleus using two different nuclear PDF schemes. Depending on the scheme used, we find the reduction of the flux due to nuclear effects varies from $10\%$ to $50 \%$ at the highest energies. Finally, we compare our results with the IceCube limit on the prompt neutrino flux, which is already providing valuable information about some of the QCD models.

We re-evaluate the prompt atmospheric neutrino flux, using the measured charm cross sections at RHIC and the Large Hadron Collider to constrain perturbative QCD parameters such as the factorization and renormalization scales, as well as modern parton distribution functions and recent estimates of the cosmic-ray spectra. We find that our result for the prompt neutrino flux is lower than previous perturbative QCD estimates and, consequently, alters the signal-to-background statistics of the recent IceCube measurements at high energies.

We investigate whether the recent ultra-high energy (UHE) neutrino events detected at the IceCube neutrino observatory could come from the decay of charmed mesons produced within the mildly relativistic jets of supernova-like astrophysical sources. We demonstrate that the $5.7\sigma$ excess of neutrinos observed by IceCube in the energy range 30 TeV--2 PeV can be explained by a diffuse flux of neutrinos produced in such slow-jet supernovae, using the values of astrophysical and QCD parameters within the theoretical uncertainties associated with neutrino production from charmed meson decay in astrophysical sources. We discuss the theoretical uncertainties inherent in the evaluation of charm production in high energy hadronic collisions, as well as the astrophysical uncertainties associated with slow-jet supernova sources. The proton flux within the source, and therefore also the produced neutrino flux, is cut off at around a few PeV, when proton cooling processes become dominant over proton acceleration. This directly explains the sudden drop in event rates at energies above a few PeV. We incorporate the effect of energy dependence in the spectrum-weighted charm production cross-section and show that this has a very significant effect on the shape, magnitude and cut-off energies for the neutrino flux.

We calculate the number of dark matter particles that a neutron star accumulates over its lifetime as it rotates around the center of a galaxy, when the dark matter particle is a self-interacting boson but does not self-annihilate. We take into account dark matter interactions with baryonic matter and the time evolution of the dark matter sphere as it collapses within the neutron star. We show that dark matter self-interactions play an important role in the rapid accumulation of dark matter in the core of the neutron star. We consider the possibility of determining an exclusion region of the parameter space for dark matter mass and dark matter interaction cross sections based on the observation of old neutron stars with strong dark matter self-interactions. We show that for a dark matter density of $~10^3$ GeV/cm$^3$ and dark matter mass $m_\chi$ less than approximately 10 GeV, there is a potential exclusion region for dark matter interactions with nucleons that is three orders of magnitude more stringent than without self-interactions. The potential exclusion region for dark matter self-interaction cross sections is many orders of magnitude stronger than the current Bullet Cluster limit. For example, for high dark matter density regions, we find that for $m_\chi\sim 10$ GeV when the dark matter interaction cross section with the nucleons ranges from $\sigma_{\chi n}\sim 10^{-52}$ cm$^2$ to $10^{-57}$ cm$^2$, the dark matter self-interaction cross section limit is $\sigma_{\chi\chi}< 10^{-33}$ cm$^2$, which is about ten orders of magnitude stronger than the Bullet Cluster limit.

Neutrino production of radio Cherenkov signals in the Moon is the object of radio telescope observations. Depending on the energy range and detection parameters, the dominant contribution to the neutrino signal may come from interactions of the neutrino on the Moon facing the telescope, rather than neutrinos that have traversed a portion of the Moon. Using the approximate analytic expression of the effective lunar aperture from a recent paper by Gayley, Mutel and Jaeger, we evaluate the background from cosmic ray interactions in the lunar regolith. We also consider the modifications to the effective lunar aperture from generic non-standard model neutrino interactions. A background to neutrino signals are radio Cherenkov signals from cosmic ray interactions. For cosmogenic neutrino fluxes, neutrino signals will be difficult to observe because of low neutrino flux at the high energy end and large cosmic ray background in the lower energy range considered here. We show that lunar radio detection of neutrino interactions is best suited to constrain or measure neutrinos from astrophysical sources and probe non-standard neutrino-nucleon interactions such as microscopic black hole production.

We calculate the contained and upward muon and shower fluxes due to neutrinos produced via dark matter annihilation or decay in the Galactic center. We consider dark matter models in which the dark matter particle is a gravitino, a Kaluza-Klein particle and a particle in leptophilic models. The Navarro-Frenk-White profile for the dark matter density distribution in the Galaxy is used. We incorporate neutrino oscillations by assuming maximal mixing and parametrize our results for muon and shower distributions. The muon and shower event rates and the minimum observation times in order to reach 2$\sigma$ detection significance are evaluated. We illustrate how observation times vary with the cone half angle chosen about the Galactic center, with the result that the optimum angles are about 10$^\circ$ and 50$^\circ$ for the muon events and shower events, respectively. We find that for the annihilating dark matter models such as the leptophilic and Kaluza-Klein models, upward and contained muon as well as showers are promising signals for dark matter detection in just a few years of observation, whereas for decaying dark matter models, the same observation times can only be reached with showers. We also illustrate for each model the parameter space probed with the 2$\sigma$ signal detection in five years. We discuss how the shape of the parameter space probed change with significance and the observation time.

We calculate contained and upward muon flux and contained shower event rates from neutrino interactions, when neutrinos are produced from annihilation of the dark matter in the Galactic Center. We consider model-independent direct neutrino production and secondary neutrino production from the decay of taus, W bosons and bottom quarks produced in the annihilation of dark matter. We illustrate how muon flux from dark matter annihilation has a very different shape than the muon flux from atmospheric neutrinos. We also discuss the dependence of the muon fluxes on the dark matter density profile and on the dark matter mass and of the total muon rates on the detector threshold. We consider both the upward muon flux, when muons are created in the rock below the detector, and the contained flux when muons are created in the (ice) detector. We also calculate the event rates for showers from neutrino interactions in the detector and show that the signal dominates over the background for $150 {\rm GeV} <m_\chi < 1$ TeV for $E_{sh}^{th} = 100$ GeV.

Charm production gives rise to a flux of very high energy neutrinos from astrophysical sources with jets driven by central engines, such as gamma ray bursts or supernovae with jets. The neutrino flux from semi-leptonic decays of charmed mesons is subject to much less hadronic and radiative cooling than the conventional flux from pion and kaon decays and therefore has a dominant contribution at higher energies, of relevance to future ultrahigh energy neutrino experiments.

We calculate the prompt neutrino flux from atmospheric charm production by cosmic rays, using the dipole picture in a perturbative QCD framework, which incorporates the parton saturation effects present at high energies. We compare our results with the next-to-leading order perturbative QCD result and find that saturation effects are large for neutrino energies above 10^6 GeV, leading to a substantial suppression of the prompt neutrino flux. We comment on the range of prompt neutrino fluxes due to theoretical uncertainties.

We study the signal for the detection of quasi-stable supersymmetric particles produced in interactions of cosmogenic neutrinos. We consider energy loss of high energy staus due to photonuclear and weak interactions. We show that there are optimal nadir angles for which the stau signal is a factor of several hundred larager than muons. We discuss how one could potentially eliminate the muon background by considering the energy loss of muons in the detector. We also show results for the showers produced by weak interactions of staus that reach the detector.

Models of late-time neutrino mass generation contain new interactions of the cosmic background neutrinos with supernova relic neutrinos (SRNs) through exchange of the on-shell light boson, leading to significant modification of the differential SRN flux observed at earth. We consider Abelian U(1) model for generating neutrino masses at low scales and we show that there is a large parameter space in this model for which the changes induced in the flux by the exchange of the light bosons might allow one to distinguish between neutrinos being Majorana or Dirac particles, the type of neutrino mass hierarchy (normal or inverted or quasi-degenerate), and could also possibly determine the absolute values of the neutrino masses. Measurements of the presence of these effects would be possible at the next-generation water Cerenkov detectors enriched with Gadolinium, or a large 100 kton liquid argon detector.

In models in which neutrinos are light, due to a low scale of symmetry breaking, additional light bosons are generically present. We show that the interaction between diffuse relic supernova neutrinos (RSN) and the cosmic background neutrinos, via exchange of these light scalars, can result in a dramatic change of the supernova (SN) neutrinos flux. Measurement of this effect with current or future experiments can provide a spectacular direct evidence for the low scale models. We demonstrate how the observation of neutrinos from SN1987A constrains the symmetry breaking scale of the above models. We also discuss how current and future experiments may confirm or further constrain the above models, either by detecting the ``accumulative resonance'' that diffuse RSN go through or via a large suppression of the flux of neutrinos from nearby < O(Mpc) SN bursts.

The potential for neutrino telescopes to discover charged stau production in neutrino-nucleon interactions in Earth depends in part on the stau lifetime and range. In some supersymmetric scenarios, the next lightest supersymmetric particle is a stau with a decay length on the scale of 10 km. We evaluate the electromagnetic energy loss as a function of energy and stau mass. The energy loss parameter $\beta$ scales as the inverse stau mass for the dominating electromagnetic processes, photonuclear and $e^+e^-$ pair production. The range can be parameterized as a function of stau mass, initial energy and minimum final energy. In comparison to earlier estimates of the stau range, our results are as much as a factor of two larger, improving the potential for stau discovery in neutrino telescopes.

Astrophysical sources of ultrahigh energy neutrinos yield tau neutrino fluxes due to neutrino oscillations. We study in detail the contribution of tau neutrinos with energies above PeV relative to the contribution of the other flavors. We consider several different initial neutrino fluxes and include tau neutrino regeneration in transit through the Earth and energy loss of charged leptons. We discuss signals of tau neutrinos in detectors such as IceCube, RICE and ANITA.

In scenarios with large extra dimensions and TeV scale gravity ultrahigh energy neutrinos produce black holes in their interactions with the nucleons. We show that ICECUBE and OWL may observe large number of black hole events and provide valuable information about the fundamental Planck scale and the number of extra dimensions. OWL is especially well suited to observe black hole events produced by neutrinos from the interactions of cosmic rays with the 3 K background radiation. Depending on the parameters of the scenario of large extra dimensions and on the flux model, as many as 28 events per year are expected for a Planck scale of 3 TeV.

We evaluate the upward shower and muon event rates for two characteristic four neutrino mixing models for extragalactic neutrinos, as well as for the atmospheric neutrinos, with energy thresholds of 1 TeV, 10 TeV and 100 TeV. We show that by comparing the shower to muon event rates, one can distinguish between oscillation and no-oscillation models. By measuring shower and muon event rates for energy thresholds of 10 TeV and 100 TeV, and by considering their ratio, it is possible to use extragalactic neutrino sources to determine the type of four-flavor mixing pattern. We find that one to ten years of data taking with kilometer-size detector has a very good chance of providing valuable information about the physics beyond the Standard Model.

The photonuclear contribution to charged lepton energy loss has been re-evaluated taking into account HERA results on real and virtual photon interactions with nucleons. With large $Q^2$ processes incorporated, the average muon range in rock for muon energies of $10^9$ GeV is reduced by only 5% as compared with the standard treatment. We have calculated the tau energy loss for energies up to $10^9$ GeV taking into consideration the decay of the tau. A Monte Carlo evaluation of tau survival probability and range show that at energies below $10^7-10^8$ GeV, depending on the material, only tau decays are important. At higher energies the tau energy losses are significant, reducing the survival probability of the tau. We show that the average range for tau is shorter than its decay length and reduce to 17 km in water for an incident tau energy of $10^9$ GeV, as compared with its decay length of 49 km at that energy. In iron, the average tau range is 4.7 km for the same incident energy.

The appearance of high energy tau neutrinos due to $\nu_\mu \to \nu_\tau$ oscillations of extragalactic neutrinos can be observed by measuring the neutrino induced upward hadronic and electromagnetic showers and upward muons. We evaluate quantitatively the tau neutrino regeneration in the Earth for a variety of extragalactic neutrino fluxes. Charged-current interactions of the upward tau neutrinos below and in the detector, and the subsequent tau decay create muons or hadronic and electromagnetic showers. The background for these events are muon neutrino and electron neutrino charged-current and neutral-current interactions, where in addition to extragalactic neutrinos, we consider atmospheric neutrinos. We find significant signal to background ratios for the hadronic/electromagnetic showers with energies above 10 TeV to 100 TeV initiated by the extragalactic neutrinos. We show that the tau neutrinos from point sources also have the potential for discovery above a 1 TeV threshold. A kilometer-size neutrino telescope has a very good chance of detecting the appearance of tau neutrinos when both muon and hadronic/electromagnetic showers are detected.

We propose a novel approach for studying $\nu_\mu \to \nu_\tau$ oscillations with extragalactic neutrinos. Active Galactic Nuclei and Gamma Ray Bursts are believed to be sources of ultrahigh energy muon neutrinos. With distances of 100 Mpc or more, they provide an unusually long baseline for possible detection of $\nu_\mu \to \nu_\tau$ with mixing parameters $\Delta m^2$ down to $10^{-17}$eV$^2$, many orders of magnitude below the current accelerator experiments. By solving the coupled transport equations, we show that high-energy $\nu_\tau$'s, as they propagate through the earth, cascade down in energy, producing the enhancement of the incoming $\nu_\tau$ flux in the low energy region, in contrast to the high-energy $\nu_\mu$'s, which get absorbed. For an AGN quasar model we find the $\nu_\tau$ flux to be a factor of 2 to 2.5 larger than the incoming flux in the energy range between $10^2$ GeV and $10^4$ GeV, while for a GRB fireball model, the enhancement is 10%-27% in the same energy range and for zero nadir angle. This enhancement decreases with larger nadir angle, thus providing a novel way to search for $\nu_\tau$ appearance by measuring the angular dependence of the muons. To illustrate how the cascade effect and the $\nu_\tau$ final flux depend on the steepness of the incoming $\nu_\tau$, we show the energy and angular distributions for several generic cases of the incoming tau neutrino flux, $F_\nu^0 \sim E^{-n}$ for n=1,2 and 3.6. We show that for the incoming flux that is not too steep, the signal for the appearance of high-energy $\nu_\tau$ is the enhanced production of lower energy $\mu$ and their distinctive angular dependence, due to the contribution from the $\tau$ decay into $\mu$ just below the detector.

Using next-to-leading order perturbative QCD, we find that the charm contribution to the atmospheric lepton fluxes dominates over the conventional ones from pi and K decays for energies higher than 10^5 GeV. We also discuss theoretical uncertainties involved in the calculations and compare our results with previous evaluations.

We report new calculations of the cross sections for deeply inelastic neutrino-nucleon scattering at neutrino energies between $10^{9}\ev$ and $10^{21}\ev$. We compare with results in the literature and assess the reliability of our predictions. For completeness, we briefly review the cross sections for neutrino interactions with atomic electrons, emphasizing the role of the $W$-boson resonance in $\bar{\nu}_{e}e$ interactions for neutrino energies in the neighborhood of $6.3\pev$. Adopting model predictions for extraterrestrial neutrino fluxes from active galactic nuclei, gamma-ray bursters, and the collapse of topological defects, we estimate event rates in large-volume water \v{C}erenkov detectors and large-area ground arrays.

The charm contribution to the atmospheric fluxes of muons and muon neutrinos may be enhanced by as much as a factor of 10 when one includes the contributions of D->pi,K->leptons and folds in uncertainties in the charm cross section and energy distribution. In the energy range considered here, from 100 GeV to 10 TeV, the charm contribution is small compared to the conventional flux of muons and muon neutrinos.

We present a calculation of the fluxes of muons and muon neutrinos from the decays of pions and kaons that are themselves the decay products of charmed particles produced in the atmosphere by cosmic ray-air collisions. Using the perturbative cross section for charm production, these lepton fluxes are two to three orders of magnitude smaller than the fluxes from the decays of pions and kaons directly produced in cosmic ray-air collisions. Intrinsic charm models do not significantly alter our conclusions, nor do models with a charm cross section enhanced in the region above an incident cosmic ray energy of 1 TeV.

The recent detection of a gamma-ray flux from the direction of the Galactic center by EGRET on the Compton GRO raises the question of whether this is a point source (possibly coincident with the massive black hole candidate Sgr A*) or a diffuse emitter. Using the latest experimental particle physics data and theoretical models, we examine in detail the gamma-ray spectrum produced by synchrotron, inverse Compton scattering and mesonic decay resulting from the interaction of relativistic protons with hydrogen accreting onto a point-like object. Such a distribution of high-energy baryons may be expected to form within an accretion shock as the inflowing gas becomes supersonic. This scenario is motivated by hydrodynamic studies of Bondi-Hoyle accretion onto Sgr A*, which indicate that many of its radiative characteristics may ultimately be associated with energy liberated as this plasma descends down into the deep potential well. Earlier attempts at analyzing this process concluded that the EGRET data are inconsistent with a massive point-like object. Here, we demonstrate that a more careful treatment of the physics of p-p scattering suggests that a ~10^6 solar mass black hole may be contributing to this high-energy emission.

We present results for neutrino-nucleon cross sections for energies up to 10$^21$eV, of relevance to the detection of ultrahigh energy galactic and extragalactic neutrinos. At the highest energies, our results are about 2.4 times larger than previous estimates. Using these new cross sections, we predict neutrino telescope event rates for the upward moving muons initiated by the neutrino interactions in the Earth and for contained-vertex events in the PeV range due to neutrino-electron interactions. We show that future neutrino detectors, such as AMANDA, BAIKAL, DUMAND and NESTOR have a very good chance of detecting neutrinos which originate in the Active Galactic Nuclei.

Cross sections for the interactions of ultrahigh-energy neutrinos with nucleons are evaluated in light of new information about nucleon structure functions. For $10^{20}$--eV neutrinos, the cross section is about 2.4 times previous estimates. We also review the cross sections for neutrino interactions with atomic electrons. Some consequences for interaction rates in the Earth and for event rates from generic astrophysical sources in large-scale detectors are noted.

In light of recent measurements of the nucleon structure function in the small-$x$ deep-inelastic regime at HERA and the consequently improved theoretical understanding of the quark distributions in this range of parton fractional momentum, we present new results for the neutrino-nucleon cross section at ultra high energies, up to $ E_{\nu} \sim 10^{21}$ eV. The results are relevant to deep underground muon detectors looking for such neutrinos. For $\simeq 10^{20}$ eV neutrinos, our cross sections are about a factor of $4$ to $10$ times larger than the previously reported results by Reno, Quigg and Walker. We discuss implications of this new neutrino-nucleon cross section for a variety of current and future neutrino detectors.