The magnetic field in intergalactic space gives important information about magnetogenesis in the early universe. The properties of this field can be probed by searching for radiation of secondary e$^+ $e$^-$ pairs created by TeV photons, that produce GeV range radiation by Compton-scattering cosmic microwave background (CMB) photons. The arrival times of the GeV "echo" photons depend strongly on the magnetic field strength and coherence length. A Monte Carlo code that accurately treats pair creation is developed to simulate the spectrum and time-dependence of the echo radiation. The extrapolation of the spectrum of powerful gamma-ray bursts (GRBs) like GRB 130427A to TeV energies is used to demonstrate how the IGMF can be constrained if it falls in the $10^{-21}$ - $10^{-17}$ G range for 1 Mpc coherence length.

The {\gamma}-ray flares from the Crab nebula observed by AGILE and Fermi-LAT reaching GeV energies and lasting several days challenge the standard models for particle acceleration in pulsar wind nebulae, because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Previous modeling has suggested that the synchrotron limit can be exceeded if the electrons experience electrostatic acceleration, but the resulting spectra do not agree very well with the data. As a result, there are still some important unanswered questions about the detailed particle acceleration and emission processes occurring during the flares. We revisit the problem using a new analytical approach based on an electron transport equation that includes terms describing electrostatic acceleration, stochastic wave-particle acceleration, shock acceleration, synchrotron losses, and particle escape. An exact solution is obtained for the electron distribution, which is used to compute the associated {\gamma}-ray synchrotron spectrum. We find that in our model the {\gamma}-ray flares are mainly powered by electrostatic acceleration, but the contributions from stochastic and shock acceleration play an important role in producing the observed spectral shapes. Our model can reproduce the spectra of all the Fermi-LAT and AGILE flares from the Crab nebula, using magnetic field strengths in agreement with the multi-wavelength observational constraints. We also compute the spectrum and duration of the synchrotron afterglow created by the accelerated electrons, after they escape into the region on the downstream side of the pulsar wind termination shock. The afterglow is expected to fade over a maximum period of about three weeks after the {\gamma}-ray flare.

Knowledge of the distant high-energy universe comes from photons, ultra-high energy cosmic rays (UHECRs), high-energy neutrinos, and gravitational waves. The theory of high-energy messengers reviewed here focuses on the extragalactic background light at all wavelengths, cosmic rays and magnetic fields in intergalactic space, and neutrinos of extragalactic origin. Comparisons are drawn between the intensities of photons and UHECRs in intergalactic space, and the high-energy neutrinos recently detected with IceCube at about the Waxman-Bahcall flux. Source candidates for UHECRs and high-energy neutrinos are reviewed, focusing on star-forming and radio-loud active galaxies. HAWC and Advanced LIGO are just underway, with much anticipation.

Active Galactic Nuclei can be copious extragalactic emitters of MeV-GeV-TeV gamma rays, a phenomenon linked to the presence of relativistic jets powered by a super-massive black hole in the center of the host galaxy. Most of gamma-ray emitting active galactic nuclei, with more than 1500 known at GeV energies, and more than 60 at TeV energies, are called "blazars". The standard blazar paradigm features a jet of relativistic magnetized plasma ejected from the neighborhood of a spinning and accreting super-massive black hole, close to the observer direction. Two classes of blazars are distinguished from observations: the flat-spectrum radio-quasar class (FSRQ) is characterized by strong external radiation fields, emission of broad optical lines, and dust tori. The BL Lac class (from the name of one of its members, BL Lacertae) corresponds to weaker advection-dominated flows with gamma-ray spectra dominated by the inverse Compton effect on synchrotron photons. This paradigm has been very successful for modeling the broadband spectral energy distributions of blazars. However, many fundamental issues remain, including the role of hadronic processes and the rapid variability of those BL Lac objects whose synchrotron spectrum peaks at UV or X-ray frequencies. A class of gamma-ray--emitting radio galaxies, which are thought to be the misaligned counterparts of blazars, has emerged from the results of the Fermi-Large Area Telescope and of ground-based Cherenkov telescopes. Blazars and their misaligned ounterparts make up most of the >100 MeV extragalactic gamma ray background (EGB), and are uspected of being the sources of ultra-high energy cosmic rays. The future "Cherenkov Telescope Array", in synergy with the Fermi-Large Area Telescope and a wide range of telescopes in space and on he ground, will write the next chapter of blazar physics.

Anisotropy in the arrival direction distribution of ultrahigh-energy cosmic rays (UHECRs) produced by powerful sources is numerically evaluated. We show that, taking account of the Galactic magnetic field, nondetection of significant anisotropy at $\approx 10^{19}$ eV at present and in future experiments imposes general upper limits on UHECR proton luminosity of steady sources as a function of source redshifts. The upper limits constrain the existence of typical steady sources in the local universe and limit the local density of $10^{19}$ eV UHECR sources to be $\gtrsim 10^{-3}$ Mpc$^{-3},$ assuming average intergalactic magnetic fields less than $10^{-9}$ G. This isotropy, which is stronger than measured at the highest energies, may indicate the transient generation of UHECRs. Our anisotropy calculations are applied for extreme high-frequency-peaked BL Lac objects 1ES 0229+200, 1ES 1101-232, and 1ES 0347-121, to test the UHECR-induced cascade model, in which beamed UHECR protons generate TeV radiation in transit from sources. While the magnetic-field structure surrounding the sources affects the required absolute cosmic-ray luminosity of the blazars, the magnetic-field structure surrounding the Milky Way directly affects the observed anisotropy. If both of the magnetic fields are weak enough, significant UHECR anisotropy from these blazars should be detectable by the Pierre Auger Observatory unless the maximum energy of UHECR protons is well below $10^{19}$ eV. Furthermore, if these are the sources of UHECRs above $10^{19}$ eV, a local magnetic structure surrounding the Milky Way is needed to explain the observed isotropy at $\sim 10^{19}$ eV, which may be incompatible with large magnetic structures around all galaxies for the UHECR-induced cascade model to work with reasonable jet powers.

We search for evidence of dark matter (DM) annihilation in the isotropic gamma-ray background (IGRB) measured with 50 months of Fermi Large Area Telescope (LAT) observations. An improved theoretical description of the cosmological DM annihilation signal, based on two complementary techniques and assuming generic weakly interacting massive particle (WIMP) properties, renders more precise predictions compared to previous work. More specifically, we estimate the cosmologically-induced gamma-ray intensity to have an uncertainty of a factor ~20 in canonical setups. We consistently include both the Galactic and extragalactic signals under the same theoretical framework, and study the impact of the former on the IGRB spectrum derivation. We find no evidence for a DM signal and we set limits on the DM-induced isotropic gamma-ray signal. Our limits are competitive for DM particle masses up to tens of TeV and, indeed, are the strongest limits derived from Fermi LAT data at TeV energies. This is possible thanks to the new Fermi LAT IGRB measurement, which now extends up to an energy of 820 GeV. We quantify uncertainties in detail and show the potential this type of search offers for testing the WIMP paradigm with a complementary and truly cosmological probe of DM particle signals.

Fermi-LAT analyses show that the gamma-ray photon spectral indices Gamma_gamma of a large sample of blazars correlate with the vFv peak synchrotron frequency v_s according to the relation Gamma_gamma = d - k log v_s. The same function, with different constants d and k, also describes the relationship between Gamma_gamma and peak Compton frequency v_C. This behavior is derived analytically using an equipartition blazar model with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, k = k_EC = 3b/4 for external Compton processes and k = k_SSC = 9b/16 for synchrotron self-Compton (SSC) processes, where b is the log-parabola width parameter of the EED. The BL Lac object Mrk 501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from equipartition. Corrections are made to the spectral-index diagrams for a low-energy power-law EED and departures from equipartition, as constrained by absolute jet power. Analytic expressions are compared with numerical values derived from self-Compton and external Compton scattered gamma-ray spectra from Ly alpha broad-line region and IR target photons. The Gamma_gamma vs. v_s behavior in the model depends strongly on b, with progressively and predictably weaker dependences on gamma-ray detection range, variability time, and isotropic gamma-ray luminosity. Implications for blazar unification and blazars as ultra-high energy cosmic-ray sources are presented. Arguments by Ghisellini et al. (2014) that the jet power exceeds the accretion luminosity depend on the doubtful assumption that we are viewing at the Doppler angle.

This article reviews a few topics relevant to Galactic cosmic-ray astrophysics, focusing on the recent AMS-02 data release and Fermi Large Area Telescope data on the diffuse Galactic gamma-ray emissivity. Calculations are made of the diffuse cosmic-ray induced p+p --> pi^0 --> 2 gamma spectra, normalized to the AMS-02 cosmic-ray proton spectrum at ~ 10 - 100 GV, with and without a hardening in the cosmic-ray proton spectrum at rigidities R >~ 300 GV. A single power-law momentum "shock" spectrum for the local interstellar medium cosmic-ray proton spectrum cannot be ruled out from the gamma-ray emissivity data alone without considering the additional contribution of electron bremsstrahlung. Metallicity corrections are discussed, and a maximal range of nuclear enhancement factors from 1.52 to 1.92 is estimated.Origins of the 300 GV cosmic-ray proton and alpha-particle hardening are discussed.

In this paper, we use previously developed exact solutions to present some of the curious features of a force-free magnetosphere in a Kerr background. More precisely, we obtain a hitherto unseen timelike current in the force-free magnetosphere that does not flow along a geodesic. The electromagnetic field in this case happens to be magnetically dominated. This too is a feature that has entered the literature for the first time. Changing the sign of a single parameter in our solutions generates a spacelike current that creates an electromagnetic field that is electrically dominated.

