###### Abstract

TeV gravity models provide a scenario for black hole formation at energies much smaller than 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 and 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 of neutrinos, a of electrons, a of photons and a of protons. Emission into the bulk would range from a of the total energy for to a for .

Signatures of TeV gravity from the evaporation
of cosmogenic black holes^{1}^{1}1Talk presented by I.M. at the Workshop on Black Holes in General Relativity and String Theory,
August 24-30 2008, Veli Lošinj, Croatia.

Iacopo Mastromatteo, Petros Draggiotis, Manuel Masip

Dipartimento di Fisica Teorica

Università degli Studi di Trieste, I-34014 Trieste, Italy

CAFPE and Departamento de Física Teórica y del Cosmos

Universidad de Granada, E-18071 Granada, Spain

mastroma@sissa.it, pdrangiotis@ugr.es, masip@ugr.es

## 1 Introduction

Models with extra dimensions [1] provide one of the most promising solutions to the hierarchy problem, namely, the huge difference between the scale of gravity GeV and the electroweak (EW) scale GeV. In these models appears as an effective scale related with the fundamental one, 1–10 TeV, by the volume of the compact space or by an exponential warp factor. The difference between and would then just define a little hierarchy problem that should be easier to solve consistenly with all collider data. The phenomenological consequencies of this framework are quite intriguing: the fundamental scale would be at accessible energies, and processes with would probe a transplanckian regime where gravity is expected to dominate over the other interactions [2]. The spin two of the graviton implies then gravitational cross sections that grow fast with and become long distance interactions. As a consequence, quantum gravity or other short distance effects become irrelevant as they are screened by black hole (BH) horizons [3].

One of the scenarios in which TeV gravity effects could play a significant role is provided by cosmic rays physics. The Earth is constantly hit by a flux of protons with energy of up to GeV and, associated to that flux, it is also expected a flux of cosmogenic neutrinos (still unobserved) with a typical energy peaked around GeV [4]. These are energies much larger than the ones to be explored at the LHC, where there would be no evidence for gravitational interactions if the scale is above a few TeV. In addition, notice that the new physics should be more relevant in collisions of particles with a small SM cross section, as it is expected for the interaction of a proton with a dark matter particle if it is taken to be a weakly interacting massive particle (WIMP). We will discuss here the interaction of ultra high energy cosmic rays (UHECR) with dark matter particles in our galactic halo. No detail about the nature of other than its mass, which defines the center-of-mass energy in the collision, is going to be significant to the present analysis. We also consider collisions of UHECR with other cosmic rays. These are arguably the most energetic elementary processes that we know that occur in nature at the present time, and would produce mini BHs significantly colder and longer-lived than the ones usually considered in the literature. We will focus just on BH production and evaporation, being this analysis a necessary first step in order to understand the full effects of TeV gravity on UHECR phenomenology.

## 2 Cosmogenic black hole production

BH production processes are the most widely and detailfully discussed aspect of TeV-gravity phenomenology [5], and they have been considered both in the LHC [6] and in the UHECR context [7]. Here we will assume a scenario with flat extra dimensions of common lenght where gravity is free to propagate, while matter fields are trapped on a (non-compact) four-dimensional brane. We will use the basic estimate that the collision of two pointlike particles at impact parameters smaller than the Schwarzschild radius of the system leads to the production of a BH whose mass is given by . The BHs that we are considering ( GeV) will be described by a ()-dimensional metric (they are smaller than the volume of the compact space), being their radius

(1) |

For two pointlike particles, the cross section to produce a BH is then written as

(2) |

If the collision involves non-elementary (at the scale ) protons, then its partonic structure has to be included in order to find the total cross section, as it is usually done for analyses of BH production at LHC [6]. The p– (or p–) cross section may therefore be written as

(3) |

