a multi - messenger approach is necessary to study the most energetic particles in nature .
we expect the sources to produce cosmic rays ( i.e. charged particles ) , photons and neutrinos .
their detection requires multiple techniques .
cosmic rays are detected in ground arrays and atmospheric fluorescence telescopes ( i.e. the pierre auger observatory ) , detectors in balloons ( e.g. heat , act ) and satellites ( e.g. pamela , fermi space telescope ) and the space station ( ams ) .
high energy gamma rays are studied with air cherenkov telescopes ( act , such as hess , veritas , cangaroo , magic ) , detectors in satellites ( e.g. integral and fermi ) , ground arrays and atmospheric fluorescence telescopes ( auger ) .
high energy neutrinos require under ice or water kilometer cube detectors , such as icecube in antarctica or the km3net in the mediterranean see , and radio telescopes in balloons , such as anita .
i will concentrate here on cosmic rays and will mention photons ( see the talk of francis halzen on neutrinos in this conference ) .
cosmic rays have a non - thermal spectrum well approximated by a power law close to @xmath0 for 10 orders of magnitude in energy , from @xmath1 ev to above @xmath2 ev , and 30 orders of magnitude in the flux .
the flux goes from 1 particle/(m@xmath3 sec ) at 10@xmath4 ev , to 1 particle/(m@xmath3yr ) at 10@xmath5 ev , to 1 particle/(km@xmath3yr ) at 10@xmath6 ev ( see fig .  1 ) .
the non - thermal spectrum of cosmic rays indicates that their acceleration results from stochastic processes in the presence of magnetic fields , as first proposed by fermi .
the diffusive acceleration of charged particles in the shock waves of supernova explosions may give origin to most cosmic rays .
but the highest energies of cosmic rays are so extreme that it is not understood if a similar mechanism can accelerate these particles even in the most extreme sources , which may or may not be large enough and have sufficiently large magnetic fields .
each energy range of cosmic rays addresses different physics . from 10@xmath7 ev to 10@xmath8 ev the modulation of the incoming charged particles by the solar wind can be studied .
the expanding plasma generated by the sun decelerates and partially sweeps - off the lower energy galactic cosmic rays from the vicinity of the earth .
this effect depends on the level of solar activity , which has an 11 year cycle . from 10@xmath8 ev to 10@xmath9 or 10@xmath10 ev , cosmic rays allow us to study the galactic sources that produce them and the propagation of charged particles within the galaxy . above this energy range and up to
at least a few 10@xmath11 ev the cosmic rays must be of extragalactic origin although the characteristics of the transition is still not understood .
thus the so called ultra - high energy cosmic rays ( uhecr ) , those above 10@xmath10 ev , allow us to study their extragalactic sources and the propagation of charged particles in the extragalactic medium .
cosmic ray spectrum @xmath0 over 10 orders of magnitude from @xmath1 ev to above @xmath2 ev(from ref .
the pamela satellite ( payload for antimatter exploration and light - nuclei astrophysics ) , consists of a magnetic spectrometer in orbit , launched in june 2006 , and has produced its first data releases in 2008 .
of particular interest has been the positron fraction @xmath13 data the pamela collaboration announced first in july 2008  @xcite ( see fig .  1 ) , which shows an excess in the 10 -100 gev range with respect to what is expected from secondary cosmic rays , i.e. a fraction diminishing as the energy increases ( secondary electrons and some positrons are produced in the interactions of cosmic ray protons and heavier nuclei ) .
this result confirms what was earlier called theheat excess "  @xcite , also observed by ams-01
@xcite . fig .
1 shows the pamela data ( in red ) and the data of several prior experiments ( in black ) .
the large dispersion in the data seen at energies below a few gev is expected , because the different experiments were done in different moments of the solar cycle . above 10 gev
the solar modulation is not important and this is the range where the excess is significant . as seen in fig .  1
, the positron fraction is about 0.05 at 10 gev and grows steadily with energy to reach 0.1 - 0.2 at 100 gev .    soon after the publication of the pamela positron data , in november 2008 , the atic collaboration announced a 6@xmath14 excess in the @xmath15 flux in the range 300 to 800 gev  @xcite , with respect to what is expected from secondary cosmic rays .
if the pamela positron fraction excess is due to an excess of positrons ( as opposed to a defect of electrons , which are much more abundant than positrons in secondary cosmic rays ) it may mean that positrons and electrons are been created by an unknown source .
the excess electron flux in the pamela energy range would not be noticed until it becomes comparable to , and then larger than , the secondary electron flux as the energy increases .
if this is so , atic may have measured the continuation to higher energies of the excess found by pamela .
energy spectrum of @xmath15 multiplied by @xmath16 measured by hess  @xcite , atic  @xcite and other experiments .
the double arrow shaded area indicate the approximate range of the systematic error in the hess spctrum ( from ref .
@xcite ) ]     energy spectrum of @xmath15 multiplied by @xmath16 measured by hess  @xcite , atic  @xcite and other experiments .
the double arrow shaded area indicate the approximate range of the systematic error in the hess spctrum ( from ref .
@xcite ) ]    atic ( the advanced thin ionization calorimeter instrument ) is a balloon born calorimeter launched from the mcmurdo base in antarctica , which has flown several times .
the data presented are from the 2000 - 2001 and 2002 - 2003 flights .
a few days after the announcement of atic , the hess collaboration published its @xmath15 flux measurement above 600 gev , which is steeply falling with increasing energy  @xcite .
fig .  2 ( from ref .
@xcite ) shows the data of hess , atic and other earlier experiments .
the @xmath15 flux is multiplied by @xmath16 .
the double arrow shows the possible shift of the hess spectrum due to @xmath17% energy scale uncertainty and the shaded area indicates the approximate range of the systematic error in the hess spectrum due to other reasons .
the hess data show that , if the atic data are correct , there is a sharp cutoff of any excess above about 800 gev , which is compatible with the sharp cutoff expected from dark matter particle annihilation at the parent particle mass .
since their appearance , at the moment of writing these notes the data of pamela and atic have generated more than 50 papers , at the rate of several per week , providing different explanations for their excesses .
the main two explanations are the existence of an astrophysical source nearby , or dark matter particles annihilating nearby , within a volume around earth much smaller than the volume of the galaxy .
electron and positrons rapidly loose energy through synchrotron and inverse compton processes , thus at the energies observed they must come from less than 1 kpc away ( for comparison , the distance between the sun and the center of the galaxy is about 8 kpc ) .
known possible astrophysical sources are pulsars , of which there are some , such as geminga and vela , within a few 100 pc from earth  @xcite .
it is not known the mechanism through which they would produce such energetic electrons and positrons , thus the spectrum they could produce is unknown too .
the maximum energy cutoff of the spectrum observed would in this case correspond to the maximum energy of emission of positrons and electrons .    in the case of annihilating dark matter particles
, the cutoff in energy would correspond to the mass of the parent particle .
the cutoff could be very sharp if the particles annihilate directly into @xmath18 pairs , and this would be a clear indication of the dark matter origin of the signal  @xcite ( see fig .
if instead @xmath18 are produced after a chain of decays of annihilation products , the cutoff would be less sharp , similar to what would be expected from astrophysical sources .
fig .  4 ( taken from ref .
@xcite ) shows three fits to the atic data , assuming an astrophysical source , a dark matter particle annihilating directly to @xmath18 and another producing @xmath18 as secondaries of the annihilation products .    in order to account for the pamela and/or atic data , the annihilation rate of most dark matter candidates needs to be enhanced by a factor , called
 boost factor " @xmath19 , which may reach @xmath20 in the case of many usual wimp dark matter candidates  @xcite and may be much smaller , close to @xmath21 ( see fig .
4 ) , for some candidates that decay preferentially into electrons and positrons ( the atic collaboration finds that a kaluza - klein particle decaying mostly into electrons and positrons with mass 620 gev fits well their electron - positron data  @xcite ) .
examples of fits to the atic data assuming a pulsar nearby ( red ) , the annihilation of 800 gev dark matter particles producing @xmath18 as secondaries with @xmath22 ( blue ) and the annihilation of 620 gev kaluza - klein dark matter particle decaying directly into @xmath18 with @xmath23 ( black ) . from ref .
@xcite ]    several ways of enhancing the local annihilation rate have been proposed : there could be a dense dark matter clamp nearby ( however , the probability found in structure formation simulations for this to happen is about 10@xmath24  @xcite ) , which would enhance the rate , which depends on the square of the dark matter density , or the annihilation cross section itself could be enhanced . for example
a  sommerfeld resonance " results from the almost formation of a bound state by the annihilating particles at very small relative velocities .
this effect would enhance the annihilation cross section in the dark halo ( but not in the early universe , where the particles move at larger velocities )  @xcite .
whatever produces the positron excess seen by pamela should not produce an excess of antiprotons , since the antiproton to proton ratio @xmath25 was also measured by pamela  @xcite and is compatible with what is expected from secondary cosmic rays .
this could be partially explained because there are many more protons than electrons in the secondary cosmic rays ( pamela requires a 10@xmath26 discrimination of protons to measure positrons ) and because @xmath27 and @xmath28 reach us from a larger volume than @xmath29 and @xmath30 , from a large fraction of the total volume of the galaxy .
the fermi space telescope ( ex glast - gamma - ray large area space telescope ) can also detect @xmath29 and @xmath30 ( it does discriminate between the two ) , besides photons .
fermi was launched in june 2008 and is providing @xmath31-ray spectroscopic data of unprecedented quality .
it has two detectors ( the glast burst monitor , gmb , and the large area telescope , lat ) which detect photons between 10 kev and 300 gev , and observe the whole sky every 3 hours .
photons of even higher energy will be detected by air cherenkov telescopes ( acts ) , cangaroo iii in australia , hess in namibia , magic in las palmas and veritas in the us ( see the talk of manel martinez in this conference ) .
these can also detect @xmath32 , as we already mentioned in the case of hess .
the universe is totally transparent to photons below 100 gev . at higher energies photons
interact with infra - red and optical backgrounds but still at energies below 10 s of tev arrive to us from cosmological distances ( at higher energies they are absorbed by the cosmic microwave background radiation ) .
photons can again reach earth from cosmological distance at energies above @xmath1 gev , the range of uhecr .
the photons observed by fermi reveal , therefore , the spacial distribution of their sources .
in particular , if fermi observes the annihilation of dark matter particles into photons , it is expected to detect not only the center of our galaxy but also a large amount of dark matter clumps . as mentioned above , since the annihilation rate depends on the square of the density , high density regions boost the rate . a large amount of dark matter clumps
are expected to remain within the dark halo of our galaxy .
haloes grow hierarchically , incorporating lumps and tidal streams from earlier phases of structure formation .
the best simulations of structure formation at present ( the via lactea ii  @xcite and the aquarius project  @xcite ) show that 10 s of subhaloes could be discovered by fermi .
however , subhaloes are more effectively destroyed by tidal effects near the center of the galaxy , thus most of them are far from the sun . the chance that a random point close to the sun is lying within a clump is smaller than 10@xmath24  @xcite .
dwarf galaxies and other objects which fermi might observe as a whole through dark matter annihilation into photons will have the annihilation rate increased by a boost factor @xmath19 , may be as large of 100 , due to the existence of substructure within them ( with respect to the rate they would have without dark matter lumps ) . but
the concept of a boost factor can not be applied to the possible enhancement of a signal in electron and positrons as seen by pamela or atic , since at most only one or very few dark matter clump could be present within the volume these particles come from .
the fermi space telescope can measure the @xmath15 flux in the range explored by atic , possibly to 1tev , with better accuracy and its results are expected to be made public soon , in the april 2009 meeting of the american physical society ( which will take place in the first week of may 2009 ) .
all - particle energy spectrum of cosmic rays , the flux is multiplied by @xmath16 ( from ref .
the lines represent spectra for elemental groups @xmath34 as indicated according to a particular model .
the sum of all elements ( galactic ) and a presumably extragalactic component are shown as well ( see ref .
@xcite for details ) . ]    in fig .
1 two features of the cosmic ray spectrum are indicated : the  knee " above 10@xmath35 ev and the  ankle " above 10@xmath10 ev .
these are shown also in fig .  5 ( taken from ref .
@xcite ) in which the flux has been multiplied by a factor @xmath16 . at the knee
there is a slight steepening of the spectrum , the spectral index changes from @xmath36 at low energies to about -3.1 at higher energies .
second knee " feature is sometimes mentioned , at about 4@xmath37 ev , where the spectrum seems to become again steeper at higher energies with a spectral index approximately @xmath38 ( see fig .
5 ) . at the ankle
the spectrum seems to flatten again to a spectral index of about @xmath39 .
the  ankle " is sometimes seen as just part of a  dip " in the spectrum between 10@xmath10 and 10@xmath6 ev .    at energies below 10@xmath40 ev the elemental composition of cosmic rays
is well known because cosmic ray particles are detected directly .
the composition and flux of cosmic rays at these energies are in agreement with their production in galactic supernova explosions  @xcite .
direct measurements above the atmosphere to energies slightly above 10@xmath40 ev are performed by various instruments , such as atic , cream , bess and tracer . at higher energies ,
the flux of cosmic rays becomes smaller than several particles /(m@xmath3 day ) so their detection with instruments carried in balloons or spacecrafts becomes impractical and cosmic rays are detected from the ground through the shower of secondary particles they produced in the atmosphere .
these are called extensive air showers ( eas ) .
arrays of detectors covering a large area , usually from a fraction to many km@xmath3 , operate for years to obtain enough statistics .
because of the sparse sampling of the shower , the uncertainties in the modeling of the shower development and the large statistical fluctuations in the interactions in the atmosphere , the inferred energy and characteristics of the shower can only determine the relative contribution of groups of elements to the cosmic ray flux .
for example , a heavy nucleus ( usually referred to just as  iron " ) generates a shower with a relatively higher ratio of muons to electrons that a proton initiated shower of the same energy .
a good understanding of the composition around the knee would help determining what is the origin of this feature , but this is difficult .
the best measurements of the composition of cosmic rays in the knee region come from the kascade ( karlsruhe shower core and array detector ) experiment , which observes a steepening of the spectrum first of light elements and successively at higher energies for heavier ones  @xcite ( see fig .  5 ) .
the knee can be explained by the cutoff of the galactic proton spectrum due to the leakage of protons from the galaxy .
the galactic magnetic fields are not strong enough to keep protons of this energy bound to the galaxy .
heavier elements are bound to the galaxy up to higher energies , so their cutoff energy increases approximately as their charge @xmath34 .
however , the data may also be explained if nuclei are accelerated in the supernova shock fronts to maximum energies proportional to @xmath34  @xcite .
two possibilities for the crossing of galactic to extragalactic cosmic rays : 3a .
( left ) shows the crossing at the  second knee " @xmath41 gev ( the ankle is part of a dip due to the interaction of extragalactic protons with the cmb ) ; 3b .
( right ) shows the crossing at the ankle , above 10@xmath42 gev .
figure taken from ref .
@xcite ]    the interpretation of the second knee is linked to that of the ankle .
the second knee is sometimes explained by the fall off of the heaviest elements ( including elements heavier than iron ) in the galactic cosmic rays  @xcite ( see fig .
the second knee would then be the place in the spectrum where an extragalactic component becomes dominant , as shown in fig .
compilation by r. engel of the main cosmic rays data up to 2007 .
the fluxes are multiplied by @xmath43 .
the agasa spectrum does not show a gzk feature , which is present in the hires and 2007 auger spectra . ]    in this case , the ankle is an absorption  dip " in the spectrum between 10@xmath10 and 10@xmath6 ev proposed to be caused by electron - positron pair production of extragalactic protons on cosmic microwave background ( cmb ) photons over large distances  @xcite .
this is one of the possible explanations for the ankle .
historically the ankle was interpreted as the transition from a rapidly falling galactic ( iron dominated ) component to a flatter spectrum of extragalactic origin subdominant at lower energies .
both interpretations are shown in fig .
6b  @xcite .
the cosmic rays above 10@xmath10 ev are called ultra - high energy cosmic rays ( uhecr ) .
pierre auger observed cosmic rays with energies above 10@xmath35 ev already in the 1930 s , and 80 years later their origin is still uncertain .
he was first in understanding that extensive air showers are produced by high energy cosmic rays .
giuseppe cocconi argued in 1956 that cosmic rays above the ankle are extragalactic in origin because they can not be confined by the galactic magnetic field . in 1962
john linsley at volcano ranch in new mexico observed the first 10@xmath11 ev shower .    immediately after the discovery of the microwave background radiation by penzias and wilson in 1965 , in two separate papers in 1966 , greisen , and independently zatsepin and kuzmin  @xcite pointed out that if the uhecr are of extragalactic origin there should be an absorption feature in their spectrum at about 4@xmath44 ev , the so called  gzk cut - off " . in fact , if the uhecr are extragalactic their energy is degraded by inelastic collisions with background photons as they propagate .
protons would loose energy mainly through photopion production in the cmb background @xmath45 , where the pion carries away @xmath46 of the original nucleon energy .
the mean free path for this reaction is only @xmath47  mpc .
heavier nuclei mostly loose energy through photodisintegration .
for both , protons and heavy nuclei , the energy losses become important above about 4@xmath44 ev , thus they could not reach earth from a distance beyond 50 to 100 mpc .
hillas plot " showing the magnetic field versus size of potential uhecr sources .
acceleration of protons to 1@xmath48zev@xmath49 ev , or protons or fe nuclei to 100 eev @xmath50 ev require conditions above the respective line . ]    until recently the existence of the gzk absorption feature was in question . the results of the main uhecr experiments until a few years ago , agasa and hires , were contradictory ( see fig.7 ) .
the agasa experiment in japan , an array of 111 surface detectors ( scintillator counters ) distributed in 100 km@xmath3 with 1 km spacing between them , took data from 1984 to 2003 and found no gzk cutoff  @xcite .
this gave rise to several models to explain the  agasa excess " above 10@xmath11 ev , in which the uhecr were produced locally so the cosmological energy absorption would not be important .
these so called  top - down " models invoked new physics to produce uhecr directly with the high energies required , as opposed to  bottom - up " models in which uhecr are accelerated in astrophysical sites .    in 1993 , the fly s eye experiment , consisting of a fluorescence telescope in utah , usa , detected a 3@xmath5110@xmath11 ev event .
the upgrade of this experiment , the high resolution fly s eye ( hires ) , with two fluorescence telescopes took data between 1997 and 2006 and found a spectrum compatible with the gzk cutoff  @xcite .    the pierre auger observatory in argentina  @xcite , a hybrid experiment which combines the two detection methods of agasa and hires , was designed to prove or disprove the existence of the gzk cutoff and to elucidate the mystery of the origin of the highest energy cosmic rays .
it consists of both surface detectors ( water tanks ) , sd , and fluorescence detectors , fd .
the spectrum of agasa showing no gzk feature , and those of hires and auger in 2007 compatible with a gzk feature , are shown in fig .
top - down models were introduced not only to produce the highest energy cosmic rays locally and avoid the gzk absorption feature , but also as an alternative to acceleration models to explain the highest energy cosmic rays , which the latter models have difficulty explaining .
only the most extreme sources in the universe are possible sites for the acceleration to protons or other nuclei at those energies : active galactic nuclei ( agn ) , gamma ray bursts ( grb ) , magnetars ( i.e. neutron stars with enormous magnetic fields @xmath52 g ) , huge lobes and hot spots of radio galaxies , super nova remnants ( snr ) , colliding galaxies , large shocks in galaxy clusters ( for a review see for example ref .
@xcite , from which fig .  8 was taken ) .
the maximum energy at which particles can be accelerated in astrophysical sources is limited by the acceleration rate and the energy loss rate due to inverse compton scattering or inelastic scattering .
the  hillas plot "  @xcite ( see fig .
8) shows astrophysical sites for which the gyroradius @xmath53 of a particle of charge @xmath34 and the desired energy @xmath54 , in the estimated magnetic field @xmath19 of the site , is smaller than the size @xmath55 of the site .
the idea is that only in sites where the particles can be confined for a while , the acceleration could take place .
8 shows a modern version of the hillas plot  @xcite , including many more potential production sites than the original version of the plot ( from 1984  @xcite ) . as the figure shows , the possible sites where protons could be accelerated diminishes fast as the energy increases above 10@xmath11 ev .
top - down models circumvented the need for acceleration of particles by producing them at very high energies . in super heavy dark matter ( shdm )
models , super heavy metastable particles were produced in early universe and remained at present as part of the dark matter of the universe .
these particles ( with colorful names such as ` cryptons ' or ` wimpzillas ' ) could decay or annihilate  @xcite into quarks and leptons etc which in turn produced the observed uhecr .
these would come primordially from within the dark halo of our galaxy , thus shdm models predicted an excess of events towards the galactic center , which has not been found . in these models
the maximum energy of cosmic rays was given by the parent particle mass .
in other models , the parent particles instead of being metastable were instead produced by topological defects , such as cosmic strings or necklaces ( for a review see ref .
@xcite ) . in the  z - burst " model  @xcite
, ultra - high energy neutrinos from remote sources would annihilate at the z - resonance with relic background neutrinos near earth ( this model encountered problems with the large flux of ultrahigh energy neutrinos required ) .
the z bosons would then decay , producing secondary protons , neutrinos and photons .
the z - resonance , which is the uhecr energy cutoff , occurs when the energy of the incoming @xmath56 is @xmath57 .
the spectra of the uhecr produced in top - down models are determined by the elementary particle physics of z - boson decays and of qcd fragmentation , which predict photon domination of the spectrum at high energies .
9a ( left panel ) fluxes measured by yakutsk , agasa , hires and auger , multiplied by @xmath16 .
9b ( righ panel ) the  dip " has been used to calibrate the energy of the experiments by factors within their systematic errors .
the spectra now coincide , except at the highest energies .
see ref .
@xcite for details . ]
the pierre auger observatory consists of 1600 water tanks with 1.5 km spacing , distributed in a 3,000 km@xmath3 surface , the surface detector ( sd ) and four fluorescence detector sites , each with 6 telescopes that overlook the surface array , the fluorescence detector ( fd ) .
the tanks detect the cherenkov light produced by the electromagnetic and muon component of the air showers .
the telescopes detect the near ultra - violate light from the flourescent emission of nitrogen induced by the shower .
auger started operating in 2004 and has been completed a few months ago , in june 2008 .
the auger exposure now is of the order of 10@xmath58 km@xmath3 sr yr .
the fd measures the longitudinal shower distribution and provides a calorimetric measurements of the shower energy , but can operate only in moonless nights , which gives it a 10% duty cycle . since the fd telescopes observe the light emitted by the air molecules and not the energy deposited by the shower particles in the air , the energy measurement inferred with the fd depends sensitively on the value adopted for the fluorescence yield , which is still not well known .
this is an important source of systematic error .
for example , auger and hires use values of the fluorescence yield that differ by about 10%
. a better determination of the fluorescence yield is expected soon .
the sd measures the transverse shower distribution and observes the sky all the time , thus it provides most of the statistics at the highest energies .
about 15% of the showers observed by auger are hybrid events , observed both by the sd and the fd .
the hybrid events are used to calibrate the energy of the much larger number of sd only events .
the systematic uncertainly in the energy determination of auger is at present 22% .
the fluxes measured by yakutsk , agasa , hires and auger are displayed in fig .
9a ( left panel ) .
the multiplication of the fluxes by @xmath16 exacerbates the importance of the systematic errors in the energy calibration .
fig.9b ( right panel ) shows how the spectra coincide , except at the highest energies , if the dip feature that appears in all of them is used to calibrate the energy by factors within the systematic errors of the different experiments .
the energy calibration factors used to produce fig .
9b are 1.0 , 1.2 , 0.75 , 0.83 and 0.625 for hires , auger , agasa , akeno and yakutsk respectively ( see ref .
@xcite for details ) .
the pierre auger observatory was optimized to observe in the energy range between 10@xmath10 ev and 10@xmath11 ev , because this is the window in which cosmic rays could reach us from extragalactic sources .
the extragalactic magnetic fields are not expected to be larger than 10@xmath59 g. recent realistic simulations of the expected large scale extragalactic magnetic fields , show that except in the regions close to the virgo , perseus and coma clusters the obtained magnetic fields are not larger than 3@xmath60  g  @xcite ( however see ref .
we can see then from eq .  1 that with @xmath61 g at energies @xmath62 ev protons would only be weakly deflected over distances of mpc . at energies above 10@xmath11 ev the gzk cutoff
will considerably reduce the flux of uhecr .
thus there is a window between 10@xmath10 ev and 10@xmath11 ev to possibly do extragalactic astronomy with charged particles .
the main properties of uhecr auger is expected to determine are four : 1 ) the spectrum of the highest energy cosmic rays , in particular to prove or disprove the existence of a gzk cutoff ; 2 ) the composition of those high energy primaries ( are they protons , heavy nuclei , photons ? ) ; 3 ) the arrival direction distribution ( is there a large scale anisotropy , a correlation with particular sources ? ) ; 4 ) the sources . in the following we address these issues in turn
.     fits to the auger spectrum with different compositions and injection spectra of primary uhecr  @xcite . ]
fits to the auger spectrum with different compositions and injection spectra of primary uhecr  @xcite . ]
the about 2@xmath63 events above 2.5@xmath64 ev collected by auger up to august 2007 , with about 4 times the exposure of agasa and 2 times that of hires , have a spectrum compatible with a gzk cutoff at 4@xmath44 ev , as shown in fig .
10  @xcite .
the lower panel of fig .
10 shows the fractional differences between the hires and auger data at the highest energies and a power law @xmath65 which fits the spectrum below 4@xmath44 ev .
the number of events expected if this power law spectrum were to continue to energies above 4@xmath44 ev or above 1@xmath66 ev would be 167@xmath673 and 35@xmath671 respectively , while auger observed 69 and 1 respectively .
this amounts to a 6@xmath14 significance of the flux suppression at the gzk energy , which is compatible with what was found by hires .
the spectrum can not be used to determine the composition of uhecr , since models with pure proton primaries and mixed proton - iron primaries all produce good fits to the spectrum .
fig.11 shows several fits to the spectrum of auger with proton primaries and mixed composition primaries and several primary spectral indices @xmath31 and maximum energies @xmath68 characterizing the spectrum emitted by the sources  @xcite . even the extreme unrealistic assumption of having pure iron emitted at the sources provides good fits to the auger spectrum  @xcite .
it is also not possible either to determine from the spectrum if the observed suppression is due to the gzk effect or to the energy cut of the spectra of the cosmic rays emitted at the sources .
illustration of the observation of the longitudinal shower development with the auger fd .
the maximum defines @xmath69 . from ref .
@xcite ]    the main tool to study the mass composition of cosmic rays is the measurement of the longitudinal profile of the showers .
12 shows an illustration of the measurement of the longitudinal shower development observed with the fluorescence detectors of auger . for a given energy ,
heavy nuclei interact earlier in the atmosphere than protons ( nuclei interact as a collection of nucleons ) .
photon initiated showers are much more penetrating in the atmosphere than hadrons , up to energies higher than 10@xmath11 ev at which photons produce a  preshower " ( through interaction with the magnetic field of the earth ) before entering into the atmosphere .
thus the column depth of air at which the shower profile is maximum , @xmath69 , as measured with the fd , is well correlated with the primary particle mass .
this can be seen in fig .
13 where the expected @xmath69 as function of energy for photon , proton and iron initiated showers are shown , as predicted by several simulation programs .
the value of the @xmath69 for protons is about 100 g @xmath70 larger than for iron .
the fd detector of auger can measure @xmath69 with a systematic uncertainty of about 15 g / cm@xmath3 .
13 also shows a compilation of the earlier data on @xmath69 as function of the energy , for energies above 10@xmath40 ev .
it shows a composition compatible with proton dominance at the highest energies . except for one data point of hires ,
the measurements reached 2@xmath44 ev .
compilation of earlier measurements of the mean depth of shower maximum @xmath69 as function of energy , compared with the predictions of various simulations programs for photon , proton and iron initiated showers ( see for example ref .
@xcite ) . ]
auger data on @xmath69  @xcite are presented in fig .
14 together with the predictions for proton and iron initiated showers from several shower simulation programs . although with low statistics the auger data show a change of regime ( a change in slope , called  elongation rate " ) at 2@xmath44 ev towards a heavier composition . at the highest energies
the composition is intermediate between protons and iron , compatible with a mean mass number of about 5 .
the errors shown in the figure are statistical .
the systematic uncertainty is @xmath67 15 g / cm@xmath3 .
the hires collaboration did not find the change in slope at 2@xmath44 ev present in the auger data , and its composition measurements seem compatible with proton dominance to the highest energies ( see fig.16 )
.     mean depth of shower maximum @xmath69 as function of energy measured by auger  @xcite compared with what is expected from proton and iron initiated showers .
the last data point seems to indicate a heavier composition above 2@xmath44 ev . ]
mean depth of shower maximum @xmath69 as function of energy measured by hires  @xcite ( fig . from ref .
@xcite ) compared with what is expected from proton and iron initiated showers .
no change in the slope is observed at 2@xmath44 ev . ]
auger has obtained stringent upper bounds on photons as uhecr  @xcite .
uhecr photons with energy above 10@xmath6ev are one of the key observables to distinguish top - down from bottom - up models ( see for example ref .
because in the process of hadronization pions are produced more copiously than protons , and neutral pions decay into photons , all top - down models produce more photons than protons ( and no heavy nuclei ) . also
gzk " photons , those produced in the decay of photo - produced @xmath72 in the gzk process must always be present at some level , although the level may be extremely low
. the photon fraction could be as small at 10@xmath24  @xcite .
finding gzk photons among the uhecr would help understanding the primary hadronic spectrum which produced them and the intervening extragalactic backgrounds with which they interact , the radio background and magnetic fields .
the flux of gzk photons is related to that of  cosmogenic " neutrinos , those neutrinos produced in the decay of the photo - produced charged pions in the gzk process .
neutrinos and photons produced as secondaries of the gzk process would give complementary information since neutrinos come without absorption from the whole production volume , while photons can reach us only from less than 100 mpc .
photons with energy @xmath73 ev pair - produce electrons and positrons on the radio background and can not reach earth from beyond 10 to 40 mpc ( although the photon energy - attenuation length is not well known , due to the uncertainties in the spectrum of the absorbing radio background ) .
same as fig .  13 , but with several characteristic top - down models and other experimental upper bounds .
see ref .
@xcite for details . ]
same as fig .  13 , but with several characteristic top - down models and other experimental upper bounds .
see ref .
@xcite for details . ]    at energies below 10@xmath11 ev photons produce in the atmosphere deeper showers than protons ( due to the landau - pomeranchuk - migdal ( lpm ) effect ) . at higher energies
the photons interact with the magnetic field of the earth before entering the atmosphere and produce several electrons below the lpm threshold ( what is called  preshowering " ) but the statistics of auger is too low at those energies .
bounds based on @xmath69 measurements using fd data were initially used to set an upper bound on photons  @xcite .
better bounds were obtained later  @xcite by using two parameters , the shower - front radius and the rise time of the signal , that can be measured with the sd .
late developing showers have smaller shower - front radius ( the shower - front flattens out as the shower develops ) .
photon initiated showers have a lower muon content , more low energy particles , and thus longer rise - time of the signal . with these two parameters ,
the high statistics sd could be used to search for photons .
auger found no events compatible with what is expected from a photon primary . as a result
it has obtained stringent upper bounds shown in figs .
16 and 17 .
the range of photon fractions in the integrated flux for some top - down models and for gzk photons from ref .
@xcite together with the 2006 and 2007 auger and agasa - yakutzk ( ay )  @xcite upper bounds are shown in fig .
. the auger bounds together with other upper bounds on the photon fraction and characteristic top - down models are shown in fig .
17  @xcite .
these photon bounds provide an argument independent of possible correlation with astrophysical sources against most top - down models .
if the suppression of the cosmic rays spectrum at 4@xmath44 ev is due to the gzk effect , namely to the interactions with the cmb photons , then the cosmic rays with higher energies should come from nearby sources .
for example , above 6@xmath44 ev 90% of the protons should come from less than about 200 mpc away , while 50% should come from less than 90 mpc away .
thus , the arrival directions of the highest energy cosmic rays should correlate with the distribution of visible matter nearby , which is very inhomogeneus .
the angular resolution of auger above 1@xmath44 ev is less than one degree , while the deflection of protons or nuclei in the @xmath74 g magnetic field of our galaxy , at those energies is expected to be about 10@xmath75 .
auger has , in fact , found a correlation with extragalactic astrophysical sources  @xcite .
auger reported a correlation between the arrival directions of uhecr with @xmath76 ev and the positions of nearby active galactic nuclei ( agn ) with redshift @xmath77 ( i.e. distance @xmath7871 mpc ) : 20 out of 27 cosmic ray events were found to correlate with at least one of 442 agn selected form a particular catalog , the veron - veron - cetty catalog ( 292 within the field of view of auger ) located within a circle of 3.2@xmath79 around the event arrival direction ( only 5.6 are expected to do so if the flux were isotropic ) .
the auger collaboration found a statistically significant correlation with a smaller amount of data in may 2006 with very similar energy threshold and angular distance . fixing these parameter a - priori it studied the subsequent data until august 2007 when the original correlation was confirmed to a pre - determined confidence level in the additional data set alone .
map in galactic coordinates of the positions of the agn within 71 mpc used in the correlation search of auger ( red stars ) and circles of 3.2@xmath79 centered at the arrival directions of the 27 events with @xmath76 ev observed by auger until august 2007 .
the shading indicate levels of exposure  @xcite . ]
the map of the arrival directions and of the agn positions is shown in fig .
there is a clear alignment of several of the correlating uhecr events close to the supergalactic plane , indicated in the figure with a dashed line .
two of the events have arrival directions less than 3@xmath79 away from the nucleus of centaurus a , one of the closest agn , at 3.4 mpc away .
cen a displays jets and radio lobes which extend over a scale of about 10@xmath79 along the supergalactic plane , and a variable compact radio nucleus .
the two events mentioned correlate with the nucleus position while several lie in the vecinity of the radio lobe extension along the supergalactic plane . in general , agn with prominent radio lobes are rare and do not follow the observed spatial distribution .
so far there are only three other doublet events ( i.e. the same agn associated with two uhecr ) , one of which is close to and other two which are outside the supergalactic plane
.    the most prominent radio galaxy within 70 mpc , m87 in the virgo cluster , does not correlate with any observed event above 5.7 @xmath44 ev .
there is a lack of events in the region of the virgo cluster which seems to be present also in the data of the hires collaboration
. may be large magnetic fields in the direction of virgo could explain the lack of events in that direction .
the hires collaboration did not confirm a correlation of uhecr with agns  @xcite ) : out of 13 events with @xmath80 only two correlated with an agn within 3.2@xmath79 ( and 3.2 events were expected just by chance ) .
hires observed mostly the northern sky and auger the southern , could they be significantly different ?    given the small statistics of the sample and the negative results of hires , the correlation found by auger does not prove that the sources of uhecr are agns .
the sources may be any astrophysical objects with a space distribution similar to agns .
are agn the sources or just tracers of the sources ?
nuclei are attenuated by photodisintegration in the intergalactic medium and intermediate mass nuclei ( a @xmath81 ) would have to reach us from distances shorter than 70 mpc .
heavier nuclei , close to fe have horizons comparable to that of protons but they would be deflected in the magnetic field of our galaxy by more than 3@xmath79 .
thus the correlation indicates that the primary uhecr are most likely protons .
however , in the case the primaries are protons , the actual sources of the uhecr may be further out than 70 mpc .
we mentioned above that the  gzk horizon " ( defined as the distance from earth which contains the sources that produce 90% of the protons which arrive with energies above a given threshold ) should be about 200 mpc for protons of energy above 6 @xmath44 gev , while it would be 70 mpc for about 9 @xmath44 gev .
this may indicate that the energy calibration of auger should be shifted upwards by 30% or that there are accidental correlations with agns that are closer than the actual uhecr sources .
the correlation with agns found by auger seems to indicate that the highest energy cosmic rays are protons , while the composition data of auger suggest a change towards a heavier composition above @xmath82 ev . which of these two pieces of information is correct
? there may also be a problem with the programs used to simulates the showers and determine if they are proton or iron like .
data from lhcf are needed to improve them .
clearly more data are needed to answer all these questions , and they are coming .
there is an infield with extra muon counters in construction in the pierre auger observatory in argentina .
also a new experiment , the telescope array ( with a low energy extension called tale ) is being built in utah ( usa ) .
this auger observatory in argentina is the south part of the two - part project initially proposed for auger .
there is a current proposal for an auger - north observatory in colorado ( usa ) with 4@xmath83 surface detectors occupying an area of 2@xmath63 km@xmath3 .
there is also the jem - euso mission in preparation , the extreme universe space observatory that will be hosted in the japanese experimental module of the international space station , that will study very energetic cosmic ray showers from space .
many advances have occurred in recent years in our understanding of the origin and propagation of cosmic rays but many fundamental open question remain . in the gev to tev energy range , the solar modulation of cosmic rays and the solar wind need to be better understood , and we need to understand if the pamela/ atic excess is in fact there and it is not an instrumental issue , and if so if it is due to a pulsar nearby or the dark matter annihilation .
the fermi space telescope data expected soon will greatly help in this respect . in the 10@xmath58 - 10@xmath7 gev energy range , the cosmic ray composition , the galactic sources , the propagation of cosmic rays within the galaxy and the transition from the galactic to extra - galactic cosmic rays ( at the second knee or at the ankle ? ) need to be better understood . at energies above 10@xmath42 gev charged particle extragalactic astronomy
has started ! but we need yet to identify the sources , understand the composition of the cosmic rays , the mechanisms of acceleration at the sources and the intervening backgrounds and magnetic fields .
9              aharonian f a , atoyan a and volk h j 1995 _ astron .
astrophys . _ * 294 * l41 ; hooper d , blasi p and serpico pd 2008 _ jcap _ * 0901 * 025 [ arxiv:0810.1527 [ astro - ph ] ] ; yuksel h , kistler m d and stanev t 2008 arxiv:0810.2784 [ astro - ph ] ; profumo s 2008 arxiv:0812.4457 [ astro - ph ] .                  berezinsky v , gazizov az , grigorieva s i 2004 _ nucl .
suppl . _ * 136 * 147 [ astro - ph/0410650 ] ; berezinsky v , gazizov az , grigorieva s i 2005 _ phys .
_ b * 612 * 147 [ astro - ph/0502550 ] ; berezinsky v , gazizov az , grigorieva s i 2006 _ phys .
_ d * 74 * 043005                  berezinsky v , kachelriess m and vilenkin a 1997 _ phys .
lett . _   * 79 * 4302 ; kuzmin v a and rubakov v a 1998 _ phys .  atom
.  nucl . _   * 61 * 1028 [ _ yad .  fiz . _   * 61 * 1122 ] ; birkel m and sarkar s 1998 _ astropart .  phys . _
* 9 * 297 ; blasi p , dick r and kolb e w 2002 _ astropart .
phys . _   * 18 * 57 ; aloisio r , berezinsky v and kachelriess m 2004 _ phys .  rev . _
d * 69 * 094023 .
dolag k , grasso d , springel v and tkachev i 2004 _ jetp lett .
_ * 79 * , 583 [ _ pisma zh .  eksp .
fiz . _ * 79 * , 719 ] [ arxiv : astro - ph/0310902 ] ; dolag k , grasso d , springel v and tkachev i 2005 _ jcap _ * 0501 * 009 [ arxiv : astro - ph/0410419 ] .