The origin of the extragalactic $\gamma$-ray background (EGB) has been debated for some time. { The EGB comprises the $\gamma$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

The IceCube collaboration has reported neutrinos with energies between ~30 TeV and a few PeV that are significantly enhanced over the cosmic-ray induced atmospheric background. Viable high-energy neutrino sources must contain very high-energy and ultra-high energy cosmic rays while efficiently making PeV neutrinos. Gamma-ray Bursts (GRBs) and blazars have been considered as candidate cosmic-ray accelerators. GRBs, including low-luminosity GRBs, can be efficient PeV neutrino emitters for low bulk Lorentz factor outflows, although the photopion production efficiency needs to be tuned to simultaneously explain ultra-high-energy cosmic rays. Photopion production efficiency of cosmic-rays accelerated in the inner jets of flat spectrum radio quasars (FSRQs) is ~1-10% due to interactions with photons of the broad-line region (BLR), whereas BL Lac objects are not effective PeV neutrino sources due to the lack of external radiation fields. Photopion threshold effects with BLR photons suppress neutrino production below ~1 PeV, which imply that neutrinos from other sources would dominate over the diffuse neutrino intensity at sub-PeV energies. Reduction of the >> PeV neutrino flux can be expected when curving cosmic-ray proton distributions are employed. Considering a log-parabolic function to describe the cosmic-ray distribution, we discuss possible implications for particle acceleration in black-hole jets. Our results encourage a search for IceCube PeV neutrino events associated with gamma-ray loud FSRQs using Fermi-LAT data. In our model, as found in our previous work, the neutrino flux is suppressed below 1 PeV, which can be tested with increased IceCube exposure.

The standard blazar paradigm consists of a supermassive black hole that expels relativistic jets of magnetized plasma in our direction. This plasma entrains nonthermal synchrotron-emitting electrons that furthermore scatter internally-produced synchrotron photons as well as externally-produced photons from the accretion-disk, the broad-line region, and the infrared-emitting torus. This picture has been very successful in reproducing the two-humped blazar spectral energy distribution. Yet various discontents persist, including (1) ultra-short variability at TeV energies that is much shorter than the black hole's dynamical timescale; (2) very-high energy (VHE; >100 GeV) radiation from FSRQs; (3) evidence for a hard spectral component in high-synchrotron peaked objects, found when deabsorbing the measured VHE spectrum using conventional extragalactic background light (EBL) models; (4) an unusual slowly varying class of BL Lac objects; (5) location of the gamma-ray emission region. Some of these problems can be resolved by introducing a hadronic component into the blazar paradigm, consistent with the hypothesis that blazars are sources of ultra-high energy cosmic rays. Gamma-ray observations with the Cherenkov Telescope Array and neutrino observations with IceCube will be important to test this hypothesis.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We study high-energy neutrino production in inner jets of radio-loud active galactic nuclei (AGN), taking into account effects of external photon fields and the blazar sequence. We show that the resulting diffuse neutrino intensity is dominated by quasar-hosted blazars, in particular, flat spectrum radio quasars, and that PeV-EeV neutrino production due to photohadronic interactions with broadline and dust radiation is unavoidable if the AGN inner jets are ultrahigh-energy cosmic-ray (UHECR) sources. Their neutrino spectrum has a cutoff feature around PeV energies since target photons are due to Ly$\alpha$ emission. Because of infrared photons provided by the dust torus, neutrino spectra above PeV energies are too hard to be consistent with the IceCube data unless the proton spectral index is steeper than 2.5, or the maximum proton energy is $\lesssim100$ PeV. Thus, the simple model has difficulty in explaining the IceCube data. For the cumulative neutrino intensity from blazars to exceed $\sim{10}^{-8}~{\rm GeV}~{\rm cm}^{-2}~{\rm s}^{-1}~{\rm sr}^{-1}$, their local cosmic-ray energy generation rate would be $\sim10-100$ times larger than the local UHECR emissivity, but is comparable to the averaged gamma-ray blazar emissivity. Interestingly, future detectors such as the Askaryan Radio Array can detect $\sim0.1-1$ EeV neutrinos even in more conservative cases, allowing us to indirectly test the hypothesis that UHECRs are produced in the inner jets. We find that the diffuse neutrino intensity from radio-loud AGN is dominated by blazars with gamma-ray luminosity of $\gtrsim10^{48}~{\rm erg}~{\rm s}^{-1}$, and the arrival directions of their $\sim1-100$ PeV neutrinos correlate with the luminous blazars detected by Fermi.

We report on our study of high-energy properties of two peculiar TeV emitters: the "extreme blazar" 1ES 0347-121 and the "extreme blazar candidate" HESS J1943+213 located near the Galactic Plane. Both objects are characterized by quiescent synchrotron emission with flat spectra extending up to the hard X-ray range, and both were reported to be missing GeV counterparts in the Fermi-LAT 2-year Source Catalog. We analyze a 4.5 year accumulation of the Fermi-LAT data, resulting in the detection of 1ES 0347-121 in the GeV band, as well as in improved upper limits for HESS J1943+213. We also present the analysis results of newly acquired Suzaku data for HESS J1943+213. The X-ray spectrum is well represented by a single power law extending up to 25 keV with photon index 2.00+/-0.02 and a moderate absorption in excess of the Galactic value, in agreement with previous X-ray observations. No short-term X-ray variability was found over the 80 ks duration of the Suzaku exposure. Under the blazar hypothesis, we modeled the spectral energy distributions of 1ES 0347-121 and HESS J1943+213, and derived constraints on the intergalactic magnetic field strength and source energetics. We conclude that although the classification of HESS J1943+213 has not yet been determined, the blazar hypothesis remains the most plausible option, since in particular the broad-band spectra of the two analyzed sources along with the source model parameters closely resemble each other, and the newly available WISE and UKIDSS data for HESS J1943+213 are consistent with the presence of an elliptical host at the distance of approximately ~600 Mpc.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

We present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hours of VERITAS observations from three years, a multiwavelength light curve and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured differential very high energy (VHE; $E\ge100$ GeV) spectral indices are $\Gamma=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than $\tau=2$, where it is postulated that any variability would be small and occur on longer than year timescales if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

Blazar spectral models generally have numerous unconstrained parameters, leading to ambiguous values for physical properties like Doppler factor delta or fluid magnetic field B'. To help remedy this problem, a few modifications of the standard leptonic blazar jet scenario are considered. First, a log-parabola function for the electron distribution is used. Second, analytic expressions relating energy loss and kinematics to blazar luminosity and variability, written in terms of equipartition parameters, imply delta, B', and the principal electron Lorentz factor gamma'_pk. The external radiation field in a blazar is approximated by Ly alpha radiation from the broad line region (BLR) and ~0.1 eV infrared radiation from a dusty torus. When used to model 3C 279 SEDs from 2008 and 2009 reported by Hayashida et al. (2012), we derive delta ~ 20-30, B' ~ few G, and total (IR + BLR) external radiation field energy densities u ~ 0.01 - 0.001 erg/cm^3, implying an origin of the gamma-ray emission site in 3C 279 at the outer edges of the BLR. This is consistent with the gamma-ray emission site being located at a distance R <~ Gamma^2 c t_{var} ~ 0.1 (\Gamma/30)^2 (t_{var}/10^4 s) pc from the black hole powering 3C 279's jets, where t_{var} is the variability time scale of the radiation in the source frame, and at farther distances for narrow-jet and magnetic_reconnection models. Excess >~ 5 GeV gamma-ray emission observed with Fermi LAT from 3C 279 challenge the model, opening the possibility of hadronic origins of the emission. For low hadronic content, absolute jet powers of ~10% of the Eddington luminosity are calculated.

Time-resolved spectroscopy is performed on eight bright, long gamma-ray bursts (GRBs) dominated by single emission pulses that were observed with the {\it Fermi Gamma-ray Space Telescope}. Fitting the prompt radiation of GRBs by empirical spectral forms such as the Band function leads to ambiguous conclusions about the physical model for the prompt radiation. Moreover, the Band function is often inadequate to fit the data. The GRB spectrum is therefore modeled with two emission components consisting of optically thin nonthermal synchrotron radiation from relativistic electrons and, when significant, thermal emission from a jet photosphere, which is represented by a blackbody spectrum. To produce an acceptable fit, the addition of a blackbody component is required in 5 out of the 8 cases. We also find that the low-energy spectral index \alpha is consistent with a synchrotron component with \alpha = -0.81\pm 0.1. This value lies between the limiting values of \alpha = -2/3 and \alpha = -3/2 for electrons in the slow and fast-cooling regimes, respectively, suggesting ongoing acceleration at the emission site. The blackbody component can be more significant when using a physical synchrotron model instead of the Band function, illustrating that the Band function does not serve as a good proxy for a nonthermal synchrotron emission component. The temperature and characteristic emission-region size of the blackbody component are found to, respectively, decrease and increase as power laws with time during the prompt phase. In addition, we find that the blackbody and nonthermal components have separate temporal behaviors.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

GRB 130427A occurred in a relatively nearby galaxy; its prompt emission had the largest GRB fluence ever recorded. The afterglow of GRB 130427A was bright enough for the Nuclear Spectroscopic Telescope ARray (NuSTAR) to observe it in the 3-79 keV energy range long after its prompt emission (~1.5 and 5 days). This range, where afterglow observations were previously not possible, bridges an important spectral gap. Combined with Swift, Fermi and ground-based optical data, NuSTAR observations unambiguously establish a single afterglow spectral component from optical to multi-GeV energies a day after the event, which is almost certainly synchrotron radiation. Such an origin of the late-time Fermi/LAT > 10 GeV photons requires revisions in our understanding of collisionless relativistic shock physics.

We report the Fermi Large Area Telescope (LAT) detection of two very-high-energy (VHE, E>100 GeV) gamma-ray photons from the directional vicinity of the distant (redshift, z = 1.1) blazar PKS 0426-380. The null hypothesis that both the 134 and 122 GeV photons originate from unrelated sources can be rejected at the 5.5 sigma confidence level. We therefore claim that at least one of the two VHE photons is securely associated with the blazar, making PKS 0426-380 the most distant VHE emitter known to date. The results are in agreement with the most recent Fermi-LAT constraints on the Extragalactic Background Light (EBL) intensity, which imply a $z \simeq 1$ horizon for $\simeq$ 100 GeV photons. The LAT detection of the two VHE gamma-rays coincided roughly with flaring states of the source, although we did not find an exact correspondence between the VHE photon arrival times and the flux maxima at lower gamma-ray energies. Modeling the gamma-ray continuum of PKS 0426-380 with daily bins revealed a significant spectral hardening around the time of detection of the first VHE event (LAT photon index \Gamma\ $\simeq$ 1.4) but on the other hand no pronounced spectral changes near the detection time of the second one. This combination implies a rather complex variability pattern of the source in gamma rays during the flaring epochs. An additional flat component is possibly present above several tens of GeV in the EBL-corrected Fermi-LAT spectrum accumulated over the ~8-month high state.

Curvature effects in Gamma-ray bursts (GRBs) have long been a source of considerable interest. In a collimated relativistic GRB jet, photons that are off-axis relative to the observer arrive at later times than on-axis photons and are also expected to be spectrally softer. In this work, we invoke a relatively simple kinematic two-shell collision model for a uniform jet profile and compare its predictions to GRB prompt-emission data for observations that have been attributed to curvature effects such as the peak-flux--peak-frequency relation, i.e., the relation between the $\nu$F$_\nu$ flux and the spectral peak, E$_{pk}$ in the decay phase of a GRB pulse, and spectral lags. In addition, we explore the behavior of pulse widths with energy. We present the case of the single-pulse Fermi GRB 110920, as a test for the predictions of the model against observations.

More than 90% of the Galactic gas-related gamma-ray emissivity above 1 GeV is attributed to the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. Given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~4 GeV obtained with the Large Area Telescope on board the Fermi space observatory, it is timely to re-investigate this topic in detail, including the hadronic production mechanisms. The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

Recent data from the \emph{Fermi} Large Area Telescope have revealed about a dozen distant hard-spectrum blazars that have very-high-energy (VHE; $\gtrsim 100$ GeV) photons associated with them, but most of them have not yet been detected by imaging atmospheric Cherenkov telescopes. Most of these high-energy gamma-ray spectra, like those of other extreme high-frequency peaked BL Lac objects, can be well explained either by gamma rays emitted at the source or by cascades induced by ultra-high-energy cosmic rays, as we show specifically for KUV 00311$-$1938. We consider the prospects for detection of the VHE sources by the planned Cherenkov Telescope Array (CTA) and show how it can distinguish the two scenarios by measuring the integrated flux above $\sim 500$ GeV (depending on source redshift) for several luminous sources with $z \lesssim 1$ in the sample. Strong evidence for the origin of ultra-high-energy cosmic rays could be obtained from VHE observations with CTA. Depending on redshift, if the often quoted redshift of KUV 00311-1938 ($z = 0.61$) is believed, preliminary H.E.S.S. data favor cascades induced by ultra-high-energy cosmic rays. Accurate redshift measurements of hard-spectrum blazars are essential for this study.

Observations performed with the Fermi-LAT telescope have revealed the presence of a spectral break in the GeV spectrum of flat-spectrum radio quasars (FSRQs) and other low- and intermediate-synchrotron peaked blazars. We propose that this feature can be explained by Compton scattering of broad-line region (BLR) photons by a non-thermal population of electrons described by a log-parabolic function. We consider in particular a scenario in which the energy densities of particles, magnetic field, and soft photons in the emitting region are close to equipartition. We show that this model can satisfactorily account for the overall spectral energy distribution of the FSRQ 3C 454.3, reproducing the GeV spectal cutoff due to Klein-Nishina effects and a curving electron distribution.

Context: Cosmic rays are thought to be accelerated at supernova remnant (SNR) shocks, but conclusive evidence is lacking. Aims: New data from ground-based gamma-ray telescopes and the Large Area Telescope on the Fermi Gamma-ray Space Telescope are used to test this hypothesis. A simple model for gamma-ray production efficiency is compared with measured gamma-ray luminosities of SNRs, and the GeV to TeV fluxes ratios of SNRs are examined for correlations with SNR ages. Methods: The supernova explosion is modeled as an expanding spherical shell of material that sweeps up matter from the surrounding interstellar medium (ISM). The accumulated kinetic energy of the shell, which provides the energy available for nonthermal particle acceleration, changes when matter is swept up from the ISM and the SNR shell decelerates. A fraction of this energy is assumed to be converted into the energy of cosmic-ray electrons or protons. Three different particle radiation processes---nuclear pion-production interactions, nonthermal electron bremsstrahlung, and Compton scattering---are considered. Results: The efficiencies for gamma-ray production by these three processes are compared with gamma-ray luminosities of SNRs. Our results suggest that SNRs become less gamma-ray luminous at >~ 10^4 yr, and are consistent with the hypothesis that supernova remnants accelerate cosmic rays with an efficiency of ~10% for the dissipation of kinetic energy into nonthermal cosmic rays. Weak evidence for an increasing GeV to TeV flux ratio with SNR age is found.

Secondary nuclear production physics is receiving increased attention given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~40 GeV, obtained with the Large Area Telescope on board the Fermi space observatory. More than 90% of the gas-related emissivity above 1 GeV is attributed to gamma-rays from the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. The elementary kinematics of neutral pion production and decay are re-examined in light of two physics questions: does isobaric production follow a scaling behavior? and what is the minimum proton kinetic energy needed to make a gamma-ray of a certain energy formed through intermediate pi0 production? The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

A shock-accelerated particle flux \propto p^-s, where p is the particle momentum, follows from simple theoretical considerations of cosmic-ray acceleration at nonrelativistic shocks followed by rigidity-dependent escape into the Galactic halo. A flux of shock-accelerated cosmic-ray protons with s ~ 2.8 provides an adequate fit to the Fermi-LAT gamma-ray emission spectra of high-latitude and molecular cloud gas when uncertainties in nuclear production models are considered. A break in the spectrum of cosmic-ray protons claimed by Neronov, Semikoz, & Taylor (PRL, 108, 051105, 2012) when fitting the gamma-ray spectra of high-latitude molecular clouds is a consequence of using a cosmic-ray proton flux described by a power law in kinetic energy.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

The flat spectrum radio quasar (FSRQ) 4C +21.35 (PKS 1222+216) displays prominent nuclear infrared emission from ~1200 K dust. A 70 -- 400 GeV flare with ~10 min variations during half an hour of observations was found by the MAGIC telescopes, and GeV variability was observed on sub-day timescales with the Large Area Telescope on Fermi. We examine 4C +21.35, assuming that it is a source of ultra-high energy cosmic rays (UHECRs). UHECR proton acceleration in the inner jet powers a neutral beam of neutrinos, neutrons and gamma rays from photopion production. The radiative efficiency and production spectra of neutrals formed through photohadronic processes with isotropic external target photons of the broad line region and torus are calculated. Secondary radiations made by this process have a beaming factor ~\delta^5, where \delta is the Doppler factor. The pair-production optical depth for gamma rays and the photopion efficiency for UHECR neutrons as they pass through external isotropic radiation fields are calculated. If target photons come from the broad line region and dust torus, large Doppler factors, \delta >~100 are required to produce rapidly variable secondary radiation with isotropic luminosity >~1e47 erg/s at the pc scale. The \gamma-ray spectra from leptonic secondaries are calculated from cascades initiated by the UHECR neutron beam at the pc-scale region and fit to the flaring spectrum of 4C +21.35. Detection of >~100 TeV neutrinos from 4C +21.35 or other VHE blazars with IceCube or KM3NeT would confirm this scenario.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

The Second Catalog of Blazars and other Active Galactic Nuclei detected by the Fermi/LAT (2LAC) includes about 1100 sources, 886 of which comprise the Clean Sample. The general properties of the different populations of sources classified according to the strength of their emission lines (FSRQs, BL Lacs) or the estimated position of the synchrotron peak are reviewed.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

A simple formalism to describe nonthermal electron acceleration, evolution, and radiation in supernova remnants (SNRs) is presented. The electron continuity equation is analytically solved assuming that the nonthermal electron injection power is proportional to the rate at which the kinetic energy of matter swept up in an adiabatically expanding SNR shell. We apply this model to \fermi\ and HESS data from the SNR \rxj, and find that a one-zone leptonic model with Compton-scattered cosmic microwave background (CMB) and interstellar infrared photons has difficulty providing a good fit to its spectral energy distribution, provided the source is at a distance $\sim 1\ \kpc$ from the Earth. However, the inclusion of multiple zones, as hinted at by recent {\em Chandra} observations, does provide a good fit, but requires a second zone of compact knots with magnetic fields $B\sim 16\ \mu$G, comparable to shock-compressed fields found in the bulk of the remnant.

This is a pedagological review of some astrophysical highlights of the Fermi Gamma ray Observatory, including theoretical studies related mainly to extragalactic Fermi science.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

The spectra of BL Lac objects and Fanaroff-Riley I radio galaxies are commonly explained by the one-zone leptonic synchrotron self-Compton (SSC) model. Spectral modeling of correlated multiwavelength data gives the comoving magnetic field strength, the bulk outflow Lorentz factor and the emission region size. Assuming the validity of the SSC model, the Hillas condition shows that only in rare cases can such sources accelerate protons to much above 10^19 eV, so > 10^20 eV ultra-high-energy cosmic rays (UHECRs) are likely to be heavy ions if powered by this type of radio-loud active galactic nuclei (AGN). Survival of nuclei is shown to be possible in TeV BL Lacs and misaligned counterparts with weak photohadronic emissions. Another signature of hadronic production is intergalactic UHECR-induced cascade emission, which is an alternative explanation of the TeV spectra of some extreme non-variable blazars such as 1ES 0229+200 or 1ES 1101-232. We study this kind of cascade signal, taking into account effects of the structured extragalactic magnetic fields in which the sources should be embedded. We demonstrate the importance of cosmic-ray deflections on the gamma-ray flux, and show that required absolute cosmic-ray luminosities are larger than the average UHECR luminosity inferred from UHECR observations and can even be comparable to the Eddington luminosity of supermassive black holes. Future TeV gamma-ray observations using the Cherenkov Telescope Array and the High Altitude Water Cherenkov detector array can test for UHECR acceleration by observing >25 TeV photons from relatively low-redshift sources such as 1ES 0229+200, and > TeV photons from more distant radio-loud AGN.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

Lower limits on the power emitted in ultrahigh-energy cosmic rays (UHECRs), which are assumed to be protons with energy 1e17-1e20 eV, are derived for TeV blazars with the assumption that the observed TeV gamma rays are generated due to interactions of these protons with cosmic microwave photons. The limits depend on the spectrum of the injected UHECR protons. While for a -2.2 injection spectrum, the lower limits on the powers emitted in UHECRs by 1ES 0229+200, 1ES 1101-232 and 1ES 0347-121 are lower than their respective synchrotron luminosities (1e46 erg/s); in the case of 1ES 1426+428 it exceeds the corresponding synchrotron luminosity by up to an order of magnitude. The proposed Auger North Observatory should be able to detect 4e19 eV cosmic ray protons from 1ES 1426+428 within a few years of operation and test the TeV gamma-ray production model by UHECR energy losses while propagating along the line-of-sight, or constrain the intergalactic magnetic field to be larger than 1e-16 G in case of no detection. The lower limits on the apparent-isotropic jet power from accelerated 1e10-1e20 eV proton spectra in the blazar jet is of the order of the Eddington luminosity of a 1e9 solar mass black hole for a cosmic-ray injection spectrum -2.2 or harder for all blazars considered except for 1ES 1426+428. In the case of the latter the apparent-isotropic jet power exceeds the Eddington luminosity by an order of magnitude. For an injection spectrum softer than -2.2, as is required to fit the observed cosmic-ray data above 1e17-1e18 eV, the Eddington luminosity is exceeded by the lower limits on the jet power for all blazars considered.

Recent Fermi and VERITAS observations of the prototypical Type Ia supernova remnant (SNR) Tycho have discovered gamma-rays with energies E in the range from ~0.4 GeV to 10 TeV. Crucial for the theory of Galactic cosmic-ray origin is whether the gamma-rays from SNRs are produced by accelerated hadrons (protons and ions), or by relativistic electrons. Here we show that the broadband radiation spectrum of Tycho can be explained within the framework of a two-zone leptonic model, which is likely to apply to every SNR. A model with hadrons can also fit the radiation spectrum. The hadronic origin of gamma-rays can be confirmed by Fermi spectral measurements of Tycho and other SNRs at energies below ~300 MeV.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

Fermi has provided the largest sample of {\gamma}-ray selected blazars to date. In this work we use a complete sample of FSRQs detected during the first year of operation to determine the luminosity function (LF) and its evolution with cosmic time. The number density of FSRQs grows dramatically up to redshift \sim0.5-2.0 and declines thereafter. The redshift of the peak in the density is luminosity dependent, with more luminous sources peaking at earlier times; thus the LF of {\gamma}-ray FSRQs follows a luminosity-dependent density evolution similarly to that of radio-quiet AGN. Also using data from the Swift Burst Alert Telescope we derive the average spectral energy distribution of FSRQs in the 10 keV-100 GeV band and show that there is no correlation of the peak {\gamma}-ray luminosity with {\gamma}-ray peak frequency. The coupling of the SED and LF allows us to predict that the contribution of FSRQs to the Fermi isotropic {\gamma}-ray background is 9.3(+1.6/-1.0) (\pm3% systematic uncertainty) in the 0.1-100GeV band. Finally we determine the LF of unbeamed FSRQs, finding that FSRQs have an average Lorentz factor of {\gamma} = 11.7(+3.3/-2.2), that most are seen within 5\circ of the jet axis, and that they represent only ~0.1 % of the parent population.

Some questions raised by Fermi-LAT data about blazars are summarized, along with attempts at solutions within the context of leptonic models. These include both spectral and statistical questions, including the origin of the GeV breaks in low-synchrotron peaked blazars, the location of the gamma-ray emission sites, the correlations in the spectral energy distributions with luminosity, and the difficulty of synchrotron/SSC models to fit the spectra of some TeV blazars.

Spectral lag, which is defined as the difference in time of arrival of high and low energy photons, is a common feature in Gamma-ray Bursts (GRBs). Previous investigations have shown a correlation between this lag and the isotropic peak luminosity for long duration bursts. However, most of the previous investigations used lags extracted in the observer-frame only. In this work (based on a sample of 43 Swift long GRBs with known redshifts), we present an analysis of the lag-luminosity relation in the GRB source-frame. Our analysis indicates a higher degree of correlation -0.82 +/- 0.05 (chance probability of ~ 5.5 x 10^-5) between the spectral lag and the isotropic peak luminosity, Liso, with a best-fit power-law index of -1.2 +/- 0.2, such that Liso proportional to lag^-1.2. In addition, there is an anti-correlation between the source-frame spectral lag and the source-frame peak energy of the burst spectrum, E_pk(1+z).

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

NGC 6251 is a luminous radio galaxy ~104 Mpc away that was detected significantly with the Fermi Gamma-ray Space Telescope, and before that with EGRET (onboard the Compton Gamma-ray Observatory). Different observational constraints favor a nuclear origin for the gamma-ray emission. Here we present a study of the spectral energy distribution (SED) of the core of NGC 6251, and give results of modeling in the one-zone synchrotron/SSC framework. The SSC model provides a good description of the radio to gamma-ray emission but, as for other misaligned sources, predicts a lower Lorentz factor (\Gamma~2.4) than typically found when modeling blazars. If the blazar unification scenario is correct, this seems to point to the presence of at least two emitting regions in these objects, one with a higher and one with a lower Lorentz factor. The solution of a structured jet, with a fast moving spine surrounded by a slow layer, is explored and the consequences of the two models for the jet energetics and evolution are discussed.

In this paper, we provide exact solutions for the extraction of energy from a rotating black hole via both the electromagnetic Poynting flux and matter currents. By appropriate choice of a radially independent poloidal function $\Lambda(\theta)$, we find solutions where the dominant outward energy flux is along the polar axis, consistent with a jet-like collimated outflow, but also with a weaker flux of energy along the equatorial plane. Unlike all the previously obtained solutions (Blandford & Znajek (1977), Menon & Dermer (2005)), the magnetosphere is free of magnetic monopoles everywhere.

Recent claims that the strength B_IGMF of the intergalactic magnetic field (IGMF) is >~ 1e-15 G are based on upper limits to the expected cascade flux in the GeV band produced by blazar TeV photons absorbed by the extragalactic background light. This limit depends on an assumption that the mean blazar TeV flux remains constant on timescales >~2 (B_ IGMF/1e-18 G)^2 / (E/{10 GeV})^2 yr for an IGMF coherence length ~ 1 Mpc, where E is the measured photon energy. Restricting TeV activity of 1ES 0229+200 to ~3 -- 4 years during which the source has been observed leads to a more robust lower limit of B_IGMF >~ 1e-18 G, which can be larger by an order of magnitude if the intrinsic source flux above ~5 -- 10 TeV from 1ES 0229+200 is strong.

A checklist of criteria is presented to help establish the sources of ultra-high energy cosmic rays (UHECRs). These criteria are applied to long-duration GRBs in order to determine if they are UHECR sources. The evidence seems to favor blazars and radio galaxies (or other sources) over GRBs.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

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 gamma-ray observations of radio galaxies on blazar unification scenarios, variability in colliding shells, and whether blazars are sources of ultra-high energy cosmic rays.

Fermi Gamma ray Space Telescope measurements of spectra, variability time scale, and maximum photon energy give lower limits to the apparent jet powers and, through gammagamma opacity arguments, the bulk Lorentz factors of relativistic jets. The maximum cosmic-ray particle energy is limited by these two quantities in Fermi acceleration scenarios. Recent data are used to constrain the maximum energies of cosmic-ray protons and Fe nuclei accelerated in colliding shells of GRBs and blazars. The Fermi results indicate that Fe rather than protons are more likely to be accelerated to ultra-high energies in AGNs, whereas powerful GRBs can accelerate both protons and Fe to >~ 10^{20} eV. Emissivity of nonthermal radiation from radio galaxies and blazars is estimated from the First Fermi AGN Catalog, and shown to favor BL Lac objects and FR1 radio galaxies over flat spectrum radio quasars, FR2 radio galaxies, and long-duration GRBs as the sources of UHECRs.

Fermi gamma-ray telescope data of GRB 080916C with ~1e55 erg in apparent isotropic gamma-ray energy, show a several second delay between the rise of 100 MeV - GeV radiation compared with keV - MeV radiation. Here we show that synchrotron radiation from cosmic ray protons accelerated in GRBs, delayed by the proton synchrotron cooling timescale in a jet of magnetically-dominated shocked plasma moving at highly relativistic speeds with bulk Lorentz factor Gamma ~ 500, could explain this result. A second generation electron synchrotron component from attenuated proton synchrotron radiation makes enhanced soft X-ray to MeV gamma-ray emission. Long GRBs with narrow, energetic jets accelerating particles to ultra-high energies could explain the Auger observations of UHE cosmic rays from sources within 100 Mpc for nano-Gauss intergalactic magnetic fields. The total energy requirements in a proton synchrotron model are proportional to Gamma^(16/3). This model for GRB 080916C is only plausible if Gamma ~< 500 and the jet opening angle is ~ 1 degree.

Fermi LAT (Large Area Telescope) and GBM (Gamma ray Burst Monitor) observations of GRBs are briefly reviewed, keeping in mind EGRET expectations. Using gamma\gamma constraints on outflow Lorentz factors, leptonic models are pitted against hadronic models, and found to be energetically favored. Interpretation of the Fermi data on GRBs helps establish whether GRBs accelerate cosmic rays, including those reaching $\approx 10^{20}$ eV.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

Fermi Gamma ray Space Telescope observations of the flat spectrum radio quasar 3C~454.3 show a spectral-index change $\Delta \Gamma \cong 1.2\pm 0.3$ at break energy $E_{br} \approx 2.4\pm0.3$ GeV. Such a sharp break is inconsistent with a cooling electron distribution and is poorly fit with a synchrotron self-Compton model. We show that a combination of two components, namely the Compton-scattered disk and broad-line region (BLR) radiations, explains this spectral break and gives a good fit to the quasi-simultaneous radio, optical/UV, X-ray, and $\gamma$-ray spectral energy distribution observed in 2008 August. A sharp break can be produced independent of the emitting region's distance from the central black hole if the BLR has a gradient in density $\propto R^{-2}$, consistent with a wind model for the BLR.

If gamma-ray bursts are sources of ultra-high energy cosmic rays, then radiative signatures of hadronic acceleration are expected in GRB data. Observations with the Fermi Gamma-ray Space Telescope offer the best means to search for evidence of UHECRs in GRBs through electromagnetic channels. Various issues related to UHECR acceleration in GRBs are reviewed, with a focus on the question of energetics.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

The extragalactic background light (EBL) from the far infrared through the visible and extending into the ultraviolet is thought to be dominated by starlight, either through direct emission or through absorption and reradiation by dust. This is the most important energy range for absorbing $\g$-rays from distant sources such as blazars and gamma-ray bursts and producing electron positron pairs. In previous work we presented EBL models in the optical through ultraviolet by consistently taking into account the star formation rate (SFR), initial mass function (IMF) and dust extinction, and treating stars on the main sequence as blackbodies. This technique is extended to include post-main sequence stars and reprocessing of starlight by dust. In our simple model, the total energy absorbed by dust is assumed to be re-emitted as three blackbodies in the infrared, one at 40 K representing warm, large dust grains, one at 70 K representing hot, small dust grains, and one at 450 K representing polycyclic aromatic hydrocarbons. We find our best fit model combining the Hopkins and Beacom SFR using the Cole et al. parameterization with the Baldry and Glazebrook IMF agrees with available luminosity density data at a variety of redshifts. Our resulting EBL energy density is quite close to the lower limits from galaxy counts and in some cases below the lower limits, and agrees fairly well with other recent EBL models shortward of about 5 $\mu$m. Deabsorbing TeV $\g$-ray spectra of various blazars with our EBL model gives results consistent with simple shock acceleration theory. We also find that the universe should be optically thin to $\g$-rays with energies less than 20 GeV.

Bright gamma-ray flares observed from sources far beyond our Galaxy are best explained if enormous amounts of energy are liberated by black holes. The highest-energy particles in nature--the ultra-high energy cosmic rays--cannot be confined by the Milky Way's magnetic field, and must originate from sources outside our Galaxy. Here we summarize the themes of our book, "High Energy Radiation from Black Holes: Gamma Rays, Cosmic Rays, and Neutrinos", just published by Princeton University Press. In this book, we develop a mathematical framework that can be used to help establish the nature of gamma-ray sources, to evaluate evidence for cosmic-ray acceleration in blazars, GRBs and microquasars, to decide whether black holes accelerate the ultra-high energy cosmic rays, and to determine whether the Blandford-Znajek mechanism for energy extraction from rotating black holes can explain the differences between gamma-ray blazars and radio-quiet AGNs.

We investigate the spectral and timing signatures of the internal-shock model for blazars. For this purpose, we develop a semi-analytical model for the time-dependent radiative output from internal shocks arising from colliding relativistic shells in a blazar jet. The emission through synchrotron and synchrotron-self Compton (SSC) radiation as well as Comptonization of an isotropic external radiation field are taken into account. We evaluate the discrete correlation function (DCF) of the model light curves in order to evaluate features of photon-energy dependent time lags and the quality of the correlation, represented by the peak value of the DCF. The almost completely analytic nature of our approach allows us to study in detail the influence of various model parameters on the resulting spectral and timing features. This paper focuses on a range of parameters in which the gamma-ray production is dominated by Comptonization of external radiation, most likely appropriate for gamma-ray bright flat-spectrum radio quasars (FSRQs) or low-frequency peaked BL Lac objects (LBLs). In most cases relevant for FSRQs and LBLs, the variability of the optical emission is highly correlated with the X-ray and high-energy (HE: > 100 MeV) gamma-ray emission. Our baseline model predicts a lead of the optical variability with respect to the higher-energy bands by 1 - 2 hours and of the HE gamma-rays before the X-rays by about 1 hour. We show that variations of certain parameters may lead to changing signs of inter-band time lags, potentially explaining the lack of persistent trends of time lags in most blazars.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

Results of a phenomenological model to estimate the GRB detection rate by the Fermi Gamma ray Space Telescope are reported. This estimate is based on the BATSE 4B GRB fluence distribution, the mean ratio of fluences measured at 100 MeV - 5 GeV with EGRET and at 20 keV - 2 MeV with BATSE, and the mean EGRET GRB spectrum for the 5 EGRET spark-chamber GRBs. For a 10% fluence ratio and a number spectral index alpha_1 = -2 at 100 MeV - 5 GeV energies, we estimate a rate of ~ 20 and 4 GRBs per yr in the Fermi Large Area Telescope field of view with at least 5 photons with energy E > 100 MeV and E > 1 GeV, respectively. We also estimate ~ 1.5 GRBs per yr in the Fermi FoV where at least 1 photon with energy E > 10 GeV is detected. For these parameters, we estimate = 1 - 2 GRBs per year detected with the Fermi telescope with more than 100 gamma rays with E >~ 100$ MeV. Comparison predictions for alpha_1 = -2.2, different fluence ratios, and the AGILE gamma-ray satellite are made. Searches for different classes of GRBs using a diagram plotting 100 MeV - 10 GeV fluence vs. 20 keV - 20 MeV fluence is considered as a way to search for separate classes of GRBs and, specifically, spectral differences between the short-hard and long duration GRB classes, and for hard components in GRBs.

The Auger Collaboration reports that the arrival directions of >60 EeV ultra-high energy cosmic rays (UHECRs) cluster along the supergalactic plane and correlate with active galactic nuclei (AGN) within ~100 Mpc. The association of several events with the nearby radio galaxy Centaurus A supports the paradigm that UHECRs are powered by supermassive black-hole engines and accelerated to ultra-high energies in the shocks formed by variable plasma winds in the inner jets of radio galaxies. The GZK horizon length of 75 EeV UHECR protons is ~100 Mpc, so that the Auger results are consistent with an assumed proton composition of the UHECRs. In this scenario, the sources of UHECRs are FR II radio galaxies and FR I galaxies like Cen A with scattered radiation fields that enhance UHECR neutral-beam production. Radio galaxies with jets pointed away from us can still be observed as UHECR sources due to deflection of UHECRs by magnetic fields in the radio lobes of these galaxies. A broadband ~1 MeV -- 10 EeV radiation component in the spectra of blazar AGN is formed by UHECR-induced cascade radiation in the extragalactic background light (EBL). This emission is too faint to be seen from Cen A, but could be detected from more luminous blazars.

A fundamental question that can be answered in the next decade is: WHAT IS THE ORIGIN OF THE HIGHEST ENERGY COSMIC PARTICLES? The discovery of the sources of the highest energy cosmic rays will reveal the workings of the most energetic astrophysical environments in the recent universe. Candidate sources range from the birth of compact objects to explosions related to gamma-ray bursts or generated around supermassive black holes in active galactic nuclei. In addition to beginning a new era of high-energy astrophysics, the study of ultra-high energy cosmic rays will constrain the structure of the Galactic and extragalactic magnetic fields. The propagation of these particles from source to Earth also probes the cosmic background radiation and gives insight into particle interactions at orders of magnitude higher energy than can be achieved in terrestrial laboratories. Next generation observatories designed to study the highest energy cosmic rays will have unprecedented sensitivity to ultra-high energy photons and neutrinos, which will further illuminate the workings of the universe at the most extreme energies. For this challenge to be met during the 2010-2020 decade, a significant increase in the integrated exposure to cosmic rays above 6 1019 eV will be necessary. The technical capabilities for answering this open question are at hand and the time is ripe for exploring Charged Particle Astronomy.

Gamma-ray bursts (GRBs) have been an enigma since their discovery forty years ago. However, considerable progress unraveling their mysteries has been made in recent years. Developments in observations, theory, and instrumentation have prepared the way so that the next decade can be the one in which we finally answer the question, "What are gamma-ray bursts?" This question encompasses not only what the progenitors are that produce the GRBs, but also how the enormous luminosity of the GRBs, concentrated in gamma rays, is achieved. Observations across the electromagnetic spectrum, from both the ground and space, will be required to fully tackle this important question. This white paper, mostly distilled from a recent study commissioned by the Division of Astrophysics of the American Physical Society, focuses on what very high energy (~100 GeV and above) gamma-ray observations can contribute. Very high energy gamma rays probe the most extreme high energy particle populations in the burst environment, testing models of lepton and proton acceleration in GRBs and constraining the bulk Lorentz factor and opacity of the outflow. Sensitivity improvements of more than an order of magnitude in the very high energy gamma-ray band can be achieved early in the next decade, in order to contribute to this science.

TeV gamma rays from distant astrophysical sources are attenuated due to electron-positron pair creation by interacting with ultraviolet/optical to infrared photons which fill the universe and are collectively known as the extra-galactic background light (EBL). We model the ~0.1-10 eV starlight component of the EBL derived from expressions for the stellar initial mass function, star formation history of the universe, and wavelength-dependent absorption of a large sample of galaxies in the local universe. These models are simultaneously fitted to the EBL data as well as to the data on the stellar luminosity density in our local universe. We find that the models with modified Salpeter A initial mass function together with Cole et al. (2001) or Hopkins and Beacom (2006) star formation history best represent available data. Since no dust emission is included, our calculated EBL models can be interpreted as the lower limits in the ~0.1-1 eV range. We present simple analytic fits to the best-fit EBL model evolving with redshift. We then proceed to calculate gamma-ray opacities, and absorption of ~10-300 GeV gamma-rays coming from different redshifts. We discuss implications of our results for the Fermi Gamma Ray Space Telescope and ground-based Air Cherenkov Telescopes.

Gamma-ray bursts (GRBs) have long been considered as candidates of ultrahigh-energy cosmic rays (UHECRs). We investigate the signatures of CR proton acceleration in the GRBs by consistently taking into account their hadronic and electromagnetic interactions. We discuss the implications of our findings for high-energy gamma ray observations with the recently launched Fermi Gamma-ray Space Telescope.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

Protons accelerated in the internal shocks of a long duration gamma ray burst can escape the fireball as cosmic rays by converting to neutrons. Hadronic interactions of these neutrons inside a stellar wind bubble created by the progenitor star will produce TeV gamma rays via neutral meson decay and synchrotron radiation by charged pion-decay electrons in the wind magnetic field. Such gamma rays should be observable from nearby gamma ray bursts by currently running and upcoming ground-based detectors.

We extend a method for modeling synchrotron and synchrotron self-Compton radiations in blazar jets to include external Compton processes. The basic model assumption is that the blazar radio through soft X-ray flux is nonthermal synchrotron radiation emitted by isotropically-distributed electrons in the randomly directed magnetic field of outflowing relativistic blazar jet plasma. Thus the electron distribution is given by the synchrotron spectrum, depending only on the Doppler factor $\delta_{\rm D}$ and mean magnetic field $B$, given that the comoving emission region size scale $R_b^\prime \lesssim c \dD t_v/(1+z)$, where $t_v$ is variability time and $z$ is source redshift. Generalizing the approach of Georganopoulos, Kirk, and Mastichiadis (2001) to arbitrary anisotropic target radiation fields, we use the electron spectrum implied by the synchrotron component to derive accurate Compton-scattered $\gamma$-ray spectra throughout the Thomson and Klein-Nishina regimes for external Compton scattering processes. We derive and calculate accurate $\gamma$-ray spectra produced by relativistic electrons that Compton-scatter (i) a point source of radiation located radially behind the jet, (ii) photons from a thermal Shakura-Sunyaev accretion disk and (iii) target photons from the central source scattered by a spherically-symmetric shell of broad line region (BLR) gas. Calculations of broadband spectral energy distributions from the radio through $\gamma$-ray regimes are presented, which include self-consistent $\gamma\gamma$ absorption on the same radiation fields that provide target photons for Compton scattering. Application of this baseline flat spectrum radio/$\gamma$-ray quasar model is considered in view of data from $\gamma$-ray telescopes and contemporaneous multi-wavelength campaigns.

We report the results from our analysis of {\it Suzaku} XIS (0.5-10 keV) and HXD/PIN (15-40 keV) observations of five well-known local ULIRGs: {\em IRAS} F05189-2524, {\em IRAS} F08572+3915, Mrk 273, PKS 1345+12, and Arp 220. The XIS observations of F05189-2524 and Mrk 273 reveal strong iron lines consistent with Fe K$\alpha$ and changes in spectral shapes with respect to previous {\it Chandra} and {\it XMM-Newton} observations. Mrk 273 is also detected by the HXD/PIN at $\sim$1.8-$\sigma$. For F05189-2524, modeling of the data from the different epochs suggests that the change in spectral shape is likely due to the central source switching off, leaving behind a residual reflection spectrum, or an increase in the absorbing column. An increase in the covering fraction of the absorber can describe the spectral variations seen in Mrk 273, although a reduction in the intrinsic AGN luminosity cannot be formally ruled out. The {\it Suzaku} spectra of Mrk 273 are well fit by a ~94% covering fraction model with a column density of $\sim10^{24}$ cm$^{-2}$. The absorption-corrected log[$L_{\rm 2-10 keV}$ / $L_{\rm IR}$] ratio is consistent with those found in PG Quasars. The 0.5-10 keV spectrum of PKS 1345+12 and Arp 220 seem unchanged from previous observations and their hard X-ray emission is not convincingly detected by the HXD/PIN. The large column density derived from CO observations and the large equivalent width of an ionized Fe line in Arp 220 can be reconciled by an ionized reflection model. F08572+3915 is undetected in both the XIS and HXD/PIN, but the analysis of unpublished {\em Chandra} data provides a new measurement at low energies.

Kinematical and luminosity relations for black-hole jet sources are reviewed. If the TeV flares observed from PKS 2155-304 in 2006 July are assumed to originate from a black hole with mass $\approx 10^8 M_8 M_\odot$, then the $\sim 5$ minute variability timescale is consistent with the light-travel time across the Schwarzschild radius of the black hole if $M_8\sim 1$. The absolute jet power in a synchrotron/SSC model exceeds, however, the Eddington luminosity for a black hole with $M_8\sim 1$ unless the jet is highly efficient. The maximum Blandford-Znajek power is $\sim 10^{46}M_8$ ergs s$^{-1}$ if the magnetic-field energy density threading the horizon is equated with the luminous energy density in the vicinity of the black hole. An external Compton component can relax power requirements, so a black hole with mass $\sim 10^8 M_\odot$ could explain the observed flaring behavior. For the Swift and HESS data taken in 2006 July, relativistic outflows with bulk Lorentz factor $\Gamma \gtrsim 30$ satisfy $\gamma$-$\gamma$ attenuation limits. If this system harbors a binary black hole, then the accretion disk from a more massive, $\sim 10^9 M_\odot$ black-hole primary would make an additional external radiation component. Dual thermal accretion disk signatures would confirm this scenario.

This is a report on the findings of the gamma ray burst working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe gamma ray bursts at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that major advances are possible and that the detection of very high energy emission would have strong implications for GRB models, as well as cosmic ray origin.

Ultra-high energy cosmic rays (UHECRs) accelerated in the jets of active galactic nuclei can accumulate in high magnetic field, ~100 kpc-scale regions surrounding powerful radio galaxies. Photohadronic processes involving UHECRs and photons of the extragalactic background light make ultra-relativistic electrons and positrons that initiate electromagnetic cascades, leading to the production of a gamma-ray synchrotron halo. We calculate the halo emission in the case of Cygnus A and show that it should be detectable with the Fermi Gamma ray Space Telescope and possibly detectable with ground-based gamma-ray telescopes if radio galaxies are the sources of UHECRs.

Gamma-ray bursts (GRBs) are a mixed class of sources consisting of, at least, the long duration and short-hard subclasses, the X-ray flashes, and the low-luminosity GRBs. In all cases, the release of enormous amounts of energy on a short timescale makes an energetic, relativistic or mildly relativistic fireball that expands until it reaches a coasting Lorentz factor determined by the amount of baryons mixed into the fireball. Radiation is produced when the blast wave interacts with the surrounding medium at an external shock, or when shell collisions dissipate kinetic energy at internal shocks. This series of notes is organized as follows: (1) The observational situation of GRBs is summarized; (2) Progenitor models of GRBs are described; (3) An overview of the the blast-wave physics used to model leptonic emissions is given; (4) GRB physics is applied to hadronic acceleration and ultra-high energy cosmic ray production; (5) Prospects for GRB physics and gamma-ray astronomy with the Fermi Gamma-ray Space Telescope (FGST, formerly GLAST), and space-based and ground-based observatories are considered. Also included are exercises and problems.

We introduce a methodology for analysis of multiwavelength data from X-ray selected BL Lac (XBL) objects detected in the TeV regime. By assuming that the radio--through--X-ray flux from XBLs is nonthermal synchrotron radiation emitted by isotropically-distributed electrons in the randomly oriented magnetic field of a relativistic blazar jet, we obtain the electron spectrum. This spectrum is then used to deduce the synchrotron self-Compton (SSC) spectrum as a function of the Doppler factor, magnetic field, and variability timescale. The variability timescale is used to infer the comoving blob radius from light travel-time arguments, leaving only two parameters. With this approach, we accurately simulate the synchrotron and SSC spectrum of flaring XBLs in the Thomson through Klein-Nishina regimes. Photoabsorption by interactions with internal jet radiation and the intergalactic background light (IBL) is included. Doppler factors, magnetic fields, and absolute jet powers are obtained by fitting the {\em HESS} and {\em Swift} data of the recent giant TeV flare observed from \object{PKS 2155--304}. For the contemporaneous {\em Swift} and {\em HESS} data from 28 and 30 July 2006, respectively, Doppler factors $\gtrsim 60$ and absolute jet powers $\gtrsim 10^{46}$ ergs s$^{-1}$ are required for a synchrotron/SSC model to give a good fit to the data, for a low intensity of the IBL and a ratio of 10 times more energy in hadrons than nonthermal electrons. Fits are also made to a TeV flare observed in 2001 from Mkn 421 which require Doppler factors $\gtrsim 30$ and jet powers $\gtrsim 10^{45}$ erg s$^{-1}$.

An analysis of the interaction between a spherical relativistic blast-wave shell and a stationary cloud with a spherical cap geometry is performed assuming that the cloud width << x, where x is the distance of the cloud from the GRB explosion center. The interaction is divided into three phases: (1) a collision phase with both forward and reverse shocks; (2) a penetration phase when either the reverse shock has crossed the shell while the forward shock continues to cross the cloud, or vice versa; and (3) an expansion phase when, both shocks having crossed the cloud and shell, the shocked fluid expands. Temporally evolving spectral energy distributions (SEDs) are calculated for the problem of the interaction of a blast-wave shell with clouds that subtend large and small angles compared with the Doppler(-cone) angle 1/Gamma_0, where Gamma_0 is the coasting Lorentz factor. The Lorentz factor evolution of the shell/cloud collision is treated in the adiabatic limit. Behavior of the light curves and SEDs on, e.g., Gamma_0, shell-width parameter eta, and properties and locations of the cloud is examined. Short timescale variability (STV) in GRB light curves, including ~100 keV gamma-ray pulses observed with BATSE and delayed ~1 keV X-ray flares found with Swift, can be explained by emissions from an external shock formed by the GRB blast wave colliding with small density inhomogeneities in the "frozen pulse" approximation (eta -> 0, where Delta_0 + eta x/Gamma_0^2 is the blast-wave shell width), and perhaps in the thin-shell approximation (eta ~ 1/Gamma_0), but not when eta ~ 1. If this approximation is valid, then external shock processes could make the dominant prompt and afterglow emissions in GRB light curves, consistent with short delay two-step collapse models for GRBs.

Observations of very-high-energy (VHE, E > 250 GeV) gamma-ray emission from several blazars at z > 0.1 have placed stringent constraints on the elusive spectrum and intensity of the intergalactic infrared background radiation (IIBR). Correcting their observed VHE spectrum for gamma-gamma absorption even by the lowest plausible level of the IIBR provided evidence for a very hard (photon spectral index Gamma_{ph} < 2) intrinsic source spectrum out to TeV energies. Such a hard VHE gamma-ray spectrum poses a serious challenge to the conventional synchrotron-self-Compton interpretation of the VHE emission of TeV blazars and suggests the emergence of a separate emission component beyond a few hundred GeV. Here we propose that such a very hard, slowly variable VHE emission component in TeV blazars may be produced via Compton upscattering of Cosmic Microwave Background (CMB) photons by shock-accelerated electrons in an extended jet. For the case of 1ES 1101-232, this component could dominate the bolometric luminosity of the extended jet if the magnetic fields are of the order of typical intergalactic magnetic fields B ~ 10 micro-Gauss and electrons are still being accelerated out to TeV energies gamma > 4 X 10^6) on kiloparsec scales along the jet.

