The impact of the Fermi Gamma-ray Space Telescope on blazar research is reviewed. This includes a brief description of the Fermi Large Area Telescope, a summary of the various classes of extragalactic sources found in the First Large Area Telescope AGN Catalog, and more detailed discussion of the flat spectrum radio quasar 3C454.3 and the BL Lac object PKS 2155-304. Some theoretical studies related to ongoing blazar research with Fermi are mentioned, including implications of -ray observations of radio galaxies on blazar unification scenarios, variability in colliding shells, and whether blazars are sources of ultra-high energy cosmic rays.
IAU 275. Blazar jet physics in the age of Fermi] Blazar jet physics in the age of Fermi Charles D. Dermer] Charles D. Dermer, on behalf of the Fermi LAT Collaboration 2010 \volume275 \pagerangexxx–yyy \jnameJets at all Scales \editorsG.E. Romero, R.A. Sunyaev and T. Belloni, eds.
Even without the Fermi Gamma-ray Space Telescope, the study of blazar jets would be rapidly advancing from discoveries made with ground-based -ray Cherenkov telescopes and AGILE. With Fermi, we have a telescope that is giving the most complete census yet of the MeV sky. Fermi makes sensitive measurements of the -ray spectra of blazars where the energy output often dominates. Unexpected puzzles have arisen. This includes the origin of the multi-GeV cutoffs in 3C 454.3 and other FSRQs. Another is the relation of blazars to radio galaxies detected at -ray energies. The blazar sequence ([Fossati et al.(1998), Fossati et al. 1998)], if it is that, remains sketchily understood, even with our new knowledge from Fermi.
Here we briefly review the results of blazar studies in the age of Fermi, including the 1 Fermi LAT AGN Catalog (1LAC; [Abdo et al. (2010a), Abdo et al. 2010a]) and the misaligned AGN (MAGN) sample of radio galaxies detected with Fermi. An overview of various classes of galaxies and AGNs detected with Fermi is given, including some detail on individual sources. 3C 454.3 and PKS 2155-304 Theoretical questions arising from these observations are raised, and some speculations in terms of colliding shells ejected from the central supermassive black hole are made. The likelihood that blazars accelerate ultra-high energy cosmic rays is considered.
2 The Large Area Telescope on Fermi
The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope is a pair conversion telescope that tracks the direction of an electron-positron pair formed by an incident ray. The converter-tracker is made of a 44 array of 16 modules. Each module has 18 tracker planes that uses Si strips to follow the direction of the converted pair. The upper 16 tracker planes have tungsten converters, with the uppermost 12 “front” planes using thin layers of tungsten to minimize multiple scattering and provide a narrow point spread function (psf), and the lower 4 “back” planes using thick tungsten layers to increase the effective area at the expense of the psf. The bottom two planes have no converters, only Si tracking planes. The energy is deposited in the calorimeter, which is made of 96 CsI(Tl) crystals. The entire assembly is surrounded by a segmented anti-coincidence detector used to reject charged particles. See [Atwood et al.(2009), Atwood et al. (2009)] for more detail.
The LAT has a field-of-view, roughly defined by the solid angle within which the effective area is more than one-half of the on-axis effective area, of 2.4 sr, which corresponds to of the full sky. For the most stringent diffuse event class used to identify rays, the on-axis effective area increases from cm at 100 MeV to cm at 1 GeV, and is roughly constant at higher energies. The energy resolution is better than 10% between MeV and GeV. For front-converting events, the 68% acceptance cone is and at photon energies MeV and GeV, respectively. Updated instrument performance that improves on pre-flight Monte Carlo and muon calibrations by using inflight corrections to the instrument response functions are found at http://www-glast.slac.stanford.edu/software/IS/glastlatperformance.htm.
Fermi nominally operates in a zenith-pointing rocking mode that allows the LAT to survey the entire sky every two orbits, or every 3 hours. It typically rocks 35 north and south of the zenith on alternate orbits. Autonomous repoints to gamma-ray bursts and two dedicated pointings, one to 3C 454.3 and a second to the flaring Crab, have been executed by Fermi.
3 LBAS and 1LAC
Two lists of AGNs detected with Fermi have now been published by the Fermi Collaboration. After 3 months of science observations between 2008 August 4 and 2008 October 30, the LAT Bright AGN Sample (LBAS) was released ([Abdo et al.(2009a), Abdo et al. 2009a]). It consists of 106 high-latitude () sources associated with AGNs. These sources have a test statistic , corresponding to significance, and are a subset of the 205 sources listed in the Fermi LAT bright source list ([Abdo et al.(2009b), Abdo et al. 2009b]). By comparison, the 3 Egret Catalog (3EG; [Hartman et al.(1999), Hartman et al. 1999]) of -ray sources and the revised EGRET catalog (EGR; [Casandjian & Grenier(2008), Casandjian & Grenier 2008]) list 31 sources with significance , of which 10 are at high latitude. Remarkably, 5 of the EGRET sources are not found in the LAT bright source list. These are the flaring blazars NRAO 190, NRAO 530, 1611+343, 1406-076, and PKS 1622-297, the most luminous EGRET blazar ([Mattox et al.(1997), Mattox et al. 1997)].
The 1LAC ([Abdo et al. (2010a), Abdo et al. 2010a]) is a subset of the 1451 sources in the First LAT Fermi Source Catalog (1FGL; [Abdo et al. (2010b), Abdo et al. 2010b]) derived from analysis of data taken during 11 months of observation between 2008 August 4 and 2009 July 4. There are 1043 1FGL sources at high latitudes, of which 671 are associated with 709 AGNs, with the larger number of AGNs than 1LAC sources due to multiple associations. Associations are made by comparing the localization contours with counterparts in various source catalogs, for example, the flat-spectrum 8.4 GHz CRATES (Combined Radio All-Sky Targeted Eight GHz Survey; [Healey et al.(2007), Healey et al. 2007]) and the Roma BZCAT blazar catalog ([Massaro et al.(2009), Massaro et al. 2009]). The probability of association is calculated by comparing the likelihood of chance associations with catalog sources if randomly distributed. Positional coincidence only affords an association; correlated flux variability between different wavebands is required for a firm identification.
|BL Lac objects||275||weak emission lines||AO 0235+164|
|…ISP||44||Hz Hz||3C 66A, W Comae|
|…HSP||114||Hz||PKS 2155-304, Mrk 501|
|FSRQs||248||strong emission lines||3C 279, 3C 354.3|
|…Starburst||3||active star formation||M82, NGC 253|
|…MAGN||7||steep radio spectrum AGNs||M87, Cen A, NGC 6251|
|…RL-NLS1s||4||strong FeII, narrow permitted lines||PMN J0948+0022|
|…NLRGs||4||narrow line radio galaxy||4C+15.05|
Total adds to 27, because the RL-NLS1 source PMN J0948+0022 is also classified as an FSRQ in the 1LAC
Includes PKS 0336-177, BZU J0645+6024, B3 0920+416, CRATES J1203+6031, CRATES J1640+1144, CGRaBS J1647+4950, B2 1722+40, 3C 407, and 4C +04.77
Of the 671 associations, 663 are considered “high-confidence” associations due to more secure positional coincidences. The “clean” sample is a subset of the high-confidence associations consisting of 599 AGNs with no multiple associations or other peculiarities. As listed in Table 1, these subdivide into 275 BL Lac objects, 248 flat spectrum radio quasars, 26 other AGNs, and 50 AGNs of unknown types. The “New Classes” category contains non-blazar AGNs, including starburst galaxies and radio galaxies of various types (e.g., narrow line and broad line). An AGN is classified as an “unknown” type either because it lacks an optical spectrum, or the optical spectrum has insufficient statistics to determine if it is a BL Lac objects or a FSRQ. In comparison with the 671 AGNs in the 1LAC taken with 11 months of Fermi data, EGRET found 66 high-confidence () detections of blazars out of total of 271 sources in the 3EG, with another 27 lower-confidence detections with significance between and . Thus the 1LAC already represents an order-of-magnitude increase in the number of AGNs over EGRET. There are – 400 unassociated and therefore unidentified high-latitude Fermi sources in the LBAS.
4 Classification of radio-emitting AGNs and unification
Different classes of extragalactic AGNs are defined according to observing frequency. We already noted the association of Fermi sources with BL Lac objects and flat spectrum radio quasars (FSRQs). This represents an optical classification. The precise definition used by the Fermi team is that an AGN is a BL Lac object if the equivalent width of the strongest optical emission line is Å, and the optical spectrum shows a Ca II H/K break ratio in order to ensure that the radiation is predominantly nonthermal (the Ca II break arises from old stars in elliptical galaxies). The wavelength coverage of the spectrum must satisfy in order that at least one strong emission line would have been detected if present. This helps guard against biasing the classification for AGNs at different redshifts where the emission lines could be redshifted out of the relevant wavelength range. For sources exhibiting BL Lac or FSRQ characteristics at different times, the criterion adopted is that if the optical spectrum conforms to BL Lac properties at any time, then it is classified as a BL Lac object.
The criterion for classification of radio galaxies according to their radio properties stems from the remarkable correlation between radio morphology and radio luminosity ([Fanaroff & Riley(1974), Fanaroff & Riley 1974)]. The twin-jet morphology of radio galaxies is seen in low-power radio galaxies, whereas the lobe and edge-brightened morphology is found in high-power radio galaxies, with a dividing line at W/Hz at 178 MHz, or at a bolometric radio luminosity of erg s. Besides a radio-morphology/radio-power classification, radio spectral hardness can also be used to characterize sources as flat-spectrum and steep-spectrum sources. The misaligned AGNs (MAGNs) are 1LAC sources associated with steep (, with ) radio-spectrum objects typically at 178 MHz in the Third Cambridge and Molonglo radio catalogs that show extended radio structures in radio maps. By contrast, some radio galaxies show hard radio spectra in the 1 – 10 GeV range, and can be subdivided according to the widths of the optical emission lines into broad- and narrow-line radio galaxies ([Perlman et al.(1998), Perlman et al. 1998]). Because of the close relation between core dominance and -ray spectral properties, core dominance can be used to infer the alignment of a radio-emitting AGN (e.g., [Lister & Homan (2005), Lister & Homan 2005)].
Blazars and radio galaxies can also be classified according to their broadband spectral energy distribution (SED) when there is sufficient multiwavelength coverage to reconstruct a spectrum from the radio through the optical and X-ray bands. When the peak frequency of the synchrotron component of the spectrum is Hz, then a source is called low synchrotron-peaked (LSP), whereas if the SED has Hz, then it is referred to as high synchrotron-peaked (HSP). Intermediate synchrotron-peaked (ISP) objects have Hz Hz. SEDs of the bright Fermi LBAS sources are constructed in [Abdo et al. (2010c), Abdo et al. (2010c)]. Essentially all FSRQs are LSP blazars, whereas BL Lac objects sample the LSP, ISP, and HSP range.
According to the standard unification scenario for radio-loud AGNs ([Urry & Padovani(1995), Urry & Padovani 1995)], radio galaxies are misaligned blazars, and FR1 and FR2 radio galaxies are the parent populations of BL Lac objects and FSRQs, respectively. To establish this relationship requires a census of the various classes of sources that takes into account the different beaming properties for the Doppler-boosted radiation of blazars. Even if analysis of data of radio galaxies and blazars supports the unification hypothesis, this paradigm still does not explain the reasons for the differences between radio-quiet and radio-loud AGNs, or between BL Lac objects and FSRQs.
Other classes of extragalactic Fermi sources found in the 1LAC include starburst galaxies, narrow line radio galaxies (NLRGs), radio-loud narrow line Sy 1s (RL-NLS1s), and radio-quiet AGNs. The recent -ray detections of the starburst galaxies M82, NGC 253, and NGC 4945 (though the latter source has a Sy 2 nucleus that might contribute to the -ray emission) with Fermi, VERITAS, and HESS opens up a new avenue of research in cosmic-ray studies. Five NLRGs are reported in the 1LAC. These objects have narrow emission lines in their optical spectrum, suggesting that they are observed at large angles with respect to the jet direction, with the surrounding dust torus obscuring the broad line region (BLR).
RL-NLS1s have also been recently established as a -ray source class ([Abdo et al. 2009c, Abdo et al. 2009c]). These objects show narrow H lines with FWHM line widths km s, weak forbidden lines () and a strong Fe II bump, and are therefore classified as narrow-line type I Seyferts ([?, Pogge 2000)]. By comparison with the black holes in blazars, the host galaxies of RL-NLS1s are spirals that host nuclear black holes with relatively small ( – ) mass that accrete at a high Eddington rate. The detection of these objects challenges scenarios where radio-loud AGNs are hosted by elliptical galaxies that form as a consequence of galaxy mergers.
The 1LAC includes 10 associations with radio-quiet AGNs. Radio-quiet AGNs have not yet been established as a -ray source class. In 8 of these cases, at least one blazar, radio galaxy, or CRATES source is also found close to the -ray source. In the remaining two cases, the association probabilities are weak. Thus none appear in the 1LAC “clean” sample. It remains unclear whether any radio-quiet sources, including Sy 2 galaxies such as NGC 4945 or NGC 1068 ([Lenain et al.(2010), Lenain et al. 2010]), produce rays from Sy nuclei rather than cosmic-ray processes.
5 Properties of Fermi AGNs
Various correlations are found by comparing -ray properties of Fermi AGNs according to their radio, optical, or SED classification. Probably the most pronounced correlation is between the -ray spectral index measured between 100 MeV and 100 GeV, and optical AGN type. FSRQs have significantly softer spectra than BL Lac objects, with for FSRQs and for BL Lac objects in the LBAS ([Abdo et al.(2009a), Abdo et al. 2009a]). The SED classification shows that the -ray spectral index progressively hardens from and 1.96 when the class varies from FSRQs to LSP-BL Lacs, ISP-BL Lacs, and HSP-BL Lacs, respectively.
Fermi data reveal complex -ray blazar spectra. All FSRQs and LSP-BL Lac objects, and most ISP blazars show breaks in the – 10 GeV range ([Abdo et al.(2010d), Abdo et al. 2010d)]. Such a break was apparent from the first observation of the bright blazar 3C 454.3 ([Abdo et al. (2009d), Abdo et al. 2009d)], and will be discussed in more detail below. This has the unfortunate consequence, however, to reduce the utility of the FSRQs for EBL studies. The HSP blazars, though, are generally well-described by a flat or rising SED in the GeV range, with spectral breaks between – 100 GeV implied from GeV measurements with air Cherenkov telescopes.
There is a much higher ratio, close to unity, of BL Lacs to FSRQs in the 1LAC compared to the 3EG, where the ratio was %. This is in large part due to the better multi-GeV sensitivity of Fermi over EGRET which finds far more HSP BL Lac objects than with EGRET, though stronger redshifting effects on the higher-redshift FSRQs than BL Lacs could also contribute. Only 121 out of 291 BL Lac objects had measured redshifts at the time of publication of the 1LAC. For sources with measured redshift , BL Lac objects are mostly found at low () redshifts, with only a few HSP BL Lac objects at higher redshifts. By contrast, the FSRQs span a wide range from to the highest redshift 1LAC blazar with .
This significant redshift incompleteness hampers interpretation of AGN properties, and the overall behavior and explanation of the blazar sequence (see discussion in Abdo et al. 2010a), and the role of RL NLS1s in understanding the sequence [Foschini et al.(2009), Foschini et al.(2009)], if it is not due to selection biases ([Padovani et al.(2003), Padovani et al. 2003)]). Nevertheless, a plot of as a function of apparent isotropic -ray luminosity in the 100 MeV – 10 GeV band for those sources with measured redshift can be constructed (Fig. 1; from [Abdo et al.(2010e), Abdo et al. 2010e]). It shows that the hard-spectrum BL Lac objects typically have much lower than the FSRQs. This divide has been interpreted as a change in the accretion regime at approximately 1% of the Eddington luminosity ([Ghisellini et al. 2009, Ghisellini et al. 2009]). In addition, the nearby radio galaxies with inhabit a separate portion of the vs. plane, and are characterized by lower -ray luminosities than their parent populations. The values of and of the two more distant steep spectrum radio sources, 3C 207 () and 3C 380 (), and the FR2 radio galaxy PKS 0943-76 () fall, however, within the range of -ray luminosities measured from FSRQs.
Both core and lobes can significantly contribute to the measured -ray luminosities. In the case of Centaurus A, the values of of the core and lobes are comparable, with the lobe emission primarily attributed to Compton-scattered CMBR ([?, Abdo et al. 2010f)]. The significant or dominant lobe component means that the core luminosity of misaligned AGNs can be less than the measured unless the lobe emission is resolved.
5.1 3c 454.3
The FSRQ 3C 454.3, at , underwent giant flares and became the brightest -ray source in the sky for a week in 2009 December and again in 2010 April ([Ackermann et al.(2010), Ackermann et al. 2010)]. The latter outburst triggered a pointed-mode observation by Fermi. During the December outburst, its daily flux reached ph( MeV) cm s, corresponding to an apparent isotropic luminosity of erg s, making it the most luminous blazar yet detected with Fermi. Using the measured flux and a one-day variability timescale at the time that the most energetic photon (with energy GeV) was detected, implies a minimum Doppler factor of . Assuming that the outflow Lorentz factor , consistent with the inferred value of and with radio observations at a different epoch ([Jorstad et al. 2005, Jorstad et al. 2005]), then simple arguments suggests a location pc, which is at the outer boundary of the BLR. Flux variability on time scales as short as 3 hr was measured at another bright flux state, which suggests that the -ray emission site would be even deeper in the BLR. This stands in contrast to inferences based on coherent optical polarization changes over timescales of weeks prior to a -ray flare in 3C 279 ([Abdo et al. (2010g), Abdo et al. 2010g)], which seems to place the emission site at much larger distances.
Figure 2 shows the light curve of 3C 454.3 taken from the public website on Fermi monitored sources
As noted already in the initial Fermi report ([Abdo et al. (2009d), Abdo et al. 2009d)], the spectrum of 3C 454.3 breaks strongly by units at GeV. Such a break is inconsistent with simple radiative cooling scenarios, which predict a break by 0.5 units. The more recent analysis of 3C 454.3 data ([Ackermann et al.(2010), Ackermann et al. 2010)], including the major outbursts, confirms the strong spectral break and finds that is very weakly dependent on the flux state, even when the flux changes by more than an order of magnitude. No obvious pattern is found in the spectral index/flux plane, as might be expected in simple acceleration and cooling scenarios ([Kirk et al.(1998), e.g., Kirk et al. 1998)].
The origin of the spectral break in 3C 454.3 bears on several important issues in FSRQs: the location of the -ray emission site; the source of soft target photons in Compton-scattering models; and the relation of FSRQs and BL Lacs in view of the disappearance of such breaks in ISP and HSP blazars. Such a break would be readily understood if the target field was sufficiently intense to attenuate the blazar radiation by absorption processes, but the intense line Ly radiation field at 10.2 eV observed in GALEX measurements of 3C 454.3 ([Bonnoli et al.(2010), Bonnoli et al. 2010)] implies GeV ([Reimer (2007), Reimer 2007)]. Photon attenuation deep within the BLR by He II recombination and Ly radiation with eV has been proposed ([Poutanen & Stern(2010), Poutanen & Stern 2010)], but the model lacks a consistent treatment of synchrotron and Compton-scattered blazar flare emission. A full spectral model for the SED of 3C 454.3 can fit the break with a distribution of nonthermal electrons that scatters accretion-disk and BLR radiation ([Finke & Dermer(2010), Finke & Dermer 2010)], but a solution insensitive to changes in the location of the emission site requires a BLR with a wind-like density.
The spectral break could also be due to Klein-Nishina effects in scattering, as has been proposed to explain the SED of PKS 1510-089 ([Abdo et al.(2010h), Abdo et al. 2010h)]. The KN break due to upscattered Ly radiation occurs at a few GeV, and the observed break energy is insensitive to the Doppler factor. But the break is not sharp enough to fit the spectrum of 3C 454.3 spectrum for a power-law electron distribution, and the electron spectrum resulting from cooling is likely to exhibit a hardening due to competition between synchrotron and Compton losses ([Dermer & Atoyan(2002), Dermer & Atoyan 2002,][Moderski et al. (2005), Moderski et al. 2005)].
5.2 PKS 2155-304, HSP BL Lac objects, and FR1 radio galaxies
PKS 2155-304, an X-ray selected BL Lac object at , is one of the most prominent representatives of the HSP blazar population. Its two-peaked SED was measured for 11 days between 2008 August 25 and 2008 September 6 at optical (ATOM), X-ray (RXTE and Swift), and -ray (Fermi and HESS) frequencies. This low-state SED is well fit by a one-zone synchrotron/SSC model with Doppler factor , magnetic field G, and comoving size cm (corresponding to a variability time of 2 d; [Aharonian et al.(2009), Aharonian et al. 2009]). During a period of extraordinary flaring on 2006 July 28, PKS 2155-304 exhibited a sequence of -ray flares varying on time scales as short as min ([Aharonian et al.(2007), Aharonian et al. 2007)]. Detailed one-zone synchrotron SSC model fits using a low EBL require in order to avoid attenuation of the TeV rays and provide a good fit to the broadband SED ([Finke et al.(2008), Finke et al. 2008)].
One-zone synchrotron/SSC models with give good fits to other HSP BL Lac objects such as Mrk 421 and Mrk 501. The SEDs of radio galaxies, in contrast, are fit with much smaller Doppler factors. The SED of the core of Cen A, for instance, can be fit with and bulk Lorentz factors few ([Abdo et al.(2010i), Abdo et al. 2010i)]. Likewise, the SEDs of the FR1 radio galaxies NGC 1275 ([Abdo et al.(2009e), Abdo et al. 2009e]) and M87 ([Abdo et al. (2009f), Abdo et al. 2009f]) are well fit with and .
The much larger values of for BL Lac objects than for their putative parent population, the misaligned FR1 radio galaxies, is contrary to simple unification expectations. Moreover, even though the -ray luminosities from FR1 radio galaxies are much smaller than that of BL Lac objects (Fig. 1), they are still larger than expected by debeaming the radiation of BL Lac objects with . Additional soft target photons that can be Compton scattered to high energies result in a reduction of the value of compared to those implied by the one-zone synchrotron/SSC model. These target photons can be produced in a structured jet, as in the spine and sheath model ([Chiaberge et al.(2000), Chiaberge et al. 2000)]. Another soft photon source arises if blazar flows decelerate from the inner jet to the pc scale ([Georganopoulos & Kazanas(2003), Georganopoulos & Kazanas 2003]), in accord with the mildly relativistic flows at the sub-pc scale found in radio observations of Mrk 421 and Mrk 501.
Given the results of synchrotron/SSC modeling, the Fermi observations suggest that the -ray emission from the core of a radio galaxy originates from a slower region than the emission of BL Lac objects. This can be understood in a colliding shell model if shells with large opening angles have lower and are less energetic than shells with narrow opening angles. Such a circumstance might also explain the short variability timescale of PKS 2155-304, as we now show.
Consider a simple colliding shell model where the second shell is much more powerful and has a much larger Lorentz factor than the first shell. Denote the Lorentz factors of the first and second shells by , with . The collision radius takes place at , where , and ejection time scales are required to be greater than /c, where is the Schwarzschild radius. If the second shell is much more powerful than the first, then a strong forward shock is formed that travels through shell with Lorentz factor given by the relative Lorentz factor of the two shells, so that (see [Sari & Piran(1995), Sari & Piran 1995)].
For short timescale variability, both the radial and angular timescales must be much shorter than . The radial timescale
for low-redshift sources, noting that is the width of shell in the frame of shell , and shell widths in the engine frame are required to be . Here the Doppler factor , where is the shocked fluid Lorentz factor. The angular timescale
In both cases, the variability timescale can be much shorter than when .
The criterion for a strong forward shock is that the ratio of apparent luminosities . This may seem an extreme requirement if the opening angles of the slow and fast shells are the same. But if , , a condition arising in simulations of relativistic jets (e.g., [Komissarov et al.(2009), Komissarov et al. 2009], and references therein), then excessive energy requirements for the fast shell are mitigated. The maximum radiative efficiency is ([Dermer & Razzaque (2010), Dermer & Razzaque 2010]), so this system is necessarily very inefficient. This suggestion may not only resolve the short variability measured in PKS 2155-304, but also, as noted above, possibly relieve the Doppler-factor conflict in unification schemes of blazars and radio galaxies. A more detailed treatment is in preparation.
6 Blazars and ultra-high energy cosmic rays
In Fermi acceleration scenarios, two conditions are required to accelerate particles to ultra-high energies eV. The first is that the isotropic power must exceed erg s, where is the shocked fluid Lorentz factor, and is the charge. The second is that the sources have luminosity density erg Mpc yr within the GZK radius to power UHECRs against photohadronic losses on the CMBR. Consider the FR1 radio galaxy NGC 1275, at Mpc. Its -ray luminosity from the 1LAC is erg s. For the mildly relativistic outflow speeds from synchrotron/SSC models, it has adequate power to accelerate Fe nuclei. Moreover, it is radiating nonthermal -ray power within this volume with a luminosity density erg Mpc yr. If a small fraction of this power is channeled into UHECRs, then FR1 radio galaxies and BL Lac objects are favored to be the sources of UHECRs, provided that UHECRs are high-Z ions ([Dermer & Razzaque (2010), Dermer & Razzaque 2010]).
The Fermi LAT is providing a uniform spectral and temporal GeV database with much better sensitivity than either the EGRET or AGILE missions. Besides BL Lac objects and Flat Spectrum Radio Quasars, several new classes of -ray galaxies have been established, including radio galaxies, radio-loud narrow line Sy 1 galaxies, and star-forming galaxies. The -ray spectral slope of blazars is strongly correlated with -ray luminosity and whether they are BL Lac objects or FSRQs, though a large fraction of BL Lac objects still lack redshift measurements. FSRQs exhibit spectral cutoffs between 1 – 10 GeV, the reason for which is not well understood. Gamma-ray emitting FR1 radio galaxies have much lower -ray luminosities than BL Lac objects, but not as low as expected for one-zone synchrotron/SSC models. Colliding shells with a range of collision Lorentz factors might account for the different -ray luminosities of radio galaxies and blazars, and the short variability timescale measured in PKS 2155-304. BL Lacs and FR1 galaxies have sufficient power and luminosity density to account for the UHECRs.
Acknowledgements.This work is supported by the Office of Naval Research and NASA. I thank L. Foschini for comments and T. Piran for criticism. The LAT Collaboration acknowledges support from a number of agencies and institutes for both development and the operation of the LAT as well as scientific data analysis. These include NASA and DOE in the United States, CEA/Irfu and IN2P3/CNRS in France, ASI and INFN in Italy, MEXT, KEK, and JAXA in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the National Space Board in Sweden. Additional support from INAF in Italy and CNES in France for science analysis during the operations phase is also gratefully acknowledged.
David MeierEarly in your talk you mentioned 10 possible RQ objects that have been detected, but I did not see them on your vs. luminosity plot. Where do they lie on that plot?
Chuck DermerThe association of radio-quiet AGN with -ray sources is very tentative, as we discuss in the 1LAC paper. Moreover, none of them survive in the “clean” sample shown in the plot. Because the detection of -ray emission from radio-quiet AGNs could also originate from cosmic-ray processes, as in the case of starburst galaxies, and so might require the detection of -ray flux variability to confirm, the Fermi Collaboration is not yet ready to make a definitive statement. The analysis is ongoing.
Tsvi PiranIn internal shocks the time scale, because of angular spreading, is the variability time of the inner engine. I don’t see how you can avoid this in your model.
Chuck Dermer The collision radius is set by the speed of the slower shell. The angular spreading is set by the speed of the shocked fluid which, for a strong forward shock, can be much greater than the speed of the slower shell. The smaller opening angle of the high speed, large luminosity shell also reduces the angular spreading time scale. I agree that this process cannot be very efficient, but if % of the Eddington luminosity is channeled into a narrow opening angle jet for during bright flares, only % need be converted to yield erg s from PKS 2155-304.
Elisabete M. de Gouveia Dal Pino Regarding the production of UHECRs by FR1/BL Lac sources, surely they have the “power” to produce them, but how do the CRs escape from these compact sources?
Chuck Dermer The escape of UHECRs without suffering photodisintegration in the AGN radiation fields poses an important constraint on the location of UHECR acceleration sites. In FR1 radio galaxies and BL Lac objects, the broad line region is absent, and the accretion-disk radiation is weak. Whether the host galaxy radiation field poses a challenge to this model will require more detailed calculations.
- Aharonian, F., et al. 2007, ApJL, 664, L71
- Aharonian, F., et al. 2009, ApJL, 696, L150
- Abdo, A. A., et al. 2009a, ApJ, 700, 597 (LBAS)
- Abdo, A. A., et al. 2009b, ApJS, 183, 46 (0FGL)
- Abdo, A. A., et al. 2009c, ApJL, 707, L142 (RL-NLS1s)
- Abdo, A. A., et al. 2009d, ApJ, 699, 817 (3C 454.3)
- Abdo, A. A., et al. 2009e, ApJ, 699, 31 (NGC 1275)
- Abdo, A. A., et al. 2009f, ApJ, 707, 55 (M87)
- Abdo, A. A., et al. 2010a, ApJ, 715, 429 (1LAC)
- Abdo, A. A., et al. 2010b, ApJS, 188, 405 (1FGL)
- Abdo, A. A., et al. 2010c, ApJ, 716, 30 (LBAS broadband SED)
- Abdo, A. A., et al. 2010d, ApJ, 710, 1271 (LBAS GeV SED)
- Abdo, A. A., et al. 2010e, ApJ, 720, 912 (MAGN)
- Abdo, A. A., et al. 2010f, Science, 328, 725 (Cen A lobes)
- Abdo, A. A., et al. 2010g, Nature, 463, 919 (3C 379)
- Abdo, A. A., et al. 2010h, ApJ, 721, 1425 (PKS 1510-089)
- Abdo, A. A., et al. 2010i, ApJ, 719, 1433 (Cen A core)
- Ackermann, M., et al. 2010, ApJ, 721, 1383 (3C 454.3 flares)
- Atwood, W. B., et al. 2009, ApJ, 697, 1071
- Bonnoli, G., Ghisellini, G., Foschini, L., Tavecchio, F., & Ghirlanda, G. 2010, MNRAS, in press (arXiv:1003.3476)
- Casandjian, J.-M., & Grenier, I. A. 2008, A&A, 489, 849
- Chiaberge, M., Celotti, A., Capetti, A., & Ghisellini, G. 2000, A&A, 358, 104
- Dermer, C. D., & Atoyan, A. M. 2002, ApJL, 568, L81
- Dermer, C. D., & Razzaque, S. 2010, ApJ, 724, 1366
- Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P
- Finke, J. D., & Dermer, C. D. 2010, ApJL, 714, L303
- Finke, J. D., Dermer, C. D., Böttcher, M. 2008, ApJ, 686, 181
- Foschini, L., for the Fermi/LAT Collaboration, Ghisellini, G., Maraschi, L., Tavecchio, F., & Angelakis, E. 2009, arXiv:0908.3313
- Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 1998, MNRAS, 299, 433
- Georganopoulos, M., & Kazanas, D. 2003, ApJL, 594, L27
- Ghisellini, G., Maraschi, L., & Tavecchio, F. 2009, MNRAS, 396, L105
- Healey, S. E., Romani, R. W., Taylor, G. B., Sadler, E. M., Ricci, R., Murphy, T., Ulvestad, J. S., & Winn, J. N. 2007, ApJS, 171, 61
- Hartman, R. C., et al. 1999, ApJS, 123, 79
- Jorstad, S. G., et al. 2005, AJ, 130, 1418
- Kirk, J. G., Rieger, F. M., & Mastichiadis, A. 1998, A&A, 333, 452
- Komissarov, S. S., Vlahakis, N., Königl, A., & Barkov, M. V. 2009, MNRAS, 394, 1182
- Lenain, J.-P., Ricci, C., Türler, M., Dorner, D., & Walter, R. 2010, arXiv:1008.5164
- Lister, M. L., & Homan, D. C. 2005, AJ, 130, 1389
- Massaro, E., Giommi, P., Leto, C., Marchegiani, P., Maselli, A., Perri, M., Piranomonte, S., & Sclavi, S. 2009, A&A, 495, 691
- Mattox, J. R., Wagner, S. J., Malkan, M., McGlynn, T. A., Schachter, J. F., Grove, J. E., Johnson, W. N., & Kurfess, J. D. 1997, ApJ, 476, 692
- Moderski, R., Sikora, M., Coppi, P. S., & Aharonian, F. 2005, MNRAS, 363, 954
- Padovani, P., Perlman, E. S., Landt, H., Giommi, P., & Perri, M. 2003, ApJ, 588, 128
- Perlman, E. S., Padovani, P., Giommi, P., Sambruna, R., Jones, L. R., Tzioumis, A., & Reynolds, J. 1998, AJ, 115, 1253
- Pogge, R. W. 2000, New Astronomy Reviews, 44, 381
- Poutanen, J., & Stern, B. 2010, ApJL, 717, L118
- Reimer, A. 2007, ApJ, 665, 1023
- Sari, R., & Piran, T. 1995, ApJL, 455, L143
- Urry, C. M., & Padovani, P. 1995, PASP, 107, 803