This formula expresses the cross section as the sum of
partial contributions to produce a BH of mass
resulting from the collision of a parton
that carries a fraction of momentum with a pointlike target.
It is crucial to notice that the scale in the collision
is fixed by the inverse Scwarzschild radius, rather than by the
BH mass [3] [8], since the
scattering is probing a
lenght scale that grows (not decreases!) with . Actually, we
expect that for large enough the scale that we are exploring goes
above its radius and a pointlike behaviour for the proton
will emerge. In contrast with a QED scattering, here at lower
energies
( GeV) we can see the composite structure
of the proton, while at higher energies ( GeV)
the proton will scatter coherently as a whole.
Since Eq. 3 does not reproduce this behaviour,
it is necessary to include matching corrections between the two
energy regions. The cross section in Eq. 3 describes
the low-energy regime, and it is dominated by the large number of
partons of low that may produce a BH of mass near the threshold
. This scheme explains why
.
When the cross section approaches the proton size
( mbarn), then the
density of partons with enough energy to produce a BH is so large
that the parton cross sections overlap, and
the BHs produced are big enough to trap
other spectator partons. This overlapping reduces
the total cross section and increases the average
mass of the produced BH. In this regime
is basically constant with until it matches
the pointlike behaviour in .
A similar behavior is also expected in – collisions,
where the partonic
enhancement of the cross section is even more important
at lower energies (in this regime
) and the
intermadiate regime of constant total cross section is reached
at lower energies.
The smooth transition from these regimes can be
modelled numerically discounting the
contributions from spectator partons, and are
summarized in Fig. 1. There we plot
the BH production cross section for different kind of
particles^{2}^{2}2We assumed a CTEQ6M
set of PDF [9].

We will analyze two processes that can lead to BH production (see [10] for the fluxes of proton, cosmogenic neutrinos and for the dark matter density).

(i) A cosmic ray of energy colliding with a dark matter particle at rest in the frame of reference of our galaxy. The average number of BHs produced per unit time and volume depends on the density , the cross section and the differential flux of cosmic rays (with ):

(4) |

Here the center of mass energy can run from to GeV.

(ii) A cosmic ray of energy colliding with a cosmic ray of energy . In this case the center of mass energy depends upon the relative angle , and results into . The interaction rate per unit time and volume is expressed by:

(5) |

These processes generate BH masses that can reach GeV.

In Fig. 2 we plot the production rate of BHs from both types of collisions.

## 3 Black hole evaporation

To understand what kind of signal one could observe from such an event, it is necessary to estimate how the BH evolves after its production. It is expected that initially the BH undergoes a quick balding phase, in which it loses its gauge hair and asymmetries. Then it experiences a spin down phase, where its angular momentum is radiated while losing just a small fraction of its mass [11]. Finally, during most of its life the BH is in a Schwarzschild phase, losing mass through spherically symmetric Hawking radiation [12]. The spectrum is, in a first approximation, that of a black body of temperature [13]

(6) |

This means that the scale of emission is fixed by the inverse Schwarzshild radius. This formula has important corrections arising from the gravitational barrier that the particles have to cross once emitted. These corrections are usually expressed in terms of the so called greybody factors, effective emission areas that depend on the dimensionality of the space-time, the spin of the particle emitted, and its energy [14]. These factors give corrections of order 1 to the black-body emission rates for all particles species except for the graviton, which can have a stronger correction depending upon the number of extra dimensions. We will assume here the numerical greybody factors given in [15].

The number of particles of the species emitted with ()-dimensional momenta between and in a time interval can be written as

(7) |

while the radiated energy is given by

(8) |

Some remarks are here in order.

(i) On dimensional grounds and , so each degree of freedom should contribute equally (up to order one geometric and greybody factors) to the total emission independently from its bulk or brane localization [16].

(ii) We are considering BH temperatures above ( GeV leads to GeV), so QCD degrees of freedom (quarks and gluons) are also radiated and dominate the total emission.

Once the instant spectrum is known, we integrate it over time to get the BH lifetime. On dimensional grounds , although the dependence upon the number of the radiated degrees of freedom at different temperatures may be significant. In Fig. 3 we plot the correlation between lifetime, mass and initial temperatures for BHs of mass ranging from TeV to GeV, and TeV; it is there shown that lifetimes go from a maximum of s for the most heavy BHs to a minimum around s for LHC-like BHs.

## 4 Thermal properties of the radiation

An important issue about the evolution of the radiation from a BH is the debated question relative to its thermalization. It has been argued [17] that the emitted particles should produce a thick shell of almost-thermal plasma of QED (QCD) particles usually called photosphere (chromosphere). This would occur for BHs above a critical temperature (), and would change the average energy of the emitted particles from to (or ). The argument leading to these shells is based on the average number of interactions of the particles exiting the BH, so should suffice to confirm the presence of the plasma shell. Initial estimates [17] used the expression