Auger results on clustering of > 60 EeV ultra-high energy cosmic ray (UHECR) ions and the interpretation of the gamma-ray spectra of TeV blazars are connected by effects from the extragalactic background light (EBL). The EBL acts as an obscuring medium for gamma rays and a reprocessing medium for UHECR ions and protons, causing the GZK cutoff. The study of the physics underlying the coincidence between the GZK energy and the clustering energy of UHECR ions favors a composition of > 60 EeV UHECRs in CNO group nucleons. This has interesting implications for the sources of UHECRs. We also comment on the Auger analysis.

This is a short report on the preliminary findings of the gamma ray burst (GRB) working group for the white paper on the status and future of very high energy (VHE; >50 GeV) gamma-ray astronomy. The white paper discusses the status of past and current attempts to observe GRBs at GeV-TeV energies, including a handful of low-significance, possible detections. The white paper concentrates on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that the detection of VHE emission would have strong implications for GRB models, as well as cosmic ray origin. In particular, the extended emission phase (including both afterglow emission and possible flaring) of nearby long GRBs could provide the best possibility for detection. The difficult-to-obtain observations during the prompt phase of nearby long GRBs and short GRBs could also provide particularly strong constraints on the opacity and bulk Lorentz factors surrounding the acceleration site. The synergy with upcoming and existing observatories will, of course, be critical for both identification of GRBs and for multiwavelength/multimessenger studies.

