The Antarctic Impulsive Transient Antenna (ANITA) experiment has observed two air shower events with energy $\sim 500~{\rm PeV}$ emerging from the Earth with exit angles $\sim 30^\circ$ above the horizon. As was immediately noted by the ANITA Collaboration, these events (in principle) could originate in the atmospheric decay of an upgoing $\tau$-lepton produced through a charged current interaction of a $\nu_\tau$ inside the Earth. However, the relatively steep arrival angles of these perplexing events create tension with the standard model (SM) neutrino-nucleon interaction cross section. Deepening the conundrum, the IceCube neutrino telescope and the Pierre Auger Observatory with substantially larger exposures to cosmic $\nu_\tau$'s in this energy range have not observed any events. This lack of observation implies that the messenger particle (MP) giving rise to ANITA events must produce an air shower event rate at least a factor of 40 larger than that produced by a flux of $\tau$-neutrinos to avoid conflicts with the upper limits reported by the IceCube and the Pierre Auger collaborations. In addition, the sensitivity of ANITA to MP-induced events must be comparable to or larger than those of IceCube and Auger to avoid conflict with the non-observation of any signal at these facilities. Beyond SM interpretations of ANITA events can be classified according to whether the MPs: (i) live inside the Earth, (ii) originate in neutrino-nucleon collisions inside the Earth, (iii) come from cosmological distances. In this communication we investigate the positive and negative facets of these three classes of models.

We discuss the spectrum of the different components in the astrophysical neutrino flux reaching the Earth and the possible contribution of each component to the high-energy IceCube data. We show that the diffuse flux from cosmic ray interactions with gas in our galaxy implies just 2 events among the 54 event sample. We argue that the neutrino flux from cosmic ray interactions in the intergalactic (intracluster) space depends critically on the transport parameter $\delta$ describing the energy dependence in the diffusion coefficient of galactic cosmic rays. Our analysis motivates a E^{-2.1} neutrino spectrum with a drop at PeV energies that fits well the data, including the non-observation of the Glashow resonance at 6.3 PeV. We also show that a cosmic ray flux described by an unbroken power law may produce a neutrino flux with interesting spectral features (bumps and breaks) related to changes in the cosmic ray composition.

We find that the ratio r_{mu e} of the muon to the electromagnetic component of an extended air shower at the ground level provides an indirect measure of the depth X_{max} of the shower maximum. This result, obtained with the air-shower code AIRES, is independent of the hadronic model used in the simulation. We show that the value of r_{mu e} in a particular shower discriminates its proton or iron nature with a 98% efficiency. We also show that the eventual production of forward heavy quarks inside the shower may introduce anomalous values of r_{mu e} in isolated events.

IceCube has recently observed 37 events of TeV-PeV energies. The angular distribution, with a strong preference for downgoing directions, the spectrum, and the small muon to shower ratio in the data can not be accommodated assuming standard interactions of atmospheric neutrinos. We obtain an excellent fit, however, if a diffuse flux of ultrahigh energy (cosmogenic) neutrinos experiences collisions where only a small fraction of the energy is transferred to the target nucleon. We show that consistent models of TeV gravity or other non-Wilsonian completions of the standard model provide cross sections with these precise features. An increased statistics could clearly distinguish our scenario from the one assumed by IceCube (a diffuse flux of astrophysical neutrinos with a E^{-2} spectrum) and establish the need for new physics in the interpretation of the data.

Very energetic charm and bottom hadrons may be produced in the upper atmosphere when a primary cosmic ray or the leading hadron in an extensive air shower collide with a nucleon. At $E\approx 10^8$ GeV their decay length becomes of the order of 10 km, implying that they tend to interact in the air instead of decaying. Since the inelasticity in these collisions is much smaller than the one in proton and pion collisions, there could be rare events where a heavy-hadron component transports a significant amount of energy deep into the atmosphere. We have developed a module for the detailed simulation of these processes and have included it in a new version of the air shower simulator AIRES. We study the frequency, the energy distribution and the depth of charm and bottom production, as well as the depth and the energy distribution of these quarks when they decay. As an illustration, we consider the production and decay of tau leptons (from $D_s$ decays) and the lepton flux at PeV energies from a 30 EeV proton primary. The proper inclusion of charm and bottom hadrons in AIRES opens the possibility to search for air-shower observables that are sensitive to heavy quark effects.

Collisions of cosmic rays in the atmosphere may produce heavy hadrons of very high energy. The decay length of a $B$ meson of energy above $10^7$ GeV is larger than 1 km, implying that such a particle tends to interact in the air before it decays. We show that the fraction of energy deposited in these interactions is much smaller than in proton and pion collisions. We parameterize their elasticity and determine the average number of interactions and the atmospheric depth at the decay point for different initial energies. We find that the profile of a $3\times 10^9$ GeV bottom shower may be very different from the profile of a proton shower of the same energy, defining either a very deep maximum, or two maxima, or other features that cannot be parameterized with a single Gaisser-Hillas function. Finally, we discuss under what conditions a bottom hadron inside the parent air shower may provide observable effects.

Charmed mesons may be produced when a primary cosmic ray or the leading hadron in an air shower collide with an atmospheric nucleon. At energies \ge 10^8 GeV their decay length becomes larger than 10 km, which implies that they tend to interact in the air instead of decaying. We study the collisions of long-lived charmed hadrons in the atmosphere. We show that (\Lambda_c,D)-proton diffractive processes and partonic collisions of any q^2 where the charm quark is an spectator have lower inelasticity than (p,\pi)-proton collisions. In particular, we find that a D meson deposits in each interaction just around 55% of the energy deposited by a pion. On the other hand, collisions involving the valence c quark (its annihilation with a sea cbar quark in the target or c-quark exchange in the t channel) may deposit most of D meson energy, but their frequency is low (below 0.1% of inelastic interactions). As a consequence, very energetic charmed hadrons may keep a significant fraction of their initial energy after several hadronic interactions, reaching much deeper in the atmosphere than pions or protons of similar energy.