(9) |

which describes the case of particles scattering against a fixed target.
Recently, however, it has been noticed that the kinematic
differences between that case and the case of particles
exiting radially from a BH are so significant that lead to a
complete suppression of the interaction rate [18].
We will show, following the approach of Carr, McGibbon and Page,
that their arguments^{3}^{3}3These arguments are supported
by the numerical analysis in [19]. (formulated
for ordinary BHs) hold also for BHs in TeV gravity models.

The first kinematic effect is due to causality, and depends on the fact that the scattered particles do not come from infinity (as in a regular collision), they are created in definite points of space-time. This introduces a minimal separation between particles successively emitted, both in time and lenght, and induces via Heisenberg’s indetermination principle an UV cutoff in the scale of the exchanged momenta. The scattering cross section is reduced because not all the energies can be interchanged. In particular, in QED (QCD) Bremsstrahlung and pair production the momenta dominating the collision are of order (or right above ). If the particle wave functions do not overlap, and their minimum distance (in units of their Compton wavelength) is larger than the dominant inverse momenta, then the process will be suppressed. Checking the parameter is sufficient to decide about the effective connection between emitted particles, and eventually exclude thermal interactions. In [10] We have shown that this argument excludes the presence of a photosphere for any number of extra dimensions, but not of a chromosphere when .

The second suppression effect is based on the fact that the interaction between two particles is not instantaneous, it takes a finite time to complete. It is easy to see that when this occurs the particles are already far away from each other, so that they can not interact again. In particular, after completing a QCD interaction partons will be at distances larger than , where QCD is already ineffective. To fully understand this point, one first has to notice that the interaction between particles moving radially in the same direction (within the exclusion cone) is negligible, as the density in such a region is low. Also, that particles moving radially keep on moving radially, as the average angular deviation due to Bremsstrahlung-like processes is small. This implies that the distance of a particle to the particles out of the exclusion cone will always increase (they never approach to each other), and when it reaches a radius this distance will be bigger than (or ) and the particle is no longer able to interact. If the BH temperature is above , as it is the case for the BHs under study here, it is easy to see that after the particle has completed one interaction it will have already crossed .

## 5 Stable particle spectrum

Once the greybody spectrum of emission has been established, it is necessary to study how it evolves at astrophysical distances: unstable particles will decay, and colored particles (which dominate the spectrum) will fragment into hadrons and then shower into stable species. We present our results following the approach of [20], who first studied this issue for primordial BHs. The main difference with their analysis is that while the authors in [20] compute the stationary spectrum at a given (which only changes on astrophysical time scales), we need here to evaluate the spectrum integrated over the whole (very short) BH lifetime. In any case, our results will be analogous, since the temperature of a BH variates little for most of its lifetime. Of course, our framework also deals with a different scale of gravity and extra dimensions where gravitons propagate. This implies emission into the bulk and different greybody factors for all the species. Notice finally that the spectrum that we are discussing is in the BH rest frame, it is not the one to be observed at the Earth as the BHs produced in cosmic ray collisions will be highly boosted.

We will assume that the evolution of the species emitted by a BH at rest coincides with the one in in the center of mass frame, so we will use the MonteCarlo jet code HERWIG6 [21] to evolve the greybody spectrum described before. Namely, we compute the convolution

(10) |

to obtain the number of stable particles of species with energy between and emitted in a time . The first term in parenthesis stands for the greybody spectrum of emission for particle species , while the second encodes the probability for the species of energy to give a of energy . For a given , this has been implemented via MonteCarlo including all particles of mass (leptons, quarks and gauge bosons, neglecting the Higgs or the dark matter particle) and has resulted in a final spectrum of neutrinos, electrons, photons and protons. The spectrum includes the same number of particles and antiparticles (they are generated at the same rate), and the three families of neutrinos (their flavor oscillates at astrophysical lenght scales). In Fig. 4 we plot the spectrum at fixed temperature GeV, whereas in Fig. 5 we give the complete spectrum for initial masses of GeV and and .

The results can be summarized as follows.