Results of a leptonic jet model for the prompt emission and early afterglows of GRBs are presented. The synchrotron component is modeled with the canonical Band spectrum and the synchrotron self-Compton component is calculated from the implied synchrotron-emitting electron spectrum in a relativistic plasma blob. In the comoving frame the magnetic field is assumed to be tangled and the electron and photon distributions are assumed to be isotropic. The Compton-scattered spectrum is calculated using the full Compton cross-section in the Thomson through Klein-Nishina using the Jones formula. Pair production photoabsorption, both from ambient radiation in the jet and from the extragalactic background light (EBL), is taken into account. Results are presented as a function of a small set of parameters: the Doppler factor, the observed variability timescale, the comoving magnetic field, the peak synchrotron flux, and the redshift of the burst. Model predictions will be tested by multiwavelength observations, including the {\em Swift} and {\em GLAST} satellites, which will provide unprecedented coverage of GRBs.

The discoveries of the GZK cutoff with the HiRes and Auger Observatories and the discovery by Auger of clustering of >~60 EeV ultra-high energy cosmic rays (UHECRs) towards nearby <~75 Mpc) AGNs along the supergalactic plane establishes the astrophysical origin of the UHECRs. The likely sources of the UHECRs are gamma-ray bursts and radio-loud AGNs because: (1) they are extragalactic; (2) they are sufficiently powerful; (3) acceleration to ultra-high energies can be achieved in their relativistic ejecta; (4) anomalous X-ray and gamma-ray features can be explained by nonthermal hadron acceleration in relativistic blast waves; and (5) sources reside within the GZK radius. Two arguments for acceleration to UHE are presented, and limits on UHECR ion acceleration are set. UHECR ions are shown to be able to survive without photodisintegrating while passing through the AGN scattered radiation field, even if launched deep in the broad line region. UHECR injection throughout cosmic time fits the measured energy spectrum of UHECRs, at least for protons. Local UHECR proton and ion interaction and energy-loss mean free paths are calculated using an empirical fit to the extragalactic background light (EBL) at IR and optical energies. Minimum intergalactic magnetic (IGM) fields ~1e-11 G are derived from clustering assuming specific source origins, e.g., Cen A, nearby AGNs, or GRBs for the super-GZK CRs seen with Auger. Besides distinct cosmic-ray induced gamma-ray signatures that should be observed with the Gamma ray Large Area Space Telescope (GLAST), source and GZK neutrino detections and the arrival distribution of UHECR in direction and time can finally decide the sources of cosmic rays at the highest energies.

In recent years, ground-based gamma-ray observatories have made a number of important astrophysical discoveries which have attracted the attention of the wider scientific community. The Division of Astrophysics of the American Physical Society has requested the preparation of a white paper on the status and future of ground-based gamma-ray astronomy to define the science goals of the future observatory, to determine the performance specifications, and to identify the areas of necessary technology development. In this contribution we give a brief overview of the activities of the current white paper team and invite the international community to contribute to the white paper.

Relativistic black-hole jet sources are leading candidates for high energy (>> TeV) neutrino production. The relations defining (a) efficient photopion losses of cosmic-ray protons on target photons and (b) gamma-gamma opacity of gamma rays through that same target photon field imply clear multiwavelength predictions for when and at what energies blazars and GRBs should be most neutrino bright and gamma-ray dim. The use of multiwavelength observations to test the standard relativistic jet model for these source is illustrated.

One way to understand the nonthermal history of the universe is by establishing the origins of the unresolved and truly diffuse extragalactic gamma rays. Dim blazars and radio/gamma galaxies certainly make an important contribution to the galactic gamma-ray background given the EGRET discoveries, and previous treatments are reviewed and compared with a new analysis. Studies of the gamma-ray intensity from cosmic rays in star-forming galaxies and from structure formation shocks, as well as from dim GRBs, are briefly reviewed. A new hard gamma-ray source class seems required from the predicted aggregate intensity compared with the measured intensity.

GRB X-ray light curves display rapid declines followed by a gradual steepening or plateau phase in >~ 30% of GRBs in the Swift sample. Treating the standard relativistic blastwave model in a uniform circumburst medium, it is shown that if GRBs accelerate ultra-high energy cosmic rays through a Fermi mechanism, then the hadronic component can be rapidly depleted by means of photopion processes on time scales ~100 - 10^4 s after the GRB explosion. While discharging the hadronic energy in the form of ultra-high energy cosmic ray neutrals and escaping cosmic-ray ions, the blast wave goes through a strongly radiative phase, causing the steep declines observed with Swift. Following the discharge, the blast wave recovers its adiabatic behavior, forming the observed plateaus or slow declines. These effects are illustrated by calculations of model bolometric light curves. The results show that steep X-ray declines and plateau features occur when GRB sources take place in rather dense media, with n >~ 100 cc out to >~ 10^17 cm.

We summarize our model for long-duration gamma ray bursts (GRBs) that fits the redshift (z) distributions measured with Swift and missions before Swift, and the pre-Swift GRB jet opening-angle distribution inferred from achromatic breaks in the optical light curves. We find that the comoving rate density of GRB sources exhibits positive evolution to $z \gtrsim 3$ -- 5, whereas the star formation rate inferred from measurements of the blue and UV luminosity density peaks at $z\sim 1$ -- 3. The mean intrinsic beaming factor of GRBs is found to be $\approx 34$ -- 42, and we predict that the mean GRB optical jet opening half-angle measured with Swift is $\approx 10^\circ$. We estimate the number of GRBs per year that GLAST is expected to observe based on ratios of BATSE and EGRET GRB fluences.

A simple physical model for long-duration gamma ray bursts (GRBs) is used to fit the redshift (z) and the jet opening-angle distributions measured with earlier GRB missions and with Swift. The effect of different sensitivities for GRB triggering is sufficient to explain the difference in the z distributions of the pre-Swift and Swift samples, with mean redshifts of <z> ~1.5 and <z>~2.7, respectively. Assuming that the emission properties of GRBs do not change with time, we find that the data can only be fitted if the comoving rate-density of GRB sources exhibits positive evolution to z >~ 3-5. The mean intrinsic beaming factor of GRBs is found to range from ~34-42, with the Swift average opening half-angle <\theta_j> ~10 degree, compared to the pre-Swift average of <\theta_j> ~7 degree. Within the uniform jet model, the GRB luminosity function is proportional to L^{-3.25}_*, as inferred from our best fit to the opening angle distribution. Because of the unlikely detection of several GRBs with z <~ 0.25, our analysis indicates that low redshift GRBs represent a different population of GRBs than those detected at higher redshifts. Neglecting possible metallicity effects on GRB host galaxies, we find that ~1 GRB occurs every 600,000 yrs in a local L_* spiral galaxy like the Milky Way. The fraction of high-redshift GRBs is estimated at 8-12% and 2.5-6% at z >= 5 and z >= 7, respectively, assuming continued positive evolution of the GRB rate density to high redshifts.

A study of the statistics of cosmological black-hole jet sources is applied to EGRET blazar data, and predictions are made for GLAST. Black-hole jet sources are modeled as collimated relativistic plasma outflows with radiation beamed along the jet axis due to strong Doppler boosting. The comoving rate density of blazar flares is assumed to follow a blazar formation rate (BFR), modeled by analytic functions based on astronomical observations and fits to EGRET data. The redshift and size distributions of gamma-ray blazars observed with EGRET, separated into BL Lac object (BL) and flat spectrum radio quasar (FSRQ) distributions, are fit with monoparametric functions for the distributions of the jet Lorentz factor \Gamma, comoving directional power l'_e, and spectral slope. A BFR factor ~10 x greater at z ~ 1 than at present is found to fit the FSRQ data. A smaller comoving rate density and greater luminosity of BL flares at early times compared to the present epoch fits the BL data. Based on the EGRET observations, ~1000 blazars consisting of ~800 FSRQs and FR2 radio galaxies and ~200 BL Lacs and FR1 radio galaxies will be detected with GLAST during the first year of the mission. Additional AGN classes, such as hard-spectrum BL Lacs that were mostly missed with EGRET, could add more GLAST sources. The FSRQ and BL contributions to the EGRET gamma-ray background at 1 GeV are estimated at the level of ~10 - 15% and ~2 - 4%, respectively. EGRET and GLAST sensitivities to blazar flares are considered in the optimal case, and a GLAST analysis method for blazar detection is outlined.

Our knowledge of the high-energy universe is undergoing a period of rapid change as new astronomical detectors of high-energy radiation start to operate at their design sensitivities. Now is a boomtime for high-energy astrophysics, with new discoveries from Swift and HESS, results from MAGIC and VERITAS starting to be reported, the upcoming launches of the gamma-ray space telescopes GLAST and AGILE, and anticipated data releases from IceCube and Auger. A formalism for calculating statistical properties of cosmological gamma-ray sources is presented. Application is made to model calculations of the statistical distributions of gamma-ray and neutrino emission from (i) beamed sources, specifically, long-duration GRBs, blazars, and extagalactic microquasars, and (ii) unbeamed sources, including normal galaxies, starburst galaxies and clusters. Expressions for the integrated intensities of faint beamed and unbeamed high-energy radiation sources are also derived. A toy model for the background intensity of radiation from dark-matter annihilation taking place in the early universe is constructed. Estimates for the gamma-ray fluxes of local group galaxies, starburst, and infrared luminous galaxies are briefly reviewed. Because the brightest extragalactic gamma-ray sources are flaring sources, and these are the best targets for sources of PeV -- EeV neutrinos and ultra-high energy cosmic rays, rapidly slewing all-sky telescopes like MAGIC and an all-sky gamma-ray observatory beyond Milagro will be crucial for optimal science return in the multi-messenger age.

Likely astrophysical sources of detectable high-energy (>> TeV) neutrinos are considered. Based on gamma-ray emission properties, the most probable sources of neutrinos are argued to be GRBs, blazars, microquasars, and supernova remnants. Diffuse neutrino sources are also briefly considered.

Highly variable gamma-ray pulses and X-ray flares in GRB light curves can result from external shocks rather than central engine activity under the assumption that the GRB blast-wave shell does not spread. Acceleration of cosmic rays to 10^{20} eV energies can take place in the external shocks of GRBs. Escape of hadronic energy in the form of UHECRs leads to a rapidly decelerating GRB blast wave, which may account for the rapid X-ray declines observed in Swift GRBs.

The Advanced Compton Telescope (ACT), the next major step in gamma-ray astronomy, will probe the fires where chemical elements are formed by enabling high-resolution spectroscopy of nuclear emission from supernova explosions. During the past two years, our collaboration has been undertaking a NASA mission concept study for ACT. This study was designed to (1) transform the key scientific objectives into specific instrument requirements, (2) to identify the most promising technologies to meet those requirements, and (3) to design a viable mission concept for this instrument. We present the results of this study, including scientific goals and expected performance, mission design, and technology recommendations.

Gamma-ray bursts (GRBs) are sources of energetic, highly variable fluxes of gamma rays, which demonstrates that they are powerful particle accelerators. Besides relativistic electrons, GRBs should also accelerate high-energy hadrons, some of which could escape cooling to produce ultra-high energy cosmic rays (UHECRs). Acceleration of high-energy hadrons in GRB blast waves will be established if high-energy neutrinos produced through photopion interactions in the blast wave are detected from GRBs. Limitations on the energy in nonthermal hadrons and the number of expected neutrinos are imposed by the fluxes from pair-photon cascades initiated in the same processes that produce neutrinos. Only the most powerful bursts at fluence levels >~ 3e-4 erg/cm^2 offer a realistic prospect for detection of >> TeV neutrinos. Detection of high-energy neutrinos is likely if GRB blast waves have large baryon loads and Doppler factors <~ 200. Cascade gamma rays will accompany neutrino production and might already have been detected as anomalous emission components in the spectra of some GRBs. Prospects for detection of GRBs in the Milky Way are also considered.

Recent HESS observations show that microquasars in high-mass systems are sources of VHE gamma-rays. A leptonic jet model for microquasar gamma-ray emission is developed. Using the head-on approximation for the Compton cross section and taking into account angular effects from the star's orbital motion, we derive expressions to calculate the spectrum of gamma rays when nonthermal jet electrons Compton-scatter photons of the stellar radiation field. Calculations are presented for power-law distributions of nonthermal electrons that are assumed to be isotropically distributed in the comoving jet frame, and applied to $\gamma$-ray observations of LS 5039. We conclude that (1) the TeV emission measured with HESS cannot result only from Compton-scattered stellar radiation (CSSR), but could be synchrotron self-Compton (SSC) emission or a combination of CSSR and SSC; (2) fitting both the HESS data and the EGRET data associated with LS 5039 requires a very improbable leptonic model with a very hard electron spectrum. Because the gamma rays would be variable in a leptonic jet model, the data sets are unlikely to be representative of a simultaneously measured gamma-ray spectrum. We therefore attribute EGRET gamma rays primarily to CSSR emission, and HESS gamma rays to SSC emission. Detection of periodic modulation of the TeV emission from LS 5039 would favor a leptonic SSC or cascade hadron origin of the emission in the inner jet, whereas stochastic variability alone would support a more extended leptonic model. The puzzle of the EGRET gamma rays from LS 5039 will be quickly solved with GLAST. (Abridged)

Stochastic acceleration of charged particles due to interactions with magnetohydrodynamic (MHD) plasma waves is the dominant process leading to the formation of the high-energy electron and ion distributions in a variety of astrophysical systems. Collisions with the waves influence both the energization and the spatial transport of the particles, and therefore it is important to treat these two aspects of the problem in a self-consistent manner. We solve the representative Fokker-Planck equation to obtain a new, closed-form solution for the time-dependent Green's function describing the acceleration and escape of relativistic ions interacting with Alfven or fast-mode waves characterized by momentum diffusion coefficient $D(p)\propto p^q$ and mean particle escape timescale $t_esc(p) \propto p^{q-2}$, where $p$ is the particle momentum and $q$ is the power-law index of the MHD wave spectrum. In particular, we obtain solutions for the momentum distribution of the ions in the plasma and also for the momentum distribution of the escaping particles, which may form an energetic outflow. The general features of the solutions are illustrated via examples based on either a Kolmogorov or Kraichnan wave spectrum. The new expressions complement the results obtained by Park and Petrosian, who presented exact solutions for the hard-sphere scattering case ($q=2$) in addition to other scenarios in which the escape timescale has a power-law dependence on the momentum. Our results have direct relevance for models of high-energy radiation and cosmic-ray production in astrophysical environments such as $\gamma$-ray bursts, active galaxies, and magnetized coronae around black holes.