The observation of astrophysical neutrinos requires a detailed understanding of the atmospheric neutrino background. Since neutrinos are produced in meson decays together with a charged lepton, important constraints on this background can be obtained from the measurement of the atmospheric muon flux. Muons, however, can also be produced as mu+ mu- pairs by purely electromagnetic processes. We use the Z-moment method to study and compare the contributions to the atmospheric muon and neutrino fluxes from different sources (pi/K decay, charmed and unflavored hadron decay, and photon conversion into a muon pair). We pay special attention to the contribution from unflavored mesons (eta, eta', rho0, omega and phi). These mesons are abundant in air showers, their lifetimes are much shorter than those of charged pions or kaons, and they have decay branching ratios of order 10^-4 into final states containing a muon pair. We show that they may be the dominant source of muons at E_mu >10^3 TeV.

In order to estimate the possibility to observe exotic physics in a neutrino telescope, it is essential to first understand the flux of atmospheric neutrinos, muons and dimuons. We study the production of these leptons by high-energy cosmic rays. We identify three main sources of muons of energy E > 10^6 GeV: the weak decay of charm and bottom mesons and the electromagnetic decay of unflavored mesons. Contrary to the standard assumption, we find that eta mesons, not the prompt decay of charm hadrons, are the dominant source of atmospheric muons at these energies. We show that, as a consequence, the ratio between the neutrino and muon fluxes is significantly reduced. For dimuons, which may be a background for long-lived staus produced near a neutrino telescope, we find that pairs of E ~ 10^7 GeV forming an angle above 10^-6 rad are produced through D (80%) or B (10%) meson decay and through Drell-Yan proceses (10%). The frequency of all these processes has been evaluated using the jet code PYTHIA.

The collision of a high energy cosmic ray with a nucleon in the upper atmosphere could produce long-lived heavy particles. Such particles would be very penetrating, since the energy loss in matter scales as the inverse mass, and could reach a neutrino telescope like IceCube from large zenith angles. Here we study this possibility and focus on the long-lived stau of SUSY models with a gravitino LSP. The signal would be a pair of muon-like parallel tracks separated by 50 meters along the detector. We evaluate the background of muon pairs and show that any events from zenith angles above 80 degrees could be explained by the production of these heavy particles by cosmic rays.

Cosmic rays from outer space enter the atmosphere with energies of up to 10^{11} GeV. The initial particle or a secondary hadron inside the shower may then interact with an air nucleon to produce nonstandard particles. In this article we study the production of new physics by high energy cosmic rays, focusing on the long-lived gluino of split-SUSY models and a WIMP working as dark matter. We first deduce the total flux of hadron events at any depth in the atmosphere, showing that secondary hadrons can not be neglected. Then we use these results to find the flux of gluinos and WIMPs that reach the ground after being produced inside air showers. We also evaluate the probability of producing these exotic particles in a single proton shower of ultrahigh energy. Finally we discuss the possible signal in current and projected experiments. While the tiny flux of WIMPs does not seem to have any phenomenological consequences, we show that the gluinos could modify substantially the profile of a small fraction of extensive air showers. In particular, they could produce a distinct signal observable at AUGER in showers of large zenith angle.

Cosmogenic neutrinos reach the Earth with energies around 10^9 GeV, and their interactions with matter will be measured in upcoming experiments (Auger, IceCube). Models with extra dimensions and the fundamental scale at the TeV could imply signals in these experiments. In particular, the production of microscopic black holes by cosmogenic neutrinos has been extensively studied in the literature. Here we make a complete analysis of gravity-mediated interactions at larger distances, where they can be calculated in the eikonal approximation. In these processes a neutrino of energy E_\nu interacts elastically with a parton inside a nucleon, loses a small fraction y of its energy, and starts a hadronic shower of energy y E_\nu << E_\nu. We analyze the ultraviolet dependence and the relevance of graviton emission in these processes, and show that they are negligible. We also study the energy distribution of cosmogenic events in AMANDA and IceCube and the possibility of multiple-bang events. For any neutrino flux, the observation of an enhanced rate of neutral current events above 100 TeV in neutrino telescopes could be explained by TeV-gravity interactions. The values of the fundamental scale of gravity that IceCube could reach are comparable to those to be explored at the LHC.

The existence of extra dimensions allows the possibility that the fundamental scale of gravity is at the TeV. If that is the case, gravity could dominate the interactions of ultra-high energy cosmic rays. In particular, the production of microscopic black holes by cosmogenic neutrinos has been estimated in a number of papers. We consider here gravity-mediated interactions at larger distances, where they can be calculated in the eikonal approximation. We show that for the expected flux of cosmogenic neutrinos these elastic processes give a stronger signal than black hole production in neutrino telescopes. Taking the bounds on the higher dimensional Planck mass M_D (D=4+n) from current air shower experiments, for n=2 (6) elastic collisions could produce up to 118 (34) events per year at IceCube. On the other hand, the absence of any signal would imply a bound of M_D>~5 TeV.

In scenarios with the fundamental unification scale at the TeV one expects string excitations of the standard model fields at accessible energies. We study the neutrino-nucleon cross section in these models. We show that duality of the scattering amplitude forces the existence of a tower of massive leptoquarks that mediate the process in the s-channel. Using the narrow-width approximation we find a sum rule for the production rate of resonances with different spin at each mass level. We show that these contributions can increase substantially the standard model neutrino-nucleon cross section, although seem insufficient in order to explain the cosmic ray events above the GZK cutoff energy.