(i) The main product of the emission is constituted by particles resulting from the showering of QCD species; this explains the primary peak at GeV observed in the spectrum. It is also possible to detect at the direct greybody emission as a secondary peak. Gravitons decouple, since they are not produced by decay of unstable species.

(ii) The relative emissivities of Standard Model particles are an approximate of neutrinos, a of electrons, a of photons and a of protons. This is only mildly sensitive to the BH mass or , as it is determined by the showering of colored particles.

(iii) Emission into the bulk goes from the 0.4% of the total number of particles (16% of the total energy) emitted if and GeV to the 0.02% of the particles (1.4% of the energy) emitted for and GeV.

## 6 Summary and outlook

The head to head collision of two cosmic rays provides center-of-mass energies of up to GeV. In models with extra dimensions and a fundamental scale of gravity at the TeV such collision should result in the formation of a mini BH. Its evaporation and showering into stable particles could provide an observable signal.

We have estimated the production rate of these BHs (Fig. 2) via collisions of two cosmic rays or, more frequently, in the collision of a cosmic ray and a dark matter particle. In particular, it seems worth to analyze the possibility that (i) extragalactic cosmic rays crossing the galactic DM halo produce a flux of secondary particles with a characteristic shape and strongly dependent upon galactic latitude; (ii) a fraction of the flux of cosmic rays with energy up to GeV trapped in our galaxy by magnetic fields can be processed by TeV interactions into a secondary flux peaked at smaller energies. Notice that the physics proposed in this talk is expected to become relevant just at center of mass energies above TeV, i.e., at cosmic ray energies around the cosmic ray knee. These considerations will be worked out in [22], where the additional effects of gravitational elastic interactions will also be included.

Here we have discussed the properties of BHs with masses between and GeV. Such objects have a proper lifetime between and s (Fig. 3), and their desintegration products are mainly determined by the fragmentation of QCD species produced via Hawking radiation. Interactions among emitted particles are not able to produce a thermal shell of radiation, so the spectrum of fundamental species exit the BH with basically the black body spectrum described by Hawking. The final spectrum of stable particles at large distances, however, is peaked around , and exhibits features weakly dependent upon number of extra dimensions or the BHs mass. Standard Model modes are constituted by an approximate of neutrinos, a of electrons, a of photons and a of protons. The gravitons produced are a fraction that goes from the of the total number of particles ( of the energy) for GeV and to the ( of the energy) for TeV and .

This work is a preliminary analysis, with results that can be useful for future search for effects of TeV gravity on cosmic ray physics.

## References

