High Energy Neutrinos with a Mediterranean Neutrino Telescope
E. Borriello, A. Cuoco, G. Mangano, G. Miele, S. Pastor, O. Pisanti, P. D. Serpico
aa r X i v : . [ a s t r o - ph ] S e p TH I NTERNATIONAL C OSMIC R AY C ONFERENCE
High Energy Neutrinos with a Mediterranean Neutrino Telescope
E. B
ORRIELLO , , A. C UOCO , G. M ANGANO , G. M IELE , , S. P ASTOR , O. P ISANTI ,P. D. S ERPICO Dipartimento di Scienze Fisiche, Universit`a di Napoli “Federico II” and INFN Sezione di Napoli, Com-plesso Universitario di Monte S. Angelo, Via Cintia, Napoli, 80126, Italy Instituto de F´ısica Corpuscular (CSIC-Universitat de Val`encia), Ed. Institutos de Investigaci´on, Apdo.22085, E-46071 Valencia, Spain Institut for Fysik og Astronomi, Aarhus Universitet Ny Munkegade, Bygn. 1520 8000 Aarhus Denmark Center for Particle Astrophysics, Fermi National Accelerator Laboratory, Batavia, IL 60510-0500, USA [email protected]
Abstract:
The high energy neutrino detection by a km Neutrino Telescope placed in the Mediterraneansea provides a unique tool to both determine the diffuse astrophysical neutrino flux and the neutrino-nucleon cross section in the extreme kinematical region, which could unveil the presence of new physics.Here is performed a brief analysis of possible NEMO site performances.
Introduction
Neutrinos are one of the main components of thecosmic radiation in the high energy regime. Al-though their fluxes are uncertain and depend on theproduction mechanism, their detection can provideinformation on the sources and origin of the highenergy cosmic rays.From the experimental point of view the detectionperspectives are stimulated by the Neutrino Tele-scopes (NT’s) constructed, like
Baikal [1] and
AMANDA [2], or under construction like
IceCube [3] under the ice and
ANTARES [4] in the deep wa-ter of the Mediterranean sea. Here, also the exper-iments
NESTOR [5] and
NEMO [6] are in the R&Dphase and, together with
ANTARES , in the futurecould lead to the construction of a km telescopeas pursued by the KM3NeT project [7].Although NT’s were originally thought as ν µ de-tectors, their capability as ν τ detectors has becomea hot topic in view of the fact that flavor neu-trino oscillations lead to nearly equal astrophysi-cal fluxes for the three neutrino flavors. Despitethe different behavior of the produced tau leptonswith respect to muons in terms of energy loss anddecay length, both ν µ and ν τ event detection rates are sensitive to the matter distribution near the NTarea. In principle, the elevation profile of the Earthsurface around the detector may be relevant. InRef. [8], some of the present authors calculatedthe aperture of the Pierre Auger Observatory [9]for Earth-skimming UHE ν τ ’s, by using the Digi-tal Elevation Map (DEM) of the site (GTOPO30)[10]. In Ref. [11] the DEM’s of the under-waterEarth surface, provided by the Global Relief Datasurvey (ETOPO2) [12], was used to estimate theeffective aperture for ν τ and ν µ detection of a km NT in the Mediterranean sea placed at any of thethree locations proposed by the
ANTARES , NEMO and
NESTOR collaborations. In the present pa-per we further develop the approach of Ref. [11]to evaluate the performances of a MediterraneanNT in the simultaneous determination of the neu-trino flux and the ν -Nucleus cross section in ex-treme kinematical regions (which may probe newphysics, see e.g. [13]). Since the three differentproposed sites for the under-water km telescopeshow event rate differences of the order of 20%, forthe sake of brevity we report the results our analy-sis for the NEMO site only, which presents interme-diate performances.
IGH E NERGY N EUTRINOS WITH A M EDITERRANEAN N EUTRINO T ELESCOPE
Formalism and results
Following the formalism developed in [8, 11] wedefine the km NT fiducial volume as that boundedby the six lateral surfaces Σ a (the index a =D, U, S,N, W, and E labels each surface through its orien-tation: Down, Up, South, North, West, and East),and indicate with Ω a ≡ ( θ a , φ a ) the generic direc-tion of a track entering the surface Σ a (see Figure4 of Ref. [11] for notations). We introduce all rele-vant quantities with reference to ν τ events, the caseof ν µ being completely analogous.Let dΦ ν / (d E ν dΩ a ) be the differential flux ofUHE ν τ + ¯ ν τ . The number per unit time of τ lep-tons emerging from the Earth surface and enteringthe NT through Σ a with energy E τ is given by (cid:18) d N τ d t (cid:19) a = Z dΩ a Z d S a Z d E ν dΦ ν ( E ν , Ω a )d E ν dΩ a Z d E τ cos ( θ a ) k τa ( E ν , E τ ; ~r a , Ω a ) . (1)The kernel k τa ( E ν , E τ ; ~r a , Ω a ) represents the pro-bability that an incoming ν τ crossing the Earth,with energy E ν and direction Ω a , produces a τ -lepton which enters the NT fiducial volume th-rough the lateral surface d S a at the position ~r a withenergy E τ . For an isotropic flux and an exposuretime T , the total number of τ leptons (and similarlyfor muons) crossing the NT is N τ = T X a Z dΩ a Z d S a Z d E ν Z d E τ (cid:18) π dΦ ν ( E ν )d E ν (cid:19) cos ( θ a ) k τa ( E ν , E τ ; ~r a , Ω a ) . (2)Although the exact dependence of Eq. (2) onthe neutrino flux and the neutrino-nucleon chargedcurrent cross section σ νNCC may be quite compli-cated, basic physical considerations show that evena rough binning of the events for energy loss andarrival direction may be used to obtain informationon both these quantities (see e.g. [14, 15]). In par-ticular, in the following we shall consider the sumof the µ and τ contributions as the experimentalobservable, namely the energy deposited in the de-tector and not the energy and/or the nature of thecharged lepton crossing the NT. In fact, only for aminor fraction of the detected events the nature ofthe charged lepton can be reliably established. According to Ref.s [16, 17], the differential energyloss of the τ leptons per unit of length in an un-derwater NT can be simply taken as d E τ / d λ = − β τ E τ ̺ w , with β τ = 0 . × − cm g − and ̺ w denoting the water density. Analogously, formuons one just needs to replace β τ with the corre-sponding value β µ = 0 . × − cm g − . As-suming that the lepton energy loss in the NT bye.m. interactions, ∆ E l , is just a small fraction ofits energy at the entrance, E l , we simply obtain ∆ E l = λ ( ~r a , Ω a ) β l E l ̺ w , where λ ( ~r a , Ω a ) is thelength crossed in the NT by the lepton whose trackis defined the geometrical quantities ~r a , Ω a .Using these relations one can derive the spectrumof leptons detected in the NT as a function of theirdeposited energy, ∆ E , and their arrival direction, Ω ≡ ( θ, φ ) , measured in the zenith-azimuth refer-ence frame d N d(∆ E )dΩ = T X α = µ,τ X a Z d S a Z d E ν π dΦ ν ( E ν )d E ν cos ( θ a ) k αa λ ( ~r a , Ω a ) β α ̺ w . (3)By denoting with X i a given bin in energy loss,and with Y j the one for the zenithal angle, we canintegrate the expression (3) to get the number ofexpected events in X i × Y j , N ij = T X α = µ,τ X a Z X i d(∆ E ) Z Y j dΩ (4) Z d S a Z d E ν π dΦ ν ( E ν )d E ν cos ( θ a ) k αa λ ( ~r a , Ω a ) β α ̺ w . To take into account the underwater surface profileone can numerically compute the above integral asdescribed in Ref. [8]: by using the available DEMof the area near the
NEMO site, one can isotrop-ically generate a large number of oriented trackswhich cross the
NEMO fiducial volume (see Figure4 of Ref. [11]) and sample the above integrand.This technique allows also to account for the radialdensity profile of the Earth (we use the formula re-ported in [18]).In order to study the sensitivity to both neutrinoflux and σ νNCC it is necessary to parameterize theirstandard expressions and the possible departuresfrom them. In particular, we parameterize the fluxas dΦ ν / d E ν dΩ a = C · . · − ( E ν / GeV ) − D TH I NTERNATIONAL C OSMIC R AY C ONFERENCE Π€€€€€€ Π€€€€€€
Π€€€€€€€€€€€ ΠΘ dN i €€€€€€€€€€€€€ d Θ N = GeV < D E l < GeV N = N = N = D E l > GeV N = N = j = = Figure 1: Angular distributions of ( µ + τ ) eventscollected in five years from a km NT placed at theNEMO site (see text).GeV − cm − s − sr − , which gives a standardWaxman-Bahcall flux [19] for C = D = 1 . Forthe neutrino-nucleon cross section we use: σ νNCC . nb = (cid:16) E ν E (cid:17) . A E ν ≤ E (cid:16) E E (cid:17) . A (cid:16) E ν E (cid:17) . B E ν > E where E = 10 . GeV is the energy below whichthe atmospheric flux is expected to dominate (sowe consider only the region
E > E ) and E =10 . GeV. In the low-energy bin this cross-sectionmatches the standard expression [20] for A = 1 .A value of B significantly larger than 1 may beassociated with new physics. Note that the factor C only enters via the product CT as a normalizationand can be fixed to C = 1 , considering instead theexposure time T as the effective variable.For illustrative purposes, in Figure 1 we report theevent angular distribution, for a km NT placed atthe
NEMO site in five years of operations. The solidand dashed lines correspond to events whose en-ergy loss in the detector belongs to the intervals . - GeV or > GeV, respectively. Thepredictions are obtained for standard flux and crosssection ( A = B = C = 1 ). In the plot are also re-ported the number of events N ij (see Eq. (4)) whenwe consider i = 1 , for the previous two energy D - C T H y r L D < N j = N j = æ D - C T H y r L Figure 2: ( CT, D ) region corresponding to the ob-servation of at least one event in each bin (for stan-dard cross section).bins and j = 1 , when the zenith arrival directionis between 0 ◦ and 90 ◦ or 90 ◦ and 180 ◦ .Clearly, for a very steep flux power-law index D ,the number of events decreases. We shall re-quire that at least one event falls in each bin, in T years of running, in the case of standard cross-section; this rough criterion constrains the param-eter range that one experiment is able to explore tothe brighter region of Fig. 2, corresponding to theintersection of the regions where N ij ≥ , for all i, j .As a preliminary result, in Fig. 3 we show theconstraints (contours at the 68 % and 95 % CL)which can be obtained on the physical parameters A and D after the marginalization over C is made.Here we are assuming B = A , so that the plotrepresents the capability of the telescope to dis-entangle the energy dependence of the flux fromthe energy dependence of the cross-section (in thetoy model where both are described by a singleparameter). We performed a multi-Poisson like-lihood analysis [21], in which the likelihood func-tion, L = exp( − χ / , is defined using the fol-lowing expression for the χ ( N ij being the eventnumbers of the reference model): χ = 2 X ij (cid:2) ( N ij − N ij ) + N ij ln( N ij /N ij ) (cid:3) . (5) IGH E NERGY N EUTRINOS WITH A M EDITERRANEAN N EUTRINO T ELESCOPE A D
95% 68%
Figure 3: Marginalized contour levels in the ( A, D ) plane (for A = B ) (see text for details). Conclusions
We have performed an analysis of the capability ofa km NT in the Mediterranean to disentangle thehigh energy neutrino flux and neutrino-nucleoncross section in an unexplored kinematical region.Our statistical analysis exploits the dependence ofobservables on energy and arrival direction (underthe hypothesis of an isotropic diffuse flux). Usinga simplified toy model to parameterize fluxesand cross-sections, preliminary results confirmthat this approach is very promising, and couldpotentially detect hints of new physics. Of coursethe real feasibility of such measurements willdepend crucially on the size of the neutrino fluxwhich fixes the time required to reach a reasonablestatistics. A complete account of this research willbe reported in a forthcoming publication.
Acknowledgements:
P.S. acknowledges sup-port by the US Department of Energy and byNASA grant NAG5-10842. G. Miele acknowl-edges support by Generalitat Valenciana (ref.AINV/2007/080 CSIC).
References [1] V. A. Balkanov et al. [Baikal Collabora-tion], Phys. Atom. Nucl. , 951 (2000)[Yad. Fiz. , 1027 (2000)] [arXiv:astro-ph/0001151]. [2] E. Andres et al. , Astropart. Phys. , 1 (2000)[arXiv:astro-ph/9906203].[3] J. Ahrens et al. [IceCube Collaboration], As-tropart. Phys. , 507 (2004) [arXiv:astro-ph/0305196].[4] M. Spurio [ANTARES Collaboration],arXiv:hep-ph/0611032.[5] G. Aggouras et al. [NESTOR Collaboration],Nucl. Instrum. Meth. A , 452 (2006).[6] P. Piattelli [NEMO Collaboration], Nucl.Phys. Proc. Suppl. , 172 (2007).[7] U. F. Katz, Nucl. Instrum. Meth. A , 457(2006) [arXiv:astro-ph/0606068].[8] G. Miele, S. Pastor and O. Pisanti, Phys. Lett.B , 137 (2006) [arXiv:astro-ph/0508038].[9] J. Abraham et al. [Pierre Auger Collabora-tion], Nucl. Instrum. Meth. A , 50 (2004).[10] U.S. Geological Survey’s Center for EarthResources Observation and Science (EROS),1996, http://asterweb.jpl.nasa.gov [11] A. Cuoco et al. , JCAP , 007 (2007)[arXiv:astro-ph/0609241].[12] National Geophysical Data Center, 2001, [13] J. Alvarez-Muniz, F. Halzen, T. Han andD. Hooper, Phys. Rev. Lett. , 021301(2002) [arXiv:hep-ph/0107057].[14] D. Hooper, Phys. Rev. D , 097303 (2002)[arXiv:hep-ph/0203239].[15] S. Hussain, D. Marfatia, D. W. McKay andD. Seckel, Phys. Rev. Lett. , 161101 (2006)[arXiv:hep-ph/0606246].[16] C. Aramo et al. , Astropart. Phys. , 65(2005) [arXiv:astro-ph/0407638].[17] S. I. Dutta, Y. Huang and M. H. Reno,Phys. Rev. D , 013005 (2005) [arXiv:hep-ph/0504208].[18] R. Gandhi, C. Quigg, M. H. Reno andI. Sarcevic, Astropart. Phys. , 81 (1996)[arXiv:hep-ph/9512364].[19] E. Waxman and J. N. Bahcall, Phys. Rev. D , 023002 (1999) [arXiv:hep-ph/9807282].[20] R. Gandhi, C. Quigg, M. H. Reno andI. Sarcevic, Phys. Rev. D , 093009 (1998)[arXiv:hep-ph/9807264].[21] S. Baker and R. D. Cousins, Nucl. Instrum.Meth. A221