arXiv:0711.3635v2 [astro-ph] 15 Oct 2008 i’ hTF ZhET Pis’m v a detected s R C per source,it years any take ay m posures until m ill w inevitably exposuresky is increase yet not - large E yet ould in H this ed w be U hile of enough. ber num scenario the W larger ex- that the accum ulated ofthe order sam be to e agnitude ofm [6]. ections e D alactic the in G agnetic m eld are expected ACS:98.70.Sa C PA f Gis typically nG of1 angularthe agnetic spread m elds random caused by error in the energy determsources from ination. com e closer observed EeV an tons with ver energy pc,even 80 such O of than distances, 100 M accounting ogeneouslyhom for 20% a throughoutthetypical icrow universe,70% energy m ic spectram instance, background ). ave For B sources (C ofthe and M pro- for distributedreduced their by interactions roughly ith w photons from the cos- should particularly be ni cantly ected de their [5]at cuto from true agnetic initial above sig- ZK arem not the G elds s is allso ECR sm anticipate H that U ay direction. m one his T that in the uence extragalactic of although active this and s correlations galactic ECR tion H ofU nuclei(AG [3],(A uger) ) N has been contested [2]. Second,the ellasw [4]. uger the in A data data ierre from the hint P hus for T correla- bservatory uger a O A ing all available experimanisotropies \old" of R data ents C [1] \charged of as particle daw n the at edium scales m on astronom have First, been foundthere y." com bin- pieces are tw o of evidencepossibility indicatingprogress identify to are that e cosm physics w in individual ic ) (C ray R im the been - has ow ever, sources. H eGZ oio n constraints source horizon spec ic ray cosm and the on ZK G he T 5 emanrao yn sources ain reason no m hy he w have identi-T been Introduction| 10 index. particular, In typical that for show probability a e w events source two source per or densi applied only one obse are be hen for given w even a source populate to ultrahigh dierent energy supposedly s) from energy cosm R C ic rays E b (U H iden ithw above 6 19 V,hntherangeofUeV ,when sissigni ECR H cantly ics h Khorizon c ethod to protons present discuss of m e and ZK a G the W E 2 : 7 10 t9%condence 95% at level. eo h i obstacles ain m fast the of to ne O 19 < Valwsoei 7 of allcases allow 97% eV in one distinguish to s source a spect 2 3 C,0reAlc mo tLoi uquet,F{75205 Leonie et D lice Paris on A ,10,rue ede om C PC D A 4 for such high-energy protons. NRRAS6t ctober nniversary O S,60th A A R prospectIN R osco 7a,117312 M r mosrain nteGZ regim e observations ZK from G the in Kachelrie .K M 3 1 er ivision,C heory D T eneva itzerland 23,Sw {1211 G N H R E C Institutt Trondheim {7491 for orw fysikk,N ay ,N ,N U N T 1 ,E.Parizot 1 2 adD. Smikoz .Sem .V D ,and s p e c t r ai n d icannotbe in v know i advance n d u a (i.e. l lthat,even y ) before , t h er though e easuring l the constrain m to relative a t i their v i ewe the spectral g h eightfor to f d w i the near er e future ofdi index n te n erent - of individual sourceconstraints ofproton sources spectrum R astronom C E H suitable sources U the , basic he y. T on idea is hand, pinch accelerationgion in serve exam the usualshock ple for ay an as m accelerationcount large a photon background [12]. the other in the acceleration n O re- notherA possibilitythe observed obtain properties to black holes of suggested largeacceleration scale ef. R in [11]. [9, jets N 10] that G in A explains electric the in also eld around superm assive r e s u l t si n determ be can ined best.source are detected,lation the injection ask how to studies ediatem typically spectrum but phase sources hen w identi be ay e m events s only tw corre- or o one s by ed e per n t ithe a source l l the k and di n o yun use wn.I ti spectra st h e r e at higher f o r energies et i mel y ,i are ei nth n t e r - x mpeframodelwithexam m ple for a that predict either harder uch softer m or spectra. n A distributionwi axim t h of m al energies below dividualspectrumtense di he use T energy . source spectrum in s theenough R ofC sky events to allow are one draw accum ing ulated decent a ost from even in- in- the m racy requires and generation a spectrum 2 = ins,depending the injection actualC on R spectral In this ork,w w presentan e alternative ethod m to set hile rst-orderW Ferm shock i acceleration typically yof4 ty : tied extragalactic sources. ethod can his T m 7[ 1 3 ] . E < v dadi ae nteanalysis the is based on rved and the of 2 4 : around 2.1 [8], there exist various odels m for identical7 sources| appropriate an or 10 onstrain injection the spectrum of 10 5 19 pc M 2 u rum Vi nwnwith w is reasonable know n eV accu- ; 3 13,Francex Russia ,R w ; 4 3 N=E =dE dN adt e f10events data set 100 of ,a / pt 10 to up 1 E E ax m is take to into1 ac- 1 : 1 [7]| hile w both r monefrom N=dE dN 20 Vis theeV r mtrum / E 2

400 1000 10 % 90 % Kalashev et al 30 % 350 90 % Harari et al 50 % CEL 70 % CEL no redshift 90 % 300 SOPHIA ∆E/E=0% SOPHIA ∆E/E=20% 100 250

200 R [Mpc] R [Mpc] 150 10 100

50

0 1 60 80 100 120 140 1e+19 1e+20 1e+21 E [EeV] E [eV]

Fig.1:D istance R inM pcfrom w hich 90% ofU H EC R s Fig.2:T he distance R inM pc forw hich a certainfrac- arrive w ith energy >E asfunction ofthe thresholden- tion f ofU H EC R sarrivesw ith energy >E asfunction 21 ergy E for Em ax = 10 eV and = 2:7.T he thinsolid ofthe energy threshold E for = 2:7.From top tobot- red line usesC EL ina staticU niverseas[19],the green tom , f = 90% asred line, f = 70% aspink, f = 50% as line usesC EL inan expanding U niverse.T he blue line m agenta, f = 30% asblue and f = 10% as green line. labeled \SO PH IA" hastobe com pared to[20].T he red line takesintoaccountadditionallyan experimentalen- ergy resolution E =E = 20 % . ityofthecorresponding energy losses(dotted,blueline) or applying the continuous energy loss(CEL) approx- + ergy bins for a given source isa directconsequenceof imation to itsresults(dashed,green line). The e e the sourcespectrum . N ow supposethata m inimalen- pairproduction lossesweretaken from R ef.[17]. ergy Em in can be identied,above which we can trust The f = 90% horizon com puted withinthe CEL ap- that the observed CR s com e roughlyin straight lines proximation underestimatesconsiderablythefullM onte from theirsourceand,m ostimportantly,sourcesinside Carloresult.The dierenceincreasesfora larger\hori- the horizon appearwitha smallenough angularspread zon fraction", f ! 1, and as function of energy for on the sky that they do not overlap. The energy dis- E ! Em ax.Therearetwo reasonsforthelatterdiscrep- tribution ofCR s seen above Em in from a given source ancy. First,the energy transferperinteraction, y,in- should then re ect the source spectrum (modied by creaseswithenergyand violatesm oreand m orestrongly the usualpropagation e ects),and even ifone observes the formalrequirem ent y 1 needed forthe applicabil- onlyone ofthem ,itsenergy containssom e information ityofthe CEL approximation. Second,the ux taking about the source spectrum . W e show how thissimple intoaccountthestochasticnatureoftheenergylossesin argum entcan beimplem entedquantitativelyfora given pion production rem ains nitefor E ! Em ax,whilein dataset,taking intoaccountU H ECR energylossesfrom theCEL approximation no particleswith E = Em ax can pureproton sourceswithsupposedlyidentieddistances reachtheobserverfrom a sourceata nitedistance[18]. and identicalm aximum energy.W e usethistoy m odel In a realisticexperiment,the primary energy can toillustratethe basicfeaturesofthe m ethod and toex- onlybe reconstructed with a niteprecision. A ssum - ploreitspotentialpower,leaving necessaryre nem ents ing a G aussian (in log E) experimentaluncertainty of forfuturework. E =E = 20% ,we com puted the 90% horizon asa func- Propagation and horizon scaleofU H E protons| In tion ofthe m easured CR energy,forthesam econditions Fig.1,we show the \90% horizon" { i.e. the distance as above. The two resulting curvesarealso shown in R90 from which 90% ofthe U H ECR s observed above a Fig.1. Since the CR spectrum isfalling steeply,the given energy, E,originate| as function ofenergy. W e m isinterpretation oflowerenergy eventsashigh energy assum e a uniform source distribution with a density oneshasa largerimpactthan the reverse,whichinturn 5 3 ns = 410 =M pc (cf.e.g.Refs.[14,15])and a power- leads toan increaseofthe estimated horizon scale.At 19 law source spectrum dN =dE / E with = 2:7 up low energies,say <5 10 eV,the observed spectrum 21 to the m aximalenergy Em ax = 10 eV. W e used for approximateswellto a power-law and the energy res- the calculation ofphoto-pion production the program olution onlya ectsthe absolute ux,not the relative SO PH IA [16],eithertaking intoaccountthe stochastic- uxesrelevantfor R90(E). G ZK horizon and constraintson the cosmicray sourcespectrum from observationsinthe G ZK regime 3

The horizon scalefor U H E protons and nucleiwas energy rangefrom a few1017 eV ora few1018 eV up recentlydiscussed also inR efs.[19,20]. In Fig.1 we to the G ZK cuto , depending on the assum ed source com pareourcalculationstothoseofR efs.[19]and [20] com position [21,22]. for proton primaries. In R ef.[19],H arari etal. pre- W e rstproduceaM onteCarlo(M C)sam pleby gen- sented results(shown asorangeline)using the CEL ap- erating sources with constant com oving density ns = proximation and assum ing a staticU niverse,ourresult 4 105 M pc3 up to a m aximalredshiftof z = 0:1. forthe sam e assum ptionsisshown witha thinsolidred Then we choose a source i according to the declina- line.Both calculations agreewellatm oderateenergies tion dependent exposureofA uger,with an additional E = 80 100 EeV,whilethereissom e disagreem ent weight chosen according to the source distance. Fi- bothathigh and low energies.H owever,the dierences nally,we generatea CR with an initialenergy drawn atlow energiesbetween the two calculations arem uch random lyaccording tothe assum ed injection spectrum , smallerthan the dierencesbetween thosecalculations dN =dE / E 0 ,and propagateituntiliteitherreaches and the m ore correct CEL calculation in the CD M the Earthdistanceorlosesenergy down tobelow Em in. m odelfor the expanding U niverse, presented with a Inthe formercase,we then applyan energy-dependent green line. angular de ection to m imicthe e ect ofthe G alactic A llresultsusing the CEL approximation dier in m agnetic eld,witha shiftperpendiculartothe G alac- shape asa function ofenergy from the calculationssus- ticplane equal to b = 2(E =1020eV)1 , where this ing SO PH IA for pion production either directly(blue m agnitude ism otivated by the resultsofR ef.[6].The line),orusing the SO PH IA resultsina kineticequation chosenm agnetic eldlikelyoverestimatesde ectionsfar approachasinR ef.[20](magentaline).The agreem ent away from the galacticplane inm ostofm odels.H ow- between the lattertwo resultsisalmostperfectat all ever,we consider thischoice as a conservative upper energies. limit. Finally,we de ectthe CR direction to account A san illustration,weshow inFig.2 the horizon dis- fora niteexperimentalangularresolution,taking the tancecorresponding todierent CR fractions.Speci- A ugersurfacedetectorasa reference[23],witha spher- 2 2 cally,we plotthe distance Rf below which a given frac- icalG aussian density / exp(‘ =(2l ))sin(‘)d‘,where tion f ofthe U H ECR s reach the Earth with an energy ‘ = 0:85 and ‘ isthe angulardistance. largerthan E,asa function ofthatenergy,for f = 10% , A fterhaving generated N cosmicrays,we perform 30% , 50% , 70% and 90% (using always SO PH IA and a correlation analysisbetween the CR sand the sources. E =E = 20% ). First,we identify as \the source" of a given CR the Estimation of the spectral index| Since the angu- source with the smallestangular distance ‘ to the ob- larresolution ofcosmicray experimentsispoorby as- served CR arrivaldirection and m aximaldistance R = tronom ical standards,the identication of individual 100M pc. Inside thisregion, there are around 160 5 3 sourcesrequiresa relativelylarge angular distance be- sourcesforchosen density ns = 4 10 M pc . Such tween them . Thiscan only hold for su ciently high a smallnum ber m akes the probability negligiblethat energiessuch thatthe horizon scaleissmall,say ofthe sourcesoverlap,ifthey areuniformlydistributed. This orderof100 M pc,leaving a limited num berofsources probability increases,ifsourcesfollow| as expected| overthe sky.D e ning ashorizon,withinwhich 90% of the large-scalestructureofm atterand m ay constitutea allCR sobservedabovea givenenergy wereem itted,we reallimitation toresolvesinglesourcesinclustercores. nd from Fig.2thata horizon of100 M pc corresponds A dditionally,we requirethatthe angulardistance ‘ 20 to a thresholdenergy of E = 1 10 eV.At present, be smallerthan a prescribed value, ‘m ax. N ext,having the importanceofde ections inextragalacticm agnetic pre-de ned an energy E2 thatdividesthe wholeenergy elds above thisenergy isunclear.A s soon as sources rangeintotwo largebins,wecountforeachsource i the aredetected,onewillbeabletosetan upperlimitand to num bers Ni;1 and Ni;2 ofhigh energy (E E2)and low a certainextentreconstructthe extragalacticm agnetic energy events(Em in E

thatthe observed num bers Ni;j areconsistentwith the tively,or the m aximalenergy ofsourcesisdistributed 3:6 value 0 used inthe M C.Considered as a function of as dn=dEm ax / Em ax [7].M oreover,we consideronly ,thisprobabilitydistribution hasthe true value 0 as two extrem e cases,nam elya power-law with 0 = 1:1 itsexpectation value,ifour procedureisunbiased,and and 0 = 2:7. m easures how stronglythe data disfavor a dierently InFig.3 we com parethe distributions ofprobabili- assum ed value 6= 0. tiesobtained from Eq.(2)choosing astrue value 0 = 5 3 Since the dierent sourcesem itCR s independently 2:7,as source density as always ns = 4 10 =M pc , from one another,we can simply m ultiply the single as num berofCR s N = 100,and ‘m ax = 4 . The red sourceprobabilities pi(Ni;1;Ni;2j)toobtainthe global solidlineisthedistribution ofprobabilitiesobtainedas- probabilityofagivendatasetwith Ns identiedsources: sum ing = 1:1,whilethe blue dashed line corresponds = 2:7.Thetwo curveshaveonlya smalloverlap,since YNs theprobabilitiesusing thecorrect arerathernarrowly p(fNi;1gi= 1;Ns j)= pi(fi;1;fi;2j): (2) i= 1 concentrated around p = 1,whilethe probabilitydistri- bution using the wrong extends from extrem elylow The basicoutcom e ofa sam pleofM C simulations valuesup toone.Thus an experimentaldierentiation for xedparam eters = 0 ::: isthusa binneddistribu- between dierentinjection spectraseem spossible,even tion, f(pj),giving the fraction f ofM Csproducing the ifonlyone or,infew cases,two eventspersource are value p. W ith how m uch con dencecan we distinguish detected,asitisthe caseforthe chosen param etersin thesedistributionsfortwo dierent 1 and 2? Clearly, Fig.3.Thisconstitutesthe m ainresultofourwork. thesmallertheoverlap ofthetwo distributions,theeas- W e quantify the chances to distinguish two dier- ierthe two param etersets can be distinguished. i ent spectralindicesinthe following way: W e calculate

1 the area A corresponding tothe desired con dencelevel 1.1 real 2.7 2.7 real 2.7 (C.L.), A, starting from 1 to the leftusing the best- p 95 t distribution (e.g.the blue line inFig.3)and obtain thereby asitslowerboundary pA. Thus onlyin1 pA 0.1 caseswe willobtainby chancea lowerprobabilityusing the correcttesthypothesis.N extwe counthow largeis

P(p) the area B ofthe wrong testhypothesison the leftof 0.01 pA.A s nalanswerwe obtainthatinthe fraction B of allcaseswecan distinguishbetween the two hypotheses with C.L. A.

0.001 Letus illustratethisprocedureforthe caseconsid- ered above,choosing ascon dencelevel A = 95% . The 1e-10 1e-08 1e-06 1e-04 0.01 1 greendashed-dottedverticallineinFig.3 enclosing 95% p ofthe area ofthe true (blue) distribution determines Fig. 3: D istribution of probability of reconstructed p95 = 0:056. The area ofthe red curve on the leftof power law spectrum ifreal power law spectrum is p95 = 0:056 is B = 0:971. H ence one can exclude in = 2:7,angle ‘ m ax = 4 . In allcases Em in = 60 EeV. B = 97:1% ofcaseswith at least95% C.L.the expo- T he red line isfor = 1:1 and the blue line for = 2:7. nent = 1:1 for the spectrum ,ifthe true exponent is 0 = 2:7. W e study now the possibility to distinguish dier- Inaddition tothe ratherextrem e casesofthe spec- ent values of the injection spectrum of CR s in m ore tral indices above,we investigated the ability of the detail.A s simplifying assum ption we assum e that the m ethod todistinguish between any ofthem and an in- injection spectrum ofallsourcesisthe sam e,i.e.inpar- termediatevalue of 0 = 2:0,often considered in the ticular that the m aximalenergy ofallsourcesisiden- contextofastrophysicalparticleacceleration. A s an il- tical.Thisassum ption allows us to study the spectra lustration,we found thatwith a datasetof100 cosmic onlyabove 4 1019 eV,becauseat lowerenergiesa raysabove6 1019 eV,itispossiblein50% ofthe cases spectralindex < 2:6 requireseitheradditionalG alac- todiscriminate 0 = 2:0 from a valueofeither1.1or2.7 ticsourcesora non-uniform sourcedistribution.Inthe with a C.L.of95% .Likewise,fora datasetof200 cos- latter case,either the source density or the lum inos- m icraysabove4 1019 eV (i.e.foressentiallythe sam e ityofsinglesourcesshouldincreaseasfunction ofred- exposureofthe sky,butwitha lowerenergy threshold), m m shift, n(z)= n0(1+ z) and L(z)= L0(1+ z) respec- an injection spectralindex of1.1can be discriminated G ZK horizon and constraintson the cosmicray sourcespectrum from observationsinthe G ZK regime 5

against 0 = 2:0 witha C.L.of95% in90% ofthecases, arXiv:0804.1088 [astro-ph]; I. V . M oskalenko et al., whilean injection spectralindex of2.7can be discrim- arXiv:0805.1260 [astro-ph]. inated against 0 = 2:0 with a C.L.of95% in70% of 5. K .G reisen,Phys.R ev.Lett. 16,748 (1966);G .T .Zat- the cases. sepin and V . A . K uzm in, JET P Lett. 4, 78 (1966) Sum m ary|W e haveproposed a m ethod toestimate [Pisma Zh.Eksp.Teor.Fiz. 4,114 (1966)]. the generation spectrum ofindividualextragalacticCR 6. M . K achelrie, P. D . Serpico and M . Teshima, A s- sourcesthat iswell-suited for the time when onlyone tropart.Phys. 26,378 (2006) [astro-ph/0510444]. ortwo eventspersourcearedetected.A n importantin- 7. M .K achelrie and D .V .Sem ikoz,Phys.Lett.B 634, gredientofthism ethod isthe relativefraction ofevents 143 (2006)[astro-ph/0510188]. contained ina prescribed energy interval.Thereforewe 8. V .S.Berezinskii etal,A strophysicsofcosmicrays,A m - haverecalculated the horizon scaleofultra-high energy sterdam : N orth-Holland 1990.T .G aisser,C osmicR ays and ParticlePhysics,C am bridge U niversityPress1991. protons,taking into account a reasonableenergy reso- R .J.Protheroe and R .W .C lay,Publ.A stron.Soc.of lution,similartothatofA uger. A ustralia 21,1 (2004) [astro-ph/0311466]. W e havedem onstratedfora toy-modelthe potential 9. A .N eronov and D .Sem ikoz,N ew A stronom y R eviews, ofthism ethod, nding that around 100 eventsabove 47,693 (2003),A .N eronov,D .Sem ikozand I.T kachev, 19 6 10 eV arerequired todistinguish with 97% prob- arXiv:0712.1737 [astro-ph] abilityat leastat the 95% C.L.the two extrem e cases 10. A .N eronov,P.T inyakov and I.T kachev,J.Exp.T heor. = 1:1 and 2.7.A dierentiation between ’sthatare Phys. 100, 656 (2005) [Zh.Eksp.Teor.Fiz. 100, 744 m oresimilarwillbe clearlym orechallenging.A n injec- (2005)][astro-ph/0402132]. tion spectralindex of2.0can stillbe distinguished from 11. A .N eronov,D .Sem ikoz,F.A haronian and O .K alashev, the two above valueswith a 95% C.L.inthe m ajority Phys.R ev.Lett. 89,051101 (2002) [astro-ph/0201410]. ofcases(with the sam e statistics). 12. E.V .D erishev etal.,Phys.R ev.D 68,043003 (2003) Severaloftheissueswehaveneglected,likethee ect [astro-ph/0301263]. ofa possible Em ax distribution,shouldbe included ina 13. V .V .V lasov,S.K .Zhdanov and B.A .Trubnikov,Fiz. m orecom pletestudy assoon asexperimentaldatawill Plasmy 16,1457 (1990). be available.A properestimation of alsorequiresto 14. P. Blasi and D . de M arco, A stropart. Phys. quantifythebiasintroducede.g.by m isidentiedevents. 20, 559 (2004) [astro-ph/0307067]; M . K achelrie In general,itprovesm oree cient torem ove from the and D . Sem ikoz, A stropart. Phys. 23, 486 (2005) data set the doubtfulevents(e.g. in regions where a [astro-ph/0405258]. 676 given catalogue used to identify sources isknown to 15. A . C uoco et al., A strophys. J. , 807 (2008) be incom plete,orwhen severalsourcesatdierentdis- [0709.2712 [astro-ph]]. 124 tancesareidentied overa smallregion ofthe sky,with 16. A . M ucke et al., C om put.Phys.C om m un. , 290 (2000) [astro-ph/9903478]. possibleoverlap due tom agneticde ection orpooran- gular resolution),and applythe m ethod with a corre- 17. V .Berezinsky,A .Z.G azizov and S.I.G rigorieva,Phys. R ev.D 74,043005 (2006) [hep-ph/0204357]. spondinglysmaller statistics. Sources physicallyclus- 18. See also R ef.[17] and V . Berezinsky,A . G azizov and tered inthe universearenota problem here,sincethey M . K achelrie, Phys. R ev. Lett. 97, 231101 (2006) are located essentiallyat the sam e distance from the [astro-ph/0612247]. Earthand thussu erfrom the sam eattenuation during 19. D .H arari,S.M ollerach and E.R oulet,JCA P 0611,012 propagation. (2006) [astro-ph/0609294]. 20. O .E.K alashev etal.,arXiv:0710.1382 [astro-ph]. 21. V . Berezinsky, A . Z. G azizov and S. I. G rigorieva, REFERENCES astro-ph/0210095; Phys. Lett. B 612 (2005) 147 [astro-ph/0502550]. 22. D .A llard,E.Parizotand A .V .O linto,A stropart.Phys. 26 1. M .K achelrie and D .V .Sem ikoz,A stropart.Phys. , 27,61 (2007)[astro-ph/0512345]. 10 (2006) [astro-ph/0512498]. 23. M .A ve etal.,in Proc.\30thInternationalC osmicRay 2. S. M ollerach et al., to appear in Proc. M exico, 2007 , C onference",M erida,M exico,2007 ,# 0297. arXiv:0706.1749 [astro-ph]. 3. J. A braham et al. A stropart. Phys. 29, 188 (2008) [astro-ph/0712.2843]. 4. R . U . A bbasi et al. [HiRes C ollaboration],, arXiv:0804.0382 [astro-ph]; D . S. G orbunov et al.,