- [1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429 (1998) 263; N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Rev. D 59 (1999) 086004; I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436 (1998) 257; L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370.
- [2] R. Emparan, Phys. Rev. D 64 (2001) 024025.
- [3] S. B. Giddings and S. D. Thomas, Phys. Rev. D 65 (2002) 056010.
- [4] D. V. Semikoz and G. Sigl, JCAP 0404 (2004) 003.
- [5] P. D. D’Eath and P. N. Payne, Phys. Rev. D 46, 658 (1992); P. D. D’Eath and P. N. Payne, Phys. Rev. D 46, 675 (1992); P. D. D’Eath and P. N. Payne, Phys. Rev. D 46, 694 (1992); T. Banks and W. Fischler, [hep-th/9906038]; D. M. Eardley and S. B. Giddings, Phys. Rev. D 66, 044011 (2002); H. Yoshino and Y. Nambu, Phys. Rev. D 66, 065004 (2002); H. Yoshino and Y. Nambu, Phys. Rev. D 67, 024009 (2003).
- [6] S. Dimopoulos and G. L. Landsberg, Phys. Rev. Lett. 87 (2001) 161602; K. Cheung, Phys. Rev. D 66 (2002) 036007; T. G. Rizzo, JHEP 0202, 011 (2002); A. Chamblin and G. C. Nayak, Phys. Rev. D 66, 091901 (2002); R. Casadio and B. Harms, Int. J. Mod. Phys. A 17 (2002) 4635; M. Cavaglia, Phys. Lett. B 569 (2003) 7; I. Mocioiu, Y. Nara and I. Sarcevic, Phys. Lett. B 557, 87 (2003); M. Cavaglia, S. Das and R. Maartens, Class. Quant. Grav. 20, L205 (2003); C. M. Harris, M. J. Palmer, M. A. Parker, P. Richardson, A. Sabetfakhri and B. R. Webber, JHEP 0505, 053 (2005); J. Tanaka, T. Yamamura, S. Asai and J. Kanzaki, Eur. Phys. J. C 41, 19 (2005); J. L. Hewett, B. Lillie and T. G. Rizzo, Phys. Rev. Lett. 95, 261603 (2005); B. Koch, M. Bleicher and S. Hossenfelder, JHEP 0510, 053 (2005); D. M. Gingrich, Int. J. Mod. Phys. A 21 (2006) 6653; T. J. Humanic, B. Koch and H. Stoecker, Int. J. Mod. Phys. E 16, 841 (2007); T. G. Rizzo, Phys. Lett. B 647, 43 (2007); M. Cavaglia, R. Godang, L. M. Cremaldi and D. J. Summers, JHEP 0706, 055 (2007); G. Dvali and S. Sibiryakov, JHEP 0803 (2008) 007; P. Meade and L. Randall, arXiv:0708.3017 [hep-ph]; M. Cavaglia, Int. J. Mod. Phys. A 18 (2003) 1843; G. L. Alberghi, R. Casadio, D. Galli, D. Gregori, A. Tronconi and V. Vagnoni, hep-ph/0601243.
- [7] J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88 (2002) 021303; L. Anchordoqui and H. Goldberg, Phys. Rev. D 65 (2002) 047502; L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 65 (2002) 124027; A. Ringwald and H. Tu, Phys. Lett. B 525 (2002) 135; M. Kowalski, A. Ringwald and H. Tu, Phys. Lett. B 529 (2002) 1; S. I. Dutta, M. H. Reno and I. Sarcevic, Phys. Rev. D 66 (2002) 033002; J. Alvarez-Muniz, J. L. Feng, F. Halzen, T. Han and D. Hooper, Phys. Rev. D 65 (2002) 124015; L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 68 (2003) 104025; J. I. Illana, M. Masip and D. Meloni, Phys. Rev. Lett. 93 (2004) 151102; E. J. Ahn, M. Cavaglia and A. V. Olinto, Astropart. Phys. 22 (2005) 377; J. I. Illana, M. Masip and D. Meloni, Phys. Rev. D 72 (2005) 024003.
- [8] R. Emparan, M. Masip and R. Rattazzi, Phys. Rev. D 65, 064023 (2002).
- [9] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 012 (2002).
- [10] P. Draggiotis, M. Masip and I. Mastromatteo, JCAP 0807, 014 (2008).
- [11] E. Winstanley, arXiv:0708.2656 [hep-th].
- [12] S. W. Hawking, Nature 248 (1974) 30; S. W. Hawking, Commun. Math. Phys. 43 (1975) 199, [Erratum-ibid. 46 (1976) 206].
- [13] R. C. Myers and M. J. Perry, Annals Phys. 172, 304 (1986).
- [14] D. N. Page, Phys. Rev. D 13 (1976) 198; Phys. Rev. D 14 (1976) 3260; Phys. Rev. D 16 (1977) 2402.
- [15] C. M. Harris and P. Kanti, JHEP 0310 (2003) 014; P. Kanti and J. March-Russell, Phys. Rev. D 66 (2002) 024023; P. Kanti and J. March-Russell, Phys. Rev. D 67 (2003) 104019; V. Cardoso, M. Cavaglia and L. Gualtieri, Phys. Rev. Lett. 96 (2006) 071301, [Erratum-ibid. 96 (2006) 219902].
- [16] R. Emparan, G. T. Horowitz and R. C. Myers, Phys. Rev. Lett. 85 (2000) 499.
- [17] A. F. Heckler, Phys. Rev. D 55 (1997) 480.
- [18] J. H. MacGibbon, B. J. Carr and D. N. Page, arXiv:0709.2380 [astro-ph]; D. N. Page, B. J. Carr and J. H. MacGibbon, arXiv:0709.2381 [astro-ph].
- [19] C. Alig, M. Drees and K. y. Oda, JHEP 0612, 049 (2006).
- [20] J. H. MacGibbon and B. R. Webber, Phys. Rev. D 41 (1990) 3052.
- [21] G. Corcella et al., JHEP 0101 (2001) 010.
- [22] M. Masip and I. Mastromatteo, JCAP 0812, 003 (2008).