The accretion of ~ 0.1 -- 1 M_o of material by a neutron star through Roche lobe overflow of its companion or through white-dwarf/neutron-star coalescence in a low mass binary system could be enough to exceed the critical mass of a neutron star and trigger its collapse to a black hole, leading to the production of a short gamma-ray burst (SGRB). In this model, SGRBs would often be found in early-type galaxies or in globular cluster environments, though they could also be formed through stellar-wind accretion in high-mass binary systems of star-forming galaxies. Although the collapse event is likely to be of short (<< 1 s) duration, afterglow complexity could be produced by neutron decoupling in the fireball and subsequent accretion of the companion remnant, for example, the disrupted white-dwarf core.

Jet models for the high-energy emission of Galactic X-ray binary sources have regained significant interest with detailed spectral and timing studies of the X-ray emission from microquasars, the recent detection by the HESS collaboration of very-high-energy gamma-rays from the microquasar LS~5039, and the earlier suggestion of jet models for ultraluminous X-ray sources observed in many nearby galaxies. Here we study the synchrotron and Compton signatures of time-dependent electron injection and acceleration, adiabatic and radiative cooling, and different jet geometries in the jets of Galactic microquasars. Synchrotron, synchrotron-self-Compton, and external-Compton radiation processes with soft photons provided by the companion star and the accretion disk are treated. An analytical solution is presented to the electron kinetic equation for general power-law geometries of the jets for Compton scattering in the Thomson regime. We pay particular attention to predictions concerning the rapid flux and spectral variability signatures expected in a variety of scenarios, making specific predictions concerning possible spectral hysteresis, similar to what has been observed in several TeV blazars. Such predictions should be testable with dedicated monitoring observations of Galactic microquasars and ultraluminous X-ray sources using Chandra and/or XMM-Newton.

A numerical simulation model of the injection and cooling of nonthermal particles energized by shocks formed in merging clusters of galaxies is used to fit radio and X-ray data observed from the Coma cluster of galaxies. The results are consistent with a primary merger-shock origin for both the diffuse radio halo emission and the hard X-ray excess measured with {\it Beppo-SAX} and {\it RXTE}. For equal (1%) efficiency of power injected in nonthermal protons and electrons, we predict that the Coma cluster of galaxies will be significantly detected with the space-based observatory {\it GLAST}, and marginally detectable with the ground-based $\gamma$-ray observatories {\it VERITAS} and {\it HESS}. Significant TeV detections are possible if the nonthermal proton intensity is greater due to a larger efficiency of nonthermal hadron acceleration, or to past merger events. The nonthermal hadronic content in Coma is also revealed by a weak, hard secondary emission component at $\sim 10$ -- 100 GHz. The spectral softening of the radio emission at large radii from the Coma cluster core derives in this scenario from the decreasing shock speed away from cluster center for an on-axis merger event. We discuss differences between merger and delayed turbulence models.

Gamma-ray bursts are known to be sources of high-energy gamma rays, and are likely to be sources of high-energy cosmic rays and neutrinos. Following a short review of observations of GRBs at multi-MeV energies and above, the physics of leptonic and hadronic models of GRBs is summarized. Evidence for two components in BATSE and EGRET/TASC data suggest that GRBs are sources of high-energy cosmic rays. GLAST observations will reveal the high-energy gamma-ray power and energy releases from GRBs, and will provide detailed knowledge of anomalous high-energy emission components, but confirmation of cosmic ray acceleration must await 100 TeV -- PeV neutrino detection from GRBs.

We analyze the constraint equation giving allowed solutions describing fields and currents in a force-free magnetosphere around a rotating black hole. Utilizing the divergence properties of the energy and angular-momentum fluxes for physically allowed solutions, we conclude that poloidal surfaces are independent of the radial coordinate for large values of $r$. Using this fact and the Znajek regularity condition, we explicitly derive all possible exact solutions admitted by the constraint equation for $r$ independent poloidal surfaces, which are given in terms of the electromagnetic angular velocity function $\Omega = 1/a\sin^2\theta$, where $a$ is the angular momentum per unit mass of the black hole.

Very high energy (VHE) gamma-rays have recently been detected from the Galactic black-hole candidate and microquasar LS 5039. A plausible site for the production of these VHE gamma-rays is the region close to the mildly relativistic outflow. However, at distances comparable to the binary separation, the intense photon field of the stellar companion will lead to substantial gamma-gamma absorption of VHE gamma-rays. If the system is viewed at a substantial inclination (i > 0), this absorption feature will be modulated on the orbital period of the binary as a result of a phase-dependent stellar-radiation intensity and pair-production threshold. We apply our results to LS 5039 and find that (1) gamma-gamma absorption effects will be substantial if the photon production site is located at a distance from the central compact object of the order of the binary separation (~ 2.5e12 cm) or less; (2) the gamma-gamma absorption depth will be largest at a few hundred GeV, leading to a characteristic absorption trough; (3) the gamma-gamma absorption feature will be strongly modulated on the orbital period of the binary, characterized by a spectral hardening accompanying periodic dips of the VHE gamma-ray flux; and (4) gamma rays can escape virtually unabsorbed, even from within ~ 10^{12} cm, when the star is located behind the production site as seen by the observer.

The Blandford-Znajek constraint equation for a stationary, axisymmetric black-hole force-free magnetosphere is cast in a 3+1 absolute space and time formulation, following Komissarov (2004). We derive an analytic solution for fields and currents to the constraint equation in the far-field limit that satisfies the Znajek condition at the event horizon. This solution generalizes the Blandford-Znajek monopole solution for a slowly rotating black hole to black holes with arbitrary angular momentum. Energy and angular momentum extraction through this solution occurs mostly along the equatorial plane. We also present a nonphysical, reverse jet-like solution.

We have developed a model that explains cosmic rays with energies E between \~0.3 PeV and the energy of the second knee at E_2 ~ 3*10^{17} eV as originating from a recent Galactic gamma-ray burst (GRB) that occurred ~1 Myr ago within 1 kpc from Earth. Relativistic shocks from GRBs are assumed to inject power-law distributions of cosmic rays (CRs) to the highest energies. Diffusive propagation of CRs from the local GRB explains the CR spectrum near and above the first knee at E_1 ~ 3*10^{15} eV. The first and the second knees are explained as being directly connected with the injection of plasma turbulence in the interstellar medium on a ~1 pc and ~100 pc scales, respectively. Transition to CRs from extragalactic GRBs occurs at E > E_2. The origin of the ankle in the CR spectrum at E ~ 4*10^{18} eV is due to photopair energy losses of UHECRs on cosmological timescales, as also suggested by Berezinsky and collaborators. Any significant excess flux of extremely high energy CRs deviating from the exponential cutoff behavior at E> E_{GZK} = 6*10^{19} eV would imply a significant contribution due to recent GRB activity on timescales t < 10^8 yrs from local extragalactic sources within ~10 Mpc.

A complete model for the origin of high-energy >~10^{14} eV) cosmic rays from gamma-ray bursts (GRBs) and implications of this hypothesis are described. Detection of high-energy neutrinos from GRBs provide an unambiguous test of the model. Evidence for cosmic-ray acceleration in GRBs is suggested by the detection of anomalous gamma-ray components such as that observed from GRB 941017. Neutron beta-decay halos around star-forming galaxies such as the Milky Way are formed as a consequence of this model. Cosmic rays from GRBs in the Galaxy are unlikely to account for the ~10^{18} eV cosmic-ray excess reported by the Sydney University Giant Air Shower Recorder (SUGAR), but could contribute to past extinction events.

The rate of terrestrial irradiation events by galactic gamma-ray bursts (GRBs) is estimated using recent standard-energy results. We assume that GRBs accelerate high-energy cosmic rays, and present results of three-dimensional simulations of cosmic rays moving in the Galactic magnetic field and diffusing through pitch-angle scattering. An on-axis GRB extinction event begins with a powerful prompt gamma-ray and neutron pulse, followed by a longer-lived phase from cosmic-ray protons and neutron-decay protons that diffuse towards Earth. Our results force a reinterpretation of reported ~ 10^{18} eV cosmic-ray anisotropies and offer a rigorous test of the model where high-energy cosmic rays originate from GRBs, which will soon be tested with the Auger Observatory.

The HESS collaboration recently reported highly significant detection of TeV gamma-rays coincident with Sgr A*. In the context of other Galactic Center (GC) observations, this points to the following scenario: In the extreme advection-dominated accretion flow (ADAF) regime of the GC black hole (BH), synchrotron radio/sub-mm emission of 100 MeV electrons emanates from an inefficiently radiating turbulent magnetized corona within 20 Schwarzschild radii of the GCBH. These electrons are accelerated through second-order Fermi processes by MHD turbulence. Closer to the innermost stable orbit of the ADAF, instabilities and shocks within the flow inject power-law electrons through first-order Fermi acceleration to make synchrotron X-ray flares observed with Chandra, XMM, and INTEGRAL. A subrelativistic MHD wind subtending a 1sr cone with power > 10^{37} erg/s is driven by the ADAF from the vicinity of the GCBH. As in pulsar powered plerions, electrons are accelerated at the wind termination shock, at > 10^{16.5} cm from the GCBH, and Compton-scatter the ADAF and the far infra-red (FIR) dust radiation to TeV energies. The synchrotron radiation of these electrons forms the quiescent X-ray source resolved by Chandra. The radio counterpart of this TeV/X-ray plerion, formed when the injected electrons cool on timescales 10^4 yrs, could explain the origin of nonthermal radio emission in the pc-scale bar of the radio nebula Sgr A West.

Analytic approximations for synchrotron, synchrotron self-Compton (SSC), and external Compton (EC) processes are used to constrain model parameters for knot and hot-spot emission in extended jets of radio galaxies. Equipartition formulas are derived that relate the Doppler factor $\delta$ and comoving magnetic field $B$ assuming a nonthermal synchrotron origin of the radio emission, and synchrotron, SSC and EC origins of the X-ray emission. Expressions are also derived for $\delta$ and $B$ that minimize the total jet powers of the emitting region in synchrotron, SSC and EC models for the X-ray emission. The results are applied to knot WK7.8 of PKS 0637-752. Predictions to test two-component synchrotron and EC models are made for {\it Chandra} and {\it GLAST}.

An elementary kinematic model for emission produced by relativistic spherical colliding shells is studied. The case of a uniform blast-wave shell with jet opening angle $\theta_j \gg 1/\Gamma$ is considered, where $\Gamma$ is the Lorentz factor of the emitting shell. The shell, with comoving width $\Delta r^\prime$, is assumed to be illuminated for a comoving time $\Delta t^\prime$ and to radiate a broken power-law $\nu L_\nu$ spectrum peaking at comoving photon energy $\e_{pk,0}^{\prime}$. Synthetic GRB pulses are calculated, and the relation between energy flux and internal comoving energy density is quantified. Curvature effects dictate that the measured $\nu F_\nu$ flux at the measured peak photon energy $\e_{pk}$ is proportional to $\e^3_{pk}$ in the declining phase of a GRB pulse. Possible reasons for discrepancy with observations are discussed, including adiabatic and radiative cooling processes that extend the decay timescale, a nonuniform jet, or the formation of pulses by external shock processes. A prediction of a correlation between prompt emission properties and times of the optical afterglow beaming breaks is made for a cooling model, which can be tested with Swift.

A widely discussed explanation for the origin of the X-ray emission observed from knots in extended quasar jets with the Chandra X-ray Obseratory is Compton-scattered CMBR by electrons with Lorentz factors $\gamma^\prime \sim 10^2$. This model faces difficulties in terms of total energy requirements, and in explaining the spatial profiles of the radio, optical, and X-ray knots in sources such as PKS 0637-752, 3C 273, or PKS 1127-145. These difficulties can be resolved in the framework of one- and two-component synchrotron models. We propose a model where the broad band radio to X-ray synchrotron emission in quasar jets is powered by collimated beams of ultra-high energy neutrons and gamma-rays formed in the sub-parsec scale jets. The decay of the neutral beam in the intergalactic medium drives relativistic shocks to accelerate nonthermal electrons out of the ambient medium. A second synchrotron component arises from the injection of leptons with Lorentz factors $\gg 10^7$ that appear in the extended jet in the process of decay of ultra-high energy gamma rays. This approach could account for qualitative differences in the extended X-ray jets of FR1 and FR2 galaxies. Detection of high-energy neutrinos from blazars and core-dominated quasars will provide strong evidence for this model.

We have developed a numerical model for the temporal evolution of particle and photon spectra resulting from nonthermal processes at the shock fronts formed in merging clusters of galaxies. Fermi acceleration is approximated by injecting power-law distributions of particles during a merger event, subject to constraints on maximum particle energies. We consider synchrotron, bremsstrahlung, Compton, and Coulomb processes for the electrons, nuclear, photomeson, and Coulomb processes for the protons, and knock-on electron production during the merging process. The broadband radio through $\gamma$-ray emission radiated by nonthermal protons and primary and secondary electrons is calculated both during and after the merger event. To test the ability of the computer model to accurately calculate the nonthermal emission expected from a cluster merger event, we apply the model to the Coma cluster of galaxies, and show that the centrally located radio emission and the Hard X-ray excess observed at 40-80$\kev$ is well fit by our model. If our model is correct, then the Coma cluster will be significantly detected with GLAST and ground-based air Cherenkov telescopes.

This rapporteur review summarizes results presented in Parallel Session GBT2 (Gamma Ray Burst Theory 2) on the Collapsar and Supranova Models held 25 July 2003 at the 10th Marcel Grossmann Meeting on General Relativity, Rio de Janeiro, Brazil. A central issue in GRB studies is the process whereby energy is released from the GRB engine. One scenario is the collapsar model, where the evolved stellar core promptly collapses to a black hole surrounded by a massive, intermittently-accreting torus of nuclear density material. A second scenario is the supranova model, where the first step of a two-step collapse process leaves behind a rapidly rotating neutron star stabilized by rotation, which later collapses to a black hole while making the GRB. In the supranova model, a powerful pulsar wind lasting days to weeks after the supernova makes distinctive signatures from the heating of the supernova remnant (SNR) shell, which should be discovered with Swift. This model also predicts nonstandard reddened excesses from the cooling SNR due to the range of delays between the two collapse events, contrary to observations of SN emissions from low redshift GRBs which favor the collapsar model. The observational basis for both models is critically reviewed. An outline of a GRB model is proposed where the second collapse takes place within minutes to hours after the primary SN event, during which time loss of centrifugal support along the rotation axis of the neutron star provides a relatively baryon-clean polar environment along which a newly formed, rapidly spinning black hole drives collimated baryon-dilute outflows.

Gonz\'alez et al. (2003) have reported the discovery of an anomalous radiation component from ~ 1 -- 200 MeV in GRB 941017. This component varies independently of and contains > 3 times the energy found in the prompt ~ 50 keV -- 1 MeV radiation component that is well described by the relativistic synchrotron-shock model. Acceleration of hadrons to very high energies can give rise to two additional emission components, one produced inside the GRB blast wave and one associated with an escaping beam of ultra-high energy (UHE; > 10^{14} eV) neutrons, gamma rays, and neutrinos. The first component extending to ~ 100 MeV is from a pair-photon cascade induced by photomeson processes with the internal synchrotron photons coincident with the prompt radiation. The outflowing UHE neutral beam can undergo further interactions with external photons from the backscattered photon field to produce a beam of hyper-relativistic electrons that lose most of their energy during a fraction of a gyroperiod in the assumed Gauss-strength magnetic fields of the circumburst medium. The synchrotron radiation of these electrons has a spectrum with vF_v index equal to +1 that can explain the anomalous component in GRB 941017. This interpretation of the spectrum of GRB 941017 requires a high baryon load of the accelerated particles in GRB blast waves. It implies that most of the radiation associated with the anomalous component is released at > 500 MeV, suitable for observations with GLAST, and with a comparable energy fluence in ~100 TeV neutrinos that could be detected with a km-scale neutrino telescope like IceCube.

This rapporteur review covers selected results presented in the Parallel Session HEA2 (High Energy Astrophysics 2) of the 10th Marcel Grossmann Meeting on General Relativity, held in Rio de Janeiro, Brazil, July 2003. The subtopics are: ultra high energy cosmic ray anisotropies, the possible connection of these energetic particles with powerful gamma ray bursts, and new exciting scenarios with a strong neutrino-nucleon interaction in the atmosphere.

We discuss the fluxes of high energy neutrinos and gamma-rays expected from AGNs if hadrons can be effectively accelerated to ultra-high energies by their relativistic jets, as currently believed. Fluxes of multi-TeV neutrinos detectable by $km$-scale detectors like IceCube could be expected from powerful blazars where strong accretion-disk radiation is present in the AGN cores. Gamma-ray fluxes of hadronic origin can be important for flares in the compact jets of these sources up to GeV energies, but they will be insignificant for BL Lac objects. Production of UHE neutral beams composed of neutrons and gamma-rays can drive straight collimated jets in the intergalactic medium on multi-kpc scales, which could be resolved by the Chandra X-ray observatory. While we do not expect any significant neutrino flux from these large-scale jets, we predict gamma-rays of synchrotron origin in the energy range from sub-GeV up to TeV energies, which would be detectable by GLAST and ground-based gamma-ray telescopes.

A derivation of the \gamma\gamma --> e^+ e^- optical depth for \gamma rays produced in a comoving spherical emitting region is presented. Employing a simplified expression for the \gamma\gamma absorption cross section, analytic expressions for the minimum Doppler factor implied by the requirement of gamma-ray transparency are derived for a broken power-law spectrum of target photons which are isotropically distributed in the comoving frame. Application to specific systems is illustrated.

A model is proposed for the origin of cosmic rays (CRs) from ~10^14 eV to the highest energies, >10^20 eV. Gamma-Ray Bursts (GRBs) are assumed to inject CR protons and ions into the interstellar medium of star-forming galaxies--including the Milky Way--with a power law spectrum extending to a maximum energy ~10^20 eV. The CR spectrum near the knee is fit with CRs trapped in the Galactic halo that were accelerated and injected by an earlier Galactic GRB. These CRs diffuse in the disk and halo of the Galaxy due to gyroresonant pitch-angle scattering with MHD turbulence in the Galaxy's magnetic field. The preliminary (2001) KASCADE data through the knee of the CR spectrum are fit by a model with energy-dependent propagation of CR ions from a single Galactic GRB. Ultra-high energy CRs (UHECRs), with energies above the ankle are assumed to propagate rectilinearly with their spectrum modified by photo-pion, photo-pair, and expansion losses. We fit the measured UHECR spectrum assuming comoving luminosity densities of GRB sources consitent with possible star formation rate histories of the universe. For power-law CR proton injection p>2 this model implies that the nonthermal content in the GRB blast waves is hadronically dominated by a factor ~60-200. Calculations show that 100 TeV-100 PeV neutrinos could be detected several times per year from all GRBs in kilometer-scale neutrino detectors such as IceCube, for GRB blast-wave Doppler factors <~200. GLAST measurements of gamma-ray components and cutoffs will constrain the product of nonthermal baryon loading and radiative efficiency, limit the Doppler factor, and test this senario.

