%Format: Latex Orientation: Portrait
\documentclass[12pt,dvips]{article}
%\usepackage{dina4p}
\usepackage[dvips]{graphicx}
\usepackage{graphicx}

\DeclareGraphicsRule{.eps.gz}{eps}{.eps.bb}{`gunzip -cd #1}
\graphicspath{{plots/}}
\DeclareGraphicsExtensions{.eps.gz,.eps,.ps,.ps.gz}

%
\oddsidemargin-5mm
\evensidemargin-5mm

\setcounter{secnumdepth}{7}
\setcounter{tocdepth}{7} 
\parskip=\itemsep               %?
\setlength{\itemsep}{0pt}       %?
\setlength{\partopsep}{0pt}     %?
\setlength{\topsep}{0pt}        %?
%%---layout fuer eine dina4 seite-------------------
\setlength{\textheight}{22.5cm}
\setlength{\textwidth}{174mm}
\setlength{\topmargin}{-1.5cm}
%
\setlength{\marginparwidth}{0cm}
\setlength{\marginparsep}{0cm}
\setlength{\parindent}{0cm}
\renewcommand{\textfraction}{0}  
\renewcommand{\topfraction}{1}
\renewcommand{\bottomfraction}{1}
\renewcommand{\baselinestretch}{1}

\newcommand{\lwig}{\mbox{\,\raisebox{.3ex}
    {$<$}$\!\!\!\!\!$\raisebox{-.9ex}{$\sim$}\,}}
\newcommand{\gwig}{\mbox{\,\raisebox{.3ex}
    {$>$}$\!\!\!\!\!$\raisebox{-.9ex}{$\sim$}}\,}
\newcommand{\lambdabar}{{\hbox{$\lambda_e$\kern-1.9ex\raise+0.45ex\hbox{--}
\kern+0.2ex}}}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%BEGINNING OF TEXT                           
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\date{\empty}

\title{{\normalsize\rightline{DESY 02-001}\rightline}
\vskip 1cm 
\bf Black Holes at Neutrino Telescopes
       \vspace{21mm}} 
\author{
M. Kowalski\\[4mm]
Deutsches Elektronen-Synchrotron DESY, Zeuthen, Germany\\[4mm]
A. Ringwald, H. Tu\\[4mm] 
Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany}

\begin{document}
\begin{titlepage} 
  \maketitle
% declarations for front matter
%\vspace{1.cm}
\begin{abstract}
In scenarios with extra dimensions and TeV-scale quantum gravity,
black holes are expected to be produced in the collision of light particles 
at center-of-mass energies above the fundamental Planck scale with small impact parameters.
Black hole production and evaporation may thus be studied in detail at the  
Large Hadron Collider (LHC).
But even before the LHC starts operating, 
neutrino telescopes such as  AMANDA/IceCube,  ANTARES,  
Baikal, and RICE have an opportunity to search for 
black hole signatures. Black hole production in the scattering of 
ultrahigh energy cosmic neutrinos on nucleons in the ice or water  may initiate 
cascades and through-going muons with distinct characteristics above the Standard Model rate. 
In this Letter, we investigate the sensitivity of neutrino telescopes  
to black hole production and compare it to the one expected at the Pierre Auger Observatory, 
an air shower array currently under construction, and at the LHC. 
We find that, already with the currently available data, AMANDA and RICE should be able
to place sensible constraints in black hole production parameter space, which are 
competitive with the present ones from the air shower facilities Fly's Eye and AGASA.
In the optimistic case that a ultrahigh energy cosmic neutrino flux significantly higher than the 
one expected from cosmic ray interactions with the cosmic microwave background radiation
is realized in nature, one even has discovery potential for black holes at neutrino telescopes 
beyond the reach of LHC.     
\end{abstract}

% typeset front matter (including abstract)

\thispagestyle{empty}
\end{titlepage}
\newpage \setcounter{page}{2}

%\section{Introduction}
{\em 1.}
The possibility that black holes will be produced in the 
collision of two light particles at center-of-mass (cm) energies above the Planck scale
with small impact parameters~\cite{'tHooft:1987rb}  
seems now within reach  in the context of theories with
$\delta = D-4\geq 1$ large compact~\cite{Arkani-Hamed:1998rs} 
or warped~\cite{Randall:1999ee} extra dimensions 
and a low fundamental Planck scale 
$M_D\,\gwig$ 1 TeV characterizing quantum 
gravity.
In these theories one might expect the copious production of black holes 
in high energy collisions
at cm energies above 
$M_D$~\cite{Argyres:1998qn}.
%
Correspondingly, the Large Hadron Collider (LHC)~\cite{Evans:2001mn}, 
expected  to have a first physics run in  2006, may turn into a factory of black holes
at which their production and evaporation may be studied in 
detail~\cite{Giddings:2001bu,Dimopoulos:2001hw,Dimopoulos:2001qe,Hossenfelder:2001dn}. 

Black hole production and subsequent decay in the scattering of ultrahigh energy cosmic neutrinos
on nucleons in the atmosphere may initiate quasi-horizontal air showers 
far above the Standard Model rate~\cite{Feng:2002ib,Emparan:2001kf} and with distinct 
characteristics~\cite{Anchordoqui:2001ei}. The reach 
of the Pierre Auger Observatory~\cite{Zavrtanik:2000zi} for extensive air showers, 
expected to be completed by the end of 2003,
in the black hole production parameter 
range has been thoroughly investigated in Refs.~\cite{Feng:2002ib,Ringwald:2001vk,Anchordoqui:2001cg}.  
It was found that, depending on the size of the flux of ultrahigh energy cosmic neutrinos, 
the first signs of black hole production may be observed at Auger already before
the start of the LHC. Furthermore, it was shown that already now sensible constraints
on black hole production can be obtained from the non-observation of horizontal 
showers~\cite{Baltrusaitis:1985mt,Yoshida:2001}
by the Fly's Eye collaboration~\cite{Baltrusaitis:1985mx} and the 
Akeno Giant Air Shower Array (AGASA) collaboration~\cite{Chiba:1992nf}, 
respectively.  These constraints turn out to be 
competitive with other currently available constraints on TeV-scale gravity which are mainly based on interactions associated with
Kaluza-Klein gravitons, 
according to which a fundamental Planck scale as low as $M_D = {\mathcal O}(1)$ TeV is still allowed for $\delta\geq 4$ flat 
or $\delta\geq 1$ warped extra dimensions~\cite{Peskin:2000ti}\footnote{For an exhaustive list of references in this 
context, see also Ref.~\cite{Ringwald:2001vk}.}. 

In this Letter, we want to consider the sensitivity of neutrino telescopes such as 
AMANDA/Ice\-Cube~\cite{Halzen:1999jy}, 
ANTARES~\cite{Montaruli:2002gh},  Baikal~\cite{Domogatsky:2001br}, NESTOR~\cite{Grieder:2001kr}, 
and RICE~\cite{Kravchenko:2001id}
to black hole production in the scattering of ultrahigh energy cosmic neutrinos on 
nucleons in ice or water, considerably extending  a first exploratory
study~\cite{Uehara:2001yk}.
After briefly reviewing the phenomenological model of black hole 
production adopted, we calculate the rate expectation for events contained in the 
fiducial volume of a neutrino telescope.   
We then extend our calculations to rates from through-going muons.
Finally, we discuss the achievable sensitivity for detection of 
black hole production and compare it with that from air-shower 
facilities and the LHC.

%\section{Black hole phenomenology}
{\em 2.}
Our phenomenological model of black hole production and decay will be based on the 
following working hypothesis which summarizes the current understanding of the 
production and decay of black holes in TeV-scale gravity 
scenarios~\cite{Argyres:1998qn,Emparan:2000rs,Giddings:2000ay,Giddings:2001bu,Dimopoulos:2001hw,%
Dimopoulos:2001qe,Hossenfelder:2001dn}.

It is assumed that 
at trans-Planckian parton-parton cm energies squared, $\hat s\gg M_D^2$,  
and at small parton-parton impact parameters, $b\ll r_S\,(M_{\rm bh}=\sqrt{\hat s})$, i.\,e. 
at impact parameters much smaller than the Schwarzschild radius $r_S$ of a $(4+\delta )$-dimensional black 
hole with mass $M_{\rm bh}=\sqrt{\hat s}$~\cite{Myers:1986un}\footnote{We define $M_D$ as in 
Ref.~\cite{Giudice:1999ck}. 
Equation~(\ref{schwarzsch}) is valid as long as $r_S\ll R_c$, with $R_c$ being the compactification 
or curvature radii in the flat or warped scenario, respectively.}, 
\begin{equation}
\label{schwarzsch}
r_S =\frac{1}{M_D}
\left[
\frac{M_{\rm bh}}{M_D}
\left(
\frac{2^\delta \pi^{\frac{\delta -3}{2}}\,\Gamma\left( \frac{3+\delta}{2}\right)}{2+\delta}
\right)
\right]^{\frac{1}{1+\delta}}
\,,
\end{equation} 
a black hole forms with a cross section%\footnote{}
\begin{equation}
\label{sig_bh_geom}
\hat\sigma (ij\to {\rm bh})\equiv  
\hat\sigma^{\rm bh}_{ij} (\hat s ) \approx \pi\,r_S^2 
\left( M_{\rm bh}=\sqrt{\hat s} \right)\,
\theta\left( \sqrt{\hat s} -M_{\rm bh}^{\rm min}\right) 
\,.
\end{equation}
Here, $M_{\rm bh}^{\rm min}\gg M_D$ parametrizes the cm energy above which the
semiclassical reasoning mentioned above is assumed to be valid.
%
The contribution of black hole production to the neutrino-nucleon  cross section is then
obtained as the convolution 
\begin{equation}
\label{sig_nuN_bh_pdf}
\sigma_{\nu N}^{\rm bh} (s) =
\sum_i \int^1_0%_{M_{\rm bh}^{{\rm min}\,2}/s}^1
{\rm d}x\,f_i (x,\mu )\,\hat \sigma_{\nu i}^{\rm bh} (xs)
\, .
\end{equation}
Here, $s=2\,m_N\,E_\nu$, with $m_N=(m_p+m_n)/2$ denoting an isoscalar nucleon mass and $E_\nu$
the neutrino energy in the laboratory frame, denotes the 
neutrino-nucleon cm  energy squared. The sum 
extends over all partons in the nucleon, with parton distribution functions $f_i(x,\mu )$ and 
factorization scale $\mu$. As in Ref.~\cite{Ringwald:2001vk}, we have used in the calculation of the
cross section~(\ref{sig_nuN_bh_pdf}) various sets of parton 
distributions as they are implemented in the parton distribution library 
PDFLIB~\cite{Plothow-Besch:1993qj}. In the following we display only the results obtained with 
the CTEQ3D~\cite{Lai:1995bb} parton distributions with $\mu = {\rm min}\,(\sqrt{\hat s},10\ {\rm TeV})$, 
keeping in mind that uncertainties associated with different parton distribution sets 
are in the ${\mathcal O}(20)\, \%$ range, whereas a choice of the other natural factorization scale 
$\mu = r_S^{-1}$~\cite{Giddings:2001bu,Emparan:2001kf} decreases the predicted rates
by a factor of ${\mathcal O}(2)$.

It is fair to say that there remain nevertheless substantial uncertainties in the size of the cross section.
As noted in Ref.~\cite{Voloshin:2001vs}, the semiclassical production of black
holes resembles largely the problem of baryon and lepton number violating processes 
(``sphaleron~\cite{Klinkhamer:1984di} production'') in multi-TeV ($\sqrt{\hat s}\gg m_W/\alpha_W$) 
particle collisions in the standard electroweak theory~\cite{Ringwald:1990ee,Morris:1991bb,Gibbs:1995cw} 
and more generally the problem of multi-particle production in weakly coupled theories~\cite{Cornwall:1990hh}, 
which is not yet completely understood~\cite{Mattis:1992bj}.  
There might be an additional exponential suppression factor rendering semiclassical 
sphaleron, multi-particle or black hole production unobservable in the TeV range~\cite{Voloshin:2001vs} 
(see, however, Ref.~\cite{Dimopoulos:2001qe}). 
%
Nevertheless, in this Letter we shall confine ourselves  to the simple geometric cross section 
estimate~(\ref{sig_nuN_bh_pdf}) and postpone an inclusion of a possible exponential suppression factor 
to a longer write-up~\cite{Kowalski:inprep}. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%FIGURE%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\begin{center}
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=570pt,bbury=608pt,%
width=12cm]{mul_Md1_Mmin5.eps}
\caption[dum]{\label{multiplicity}
The average number $\langle n\rangle$ of Standard Model particles~(\ref{eq:av_mult}), into which a 
black hole produced in neutrino-nucleon scattering decays, as a function of the
incident neutrino energy in the laboratory frame, $E_\nu$. 
}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

After production, the black holes decay primarily 
via Hawking radiation~\cite{Hawking:1974rv}, corresponding to a temperature 
\begin{equation}
\label{eq:hawk_temp}
T_{\rm H}= \frac{\delta +1}{4\pi r_S}\,,
\end{equation}
into a large number~\cite{Dimopoulos:2001hw},   
\begin{equation}
\label{eq:av_mult}
\langle n\rangle \approx \frac{1}{2}\,\frac{M_{\rm bh}}{T_{\rm H}}\,,
\end{equation}
of hard quanta (cf. Fig.~\ref{multiplicity}), with energies approaching several hundreds of GeV in the cm system. 
We shall assume that the average multiplicity~(\ref{eq:av_mult}) refers to 
Standard Model particles~\cite{Emparan:2000rs} -- an assumption which is still under
debate~\cite{Casadio:2001wh} --  in a ``flavor"-democratic fashion. 



%\section{Contained events}
{\em 3.}
Let us consider now the rate of black hole events taking place within the fiducial volume of a neutrino telescope. 
The number of black hole events initiated by neutrino-nucleon scattering occurring per unit time $t$ 
and solid angle $\Omega$ with 
energy larger than a threshold energy $E_{\rm th}$ inside a subsurface detector volume 
$V$ is (see e.g. Ref.~\cite{Morris:1991bb}) 
%
\begin{equation}
\label{rate_contained}
\frac{{\rm d}^2 N_{\rm casc}^{\rm bh}}{{\rm d}t\,{\rm d}\Omega} 
= \frac{\rho\,V}{m_N}
\int^\infty_{E_{\rm th}} {\rm d}E_\nu\,
F_\nu (E_\nu )\,\sigma_{\nu N}^{\rm bh}(E_\nu )\,
\exp\left[ -\sigma^{\rm tot}_{\nu N}(E_\nu )\,X(\theta )/m_N\right]
\,,
\end{equation}
where 
$\rho$ is the density of the material in which the neutrino interaction occurs,  
$F_\nu = \sum_i (F_{\nu_i}+F_{\bar\nu_i})$ is the sum of the differential diffuse neutrino fluxes, 
$\sigma_{\nu N}^{\rm tot}$ denotes the total interaction cross section as a sum of the usual
Standard Model charged current cross section $\sigma_{\nu N}^{\rm cc}$ and the black hole production 
cross section~(\ref{sig_nuN_bh_pdf}),
\begin{equation}
\label{eq:ccpbh}
\sigma_{\nu N}^{\rm tot} = \sigma_{\nu N}^{\rm cc} + \sigma_{\nu N}^{\rm bh}\,,
\end{equation}
and, finally,  
$X(\theta  )$ is the column density of material between the detector and the upper atmosphere, as
a function of the zenith angle $\theta$. It can be approximated by~\cite{Morris:1991bb}
\begin{equation}
\label{eq:overburden}
X (\theta )  = \rho \biggl[
\sqrt{ (R_\oplus - D)^2 \cos^2\theta + 2DR_\oplus - D^2} -
(R_\oplus-D)\cos\theta
\biggr]\,, 
%+ X_{\rm atm}(\theta ),
\end{equation}
where $R_\oplus$ denotes the radius of the earth and $D$ is the vertical depth 
of the detector. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%FIGURE%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\begin{center}
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=570pt,bbury=608pt,%
width=12cm]{e2flux_wost}
\caption[dum]{\label{e2flux_cosmogenic}
Predictions of the cosmogenic neutrino flux, $F_\nu =\sum_i\left[ F_{\nu_i}+
F_{\bar\nu_i}\right]$. 
Long-dashed (long-dashed--dotted) line: Flux from 
Ref.~\cite{Yoshida:1993pt} for cosmological evolution parameters 
$m=2$, $z_{\rm max}=2$ ($m=4$, $z_{\rm max}=4$).
Solid (dotted) line: Flux from Ref.~\cite{Protheroe:1996ft}, 
assuming a maximum energy of $E_{\rm max}=3\cdot 10^{20(21)}$ eV for the ultrahigh energy cosmic rays.
Shaded: Theoretical upper limit of the neutrino flux from ``hidden'' %hadronic 
astrophysical sources that are
non-transparent to ultrahigh energy nucleons~\cite{Mannheim:2001wp}.
}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


As in Ref.~\cite{Ringwald:2001vk}, we shall 
exploit both conservative lower and upper limits~\cite{Waxman:1999yy,Mannheim:2001wp} on the 
presently unknown differential flux $F_\nu$ of ultrahigh energy cosmic neutrinos  (for recent reviews, see  
Ref.~\cite{Protheroe:1999ei}) entering into
event rates at neutrino telescopes such as Eq.~(\ref{rate_contained}).
% 
More or less guaranteed, and therefore comprising a reasonable lower bound, are the so-called 
cosmogenic neutrinos which are produced when ultrahigh energy cosmic protons
inelastically scatter off the cosmic microwave background radiation~\cite{Greisen:1966jv}
in processes such as $p\gamma\to \Delta\to n\pi^+$, where the produced pion subsequently 
decays~\cite{Beresinsky:1969qj}.
Recent estimates of these fluxes  can be found in 
Refs.~\cite{Yoshida:1993pt,Protheroe:1996ft,Yoshida:1997ie,Engel:2001hd}, some of which 
are shown in Fig.~\ref{e2flux_cosmogenic}. 
%
The upper limit from ``hidden'' hadronic astrophysical sources\footnote{For an early determination of such an upper limit, see 
Ref.~\cite{Berezinsky:1979pd}. Fluxes of this size are required in the context of the Z-burst 
scenario~\cite{Fargion:1999ft} for the highest
energy cosmic rays.}, i.\,e. from those sources from which only 
photons and neutrinos can escape, is much larger~\cite{Mannheim:2001wp} and  
also shown in Fig.~\ref{e2flux_cosmogenic}. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%FIGURE%%%%%%%%%%%%%%%%%%%%%%%%%% 
\begin{figure}
\begin{center}
%rev
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=585pt,bbury=608pt,%
width=13cm]{bh_par_casc_km3_PJ_d6_MD_att_yr.eps}
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=585pt,bbury=608pt,%
width=13cm]{bh_par_casc_km3_hid_d6_MD_att_yr.eps}
\caption[dum]{\label{bh_par_cascades}
Reach in the black hole production parameters for $\delta = 6$ extra dimensions, 
for contained events in an under-ice detector at a depth of 2 km and with an 1 km$^3$ 
ficudial volume ($1\ {\rm yr}=3.16\cdot 10^7$~s). 
The shaded solid lines labeled ``FE'' indicates the constraint~\cite{Ringwald:2001vk} 
arising from the non-observation of 
horizontal showers by the Fly's Eye collaboration~\cite{Baltrusaitis:1985mt}.  
The shaded dotted, $M_D = M_{\rm bh}^{\rm min}$, and shaded solid, 
$M_D=(1/5)\,M_{\rm bh}^{\rm min}$, lines give a rough indication of the 
boundary of applicability of  the semiclassical picture~\cite{Giddings:2001bu}.
{\em Top:} Exploiting the 
cosmogenic neutrino flux from Ref.~\cite{Protheroe:1996ft} (cf. Fig.~\ref{e2flux_cosmogenic}).
{\em Bottom:} Exploiting 
the upper limit on the neutrino flux from ``hidden'' hadronic astrophysical sources from
Ref.~\cite{Mannheim:2001wp} (cf. Fig.~\ref{e2flux_cosmogenic}). 
}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The projected reach in the black hole production parameter space for $\delta = 6$ extra dimensions 
for contained events in an under-ice detector corresponding to the IceCube proposal~\cite{Halzen:1999jy}
(2 km depth, fiducial volume 1 km$^3$) is shown in Fig.~\ref{bh_par_cascades}, for two of the
discussed ultrahigh energy cosmic neutrino fluxes. 
Table~\ref{tb:contained_sm} summarizes the corresponding expected Standard Model background rates 
for contained events due to charged current interactions.
Taking into account the small effective volume, $V\approx 0.001\div 0.01$ km$^3$, of the currently operating 
AMANDA and Baikal neutrino telescopes, it seems that, on the basis of contained events, these facilities
cannot compete with the already established limits from air shower observatories such as 
Fly's Eye and AGASA, and
the projected ones from Auger~\cite{Feng:2002ib,Ringwald:2001vk,Anchordoqui:2001cg}. 
Even IceCube does not seem to be really competitive, as far as the contained events are concerned: 
According to the time schedule of IceCube, the construction starts 
in 2004 and the final effective volume, $V\approx 1$ km$^3$, is reached only 
after the LHC starts operating and Auger has taken data for already a 
few years.
For the under-water/ice neutrino telescopes considered so far, the best perspective for detection of 
black hole events on the basis of the contained events might thus be assigned to 
RICE~\cite{Kravchenko:2001id}, a currently operating radio-Cherenkov telescope with an 
effective volume of about 1~km$^3$ for 10$^8$~GeV electromagnetic cascades~\cite{Seckel:2001}.

%%%%%%%%%%%%%%%%%%%table%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{table}
\begin{center}
\begin{tabular}{|c||c|c|c|c|}
\hline 
 & \multicolumn{2}{|c|}{cosmogenic neutrinos} & \multicolumn{2}{|c|}{hidden sources}\\
\hline
$E_{\rm th}$& contained & through-going & contained & through-going \\
\hline    
$10^7$~GeV & 0.47   (yr km$^3)^{-1}$ & 0.85   (yr km$^2)^{-1}$ 
                      & 116  (yr km$^3)^{-1}$ & 92  (yr km$^2)^{-1}$\\
%$10^8$~GeV & 0.38   (yr km$^3)^{-1}$ & 0.44   (yr km$^3)^{-1}$ 
%& 23  (yr km$^3)^{-1}$ & 18  (yr km$^2)^{-1}$\\
\hline
\end{tabular}
\caption{Event rates ($1\ {\rm yr}=3.16\cdot 10^7$~s) of contained reactions and 
through-going muons, respectively, for  
Standard Model like cross sections for the case of the cosmogenic neutrino flux from 
Ref.~\cite{Protheroe:1996ft}, assuming a maximum energy of $E_{\rm max}=3\cdot 10^{21}$ eV 
for the ultrahigh energy cosmic rays, and the case of the upper limit 
on the neutrino flux from ``hidden'' hadronic astrophysical sources from
Ref.~\cite{Mannheim:2001wp}, respectively. An under-ice detector at a depth of 2 km has
been assumed. For the contained events, we have considered,   
as in Eq.~(\ref{eq:ccpbh}), charged current reactions from all flavors, 
but have neglected neutral current reactions, because 
the fraction of visible energy, $\langle y\rangle\approx 0.2$, as well as 
the cross section, is lower for the latter~\cite{Gandhi:1996tf}.
For the through-going events, we have considered charged current $\nu_\mu N$
and $\bar\nu_\mu N$ interactions and an equipartion of neutrino flavors in the
neutrino flux.}
\label{tb:contained_sm}
\end{center}
\end{table}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\section{Through-going muons}
{\em 4.}
The ability to detect muons from distant neutrino reactions increases an
under-ice/water detector's effective neutrino target volume dramatically
and is the premise upon which such detectors can act as neutrino telescopes. 
Let us therefore, in addition to the rate~(\ref{rate_contained}) of fully contained
black hole events, estimate the rate of through-going muons at a neutrino telescope. 
We define through-going muon events to be events with a vertex located 
outside of the detector and which then pass through the detector.
The expected rate for such events is then proportional to the area of
the detector. The muon rate above some threshold $E_{\rm th}$ per unit area $A$ and
solid angle $\Omega$ can be estimated by 
\begin{eqnarray}
\label{eq:trate}
\lefteqn{
\frac{{\rm d}^3 N^{\rm bh}_\mu}{{\rm d}A\,{\rm d}t\,{\rm d}\Omega} = \,
\sum_{k=1}^\infty \frac{1}{m_{ N}} 
\int_{E_{\rm th}}^\infty {\rm d}E_\nu\, F_\nu ( E_\nu )\, 
\int {\rm d}n_\mu\, \frac{{\rm d}\sigma^{\rm bh}_{\nu N}}{{\rm d}n_{\mu}}\ 
%p_k (n_{\mu})
\frac{n_\mu^k}{k\,!}\,{\rm e}^{-n_\mu}
\,\times \,}
\\[1.5ex]\nonumber &&
\mbox{}\hspace{7ex}\times \int_{0}^{X(\theta )} {\rm d}X^{\prime} 
\exp\left[ -\sigma^{\rm tot}_{\nu N}(E_\nu )\,(X(\theta )-X^\prime )/m_N\right]
\left( 1- \left(1-\tilde{p}^{\rm bh}_\mu(E_\nu, \,E_{\rm th},X^{\prime})\right)^k\right)\,.
\end{eqnarray}
In words, the rate of through-going muon events from a black hole reaction is obtained by 
summing over different muon multiplicities $k$, weighted with the Poisson 
probability of mean $n_\mu$   
that $k$ muons emerge from black hole decay and that at least one can be observed in the detector.
%
The quantity $\tilde p_\mu^{\rm bh}$, which may be interpreted as the 
detection probability of a ``typical'' muon produced in a black hole
event, is defined as (cf. also Ref.~\cite{Morris:1991bb})
\begin{equation}
\label{eq:ptilde}
{\tilde p}_\mu^{\rm bh} (E_\nu , E_{\rm th}, X^\prime )
=
\int_{E_{\rm th}}^{E_\nu} {\rm d}E_\mu\,
\frac{1}{n_\mu^{\rm bh}}\frac{{\rm d}n_\mu^{\rm bh}}{{\rm d}E_\mu}\,
p_\mu (E_\mu , E_{\rm th}, X^\prime )\,
\end{equation}
with $p_\mu (E_\mu , E_{\rm th}, X^\prime )$ denoting the probability 
that a muon of definite energy $E_\mu$ traverses an amount $X^\prime$
of material with $E_\mu >E_{\rm th}$. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%FIGURE%%%%%%%%%%%%%%%%%%%%%%%%%% 
\begin{figure}
\begin{center}
%rev
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=585pt,bbury=608pt,%
width=13cm]{bh_par_through_km2_PJ_d6_MD_att_yr.eps}
\includegraphics*[bbllx=20pt,bblly=221pt,bburx=585pt,bbury=608pt,%
width=13cm]{bh_par_through_km2_hid_d6_MD_att_yr.eps}
\caption[dum]{\label{bh_par_through}
Reach in the black hole production parameters for $\delta = 6$ extra dimensions, 
for through-going muons in an under-ice detector at a depth of 2 km and with an 1 km$^2$ 
effective area ($1\ {\rm yr}=3.16\cdot 10^7$~s). 
The shaded solid lines labeled ``FE'' indicates the constraint~\cite{Ringwald:2001vk}  arising from the non-observation of 
horizontal showers by the Fly's Eye collaboration~\cite{Baltrusaitis:1985mt}.  
The shaded dotted, $M_D = M_{\rm bh}^{\rm min}$, and shaded solid, 
$M_D=(1/5)\,M_{\rm bh}^{\rm min}$, lines give a rough indication of the 
boundary of applicability of  the semiclassical picture~\cite{Giddings:2001bu}.
{\em Top:} Exploiting the 
cosmogenic neutrino flux from Ref.~\cite{Protheroe:1996ft} (cf. Fig.~\ref{e2flux_cosmogenic}).
{\em Bottom:} Exploiting 
the upper limit on the neutrino flux from ``hidden'' hadronic astrophysical sources from
Ref.~\cite{Mannheim:2001wp} (cf. Fig.~\ref{e2flux_cosmogenic}). 
}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


If one neglects the stochastic 
effects of range straggling which may become important for muon energies close
to the detection threshold,  i.e. within the approximation
of a deterministic energy-range relationship, one has the step 
function
\begin{equation}
p_\mu (E_\mu , E_{\rm th}, X^\prime ) = 
\theta (E_\mu - E_\mu^\prime (E_{\rm th},X^\prime ))\,
\end{equation}
where $E_\mu^\prime$ is the solution of 
\begin{equation}
X^\prime = \int_{E_{\rm th}}^{E_\mu^\prime}
\frac{{\rm d}E}{a (E) + b (E)\,E}\,.
\end{equation}
The ionization and radiation loss parameters,  $a (E)$ and $b (E)$, 
vary only slowly with energy and take the asymptotic values 
$a \approx 3\cdot 10^{-3}$ GeV cm$^2$/g and $b\approx 4\cdot 10^{-6}$ cm$^2$/g~\cite{Groom:2000in}.
Choosing these constant values, we get for the average muon range,
\begin{equation}
\label{eq:murange}
X^\prime \approx \frac{1}{b}\,
\ln\frac{a + b E^\prime_\mu}{ a + b E_{\rm th}}
\hspace{4ex}
\Leftrightarrow
\hspace{4ex}
E_\mu^\prime (E_{\rm th},X^\prime ) \approx 
\left[ E_{\rm th} + \left( 1 - {\rm e}^{- b\,X^\prime}\right)\,
\frac{a}{b}\right] {\rm e}^{b\,X^\prime}\,. 
\end{equation}
Since for our considerations $E_{\rm th}\sim 10^7$ GeV, we can safely neglect the contribution proportional to $a$ in the following. 

For the laboratory energy distribution of the muons produced in a 
black hole event, entering the probability~(\ref{eq:ptilde}), we have taken for 
simplicity\footnote{Alternatively, one may use a 
boosted thermal distribution, $\sim \int_{E_\mu/2\gamma}^\infty 
{\rm d}E (E/T_{\rm H})^2/(e^{E/T_{\rm H}}+1)$,  
which, however, needs an additional numerical integration.}
(cf. also Ref.~\cite{Morris:1991bb})
\begin{equation}
\label{eq:cmtolab}
\frac{1}{n_\mu^{\rm bh}}\frac{{\rm d}n_\mu^{\rm bh}}{{\rm d}E_\mu}
=\frac{\langle n\rangle}{2E_\nu}\,
\theta\left(\frac{2 E_\nu}{\langle n\rangle}-E_\mu \right)\,,
\end{equation}
with $\langle n\rangle$ being the average number~(\ref{eq:av_mult}) of particles produced in black hole decay 
(cf. Fig.~\ref{multiplicity}).
The distribution~(\ref{eq:cmtolab}) is obtained by boosting a delta-function energy spectrum, 
\begin{equation}
\label{eq:spec}
\frac{{\rm d}\sigma_{\nu N}^{\rm bh}}{{\rm d}E^{\rm cm}_\mu}
\propto \delta ( E_\mu^{\rm cm} - M_{\rm bh}/\langle n\rangle )\,,
\end{equation}
in the cm frame of the black hole, with isotropic direction of the decay particles, 
with a Lorentz factor $\gamma = E_\nu/M_{\rm bh}$ to the laboratory frame.  

Finally, we have exploited the ansatz (similar to Eq.(\ref{eq:spec})), 
\begin{equation}
\frac{d\sigma^{\rm bh}_{\nu N}}{{\rm d}n_\mu} 
= \sigma^{\rm bh}_{\nu N}\ \delta \left( n_\mu-\langle n_\mu \rangle\right)\,,
\hspace{6ex}
{\rm with\ }\ \langle n_\mu\rangle \approx \langle n\rangle /30 \,,
\end{equation}
to perform also the integration over the multiplicity $n_\mu$ in our rate
estimate~(\ref{eq:trate}).  The above relation for $\langle n_\mu \rangle$ is based on
the assumption that a black hole does not discriminate between Standard Model  particles in its 
evaporation spectrum~\cite{Giddings:2001bu}. 


The projected reach in the black hole production parameter space for $\delta = 6$ extra dimensions 
for through-going muons in an under-ice detector at a depth of 2 km and with an
effective area of 1 km$^2$ is shown in Fig.~\ref{bh_par_through}, for two of the
discussed ultrahigh energy cosmic neutrino fluxes. 
The corresponding expected Standard Model background rates 
for through-going muons due to charged current interactions is summarized in 
Table~\ref{tb:contained_sm} (see also Ref.~\cite{Alvarez-Muniz:2001es} for related
calculations).
%
An additional source of background which 
is not related to an astrophysical neutrino flux is the muon flux arising 
from charm production by cosmic rays interacting in the atmosphere.
Using two available models~\cite{Bugaev:1998bi,Gelmini:2000ve} we have 
estimated a background of 
$0.33\div 1.5$ through-going muon events from atmospheric charm production per year.
To further reduce this background, one may introduce a cut in the zenith
angle, e.g. $\theta > 60^{\,o}$.

%\section{Discussions} 
\noindent
{\em 5.} 
We have considered the sensitivity of neutrino telescopes to black hole production in the
context of extra dimension scenarios with TeV-scale gravity and compared it to
the one of air shower facilities and the LHC. 
For concreteness, we have assumed a detector with 1~km$^2$ area and 1~km$^3$ volume, 
respectively, located at a vertical depth of 2~km in ice.

We have seen that the effective volumina for contained events, $V\approx
0.001\div 0.01$~km$^{3}$, available from the currently operating AMANDA and
Baikal detectors  are too small 
to yield a large enough rate to be competitive with the already existing
constraints on black hole production from Fly's Eye and AGASA and with the 
projected ones from
Auger~\cite{Feng:2002ib,Ringwald:2001vk,Anchordoqui:2001cg}.  
The planned IceCube detector with a fiducial volume of about 1 km$^3$ 
would have an interesting perspective.
However, with the start of construction of IceCube in 2004 and full 
completion only after the physics start at the LHC, the discovery potential for
black hole production is strongly reduced.
Instead, the radio-Cherenkov telescope RICE~\cite{Kravchenko:2001id}, with an effective volume of about 
0.1~(1)~km$^3$ for 10$^{7(8)}$~GeV electromagnetic cascades~\cite{Seckel:2001}, could set sensible constraints using 
already available data.  

The prospects for a search for black hole signatures using through-going muons are also good. 
For a detection area of 1 km$^2$ we have shown that in case the neutrino flux 
is just on the level of the 
cosmogenic one, only a few ($\lwig 1$) contained or through-going background events from Standard Model charged
current interactions are expected per year. The background of prompt
muons from atmospheric charm was estimated to be of similar 
magnitude and can be further reduced by application of a cut on the 
zenith angle.
It was shown in Ref.~\cite{Hundertmark:2001}, that AMANDA achieves 
an effective area of about $0.3$ km$^2$ for down-going muons above an energy 
 of $10^{7}$~GeV and that suppression of the atmospheric muon background is possible.
With 5 years of AMANDA data available and an assumed duty-factor of $2/3$ 
this leads to an exposure of about 1 km$^2\cdot$yr. 
Therefore, strong constraints on the black hole production 
parameters are expected from AMANDA if the currently available data
show no through-going muons\footnote{It is noteworthy that in this case 
AMANDA can put also strong constraints on sphaleron production parameters~\cite{Morris:1991bb}, 
since the phenomenology of black hole and sphaleron production are quite similar.} above $10^7$~GeV. In  this case, 
the sensitivity of  AMANDA -- approximately the 3 events/yr contour in Fig.~\ref{bh_par_through} (top) for a 
90 \% confidence level limit --  
compares favorably with the one of Fly's Eye, AGASA, and Auger~\cite{Feng:2002ib,Ringwald:2001vk,Anchordoqui:2001cg}.
And even at the kinematic limit of the LHC, M$_{\rm bh}^{\rm min}\approx 10$~TeV (cf. Ref.~\cite{Ringwald:2001vk}),
one is still left with a sensitivity on $M_D$ of about $1.3$~TeV. 

In the case that the ultrahigh energy cosmic neutrino flux considerably exceeds the cosmogenic one, 
discrimination between Standard Model background and black hole events
becomes crucial. This might not be achieved by observation of through-going 
events alone, but with the additional information of contained events.
For example, the spectral information from contained and through-going events 
should be remarkably different, as the through-going muons from black hole decays 
are produced with considerably lower energy. 
Another discriminator may be based on the different hadronic components 
of Standard Model processes and black hole events. In a typical Standard Model charged 
current interaction one has an average hadronic energy fraction of 
$\langle y\rangle \approx 0.2$, whereas 
black hole events have typically a large hadronic component ($5:1$)~\cite{Giddings:2001bu} 
corresponding to a $\langle y\rangle \approx 0.8$. 
In particular radio-Cherenkov detectors such as RICE have the potential to 
discriminate between electromagnetic and hadronic showers by measuring
the characteristically different angular radiation patterns~\cite{Alvarez-Muniz:2000xx}. 

To summarize, we have shown that large scale under-water/ice neutrino detectors
have a substantial potential to discover or rule out black hole production which
is comparable to that of air shower facilities. In the optimistic case that a 
neutrino flux significantly higher than the cosmogenic one is realized in nature, 
one even has discovery potential for black holes at neutrino telescopes beyond the 
reach of LHC.

\section*{Acknowledgements}
We thank C. Spiering 
for fruitful discussions.

\begin{thebibliography}{99}

%Black hole production in 4d

%rev
%\cite{'tHooft:1987rb}
\bibitem{'tHooft:1987rb}
G.~'t Hooft,
%``Graviton Dominance In Ultrahigh-Energy Scattering,''
Phys.\ Lett.\ B {198} (1987) 61; 
%%CITATION = PHLTA,B198,61;%%
%
%\cite{'tHooft:1988wk}
%\bibitem{'tHooft:1988wk}
%G.~'t Hooft,
%``On The Factorization Of Universal Poles In A Theory Of Gravitating Point Particles,''
Nucl.\ Phys.\ B {304} (1988) 867;
%%CITATION = NUPHA,B304,867;%%
%
%\cite{Amati:1987wq}
%\bibitem{Amati:1987wq}
D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Superstring Collisions At Planckian Energies,''
Phys.\ Lett.\ B {197} (1987) 81; 
%
%%CITATION = PHLTA,B197,81;%%
%\cite{Amati:1988uf}
%\bibitem{Amati:1988uf}
%D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Classical And Quantum Gravity Effects From Planckian Energy Superstring Collisions,''
Int.\ J.\ Mod.\ Phys.\ A {3} (1988) 1615; 
%%CITATION = IMPAE,A3,1615;%%
%
%\cite{Amati:1989tn}
%\bibitem{Amati:1989tn}
%D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Can Space-Time Be Probed Below The String Size?,''
Phys.\ Lett.\ B {216} (1989) 41; 
%%CITATION = PHLTA,B216,41;%%
%
%\cite{Amati:1990xe}
%\bibitem{Amati:1990xe}
%D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Higher Order Gravitational Deflection And Soft Bremsstrahlung In Planckian Energy Superstring Collisions,''
Nucl.\ Phys.\ B {347} (1990) 550; 
%%CITATION = NUPHA,B347,550;%%
%
%\cite{Amati:1992zb}
%\bibitem{Amati:1992zb}
%D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Planckian scattering beyond the semiclassical approximation,''
Phys.\ Lett.\ B {289} (1992) 87;
%%CITATION = PHLTA,B289,87;%%
%
%\cite{Amati:1993tb}
%\bibitem{Amati:1993tb}
%D.~Amati, M.~Ciafaloni and G.~Veneziano,
%``Effective action and all order gravitational eikonal at Planckian energies,''
Nucl.\ Phys.\ B {403} (1993) 707;
%%CITATION = NUPHA,B403,707;%%
%
%\cite{Fabbrichesi:1994kd}
%\bibitem{Fabbrichesi:1994kd}
M.~Fabbrichesi, R.~Pettorino, G.~Veneziano and G.~A.~Vilkovisky,
%``Planckian energy scattering and surface terms in the gravitational action,''
Nucl.\ Phys.\ B {419} (1994) 147;
%%CITATION = NUPHA,B419,147;%%
%
%\cite{Aref'eva:1995qs}
%\bibitem{Aref'eva:1995qs}
I.~Y.~Aref'eva, K.~S.~Viswanathan and I.~V.~Volovich,
%``Planckian energy scattering, colliding plane gravitational waves and black hole creation,''
Nucl.\ Phys.\ B  {452} (1995) 346
[Erratum-ibid.\ B  {462} (1995) 613].
%.
%%CITATION = ;%%

%Large extra dimensions and low scale QG

%\cite{Arkani-Hamed:1998rs}
\bibitem{Arkani-Hamed:1998rs}
N.~Arkani-Hamed, S.~Dimopoulos and G.~R.~Dvali,
%``The hierarchy problem and new dimensions at a millimeter,''
Phys.\ Lett.\ B { 429} (1998) 263;
%.
%%CITATION = ;%%
%
%\cite{Antoniadis:1998ig}
%\bibitem{Antoniadis:1998ig}
I.~Antoniadis, N.~Arkani-Hamed, S.~Dimopoulos and G.~R.~Dvali,
%``New dimensions at a millimeter to a Fermi and superstrings at a TeV,''
Phys.\ Lett.\ B { 436} (1998) 257;
%.
%%CITATION = ;%%
%
%\cite{Arkani-Hamed:1999nn}
%\bibitem{Arkani-Hamed:1999nn}
N.~Arkani-Hamed, S.~Dimopoulos and G.~R.~Dvali,
%``Phenomenology, astrophysics and cosmology of theories with  sub-millimeter dimensions and TeV scale quantum gravity,''
Phys.\ Rev.\ D { 59} (1999) 086004.
%.
%%CITATION = ;%%

%noncompact extra dimensions and TeV scale gravity

%\cite{Randall:1999ee}
\bibitem{Randall:1999ee}
L.~J.~Randall and R.~Sundrum,
%``A large mass hierarchy from a small extra dimension,''
Phys.\ Rev.\ Lett.\  {83} (1999) 3370.
%.
%%CITATION = ;%%


%Black hole production

%\cite{Argyres:1998qn}
\bibitem{Argyres:1998qn}
P.~C.~Argyres, S.~Dimopoulos and J.~March-Russell,
%``Black holes and sub-millimeter dimensions,''
Phys.\ Lett.\ B {441} (1998) 96; 
%.
%%CITATION = ;%%
%
%\cite{Banks:1999gd}
%\bibitem{Banks:1999gd}
T.~Banks and W.~Fischler,
%``A model for high energy scattering in quantum gravity,''
; 
%%CITATION = ;%%
%
%\cite{Aref'eva:1999bm}
%\bibitem{Aref'eva:1999bm}
I.~Y.~Aref'eva,
%``High energy scattering in the brane-world and black hole production,''
%arXiv:
.
%%CITATION = ;%%

%\cite{Emparan:2000rs}
\bibitem{Emparan:2000rs}
R.~Emparan, G.~T.~Horowitz and R.~C.~Myers,
%``Black holes radiate mainly on the brane,''
Phys.\ Rev.\ Lett.\  {85} (2000) 499.
%.
%%CITATION = ;%%

%\cite{Giddings:2000ay}
\bibitem{Giddings:2000ay}
S.~B.~Giddings and E.~Katz,
%``Effective theories and black hole production in warped  compactifications,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Emparan:2001ce}
%\bibitem{Emparan:2001ce}
R.~Emparan,
%``Exact gravitational shockwaves and Planckian scattering on branes,''
Phys.\ Rev.\ D {64} (2001) 024025.
%.
%%CITATION = ;%%

%LHC

%\cite{Evans:2001mn}
\bibitem{Evans:2001mn}
L.~R.~Evans,
%``The Large Hadron Collider: Present status and prospects,''
CERN-OPEN-2001-027, presented at: 
18th International Conference on High Energy Accelerators, Tsukuba, Japan, 26 - 30 Mar 2001.

%Black hole production at the LHC

%\cite{Giddings:2001bu}
\bibitem{Giddings:2001bu}
S.~B.~Giddings and S.~Thomas,
%``High energy colliders as black hole factories: The end of short  distance physics,''
.
%%CITATION = ;%%


%\cite{Dimopoulos:2001hw}
\bibitem{Dimopoulos:2001hw}
S.~Dimopoulos and G.~Landsberg,
%``Black holes at the LHC,''
Phys.\ Rev.\ Lett.\  {87} (2001) 161602.
%.
%%CITATION = ;%%

%Black hole production at LHC continued


%\cite{Dimopoulos:2001qe}
\bibitem{Dimopoulos:2001qe}
S.~Dimopoulos and R.~Emparan,
%``String balls at the LHC and beyond,''
.
%%CITATION = ;%%

%\cite{Hossenfelder:2001dn}
\bibitem{Hossenfelder:2001dn}
S.~Hossenfelder, S.~Hofmann, M.~Bleicher and H.~St\"ocker,
%``Quasi-Stable Black Holes at LHC,''
;
%%CITATION = ;%%
%
%\cite{Giddings:2001ih}
%\bibitem{Giddings:2001ih}
S.~B.~Giddings,
%``Black hole production in TeV-scale gravity, and the future of high  energy physics,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Cheung:2001ue}
%\bibitem{Cheung:2001ue}
K.~Cheung,
%``Black hole production and large extra dimensions,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Park:2001xc}
%\bibitem{Park:2001xc}
S.~C.~Park and H.~S.~Song,
%``Production of spinning black holes at colliders,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Rizzo:2001dk}
%\bibitem{Rizzo:2001dk}
T.~G.~Rizzo,
%``Black hole production rates at the LHC: Still large,''
%in {\it Proc. of the APS/DPF/DPB Summer Study on the Future of Particle Physics (Snowmass 2001) } ed. R.~Davidson and C.~Quigg,
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Landsberg:2001sj}
%\bibitem{Landsberg:2001sj}
G.~Landsberg,
%``Discovering new physics in the decays of black holes,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Giudice:2001ce}
%\bibitem{Giudice:2001ce}
G.~F.~Giudice, R.~Rattazzi and J.~D.~Wells,
%``Transplanckian collisions at the LHC and beyond,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Bleicher:2001kh}
%\bibitem{Bleicher:2001kh}
M.~Bleicher, S.~Hofmann, S.~Hossenfelder and H.~St\"ocker,
%``Black hole production in large extra dimensions at the Tevatron: New  limit on the fundamental scale of gravity,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Eardley:2002re}
%\bibitem{Eardley:2002re}
D.~M.~Eardley and S.~B.~Giddings,
%``Classical Black Hole Production in High-Energy Collisions,''
%arXiv:
.
%%CITATION = ;%%


%Black holes in nuN scattering and cosmic rays

%\cite{Feng:2002ib}
\bibitem{Feng:2002ib}
J.~L.~Feng and A.~D.~Shapere,
%``Black hole production by cosmic rays,''
Phys.\ Rev.\ Lett.\  {88} (2002) 021303.
%.
%%CITATION = ;%%

%\cite{Emparan:2001kf}
\bibitem{Emparan:2001kf}
R.~Emparan, M.~Masip and R.~Rattazzi,
%``Cosmic rays as probes of large extra dimensions and TeV gravity,''
%arXiv:
.
%%CITATION = ;%%

%\cite{Anchordoqui:2001ei}
\bibitem{Anchordoqui:2001ei}
L.~Anchordoqui and H.~Goldberg,
%``Experimental signature for black hole production in neutrino air  showers,''
%arXiv:
.
%%CITATION = ;%%

%Auger Observatory

%\cite{Zavrtanik:2000zi}
\bibitem{Zavrtanik:2000zi}
D.~Zavrtanik  [AUGER Collaboration],
%``The Pierre Auger Observatory,''
Nucl.\ Phys.\ Proc.\ Suppl.\  {85} (2000) 324.
%%CITATION = NUPHZ,85,324;%%

%\cite{Ringwald:2001vk}
\bibitem{Ringwald:2001vk}
A.~Ringwald and H.~Tu,
%``Collider versus cosmic ray sensitivity to black hole production,''
Phys.\ Lett.\ B {525} (2002) 135.
%.
%%CITATION = ;%%

%Constraints on black holes from AGASA

%\cite{Anchordoqui:2001cg}
\bibitem{Anchordoqui:2001cg}
L.~A.~Anchordoqui, J.~L.~Feng, H.~Goldberg and A.~D.~Shapere,
%``Black Holes from Cosmic Rays: Probes of Extra Dimensions and New Limits on TeV-Scale Gravity,''
%arXiv:
.
%%CITATION = ;%%

%Fly's Eye limits

%\cite{Baltrusaitis:1985mt}
\bibitem{Baltrusaitis:1985mt}
R.~M.~Baltrusaitis {\it et al.},
%``Limits On Deeply Penetrating Particles In The > 10**17-Ev Cosmic Ray Flux,''
Phys.\ Rev.\ D {31} (1985) 2192.
%%CITATION = PHRVA,D31,2192;%%

%AGASA limits on horizontal air shower

%\cite{Yoshida:2001}
\bibitem{Yoshida:2001}
S.~Yoshida {\em et al.} [AGASA Collaboration], 
in: Proc. 27th International Cosmic Ray Conference, Hamburg, Germany, 2001,
Vol. 3, pp. 1142-1145.


%Fly's Eye 

%\cite{Baltrusaitis:1985mx}
\bibitem{Baltrusaitis:1985mx}
R.~M.~Baltrusaitis {\it et al.},
%``The Utah Fly's Eye Detector,''
Nucl.\ Instrum.\ Meth.\ A {240} (1985) 410.
%%CITATION = NUIMA,A240,410;%%

%AGASA Collaboration

%\cite{Chiba:1992nf}
\bibitem{Chiba:1992nf}
N.~Chiba {\it et al.},
%``Akeno giant air shower array (AGASA) covering 100-km**2 area,''
Nucl.\ Instrum.\ Meth.\ A {311} (1992) 338.
%%CITATION = NUIMA,A311,338;%%

%Summary of constraints on large extra dimensions

%\cite{Peskin:2000ti}
\bibitem{Peskin:2000ti}
M.~E.~Peskin,
%``Theoretical summary,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Abe:2001nq}
%\bibitem{Abe:2001nq}
T.~Abe {\it et al.}  [American Linear Collider Working Group Collaboration],
%``Linear collider physics resource book for Snowmass 2001. 3: Studies of  exotic and standard model physics,''
%arXiv:
;
%%CITATION = ;%%
%
%rev
%\cite{Pagliarone:2001ff}
%\bibitem{Pagliarone:2001ff}
C.~Pagliarone,
%``Extra dimensions and black hole production,''
%arXiv:
.
%%CITATION = ;%%


%ICECUBE

%\cite{Halzen:1999jy}
\bibitem{Halzen:1999jy}
F.~Halzen {\it et al.}  [AMANDA Collaboration],
%``From the first neutrino telescope, the antarctic muon and neutrino  detector array AMANDA, to the IceCube observatory,''
in: B.L. Dingus, D.B. Kieda, M.H. Salamon (Eds.),  
26th International Cosmic Ray Conference (ICRC 99), Salt Lake City, UT, 17-25 Aug 1999,
Melville, AIP, 2000, pp. 428-431.

%ANTARES

%\cite{Montaruli:2002gh}
\bibitem{Montaruli:2002gh}
T.~Montaruli  [ANTARES Collaboration],
%``Status report of the ANTARES project,''
%arXiv:
.
%%CITATION = ;%%

%Baikal

%\cite{Domogatsky:2001br}
\bibitem{Domogatsky:2001br}
V.~Balkanov {\em et al.} [Baikal Collaboration],
%``Baikal experiment: Status report,''
%arXiv:
.
%%CITATION = ;%%

%NESTOR

%\cite{Grieder:2001kr}
\bibitem{Grieder:2001kr}
P.~K.~Grieder  [NESTOR Collaboration],
%``NESTOR neutrino telescope status report,''
Nucl.\ Phys.\ Proc.\ Suppl.\  {97} (2001) 105.
%%CITATION = NUPHZ,97,105;%%

%RICE

%\cite{Kravchenko:2001id}
\bibitem{Kravchenko:2001id}
I.~Kravchenko {\it et al.}  [RICE Collaboration],
%``Performance and simulation of the RICE detector,''
%arXiv:
.
%%CITATION = ;%%

%black holes in ICECUBE

%\cite{Uehara:2001yk}
\bibitem{Uehara:2001yk}
Y.~Uehara,
%``Production and Detection of Black Holes at Neutrino Array,''
%arXiv:
.
%%CITATION = ;%%

%Black hole Schwarzschild radius

%\cite{Myers:1986un}
\bibitem{Myers:1986un}
R.~C.~Myers and M.~J.~Perry,
%``Black Holes In Higher Dimensional Space-Times,''
Annals Phys.\  {172} (1986) 304.
%%CITATION = APNYA,172,304;%%

%\cite{Giudice:1999ck}
\bibitem{Giudice:1999ck}
G.~F.~Giudice, R.~Rattazzi and J.~D.~Wells,
%``Quantum gravity and extra dimensions at high-energy colliders,''
Nucl.\ Phys.\ B {544} (1999) 3
and erratum to appear  v2].
%%CITATION = ;%%

%pdflib

%\cite{Plothow-Besch:1993qj}
\bibitem{Plothow-Besch:1993qj}
H.~Plothow-Besch,
%``PDFLIB: A Library of all available parton density functions of the nucleon, 
%the pion and the photon and the corresponding alpha-s calculations,''
Comput.\ Phys.\ Commun.\  {75} (1993) 396;
%%CITATION = CPHCB,75,396;%%
%
%\cite{Plothow-Besch:1995ci}
%\bibitem{Plothow-Besch:1995ci}
%H.~Plothow-Besch,
%``The Parton distribution function library,''
Int.\ J.\ Mod.\ Phys.\ A {10} (1995) 2901; 
%%CITATION = IMPAE,A10,2901;%%
%
%\cite{Plothow-Besch:1995ci}
%\bibitem{Plothow-Besch:2001}
%H.~Plothow-Besch, 
http://consult.cern.ch/writeups/pdflib/main.ps

%cteq3 pdf

%\cite{Lai:1995bb}
\bibitem{Lai:1995bb}
H.~L.~Lai {\it et al.},
%``Global QCD analysis and the CTEQ parton distributions,''
Phys.\ Rev.\ D {51} (1995) 4763.
%.
%%CITATION = ;%%


%\cite{Voloshin:2001vs}
\bibitem{Voloshin:2001vs}
M.~B.~Voloshin,
%``Semiclassical suppression of black hole production in particle  collisions,''
Phys.\ Lett.\ B {518} (2001) 137;
%.
%%CITATION = ;%%
%
%\cite{Voloshin:2002fe}
%\bibitem{Voloshin:2002fe}
%M.~B.~Voloshin,
%``More remarks on suppression of large black hole production in particle  collisions,''
Phys.\ Lett.\ B {524} (2002) 376.
%.
%%CITATION = ;%%

%production of semiclassical field configurations

%\cite{Klinkhamer:1984di}
\bibitem{Klinkhamer:1984di}
F.~R.~Klinkhamer and N.~S.~Manton,
%``A Saddle Point Solution In The Weinberg-Salam Theory,''
Phys.\ Rev.\ D {30} (1984) 2212.
%%CITATION = PHRVA,D30,2212;%%

%\cite{Ringwald:1990ee}
\bibitem{Ringwald:1990ee}
A.~Ringwald,
%``High-Energy Breakdown Of Perturbation Theory In The Electroweak Instanton Sector,''
Nucl.\ Phys.\ B {330} (1990) 1;
%%CITATION = NUPHA,B330,1;%%
%
%\cite{Espinosa:1990qn}
%\bibitem{Espinosa:1990qn}
O.~Espinosa,
%``High-Energy Behavior Of Baryon And Lepton Number Violating Scattering Amplitudes And 
%Breakdown Of Unitarity In The Standard Model,''
Nucl.\ Phys.\ B {343} (1990) 310;
%%CITATION = NUPHA,B343,310;%%
%
%\cite{McLerran:1990ab}
%\bibitem{McLerran:1990ab}
L.~D.~McLerran, A.~I.~Vainshtein and M.~B.~Voloshin,
%``Electroweak Interactions Become Strong At Energy Above Approximately 10-Tev,''
Phys.\ Rev.\ D {42} (1990) 171;
%%CITATION = PHRVA,D42,171;%%
%
%\cite{Ringwald:1991bh}
%\bibitem{Ringwald:1991bh}
A.~Ringwald and C.~Wetterich,
%``How strong are weak interactions in the multi - TeV range?,''
Nucl.\ Phys.\ B {353} (1991) 303;
%%CITATION = NUPHA,B353,303;%%
%
%\cite{Ringwald:1991qz}
%\bibitem{Ringwald:1991qz}
A.~Ringwald, F.~Schrempp and C.~Wetterich,
%``Phenomenology of geometrical flavor interactions at TeV energies,''
Nucl.\ Phys.\ B {365} (1991) 3.
%%CITATION = NUPHA,B365,3;%%

%Collider Phenomenology of sphaleron production and decsy

%\cite{Morris:1991bb}
\bibitem{Morris:1991bb}
D.~A.~Morris and R.~Rosenfeld,
%``Muon bundles from cosmic ray multi - W phenomena,''
Phys.\ Rev.\ D {44} (1991) 3530; 
%%CITATION = PHRVA,D44,3530;%%
%
%constraints for new physics from Fly's Eye + cosmogenic flux
%
%\cite{Morris:1994wg}
%\bibitem{Morris:1994wg}
D.~A.~Morris and A.~Ringwald,
%``Cosmic ray signatures of multi - W processes,''
Astropart.\ Phys.\  {2} (1994) 43.
%.
%%CITATION = ;%%

%\cite{Gibbs:1995cw}
\bibitem{Gibbs:1995cw}
M.~J.~Gibbs, A.~Ringwald, B.~R.~Webber and J.~T.~Zadrozny,
%``Monte Carlo simulation of baryon and lepton number violating processes at high-energies,''
Z.\ Phys.\ C {66} (1995) 285.
%.
%%CITATION = ;%%

%\cite{Cornwall:1990hh}
\bibitem{Cornwall:1990hh}
J.~M.~Cornwall,
%``On The High-Energy Behavior Of Weakly Coupled Gauge Theories,''
Phys.\ Lett.\ B {243} (1990) 271;
%%CITATION = PHLTA,B243,271;%%
%
%\cite{Goldberg:1990qk}
%\bibitem{Goldberg:1990qk}
H.~Goldberg,
%``Breakdown Of Perturbation Theory At Tree Level In Theories With Scalars,''
Phys.\ Lett.\ B {246} (1990) 445.
%%CITATION = PHLTA,B246,445;%%

%\cite{Mattis:1992bj}
\bibitem{Mattis:1992bj}
M.~P.~Mattis,
%``The Riddle of high-energy baryon number violation,''
Phys.\ Rept.\  {214} (1992) 159;
%%CITATION = PRPLC,214,159;%%
%
%\cite{Tinyakov:1993dr}
%\bibitem{Tinyakov:1993dr}
P.~G.~Tinyakov,
%``Instanton like transitions in high-energy collisions,''
Int.\ J.\ Mod.\ Phys.\ A {8} (1993) 1823;
%%CITATION = IMPAE,A8,1823;%%
%
%\cite{Voloshin:1994yp}
%\bibitem{Voloshin:1994yp}
M.~B.~Voloshin, 
%``Nonperturbative methods,''
;
%%CITATION = ;%%
%
%\cite{Rubakov:1996vz}
%\bibitem{Rubakov:1996vz}
V.~A.~Rubakov and M.~E.~Shaposhnikov,
%``Electroweak baryon number non-conservation in the early universe and in  high-energy collisions,''
Usp.\ Fiz.\ Nauk {166} (1996) 493
[Phys.\ Usp.\  {39} (1996) 461].
%.
%%CITATION = ;%%

%\cite{Kowalski:inprep}
\bibitem{Kowalski:inprep}
M.~Kowalski, A.~Ringwald and H.~Tu, in preparation.

%Hawking radiation

%\cite{Hawking:1974rv}
\bibitem{Hawking:1974rv}
S.~W.~Hawking,
%``Black Hole Explosions,''
Nature {248} (1974) 30; 
%%CITATION = NATUA,248,30;%%
%
%\cite{Hawking:1975sw}
%\bibitem{Hawking:1975sw}
%S.~W.~Hawking,
%``Particle Creation By Black Holes,''
Commun.\ Math.\ Phys.\  {43} (1975) 199.
%%CITATION = CMPHA,43,199;%%

%\cite{Casadio:2001wh}
\bibitem{Casadio:2001wh}
R.~Casadio and B.~Harms,
%``Can black holes and naked singularities be detected in accelerators?,''
%arXiv:
.
%%CITATION = ;%%

%UHECNU flux: upper bound

%\cite{Waxman:1999yy}
\bibitem{Waxman:1999yy}
E.~Waxman and J.~N.~Bahcall,
%``High energy neutrinos from astrophysical sources: An upper bound,''
Phys.\ Rev.\ D {59} (1999) 023002; 
%. 
%%CITATION = ;%%
%
%\cite{Bahcall:2001yr}
%\bibitem{Bahcall:2001yr}
J.~N.~Bahcall and E.~Waxman,
%``High energy astrophysical neutrinos: The upper bound is robust,''
Phys.\ Rev.\ D {64} (2001) 023002.
%.
%%CITATION = ;%%

%\cite{Mannheim:2001wp}
\bibitem{Mannheim:2001wp}
K.~Mannheim, R.~J.~Protheroe and J.~P.~Rachen,
%``On the cosmic ray bound for models of extragalactic neutrino  production,''
Phys.\ Rev.\ D {63} (2001) 023003.
%.
%%CITATION = ;%%


%uhecnu fluxes: reviews

%\cite{Protheroe:1999ei}
\bibitem{Protheroe:1999ei}
R.~J.~Protheroe,
%``High energy neutrino astrophysics,''
Nucl.\ Phys.\ Proc.\ Suppl.\  {77} (1999) 465;
%%CITATION = NUPHZ,77,465;%%
%
%\cite{Gandhi:2000kq}
%\bibitem{Gandhi:2000kq}
R.~Gandhi,
%``Ultra-high energy neutrinos: A review of theoretical and  phenomenological issues,''
Nucl.\ Phys.\ Proc.\ Suppl.\  {91} (2000) 453;
%.
%%CITATION = ;%%
%
%\cite{Learned:2000sw}
%\bibitem{Learned:2000sw}
J.~G.~Learned and K.~Mannheim,
%``High-energy neutrino astrophysics,''
Ann.\ Rev.\ Nucl.\ Part.\ Sci.\  {50} (2000) 679.
%%CITATION = ARNUA,50,679;%%

%gzk effect

%\cite{Greisen:1966jv}
\bibitem{Greisen:1966jv}
K.~Greisen,
%``End To The Cosmic Ray Spectrum?,''
Phys.\ Rev.\ Lett.\  {16} (1966) 748;
%%CITATION = PRLTA,16,748;%%
%
%\cite{Zatsepin:1966jv}
%\bibitem{Zatsepin:1966jv}
G.~T.~Zatsepin and V.~A.~Kuzmin,
%``Upper Limit Of The Spectrum Of Cosmic Rays,''
JETP Lett.\  {4} (1966) 78
[Pisma Zh.\ Eksp.\ Teor.\ Fiz.\  {4} (1966) 114].
%%CITATION = JTPLA,4,78;%%



%cosmogenic neutrino fluxes

%history

%\cite{Beresinsky:1969qj}
\bibitem{Beresinsky:1969qj}
V.~S.~Berezinsky and G.~T.~Zatsepin,
%``Cosmic Rays At Ultrahigh-Energies (Neutrino?),''
Phys.\ Lett.\ B {28} (1969) 423;
%%CITATION = PHLTA,B28,423;%%
%
%\cite{Beresinsky:1970}
%\bibitem{Beresinsky:1970}
V.~S.~Berezinsky and G.~T.~Zatsepin,
Sov.\ J.\ Nucl.\ Phys.\  {11} (1970) 111
[Yad.\ Fiz.\ {11} (1970) 200].
%%CITATION = SJNCA,11,111;%%

%recent


%\cite{Yoshida:1993pt}
\bibitem{Yoshida:1993pt}
S.~Yoshida and M.~Teshima,
%``Energy spectrum of ultrahigh-energy cosmic rays with extragalactic origin,''
Prog.\ Theor.\ Phys.\  {89} (1993) 833.
%%CITATION = PTPKA,89,833;%%

%\cite{Protheroe:1996ft}
\bibitem{Protheroe:1996ft}
R.~J.~Protheroe and P.~A.~Johnson,
%``Propagation of ultrahigh-energy protons over cosmological distances and implications for topological defect models,''
Astropart.\ Phys.\  {4} (1996) 253
[Erratum-ibid.\  {5} (1996) 215].
%.
%%CITATION = ;%%

%\cite{Yoshida:1997ie}
\bibitem{Yoshida:1997ie}
S.~Yoshida, H.~Dai, C.~C.~Jui and P.~Sommers,
%``Extremely high energy neutrinos and their detection,''
Astrophys.\ J.\  {479} (1997) 547.
%.
%%CITATION = ;%%

%\cite{Engel:2001hd}
\bibitem{Engel:2001hd}
R.~Engel and T.~Stanev,
%``Neutrinos from propagation of ultra-high energy protons,''
Phys.\ Rev.\ D {64} (2001) 093010.
%.
%%CITATION = ;%%




%\cite{Berezinsky:1979pd}
\bibitem{Berezinsky:1979pd}
V.~S.~Berezinsky,
%``High-Energy Neutrino Astronomy Versus Gamma Astronomy. (Talk),''
in: J.~G.~Learned (Ed.), DUMAND Summer Workshops,  
Khabarovsk and Lake Baikal, 22-31 Aug 1979, Hawaii DUMAND Center, 
University of Hawaii, 1980, pp. 245-261.

%Z-burst scenario

%\cite{Fargion:1999ft}
\bibitem{Fargion:1999ft}
D.~Fargion, B.~Mele and A.~Salis,
%``Ultrahigh energy neutrino scattering onto relic light neutrinos in  galactic 
%halo as a possible source of highest energy extragalactic  cosmic rays,''
Astrophys.\ J.\  {517} (1999) 725; 
%.
%%CITATION = ;%%
%
%\cite{Weiler:1999sh}
%\bibitem{Weiler:1999sh}
T.~J.~Weiler,
%``Cosmic ray neutrino annihilation on relic neutrinos revisited:  
%A mechanism for generating air showers above the  Greisen-Zatsepin-Kuzmin cut-off,''
Astropart.\ Phys.\  {11} (1999) 303; 
%.
%%CITATION = ;%%
%
%\cite{Yoshida:1998it}
%\bibitem{Yoshida:1998it}
S.~Yoshida, G.~Sigl and S.~J.~Lee,
%``Extremely high energy neutrinos, neutrino hot dark matter, and the  highest energy cosmic rays,''
Phys.\ Rev.\ Lett.\  {81} (1998) 5505; 
%.
%%CITATION = ;%%
%
%\cite{Fodor:2001qy}
%\bibitem{Fodor:2001qy}
Z.~Fodor, S.~D.~Katz and A.~Ringwald,
%``Determination of absolute neutrino masses from Z-bursts,''
%arXiv:
;
%%CITATION = ;%%
%
%\cite{Fodor:2001rg}
%\bibitem{Fodor:2001rg}
%Z.~Fodor, S.~D.~Katz and A.~Ringwald,
%``Possible detection of relic neutrinos and determination of their mass:  Quantitative analysis,''
%arXiv:
; 
%%CITATION = ;%%
%
%\cite{Ringwald:2001mx}
%\bibitem{Ringwald:2001mx}
A.~Ringwald,
%``Possible detection of relic neutrinos and their mass,''
%arXiv:
;
%%CITATION = ;%%
%\cite{Kalashev:2001sh}
%\bibitem{Kalashev:2001sh}
O.~E.~Kalashev, V.~A.~Kuzmin, D.~V.~Semikoz and G.~Sigl,
%``Ultra-High Energy Cosmic Rays from Neutrino Emitting Acceleration Sources?,''
%arXiv:
.
%%CITATION = ;%%

%\cite{Gandhi:1996tf}
\bibitem{Gandhi:1996tf}
R.~Gandhi, C.~Quigg, M.~H.~Reno and I.~Sarcevic,
%``Ultrahigh-energy neutrino interactions,''
Astropart. Phys.  {5} (1996) 81.
%.
%%CITATION = ;%%


\bibitem{Seckel:2001}
D.~Seckel {\em et al.} [RICE Collaboration], 
in: Proc. 27th International Cosmic Ray Conference, Hamburg, Germany, 2001, 
Vol. 3, pp. 1137-1140. 


%\cite{Groom:2000in}
\bibitem{Groom:2000in}
D.~E.~Groom {\it et al.}  [Particle Data Group Collaboration],
%``Review of particle physics,''
Eur. Phys. J. C {15} (2000) 1.
%%CITATION = EPHJA,C15,1;%%

%\cite{Alvarez-Muniz:2001es}
\bibitem{Alvarez-Muniz:2001es}
J.~Alvarez-Muniz and F.~Halzen,
%``10**20-eV cosmic-ray and particle physics with kilometer-scale neutrino  telescopes,''
Phys.\ Rev.\ D {63} (2001) 037302.
%.
%%CITATION = ;%%

%\cite{Bugaev:1998bi}
\bibitem{Bugaev:1998bi}
E.~V.~Bugaev, A.~Misaki, V.~A.~Naumov, T.~S.~Sinegovskaya,
S.~I.~Sinegovsky and N.~Takahashi,
%``Atmospheric muon flux at sea level, underground and underwater,''
Phys.\ Rev.\ D {58} (1998) 054001.
%.
%%CITATION = ;%%

%\cite{Gelmini:2000ve}
\bibitem{Gelmini:2000ve}
G.~Gelmini, P.~Gondolo and G.~Varieschi,
%``Prompt atmospheric neutrinos and muons: NLO vs LO {QCD} predictions,''
Phys.\ Rev.\ D {61} (2000) 036005.
%.
%%CITATION = ;%%


%\cite{Hundertmark:2001}
\bibitem{Hundertmark:2001}
S.~Hundertmark {\em et al.} [AMANDA Collaboration], 
in: Proc. 27th International Cosmic Ray Conference, Hamburg, Germany, 2001,
Vol. 3, pp. 1129-1132.





%\cite{Alvarez-Muniz:2000xx}
\bibitem{Alvarez-Muniz:2000xx}
J.~Alvarez-Muniz, R.~A.~Vazquez and E.~Zas,
%``Characterization of neutrino signals with radiopulses in dense media  through the LPM effect,''
Phys.\ Rev.\ D {61} (2000) 023001.
%.
%%CITATION = ;%%

\end{thebibliography}

\end{document}

%%% Local Variables: 
%%% mode: latex
%%% TeX-master: t
%%% End: 



