High energy cosmic rays reach the surface of the Sun and start showers with thousands of secondary particles. Most of them will be absorbed by the Sun, but a fraction of the neutral ones will escape and reach the Earth. Here we incorporate a new ingredient that is essential to understand the flux of these solar particles: the cosmic ray shadow of the Sun. We use Liouville's theorem to argue that the only effect of the solar magnetic field on the isotropic cosmic ray flux is to interrupt some of the trajectories that were aiming to the Earth and create a shadow. This shadow reveals the average solar depth crossed by cosmic rays of a given rigidity. The absorbed cosmic ray flux is then processed in the thin Solar surface and, assuming that the emission of neutral particles by low-energy charged particles is isotropic, we obtain (i) a flux of gammas that is consistent with Fermi-LAT observations, (ii) a flux of 100-300 neutrons/(year m^2) produced basically in the spallation of primary He nuclei, and (iii) a neutrino flux that is above the atmospheric neutrino background at energies above 0.1-0.6 TeV (depending on the solar phase and the zenith inclination). More precise measurements of the cosmic ray shadow and of the solar gamma flux, together with the possible discovery of the neutron and neutrino signals, would provide valuable information about the magnetic field, the cycle, and the interior of the Sun.

The production of very high energy muons inside an extensive air shower is observable at $\nu$ telescopes and sensitive to the composition of the primary cosmic ray. Here we discuss five different sources of these muons: pion and kaon decays; charmed hadron decays; rare decays of unflavored mesons; photon conversion into a muon pair; and photon conversion into a $J/\psi$ vector meson decaying into muons. We solve the cascade equations for a $10^{10.5}$ GeV proton primary and find that unflavored mesons and gamma conversions are the two main sources of $E\ge 10^{8.5}$ GeV muons, while charm decays dominate at $10^{5.5}\,{\rm GeV}< E< 10^{8.5}\,{\rm GeV}$. In inclined events one of these muons may deposite a large fraction of its energy near the surface, implying fluctuations in the longitudinal profile of the shower and in the muon to electron count at the ground level. In particular, we show that 1 out of 6 proton showers of $10^{10.5}$ GeV include an $E>10^6$ GeV deposition within 500 g/cm$^2$, while only in 1 out of 330 showers it is above $10^7$ GeV. We also show that the production of high energy muons is very different in proton, iron or photon showers ({e.g., conversions $\gamma\to \mu^+ \mu^-$ are the main source of $E\ge 10^4$ GeV muons in photon showers). Finally, we use Monte Carlo simulations to discuss the validity of our results.

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 perform spherically symmetric simulations of the core collapse of a single progenitor star of zero age main sequence mass $M_{\rm ZAMS} = 15 \, M_{\odot}$ with two models of heavy sterile neutrinos in the mass range of hundred MeV$/c^2$. According to both models, these hypothetical particles are copiously produced in the center, stream outwards a subsequently decay releasing energy into final states (including neutrinos) of the Standard Model. We find that they can lead to a successful explosion in otherwise non-exploding progenitors. Depending on their unknown parameters (e.g., mass and coupling constants with matter), we obtain either no explosion or an explosion of one of two types, i.e., through heating of gas downstream of the stalled shock wave, similarly to the standard scenario for supernova explosions or through heating of gas at higher radii that ejects matter from the outer core or the envelope while the center continues to accrete matter. In both cases, the explosion energies can be very high. We presume that this new type of explosion would produce an electromagnetic signal that significantly differs from common events because of the relative absence of heavy elements in the ejecta. The combination of core-collapse simulations and astrophysical observations may further constrain the parameters of the sterile neutrinos.

We describe several components in the diffuse flux of high energy neutrinos reaching the Earth and discuss whether they could explain IceCube's observations. Then we focus on TeV neutrinos from the Sun. We show that this solar neutrino flux is correlated with the cosmic-ray shadow of the Sun measured by HAWC, and we find that it is much larger than the flux of atmospheric neutrinos. Stars like our Sun provide neutrinos with a very steep spectrum and no associated gammas. We argue that this is the type of contribution that could solve the main puzzle presented by the high energy IceCube data.

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.

The Sun is a main source of high energy neutrinos. These neutrinos appear as secondary particles after the Sun absorbs high-energy cosmic rays, that find there a low-density environment (much thinner than our atmosphere) where most secondary pions, kaons and muons can decay before they lose energy. The main uncertainty in a calculation of the solar neutrino flux is due to the effects of the magnetic fields on the absorption rate of charged cosmic rays. We use recent data from HAWC on the cosmic-ray shadow of the Sun to estimate this rate. We evaluate the solar neutrino flux and show that at 1 TeV it is over ten times larger than the atmospheric one at zenith theta_z=30^o/150^o. The flux that we obtain has a distinct spectrum and flavor composition, and it could be detected in current and upcoming neutrino telescopes. KM3NeT, in particular, looks very promising: it will see the Sun high in the sky (the atmospheric background is lower there than near the horizon) and expects a very good angular resolution (the Sun's radius is just 0.27^o).

We calculate the fluxes of neutrinos and gamma rays from interactions of cosmic rays with interstellar matter in our galaxy. We use EPOS-LHC, SIBYLL and GHEISHA to obtain the yield of these particles in proton, helium and iron collisions at kinetic energies between 1 and 10^8 GeV, and we correlate the cosmic ray density with the mean magnetic field strength in the disk and the halo of our galaxy. We find that at E>1 PeV the fluxes depend very strongly on the cosmic-ray composition, whereas at 1-5 GeV the main source of uncertainty is the cosmic-ray spectrum out of the heliosphere. We show that the diffuse flux of galactic neutrinos becomes larger than the conventional atmospheric one at E>1 PeV, but that at all IceCube energies it is 4 times smaller than the atmospheric flux from forward-charm decays.

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.

We discuss the production of a class of heavy sterile neutrinos $\nu_h$ in proto-neutron stars. The neutrinos, of mass around $50$ MeV, have a negligible mixing with the active species but relatively large dimension-5 electromagnetic couplings. In particular, a magnetic dipole moment $\mu\approx 10^{-6}$ GeV$^{-1}$ implies that they are thermally produced through $e^+ e^-\to \bar \nu_h \nu_h$ in the early phase of the core collapse, whereas a heavy--light transition moment $\mu_{\rm tr}\approx 10^{-8}$ GeV$^{-1}$ allows their decay $\nu_h\to \nu_i \gamma$ with a lifetime around $10^{-3}$ s. This type of electromagnetic couplings has been recently proposed to explain the excess of electron-like events in baseline experiments. We show that the production and decay of these heavy neutrinos would transport energy from the central regions of the star to distances $d\approx 400$ km, providing a very efficient mechanism to enhance the supernova shock front and heat the material behind it.

Several cosmic-ray observatories have provided a high accuracy map of the sky at TeV--PeV energies. The data reveals an O(0.1%) deficit from north galactic directions that peaks at 10 TeV and then evolves with the energy, together with other anisotropies at smaller angular scales. Using Boltzmann's equation we derive expressions for the cosmic-ray flux that fit most of these features. The anisotropies depend on the local interstellar magnetic field B_{IS}, on the average galactic field B_{R} in our vicinity, and on the distribution of dominant cosmic-ray sources. We show that the initial dipole anisotropy along B_{IS} can be modulated by changes in the global cosmic ray wind, and that a variation in the dipole direction would imply a given radius of coherence for B_{IS}. We also show that small and medium-scale anisotropies may appear when the full-sky anisotropy finds a turbulence acting as a magnetic lens.

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.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

We discuss the possibility that a sterile neutrino of mass around 50 MeV slightly mixed with the muon flavor may be the origin of the MiniBooNE anomaly. We show that its production in the atmosphere in a fraction of kaon decays would imply an excess of contained showers at IceCube from down-going and near-horizontal directions.

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.

It has been recently proposed that a sterile neutrino \nu_h of mass m_h=40--80 MeV, mixing |U_{\mu h}|^2=10^{-3}--10^{-2}, lifetime \tau_h \lsim 10^{-9} s, and a dominant decay mode (\nu_h \to \nu_\mu \gamma) could be the origin of the experimental anomalies observed at LSND, KARMEN and MiniBooNE. Such a particle would be abundant inside air showers, as it can be produced in kaon decays (K -> \nu_h \mu, K_L -> \nu_h \pi \mu). We use the Z-moment method to evaluate its atmospheric flux and the frequency of its decays inside neutrino telescopes. We show that the \nu_h would imply around 10^4 contained showers per year inside a 0.03 km^3 telescope like ANTARES or the DeepCore in IceCube. These events would have a characteristic energy and zenith-angle distribution (E_\nu = 0.1--10 TeV and \theta < 90^o), which results from a balance between the reach of the heavy neutrino (that disfavors low energies) and a sizeable production rate and decay probability. The standard background from contained neutrino events (\nu_e N \to e X and neutral-current interactions of high inelasticity) is 100 times smaller. Therefore, although it may be challenging from an experimental point of view, a search at ANTARES and IceCube could confirm this heavy-neutrino possibility.

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.

Magnetic fields play a critical role in the propagation of charged cosmic rays. Particular field configurations supported by different astrophysical objects may be observable in cosmic ray maps. We consider a simple configuration, a constant azimuthal field in a disk-like object, that we identify as a cosmic magnetic lens. Such configuration is typical in most spiral galaxies, and we assume that it can also appear at smaller or larger scales. We show that the magnetic lens deflects cosmic rays in a regular geometrical pattern, very much like a gravitational lens deflects light but with some interesting differences. In particular, the lens acts effectively only in a definite region of the cosmic-ray spectrum, and it can be convergent or divergent depending on the (clockwise or counterclockwise) direction of the magnetic field and the (positive or negative) electric charge of the cosmic ray. We find that the image of a point-like monochromatic source may be one, two or four points depending on the relative position of source, observer and center of the lens. For a perfect alignment and a lens in the orthogonal plane the image becomes a ring. We also show that the presence of a lens could introduce low-scale fluctuations and matter-antimatter asymmetries in the fluxes from distant sources. The concept of cosmic magnetic lens that we introduce here may be useful in the interpretation of possible patterns observed in the cosmic ray flux at different energies.

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.

We consider the early universe at temperatures close to the fundamental scale of gravity (M_D << M_Planck) in models with extra dimensions. At such temperatures a small fraction of particles will experience transplanckian collisions that may result in microscopic black holes (BHs). BHs colder than the environment will gain mass, and as they grow their temperature drops further. We study the dynamics of a system (a black hole gas) defined by radiation at a given temperature coupled to a distribution of BHs of different mass. Our analysis includes the production of BHs in photon-photon collisions, BH evaporation, the absorption of radiation, collisions of two BHs to give a larger one, and the effects of the expansion. We show that the system may follow two different generic paths depending on the initial temperature of the plasma.

In a plasma at temperature close to the fundamental scale a small fraction of particles will experience transplanckian collisions that may result in microscopic black holes (BHs). We study the dynamics of a system (a black hole gas) defined by radiation at a given temperature coupled to a distribution of BHs of different mass. Our analysis includes the production of BHs in photon-photon collisions, BH evaporation, the absorption of radiation, collisions of two BHs to give a larger one, and the effects of the expansion. We find that the system may follow two different generic paths depending on the initial temperature of the plasma.

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 air-shower observatory Milagro has detected a large-scale anisotropy of unknown origin in the flux of TeV cosmic rays. We propose that this anisotropy is caused by galactic magnetic fields, in particular, that it results from the combined effects of the regular and the turbulent (fluctuating) magnetic fields in our vicinity. Instead of a diffusion equation, we integrate Boltzmann's equation to show that the turbulence may define a preferred direction in the cosmic-ray propagation that is orthogonal to the local regular magnetic field. The approximate dipole anisotropy that we obtain explains well Milagro's data.

We analyze the possibility that the cosmic ray knee appears at an energy threshold where the proton-dark matter cross section becomes large due to new TeV physics. It has been shown that such interactions could break the proton and produce a diffuse gamma ray flux consistent with MILAGRO observations. We argue that this hypothesis implies knees that scale with the atomic mass for the different nuclei, as KASKADE data seem to indicate. We find that to explain the change in the spectral index in the flux from E^{-2.7} to E^{-3.1} the cross section must grow like E^{0.4+\beta} above the knee, where \beta=0.3-0.6 parametrizes the energy dependence of the age (\tau \propto E^{-\beta}) of the cosmic rays reaching the Earth. The hypothesis also requires mbarn cross sections (that could be modelled with TeV gravity) and large densities of dark matter (that could be clumped around the sources of cosmic rays). We argue that neutrinos would also exhibit a threshold at E=(m_\chi/m_p)E_{knee}\approx 10^8 GeV where their interaction with a nucleon becomes strong. Therefore, the observation at ICECUBE or ANITA of standard neutrino events above this threshold would disprove the scenario.

The SHALON Cherenkov telescope has recorded over 2x10^6 extensive air showers during the past 17 years. The analysis of the signal at different zenith angles (\theta) has included observations from the sub-horizontal direction \theta=97^o. This inclination defines an Earth skimming trajectory with 7 km of air and around 1000 km of rock in front of the telescope. During a period of 324 hours of observation, after a cut of shower-like events that may be caused by chaotic sky flashes or reflections on the snow of vertical showers, we have detected 5 air showers of TeV energies. We argue that these events may be caused by the decay of a long-lived penetrating particle entering the atmosphere from the ground and decaying in front of the telescope. We show that this particle can not be a muon or a tau lepton. As a possible explanation, we discuss two scenarios with an unstable neutrino of mass m\approx 0.5 GeV and c\tau\approx 30 m. Remarkably, one of these models has been recently proposed to explain an excess of electron-like neutrino events at MiniBooNE.

TeV gravity models provide a scenario for black hole formation at energies much smaller than G_N^(-1/2) \sim 10^19 GeV. In particular, the collision of a ultrahigh energy cosmic ray with a dark matter particle in our galactic halo or with another cosmic ray could result into a black hole of mass between 10^4 and 10^11 GeV. Once produced, such object would evaporate into elementary particles via Hawking radiation. We show that the interactions among the particles exiting the black hole are not able to produce a photosphere nor a chromosphere. We then evaluate how these particles evolve using the jet-code HERWIG, and obtain a final diffuse flux of stable 4-dimensional particles peaked at 0.2 GeV. This flux consists of an approximate 43% of neutrinos, a 28% of electrons, a 16% of photons and a 13% of protons. Emission into the bulk would range from a 1.4% of the total energy for n=2 to a 16% for n=6.

In models with extra dimensions the fundamental scale of gravity M_D could be of order TeV. In that case the interaction cross section between a cosmic proton of energy E and a dark matter particle \chi will grow fast with E for center of mass energies \sqrt{2m_\chi E} above M_D, and it could reach 1 mbarn at E\approx 10^9 GeV. We show that these gravity-mediated processes would break the proton and produce a diffuse flux of particles/antiparticles, while boosting \chi with a fraction of the initial proton energy. We find that the expected cross sections and dark matter densities are not enough to produce an observable asymmetry in the flux of the most energetic (extragalactic) cosmic rays. However, we propose that unsuppressed TeV interactions may be the origin of the knee observed in the spectrum of galactic cosmic rays. The knee would appear at the energy threshold for the interaction of dark matter particles with cosmic protons trapped in the galaxy by \muG magnetic fields, and it would imply a well defined flux of secondary antiparticles and TeV gamma rays.

We study the properties of black holes of mass 10^4--10^{11} GeV in models with the fundamental scale of gravity at the TeV. These black holes could be produced in the collision of a ultrahigh energy cosmic ray with a dark matter particle in our galactic halo or with another cosmic ray. We show that QCD bremsstrahlung and pair production processes are unable to thermalize the particles exiting the black hole, so a chromosphere is never formed during Hawking evaporation. We evaluate with HERWIG the spectrum of stable 4-dim particles emitted during the Schwarzschild phase and find that in all cases it is peaked at energies around 0.2 GeV, with an approximate 43% of neutrinos, 28% of photons, 16% of electrons and 13% of protons. Bulk gravitons are peaked at higher energies, they account for a 0.4% of the particles (16% of the total energy) emitted by the most massive black holes in n=6 extra dimensions or just the 0.02% of the particles (1.4% of the energy) emitted by a 10 TeV black hole for n=2.

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.

Models where supersymmetry (SUSY) is manifest only in a sector of the low-energy spectrum have been recently proposed as an alternative to the MSSM. In these models the electroweak scale is explained by a fine-tuning between different Higgs mass contributions (split-SUSY models), or by the localization of the Higgs sector in a point of an extra dimension where all the mass parameters are suppressed by the metric (partly-SUSY models). Therefore, the presence of a good dark matter candidate becomes the main motivation for (partial) low-energy SUSY. We study this issue in minimal frameworks where the higgsinos are the only light supersymmetric particles. Whereas in split-SUSY models the higgsino should have a mass around 1 TeV, we show that in partly-SUSY models the lightest higgsino could also be found below MW.

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.

We examine the experimental signatures for the production of gluinos at colliders and in cosmic rays within the split supersymmetry scenario. Unlike in the MSSM, the gluinos in this model are relatively long-lived due to the large value of the squark masses which mediate their decay. Searches at colliders are found to be sensitive to the nature of gluino fragmentation as well as the gluino-hadron interactions with nuclei and energy deposition as it traverses the detector. We find that the worst-case scenario, where a neutral gluino-hadron passes through the detector with little energy deposition, is well described by a monojet signature. For this case, using Run I data we obtain a bound of $m_{\tilde g} > 170$ GeV; this will increase to 210(1100) GeV at Run II(LHC) if no excess events are observed. In the opposite case, where a charged gluino-hadron travels through the detector, a significantly greater reach is obtained via stable charged particle search techniques. We also examine the production of gluino pairs in the atmosphere by cosmic rays and show they are potentially observable at IceCube; this would provide a cross-check for observations at hadron colliders.

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.