We derive the minimum observing time scales to detect a blazar at a given flux level with the LAT on GLAST in the scanning and pointing modes. Based upon Phase 1 observations with EGRET, we predict the GLAST detection rate of blazar flares at different flux levels. With some uncertainty given the poor statistics of bright blazars, we predict that a blazar flare with integral flux >~ 200e-8 ph(> 100 MeV) cm^{-2} s^{-1}, which are the best candidates for Target of Opportunity pointings and extensive temporal and spectral studies, should occur every few days.

Several lines of evidence point to a relationship between gamma-ray bursts (GRBs) and the high mass stars that explode as supernovae. Arguments that GRB sources accelerate cosmic rays (CRs) are summarized. High-energy neutrino detection from individual GRBs would mean that they are hadronically dominated, that is, that the amount of energy deposited in the form of nonthermal hadrons in GRB blast waves greatly exceeds the radiated energy inferred from their photon emission spectra. Such a detection would make GRBs the favored candidate sources of ultra-high energy and super-GZK CRs. Cascade radiation induced by high-energy hadrons could be detected from GRBs with gamma-ray telescopes.

Gonz\'alez et al. (2003) have reported the discovery of an anomalous radiation component from ~ 1 -- 200 MeV in GRB 941017. This component varies independently of and contains >~ 3 times the energy found in the prompt ~ 50 keV -- 1 MeV radiation component that is well described by the relativistic synchrotron-shock model. Acceleration of hadrons to very high energies by GRBs could give rise to a separate emission component. Two models, both involving acceleration of ultra-high energy cosmic rays with subsequent photomeson interactions, are considered. The first involves a pair-photon cascade initiated by photohadronic processes in the GRB blast wave. Calculations indicate that the cascade produces a spectrum that is too soft to explain the observations. A second model is proposed where photopion interactions in the GRB blast-wave shell give rise to an escaping collimated neutron beam. The outflowing neutrons undergo further photopion interactions to produce a beam of hyper-relativistic electrons that can lose most of their energy during a fraction of a gyroperiod in the Gauss-strength magnetic fields found in the circumburst medium. This secondary electron beam produces a hard synchrotron radiation spectrum that could explain the anomalous component in GRB 941017.

Acceleration of high-energy hadrons in GRB blast waves will be established if high-energy neutrinos are detected from GRBs. Recent calculations of photomeson neutrino production are reviewed, and new calculations of high-energy neutrinos and the accompanying hadronic cascade radiation are presented. If hadrons are injected in GRB blast waves with an energy corresponding to the measured hard X-ray/soft gamma-ray emission, then only the most powerful bursts at fluence levels >~ 3e-4 erg cm^{-2} offer a realistic prospect for detection of muon neutrinos. Detection of high-energy neutrinos are likely if GRB blast waves have large baryon loads and Doppler factors <~ 200. Significant limitations on the hadronic baryon loading and the number of expected neutrinos are imposed by the fluxes from pair-photon cascades initiated in the same processes that produce neutrinos.

A model for high-energy (>10^14 eV) cosmic rays (HECRs) from galactic and extragalactic gamma-ray bursts (GRBs) is summarized. Relativistic outflows in GRBs are assumed to inject power-law distributions of CR protons and ions to the highest (>10^20 eV) energies. A diffusive propagation model for HECRs from a single recent GRB within 1 kpc from Earth explains the CR spectrum near and above the knee. The CR spectrum at energies above 10^18 eV is fit with a component from extragalactic GRBs. By normalizing the energy injection rate to that required to produce the CR flux from extragalactic sources observed locally, we determine the amount of energy a typical GRB must release in the form of nonthermal hadrons. Our interpretation of the HECR spectrum requires that GRBs are hadronically dominated, which would be confirmed by the detection of HE neutrinos from GRBs.

We treat high-energy neutrino production in GRBs. Detailed calculations of photomeson neutrino production are presented for the collapsar model, where internal nonthermal synchrotron radiation is the primary target photon field, and the supranova model, where external pulsar-wind synchrotron radiation provides important additional target photons. Detection of > 10 TeV neutrinos from GRBs with Doppler factors > 200, inferred from gamma-ray observations, would support the supranova model. Detection of < 10 TeV neutrinos is possible for neutrinos formed from nuclear production. Only the most powerful bursts at fluence levels > 3x10^{-4} erg/cm^2 offer a realistic prospect for detection of muon neutrinos.

We have developed a numerical model for the temporal evolution of particle and photon spectra resulting from nonthermal processes at the shock fronts formed in merging clusters of galaxies. Fermi acceleration is approximated by injecting power-law distributions of particles during a merger event, subject to constraints on maximum particle energies. We consider synchrotron, bremsstrahlung, Compton, and Coulomb processes for the electrons, nuclear, photomeson, and Coulomb processes for the protons, and knock-on electron production during the merging process. The broadband radio through gamma-ray emission radiated by nonthermal protons and primary and secondary electrons is calculated both during and after the merger event. Using ROSAT observations to establish typical parameters for the matter density profile of clusters of galaxies, we find that typical merger shocks are weak and accelerate particles with relatively soft spectra. We consider the prospects for detecting nonthermal radio and gamma-ray emission from clusters of galaxies and implications for the origin of ultra-high energy cosmic rays and the diffuse gamma-ray background. Our results suggest that only a few of the isotropically-distributed unidentified EGRET sources are due to shocks formed in cluster mergers, and that only a minor contribution to the diffuse extragalactic gamma-ray background can originate from cluster merger shocks. Cluster merger shocks can accelerate protons to $\lesssim 10^{19}$ eV for the standard parameters considered here. We predict that {\it GLAST} will detect several cluster mergers and, depending on the mean magnetic fields in the intracluster medium, the Low Frequency Array could detect anywhere from several to several hundred.

An external shock model for the prompt gamma-ray luminous phase of gamma-ray bursts (GRBs) is treated both analytically and numerically. A widely cited derivation claiming that an external shock model for rapidly variable GRBs must be very inefficient employs an incorrect expression for the angular timescale. Numerical results show that variable GRBs can be formed with >10% efficiency to transform the directed kinetic energy of the relativistic fireball to gamma-rays. Successes of the external shock model and difficulties with an internal wind/colliding shell model are summarized. An impulsive external shock model is consistent with the supranova model.

We treat the production of neutrons, photons, and neutrinos through photomeson interactions of relativistic protons with ambient photons in the compact inner jets of blazars. We show that the presence of the external UV radiation field makes possible strong energy losses already for protons with energies > 1 PeV, while without this component effective energy losses of protons begin only at E > 10^{18} eV. We develop a model describing the production and escape of neutrons from a comoving spherical blob, which continue to interact with the ambient external radiation field on the parsec-scale broad line region (BLR). Neutrons may carry ~10% of the overall energy of the protons accelerated beyond E ~ 1 PeV outside the BLR. Ultra-high energy (UHE) gamma rays produced in photomeson interactions of neutrons outside the blob can also escape the BLR. The escaping neutrons, gamma rays and neutrinos form a collimated neutral beam with a characteristic opening angle ~ 1/Gamma, where Gamma is the bulk Lorentz factor of the inner jet. The energy and the momentum of such beam is then mostly deposited in the extended jet due to neutron decay at distances ~ (E_n/10^{17} eV}) kpc, and through pair-production attenuation of gamma rays with energies E_g ~ 10^{15}-10^{18} eV which can propagate to distances beyond (10-100) kpc. In this scenario, neutral beams of UHE gamma rays and neutrons can be the reason for straight extended jets such as in Pictor A. Fluxes of neutrinos detectable with km-scale neutrino telescopes are predicted from flat spectrum radio quasars such as 3C 279.

We have developed a numerical model for the temporal evolution of particle and photon spectra resulting from nonthermal processes at the shock fronts formed in merging clusters of galaxies. Fermi acceleration is approximated by injecting power-law distributions of particles during a merger event, subject to constraints on maximum particle energies. We consider synchrotron, bremsstrahlung, Compton, and Coulomb processes for the electrons, nuclear, photomeson, and Coulomb processes for the protons, and knock-on electron production. Broadband radio through gamma-ray light curves radiated by nonthermal protons and primary and secondary electrons are calculated both during and after the merger event. Using ROSAT observations to establish typical parameters for the matter density profile of clusters of galaxies, we find that merger shocks are weak and accelerate particles with relatively soft spectra. Our results suggest that only a minor contribution to the diffuse extragalactic gamma-ray background can originate from cluster merger shocks.

Excess reddened optical emission has been measured in the power-law decay afterglow light curves of several GRBs some tens of days after the GRB events. This emission is claimed to be a signature of a supernova (SN) taking place at the same time as the GRB event, in support of the collapsar scenario. Observations of features in prompt and afterglow X-ray spectra of GRBs require high density, high metallicity material near GRB sources, supporting the supranova (SA) scenario. These conflicting lines of evidence are resolved in likely regimes of parameter space if the supernova remnant (SNR) shell is illuminated and heated by a pulsar wind (PW), and cools following the GRB event. The heating of the SNR shell by the PW is parameterized, and simplified expressions for the synchrotron, synchrotron self-Compton, and thermal radiation fields within the SNR shell and pulsar wind bubble are derived for monoenergetic injection of wind electrons. A cooling SNR shell could produce the excess emission detected from GRB 970228, GRB 980326, and GRB 011121. X-ray synchrotron radiation from the PW electrons provides an important source of ionizing radiation.

We treat transformation properties of external radiation fields in the proper frame of a plasma moving with constant speed. The specific spectral energy densities of external isotropic and accretion-disk radiation fields are derived in the comoving frame of relativistic outflows, such as those thought to be found near black-hole jet and gamma-ray burst sources. Nonthermal electrons and positrons Compton-scatter this radiation field, and high-energy protons and ions interact with this field through photomeson and photopair production. We revisit the problem of the Compton-scattered spectrum associated with an external accretion-disk radiation field, and clarify a past treatment by the authors. Simple expressions for energy-loss rates and Thomson-scattered spectra are given for ambient soft photon fields consisting either of a surrounding external isotropic monochromatic radiation field, or of an azimuthally symmetric, geometrically thin accretion-disk radiation field. A model for blazar emission is presented that displays a characteristic spectral and variability behavior due to the presence of a direct accretion-disk component. The disk component and distinct flaring behavior can be bright enough to be detected from flat spectrum radio quasars with {\it GLAST}. Spectral states of blazars are characterized by the relative importance of the accretion-disk and scattered radiation fields and, in the extended jet, by the accretion disk, inner jet, and cosmic microwave background radiation fields.

The hypothesis that ultra-high energy (>~ 10^19 eV) cosmic rays (UHECRs) are accelerated by gamma-ray burst (GRB) blast waves is assumed to be correct. Implications of this assumption are then derived for the external shock model of gamma-ray bursts. The evolving synchrotron radiation spectrum in GRB blast waves provides target photons for the photomeson production of neutrinos and neutrons. Decay characteristics and radiative efficiencies of the neutral particles that escape from the blast wave are calculated. The diffuse high-energy GRB neutrino background and the distribution of high-energy GRB neutrino events are calculated for specific parameter sets, and a scaling relation for the photomeson production efficiency in surroundings with different densities is derived. GRBs provide an intense flux of high-energy neutrons, with neutron-production efficiencies exceeding ~ 1% of the total energy release. The radiative characteristics of the neutron beta-decay electrons from the GRB "neutron bomb" are solved in a special case. Galaxies with GRB activity should be surrounded by radiation halos of ~ 100 kpc extent from the outflowing neutrons, consisting of a nonthermal optical/X-ray synchrotron component and a high-energy gamma-ray component from Compton-scattered microwave background radiation. The luminosity of sources of GRBs and relativistic outflows in L* galaxies such as the Milky Way is at the level of ~10^40+-1 ergs/s. This is sufficient to account for UHECR generation by GRBs. We briefly speculate on the possibility that hadronic cosmic rays originate from the subset of supernovae that collapse to form relativistic outflows and GRBs. (abridged)

In an essay published in 1997, M. J. Harris pointed out that the development of a theoretical explanation, or "synthesis," of a group of astronomical results can lead to a situation where subsequent observations are forced to fit within that framework. Moreover, new results may be reported that fit within the framework but which can no longer be reproduced following the advent of more sensitive observations that overturn the synthesis. Although Harris considered gamma-ray line transients, a similar situation regarding the collapsar model for gamma-ray bursts now holds. In two recent papers posted on astro-ph by J. S. Bloom et al., and by P. A. Price et al., the synthesis effect is illustrated with respect to the collapsar model. These authors dismiss observations that contradict the collapsar model, appeal to "simplicity" to justify their preferred model, and dismiss (or better yet, ignore) competing models with unscientific reasoning. Whether or not the collapsar model is correct, these papers serve as textbook examples of the synthesis effect in science.

Chandra observations of knots and hot spots in spatially resolved X-ray jets of radio galaxies show that the X-ray fluxes often lie above an extrapolation from the radio-to-optical continuum fluxes. We show that combined synchrotron and Compton losses on a single power-law electron injection function can produce a hardening in the electron spectrum at electron Lorentz factors \gamma \~ 2e8/[\Gamma(1+z)] due to KN energy losses on the cosmic microwave background radiation. Here \Gamma is the bulk Lorentz factor of the outflow, and z is the source redshift. This produces a flattening in the spectrum at frequencies >~ 8e16 \delta B_{\mu G}/[\Gamma^2(1+z)^3] Hz, where B_{\mu G} is the magnetic field in the comoving plasma frame in units of micro-Gauss and \delta is the Doppler factor. A single population of synchrotron-emitting electrons may therefore produce the radio-to-X-ray continuum in some radio galaxy knots, such as those in 3C 273.

This paper presents the theoretical basis of the fireball/blast wave model, and some implications of recent results on GRB source models and cosmic-ray production from GRBs. BATSE observations of the prompt gamma-ray luminous phase, and Beppo-SAX and long wavelength afterglow observations of GRBs are briefly summarized. Derivation of spectral and temporal indices of an adiabatic blast wave decelerating in a uniform surrounding medium in the limiting case of a nonrelativistic reverse shock, both for spherical and collimated outflows, is presented as an example of the general theory. External shock model fits for the afterglow lead to the conclusion that GRB outflows are jetted. The external shock model also explains the temporal duration distribution and clustering of peak energies in prompt spectra of long-duration GRBs, from which the redshift dependence of the GRB source rate density can be derived. Source models are reviewed in light of the constant energy reservoir result of Frail et al. that implies a total GRB energy of a few times 10^{51} ergs and an average beaming fraction of ~1/500 of full sky. Paczynski's isotropic hypernova model is ruled out. The Vietri-Stella model two-step collapse process is preferred over a hypernova/collapsar model in view of the X-ray observations of GRBs and the constant energy reservoir result. Second-order processes in GRB blast waves can accelerate particles to ultra-high energies. GRBs may be the sources of UHECRs and cosmic rays with energies above the knee of the cosmic ray spectrum. High-energy neutrino and gamma-ray observations with GLAST and ground-based gamma-ray telescopes will be crucial to test GRB source models.

The recent detection of a transient absorption feature in the prompt emission of GRB 990705 has sparked multiple attempts to fit this feature in terms of photoelectric absorption or resonance scattering out of the line of sight to the observer. However, the physical conditions required to reproduce the observed absorption feature turn out to be rather extreme compared to the predictions of current GRB progenitor models. In particular, strong clumping of ejecta from the GRB progenitor seems to be required. Using detailed 3D hydrodynamic simulations of supernova explosions as a guideline, we have investigated the dynamics and structure of pre-GRB ejecta predicted in various GRB progenitor models. Based on our results, combined with population synthesis studies relevant to the He-merger model, we estimate the probability of observing X-ray absorption features as seen in GRB 990705 to << 1 %. Alternatively, if the supranova model is capable of producing highly collimated long-duration GRBs, it may be a more promising candidate to produce observable, transient X-ray absorption features.

An important radiation field for photomeson neutrino production in blazars is shown to be the radiation field external to the jet. Assuming that protons are accelerated with the same power as electrons and injected with a -2 number spectrum, we predict that km^2 neutrino telescopes will detect about 1-to-several neutrinos per year from flat spectrum radio quasars (FSRQs) such as 3C 279. The escaping high-energy neutron and photon beams transport inner jet energy far from the black-hole engine, and could power synchrotron X-ray jets and FR II hot spots and lobes.

Blazar subclasses ranging from flat-spectrum radio quasars (FSRQs) through low-frequency-peaked BL Lac objects (LBLs) to high-frequency-peaked BL Lac objects (HBLs) exhibit a sequence of increasing spectral hardness with decreasing luminosity that cannot be explained solely by orientation effects. Using an analytic model for the synchrotron, synchrotron-self-Compton, and Compton-scattered external radiation from blazar jets, we propose an evolutionary scenario that links these blazar subclasses in terms of a reduction of the black-hole accretion power with time. As the circumnuclear material accretes to fuel the central engine, less gas and dust is left to scatter accretion-disk radiation and produce an external Compton-scattered component in blazar spectra. This evolutionary trend produces the sequence FSRQ -> LBL -> HBL. Such a scenario may also link radio-loud AGNs with ultraluminous infrared galaxies and optical QSOs, if the latter constitute the high Eddington-ratio epoch of supermassive black hole growth, as suggested by the observed anti-correlation between radio and soft X-ray activity in some Galactic black-hole candidates.

Many classes of gamma-ray sources, such as gamma-ray bursts, blazars, Seyfert galaxies, and galactic black hole sources are surrounded by large amounts of gas and dust. X-rays and gamma-rays that traverse this material will be attenuated by Compton scattering and photoelectric absorption. One signature of an intervening scattering cloud is radiation-hardening by electrons that have been scattered and heated by the incident radiation, as illustrated by a Monte Carlo calculation. Compton scattering provides backscattered photons that will attenuate subsequent gamma rays through \gamma\gamma pair-production processes. We calculate the pair efficiency for a cloud illuminated by gamma-ray burst radiation. An analytic calculation of the flux of X-rays and gamma rays Thomson scattered by an intervening cloud is presented. Illuminated clouds near GRBs will form relativistic plasmas containing large numbers of electron-positron pairs that can be detected within ~1-2 days of the explosion before expanding and dissipating. Localized regions of pair annihilation radiation in the Galaxy could reveal gamma-ray sources embedded in dense clouds, or sites of past GRB explosions.

