Arxiv:0811.1848V3 [Astro-Ph] 18 Dec 2008 Cieglci Uli Ore O Lr Iheeg Cosmi Energy High Biermann Ultra L
Total Page:16
File Type:pdf, Size:1020Kb
Active Galactic Nuclei: Sources for ultra high energy cosmic rays? Peter L. Biermannabcde, Julia K. Beckerf g, Laurent¸iu Carametea h, Alex Curut¸iua, Ralph Engele, Heino Falckei j, L´aszl´o A.´ Gergelyk l, P. Gina Isare h, Ioana C. Mari¸se, Athina Melim, Karl-Heinz Kampertn, Todor Stanevo, Oana Ta¸sc˘aun, Christian Ziera p aMPI for Radioastronomy, Bonn, Germany bDept. of Phys. & Astron., Univ. of Bonn, Germany cDept. of Phys. & Astr., Univ. of Alabama, Tuscaloosa, AL, USA dDept. of Phys., Univ. of Alabama at Huntsville, AL, USA eInst. Nucl. Phys. FZ, Karlsruhe Inst. of Techn. (KIT), Germany f Institution f¨or Fysik, G¨oteborgs Univ., Sweden gDept. of Phys., Univ. Dortmund, Dortmund, Germany hInstitute for Space Studies, Bucharest, Romania iDept. of Astrophys., IMAP, Radboud Univ., Nijmegen, Netherlands jASTRON, Dwingeloo, Netherlands kDept. Appl. Sci., London South Bank University, UK lDept. of Theoret. & Exp. Phys., Univ. of Szeged, Szeged, Hungary mPhysik. Inst. Univ. Erlangen-N¨urnberg, Germany nPhys. Dept., Univ. Wuppertal, Germany oBartol Research Inst., Univ. of Delaware, Newark, DE, USA arXiv:0811.1848v3 [astro-ph] 18 Dec 2008 pRaman Res. Inst., Bangalore, India The origin of ultra high energy cosmic rays promises to lead us to a deeper understanding of the structure of matter. This is possible through the study of particle collisions at center-of-mass energies in interactions far larger than anything possible with the Large Hadron Collider, albeit at the substantial cost of no control over the sources and interaction sites. For the extreme energies we have to identify and understand the sources first, before trying to use them as physics laboratories. Here we describe the current stage of this exploration. The most promising contenders as sources are radio galaxies and gamma ray bursts. The sky distribution of observed events yields a hint favoring radio galaxies. Key in this quest are the intergalactic and galactic magnetic fields, whose strength and structure are not yet fully understood. Current data and statistics do not yet allow a final judgement. We outline how we may progress in the near future. 1. Introduction In the quest of understanding the fundamental structure of matter the physics community has 1 2 built the Large Hadron Collider (LHC) at CERN in Geneva, which will have particle energies up into the TeV region. In the universe, as we know from direct observation (first Linsley 1963), we have particles up to 300 EeV (= 3 1020 eV); even in the center of mass collision with· an identical particle (proton or heavier) this yields energies in the center of mass frame of 5 1014 eV, so more. If we could identify the sources· and interaction regions of these extreme energy particles (see the books by Ginzburg & Syrovatskii 1964, Berezin- sky et al. 1990, Gaisser 1991, Stanev 2004, and recent reviews by Gaisser & Stanev 2006, Bier- mann et al. 2003, 2006, as well as Nagano & Watson 2000) we may be able to learn some of Figure 1. This overlay spectrum shows the public the physics at such energies, so perhaps go be- data, as of summer 2008, of AGASA, HiRes and yond the LHC. Auger (Abbasi et al. (HiRes) 2008; Auger 2008b). In this short review we will discuss the latest The AGASA event energy has been scaled down trends in the quest to understand where the ex- by 25 %, and the HiRes event energy has been tremely high energy particles come from, and how scaled down by 10 %. The spectra are adjusted we might be able to test our ideas. For lack of to show a common flux near 1019 eV. The main space we only give a small fraction of all refer- problem here is the energy estimate of the differ- ences. ent detectors. 2. Source candidates While very many ideas exist based on de- ies. On the other hand, the arrival directions do tailed physical models for possible sources of ul- seem compatible with the nearest radio galaxy, tra high energy cosmic ray particles, the best Cen A = NGC5128, a source long suspected to bet candidates to explain them are radio galaxies emit cosmic rays (Ginzburg & Syrovatski 1963). (Ginzburg & Syrovatskii, Blandford, Biermann, However, this radio galaxy has so little power, it et multi al.) and gamma ray bursts (M´esz´aros, presents a special challenge to understand how it Piran, Rees, Vietri, Waxman, et multi al., with a could accelerate particles to 1020 eV, and beyond. recent summary, e.g., Waxman 2006). As gamma For another nearby radio galaxy, M87 in the ray bursts are special cases of very massive star Virgo cluster (Ginzburg & Syrovatski 1963, Cun- explosions, their occurrence should correlate with ningham et al. 1980), the synchrotron spectrum galaxies, which have a current starburst, so are of the knots in the jet has been used to argue that strong in the far infrared, such as, e.g., M82 (e.g. it requires protons at 1021 eV to initiate the cas- Kronberg et al. 1985), NGC253, NGC2146 and cade in the plasma for scattering the non-thermal the like; the early models (Biermann 1976, Bier- electrons in order to yield a parameter-free cut- mann & Fricke 1977) already allowed the predic- off frequency of near 3 1014 Hz, as observed in tion of far-infrared from radio fluxes from star- many knots, hot spots and· nuclei (e.g. Rieke et burst galaxies, such as NGC2146 (Kronberg & al. 1976) of radio emitting active galactic nuclei Biermann 1981), and so gave a prediction of the (Biermann & Strittmatter 1987). This is in fact supernova rate (today perhaps equivalent to a the only argument based on observations which prediction of the gamma ray burst rate). At implies the existence of these ultra high energy present the statistics of the arrival directions do particles in the source. However, it has to be not support a correlation with starburst galax- noted, that this does not necessarily imply that 3 the particles we observe come from such sources; and it is just plausible for lack of many alternatives. F S2 7 (2) CR ∼ . ,tot 2.1. Complete samples This flux corresponds to distance attenuation In order to test the idea that radio galaxies are −2 source candidates, we have developed the jet-disk only with DL . Interestingly, the mass of the symbiosis concept (papers by Falcke et al. 1995a, black hole does not even enter here due to the b, Markoff et al. 2001, Yuan et al. 2002, Massi simplicity of the Poynting flux argument. Follow- & Kaufman Bernardo 2008, etc.). Therefore we ing the argument below we might have to multi- need a complete sample of steep spectrum radio ply the maximum particle energy by 6 - 8 or so sources (e.g. teams led by Witzel, see K¨uhr et al. to simulate the seeding with heavier nuclei from 1981). Table 1 presents such a complete list, dif- a weak starburst, indicated by a relatively large ferentiated in two sets in redshift range; the com- FIR/radio ratio as given in Table 1; this ratio is plete sample takes all extragalactic steep spec- still far below that of a pure starburst, for which trum sources down to a flux density limit, which it is of order 300. However, the numbers given are not already known as predominantly starburst for the maximal energy do not include this extra galaxies from their far-infrared/radio flux density factor. In Table 2 this is indicated by an asterisk. ratio, and that are within the redshift specified. Accretion powered jets are the other alterna- Table 2 extends the list to slightly higher redshift. tive, which works well for relatively high current accretion rates (Falcke et al. 1995a, b, Ta¸sc˘au 2.2. Particle energy and particle flux pre- 2004, Ta¸sc˘au et al. 2008): dictions E† S1/3 D2/3 M (3) The main indefinite parameter in the jet-disk max ∼ rad L BH symbiosis picture is the anchoring of the mag- netic field at the base of the jet. Spin-down and powered jets emanating from very near super- † 2/3 −2/3 F S D (4) massive black holes (Blandford & Znajek 1977, CR ∼ rad L Blandford & Koenigl 1979, Boldt & Ghosh 1999) For distances < 50 Mpc usually NGC5128 = are one possibility to do this: The jet power is Cen A, possibly NGC1316 = For A, and a group roughly proportional to the total radio luminos- around M87 = Vir A dominate in predicted ity (Enßlin et al. 1997), especially if we include UHECR flux (Ginzburg & Syrovatskii 1963). low power sources in the crude fit. If we identify The first five in flux density of the extended the jet power as an upper limit to the Poynting flux are ESO137-G006, NGC1316, NGC4261, flux, and use the relationship between Poynting NGC4486=M87, and NGC5128. flux and maximal particle energy (Lovelace 1976, An early attempt to fit older data is shown in see below), we obtain an expression for the maxi- Fig. 2. mal particle energy. Furthermore we can assume These two approaches allow to understand the that the cosmic ray flux is a fraction of the total huge range in radio to optical flux ratios from jet power, and so obtain a simple proportional- active galactic nuclei (Strittmatter & Witzel et ity. Calling the mass of the black hole MBH , the al., 1980), and the ubiquity of low flux densities observed compact radio flux density Srad or ex- of compact radio emission from basically all early tended total flux density at 2.7 GHz S2.7,tot, the Hubble type galaxies (e.g. Perez-Fournon & Bier- luminosity distance DL to the radio galaxy, the mann 1984).