Subcluster interactions within clusters of galaxies produce shocks that accelerate nonthermal particles. We treat Fermi acceleration of nonthermal electrons and protons by injecting power-law distributions of particles during the merger event, subject to constraints on maximum particle energies. The broadband nonthermal spectrum emitted by accelerated electrons and protons is calculated during and following the subcluster interaction for a standard parameter set. The intensity of gamma-ray emission from primary and secondary processes is calculated and discussed in light of detection capabilities at radio and gamma-ray energies.

Observations of annihilation radiation in the Galaxy are briefly reviewed. We summarize astrophysical mechanisms leading to positron production, and recent estimates for production rates from nova and supernova nucleosynthesis in the Galaxy. The physical processes involved in the production of annihilation radiation in the interstellar medium are described. These include positron thermalization, charge exchange, radiative recombination, and direct annihilation. Calculations of 2\gamma and 3\gamma spectra and the positronium (Ps) fraction due to the annihilation of positrons in media containing H and He at different temperatures and ionization states are presented. Quenching of Ps by high temperature plasmas or dust could account for differences between 0.511 MeV and 3\gamma Ps continuum maps. These results are presented in the context of the potential of INTEGRAL to map sites of annihilation radiation in the Galaxy. Positron production by compact objects is also considered.

The measured spectral energy distribution and variability time scale are used to determine the radiation and magnetic-field energy densities in the relativistic plasma that forms the gamma-ray emitting jet in the blazar 3C 279. Assuming that protons are accelerated as efficiently as electrons to a maximum energy determined by the size and magnetic field of the emitting region, we calculate the emissivity of neutrinos produced by protons that interact with the external radiation field intercepted by the jet. The external radiation field provides the most important target photons for photomeson production of high-energy neutrinos in flat spectrum radio quasars (FSRQs). Because of photomeson interactions with this field, km^2 neutrino telescopes are predicted to detect > 0.1-1 neutrinos per year from blazars such as 3C 279. BL Lac objects are weaker neutrino sources if, as widely thought, their gamma-ray emission is due to Compton-scattered synchrotron (SSC) radiation.

We derive the maximum accelerated particle energy from first-order and second-order Fermi acceleration at nonrelativistic and relativistic shocks for explosions taking place in a uniform surrounding medium. Second-order stochastic processes in relativistic flows are shown to be capable of accelerating cosmic rays to ultra-high energies. Cosmic rays above the knee of the cosmic ray spectrum can be accelerated by the second-order Fermi mechanism in relativistic flows, such as those occuring in gamma-ray bursts and unusual supernovae like SN 1998bw.

We treat the problem of adiabatic losses and stochastic particle acceleration in gamma-ray burst (GRB) blast waves that decelerate by sweeping up matter from an external medium. The shocked fluid is assumed to be represented by a homogeneous expanding shell. The energy lost by nonthermal particles through adiabatic expansion is converted to the bulk kinetic energy of the outflow, permitting the evolution of the bulk Lorentz factor Gamma of the blast wave to be self-consistently calculated. The behavior of the system is shown to reproduce the hydrodynamic self-similar solutions in the relativistic and nonrelativistic limits, and the formalism is applicable to scenarios that are intermediate between the adiabatic and fully radiative regimes. Nonthermal particle energization through stochastic gyroresonant acceleration with magnetic turbulence in the blast wave is treated by employing energy-gain rates and diffusive escape timescales based upon expressions derived in the quasilinear regime. If the magnetic field in the shocked fluid approaches its equipartition value, this process can accelerate escaping particles to >~10^{20} eV energies, consistent with the hypothesis that ultra-high energy cosmic rays (UHECRs) are accelerated by GRB blast waves. Due to particle trapping by the magnetic turbulence, only the highest energy particles can escape during the prompt and afterglow phases of a GRB for acceleration by a Kolmogorov spectrum of MHD turbulence. Lower energy particles begin to escape as the blast wave becomes nonrelativistic and shock Fermi acceleration becomes more important.

Simple arguments concerning power and acceleration efficiency show that ultra-high energy cosmic rays (UHECRS) with energies >~ 10^{19} eV could originate from GRBs. Neutrons formed through photo-pion production processes in GRB blast waves leave the acceleration site and travel through intergalactic space, where they decay and inject a very energetic proton and electron component into intergalactic space. The neutron-decay protons form a component of the UHECRs, whereas the neutron-decay electrons produce optical/X-ray synchrotron and gamma radiation from Compton-scattered background radiation. A significant fraction of galaxies with GRB activity should be surrounded by neutron-decay halos of characteristic size ~ 100 kpc.

Statistical arguments show that the volume- and time-averaged kinetic power of GRBs and fireball transients (FTs) into an L* galaxy like the Milky Way is at the level of 10^40 ergs/s. This number, though with wide uncertainties related to the internal or external shock efficiency, is sufficient to power hadronic cosmic rays observed locally. The release of energy by the high-mass progenitor stars of GRBs and FTs is sufficient to power the ultrahigh energy cosmic rays, as already shown by Waxman and Vietri in 1995. It is sufficient to power the cosmic rays above the knee of the cosmic-ray spectrum. Indeed, all hadronic cosmic rays could originate from the high-mass (>~ 100 M_o) stars that collapse to black holes, in the process forming GRBs and FTs. This source class represents a new solution to the problem of cosmic-ray origin. The ~10^4 - 10^7 black holes made by these stars could make their presence known by radiating as they accrete from the ISM, by microlensing background radiations, and by forming luminous binary systems. Some unidentified EGRET sources could be isolated black holes that accrete from the ISM. Better imaging and sensitivity with GLAST and TeV observatories will test this model for the unidentified gamma-ray sources, and this theory for cosmic-ray origin.

Low frequency electromagnetic turbulence is generated in relativistically outflowing plasma that sweeps up particles from the surrounding environment. Electrons are energized by stochastic gyroresonant acceleration with the turbulence produced by the isotropization of captured protons and charged dust. Protons are accelerated stochastically by the dust-induced turbulence. Analytical solutions for the proton and electron energy distributions are obtained and used to calculate broadband synchrotron emission. The solutions are compared with generic spectral behavior of blazars and gamma-ray bursts. Dust captured by a blast wave can generate turbulence that could accelerate protons to very high energies.

Sari, Piran, and Narayan have derived analytic formulas to model the spectra from gamma-ray burst blast waves that are energized by sweeping up material from the surrounding medium. We extend these expressions to apply to general radiative regimes and to include the effects of synchrotron self-absorption. Electron energy losses due to the synchrotron self-Compton process are also treated in a very approximate way. The calculated spectra are compared with detailed numerical simulation results. We find that the spectral and temporal breaks from the detailed numerical simulation are much smoother than the analytic formulas imply, and that the discrepancies between the analytic and numerical results are greatest near the breaks and endpoints of the synchrotron spectra. The expressions are most accurate (within a factor of ~ 3) in the optical/X-ray regime during the afterglow phase, and are more accurate when epsilon_e, the fraction of swept-up particle energy that is transferred to the electrons, is <~ 0.1. The analytic results provide at best order-of-magnitude accuracy in the self-absorbed radio/infrared regime, and give poor fits to the self-Compton spectra due to complications from Klein-Nishina effects and photon-photon opacity.

The blast-wave model for gamma-ray bursts (GRBs) has been called into question by observations of spectra from GRBs that are harder than can be produced through optically thin synchrotron emission. If GRBs originate from the collapse of massive stars, then circumstellar clouds near burst sources will be illuminated by intense gamma radiation, and the electrons in these clouds will be rapidly scattered to energies as large as several hundred keV. Low-energy photons that subsequently pass through the hot plasma will be scattered to higher energies, hardening the intrisic spectrum. This effect resolves the "line-of-death" objection to the synchrotron shock model. Illuminated clouds near GRBs will form relativistic plasmas containing large numbers of electron-positron pairs that can be detected within ~ 1-2 days of the explosion before expanding and dissipating. Localized regions of pair annihilation radiation in the Galaxy would reveal past GRB explosions.

The hard X-ray and gamma-ray phenomenology of gamma-ray bursts (GRBs) can be explained by an external shock model where a single relativistic blast wave interacts with the surrounding medium. Besides reproducing the generic spectral behavior of GRB profiles, the external shock model provides quantitative fits to the peak flux distribution, the > 1 s t_{50} duration distribution, and the distribution of the peaks E_pk of the vF_v spectra of GRBs measured with BATSE. The apparent paradox between a relativistic beaming scenario and the empirical finding that E_pk values are preferentially measured within the triggering range of a GRB detector is resolved by this model when blast wave physics and detector triggering criteria are taken into account. Some surprising implications follow, namely that the fireball event rate is ~1 per 10^4 years per Milky Way galaxy for unbeamed sources, and proportionally more if fireball outflows are collimated. This is ~ 3 orders of magnitude larger than normally quoted. Most of the clean and dirty fireball transients are undetected due to telescope sensitivity and design limitations. Strongly variable GRB time histories with good radiative efficiencies are possible because of the strongly enhanced emissions when a blast wave interacts with density inhomogeneities located nearly along the line-of-sight to the observer. Arguments against short timescale variability in an external shock model are answered, and difficulties in an internal shock/colliding shell model are mentioned.

Gamma-ray burst (GRB) engines are probed most intimately during the prompt gamma-ray luminous phase when the expanding blast wave is closest to the explosion center. Using GRBs 990123 and 940217 as guides, we briefly review observations of high-energy emission from GRBs and summarize some problems in GRB physics. \gamma\gamma transparency arguments imply relativistic beaming. The parameters that go into the external shock model are stated, and we show numerical simulation results of gamma-ray light curves from relativistic blast waves with different amounts of baryon loading. A distinct component due to the synchrotron self-Compton process produces significant emission at GeV and TeV energies. Predictions for spectral and temporal evolution at these energies are presented for a blast wave expanding into uniform surroundings. Observations of the slow decay of GeV-TeV radiation provide evidence for ultra-high energy cosmic ray acceleration in GRBs.

We present simple analytical expressions for the predicted spectral and temporal behavior of the early afterglow radiation from gamma-ray bursts in radiative regimes intermediate between the adiabatic and the fully radiative solutions of the blastwave hydrodynamic equations. Our expressions are valid as long as the relativistic electrons responsible for the observed synchrotron emission are in the fast cooling regime and the blast wave is relativistic. We show that even a slight deviation from a perfectly adiabatic evolution results in significant changes of the temporal characteristics of the afterglow emission.

We present detailed calculations of nonthermal synchrotron and synchrotron self-Compton (SSC) spectra radiated by blast waves that are energized by interactions with a uniform surrounding medium. Radio, optical, X-ray and gamma-ray light curves and spectral indices are calculated for a standard parameter set that yields hard GRB spectra during the prompt emission phase. Because no lateral spreading of the blast-wave is assumed, the calculated temporal breaks represent the sharpest breaks possible from collimated outflows in a uniform surrounding medium. Absence of SSC hardenings in observed GRB X-ray afterglows indicates magnetic field generation toward equipartition as the blast wave evolves. EGRET detections of 100 MeV-GeV photons observed promptly and 90 minutes after GRB 940217 are attributed to nonthermal synchrotron radiation and SSC emission from a decelerating blast wave, respectively. The SSC process will produce prompt TeV emission that could be observed from GRBs with redshifts $z \lesssim 0.1$, provided $\gamma$-$\gamma$ opacity in the source is small. Measurements of the time dependence of the 100 MeV-GeV spectral indices with the planned {\it GLAST} mission will chart the evolution of the SSC component and test the external shock scenario. Transient optical and X-ray emissions from misaligned GRBs are generally much weaker than on-axis emissions produced by dirty and clean fireballs that would themselves not trigger a GRB detector; thus detection of long wavelength transients not associated with GRBs will not unambiguously demonstrate GRB beaming.

We use a new method of analysis to determine parameters of cosmological gamma-ray bursts (GRBs), assuming that their distribution follows the star-formation history of the universe. Spectral evolution is calculated from an external shock model for fireball/blast wave evolution, and used to evaluate the measured peak flux, duration, and $\nu F_\nu$ peak photon energy for a GRB source occuring at a given redshift and with given values of total energy, baryon-loading and environmental parameters. We then fit model distributions of GRB sources to the observed peak flux, duration and $\nu F_\nu$ peak photon energy distributions. We find that the observed width of the $E_p$ and duration distributions can not be explained by cosmological redshift and time dilation effects. Rather, broad distributions of total blastwave energies and bulk Lorentz factors are necessary to explain the observed distributions simultaneously within the framework of our unifying GRB model. We discuss implications for source parameter distributions and determine a range of burst parameters consistent with the data.

In the external shock model, gamma-ray burst (GRB) emissions are produced by the energization and deceleration of a thin relativistic blast wave due to its interactions with the circumburst medium (CBM). We study the physical properties of an analytic function which describes temporally-evolving GRB spectra in the limit of a smooth CBM with density n(x)\propto x^(-\eta), where x is the radial coordinate. The hard-to-soft spectral evolution and the intensity-hardness correlation of GRB peaks are reproduced. We predict that (1) GRB peaks are aligned at high photon energies and lag at low energies according to a simple rule; that (2) temporal indices at the leading edge of a GRB peak display a well-defined shift with photon energy; and that (3) the change in the spectral index values between the leading and trailing edges of a GRB peak decreases at higher photon energies. The reason that GRBs are usually detected with vF_v peaks in the 50 keV - several MeV range for detectors which trigger on peak flux over a fixed time interval is shown to be a consequence of the inverse correlation of peak flux and duration of the radiation emitted by decelerating blast waves.

Blazar flare data show energy-dependent lags and correlated variability between optical/X-ray and GeV-TeV energies, and follow characteristic trajectories when plotted in the spectral-index/flux plane. This behavior is qualitatively explained if nonthermal electrons are injected over a finite time interval in the comoving plasma frame and cool by radiative processes. Numerical results are presented which show the importance of the effects of synchrotron self-Compton cooling and plasmoid deceleration. The use of INTEGRAL to advance our understanding of these systems is discussed.

The kinetic energy of bulk relativistic plasma ejected from the central engine of blazars is converted into nonthermal particle energy in the comoving frame through a process of sweeping up material from the surrounding medium. The resulting deceleration of the bulk plasma introduces a number of effects which must be included in blazar modeling. For example, the varying Doppler factor means that model fits must employ time integrations appropriate to the observing times of the detectors. We find that the ratio of the peak synchrotron fluxes reached at two different photon energies is largest when viewing along the jet axis, and becomes smaller at large angles to the jet axis. This effect is important in studies of the statistics of jet sources. Variability due either to bulk plasma deceleration or radiative cooling must be distinguished in order to apply recently proposed tests for beaming from correlated X-ray and TeV observations. The blast-wave physics developed to analyze these problems implies that most of the energy injected in the comoving frame is originally in the form of nonthermal hadrons. Because plasmoid deceleration can produce rapid variability due to a changing Doppler factor, arguments against hadronic blazar models related to the long radiative cooling time scale of hadrons are not compelling.

We have developed a computer model to calculate gamma ray burst (GRB) light curves and efficiencies from the interaction of a single, thin blast wave with clouds in the external medium. Large amplitude, spiky variability occurs when the clouds have radii r << R/Gamma, where R is the mean distance of a cloud from the GRB source and Gamma is the blast-wave Lorentz factor. Efficiencies >~ 10% require a large number of small clouds, each with sufficiently large column densities to extract most of the available blast-wave energy in the region of interaction. For such efficiencies, we find that erratic variability persists in the simulated GRB light curves. If GRB sources are surrounded by clouds with such properties, then short timescale variability of GRBs is possible in the external shock model.

In the cosmological blast-wave model for gamma ray bursts (GRBs), high energy (> 10 GeV) gamma-rays are produced either through Compton scattering of soft photons by ultrarelativistic electrons, or as a consequence of the acceleration of protons to ultrahigh energies. We describe the spectral and temporal characteristics of high energy gamma-rays produced by both mechanisms, and discuss how these processes can be distinguished through observations with low-threshold Cherenkov telescopes or GLAST. We propose that Compton scattering of starlight photons by blast wave electrons can produce delayed flares of GeV -- TeV radiation.

Most models for blazars and gamma-ray bursts involve relativistic plasma outflows powered by accretion processes onto black holes. The blast wave physics developed for cosmological models of GRBs is reviewed. Two points relevant for blazar modeling are made: (1) The injection of nonthermal relativistic particles in the comoving frame is simply treated though a process of energizing the plasma as it sweeps up material from the surrounding medium. (2) The primary energy source of blazar radiation derives from the bulk kinetic energy of the outflowing plasma. Thus deceleration of the plasma blast wave must be included in blazar flaring calculations, and this process will introduce temporal and spectral effects in addition to those produced by acceleration and radiative cooling.

If gamma-ray bursts are associated with dense star-forming regions in young galaxies, photoelectric absorption by the dense circumburster material (CBM) will occur. As the burst evolves, the surrounding material is photoionized, leading to fluorescence line emission and reduced photoelectric absorption opacity. We have analyzed this process in detail, accounting for the time-dependent photoelectric absorption, photoionization and fluorescence line emission from the CBM. We find that even if GRBs are hosted in dense star-forming regions, photoionization of the GRB environment leads to a constant, but very weak level of delayed fluorescence line emission on timescales of weeks to years after the burst. A temporally evolving iron~K edge absorption feature can serve as diagnostic tool to reveal the density structure of the CBM and may provide an opportunity for redshift measurements. We also investigated whether photoelectric absorption could be responsible for the spectral evolution of the low-energy slopes of some bright BATSE $\gamma$-ray bursts displaying extremely hard spectra below the peak energy, inconsistent with the optically-thin synchrotron shock model. We find that a very strong metal enrichment (~ 100 times solar-system abundances) in the gamma-ray burst environment and a rather peculiar spatial distribution of the CBM would be necessary in order to account for the observed hard spectra below a few 100 keV and their temporal evolution.

A simple function for the spectral power $P(\epsilon,t) \equiv \nu L(\nu) $ is proposed to model, with 9 parameters, the spectral and temporal evolution of the observed nonthermal synchrotron power flux from GRBs in the blast wave model. Here $\epsilon = h\nu/$m$_e$c$^2$ is the observed dimensionless photon energy and $t$ is the observing time. Assumptions and an issue of lack of self-consistency are spelled out. The spectra are found to be most sensitive to the baryon loading, expressed in terms of the initial bulk Lorentz factor $\Gamma_0$, and an equipartition term $q$ which is assumed to be constant in time and independent of $\Gamma_0$. Expressions are given for the peak spectral power $P_p(t) = P(\epsilon_p,t)$ at the photon energy $\epsilon = \epsilon_p(t)$ of the spectral power peak. A general rule is that the total fireball particle kinetic energy $E_0 \sim \Pi_0 t_d$, where $t_d \propto \Gamma_0^{-8/3}$ is the deceleration time scale and $\Pi_0 \equiv P(\epsilon_p,t_d) \propto \Gamma_0^{8/3}$ is the maximum measured bolometric power output in radiation, during which it is carried primarily by photons with energy ${\cal E}_0 = \epsilon_p(t_d) \propto q\Gamma_0^4$.

A simple model for variability in relativistic plasma outflows is studied, in which nonthermal electrons are continuously and uniformly injected in the comoving frame over a time interval dt. The evolution of the electron distribution is assumed to be dominated by synchrotron losses, and the energy- and time-dependence of the synchrotron and synchrotron self-Compton (SSC) fluxes are calculated for a power-law electron injection function with index s = 2. The mean time of a flare or pulse measured at photon energy E with respect to the onset of the injection event varies as E^{-1/2} and E^{-1/4} for synchrotron and SSC processes, respectively, until the time approaches the limiting intrinsic mean time (1+z)dt/(2 D), where z is the redshift and D is the Doppler factor. This dependence is in accord with recent analyses of blazar and GRB emissions, and suggests a method to discriminate between external Compton and SSC models of high-energy gamma radiation from blazars and GRBs. The qualititative behavior of the X-ray spectral index/flux relation observed from BL Lac objects can be explained with this model. This demonstrates that synchrotron losses are primarily responsible for the X-ray variability behavior and strengthens a new test for beaming from correlated hard X-ray/TeV observations.

It has recently been proposed that ultrahigh energy ($\gtrsim 10^{19}$ eV) cosmic rays (UHECRs) are accelerated by the blast waves associated with GRBs. We calculate the observed synchrotron radiation spectrum from protons and energetic leptons formed in the cascades initiated by photopion production, taking into account $\gamma\gamma$ attenuation at the source. Normalizing to the emission characteristics of GRB~970508, we predict $\sim 10$ MeV - 100 GeV fluxes at a level which may have been observed with EGRET from bright GRBs, and could be detected with the proposed GLAST experiment or with ground-based air \v Cerenkov telescopes having thresholds $\lesssim $ several hundred GeV. The temporal decay of the UHECR-induced high-energy $\gamma$-ray afterglows is significantly slower than that of the lower-energy burst and associated synchrotron self-Compton (SSC) radiation, which provides a direct way to test the hadronic origin of a high-energy GRB afterglow. Besides testing the UHECR origin hypothesis, the short wavelength emission and afterglows can be used to probe the level of the diffuse intergalactic infrared radiation field or constrain redshifts of GRB sources.

The synchrotron reflection scenario recently proposed to explain $\gamma$-ray flares observed from blazar jets is studied. Our analysis takes into account the angular distribution of the beamed radiation, the finite extent of the scattering region, and light travel-time effects. We compare energy densities and powers for synchrotron, SSC, reflected synchrotron (RSy), and external Compton (EC) scattering processes. If the width of the scattering layer is much larger than $\Gamma R^\prime_B$, where $\Gamma$ and $ R^\prime_B$ denote the bulk Lorentz factor and comoving-frame radius of the plasma blob, respectively, then the ratio of the RSy and synchrotron energy densities $\sim 4 \Gamma^3 n_{BLR} \sigma_{T} R^\prime_B$, where $n_{BLR}$ is the mean particle density in the broad line region (BLR). Our results imply that Thomson-thick scattering regions of narrow extent must be present for the synchrotron reflection mechanism to operate effectively. This process seems unlikely to cause flares in lineless BL Lac sources, where X-ray and TeV flares are common and the BLR is thought to be weak or absent. We sketch time profiles of flares for various scenarios, including a model where the blob is energized by sweeping up surrounding material.

We investigate the dynamics and radiation from a relativistic blast-wave which decelerates as it sweeps up ambient matter. The bulk kinetic energy of the blast-wave shell is converted into internal energy by the process of accreting external matter. If it takes the form of non-thermal electrons and magnetic fields, then this internal energy will be emitted as synchrotron and synchrotron self-Compton radiation. We perform analytic and numerical calculations for the deceleration and radiative processes and present time-resolved spectra throughout the evolution of the blast-wave. We also examine the dependence of the burst spectra and light curves on various parameters describing the magnetic field and non-thermal electron distributions. We find that for bursts such as GRB~910503, GRB~910601 and GRB~910814, the spectral shapes of the prompt gamma-ray emission at the peaks in $\nu F_\nu$ strongly constrain the magnetic fields in these bursts to be well below ($\la 10^{-2}$) the equipartition values. These calculations are also considered in the context of the afterglow emission from the recently detected gamma-ray burst counterparts.

An equation is derived to calculate the dynamics of relativistic magnetized plasma which decelerates by sweeping up matter from the ISM. Reduction to the non-radiative and radiative regimes is demonstrated. The evolving electron momentum distribution function in the comoving fluid frame is used to calculate the observed synchrotron radiation spectrum, assuming that a fixed fraction of the comoving proton power is instantaneously transformed into a power-law, "shock"-like electron distribution function. This permits an analytic solution for the case of a relativistic plasmoid with a constant internal magnetic field. Breaks in the temporal behavior of the primary burst emission and long wavelength afterglows occur on two time scales for external ISM density distributions. The first represents the time when the outflow sweeps up approximately M/G_0 of material from the ISM, where M is the mass in the outflow ejecta and G_0 is the initial Lorentz factor of the plasmoid. The second represents the time when synchrotron cooling begins strongly to regulate the number of nonthemal electrons that are producing radiation which is observed at a given energy. Results are applied to GRB and blazar variability. The slopes of the time profiles of the X-ray and optical light cuves of the afterglows of GRB 970228 and GRB 970508 are explained by injection of hard electron spectra with indices $s < 2$ in a deceleration regime near the radiative limit. We also propose that blazar flares are due to the transformation of the directed kinetic energy of the plasmoid when interacting with the external medium, as would occur if relativistic outflows pass through clouds orbiting the central supermassive black hole.

The flaring and fading radio, optical, and X-ray afterglows from GRB 970508 are modeled by a highly relativistic plasma sphere which decelerates by sweeping up ambient gas. The afterglow emission is assumed to be synchrotron radiation emitted by nonthermal electrons in the magnetized plasmoid. The temporal behavior of the delayed emission is controlled by the evolution of the Doppler factor and by adiabatic expansion losses of the nonthermal electrons in the plasmoid. Model fits to the optical data of GRB 970508 are provided, and the relative delay of the radio peak to the optical peak is found to result from a decrease in the observed self-absorption frequency as the plasmoid expands and decelerates. A variety of afterglow behaviors occurs for different observing angles and plasmoid parameters. The degree of collimation inferred from our fit to GRB 970508 implies a space density of GRB sources which exceeds the estimates from scenarios involving coalescing compact objects. This model can be verified through observations of superluminal motion in the delayed radio emission.

We present multiwavelength observations of the BL Lacertae object Markarian 501 (Mrk 501) in 1997 between April 8 and April 19. Evidence of correlated variability is seen in very high energy (VHE, E > 350 GeV) gamma-ray observations taken with the Whipple Observatory gamma-ray telescope, data from the Oriented Scintillation Spectrometer Experiment of the Compton Gamma-Ray Observatory, and quicklook results from the All-Sky Monitor of the Rossi X-ray Timing Explorer while the Energetic Gamma-Ray Experiment Telescope did not detect Mrk 501. Short term optical correlations are not conclusive but the U-band flux observed with the 1.2m telescope of the Whipple Observatory was 10% higher than in March. The average energy output of Mrk 501 appears to peak in the 2 keV to 100 keV range suggesting an extension of the synchrotron emission to at least 100 keV, the highest observed in a blazar and ~100 times higher than that seen in the other TeV-emitting BL Lac object, Mrk 421. The VHE gamma-ray flux observed during this period is the highest ever detected from this object. The VHE gamma-ray energy output is somewhat lower than the 2-100 keV range but the variability amplitude is larger. The correlations seen here do not require relativistic beaming of the emission unless the VHE spectrum extends to >5 TeV.

We report results of soft gamma-ray observations of the ultraluminous infrared galaxies Arp 220, Mrk 273, and Mrk 231 in order to test whether the infrared radiation from these sources originates from buried active galactic nuclei (AGNs). Only upper limits are measured, implying that the emergent soft gamma-ray luminosities are 1-2 orders of magnitude smaller than the infrared luminosities. Monte Carlo simulations of radiation transport through tori are used to infer the minimum column densities $N_H$ required to block transmission of soft gamma-rays from a buried AGN, assuming that spectra of AGNs in such sources are similar to those of radio-quiet quasars. Lack of measured gamma-ray emission provides no supporting evidence for the existence of buried AGNs in these galaxies, but is consistent with an origin of the infrared luminosity from starburst activity.

Two different model-independent mapping techniques have been applied to OSSE, SMM, TGRS and balloon data and reveal a feature in the 0.511 MeV electron-positron annihilation radiation pattern of our galaxy centered roughly at l=-2 deg. and b=10 deg. with a flux of 5x10^(-4) 0.511 MeV ph/cm^2/s. If near the galactic center, then positron sources are producing approximately 10^42 positrons/s which annihilate 1-2 kpc above the galactic plane. A starburst episode within the inner few hundred pc of our galaxy would drive hot pair-laden gas into the halo, with the one-sidedness pointing to the site of initial pressure release at the onset of the starburst activity. Positrons lose energy and annihilate as they are convected upward with the gas flow, and we calculate high-latitude annihilation patterns and fluxes in accord with the observations. Changes in the ionization state when the escaping gas cools could give annihilation radiation substructure. The fountain of hot (10^6-10^7 K) gas rising into the galactic halo would be seen through its enhanced dispersion measure, thermal emission, and recombination radiation.

The beaming pattern of thermal annihilation radiation is broader than the beaming pattern produced by isotropic nonthermal electrons and positrons in the jets of radio-emitting active galactic nuclei which Compton scatter photons from an external isotropic radiation field. Thus blueshifted thermal annihilation radiation can provide the dominant contribution to the high-energy radiation spectrum at observing angles theta > 1/Gamma, where Gamma is the bulk Lorentz factor of the outflowing plasma. This effect may account for the spectral features of MeV blazars discovered with the Compton Telescope on the Compton Gamma Ray Observatory. Coordinated gamma-ray observations of annihilation line radiation to infer Doppler factors and VLBI radio observations to measure transverse angular speeds of outflowing plasma blobs can be used to determine the Hubble constant.

The unusual properties of the bursting X-ray pulsar GRO J1744-28 are explained in terms of a low-mass X-ray binary system consisting of an evolved stellar companion transferring mass through Roche-lobe overflow onto a neutron star, implying that the inclination of the system is < 18 degrees. Interpretation of the QPO at frequency nu(QPO) = 40 Hz using the beat-frequency model of Alpar \& Shaham and the measured period derivative with the Ghosh \& Lamb accretion-torque model implies that the persistent X-ray luminosity of the source is approximately equal to the Eddington luminosity and that the neutron star has a surface equatorial magnetic field = 2e10 [40 Hz/nu(QPO)] G for standard neutron star parameters. This implies a distance to GRO J1744-28 of ~ 5 [nu(QPO)/40 Hz]^{1/6} b^{1/2} kpc, where b < 1 is a correction factor that depends on the orientation of the neutron star.

We show that Compton scattering of soft photons with energies near 100 eV in thermally stratified black-hole accretion plasmas with temperatures in the range 100 keV - 1 MeV can give rise to an X-ray spectral hardening near 10 keV. This could produce the hardening observed in the X-ray spectra of black holes, which is generally attributed to reflection or partial covering of the incident continuum source by cold optically thick matter. In addition, we show that the presence of very hot (kT=1 MeV) cores in plasmas leads to spectra exibiting high energy tails similar to those observed from Galactic black-hole candidates.

The broadband emission observed from radio galaxies, BL Lac objects and quasars is thought to be produced by energetic electrons confined in plasma blobs which are ejected from supermassive black holes at relativistic speeds. The lower energy radio-through-optical component is almost certainly nonthermal electron synchrotron radiation, while the higher-energy $\gamma$-ray component is probably due to Compton scattering of target photons by these same electrons. If the high-energy component is formed by jet electrons Compton-scattering radiation from outside the jet, such as the direct or rescattered accretion-disk photons, then the ratio of the power in the high-energy Compton component to the power in the low-energy synchrotron component stands in a simple relation depending on the observing angle to the jet and the plasma outflow speed. When combined with contemporaneous VLBI measurements of apparent transverse speed, we find that a simple diagram relates different classes of radio-loud active galactic nuclei (AGNs) and makes definite predictions for multiwavelength observations of these sources. If an equipartition parameter remains constant between episodes of plasma ejection, then multiple observations of a single source can in principle determine the Hubble constant.

We consider the stochastic acceleration of particles which results from resonant interactions with plasma waves in black hole magnetospheres. We calculate acceleration rates and escape time scales for protons and electrons resonating with Alfv\'en waves, and for electrons resonating with whistlers. Assuming either a Kolmogorov or Kraichnan wave spectrum, accretion at the Eddington limit, magnetic field strengths near equipartition, and turbulence energy densities $\sim 10\%$ of the total magnetic field energy density, we find that Alfv\'en waves accelerate protons to Lorentz factors $\lte 10^4$--$10^6$ before they escape from the system. Acceleration of electrons by fast mode and whistler waves can produce a nonthermal population of relativistic electons whose maximum energy is determined by a competition with radiation losses. Particle energization and outflow is not possible at lower accretion rates, magnetic field strengths, or turbulence levels due to dominant Coulomb losses. Increases in the accretion luminosity relative to the Eddington luminosity can trigger particle acceleration out of the thermal background, and this mechanism could account for the differences between radio-quiet and radio-loud active galactic nuclei. Observations of outflowing radio-emitting components following transient X-ray events in galactic X-ray novae and gamma-ray flares in blazars are in accord with this scenario.

We review models of cosmological gamma-ray bursts (GRBs). The statistical and $\gamma$-ray transparency issues are summarized. Neutron-star and black-hole merger scenarios are described and estimates of merger rates are summarized. We review the simple fireball models for GRBs and the recent work on non-simple fireballs. Alternative cosmological models, including models where GRBs are analogs of active galactic nuclei and where they are produced by high-field, short period pulsars, are also mentioned. The value of neutrino astronomy to solve the GRB puzzle is briefly reviewed.

We calculate the time- and energy-dependent opacity of high-energy gamma-rays attenuated by pair production interaction with accretion-flare photons that are scattered by gas and dust surrounding thenuclei of active galaxies. We show that the temporal behavior of the high-energy opacity cutoff can be used in conjunction with the time history of the accretion flare to determine the location of the gamma-ray emission site, as well as the column density and scale height of the material surrounding the central engine. Reverberation mapping using this technique is possible for nearby BL Lac objects and will be particularly valuable when the 10 -- 100 GeV gamma-ray window is opened, because there the attenuation of distant blazar radiation by the IIBR is negligible.

We use a Monte Carlo simulation to calculate the spectra of mildly relativistic thermal plasmas in pair balance. We use the exact integral expression for the electron-positron thermal annihilation spectrum, and provide accurate expressions for the Gaunt factors of electron-ion, electron-electron, and electron-positron thermal bremsstrahlung in the transrelativistic temperature regime. The particles are assumed to be uniformly distributed throughout a sphere, and the pair opacity is self-consistently calculated from the energy and angular distribution of scattered photons. The resultant photon spectra are compared with the nonrelativistic diffusion treatment of Sunyaev and Titarchuk, the bridging formulas of Zdziarski, and the relativistic corrections proposed by Titarchuk. We calculate allowed pair-balanced states of thermal plasmas with no pair escape which include bremsstrahlung and internal soft photons. The results are presented in the spectral index/compactness plane, and can be directly compared with observations of spectra from AGNs and Galactic black hole candidates. By comparing with X-ray spectral indices of Seyfert AGNs and compactnesses inferred from X-ray variability data, we find that the allowed solutions for pair equilibrium plasma imply that the temperatures of Seyfert galaxies are $\lesssim 300$ keV. This prediction can be tested with more sensitive gamma-ray observations of Seyfert galaxies. We find that if the X-ray variability time scale gives an accurate measure of the compactness, pair-dominated solutions are inconsistent with the data.

We propose that some of the unidentified Galactic plane (-10 degrees < b < 10 degrees) sources listed in the first EGRET catalog (Fichtel et al. 1994) are associated with nearby supernova remnants. The probability that the association of strong EGRET sources with supernova remnants is due to chance alignment is found to be < 1/2200. Three of the most convincing associations are radio-bright, nearby shell remnants and a fourth is associated with a plerion. We examine possible gamma-ray production mechanisms and conclude that pion production and decay from high-energy protons interacting with the remnant gas could produce the gamma rays. The energetics of this process are in accord with theories of cosmic ray production through supernova shock acceleration. The identification of pion-decay features in the gamma-ray spectra of supernova remnants would support this interpretation.

We show that gamma ray line emission at $\sim$0.4 MeV and $\lsim 0.2$ MeV can be produced by Compton scattering of beamed radiation in the jets of Galactic black hole candidates. This mechanism has the novel feature of not invoking the presence of $e^+$--$e^-$ pairs. To produce the two lines we employ a symmetric double sided jet with bulk flow velocity of about 0.5 c and incident beam radiation with a hard energy spectrum. We show that the two lines can be seen at viewing angle cosines relative to the jet ranging from 0.2 to 0.6. This comprises 40\% of the total solid angle. In addition, the line radiation is approximately 10\% polarized. Depending on the bulk flow and viewing angle the model can produce lines at other energies as well. In particular a broad feature near 1 MeV can be seen by viewing the jet close to its axis. Our model can also accommodate single line spectra if the beamed gamma ray emission or the jets themselves are asymmetric.

We consider whether the hard X-ray and soft gamma-ray emission from Centaurus A is beamed radiation from the active nucleus which is Compton-scattered into our line-of-sight. We derive the spectrum and degree of polarization of scattered radiation when incident beamed radiation is scattered from a cold ($kT<<m_ec^2$) electron cloud moving with bulk relativistic motion along the jet axis, and calculate results for an unpolarized, highly-beamed incident power-law photon source. We fit the OSSE data from Centaurus A with this model and find that if the scatterers are not moving relativistically, then the angle the jet makes with respect to our line-of-sight is $61^\circ\pm 5^\circ$. We predict a high degree of polarization of the scattered radiation below $\sim300$ keV. Future measurements with X-ray and gamma-ray polarimeters could be used to constrain or rule out such a scenario.

Compton scattering of external radiation by nonthermal particles in outflowing blazar jets is dominated by accretion-disk photons rather than scattered radiation to distances $\sim 0.01-0.1$ pc from the central engine for standard parameters, thus clarifying the limits of validity of the model by the present authors and the model of Sikora, Begelman, \& Rees. On the basis of contemporaneous Ginga X-ray and EGRET gamma-ray observations, we estimate the radius of 3C 279's gamma-ray photosphere to be smaller than estimated by Blandford. There is thus no need to require that the acceleration and emission sites of gamma-ray blazars to be located farther than $\sim 10^{2-3}$ gravitational radii from the central engine. We argue that lineless BL Lac objects, rather than quasars, are more likely to be detected in the TeV energy range.